[Libmesh-users] imposing boundary conditions inside the domain

[Yu Z. <eli...@li...>]

Dear all,

For a 2D problem, I'm trying to impose Dirichlet BC on a physical boundary connecting two subdomains in a full domain (containing more than 2 subdomains). With the penalty method, I tried with

c.has_side_boundary_id(ID)

to find this physical boundary, where ID is the physical curve ID defined in gmsh. However, it does not work. This method seems only working with outer boundaries instead of those within the full domain.

I went through the examples, and found the Miscellaneous Example 9 provides a solution for this issue, where two different meshes are created in two subdomains. Then the boundary connecting these two subdomains can be found as those outer boundaries for a single domain. This method would ask for heavy works on the mesh generation side when there are many such interior B.C. needed within the full domain.

I just wonder whether libmesh has any routine that can automatically find the interior edges with specific IDs, such that specific B.C. can be imposed onto these edges.

Thank you very much.

Best regards,
Eli