Hi, Folks -- I am trying to use the DofMap to impose Dirichlet-type boundary conditions within an explicit nonlinear mechanics solver which is embedded in a larger fluid-structure interaction solver. I constrain the velocity and force DOFs at Dirichlet boundaries to be zero by adding a DofConstraintRow for each of the constrained DOFs.
This seems to work well as long as I do not also have subdomains in my mesh; however, as soon as I start labeling different subdomains, the dof constraints no longer impose the boundary conditions. Is this expected behavior? Is there something special that I should be doing in the case of multiple subdomains? By the way, subdomain_id seems to be the most direct way to assign different material properties to different parts of the mesh. Are there any alternatives which I should use instead? Thanks! -- Boyce ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
