Hi,
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