Problem with visualizing material type in ParaView

Refaat_G_Hashish
4 yrs ago
6replies

Hi every one,

I have run initialization simulation to obtain initial pressure and temperature conditions that follow hydrostatic pressure gradient and geothermal temperature gradient, respectively. The model has three material types; sand, well1, and well2. I have read the "INFILE" using Python and exported it to an unstructured VTK file and visualized it with ParaView. I could see the pressure and temperature results, however I couldn't see the material type. This is the code that I used to read the input file and export it to an unstructured VTK file.

import numpy
import toughio

parameters = toughio.read_input("INFILE")
points = numpy.array([v["center"] for v in parameters["elements"].values()])
incons = numpy.array([v["values"] for v in parameters["initial_conditions"].values()])

mesh = toughio.meshmaker.triangulate(points)
mesh.add_point_data("pressure", incons[:, 0])
mesh.add_point_data("temperature", incons[:, 3])
mesh.write("INFILE.vtu")

Thanks in advance

- Keurfon_Luu
- 4 yrs ago
- Reported - view
You also have to export material data (which are stored for each element under the key "material"):

rocks = numpy.array([v["material"] for v in parameters["elements"].values()]) mesh.add_point_data("material", rocks)

Note that because the mesh is reconstructed using Delaunay triangulation from a point cloud, material properties will be associated to point data "material". But you will also notice in Paraview another property called "material" associated to cell data as toughio always create a uniform cell data array "material" (as shown in your figure).
- Refaat_G_Hashish
  
  4 yrs ago
  
  Reported - view
  Keurfon Luu I have exported the material data, but it seems that the problem still exists!. I find the data range is [1,1] for the exported material!.
- Keurfon_Luu
  
  4 yrs ago
  
  Reported - view
  Refaat G Hashish Have you appended the code in my previous answer to your script?
  
  import numpy import toughio parameters = toughio.read_input("INFILE") points = numpy.array([v["center"] for v in parameters["elements"].values()]) incons = numpy.array([v["values"] for v in parameters["initial_conditions"].values()]) rocks = numpy.array([v["material"] for v in parameters["elements"].values()]) mesh = toughio.meshmaker.triangulate(points) mesh.add_point_data("pressure", incons[:, 0]) mesh.add_point_data("temperature", incons[:, 3]) mesh.add_point_data("material", rocks) mesh.write("INFILE.vtu")
  
  Running this script, I get a new point data array named "material" in ParaView:
- Refaat_G_Hashish
  
  4 yrs ago
  
  Reported - view
  Keurfon Luu I got it. Thanks!
- Keurfon_Luu
  
  4 yrs ago
  
  Reported - view
  Refaat G Hashish Anyway, with the latest version of toughio (v1.4.4), you can accurately reconstruct a mesh generated by TOUGH's MESHMAKER using the function toughio.meshmaker.voxelize.
  
  mesh = toughio.meshmaker.voxelize(points, origin=[0.0, 0.0, -50.0]) mesh.add_cell_data("pressure", incons[:, 0]) mesh.add_cell_data("temperature", incons[:, 3]) mesh.add_cell_data("material", rocks) mesh.write("INFILE.vtu")
- Refaat_G_Hashish
  
  4 yrs ago
  
  Reported - view
  Keurfon Luu I got it. Thanks again!