Hi everyone,
I am using Moose Framework I have been working on getting a mesh that
can be read by libmesh, with 2D fractures in a 3D rock matrix.
I am using the .msh format for my mesh. The gmsh reader of Moose wasn't
able to recognize those 2D blocks in a 3D domain as it only assigned
them as sidesets.
So I changed a bit the file gmsh_io.C and it's now working for me. In
order to do that however, I had to add a keyword in my Gmsh file in
order to treat it slightly differently in the gmsh_io.C reader function.
My question to the list is therefore simple: which solution do you
recommend to make this work in the proper way?
My current hack is to add an extra keyword in the .msh file in the
$PhyisicalNames paragraph, at the end of the line defining the tag that
we want to consider as a lower dimensional block (instead of a sideset).
Since the gmsh_io.C reader function doesn't currently use that
information, it's not breaking anything. Also, Gmsh itself doesn't seem
to bother and can still read the file. But I understand that it might
not be the best solution since it changes (very slightly) the official
Gmsh format...
What do you recommend?
For info, here's what an example of my .msh section looks like
$PhysicalNames
3
2 1 "fractures" MYKEYWORD
2 2 "a_normal_sideset"
3 3 "a_normal_block"
$EndPhysicalNames
Here is the link to the wiki I wrote about how to use and install my
hack : https://github.com/pou036/redback/wiki/Gmsh%20reader
By the way, what is the proper way to compile a single file in the
libmesh subdirectory?
Cheers,
Martin Lesueur
------------------------------------------------------------------------------
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
