NAJul 25, 2011
A new numerical strategy with space-time adaptivity and error control for multi-scale streamer discharge simulationsMax Duarte, Zdenek Bonaventura, Marc Massot et al.
This paper presents a new resolution strategy for multi-scale streamer discharge simulations based on a second order time adaptive integration and space adaptive multiresolution. A classical fluid model is used to describe plasma discharges, considering drift-diffusion equations and the computation of electric field. The proposed numerical method provides a time-space accuracy control of the solution, and thus, an effective accurate resolution independent of the fastest physical time scale. An important improvement of the computational efficiency is achieved whenever the required time steps go beyond standard stability constraints associated with mesh size or source time scales for the resolution of the drift-diffusion equations, whereas the stability constraint related to the dielectric relaxation time scale is respected but with a second order precision. Numerical illustrations show that the strategy can be efficiently applied to simulate the propagation of highly nonlinear ionizing waves as streamer discharges, as well as highly multi-scale nanosecond repetitively pulsed discharges, describing consistently a broad spectrum of space and time scales as well as different physical scenarios for consecutive discharge/post-discharge phases, out of reach of standard non-adaptive methods.
APFeb 24, 2015
A numerical strategy to discretize and solve the Poisson equation on dynamically adapted multiresolution grids for time-dependent streamer discharge simulationsMax Duarte, Zdenek Bonaventura, Marc Massot et al.
We develop a numerical strategy to solve multi-dimensional Poisson equations on dynamically adapted grids for evolutionary problems disclosing propagating fronts. The method is an extension of the multiresolution finite volume scheme used to solve hyperbolic and parabolic time-dependent PDEs. Such an approach guarantees a numerical solution of the Poisson equation within a user-defined accuracy tolerance. Most adaptive meshing approaches in the literature solve elliptic PDEs level-wise and hence at uniform resolution throughout the set of adapted grids. Here we introduce a numerical procedure to represent the elliptic operators on the adapted grid, strongly coupling inter grid relations that guarantee the conservation and accuracy properties of multiresolution finite volume schemes. The discrete Poisson equation is solved at once over the entire computational domain as a completely separate process. The accuracy and numerical performance of the method is assessed in the context of streamer discharge simulations.