Error when reading Gmsh mesh with overlapping physical groups

dokken · November 3, 2024, 10:02am

Note that you can supply tuples in integration measures, ie.

dx((5,7)) will integrate over all cells marked with either 5 or 7.

An example for a more complex mesh is shown

Wells-Group/TEAM30/blob/main/generate_team30_meshes_3D.py#L66


      
          
          # Copper wires is ordered in counter clock-wise order from angle = 0, 2*np.pi/num_segments...
          _domain_map_single = {
              "Cu": (7, 8),
              "Stator": (6,),
              "Rotor": (5,),
              "Al": (4,),
              "AirGap": (2, 3),
              "Air": (1,),
          }
          _domain_map_three = {
              "Cu": (7, 8, 9, 10, 11, 12),
              "Stator": (6,),
              "Rotor": (5,),
              "Al": (4,),
              "AirGap": (2, 3),
              "Air": (1,),
          }
          
          # Currents mapping to the domain marker sof the copper
          _currents_single = {7: {"alpha": 1, "beta": 0}, 8: {"alpha": -1, "beta": 0}}

github.com

Wells-Group/TEAM30/blob/main/team30_A_phi.py#L132


      
          sigma = fem.Function(DG0)
          density = fem.Function(DG0)
          for material, domain in domains.items():
              for marker in domain:
                  cells = ct.find(marker)
                  mu_R.x.array[cells] = model_parameters["mu_r"][material]
                  sigma.x.array[cells] = model_parameters["sigma"][material]
                  density.x.array[cells] = model_parameters["densities"][material]
          
          # Split domain into conductive and non-conductive parts
          Omega_n = domains["Cu"] + domains["Stator"] + domains["Air"] + domains["AirGap"]
          Omega_c = domains["Rotor"] + domains["Al"]
          
          sub_cells = np.hstack([ct.find(tag) for tag in Omega_c])
          sort_order = np.argsort(sub_cells)
          conductive_domain, sub_to_parent, _, _ = dolfinx.mesh.create_submesh(
              mesh, mesh.topology.dim, sub_cells[sort_order]
          )
          
          parent_cell_map = mesh.topology.index_map(mesh.topology.dim)
          num_cells = parent_cell_map.size_local + parent_cell_map.num_ghosts

See: Renumbered .msh imported data giving wrong subdomains - #4 by dokken

Topic		Replies	Views
A malloc() mismatch error when importing a GMSH mesh Errors	3	533	July 8, 2023
Issue reading the mesh tags imported from gmsh General mesh , boundary-conditions , dolfinx , mpi , gmsh	7	162	July 9, 2024
Reading xdmf mesh leads to illegal topological dimension mesh gmsh , xdmf	0	13	December 17, 2024
Meshio does not read Gmsh physical groups mesh	19	7439	February 16, 2022
Want to import already defined physical groups of Microstructure in Gmsh to FEnics mesh mesh	40	907	January 13, 2024

Error when reading Gmsh mesh with overlapping physical groups

Related topics