How to read and visualize MESH file using PyVista
I have created a Cartesian grid model using Mesh Maker which is given the attached MESH file. How to read MESH file and visualize the grid model using PyVista?
12 replies
-
Hi Refaat G Hashish ,
Generally, given any MESH file, you can use the function `toughio.meshmaker.triangulate`. As you guess from the name, this function performs a Delaunay triangulation to reconstruct a mesh from a cloud of points (here, the X, Y, Z coordinates of the elements in the MESH file), so it not the actual mesh generated by MESHMAKER.
import toughio parameters = toughio.read_input("MESH") points = [v["center"] for v in parameters["elements"].values()] mesh = toughio.meshmaker.triangulate(points) mesh.plot(show_edges=True)Otherwise, you can also use the function `toughio.meshmaker.from_meshmaker` to generate a mesh from your TOUGH input file that you used to generate your mesh with MESHMAKER (a file that contains a MESHMAKER block).
mesh = toughio.meshmaker.from_meshmaker("INFILE") mesh.plot(show_edges=True)Alternatively, why not simply use the function `toughio.meshmaker.structured_grid` or `toughio.meshmaker.cylindric_grid` to generate your mesh (and MESH file via the method `mesh.write_tough`)?
import numpy as np dx = np.arange(5) dy = np.arange(6) dz = np.arange(7) mesh = toughio.meshmaker.structured_grid(dx, dy, dz) mesh.write_tough() mesh.plot(show_edges=True)
-
Keurfon Luu Thank you so much. All of the methods are working fine. I wonder why only a single material is shown on the grid model while the MESH file includes grid cells with different materials. Is there a way to show all the materials on the grid model?
-
Refaat G Hashish I am not sure how you got material defined in your MESH file since TOUGH built-in MESHMAKER does not have such feature.
Anyway, none of the previous sample codes have imported the material properties from your MESH file.
import toughio parameters = toughio.read_input("MESH") points = [v["center"] for v in parameters["elements"].values()] materials = [v["material"] for v in parameters["elements"].values()] mesh = toughio.meshmaker.triangulate(points) mesh.add_point_data("material", materials) mesh.plot(scalars="material", show_edges=True)If you use triangulation, you should use the method `add_point_data`, otherwise you should use `set_material`.
-
Keurfon Luu Thank you. It now shows the material types on the grid model. By the way, I have assigned the material types to the MESH file after creating it using MESHMAKER using the following python code:
------------------------------------------------------------------------------------------------------
import toughio
# Read file with block ELEME (e.g. MESH)
parameters = toughio.read_input("MESH")# Loop over elements
for label, element in parameters["elements"].items():
x, y, z = element["center"]# Assign material name depending on element's depth
if 0.0 >= z > -10.0:
element["material"] = "mat01"
elif -10.0 >= z > -20.0:
element["material"] = "mat02"
elif -20.0 >= z > -30.0:
element["material"] = "mat03"
# and so on# Write your file back with modified materials
toughio.write_input("MESH", parameters)--------------------------------------------------------------------------------------------------------
-
Refaat G Hashish Since you are using toughio to set the material, why not directly use toughio built-in mesh generation functions to generate your structured mesh?
dx = 10 * [10.0] mesh = toughio.meshmaker.structured_grid(dx, dx, dx, material="MID") mesh.set_material("BOT", mesh.filter.box(z0=-100, dz=20)) mesh.set_material("TOP", mesh.filter.box(z0=-20, dz=20)) mesh.plot(show_edges=True)I guess that's on me and the missing documentation about all these features.
-
Keurfon Luu Thank you. I have a final question. Is there a way to modify the background color?
-
Refaat G Hashish The method `plot` relies on PyVista, so you can refer to PyVista's documentation to customize your plots. For instance, to change the color of your background, see here: https://docs.pyvista.org/api/plotting/_autosummary/pyvista.themes.ParaViewTheme.background.html?highlight=background+color
-
Keurfon Luu I got the following error message.
Traceback (most recent call last):
File "<ipython-input-11-4d803777092a>", line 1, in <module>
pyvista.global_theme.background = 'white'AttributeError: module 'pyvista' has no attribute 'global_theme'
-
Keurfon Luu Got it. Thank you.
-
Keurfon Luu Is it possible to use the function "toughio.meshmaker" to reconstruct a mesh from a cloud of points such that grid cells are rectangular?
-
Refaat G Hashish Generally speaking, no. To reconstruct an actual mesh for visualization, you need the corner points coordinates of the cells which are lacking in TOUGH generated MESH files (X, Y, Z are the coordinates of cell centers).
However, in case your mesh is actually a structured grid, you can try to voxelize it using `toughio.meshmaker.voxelize` (https://keurfonluu.github.io/toughio/api/mesh.html?highlight=voxelize#toughio.meshmaker.voxelize) which will try to reconstruct the original mesh if you provide the coordinate of the lower left corner point. The only issue is that the X, Y, Z coordinates in MESH files generated by TOUGH's built-in MESHMAKER usually are not accurate enough (only 3 decimal places) and the algorithm will fail because some elements have the same coordinates.
-
If I use this method to generate MESH, how can I write that only one center column of cells (say the center column in an 11*11*11 model) is assigned to wellbore cells in T2well? So basically, I need to add those cells as another material called (for example) wella?
-
