1-D Geometry and Meshing

Hello,

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)

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

Hi,

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.