Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNoiseless Vlasov-Poisson simulations with linearly transformed particles
126
0
0.0
(
0
)
Added by
Martin Campos Pinto
Publication date
2012
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We introduce a deterministic discrete-particle simulation approach, the Linearly-Transformed Particle-In-Cell (LTPIC) method, that employs linear deformations of the particles to reduce the noise traditionally associated with particle schemes. Formally, transforming the particles is justified by local first order expansions of the characteristic flow in phase space. In practice the method amounts to using deformation matrices within the particle shape functions; these matrices are updated via local evaluations of the forward numerical flow. Because it is necessary to periodically remap the particles on a regular grid to avoid excessively deforming their shapes, the method can be seen as a development of Denavits Forward Semi-Lagrangian (FSL) scheme [J. Denavit, J. Comp. Physics 9, 75 (1972)]. However, it has recently been established [M. Campos Pinto, Smooth particle methods without smoothing, arXiv:1112.1859 (2012)] that the underlying Linearly-Transformed Particle scheme converges for abstract transport problems, with no need to remap the particles; deforming the particles can thus be seen as a way to significantly lower the remapping frequency needed in the FSL schemes, and hence the associated numerical diffusion. To couple the method with electrostatic field solvers, two specific charge deposition schemes are examined, and their performance compared with that of the standard deposition method. Finally, numerical 1d1v simulations involving benchmark test cases and halo formation in an initially mismatched thermal sheet beam demonstrate some advantages of our LTPIC scheme over the classical PIC and FSL methods. Benchmarked test cases also indicate that, for numerical choices involving similar computational effort, the LTPIC method is capable of accuracy comparable to or exceeding that of state-of-the-art, high-resolution Vlasov schemes.
rate research
Read More
In the paper, gridless particle techniques are presented in order to solve problems involving electrostatic, collisionless plasmas. The method makes use of computational particles having the shape of spherical shells or of rings, and can be used to study cases in which the plasma has spherical or axial symmetry, respectively. As a computational grid is absent, the technique is particularly suitable when the plasma occupies a rapidly changing space region.
Many phenomena in collisionless plasma physics require a kinetic description. The evolution of the phase space density can be modeled by means of the Vlasov equation, which has to be solved numerically in most of the relevant cases. One of the problems that often arise in such simulations is the violation of important physical conservation laws. Numerical diffusion in phase space translates into unphysical heating, which can increase the overall energy significantly, depending on the time scale and the plasma regime. In this paper, a general and straightforward way of improving conservation properties of Vlasov schemes is presented that can potentially be applied to a variety of different codes. The basic idea is to use fluid models with good conservation properties for correcting kinetic models. The higher moments that are missing in the fluid models are provided by the kinetic codes, so that both kinetic and fluid codes compensate the weaknesses of each other in a closed feedback loop.
Comparisons are presented between a hybrid Vlasov-Maxwell (HVM) simulation of turbulence in a collisionless plasma and fluid reductions. These include Hall-magnetohydrodynamics (HMHD) and Landau fluid (LF) or FLR-Landau fluid (FLR-LF) models that retain pressure anisotropy and low-frequency kinetic effects such as Landau damping and, for the last model, finite Larmor radius (FLR) corrections. The problem is considered in two space dimensions, when initial conditions involve moderate-amplitude perturbations of a homogeneous equilibrium plasma subject to an out-of-plane magnetic field. LF turns out to provide an accurate description of the velocity field up to the ion Larmor radius scale, and even to smaller scales for the magnetic field. Compressibility nevertheless appears significantly larger at the sub-ion scales in the fluid models than in the HVM simulation. High frequency kinetic effects, such as cyclotron resonances, not retained by fluid descriptions, could be at the origin of this discrepancy. A significant temperature anisotropy is generated, with a bias towards the perpendicular component, the more intense fluctuations being rather spread out and located in a broad vicinity of current sheets. Non-gyrotropic pressure tensor components are measured and their fluctuations are shown to reach a significant fraction of the total pressure fluctuation, with intense regions closely correlated with current sheets.
We present a way to combine Vlasov and two-fluid codes for the simulation of a collisionless plasma in large domains while keeping full information of the velocity distribution in localized areas of interest. This is made possible by solving the full Vlasov equation in one region while the remaining area is treated by a 5-moment two-fluid code. In such a treatment, the main challenge of coupling kinetic and fluid descriptions is the interchange of physically correct boundary conditions between the different plasma models. In contrast to other treatments, we do not rely on any specific form of the distribution function, e.g. a Maxwellian type. Instead, we combine an extrapolation of the distribution function and a correction of the moments based on the fluid data. Thus, throughout the simulation both codes provide the necessary boundary conditions for each other. A speed-up factor of around 20 is achieved by using GPUs for the computationally expensive solution of the Vlasov equation and an overall factor of at least 60 using the coupling strategy combined with the GPU computation. The coupled codes were then tested on the GEM reconnection challenge.
Electrostatic structures have been observed in many regions of space plasmas, including the solar wind, the magnetosphere, the auroral acceleration region. One possible theoretical description of some of these structures is the concept of Bernstein-Greene-Kruskal (BGK) modes, which are exact nonlinear steady-state solutions of the Vlasov-Poisson system of equations in collisionless kinetic theory. We generalize exact solutions of two-dimensional BGK modes in a magnetized plasma with finite magnetic field strength [Ng, Bhattacharjee, and Skiff, Phys. Plasmas {bf13}, 055903 (2006)] to cases with azimuthal magnetic fields so that these structures carry electric current as well as steady electric and magnetic fields. Such nonlinear solutions now satisfy exactly the Vlasov-Poisson-Amp`{e}re system of equations. Explicit examples with either positive or negative electric potential structure are provided.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
