.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "examples/heat_pipe_in_cylindrical_geometry/2_postprocessing.py"
.. LINE NUMBERS ARE GIVEN BELOW.

.. only:: html

    .. note::
        :class: sphx-glr-download-link-note

        :ref:`Go to the end <sphx_glr_download_examples_heat_pipe_in_cylindrical_geometry_2_postprocessing.py>`
        to download the full example code

.. rst-class:: sphx-glr-example-title

.. _sphx_glr_examples_heat_pipe_in_cylindrical_geometry_2_postprocessing.py:


Plot profiles
=============

The objective of this example is to import the output results and plot the profiles of temperature, pressure, liquid saturation and air mass fraction.

.. GENERATED FROM PYTHON SOURCE LINES 10-11

Here, we assume that the simulation's output have been written in the file "OUTPUT". To import the results, we use the function :func:`toughio.read_output`. The variable `outputs` is a list with three :class:`toughio.Output` corresponding to the three time steps requested in the preprocessing example. In this example, we want to look at the last time step (index -1).

.. GENERATED FROM PYTHON SOURCE LINES 11-23

.. code-block:: Python


    import numpy as np
    import toughio

    outputs = toughio.read_output("OUTPUT")
    output = outputs[-1]

    t = output.data["T"]
    p = output.data["P"]
    sl = output.data["SL"]
    xm = output.data["XAIRG"]








.. GENERATED FROM PYTHON SOURCE LINES 26-27

It is well known that for the stated conditions (1-D radial geometry, homogeneous medium, uniform initial conditions, and a constant-rate line source) the problem has a similarity solution: The partial differential equations for this complex two-phase flow problem can be rigorously transformed into a set of ordinary differential equations in the variable :math:`Z = R/\sqrt{t}`, which can be easily solved to any degree of accuracy desired by means of one-dimensional numerical integration (O'Sullivan, 1981). Comparison of TOUGH2 simulations with the semi-analytical similarity solution has shown excellent agreement (Doughty and Pruess, 1992). To define such variable, we first need to import the mesh 

.. GENERATED FROM PYTHON SOURCE LINES 27-31

.. code-block:: Python


    mesh = toughio.read_mesh("mesh.pickle")
    R = np.log(mesh.centers[:, 0] / (output.time) ** 0.5)








.. GENERATED FROM PYTHON SOURCE LINES 34-35

Now that the required data have been imported, we can plot the results using :mod:`matplotlib`.

.. GENERATED FROM PYTHON SOURCE LINES 35-62

.. code-block:: Python


    import matplotlib.pyplot as plt
    plt.rc("font", size=12)

    fig, ax1 = plt.subplots(figsize=(8, 5))
    ax2 = ax1.twinx()

    ax1.plot(R, t, color="black", linestyle="--", label="Temperature")
    ax1.set_ylim(0.0, 260.0)
    ax1.set_ylabel("Temperature ($\degree$C)")

    ax2.plot(R, sl, label="Liquid saturation")
    ax2.plot(R, xm, label="Air mass fraction")
    ax2.plot(R, p * 1.0e-5, label="Pressure (bar)")
    ax2.set_ylim(0.0, 1.3)
    ax2.set_ylabel("Liquid saturation, air mass fraction, pressure")

    ax1.set_xlim(R.min(), -4.0)
    ax1.set_xlabel("$ln(R/\sqrt{t})$")

    fig.legend(
        loc="lower left",
        bbox_to_anchor=(0.0, 0.0),
        bbox_transform=ax1.transAxes,
        frameon=False,
    )




.. image-sg:: /examples/heat_pipe_in_cylindrical_geometry/images/sphx_glr_2_postprocessing_001.png
   :alt: 2 postprocessing
   :srcset: /examples/heat_pipe_in_cylindrical_geometry/images/sphx_glr_2_postprocessing_001.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <matplotlib.legend.Legend object at 0x7f70a6822460>




.. rst-class:: sphx-glr-timing

   **Total running time of the script:** (0 minutes 0.212 seconds)


.. _sphx_glr_download_examples_heat_pipe_in_cylindrical_geometry_2_postprocessing.py:

.. only:: html

  .. container:: sphx-glr-footer sphx-glr-footer-example

    .. container:: sphx-glr-download sphx-glr-download-jupyter

      :download:`Download Jupyter notebook: 2_postprocessing.ipynb <2_postprocessing.ipynb>`

    .. container:: sphx-glr-download sphx-glr-download-python

      :download:`Download Python source code: 2_postprocessing.py <2_postprocessing.py>`


.. only:: html

 .. rst-class:: sphx-glr-signature

    `Gallery generated by Sphinx-Gallery <https://sphinx-gallery.github.io>`_