Multi-scale two-domain numerical modeling of stationary positive DC corona discharge/drift-region coupling

被引:7
作者
Monrolin, Nicolas [2 ]
Plouraboue, Franck [1 ]
机构
[1] Toulouse Univ, Inst Fluid Mech Toulouse IMFT, CNRS, INPT,UPS, Toulouse, France
[2] Ecole Natl Aviat Civile, Toulouse, France
关键词
Corona discharge; Multi-scale modeling; Multi-domain coupling; Lagrange multipliers; Asymptotic analysis; Ionic wind; FINITE-ELEMENT-METHOD; SIMULATION; AIR; ALGORITHM;
D O I
10.1016/j.jcp.2021.110517
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
YY Corona discharge modelingmostly relies on two, mostly distinct, approaches: high-fidelity, numerically challenging, unsteady simulations having high-computational cost or lowfidelity simulations based on empirical assumptions such as constant electric field at the emitter electrode. For the purpose of steady discharge current predictions, high-fidelity models are very costly to use whilst empirical models have limited range of validity owing the subtle use of tuned parameters. We propose an intermediate approach: an asymptotic multi-scale/two-domain numerical modeling based upon generalizing previous asymptotic axi-symmetrical analysis [1,2]. We show how the initial elliptic (electric potential), hyperbolic (charge transport), non-local (photo-ionization) problem can be formulated into two local problems coupled by matching conditions. The approach relies on a multipole expansion of the radiative photo-ionization source term (in two dimensions for cylindrical emitters). The analytical asymptotic matching conditions derived in [2] result in flux continuity conditions at the boundary of the two domains. These coupling conditions are enforced by Lagrange multipliers, within a variational formulation, leading to a hierarchy of non-linear coupled problems. The proposed approach is both monolithic and two-domains: two asymptotic regions, an inner-one associated with corona discharge, and an outer-one, the ion drift region. Numerical convergence and validations of the finite element implementation is provided. A comparison with various experimental results convincingly demonstrate the applicability of the method, which avoids tuning parameters dedicated to each specific configuration, but, on the contrary, exclusively relies on known and measurable physical quantities (e.g., ion mobilities, photo-ionization coefficient, ionization electric field, Townsend discharge coefficient, etc...). (C) 2021 Elsevier Inc. All rights reserved.
引用
收藏
页数:28
相关论文
共 54 条
[1]  
Abramowitz I. A. M., 1965, HDB MATH FUNCTIONS
[2]   Simulation of corona discharge in point-plane configuration [J].
Adamiak, K ;
Atten, P .
JOURNAL OF ELECTROSTATICS, 2004, 61 (02) :85-98
[3]  
[Anonymous], 2016, BOLSIG SOLVER VER 03
[4]  
[Anonymous], 2020, PHELPS DATABASE N2 O
[5]  
[Anonymous], 1915, ELECT GASES
[6]   FINITE-ELEMENT METHOD WITH LAGRANGIAN MULTIPLIERS [J].
BABUSKA, I .
NUMERISCHE MATHEMATIK, 1973, 20 (03) :179-192
[7]   Multiple solutions in the theory of dc glow discharges and cathodic part of arc discharges. Application of these solutions to the modeling of cathode spots and patterns: a review [J].
Benilov, M. S. .
PLASMA SOURCES SCIENCE & TECHNOLOGY, 2014, 23 (05)
[8]   Modelling cathode spots in glow discharges in the cathode boundary layer geometry [J].
Bieniek, M. S. ;
Almeida, P. G. C. ;
Benilov, M. S. .
JOURNAL OF PHYSICS D-APPLIED PHYSICS, 2016, 49 (10)
[9]   Efficient models for photoionization produced by non-thermal gas discharges in air based on radiative transfer and the Helmholtz equations [J].
Bourdon, A. ;
Pasko, V. P. ;
Liu, N. Y. ;
Celestin, S. ;
Segur, P. ;
Marode, E. .
PLASMA SOURCES SCIENCE & TECHNOLOGY, 2007, 16 (03) :656-678
[10]   Narrow-Flow-Channel-Driven EHD Gas Pump for an Advanced Thermal Management of Microelectronics [J].
Chang, Jen-Shih ;
Tsubone, Hiroaki ;
Harvel, Glenn D. ;
Urashima, Kuniko .
IEEE TRANSACTIONS ON INDUSTRY APPLICATIONS, 2010, 46 (03) :1151-1158