A fitting formula for radiative cooling based on non-local thermodynamic equilibrium population from weakly-ionized air plasma

A fitting formula for radiative cooling with collisional-radiative population for air plasma flowfield has been developed. Population number densities are calculated from rate equations in order to evaluate the effects of nonequilibrium atomic and molecular processes. Many elementary processes are integrated to be applied to optically-thin plasmas in the number density range of 10¹²/cm³ ≤ N ≤ 10¹⁹/cm³ and the temperature range of 300 K ≤ T ≤ 40,000 K. Our results of the total radiative emissivity calculated from the collisional-radiative population are fitted in terms of temperature and total number density. To validate the analytic fitting formula, numerical simulation of a laser-induced blast wave propagation with the nonequilibrium radiative cooling is conducted and successfully reproduces the shock and plasma wave front time history observed by experiments. In addition, from the comparison between numerical simulations with the radiation cooling effect based on the fitting formula and those with a gray gas radiation model that assumes local thermodynamic equilibrium, we find that the displacement of the plasma front is slightly different due to the deviation of population probabilities. By using the fitting formula, we can easily and more accurately evaluate the radiative cooling effect without solving detailed collisional-radiative rate equations.