A Fast Solver for Implicit Integration of the Vlasov--Poisson System in the Eulerian Framework
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
- Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
In this paper, we present a domain decomposition algorithm to accelerate the solution of Eulerian-type discretizations of the linear, steady-state Vlasov equation. The steady-state solver then forms a key component in the implementation of fully implicit or nearly fully implicit temporal integrators for the nonlinear Vlasov--Poisson system. The solver relies on a particular decomposition of phase space that enables the use of sweeping techniques commonly used in radiation transport applications. The original linear system for the phase space unknowns is then replaced by a smaller linear system involving only unknowns on the boundary between subdomains, which can then be solved efficiently with Krylov methods such as GMRES. Steady-state solves are combined to form an implicit Runge--Kutta time integrator, and the Vlasov equation is coupled self-consistently to the Poisson equation via a linearized procedure or a nonlinear fixed-point method for the electric field. Finally, numerical results for standard test problems demonstrate the efficiency of the domain decomposition approach when compared to the direct application of an iterative solver to the original linear system.
- Research Organization:
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
- Sponsoring Organization:
- USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR); ORNL Laboratory Directed Research and Development (LDRD) Program
- Grant/Contract Number:
- AC52-06NA25396; AC05-00OR22725
- OSTI ID:
- 1435515
- Report Number(s):
- LA-UR-17-24638
- Journal Information:
- SIAM Journal on Scientific Computing, Vol. 40, Issue 2; ISSN 1064-8275
- Publisher:
- SIAMCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Web of Science
Similar Records
Recovery Discontinuous Galerkin Jacobian-Free Newton-Krylov Method for All-Speed Flows
Terascale Optimal PDE Simulations (TOPS) Center