skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Linear-Multi-Frequency-Grey Preconditioning for Radiative Transfer SN Calculations

  1. Los Alamos National Laboratory
Publication Date:
Research Org.:
Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Sponsoring Org.:
OSTI Identifier:
Report Number(s):
DOE Contract Number:
Resource Type:
Technical Report
Country of Publication:
United States
Mathematics & Computing(97); Nuclear Physics & Radiation Physics(73)

Citation Formats

Till, Andrew T. Linear-Multi-Frequency-Grey Preconditioning for Radiative Transfer SN Calculations. United States: N. p., 2013. Web. doi:10.2172/1090639.
Till, Andrew T. Linear-Multi-Frequency-Grey Preconditioning for Radiative Transfer SN Calculations. United States. doi:10.2172/1090639.
Till, Andrew T. Thu . "Linear-Multi-Frequency-Grey Preconditioning for Radiative Transfer SN Calculations". United States. doi:10.2172/1090639.
title = {Linear-Multi-Frequency-Grey Preconditioning for Radiative Transfer SN Calculations},
author = {Till, Andrew T.},
abstractNote = {},
doi = {10.2172/1090639},
journal = {},
number = ,
volume = ,
place = {United States},
year = {Thu Aug 15 00:00:00 EDT 2013},
month = {Thu Aug 15 00:00:00 EDT 2013}

Technical Report:

Save / Share:
  • LoopFest V is sponsored by the American Linear Collider Physics Working Group on higher-order calculations (LoopVerein) at a future international linear e{sup +}e{sup -} Collider (ILC). The LoopFest aims to provide a forum to coordinate activities focused on the theoretical challenges resulting from the ultra-high experimental precision of a future ILC.
  • We present a new linear stability analysis of three time discretizations and Monte Carlo interpretations of the nonlinear, grey thermal radiative transfer (TRT) equations: the widely used 'Implicit Monte Carlo' (IMC) equations, the Carter Forest (CF) equations, and the Ahrens-Larsen or 'Semi-Analog Monte Carlo' (SMC) equations. Using a spatial Fourier analysis of the 1-D Implicit Monte Carlo (IMC) equations that are linearized about an equilibrium solution, we show that the IMC equations are unconditionally stable (undamped perturbations do not exist) if {alpha}, the IMC time-discretization parameter, satisfies 0.5 < {alpha} {<=} 1. This is consistent with conventional wisdom. However, wemore » also show that for sufficiently large time steps, unphysical damped oscillations can exist that correspond to the lowest-frequency Fourier modes. After numerically confirming this result, we develop a method to assess the stability of any time discretization of the 0-D, nonlinear, grey, thermal radiative transfer problem. Subsequent analyses of the CF and SMC methods then demonstrate that the CF method is unconditionally stable and monotonic, but the SMC method is conditionally stable and permits unphysical oscillatory solutions that can prevent it from reaching equilibrium. This stability theory provides new conditions on the time step to guarantee monotonicity of the IMC solution, although they are likely too conservative to be used in practice. Theoretical predictions are tested and confirmed with numerical experiments.« less
  • The equations of thermal radiative transfer describe the emission, absorption and transport of photons in a material. As photons travel through the material they are absorbed and re-emitted in a Planckian distribution characterized by the material temperature. As a result of these processes, the material can change resulting in a change in the Planckian emission spectrum. When the coupling between the material and radiation is strong, as occurs when the material opacity or the time step is large, standard iterative techniques converge very slowly. As a result, nested iterative algorithms have been applied to the problem. One algorithm, is tomore » use multifrequency DSA to accelerate the convergence of the multifrequency transport iteration and a grey transport acceleration (GTA) followed by a single group DSA. Here we summarize a new method which uses a grey diffusion equation (GDA) to directly solve the multifrequency transport (S{sub N}) problem. Results of Fourier analysis for both the continuous and discretized equations are discussed and the computational efficiency of GDA is compared with the DSA and GTA nested algorithms. 5 refs., 1 fig., 1 tab.« less
  • In the present study we are studying the effects of including carbon dioxide, ozone, methane, and the halocarbons in addition to water vapor in the radiating atmosphere. The study has focused on two principal issues: the effect on the spectral fluxes and cooling rates of carbon dioxide, ozone and the halocarbons at 1990 concentration levels and the change in fluxes and cooling rates as a consequence of the anticipated ten year change in the profiles of these species. For the latter study the water vapor profiles have been taken as invariant in time. The radiative line-by-line calculations using LBLRTM (Line-By-Linemore » Radiative Transfer Model) have been performed for tropical (TRP), mid-latitude winter (MLW) and mid-latitude summer (MLS) model atmospheres. The halocarbons considered in the present study are CCl{sub 4}, CFC-11, CFC-12 and CFC-22. In addition to considering the radiative effects of carbon dioxide at 355 ppM, the assumed current level, we have also obtained results for doubled carbon dioxide at 710 ppM. An important focus of the current research effort is the effect of the ozone depletion profile on atmospheric radiative effects.« less