10 documents found, page 1 of 1
Acceleration of the particle-in-cell code Osiris with graphics processing units
Lee, R. P.; Pierce, J. R.; Miller, K. G.; Almanza, M.; Tableman, A.; Decyk, V. K.; Fonseca, R. A.; Alves E. P.; Mori, W. B.Fully relativistic particle-in-cell (PIC) simulations are crucial for advancing our knowledge of plasma physics. Modern supercomputers based on graphics processing units (GPUs) offer the potential to perform PIC simulations of unprecedented scale, but require robust and feature-rich codes that can fully leverage their computational resources. In this work, this demand is addressed by adding GPU acceleration to ...
A new field solver for modeling of relativistic particle-laser interactions usi...
Li, F.; Miller, K. G.; Xu, X.; Tsung, F. S.; Decyk, V. K.; An, W.; Fonseca, R. A.; Mori, W. B.A customized finite-difference field solver for the particle-in-cell (PIC) algorithm that provides higher fidelity for wave-particle interactions in intense electromagnetic waves is presented. In many problems of interest, particles with relativistic energies interact with intense electromagnetic fields that have phase velocities near the speed of light. Numerical errors can arise due to (1) dispersion errors i...
Controlling the numerical Cerenkov instability in PIC simulations using a custo...
Li, F.; Yu, P. C.; Xu, X.; Fiuza, F.; Decyk, V. K.; Dalichaouch, T.; Davidson, A.; Tableman, A.; An, W.; Tsung, F. S.; Fonseca, R. A.; Lu, W.In this paper we present a customized finite-difference-time-domain (FDTD) Maxwell solver for the particle-in-cell (PIC) algorithm. The solver is customized to effectively eliminate the numerical Cerenkov instability (NCI) which arises when a plasma (neutral or non-neutral) relativistically drifts on a grid when using the PIC algorithm. We control the EM dispersion curve in the direction of the plasma drift of ...
Enabling Lorentz boosted frame particle-in-cell simulations of laser wakefield ...
Yu, P.; Xu, X.; Davidson, A.; Tableman, A.; Dalichaouch, T.; Li, F.; Meyers, M. D.; An, W.; Tsung, F. S.; Decyk, V. K.; Fiuza, F.; Vieira, J.When modeling laser wakefield acceleration (LWFA) using the particle-in-cell (PIC) algorithm in a Lorentz boosted frame, the plasma is drifting relativistically at beta(b)c towards the laser, which can lead to a computational speedup of similar to gamma(2)(b)=(1-beta(2)(b))-1. Meanwhile, when LWFA is modeled in the quasi-3D geometry in which the electromagnetic fields and current are decomposed into a limited n...
Mitigation of numerical Cerenkov radiation and instability using a hybrid finit...
Yu, P.; Xu, X.; Tableman, A.; Decyk, V. K.; Tsung, F. S.; Fiuza, F.; Vieira, J.; Davidson, A.; Fonseca, R. A.; Lu, W.; Silva, L. O.; Mori, W. B.A hybrid Maxwell solver for fully relativistic and electromagnetic (EM) particle-in-cell (PIC) codes is described. In this solver, the EM fields are solved in k space by performing an FFT in one direction, while using finite difference operators in the other direction(s). This solver eliminates the numerical Cerenkov radiation for particles moving in the preferred direction. Moreover, the numerical Cerenkov ins...
Elimination of the numerical Cerenkov instability for spectral EM-PIC codes
Yu, P.; Xu, X.; Decyk, V. K.; Fiuza, F.; Vieira, J.; Tsung, F. S.; Fonseca, R. A.; Lu, W.; Silva, L. O.; Mori, W. B.When using an electromagnetic particle-in-cell (EM-PIC) code to simulate a relativistically drifting plasma, a violent numerical instability known as the numerical Cerenkov instability (NCI) occurs. The NCI is due to the unphysical coupling of electromagnetic waves on a grid to wave-particle resonances, including aliased resonances, i.e., omega + 2 pi mu/Delta t = (k(1) + 2 pi v(1)/Delta x(1))nu(0), where mu an...
Modeling of laser wakefield acceleration in Lorentz boosted frame using EM-PIC ...
Yu, P.; Xu, X.; Decyk, V. K.; An, W.; Vieira, J.; Tsung, F. S.; Fonseca, R. A.; Lu, W.; Silva, L. O.; Mori, W. B.Simulating laser wakefield acceleration (LWFA) in a Lorentz boosted frame in which the plasma drifts towards the laser with nu(b) can speed up the simulation by factors of gamma(2)(b) = (1 nu(2)(b)/c(2))(-1). In these simulations the relativistic drifting plasma inevitably induces a high frequency numerical instability that contaminates the interesting physics. Various approaches have been proposed to mitigate ...
Modeling of laser wakefield acceleration in Lorentz boosted frame using a Quasi...
Yu, P.; Davidson, A.; Xu, X.; Tableman, A.; Decyk, V. K.; Tsung, F. S.; Vieira, J.; Fonseca, R. A.; Lu, W.; Silva, L. O.; Mori, W. B.Recently it was proposed in [A. F. Lifschitz, et. al., J. Comp. Phys. 228, 1803 (2009)] that laser wakefield acceleration could be modeled efficiently using a particle-in-cell code in cylindrical coordinates if the fields and currents were expanded into Fourier modes in the azimuthal angle, ?. We have implemented this algorithm into OSIRIS, including a new rigorous charge conserving deposition routine applicabl...
Numerical instability due to relativistic plasma drift in EM-PIC simulations
Xu, X.; Yu, P.; Martins, S. F.; Tsung, F. S.; Decyk, V. K.; Vieira, J.; Fonseca, R. A.; Lu, W.; Silva, L. O.; Mori, W. B.The numerical instability observed in electromagnetic particle-in-cell (EM-PIC) simulations with a plasma drifting with relativistic velocities is studied using both theory and computer simulations. We derive the numerical dispersion relation for a cold plasma drifting with a relativistic velocity, and find an instability attributed to the intersection between beam resonances and the electromagnetic modes in th...
One-to-one full-scale simulations of laser-wakefield acceleration using QuickPIC
Vieira, J.; Fiuza, F.; Fonseca, R. A.; Silva, L. O.; Huang, C. K.; Lu, W.; Tzoufras, M.; Tsung, F. S.; Decyk, V. K.; Mori, W. B.; Cooley, J.We use the quasi-static particle-in-cell (PIC) code QuickPIC to perform full-scale one-to-one laser-waketield acceleration (LWFA) numerical experiments, with parameters that closely follow current experimental conditions. The propagation of state-of-the-art laser pulses in both preformed and uniform plasma channels is examined. We show that the presence of the channel is important whenever the laser self-modula...