Braginskii coefficients

A short example of how to calculate classical transport coefficients from Bragiński’s theory.

from astropy import units as u
from plasmapy.physics.transport.braginskii import ClassicalTransport
import matplotlib.pyplot as plt
import numpy as np

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_particle = 'D+'  # a crude approximation

We now make the default ClassicalTransport object:

braginskii = ClassicalTransport(thermal_energy_per_electron,
                                electron_concentration,
                                thermal_energy_per_ion,
                                ion_concentration,
                                ion_particle)

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)

Out:

18.015542122021248
20.415575220593944
0.0
0.0

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 (). They could be made to act like variables via @Property, but I’m not sure that’s wise.

print(braginskii.resistivity())
print(braginskii.thermoelectric_conductivity())

Out:

1.1541382703097203e-09 m Ohm
0.7110839986207994

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