Dear all, I have a triangulation which has been refined along some (polygonal) patch, thus forming two non-overlapping domains (say, domain #1 and #2). My task is to extract cells for domain #1, and solve BVP for this domain (with known boundary conditions). The rest of the cells (domain #2) are not needed. An example of such situation is shown attached (with land and ocean/sea being two domains).
FWIW, the problem is posed in dim = 2, spacedim = 3. The solution is in H1. Here is what I struggle with: - I thought of using GridGenerator::create_triangulation_with_removed_cells, but this one does not accept refined non-conforming meshes. - My next thought was to detect those edges, which domains #1 and #2 share and calculate face integrals on them during assembly, thus incorporating BCs. Will DoFTools::make_flux_sparsity_pattern be the right function for this? - Next, what should I do on cells in domain #2? Logically, I should somehow treat them with FE_Nothing... But how to get this done right? Thanks! Alexander -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/6cdc5097-d5ac-4bb3-b229-dc9c74c1c3e7%40googlegroups.com. For more options, visit https://groups.google.com/d/optout.
