Comptes Rendus
Theoretical and numerical approaches for Vlasov–Maxwell equations
Coupled Particle-In-Cell and Direct Simulation Monte Carlo method for simulating reactive plasma flows
Comptes Rendus. Mécanique, Volume 342 (2014) no. 10-11, pp. 662-670.

Plasma flows with high Knudsen numbers cannot be treated with classic continuum methods, as represented for example by the Navier–Stokes or the magnetohydrodynamic equations. Instead, the more fundamental Boltzmann equation has to be solved, which is done here approximately by particle based methods that also allow for thermal and chemical non-equilibrium. The Particle-In-Cell method is used to treat the collisionless Vlasov–Maxwell system, while neutral reactive flows are treated by the Direct Simulation Monte Carlo method. In this article, a combined approach is presented that allows the simulation of reactive, partially or fully ionized plasma flows. Both particle methods are briefly outlined and the coupling and parallelization strategies are described. As an example, the results of a streamer discharge simulation are presented and discussed in order to demonstrate the capabilities of the coupled method.

Published online:
DOI: 10.1016/j.crme.2014.07.005
Keywords: Particle-In-Cell, Direct Simulation Monte Carlo, Discontinuous Galerkin, High-order, Plasma physics, Boltzmann equation

Claus-Dieter Munz 1; Monika Auweter-Kurtz 2; Stefanos Fasoulas 3; Asim Mirza 3; Philip Ortwein 1; Marcel Pfeiffer 3; Torsten Stindl 3

1 Institute of Aerodynamics and Gasdynamics (IAG), Universität Stuttgart, 70550 Stuttgart, Germany
2 German Aerospace Academy (ASA), Forum 1 am Konrad-Zuse-Platz 1, 71034 Böblingen, Germany
3 Institute of Space Systems (IRS), Universität Stuttgart, 70550 Stuttgart, Germany
     author = {Claus-Dieter Munz and Monika Auweter-Kurtz and Stefanos Fasoulas and Asim Mirza and Philip Ortwein and Marcel Pfeiffer and Torsten Stindl},
     title = {Coupled {Particle-In-Cell} and {Direct} {Simulation} {Monte} {Carlo} method for simulating reactive plasma flows},
     journal = {Comptes Rendus. M\'ecanique},
     pages = {662--670},
     publisher = {Elsevier},
     volume = {342},
     number = {10-11},
     year = {2014},
     doi = {10.1016/j.crme.2014.07.005},
     language = {en},
AU  - Claus-Dieter Munz
AU  - Monika Auweter-Kurtz
AU  - Stefanos Fasoulas
AU  - Asim Mirza
AU  - Philip Ortwein
AU  - Marcel Pfeiffer
AU  - Torsten Stindl
TI  - Coupled Particle-In-Cell and Direct Simulation Monte Carlo method for simulating reactive plasma flows
JO  - Comptes Rendus. Mécanique
PY  - 2014
SP  - 662
EP  - 670
VL  - 342
IS  - 10-11
PB  - Elsevier
DO  - 10.1016/j.crme.2014.07.005
LA  - en
ID  - CRMECA_2014__342_10-11_662_0
ER  - 
%0 Journal Article
%A Claus-Dieter Munz
%A Monika Auweter-Kurtz
%A Stefanos Fasoulas
%A Asim Mirza
%A Philip Ortwein
%A Marcel Pfeiffer
%A Torsten Stindl
%T Coupled Particle-In-Cell and Direct Simulation Monte Carlo method for simulating reactive plasma flows
%J Comptes Rendus. Mécanique
%D 2014
%P 662-670
%V 342
%N 10-11
%I Elsevier
%R 10.1016/j.crme.2014.07.005
%G en
%F CRMECA_2014__342_10-11_662_0
Claus-Dieter Munz; Monika Auweter-Kurtz; Stefanos Fasoulas; Asim Mirza; Philip Ortwein; Marcel Pfeiffer; Torsten Stindl. Coupled Particle-In-Cell and Direct Simulation Monte Carlo method for simulating reactive plasma flows. Comptes Rendus. Mécanique, Volume 342 (2014) no. 10-11, pp. 662-670. doi : 10.1016/j.crme.2014.07.005.

[1] G. Bird Molecular Gas Dynamics and the Direct Simulation of Gas Flows, Oxford University Press, Oxford, 1994

[2] R. Jahn Physics of Electric Propulsion, McGraw-Hill, New York, 1968

[3] O. Chanrion; T. Neubert A PIC–MCC code for the simulation of streamer propagation in air, J. Comput. Phys., Volume 126 (2008), pp. 434-448

[4] S. Pancheshnyi; P. Ségur; J. Capeillère; A. Bourdon Numerical simulation of filamentary discharges with parallel adaptive mesh refinement, J. Comput. Phys., Volume 227 (2008), pp. 6574-6590

[5] R.W. Hockney; J.W. Eastwood Computer Simulation Using Particles, McGraw-Hill, Inc., New York, 1988

[6] C.K. Birdsall; A.B. Langdon Plasma Physics via Computer Simulation, Hilger, Bristol, 1991

[7] T. Westermann Numerical modelling of the stationary Maxwell–Lorentz system in technical devices, Int. J. Numer. Model., Volume 7 (1994), pp. 43-67

[8] F. Hermeline; S. Layouni; P. Omnes A finite volume method for the approximation of Maxwell's equations in two space dimensions on arbitrary meshes, J. Comput. Phys., Volume 227 (2008), pp. 9365-9388

[9] F. Assous; P. Degond; E. Heintze; P. Raviart; J. Segre On a finite-element method for solving the three-dimensional Maxwell equations, J. Comput. Phys., Volume 109 (1993), pp. 222-237

[10] G.B. Jacobs; J. Hesthaven High-order nodal discontinuous Galerkin particle-in-cell method on unstructured grids, J. Comput. Phys., Volume 214 (2006), pp. 96-121

[11] F. Bergemann; I. Boyd New discrete vibrational energy model for the direct simulation Monte Carlo method, Progress in Astronautics and Aeronautics, vol. 158, 1994

[12] D. Liechty; M. Lewis Treatment of electronic energy level transition and ionization following the particle-based chemistry model, AIAA Paper, 8th Aerospace Sciences Meeting, 2010, 2010, pp. 1-12

[13] G. Bird, The DSMC method, CreateSpace Independent Publishing Platform, 2013.

[14] B.L. Haas; I.D. Boyd Models for direct Monte Carlo simulation of coupled vibration–dissociation, Phys. Fluids A, Fluid Dyn., Volume 5 (1993)

[15] G.A. Bird The Q–K model for gas-phase chemical reaction rates, Phys. Fluids, Volume 23 (2011), p. 106101

[16] A. Stock; J. Neudorfer; M. Riedlinger; G. Pirrung; G. Gassner; R. Schneider; S. Roller; C.-D. Munz Three-dimensional numerical simulation of a 30 GHz gyrotron resonator with an explicit high-order discontinuous Galerkin based parallel particle-on-cell method, IEEE Trans. Plasma Sci., Volume 40 (2012), pp. 1860-1870

[17] J. Neudorfer; A. Stock; J. Flamm; F. Hindenlang; G. Gassner; C.-D. Munz; R. Schneider; S. Roller Numerical investigation of high-order gyrotron mode propagation in launchers at 170 GHz, IEEE Trans. Plasma Sci., Volume 40 (2012), pp. 1512-1521

[18] A. Stock A high-order particle-in-cell method for low density plasma flow and the simulation of gyrotron resonator devices, University of Stuttgart, Germany, 2012 (Ph.D. thesis)

[19] C.-D. Munz; P. Omnes; R. Schneider; E. Sonnendrücker; U. Voß Divergence correction techniques for Maxwell solvers based on a hyperbolic model, J. Comput. Phys., Volume 161 (2000), pp. 484-511

[20] J. Hesthaven; T. Warburton Nodal Discontinuous Galerkin Methods: Algorithms, Analysis and Applications, Springer-Verlag, New York, 2008

[21] D.A. Kopriva Implementing Spectral Methods for Partial Differential Equations, Springer, 2009

[22] F. Hindenlang; G.J. Gassner; C. Altmann; A. Beck; M. Staudenmaier; C.-D. Munz Explicit discontinuous Galerkin methods for unsteady problems, Comput. Fluids, Volume 61 (2012), pp. 69-93

[23] T. Stindl; J. Neudorfer; A. Stock; M. Auweter-Kurtz; C.-D. Munz; S. Roller; R. Schneider Comparison of coupling techniques in a high-order discontinuous Galerkin-based particle-in-cell solver, J. Phys. D, Appl. Phys., Volume 44 (2011), p. 194004

[24] D. Baganoff; J.D. McDonald A collision selection rule for a particle simulation method suited to vector computers, Phys. Fluids A, Fluid Dyn., Volume 2 (1990), pp. 1248-1259

[25] S.E. Olson; A.J. Christlieb Gridless DSMC, J. Comput. Phys., Volume 227 (2008), pp. 8035-8064

[26] M. Pfeiffer; A. Mirza; S. Fasoulas A grid-independent particle pairing strategy for DSMC, J. Comput. Phys., Volume 246 (2013), pp. 28-36

[27] M. Laux Direkte Simulation verduennter, reagierender Stroemungen, Universitaet Stuttgart, 1995 (Ph.D. thesis)

[28] I.D. Boyd Modeling backward chemical rate processes in the direct simulation Monte Carlo method, Phys. Fluids, Volume 19 (2007), p. 126103

[29] M.H. Carpenter; C.A. Kennedy Fourth-order 2N-storage Kutta schemes Runge-schemes, NASA Technical Memorandum, 109112, 1994, pp. 1-26

[30] J. Neudorfer; A. Stock; R. Schneider; S. Roller; C.-D. Munz Efficient parallelization of a three-dimensional high-order particle-in-cell method for the simulation of a 170 GHz gyrotron resonator, IEEE Trans. Plasma Sci., Volume 41 (2013), pp. 87-98

[31] M. Pfeiffer; D. Petkow; G. Herdrich; S. Fasoulas Assessment of a numerical approach suitable for the M2P2 problem, Open Plasma Phys. J., Volume 4 (2011), pp. 24-33

Cited by Sources:

Comments - Policy