Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNumerical Methods and Comparisons for 1D and Quasi 2D Fluid Streamer Propagation Models
76
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
In this work, we propose and compare four different strategies to simulate the fluid model for streamer propagation in one-dimension (1D) and quasi two-dimension (2D), which consists of a Poissons equation for particle velocity and two continuity equations for particle transport. Each strategy involves of one method for solving Poissons equation and the other for solving continuity equations, and a total variation diminishing three-stage Runge-Kutta method in temporal discretization. The numerical methods for Poissons equation include finite volume method, discontinuous Galerkin methods, mixed finite element method and least-squared finite element method. The numerical method for continuity equations is chosen from the family of discontinuous Galerkin methods. The accuracy tests and comparisons show that all of these four strategies are suitable and competitive in streamer simulations from the aspects of accuracy and efficiency. Results show these methods are compatible. By applying any strategy in real simulations, we can study the dynamics of streamer propagations in both 1D and quasi 2D models.
rate research
Read More
We have developed an efficient algorithm for steady axisymmetrical 2D fluid equations. The algorithm employs multigrid method as well as standard implicit discretization schemes for systems of partial differential equations. Linearity of the multigrid method with respect to the number of grid points allowed us to use $256times 256$ grid, where we could achieve solutions in several minutes. Time limitations due to nonlinearity of the system are partially avoided by using multi level grids(the initial solution on $256times 256$ grid was extrapolated steady solution from $128times 128$ grid which allowed using long integration time steps). The fluid solver may be used as the basis for hybrid codes for DC discharges.
Magnetic reconnection (MR) plays a fundamental role in plasma dynamics under many different conditions, from space and astrophysical environments to laboratory devices. High-resolution in-situ measurements from space missions allow to study naturally occurring MR processes in great detail. Alongside direct measurements, numerical simulations play a key role in investigating the fundamental physics underlying MR. The choice of an adequate plasma model to be employed in numerical simulations, while also compromising with their computational cost, is crucial to efficiently address the problem. We consider a new plasma model that includes a refined electron response within the hybrid-kinetic framework (kinetic ions, fluid electrons). The extent to which this new model can reproduce a full-kinetic description of 2D MR, with particular focus on its robustness during the non-linear stage, is evaluated. We perform 2D simulations of MR with moderate guide field by means of three different plasma models: a hybrid-Vlasov-Maxwell model with isotropic, isothermal electrons, a hybrid-Vlasov-Landau-fluid (HVLF) model where an anisotropic electron fluid is equipped with a Landau-fluid closure, and a full-kinetic one. When compared to the full-kinetic case, the HVLF model effectively reproduces the main features of MR, as well as several aspects of the associated electron micro-physics and its feedback onto proton dynamics. This includes the global evolution of MR and the local physics occurring within the so-called electron-diffusion region, as well as the evolution of species pressure anisotropy. In particular, anisotropy driven instabilities (such as firehose, mirror, and cyclotron instabilities) play a relevant role in regulating electrons anisotropy during the non-linear stage of MR. As expected, the HVLF model captures all these features, except for the electron-cyclotron instability.
The propagation mechanisms of plasma streamers have been observed and investigated in a surface dielectric barrier discharge (SDBD) using 2D particle in cell simulations. The investigations are carried out under a simulated air mixture, 80% N$_2$ and 20% O$_2$, at atmospheric pressure, 100$,$kPa, under both DC conditions and a pulsed DC waveform that represent AC conditions. The simulated geometry is a simplification of the symmetric and fully exposed SDBD resulting in the simultaneous ignition of both positive and negative streamers on either side of the Al$_2$O$_3$ dielectric barrier. In order to determine the interactivity of the two streamers, the propagation behavior for the positive and negative streamers are investigated both independently and simultaneously under identical constant voltage conditions. An additional focus is implored under a fast sub nanosecond rise time square voltage pulse alternating between positive and negative voltage conditions, thus providing insight into the dynamics of the streamers under alternating polarity switches. It is shown that the simultaneous ignition of both streamers, as well as using the pulsed DC conditions, provides both an enhanced discharge and an increased surface coverage. It is also shown that additional streamer branching may occur in a cross section that is difficult to experimentally observe. The enhanced discharge and surface coverage may be beneficial to many applications such as, but are not limited to: air purification, volatile organic compound removal, and plasma enhanced catalysis.
The document describes a numerical algorithm to simulate plasmas and fluids in the 3 dimensional space by the Euler method, in which the spatial meshes are fixed to the space. The plasmas and fluids move through the spacial Euler mesh boundary. The Euler method can represent a large deformation of the plasmas and fluids. On the other hand, when the plasmas or fluids are compressed to a high density, the spatial resolution should be ensured to describe the density change precisely. The present 3D Euler code is developed to simulate a nuclear fusion fuel ignition and burning. Therefore, the 3D Euler code includes the DT fuel reactions, the alpha particle diffusion, the alpha particle deposition to heat the DT fuel and the DT fuel depletion by the DT reactions, as well as the thermal energy diffusion based on the three-temperature compressible fluid model.
To help reveal the complete picture of linear kinetic drift modes, four independent numerical approaches, based on integral equation, Euler initial value simulation, Euler matrix eigenvalue solution and Lagrangian particle simulation, respectively, are used to solve the linear gyrokinetic electrostatic drift modes equation in Z-pinch with slab simplification and in tokamak with ballooning space coordinate. We identify that these approaches can yield the same solution with the difference smaller than 1%, and the discrepancies mainly come from the numerical convergence, which is the first detailed benchmark of four independent numerical approaches for gyrokinetic linear drift modes. Using these approaches, we find that the entropy mode and interchange mode are on the same branch in Z-pinch, and the entropy mode can have both electron and ion branches. And, at strong gradient, more than one eigenstate of the ion temperature gradient mode (ITG) can be unstable and the most unstable one can be on non-ground eigenstates. The propagation of ITGs from ion to electron diamagnetic direction at strong gradient is also observed, which implies that the propagation direction is not a decisive criterion for the experimental diagnosis of turbulent mode at the edge plasmas.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
