The differential equation for the kinetic temperature profile.
Currently only takes into account adiabatic expansion.
This function is used by EnergyBalance to iterate and solve the
ODE.
- Parameters:
T (float) - The temperature at which to evaluate the differential equation
r (array/float) - The radial point(s) (cm)
v (Velocity()) - The velocity profile object
gamma (Profiler()) - The adiabatic coefficient profile as function of T
rates (Profiler()) - The heating and cooling rates, summed up (erg/s/cm3, H-C).
Includes the density factor, without the velocity component. This
speeds up iteration, limiting the radial grid evaluation where
possible. Default if only adiabatic cooling is taken into
account.
(default: None)
warn (bool) - Warn when extrapolation occurs in an interpolation object. Not
applicable to functional evaluation, or gamma in which case the
extrapolation is deemed safe (constant value at low and high end
temperatures consistent with diatomic gases)
(default: 1)
- Returns: array/float
- The derivative with respect to radius (K/cm)
|