Remeshing with triangle elements

Section author: Florian Zill (Helmholtz Centre for Environmental Research GmbH - UFZ)

This short example showcases the function remesh_with_tri which allows us to take an existing mesh and re-discretize it with triangle elements. This is useful for models, where the underlying meshing script is not available or hard to adapt.

from pathlib import Path
from tempfile import mkdtemp

import ogstools as ogs
from ogstools import examples
from ogstools.msh2vtu import msh2vtu

mesh = examples.load_meshseries_THM_2D_PVD().mesh(1)

elements.

fig = mesh.plot_contourf(ogs.variables.material_id)

Here, we do the remeshing and use msh2vtu to convert the resulting msh file to an OGS-compatible vtu file.

mesh = examples.load_meshseries_THM_2D_PVD().mesh(1)
temp_dir = Path(mkdtemp())
msh_path = temp_dir / "tri_mesh.msh"
ogs.meshlib.gmsh_meshing.remesh_with_triangle(mesh, msh_path)
msh2vtu(msh_path, temp_dir, reindex=False, log_level="ERROR")
mesh = ogs.Mesh(temp_dir / "tri_mesh_domain.vtu")
fig = mesh.plot_contourf(ogs.variables.material_id)