Dynamic Implicit 3D Adaptive Mesh Refinement for Non-Equilibrium Radiation Diffusion

Journal Article · Wed Jan 01 00:00:00 EST 2014 · Journal of Computational Physics

DOI:https://doi.org/10.1016/j.jcp.2013.12.058· OSTI ID:1116470

Philip, Bobby ^[1]; Wang, Zhen ^[1]; Berrill, Mark A ^[1]; Rodriguez Rodriguez, Manuel ^[1]; Pernice, Michael ^[2]

---

1. ORNL
2. Idaho National Laboratory (INL)

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 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 multiphysics 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 Newton Krylov methods on AMR grids for efficient nonlinear solution, and optimal multilevel preconditioner components that provide level independent linear solver convergence.

• Research Organization: Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States). Oak Ridge Leadership Computing Facility (OLCF)
• Sponsoring Organization: USDOE Office of Science (SC)
• DOE Contract Number: DE-AC05-00OR22725
• OSTI ID: 1116470
• Journal Information: Journal of Computational Physics, Vol. 262; ISSN 0021-9991
• Country of Publication: United States
• Language: English