Hello, I am implementing a sequence of two simulations using a global mesh that contains both a fluid and a structure. My first simulation (elasticity equation) involves only the structure subdomain, and so I intend to run it on the structure submesh (by extracting it from the global mesh using GridGenerator::create_triangulation_with_removed_cells for instance), and the second simulation will be done on the entire mesh.
Before starting my second simulation, I will have to transfer the adaptive refinement from the structure submesh from my first simulation to the corresponding subdomain in the global mesh. This is the part that I am not sure how to implement. I am also using the p::d::triangulation class, which further complicates things, as the refinements may be split across multiple processors. The only idea I can come up with is to associate refinement with a physical point in the domain. After solving the elasticity eq-n, I would cycle through the coarsest level of a triangulation, and store the cell->center() point with the cell's refinement history into some file. For the second simulation on the global mesh, each cell would have to find its counterpart from the first simulation by making sure their center points match, and then load in the refinement. While possible, this seems a bit like an overkill. So I am curious if there is a better alternative than what I'm thinking of. Any suggestions? Thanks, Artur -- 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 dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
