A FAMILY OF FLUX-LIMITED DIFFUSION THEORIES

We present a family of diffusion approximations to the equation of radiative transfer, parameterized functionally by a function chi-(omega, R). Here omega is the effective (treating emission as scattering) single scatter albedo, and R is roughly speaking the magnitude of the dimensionless spatial gradient of the energy density. For any member of this family, i.e., for any function chi, this diffusion theory is flux-limited, properly predicts a single asymptotic mode in a sourcefree homogeneous medium, and gives the correct weak gradient limit. If the choice chi = 1 is made, this family reduces to the flux-limited diffusion theory proposed earlier by Levermore and Pomraning. We suggest a function chi that depends only upon the single variable omega. Simple test problems indicate that this choice leads to improved accuracy over both classic and the earlier flux-limited diffusion theory. This choice also allows the radiative flux and energy density gradient to have independent directions, and this may result in increased accuracy in multidimensional problems.