Creating meshes from pyvista surfaces

Section author: Tobias Meisel (Helmholtz Centre for Environmental Research GmbH - UFZ)

For this example we create meshes from pyvista surfaces. See other examples for different meshing algorithms.

# Define a simple surface
from ogstools import meshlib as ml

bounds = (-200, 210, -200, 210)
args = {"bound2D": bounds, "amplitude": 100, "spread": 100, "n": 40}
surface1 = ml.Surface(ml.Gaussian2D(**args, height_offset=0), material_id=0)
surface2 = ml.Surface(ml.Gaussian2D(**args, height_offset=-100), material_id=1)
surface3 = ml.Surface(ml.Gaussian2D(**args, height_offset=-200), material_id=2)

ls = ml.LayerSet([ml.Layer(surface1, surface2), ml.Layer(surface2, surface3)])
mesh = ml.to_region_tetraeder(ls, 40).mesh
# interprets values as categories
mesh["regions"] = [str(m) for m in mesh["MaterialIDs"]]

Visualize the prism mesh

mesh.plot(scalars="regions", show_edges=True)

Static Scene

Interactive Scene