Creating meshes from vtu surface files

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

For this example we create meshes from surface layers.

from ogstools import meshlib as ml
from ogstools.examples import surface_paths

The loaded surfaces are defined within VTU files and adhere to properties such as non-intersecting boundaries with consistent x and y bounds. Alternatively, surfaces can also be created using PyVista with the same properties.

surface1 = ml.Surface(surface_paths[0], material_id=0)
surface2 = ml.Surface(surface_paths[1], material_id=5)
surface3 = ml.Surface(surface_paths[2], material_id=2)
surface4 = ml.Surface(surface_paths[3], material_id=3)

Create 3 layers from previously defined surfaces and add all layers to a layerset (ordererd from top to bottom)

layer1 = ml.Layer(top=surface1, bottom=surface2, num_subdivisions=2)
layer2 = ml.Layer(top=surface2, bottom=surface3, num_subdivisions=1)
layer3 = ml.Layer(top=surface3, bottom=surface4, num_subdivisions=0)

layer_set1 = ml.LayerSet(layers=[layer1, layer2, layer3])

From layerset creation of simplified meshes (sm), prism meshes (pm), voxel meshes (vm), tetraeder mesh (tm) is possible.

sm = ml.to_region_simplified(layer_set1, xy_resolution=200, rank=3).mesh
pm = ml.to_region_prism(layer_set1, resolution=200).mesh
vm = ml.to_region_voxel(layer_set1, resolution=[200, 200, 50]).mesh
tm = ml.to_region_tetraeder(layer_set1, resolution=200).mesh

Simplified mesh

%%

sm["regions"] = [str(m) for m in sm["MaterialIDs"]]
sm.scale([1, 1, 5]).plot(scalars="regions", show_edges=True)

Static Scene

Interactive Scene

Voxel mesh

%%

vm["regions"] = [str(m) for m in vm["MaterialIDs"]]
vm.scale([1, 1, 5]).plot(scalars="regions", show_edges=True)

Static Scene

Interactive Scene

Prism mesh

%%

pm["regions"] = [str(m) for m in pm["MaterialIDs"]]
pm.scale([1, 1, 5]).plot(scalars="regions", show_edges=True)

Static Scene

Interactive Scene

Tetraeder mesh

%%

tm["regions"] = [str(m) for m in tm["MaterialIDs"]]
tm.scale([1, 1, 5]).plot(scalars="regions", show_edges=True)

Static Scene

Interactive Scene