-
Multi-physics modeling of non-equilibrium phenomena in inductively coupled plasma discharges: Part II. Multi-temperature approach
Authors:
Sanjeev Kumar,
Alessandro Munafo,
Sung Min Jo,
Marco Panesi
Abstract:
This paper provides a comparison between the vibrational-specific state-to-state (StS) model for nitrogen plasma elaborated in Part I of this work and conventional two-temperature (2-T) models for simulating inductively coupled plasma (ICP) discharges under non-Local Thermodynamic Equilibrium (NLTE) conditions. Simulations are performed within the multi-physics computational framework established…
▽ More
This paper provides a comparison between the vibrational-specific state-to-state (StS) model for nitrogen plasma elaborated in Part I of this work and conventional two-temperature (2-T) models for simulating inductively coupled plasma (ICP) discharges under non-Local Thermodynamic Equilibrium (NLTE) conditions. Simulations are performed within the multi-physics computational framework established for ICP in Part I. Based on the findings of Part I, the quasi-steady-state (QSS) assumption is validated in the plasma core, thereby enabling the calculation of global rate coefficients under this assumption. This facilitates the reduction of the StS model to a "consistent" macroscopic 2-T model. Results from the StS model for nitrogen ICP torch exhibit considerable discrepancies when compared against predictions from the widely utilized Park 2-T model. On the contrary, the comparison between the newly proposed 2-T model, consistently derived from the original vibronic StS model, and the full StS results demonstrate excellent agreement in terms of plasma core location, morphology, and peak temperature distributions. This demonstrates the ability of the proposed 2-T model to capture the energy transfer and reactive processes predicted by the comprehensive StS model. Additionally, the study identifies the vibrational-translational (VT) energy transfer term in the 2-T model as the predominant factor in dictating plasma core morphology. This suggests a strong sensitivity of the ICP flow field to heavy-impact vibrational excitations and dissociative events.
△ Less
Submitted 11 October, 2023; v1 submitted 3 October, 2023;
originally announced October 2023.
-
Multi-physics modeling of non-equilibrium phenomena in inductively coupled plasma discharges: Part I. A state-to-state approach
Authors:
Sanjeev Kumar,
Alessandro Munafo,
Sung Min Jo,
Marco Panesi
Abstract:
This work presents a vibrational and electronic state-to-state model for nitrogen plasma implemented within a multi-physics modular computational framework to study non-equilibrium effects in inductively coupled plasma (ICP) discharges. Within the computational framework, the set of vibronic (i.e., vibrational and electronic) master equations are solved in a tightly coupled fashion with the flow g…
▽ More
This work presents a vibrational and electronic state-to-state model for nitrogen plasma implemented within a multi-physics modular computational framework to study non-equilibrium effects in inductively coupled plasma (ICP) discharges. Within the computational framework, the set of vibronic (i.e., vibrational and electronic) master equations are solved in a tightly coupled fashion with the flow governing equations. This tight coupling eliminates the need for invoking any simplifying assumptions when computing the state of the plasma, thereby ensuring a higher degree of physical fidelity. To mitigate computational complexity, a maximum entropy coarse-graining strategy is deployed, effectively truncating the internal state space. The efficacy of this reduced StS model is empirically substantiated through zero-dimensional isochoric simulations. In these simulations, the results obtained from the reduced-order model are rigorously compared against those obtained from the full StS model, thereby confirming the accuracy of the reduced StS framework. The developed Coarse-grained StS model was employed to study the plasma discharge within the VKI Plasmatron facility. Our results reveal pronounced discrepancies between the plasma flow fields obtained from StS simulations and those derived from Local Thermodynamic Equilibrium (LTE) models, which are conventionally used in the simulation of such facilities. The analysis demonstrates a substantial departure of the internal state populations of atoms and molecules from the Boltzmann distribution. These nonequilibrium effects have important consequences on the energy coupling dynamics, thereby impacting the overall morphology of the plasma discharge. A deeper analysis of the results demonstrates that the population distribution is in a Quasi-Steady-State in the hot plasma core.
△ Less
Submitted 11 October, 2023; v1 submitted 3 October, 2023;
originally announced October 2023.
-
Rovibrational-Specific QCT and Master Equation Study on $\text{N}_2(\text{X}^1Σ_g^+)$+$\text{O}({}^3\text{P})$ and $\text{NO}(\text{X}^2Π)$+$\text{N}({}^4\text{S})$ Systems in High-Energy Collisions
Authors:
Sung Min Jo,
Simone Venturi,
Maitreyee P. Sharma,
Alessandro Munafò,
Marco Panesi
Abstract:
This work presents a detailed investigation of the energy transfer and dissociation mechanisms in $\text{N}_2(\text{X}^1Σ_g^+)$+$\text{O}({}^3\text{P})$ and $\text{NO}(\text{X}^2Π)$+$\text{N}({}^4\text{S})$ systems using rovibrational-specific quasi-classical trajectory (QCT) and master equation analyses. The complete set of state-to-state kinetic data, obtained via QCT, allows for an in-depth inv…
▽ More
This work presents a detailed investigation of the energy transfer and dissociation mechanisms in $\text{N}_2(\text{X}^1Σ_g^+)$+$\text{O}({}^3\text{P})$ and $\text{NO}(\text{X}^2Π)$+$\text{N}({}^4\text{S})$ systems using rovibrational-specific quasi-classical trajectory (QCT) and master equation analyses. The complete set of state-to-state kinetic data, obtained via QCT, allows for an in-depth investigation of the Zel'dovich mechanism leading to the formation of $\text{NO}$ molecules at microscopic and macroscopic scales. The master equation analysis demonstrates that the low-lying vibrational states of $\text{N}_2$ and $\text{NO}$ have dominant contributions to the $\text{NO}$ formation and the corresponding extinction of $\text{N}_2$ through the exchange process. For the considered temperature range, it is found that while nearly 50% of the dissociation processes for $\text{N}_2$ and $\text{NO}$ occurs in the molecular quasi-steady-state (QSS) regime, the amount of the Zel'dovich reaction is zero. Using the QSS approximation to model the Zel'dovich mechanism leads to an overestimation of $\text{NO}$ production by more than a factor of 4 in the high-temperature range. The breakdown of this well-known approximation has profound consequences for the approaches that heavily rely on the validity of QSS assumption in hypersonic applications. The investigation of the rovibrational state population dynamics reveals substantial similarity among different chemical systems for the energy transfer and the dissociation processes, providing promising physical foundations for the use of reduced-order strategies to other chemical systems without significant loss of accuracy.
△ Less
Submitted 2 November, 2021;
originally announced November 2021.
-
Non-equilibrium plasma generation via nano-second multi-mode laser pulses
Authors:
Andrea Alberti,
Alessandro Munafó,
Munetake Nishihara,
Carlos Pantano,
Jonathan B. Freund,
Marco Panesi
Abstract:
A non-equilibrium model for laser-induced plasmas is used to describe how nano-second temporal mode-beating affects plasma kernel formation and growth in quiescent air. The chemically reactive Navier-Stokes equations describe the hydrodynamics, and non-equilibrium effects are modeled based on a two-temperature model. Inverse Bremsstrahlung and multiphoton ionization are self-consistently taken int…
▽ More
A non-equilibrium model for laser-induced plasmas is used to describe how nano-second temporal mode-beating affects plasma kernel formation and growth in quiescent air. The chemically reactive Navier-Stokes equations describe the hydrodynamics, and non-equilibrium effects are modeled based on a two-temperature model. Inverse Bremsstrahlung and multiphoton ionization are self-consistently taken into account via a coupled solution of the equations governing plasma dynamics and beam propagation and attenuation (i.e., Radiative Transfer Equation). This strategy, despite the additional challenges it may bring, allows to minimize empiricism and enables for more accurate simulations since it does not require an artificial plasma seed to trigger breakdown. The benefits of this methodology are demonstrated by the good agreement between the predicted and the experimental plasma boundary evolution and absorbed energy. The same goes for the periodic plasma kernel structures which, as suggested by experiments and confirmed by the simulations discussed here, are linked to the modulating frequency.
△ Less
Submitted 22 December, 2021; v1 submitted 16 June, 2021;
originally announced June 2021.
-
A reduced-order NLTE kinetic model for radiating plasmas of outer envelopes of stellar atmospheres
Authors:
A. Munafo,
N. Mansour,
M. Panesi
Abstract:
The present work proposes a self-consistent reduced-order NLTE kinetic model for radiating plasmas such as are found in the outer layers of stellar atmospheres. Starting from the most up-to-date set of ab-initio and experimental data, the highly complex collisional-radiative kinetic mechanism is simplified by lumping the bound energy states in groups. Different grouping strategies are investigated…
▽ More
The present work proposes a self-consistent reduced-order NLTE kinetic model for radiating plasmas such as are found in the outer layers of stellar atmospheres. Starting from the most up-to-date set of ab-initio and experimental data, the highly complex collisional-radiative kinetic mechanism is simplified by lumping the bound energy states in groups. Different grouping strategies are investigated, such as uniform and Maxwell-Boltzmann. The reduced set of governing equations for the material gas and the radiation field is obtained based on a moment method. Applications consider the steady flow across a shock wave in partially ionized hydrogen. The results show that adopting a Maxwell-Boltzmann grouping allows, on the one hand, for a substantial reduction of the number of unknowns and, on the other, to maintain accuracy for both gas and radiation quantities. It is observed that, when neglecting line radiation, the use of two groups already leads to a very accurate resolution of the photo-ionization precursor, internal relaxation and radiative cooling regions. The inclusion of line radiation requires adopting just one additional group to account for optically thin losses in the alpha, beta and gamma lines of the Balmer and Paschen series. This trend has been observed for a wide range of shock wave velocities.
△ Less
Submitted 13 December, 2016;
originally announced December 2016.
-
A Spectral-Lagrangian Boltzmann Solver for a Multi-Energy Level Gas
Authors:
Alessandro Munafo,
Jeffrey R. Haack,
Irene M. Gamba,
Thierry E. Magin
Abstract:
In this paper a spectral-Lagrangian method for the Boltzmann equation for a multi-energy level gas is proposed. Internal energy levels are treated as separate species and inelastic collisions (leading to internal energy excitation and relaxation) are accounted for. The formulation developed can also be used for the case of a mixture of monatomic gases without internal energy (where only elastic co…
▽ More
In this paper a spectral-Lagrangian method for the Boltzmann equation for a multi-energy level gas is proposed. Internal energy levels are treated as separate species and inelastic collisions (leading to internal energy excitation and relaxation) are accounted for. The formulation developed can also be used for the case of a mixture of monatomic gases without internal energy (where only elastic collisions occur). The advantage of the spectral-Lagrangian method lies in the generality of the algorithm in use for the evaluation of the elastic and inelastic collision operators. The computational procedure is based on the Fourier transform of the partial elastic and inelastic collision operators and exploits the fact that these can be written as weighted convolutions in Fourier space with no restriction on the cross- section model. The conservation of mass, momentum and energy during collisions is enforced through the solution of constrained optimization problems. Numerical solutions are obtained for both space homogeneous and space in- homogeneous problems. Computational results are compared with those obtained by means of the DSMC method in order to assess the accuracy of the proposed spectral-Lagrangian method.
△ Less
Submitted 22 May, 2013;
originally announced May 2013.
