I created a a couple of physical groups within GMSH with appropriate labels/tags. I then tried to use the same labels after importing the .msh file in firedrake. I got the following error:

Please find below the relevant section in the .py file that generates the gmsh mesh and saves it in .msh file

gmsh.option.setNumber("Mesh.MshFileVersion", 2)

# Physical curves defined with appropriate tags starting from 1000

gmsh.model.addPhysicalGroup(1, b_o_c_arcs,1001)

gmsh.model.addPhysicalGroup(1, b_i_c_arcs,1002)

gmsh.model.addPhysicalGroup(1, list(t_o_c_arcs),1003)

gmsh.model.addPhysicalGroup(1, list(t_i_c_arcs),1004)

# Physical surfaces defined with appropriate tags starting from 2000

gmsh.model.addPhysicalGroup(2, circ_pl_surf,2001)

gmsh.model.addPhysicalGroup(2, extsurfs,2002)

# Physical volumes defined with appropriate tags starting from 3000

gmsh.model.addPhysicalGroup(3, voltags,3001)

gmsh.model.mesh.generate(3)

gmsh.write("Pipe3D.msh")

Here is the snippet of the .py file executed within firedrake:

mesh=Mesh('Pipe3D.msh')

bc1=DirichletBC(V,zero,1003)

When I execute the solve command I get this error:

Error: {1003} are not a valid markers (not in {2001,2002})

I think I have a workaround that is not yet giving me errors wherein I define a class derived from DirichletBC and redefine the @utils.cached_property nodes function. However, It would be great if I could do it using the existing DirichletBC class.

I would be very grateful if you could help me out here.

Thanks

Abhishek Venketeswaran

NETL Research Associate - ORISE

National Energy Technology Laboratory

Department of Energy

abhishek.venketeswaran@netl.doe.gov

Work: 412-386-4833

Mobile: 716-507-7890

www.netl.doe.gov