A multigroup diffusion solver using pseudo transient continuation for a radiation-hydrodynamic code with patch-based AMR

DOI10.1016/J.JCP.2007.09.019MaRDI QIDQ2476878zbMATH OpenOpenAlexFDO

Authors A. I. Shestakov, Stella S. R. Offner

Publication date 12 March 2008

Published in Journal of Computational Physics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/0710.4509

zbMATH Keywords

multigroup radiation diffusionpseudo transient continuationradiation-hydrodynamic codes with adaptive mesh refinement (AMR)

Mathematics Subject Classification ID

Mesh generation, refinement, and adaptive methods for the numerical solution of initial value and initial-boundary value problems involving PDEs (65M50) Multigrid methods; domain decomposition for initial value and initial-boundary value problems involving PDEs (65M55) Magnetohydrodynamics and electrohydrodynamics (76W05) Waves and radiation in optics and electromagnetic theory (78A40)

Abstract: We present a scheme to solve the nonlinear multigroup radiation diffusion (MGD) equations. The method is incorporated into a massively parallel, multidimensional, Eulerian radiation-hydrodynamic code with adaptive mesh refinement (AMR). The patch-based AMR algorithm refines in both space and time creating a hierarchy of levels, coarsest to finest. The physics modules are time-advanced using operator splitting. On each level, separate level-solve packages advance the modules. Our multigroup level-solve adapts an implicit procedure which leads to a two-step iterative scheme that alternates between elliptic solves for each group with intra-cell group coupling. For robustness, we introduce pseudo transient continuation (PTC). We analyze the magnitude of the PTC parameter to ensure positivity of the resulting linear system, diagonal dominance and convergence of the two-step scheme. For AMR, a level defines a subdomain for refinement. For diffusive processes such as MGD, the refined level uses Dirichet boundary data at the coarse-fine interface and the data is derived from the coarse level solution. After advancing on the fine level, an additional procedure, the sync-solve (SS), is required in order to enforce conservation. The MGD SS reduces to an elliptic solve on a combined grid for a system of G equations, where G is the number of groups. We adapt the partial temperature scheme for the SS; hence, we reuse the infrastructure developed for scalar equations. Results are presented. (Abridged)

Recommendations

• Solution of Equilibrium Radiation Diffusion Problems Using Implicit Adaptive Mesh Refinement
• Radiation diffusion for multi-fluid Eulerian hydrodynamics with adaptive mesh refinement.
• A multigrid Newton-Krylov method for multimaterial equilibrium radiation diffusion
• Numerical simulations for radiation hydrodynamics. I: Diffusion limit
• Dynamic implicit 3D adaptive mesh refinement for non-equilibrium radiation diffusion

Cites work

• A conservative adaptive projection method for the variable density incompressible Navier-Stokes equations
• Adaptive mesh refinement for hyperbolic partial differential equations
• An inexact Newton method for fully coupled solution of the Navier-Stokes equations with heat and mass transport
• Convergence Analysis of Pseudo-Transient Continuation
• scientific article; zbMATH DE number 4062720 (Why is no real title available?)
• scientific article; zbMATH DE number 3215568 (Why is no real title available?)
• Physical constraints in numerical calculations of diffusion
• Radiation diffusion for multi-fluid Eulerian hydrodynamics with adaptive mesh refinement.
• Solution of the diffusion equation by finite elements in Lagrangian hydrodynamic codes
• Stability of Stationary Solutions of the Multifrequency Radiation Diffusion Equations
• Time-dependent simulations of point explosions with heat conduction

Cited in

(5)

• A pyramid scheme for three-dimensional diffusion equations on polyhedral meshes
• A monotone finite volume scheme for diffusion equations on general non-conforming meshes
• Compatible and energy conserving multi-material arbitrary Lagrangian Eulerian scheme for multi-group radiation hydrodynamics simulations
• Ray decomposition radiation transport for high performance computing
• An unstructured block-based adaptive mesh refinement approach for discontinuous Galerkin method

This page was built for publication: A multigroup diffusion solver using pseudo transient continuation for a radiation-hydrodynamic code with patch-based AMR

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2476878)