—

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

This example plots the transport data of methane 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, C, CH, CH2, CH3, CH4, C2, C2H, C2H2, C2H3, C2H4, C2H5, C2H6,

• corresponding ions,

• electrons,

• no solid carbon species.

References data

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

• [Niu2016],

• [MinPlasCalc] (& [minplascalcPaper]),

Tags: transport methane CH4 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_CH4_Wu2016_Niu2016 = ThermoTransportDataReader(
    gas_name="CH4", pressure_atm=1, source="Wu2016_Niu2016", skip_missing_values=True
)
data_CH4_minplascalc = ThermoTransportDataReader(
    gas_name="CH4", pressure_atm=1, source="minplascalc"
)

Plot the thermal conductivity vs. temperature.

fig, ax = data_CH4_Wu2016_Niu2016.plot(
    x="T",
    y="kappa",
    show=False,
    label="Wu2016_Niu2016",
    ls="-",
    lw=4,
    color="black",
)
data_CH4_minplascalc.plot(
    x="T",
    y="kappa",
    fig_ax=(fig, ax),
    show=False,
    label="MinPlasCalc",
    ls="--",
    lw=3,
    color="red",
)
ax.legend()

plt.show()

Plot the electrical conductivity vs. temperature.

fig, ax = data_CH4_Wu2016_Niu2016.plot(
    x="T",
    y="sigma",
    show=False,
    label="Wu2016_Niu2016",
    ls="-",
    lw=4,
    color="black",
)
data_CH4_minplascalc.plot(
    x="T",
    y="sigma",
    fig_ax=(fig, ax),
    show=False,
    label="MinPlasCalc",
    ls="--",
    lw=3,
    color="red",
)
ax.legend()

plt.show()

Plot the dynamic viscosity vs. temperature.

fig, ax = data_CH4_Wu2016_Niu2016.plot(
    x="T",
    y="mu",
    show=False,
    label="Wu2016_Niu2016",
    ls="-",
    lw=4,
    color="black",
)
data_CH4_minplascalc.plot(
    x="T",
    y="mu",
    fig_ax=(fig, ax),
    show=False,
    label="MinPlasCalc",
    ls="--",
    lw=3,
    color="red",
)
ax.legend()

plt.show()