When executing a solver on the boundary of the bulk mesh (e.g. the free surface solver), the Dirichlet boundary conditions imposed in the Body Force section associated to this second body are not imposed. The problem is in the SolverUtils routine where the loop is only done over NumberOfBulkElements whereas it should also include boundary elements?
For example, using the free surface solver, setting
Code: Select all
Body Force 2 Zs = Real Zs0 Zs Condition = Equals Mask End
I think it should be made consistent such that Dirichlet condition can be imposed in the body force section for solvers executed on a boundary.