Solution of equilibrium radiation diffusion problems using implicit adaptive mesh refinement

doi:10.1137/040609069

Solution of equilibrium radiation diffusion problems using implicit adaptive mesh refinement

2006-01-01

25

PRE

AI

DOI：10.1137/040609069

Original

Original request for help

Save

Abstract

En 中文

Diffusion approximations to radiation transport feature a nonlinear conduction coefficient that leads to formation of a sharp front, or Marshak wave, under suitable initial and boundary conditions. The front can vary several orders of magnitude over a very short distance. Resolving the shape of the Marshak wave is essential, but using a global. ne mesh can be prohibitively expensive. In such circumstances it is natural to consider using adaptive mesh refinement (AMR) to place a. ne mesh only in the vicinity of the propagating front. In addition, to avoid any loss of accuracy due to linearization, implicit time integration should be used to solve the equilibrium radiation diffusion equation. Implicit time integration on AMR grids introduces a new challenge, as algorithmic complexity must be controlled to fully realize the performance benefits of AMR. A Newton-Krylov method together with a multigrid preconditioner addresses this latter issue on a uniform grid. A straightforward generalization is to use a multilevel preconditioner that is tuned to the structure of the AMR grid, such as the fast adaptive composite grid (FAC) method. We describe the resulting Newton-Krylov-FAC method and demonstrate its performance on simple equilibrium radiation diffusion problems.

Keywords:

adaptive mesh refinement

implicit time integration

Newton-Krylov methods

fast adaptive composite grid

multilevel preconditioning

radiation diffusion

Solution of equilibrium radiation diffusion problems using implicit adaptive mesh refinement

Abstract

Journal

Organization