skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Noiseless Vlasov-Poisson simulations with linearly transformed particles

Journal Article · · Journal of Computational Physics
 [1];  [2];  [3];  [3];  [3]
  1. Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); CNRS, UMR 7598, Lab. Jacques-Louis Lions, Paris (France); UPMC Univ. Paris 06, UMR 7598, Lab. Jacques-Louis Lions, Paris (France)
  2. IRMA, UMR 7501, Univ. de Strasbourg & CNRS, Strasbourg Cedex (France); Project-team CALVI, Strasbourg Cedex (France)
  3. Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)

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 of 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 Denavit's Forward Semi-Lagrangian (FSL) scheme (Denavit, 1972 [8]). However, it has recently been established (Campos Pinto, 2012 [20]) 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. Lastly, 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.

Research Organization:
Lawrence Livermore National Lab., Livermore, CA (United States)
Sponsoring Organization:
USDOE
Grant/Contract Number:
AC52-07NA27344
OSTI ID:
1281673
Report Number(s):
LLNL-JRNL-696802
Journal Information:
Journal of Computational Physics, Vol. 275, Issue C; ISSN 0021-9991
Publisher:
ElsevierCopyright Statement
Country of Publication:
United States
Language:
English
Citation Metrics:
Cited by: 14 works
Citation information provided by
Web of Science

References (31)

Study of laser-plasma beat wave current drive with an Eulerian Vlasov code journal January 1992
Numerical study of plasma-wall transition using an Eulerian Vlasov code journal July 2004
Numerical study of plasma–wall transition in an oblique magnetic field journal March 2001
Simulation of heavy ion beams with a semi-Lagrangian Vlasov solver
  • Sonnendrucker, Eric; Barnard, John J.; Friedman, Alex
  • Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, Vol. 464, Issue 1-3 https://doi.org/10.1016/S0168-9002(01)00186-3
journal May 2001
Vlasov simulations of electron holes driven by particle distributions from PIC reconnection simulations with a guide field journal January 2008
Numerical simulation of plasmas with periodic smoothing in phase space journal February 1972
A forward-trajectory global semi-Lagrangian transport scheme journal September 2003
A forward semi-Lagrangian method for the numerical solution of the Vlasov equation journal October 2009
A Lagrangian Particle‐Wavelet Method journal January 2006
A Particle-in-cell Method with Adaptive Phase-space Remapping for Kinetic Plasmas journal January 2011
Wavelet-Based Poisson Solver for Use in Particle-in-Cell Simulations journal June 2005
Wavelet denoising for postprocessing of a 2D Particle - In - Cell code journal January 2007
Self-consistent kinetic simulation of plasmas journal March 2000
A high order cell-centered semi-Lagrangian scheme for multi-dimensional kinetic simulations of neutral gas flows journal April 2012
Nonlinear Study of Vlasov's Equation for a Special Class of Distribution Functions journal January 1967
Convergence of a Variable Blob Vortex Method for the Euler and Navier–Stokes Equations journal December 1990
Optimal Approximations of Transport Equations by Particle and Pseudoparticle Methods journal January 2000
Grid and Particle Hydrodynamics: journal August 1998
Fragmentation, merging, and internal dynamics for PIC simulation with finite size particles journal August 2003
A cloudy Vlasov solution journal May 2005
Sheet beam model for intense space charge: Application to Debye screening and the distribution of particle oscillation frequencies in a thermal equilibrium beam journal May 2011
A characterization of multivariate quasi-interpolation formulas and its applications journal December 1990
On the approximation power of convolution-based least squares versus interpolation journal July 1997
Cubic interpolated propagation scheme for solving the hyper-dimensional Vlasov—Poisson equation in phase space journal August 1999
Comparison of Eulerian Vlasov solvers journal February 2003
A conservative high order semi-Lagrangian WENO method for the Vlasov equation journal February 2010
A New Class of Nonlinear Finite-Volume Methods for Vlasov Simulation journal September 2010
Discontinuous Galerkin semi-Lagrangian method for Vlasov-Poisson journal October 2011
A positivity-preserving high-order semi-Lagrangian discontinuous Galerkin scheme for the Vlasov–Poisson equations journal July 2011
A discontinuous Galerkin method for the Vlasov–Poisson system journal February 2012
High-energy density physics experiments with intense heavy ion beams
  • Bieniosek, F. M.; Henestroza, E.; Leitner, M.
  • Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, Vol. 606, Issue 1-2 https://doi.org/10.1016/j.nima.2009.03.123
journal July 2009

Cited By (1)

Convergence of a linearly transformed particle method for aggregation equations preprint January 2015

Similar Records

Noiseless Vlasov–Poisson simulations with linearly transformed particles
Journal Article · Wed Oct 15 00:00:00 EDT 2014 · Journal of Computational Physics · OSTI ID:1281673

A 4th-Order Particle-in-Cell Method with Phase-Space Remapping for the Vlasov–Poisson Equation
Journal Article · Tue May 09 00:00:00 EDT 2017 · SIAM Journal on Scientific Computing · OSTI ID:1281673

Stability of nonlinear Vlasov-Poisson equilibria through spectral deformation and Fourier-Hermite expansion
Journal Article · Sun May 15 00:00:00 EDT 2011 · Physical Review. E, Statistical, Nonlinear, and Soft Matter Physics (Print) · OSTI ID:1281673