An asymptotic preserving unified gas kinetic scheme for gray radiative transfer equations
Abstract
The solutions of radiative transport equations can cover both optical thin and optical thick regimes due to the large variation of photon's meanfree path and its interaction with the material. In the small mean free path limit, the nonlinear timedependent radiative transfer equations can converge to an equilibrium diffusion equation due to the intensive interaction between radiation and material. In the optical thin limit, the photon free transport mechanism will emerge. In this paper, we are going to develop an accurate and robust asymptotic preserving unified gas kinetic scheme (APUGKS) for the gray radiative transfer equations, where the radiation transport equation is coupled with the material thermal energy equation. The current work is based on the UGKS framework for the rarefied gas dynamics [14], and is an extension of a recent work [12] from a onedimensional linear radiation transport equation to a nonlinear twodimensional gray radiative system. The newly developed scheme has the asymptotic preserving (AP) property in the optically thick regime in the capturing of diffusive solution without using a cell size being smaller than the photon's mean free path and time step being less than the photon collision time. Besides the diffusion limit, the scheme can capture themore »
 Authors:
 Institute of Applied Physics and Computational Mathematics, P.O. Box 8009, Beijing 100088 (China)
 Department of Mathematics, Hong Kong University of Science and Technology, Clear Water Bay, Kowloon (Hong Kong)
 Publication Date:
 OSTI Identifier:
 22465611
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Computational Physics; Journal Volume: 285; Other Information: Copyright (c) 2015 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; 97 MATHEMATICAL METHODS AND COMPUTING; ASYMPTOTIC SOLUTIONS; DIFFUSION EQUATIONS; EXACT SOLUTIONS; MEAN FREE PATH; NONLINEAR PROBLEMS; ONEDIMENSIONAL CALCULATIONS; PHOTON COLLISIONS; RADIANT HEAT TRANSFER; RADIATION TRANSPORT; TIME DEPENDENCE; TRANSPORT THEORY; TWODIMENSIONAL SYSTEMS
Citation Formats
Sun, Wenjun, Email: sun_wenjun@iapcm.ac.cn, Jiang, Song, Email: jiang@iapcm.ac.cn, and Xu, Kun, Email: makxu@ust.hk. An asymptotic preserving unified gas kinetic scheme for gray radiative transfer equations. United States: N. p., 2015.
Web. doi:10.1016/J.JCP.2015.01.008.
Sun, Wenjun, Email: sun_wenjun@iapcm.ac.cn, Jiang, Song, Email: jiang@iapcm.ac.cn, & Xu, Kun, Email: makxu@ust.hk. An asymptotic preserving unified gas kinetic scheme for gray radiative transfer equations. United States. doi:10.1016/J.JCP.2015.01.008.
Sun, Wenjun, Email: sun_wenjun@iapcm.ac.cn, Jiang, Song, Email: jiang@iapcm.ac.cn, and Xu, Kun, Email: makxu@ust.hk. 2015.
"An asymptotic preserving unified gas kinetic scheme for gray radiative transfer equations". United States.
doi:10.1016/J.JCP.2015.01.008.
@article{osti_22465611,
title = {An asymptotic preserving unified gas kinetic scheme for gray radiative transfer equations},
author = {Sun, Wenjun, Email: sun_wenjun@iapcm.ac.cn and Jiang, Song, Email: jiang@iapcm.ac.cn and Xu, Kun, Email: makxu@ust.hk},
abstractNote = {The solutions of radiative transport equations can cover both optical thin and optical thick regimes due to the large variation of photon's meanfree path and its interaction with the material. In the small mean free path limit, the nonlinear timedependent radiative transfer equations can converge to an equilibrium diffusion equation due to the intensive interaction between radiation and material. In the optical thin limit, the photon free transport mechanism will emerge. In this paper, we are going to develop an accurate and robust asymptotic preserving unified gas kinetic scheme (APUGKS) for the gray radiative transfer equations, where the radiation transport equation is coupled with the material thermal energy equation. The current work is based on the UGKS framework for the rarefied gas dynamics [14], and is an extension of a recent work [12] from a onedimensional linear radiation transport equation to a nonlinear twodimensional gray radiative system. The newly developed scheme has the asymptotic preserving (AP) property in the optically thick regime in the capturing of diffusive solution without using a cell size being smaller than the photon's mean free path and time step being less than the photon collision time. Besides the diffusion limit, the scheme can capture the exact solution in the optical thin regime as well. The current scheme is a finite volume method. Due to the direct modeling for the time evolution solution of the interface radiative intensity, a smooth transition of the transport physics from optical thin to optical thick can be accurately recovered. Many numerical examples are included to validate the current approach.},
doi = {10.1016/J.JCP.2015.01.008},
journal = {Journal of Computational Physics},
number = ,
volume = 285,
place = {United States},
year = 2015,
month = 3
}

This paper presents an extension of previous work (Sun et al., 2015 [22]) of the unified gas kinetic scheme (UGKS) for the gray radiative transfer equations to the frequencydependent (multigroup) radiative transfer system. Different from the gray radiative transfer equations, where the optical opacity is only a function of local material temperature, the simulation of frequencydependent radiative transfer is associated with additional difficulties from the frequencydependent opacity. For the multiple frequency radiation, the opacity depends on both the spatial location and the frequency. For example, the opacity is typically a decreasing function of frequency. At the same spatial region themore »

On the asymptotic preserving property of the unified gas kinetic scheme for the diffusion limit of linear kinetic models
The unified gas kinetic scheme (UGKS) of K. Xu et al. (2010) [37], originally developed for multiscale gas dynamics problems, is applied in this paper to a linear kinetic model of radiative transfer theory. While such problems exhibit purely diffusive behavior in the optically thick (or small Knudsen) regime, we prove that UGKS is still asymptotic preserving (AP) in this regime, but for the free transport regime as well. Moreover, this scheme is modified to include a time implicit discretization of the limit diffusion equation, and to correctly capture the solution in case of boundary layers. Contrary to many APmore » 
An asymptoticpreserving stochastic Galerkin method for the radiative heat transfer equations with random inputs and diffusive scalings
In this paper, we develop an AsymptoticPreserving (AP) stochastic Galerkin scheme for the radiative heat transfer equations with random inputs and diffusive scalings. In this problem the random inputs arise due to uncertainties in cross section, initial data or boundary data. We use the generalized polynomial chaos based stochastic Galerkin (gPCSG) method, which is combined with the micro–macro decomposition based deterministic AP framework in order to handle efficiently the diffusive regime. For linearized problem we prove the regularity of the solution in the random space and consequently the spectral accuracy of the gPCSG method. We also prove the uniform (inmore » 
A stochastic asymptoticpreserving scheme for a kineticfluid model for disperse twophase flows with uncertainty
In this paper we consider a kineticfluid model for disperse twophase flows with uncertainty. We propose a stochastic asymptoticpreserving (sAP) scheme in the generalized polynomial chaos stochastic Galerkin (gPCsG) framework, which allows the efficient computation of the problem in both kinetic and hydrodynamic regimes. The sAP property is proved by deriving the equilibrium of the gPC version of the Fokker–Planck operator. The coefficient matrices that arise in a Helmholtz equation and a Poisson equation, essential ingredients of the algorithms, are proved to be positive definite under reasonable and mild assumptions. The computation of the gPC version of a translation operatormore » 
A unified gaskinetic scheme for continuum and rarefied flows IV: Full Boltzmann and model equations
Fluid dynamic equations are valid in their respective modeling scales, such as the particle mean free path scale of the Boltzmann equation and the hydrodynamic scale of the Navier–Stokes (NS) equations. With a variation of the modeling scales, theoretically there should have a continuous spectrum of fluid dynamic equations. Even though the Boltzmann equation is claimed to be valid in all scales, many Boltzmann solvers, including direct simulation Monte Carlo method, require the cell resolution to the order of particle mean free path scale. Therefore, they are still single scale methods. In order to study multiscale flow evolution efficiently, themore »