Dynamic implicit 3D adaptive mesh refinement for non-equilibrium radiation diffusion

JOURNAL OF COMPUTATIONAL PHYSICS

卷:262

Article

Philip, B.¹  Wang, Z.²  Berrill, M. A.¹  Birke, M.²  Pernice, M.³ 

[1] Oak Ridge Natl Lab, Div Math & Comp Sci, Oak Ridge, TN 37831 USA

[2] Oak Ridge Natl Lab, Natl Ctr Computat Sci, Oak Ridge, TN 37831 USA

[3] Idaho Natl Lab, Idaho Falls, ID 83415 USA

关键词: Adaptive mesh refinement;    Jacobian Free Newton-Krylov;    Implicit methods;    Non equilibrium radiation diffusion;    Multilevel solvers;    Timestep control;   

DOI  :  10.1016/j.jcp.2013.12.058

来源: Elsevier

【 摘 要 】

The time dependent non-equilibrium radiation diffusion equations are important for solving the transport of energy through radiation in optically thick regimes and find applications in several fields including astrophysics and inertial confinement fusion. The associated initial boundary value problems that are encountered often exhibit a wide range of scales in space and time and are extremely challenging to solve. To efficiently and accurately simulate these systems we describe our research on combining techniques that will also find use more broadly for long term time integration of nonlinear multi-physics systems: implicit time integration for efficient long term time integration of stiff multiphysics systems, local control theory based step size control to minimize the required global number of time steps while controlling accuracy, dynamic 3D adaptive mesh refinement (AMR) to minimize memory and computational costs, Jacobian Free NewtonKrylov methods on AMR grids for efficient nonlinear solution, and optimal multilevel preconditioner components that provide level independent solver convergence. (C) 2014 Elsevier Inc. All rights reserved.

【 授权许可】

Free