OATAO - Open Archive Toulouse Archive Ouverte Open Access Week

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

Monrolin, Nicolas and Plouraboué, Franck Multi-scale two-domain numerical modeling of stationary positive DC corona discharge/drift-region coupling. (2021) Journal of Computational Physics, 443. 110517. ISSN 0021-9991

[img] (Document in English)

PDF (Author's version) - Depositor and staff only until 17 December 2021 - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader
1MB

Official URL: https://doi.org/10.1016/j.jcp.2021.110517

Abstract

Corona discharge modeling mostly relies on two, mostly distinct, approaches: high-fidelity, numerically challenging, unsteady simulations having high-computational cost or low-fidelity 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...).

Item Type:Article
Audience (journal):International peer-reviewed journal
Uncontrolled Keywords:
Institution:French research institutions > Centre National de la Recherche Scientifique - CNRS (FRANCE)
Université de Toulouse > Ecole Nationale de l'Aviation Civile - ENAC (FRANCE)
Université de Toulouse > Institut National Polytechnique de Toulouse - Toulouse INP (FRANCE)
Université de Toulouse > Université Toulouse III - Paul Sabatier - UT3 (FRANCE)
Laboratory name:
Funders:
Centre National d'Études Spatiales - CNES (FRANCE) - Région Occitanie (FRANCE)
Statistics:download
Deposited On:19 Jul 2021 08:40

Repository Staff Only: item control page