%matplotlib inline
A short example of how to calculate classical transport coefficients from Bragiński's theory.
from astropy import units as u
from plasmapy.formulary import ClassicalTransport
We'll use some sample ITER data, without much regard for whether the regime is even fit for classical transport theory:
thermal_energy_per_electron = 8.8 * u.keV
electron_concentration = 10.1e19 / u.m ** 3
thermal_energy_per_ion = 8.0 * u.keV
ion_concentration = electron_concentration
ion = "D+" # a crude approximation
braginskii = ClassicalTransport(
thermal_energy_per_electron,
electron_concentration,
thermal_energy_per_ion,
ion_concentration,
ion,
)
These variables are calculated during initialization and can be referred to straight away:
print(braginskii.coulomb_log_ei)
print(braginskii.coulomb_log_ii)
print(braginskii.hall_e)
print(braginskii.hall_i)
These quantities are not calculated during initialization and can be referred to via methods. To signify the need to calculate them, we call them via ().
print(braginskii.resistivity)
print(braginskii.thermoelectric_conductivity)
print(braginskii.electron_thermal_conductivity)
print(braginskii.ion_thermal_conductivity)
They also change with magnetization:
mag_braginskii = ClassicalTransport(
thermal_energy_per_electron,
electron_concentration,
thermal_energy_per_ion,
ion_concentration,
ion,
B=0.1 * u.T,
)
print(mag_braginskii.resistivity)
print(mag_braginskii.thermoelectric_conductivity)
print(mag_braginskii.electron_thermal_conductivity)
print(mag_braginskii.ion_thermal_conductivity)
all_direction_braginskii = ClassicalTransport(
thermal_energy_per_electron,
electron_concentration,
thermal_energy_per_ion,
ion_concentration,
ion,
B=0.1 * u.T,
field_orientation="all",
)
print(all_direction_braginskii.resistivity)
print(all_direction_braginskii.thermoelectric_conductivity)
print(all_direction_braginskii.electron_thermal_conductivity)
print(all_direction_braginskii.ion_thermal_conductivity)
The viscosities return arrays:
print(braginskii.electron_viscosity)
print(mag_braginskii.electron_viscosity)
print(braginskii.ion_viscosity)
print(mag_braginskii.ion_viscosity)