FEniCS / dolfinx

Next generation FEniCS problem solving environment
https://fenicsproject.org
GNU Lesser General Public License v3.0
732 stars 178 forks source link

Support periodic meshes via mesh topology/geometry #18

Open garth-wells opened 6 years ago

garth-wells commented 6 years ago

Periodic boundary conditions are currently handled via the dofmap construction. This has some drawbacks:

It would be simpler to support periodic meshes in which the topology naturally accounts for the periodicity.

jedbrown commented 4 years ago

Thanks to @haplav for linking us to this. We're having a similar discussion over at PETSc. Xdmf does not support periodic, but we're interested in extending it (and specifying an HDF5 format from which the XML can automatically be created). As a modest extension, I would propose a dataset of shape (n, 2) containing (image_node, true_node). Perhaps @knepley would prefer (image_element, face_num, true_elem, face_num, rotation), but I think that's quite a bit more complexity/specification for an existing Xdmf user.

I'd be curious to hear more details about "topology naturally accounts", because it's necessary to be able to apply some transformations when mapping true -> image point (such as rotation in the linked example).