An implicit energy-conservative 2D Fokker-Planck algorithm -- 2. Jacobian-free Newton-Krylov solver
Energy-conservative implicit integration schemes for the Fokker-Planck transport equation in multidimensional geometries require inverting a dense, non-symmetric matrix (Jacobian), which is very expensive to store and solve using standard solvers. However, these limitations can be overcome with Newton-Krylov iterative techniques, since they can be implemented Jacobian-free (the Jacobian matrix from Newton's algorithm is never formed nor stored to proceed with the iteration), and their convergence can be accelerated by preconditioning the original problem. In this document, the efficient numerical implementation of an implicit energy-conservative scheme for multidimensional Fokker-Planck problems using multigrid-preconditioned Krylov methods is discussed. Results show that multigrid preconditioning is very effective in speeding convergence and decreasing CPU requirements, particularly in fine meshes. The solver is demonstrated on grids up to 128 x 128 points in a 2D cylindrical velocity space ({upsilon}{sub r}, {upsilon}{sub p}) with implicit time steps of the order of the collisional time scale of the problem, {tau}. The method preserves particles exactly, and energy conservation is improved over alternative approaches, particularly in coarse meshes. Typical errors in the total energy over a time period of 10{tau} remain below a percent.
- Research Organization:
- Univ. of Illinois, Urbana, IL (US)
- OSTI ID:
- 20015681
- Journal Information:
- Journal of Computational Physics, Vol. 157, Issue 2; Other Information: PBD: 20 Jan 2000; ISSN 0021-9991
- Country of Publication:
- United States
- Language:
- English
Similar Records
On the Preconditioning of a Newton-Krylov Solver for a High-Order reconstructed Discontinuous Galerkin Discretization of All-Speed Compressible Flow with Phase Change for Application in Laser-Based Additive Manufacturing
Recovery Discontinuous Galerkin Jacobian-Free Newton-Krylov Method for All-Speed Flows