Example of submesh and DG

dokken · September 16, 2025, 5:34pm

Here you are hitting an issue that relates to: Discontinuity at interface using mixed domains - #2 by dokken
which has several workarounds required for coupling across the submeshes.

The support for this has been greatly improved in DOLFINx main branch, see for instance:

github.com/jorgensd/presentations

icosahom/code/emi_primal_single.py

main


      
          # Create integration measure for interface
          # Interior marker is considered as ("+") restriction
          def compute_interface_data(
              cell_tags: dolfinx.mesh.MeshTags, facet_indices: npt.NDArray[np.int32]
          ) -> npt.NDArray[np.int32]:
              """
              Compute interior facet integrals that are consistently ordered according to the `cell_tags`,
              such that the data `(cell0, facet_idx0, cell1, facet_idx1)` is ordered such that
              `cell_tags[cell0]`<`cell_tags[cell1]`, i.e the cell with the lowest cell marker is considered the
              "+" restriction".
          
              Args:
                  cell_tags: MeshTags that must contain an integer marker for all cells adjacent to the `facet_indices`
                  facet_indices: List of facets (local index) that are on the interface.
              Returns:
                  The integration data.
              """
              # Future compatibilty check
              integration_args: tuple[int] | tuple
              if Version("0.10.0") <= Version(dolfinx.__version__):

This file has been truncated. show original

where this coupling over interior submesh boundaries have been simplified greatly (copmared to my aforementioned implementation from 2024).

Topic		Replies	Views
Discontinuity at interface using mixed domains variational formulation dolfinx	12	761	September 13, 2024
Problem with nested submeshes dolfinx mesh , dolfinx	7	209	February 13, 2026
How generate polygonal meshes in DOLFIN-X? And how generate triangle submesh on each polygonal elements mesh	30	1114	May 29, 2024
Integrating over interface between two functionspaces dolfinx	11	380	December 25, 2024
How to define a Function in a submesh based on function values present in an adjoining submesh?	14	3181	December 10, 2024

Example of submesh and DG

Related topics