Stability analysis of a multilevel quasidiffusion method for thermal radiative transfer problems

Abstract In this paper we analyze a multilevel quasidiffusion (QD) method for solving time-dependent multigroup nonlinear radiative transfer problems which describe interaction of photons with matter. The multilevel method is formulated by means of the high-order radiative transfer equation and a set of low-order moment equations. The fully implicit scheme is used to discretize equations in time. The stability analysis is applied to the method in semi-continuous and discretized forms. To perform Fourier analysis, the system of equations of the multilevel method is linearized about an equilibrium solution. The effects of discretization with respect to different independent variables are studied. The multilevel method is shown to be stable and fast converging. We also consider a version of the method in which time evolution in the radiative transfer equation is treated by means of the α-approximation. The Fleck–Cummings test problem is used to demonstrate performance of the multilevel QD method and study its iterative stability.

[1]  Dmitriy Y. Anistratov,et al.  Computational transport methodology based on decomposition of a problem domain into transport and diffusive subdomains , 2012, J. Comput. Phys..

[2]  Jim E. Morel,et al.  Application of nonlinear Krylov acceleration to radiative transfer problems , 2013 .

[3]  Y. Zel’dovich,et al.  Gas Dynamics. (Book Reviews: Physics of Shock Waves and High-Temperature Hydrodynamic Phenomena. Vol. 1) , 1970 .

[4]  Edward W. Larsen.  Transport Acceleration Methods as Two-Level Multigrid Algorithms , 1991 .

[5]  Ryan G. McClarren,et al.  Semi-implicit time integration for PN thermal radiative transfer , 2008, J. Comput. Phys..

[6]  Edward W. Larsen,et al.  A grey transport acceleration method for time-dependent radiative transfer pro lems , 1988 .

[7]  Frank Graziani.  Computational Methods in Transport , 2006 .

[8]  Michael L. Norman,et al.  Implicit Adaptive-Grid Radiation Hydrodynamics , 1985 .

[9]  Rick M. Rauenzahn,et al.  A Consistent, Moment-Based, Multiscale Solution Approach for Thermal Radiative Transfer Problems , 2012 .

[10]  V.Ya. Gol'din,et al.  A quasi-diffusion method of solving the kinetic equation , 1964 .

[11]  Edward W. Larsen,et al.  A synthetic acceleration scheme for radiative diffusion calculations , 1985 .

[12]  E. Aristova.  Simulation of Radiation Transport in a Channel Based on the Quasi-Diffusion Method , 2008 .

[13]  Jim E. Morel,et al.  Linear multifrequency-grey acceleration recast for preconditioned Krylov iterations , 2007, J. Comput. Phys..

[14]  Marvin L. Adams,et al.  A hybrid transport-diffusion method for 2D transport problems with diffusive subdomains , 2014, J. Comput. Phys..

[15]  James Paul Holloway,et al.  On solutions to the Pn equations for thermal radiative transfer , 2008, J. Comput. Phys..

[16]  R. P. Drake,et al.  High-Energy-Density Physics: Fundamentals, Inertial Fusion, and Experimental Astrophysics , 2006 .

[17]  J. A. Fleck,et al.  An implicit Monte Carlo scheme for calculating time and frequency dependent nonlinear radiation transport , 1971 .

[18]  B. N. Chetverushkin,et al.  Methods of solving one-dimensional problems of radiation gas dynamics☆ , 1972 .

[19]  Frank H. Shu,et al.  The physics of astrophysics. , 1992 .

[20]  Edward W. Larsen,et al.  Fast iterative methods for discrete-ordinates particle transport calculations , 2002 .