developments for the nlte case for Argon and P1 radiation model.
Specifically, a new class called RadiativeDecay derived from Reaction has been added. This class deals with radiative decay reactions which are different from standard Arrhenius-type reactions. In RediativeDecay we calculate the get temperature and number densities as arguments, then we evaluate the number densities of the internal levels of a lumped species assuming Boltzmann distributed levels. Then, we calculate the escape factors for each transition and finally, we evaluate the effective Einstein A coefficient for the radiative reaction between lumped species. Changes also needed to be made in chemistry.cpp because the new class takes different arguments.
A new class called FieldInterpolator has been added. This is to convert the solution from a DG finite element collection (FEC) to an H1 FEC and vice versa. This is used in the solution of radiation transport and can be used also in the electromagnetics solver.
Developments for the P1 model. A new class was added in radiation.hpp. This is to handle radiation transport with the P1 model. The call for the P1 model is done in Mul2phys.cpp before the SolveStep. The used defines how after we solve the P1 equations. It does not support the evaluation of the radiation source term at internal rk substeps.
developments for the nlte case for Argon and P1 radiation model.
Specifically, a new class called RadiativeDecay derived from Reaction has been added. This class deals with radiative decay reactions which are different from standard Arrhenius-type reactions. In RediativeDecay we calculate the get temperature and number densities as arguments, then we evaluate the number densities of the internal levels of a lumped species assuming Boltzmann distributed levels. Then, we calculate the escape factors for each transition and finally, we evaluate the effective Einstein A coefficient for the radiative reaction between lumped species. Changes also needed to be made in chemistry.cpp because the new class takes different arguments.
A new class called FieldInterpolator has been added. This is to convert the solution from a DG finite element collection (FEC) to an H1 FEC and vice versa. This is used in the solution of radiation transport and can be used also in the electromagnetics solver.
Developments for the P1 model. A new class was added in radiation.hpp. This is to handle radiation transport with the P1 model. The call for the P1 model is done in Mul2phys.cpp before the SolveStep. The used defines how after we solve the P1 equations. It does not support the evaluation of the radiation source term at internal rk substeps.