—

Plot transport data vs. temperature for a plasma of hydrogen.

This example plots the transport data of hydrogen as a function of temperature.

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

Equilibrium composition

The following species are considered in the mechanism:

• H, H2,

• corresponding ions,

• electrons.

References data

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

• [Boulos2023],

• [MinPlasCalc] (& [minplascalcPaper]),

• [Gueye2017] (thermal conductivity and electrical conductivity).

Tags: transport hydrogen H2 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_H2_minplascalc = ThermoTransportDataReader(
    gas_name="H2", pressure_atm=1, source="minplascalc"
)
data_H2_Gueye2017 = ThermoTransportDataReader(
    gas_name="H2", pressure_atm=1, source="Gueye2017", skip_missing_values=True
)

Plot the thermal conductivity vs. temperature.

fig, ax = data_H2_Boulos2023.plot(
    x="T",
    y="kappa",
    show=False,
    label="Boulos2023",
    ls="-",
    lw=4,
    color="black",
)
data_H2_minplascalc.plot(
    x="T",
    y="kappa",
    fig_ax=(fig, ax),
    show=False,
    label="MinPlasCalc",
    ls="--",
    lw=3,
    color="red",
)
data_H2_Gueye2017.plot(
    x="T",
    y="kappa",
    fig_ax=(fig, ax),
    show=False,
    label="Gueye2017",
    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="Boulos2023",
    ls="-",
    lw=4,
    color="black",
)
data_H2_minplascalc.plot(
    x="T",
    y="sigma",
    fig_ax=(fig, ax),
    show=False,
    label="MinPlasCalc",
    ls="--",
    lw=3,
    color="red",
)
data_H2_Gueye2017.plot(
    x="T",
    y="sigma",
    fig_ax=(fig, ax),
    show=False,
    label="Gueye2017",
    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="Boulos2023",
    ls="-",
    lw=4,
    color="black",
)
data_H2_minplascalc.plot(
    x="T",
    y="mu",
    fig_ax=(fig, ax),
    show=False,
    label="MinPlasCalc",
    ls="--",
    lw=3,
    color="red",
)
data_H2_Gueye2017.plot(
    x="T",
    y="mu",
    fig_ax=(fig, ax),
    show=False,
    label="Gueye2017",
    ls=":",
    lw=3,
    color="blue",
)
ax.legend()

plt.show()