1

Display temperature and pressure change in 2D plot

Hi everyone,

 

I am trying to visualize pressure and temperature results for 3D model that has initial pressure and temperature values described by hydrostatic pressure and geothermal temperature gradient, respectively. Attached is a photo of temperature results at end of simulation using Paraview by exporting the excel file of the results into VTK file and open it using Paraview. I wonder if there is a way to illustrate change in temperature (T(@ time) - T (initial)) at any time.

 

Thanks in advance

3replies Oldest first
  • Oldest first
  • Newest first
  • Active threads
  • Popular
  • import toughio
    
    outputs = toughio.read_output("OUTPUT_ELEME.csv")
    T0 = 10.0 + 0.03 * outputs[0].data["Z"]
    for output in outputs:
      output.data["DTEMP"] = output.data["TEMP"] - T0
    toughio.write_output("OUTPUT_ELEME_DELTA.csv", outputs)
    

    Then, just use "toughio-export" with the new file.

    In this example, I assumed an initial geothermal gradient of 30 degrees C/km and a temperature at z = 0.0 m of 10 degrees C. It's probably different for your problem. If you have run a steady state simulation and use the SAVE file as INCON for your actual simulation, it would probably be better to import "T0" array from this SAVE file.

    Note that for the former case (linear), you can also use ParaView's built-in calculator to create a new data array.

    Like
    • Keurfon Luu My bad, the gradient should be negative since Z is negative:

      T0 = 10.0 - 0.03 * outputs[0].data["Z"]
      Like 1
  • I would like to know how to display the results in the surface plot (above) on logarithmic scale of X.

    Like
Like1 Follow
  • 1 Likes
  • 1 yr agoLast active
  • 3Replies
  • 115Views
  • 2 Following