Discretization of the Joule heating term for plasma discharge fluid models in unstructured meshes
Introduction
Fluid (continuum) models have been used since the 1960s to simulate plasma discharge phenomena. Fluid models describe the transport of electrons, ions and neutral species using moments of the species Boltzmann equation coupled to electromagnetic field equations. These models provide spatial and temporal information on averaged properties such as species density and species temperature in discharges and the electromagnetic fields in the discharge. A wide variety of plasmas such as low-pressure glow discharges [1], [2], [3], radio-frequency plasma discharges [4], [5], [6], dielectric barrier discharges [7], microdischarges [8], [9], [10], [11] and streamers [12] have been simulated using the fluid modeling framework.
A typical fluid model solves the following governing equations. The species continuity equation determines the individual species number densities in the discharge as followswhere, k represents the species index, is the species drift-diffusion flux, and is the gas-phase species generation rate through plasma chemical reactions. The self-consistent electric potential is determined using the electrostatic Poisson’s equationwhere, is the potential, e is the unit electric charge, is the permittivity of free space, is the total number of gas species, and is the charge number of species k (e.g. −1 for electrons). Electron transport coefficients and the rate coefficients of electron impact reactions are a strong function of the electron energy distribution function which may be parameterized by a local electron temperature [8], [9], [10], [11]. In case the local fluid approximation is not used, fluid models include an electron energy equation to determine the electron energy density in the dischargewhere, is the thermal conductivity of electrons, is the Boltzmann’s constant, and are the particle mass of electron and dominant background gas species, respectively, is the electron momentum transfer collision frequency with the background gas, is the energy lost per electron (in eV units) in an inelastic collision event represented by a gas-phase reaction is the rate of progress of a reaction j (in units), and is the total number of gas-phase reactions. The electrostatic field is determined from the electrostatic potential as . Finally, the gas temperature can be determined through the gas energy equationwhere, is the sensible enthalpy of species is the thermal conductivity of species is the energy lost from the thermal pool in an inelastic collision event represented by a gas reaction j, and the index h indicates the summation over all the heavy species. The above gas energy equation uses the constant pressure assumption and neglects fluid mechanical viscous dissipation.
The Joule heating source term, i.e. the first term on the RHS of (3), (4), constitutes the main source term in the electron energy balance and often times for the gas energy equation. Since reaction rates due to electron impact are exponential functions of the electron temperature, small deviations in the electron energy density can result in significant changes in the predicted discharge characteristics. Most reaction rates also depend on the number density of the background gas, and therefore on the gas temperature. The accuracy of the method used to compute the Joule heating term in the energy balance equations is therefore critical. As shown in [13], the stability of the numerical scheme also depends on the technique used to compute the source term of the electron energy equation owing to temporal stiffness introduced by this term for several types of plasma discharge problems.
Complex two-dimensional or three-dimensional geometries are now commonly encountered in the modeling of plasma discharges. Local mesh refinement is often needed to capture the steep gradients in the number density and in the electric field profiles. To optimize mesh resolution, adaptive mesh techniques can be pursued [12]. For these reasons, the use of unstructured meshes is necessary, as it already is the case for computational fluid dynamics. Robust and accurate numerical methods to simulate gas discharges on unstructured meshes must therefore be developed. Several issues specific to the self-consistent simulations of plasma phenomena using a cell-centered finite-volume approach arise in the context of unstructured meshes, the treatment of the Joule heating term being one. Here we discuss a simple and accurate approach for the numerical treatment of this term in the context of a cell-centered finite-volume discretization of the plasma discharge governing equations on generalized unstructured meshes. The technique can be used on structured meshes as well, where it reduces to a simplified numerical stencil.
Section snippets
Numerical method
The governing equations (1), (2), (3), (4) can be cast in an integral form, by integrating over an arbitrary control volume V aswhere, is the dependent variable (the species number densities , the electrostatic potential , the electron energy densities , or the gas temperature ), are the fluxes, and are the source terms. In the cell-centered finite-volume scheme for a fixed non-moving mesh, the following spatial discretization is used to approximate (5)
Test case
We consider a one-dimensional direct-current glow discharge between two infinitely large planar electrodes. The inter-electrode distance is 1 cm and the pressure is 1 Torr. A constant voltage difference of 250 V is applied between the two electrodes. The pure argon discharge gas is considered with the argon plasma modeled by three species (ground state atoms Ar, ions and electrons ). The species continuity equation is solved for the ions and the electrons, while the background argon number
Conclusion
An effective approach to numerical discretization of the Joule heating source term for fluid models of plasma discharge phenomena is presented in this paper. The method, which we call the “Reconstruction-Flux” discretization, has been developed for a cell-centered finite-volume scheme of the gas discharge governing equations on unstructured meshes. In this method, the Joule heating term is computed by evaluating the dot product of the electric field and the species flux at the cell center. The
References (17)
- et al.
Two-component plasma model for two-dimensional glow discharges in magnetic field
J. Comput. Phys.
(2004) - et al.
A numerical method to simulate radio-frequency plasma discharges
J. Comput. Phys.
(2002) - et al.
Speeding up fluid models for gas discharges by implicit treatment of the electron energy source term
J. Comput. Phys.
(2000) - et al.
Transport properties of gaseous ions over a wide energy range
At. Data Nucl. Data Tables
(1976) A two-dimensional model of dc glow discharges
J. Appl. Phys.
(1988)- S. Mahadevan, L.L. Raja, Simulations of direct-current air surface plasma discharges in supersonic flow, AIAA Paper...
- et al.
Fluid modeling of electron heating in low-pressure, high-frequency capacitively coupled plasma discharges
J. Appl. Phys.
(2004) Numerical model of rf glow discharges
Phys. Rev. A
(1987)
Cited by (7)
SOMAFOAM: An OpenFOAM based solver for continuum simulations of low-temperature plasmas
2021, Computer Physics CommunicationsAn axisymmetric unstructured finite volume method applied to the numerical modeling of an atmospheric pressure gas discharge
2015, Journal of Computational PhysicsCitation Excerpt :In this work a hybrid FVM and FEM was used together with an FCT algorithm and the Gummel scheme. One of the most recent works on unstructured meshes was presented by Deconinck et al. [20,21] and Breden et al. [22]. A cell-centered FVM was developed for the study of plasma bullets and micro-hollow cathode discharges on hybrid meshes (combination of triangles and quadrangles).
Numerical simulation of DC glow discharges for shock wave modification
2016, 54th AIAA Aerospace Sciences MeetingSimulations of direct-current air glow discharge at pressures ∼1 Torr: Discharge model validation
2010, Journal of Applied Physics