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

3 replies

null
    • Keurfon_Luu
    • 3 yrs ago
    • Reported - view
    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.

      • Keurfon_Luu
      • 3 yrs ago
      • Reported - view

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

      T0 = 10.0 - 0.03 * outputs[0].data["Z"]
    • Refaat_G_Hashish
    • 2 yrs ago
    • Reported - view

    I would like to know how to display the results in the surface plot (above) on logarithmic scale of X.

Content aside

  • 1 Likes
  • 2 yrs agoLast active
  • 3Replies
  • 156Views
  • 2 Following