1-D Geometry and Meshing


I am trying to create a similar geometry as described here (Precipitation/dissolution equilibrium reactions in a saturated column).

I am following the python files given in the video tutorial (https://www.opengeosys.org/docs/tools/getting-started/video-tutorial).

When I check the solute transport benchmark, I see two different vtu files for boundaries, upstream and downstream. After creating the geometry, I use msh2vtu.py file (given in the tutorial) to convert my gmsh file to vtu. But it only makes one single file for boundary (mesh_MICP_bondary.vtu)

I try to assign a physical group for the left boundary (Point 1) and right boundary (Point 2), but I get an error like this. How can I get two different files for the left (influent or upstream) and right boundary (effluent or downstream)?

You can also find the files in the attached zip file.

MICP.zip (33.6 KB)

I didn’t have time to go through the files yet, but maybe the developer can help? Or is it already solved? @dominik-kern

As I assumed mainly meshes in 2D and 3D, I thought OGS thinks so too, and was lazy about 1D meshes. The physical groups are cell (element) properties but for a boundary in 1D there is no difference between point (node) and cell (element) as it is just a point. So I have not yet worked this out, I will put it on my TODO list. As workaround I propose open boundary.vtu with paraview, extract one point save it as A.vtu and extract the other point, again save it as B.vtu and use these vtu-files.

Thank you for pointing out this bug.


I actually stop doing this with msh2vtu.py and started to use the generateStructuredMesh utility of OGS6.

I am glad to point out this bug accidentally. I hope you will solve it.