—

Plot transport data vs. temperature for a plasma of H2, O2 or N2 in LTE.

This example plots the transport data of H2, O2 and N2 as a function of temperature.

The data for electrical conductivity, dynamic viscosity, and thermal conductivity are compared to reference data from the literature.

References data

The reference data for the thermal conductivity and electrical conductivity at 1 atm are taken from:

• [Boulos2023]

Tags: transport hydrogen H2 O2 N2 plasma thermal conductivity electrical conductivity viscosity

Import the required libraries.

import matplotlib.pyplot as plt
import seaborn as sns

from rizer.io.thermo_transport_data_reader import ThermoTransportDataReader

# Set the style of the plots.
sns.set_theme("talk")

Load reference data.

data_H2_Boulos2023 = ThermoTransportDataReader(
    gas_name="H2", pressure_atm=1, source="Boulos2023"
)
data_N2_Boulos2023 = ThermoTransportDataReader(
    gas_name="N2", pressure_atm=1, source="Boulos2023"
)
data_O2_Boulos2023 = ThermoTransportDataReader(
    gas_name="O2", pressure_atm=1, source="Boulos2023"
)

Plot the thermal conductivity vs. temperature.

fig, ax = data_H2_Boulos2023.plot(
    x="T",
    y="kappa",
    show=False,
    label=r"$\mathregular{H_2}$",
    ls="-",
    lw=4,
    color="black",
)
data_N2_Boulos2023.plot(
    x="T",
    y="kappa",
    fig_ax=(fig, ax),
    show=False,
    label=r"$\mathregular{N_2}$",
    ls="--",
    lw=3,
    color="red",
)
data_O2_Boulos2023.plot(
    x="T",
    y="kappa",
    fig_ax=(fig, ax),
    show=False,
    label=r"$\mathregular{O_2}$",
    ls=":",
    lw=3,
    color="blue",
)
ax.legend()

plt.show()

Plot the electrical conductivity vs. temperature.

fig, ax = data_H2_Boulos2023.plot(
    x="T",
    y="sigma",
    show=False,
    label=r"$\mathregular{H_2}$",
    ls="-",
    lw=4,
    color="black",
)
data_N2_Boulos2023.plot(
    x="T",
    y="sigma",
    fig_ax=(fig, ax),
    show=False,
    label=r"$\mathregular{N_2}$",
    ls="--",
    lw=3,
    color="red",
)
data_O2_Boulos2023.plot(
    x="T",
    y="sigma",
    fig_ax=(fig, ax),
    show=False,
    label=r"$\mathregular{O_2}$",
    ls=":",
    lw=3,
    color="blue",
)
ax.legend()

plt.show()

Plot the dynamic viscosity vs. temperature.

fig, ax = data_H2_Boulos2023.plot(
    x="T",
    y="mu",
    show=False,
    label=r"$\mathregular{H_2}$",
    ls="-",
    lw=4,
    color="black",
)
data_N2_Boulos2023.plot(
    x="T",
    y="mu",
    fig_ax=(fig, ax),
    show=False,
    label=r"$\mathregular{N_2}$",
    ls="--",
    lw=3,
    color="red",
)
data_O2_Boulos2023.plot(
    x="T",
    y="mu",
    fig_ax=(fig, ax),
    show=False,
    label=r"$\mathregular{O_2}$",
    ls=":",
    lw=3,
    color="blue",
)
ax.legend()

plt.show()