Piecewise linear discretization of Symbolic Implicit Monte Carlo radiation transport in the difference formulation
Abstract
We describe a Monte Carlo solution for time dependent photon transport, in the difference formulation with the material in local thermodynamic equilibrium, that is piecewise linear in its treatment of the material state variable. Our method employs a Galerkin solution for the material energy equation while using Symbolic Implicit Monte Carlo to solve the transport equation. In constructing the scheme, one has the freedom to choose between expanding the material temperature, or the equivalent black body radiation energy density at the material temperature, in terms of finite element basis functions. The former provides a linear treatment of the material energy while the latter provides a linear treatment of the radiative coupling between zones. Subject to the conditional use of a lumped material energy in the vicinity of strong gradients, possible with a linear treatment of the material energy, our approach provides a robust solution for time dependent transport of thermally emitted radiation that can address a wide range of problems. It produces accurate results in thick media.
 Authors:
 University of California, Lawrence Livermore National Laboratory, P.O. Box 808, Livermore, CA 94550 (United States). Email: brooks3@llnl.gov
 University of California, Lawrence Livermore National Laboratory, P.O. Box 808, Livermore, CA 94550 (United States)
 Publication Date:
 OSTI Identifier:
 20840372
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Computational Physics; Journal Volume: 220; Journal Issue: 1; Other Information: DOI: 10.1016/j.jcp.2006.07.014; PII: S00219991(06)003433; Copyright (c) 2006 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; BLACKBODY RADIATION; ENERGY DENSITY; EQUATIONS; FINITE ELEMENT METHOD; LTE; MONTE CARLO METHOD; PHOTON TRANSPORT; TIME DEPENDENCE; TRANSPORT THEORY
Citation Formats
Brooks, Eugene D., Szoke, Abraham, and Peterson, Jayson D.L. Piecewise linear discretization of Symbolic Implicit Monte Carlo radiation transport in the difference formulation. United States: N. p., 2006.
Web. doi:10.1016/j.jcp.2006.07.014.
Brooks, Eugene D., Szoke, Abraham, & Peterson, Jayson D.L. Piecewise linear discretization of Symbolic Implicit Monte Carlo radiation transport in the difference formulation. United States. doi:10.1016/j.jcp.2006.07.014.
Brooks, Eugene D., Szoke, Abraham, and Peterson, Jayson D.L. Wed .
"Piecewise linear discretization of Symbolic Implicit Monte Carlo radiation transport in the difference formulation". United States.
doi:10.1016/j.jcp.2006.07.014.
@article{osti_20840372,
title = {Piecewise linear discretization of Symbolic Implicit Monte Carlo radiation transport in the difference formulation},
author = {Brooks, Eugene D. and Szoke, Abraham and Peterson, Jayson D.L.},
abstractNote = {We describe a Monte Carlo solution for time dependent photon transport, in the difference formulation with the material in local thermodynamic equilibrium, that is piecewise linear in its treatment of the material state variable. Our method employs a Galerkin solution for the material energy equation while using Symbolic Implicit Monte Carlo to solve the transport equation. In constructing the scheme, one has the freedom to choose between expanding the material temperature, or the equivalent black body radiation energy density at the material temperature, in terms of finite element basis functions. The former provides a linear treatment of the material energy while the latter provides a linear treatment of the radiative coupling between zones. Subject to the conditional use of a lumped material energy in the vicinity of strong gradients, possible with a linear treatment of the material energy, our approach provides a robust solution for time dependent transport of thermally emitted radiation that can address a wide range of problems. It produces accurate results in thick media.},
doi = {10.1016/j.jcp.2006.07.014},
journal = {Journal of Computational Physics},
number = 1,
volume = 220,
place = {United States},
year = {Wed Dec 20 00:00:00 EST 2006},
month = {Wed Dec 20 00:00:00 EST 2006}
}

We describe a Monte Carlo solution for time dependent photon transport, in the difference formulation with the material in local thermodynamic equilibrium (LTE), that is piecewise linear in its treatment of the material state variable. Our method employs a Galerkin solution for the material energy equation while using Symbolic Implicit Monte Carlo (SIMC) to solve the transport equation. In constructing the scheme, one has the freedom to choose between expanding the material temperature, or the equivalent black body radiation energy density at the material temperature, in terms of finite element basis functions. The former provides a linear treatment of themore »

Symbolic implicit Monte Carlo radiation transport in the difference formulation: a piecewise constant discretization
The equations of radiation transport for thermal photons are notoriously difficult to solve in thick media without resorting to asymptotic approximations such as the diffusion limit. One source of this difficulty is that in thick, absorbing media, thermal emission and absorption are almost completely balanced. A new formulation for thermal radiation transport, called the difference formulation, was recently introduced in order to remove the stiff balance between emission and absorption. In the new formulation, thermal emission is replaced by derivative terms that become small in thick media. It was proposed that the difficulties of solving the transport equation in thickmore » 
Accurate and efficient radiation transport in optically thick media  by means of the Symbolic Implicit Monte Carlo method in the difference formulation
The equations of radiation transport for thermal photons are notoriously difficult to solve in thick media without resorting to asymptotic approximations such as the diffusion limit. One source of this difficulty is that in thick, absorbing media thermal emission is almost completely balanced by strong absorption. In a previous publication [SB03], the photon transport equation was written in terms of the deviation of the specific intensity from the local equilibrium field. We called the new form of the equations the difference formulation. The difference formulation is rigorously equivalent to the original transport equation. It is particularly advantageous in thick media,more » 
Implicit Monte Carlo with a linear discontinuous finite element material solution and piecewise nonconstant opacity
Here, the nonlinear thermal radiativetransfer equations can be solved in various ways. One popular way is the Fleck and Cummings Implicit Monte Carlo (IMC) method. The IMC method was originally formulated with piecewiseconstant material properties. For domains with a coarse spatial grid and large temperature gradients, an error known as numerical teleportation may cause artificially noncausal energy propagation and consequently an inaccurate material temperature. Source tilting is a technique to reduce teleportation error by constructing subspatialcell (or subcell) emission profiles from which IMC particles are sampled. Several source tilting schemes exist, but some allow teleportation error to persist. We examinemore »Cited by 1 
Asymptotic analysis of the spatial discretization of radiation absorption and reemission in Implicit Monte Carlo
We perform an asymptotic analysis of the spatial discretization of radiation absorption and reemission in Implicit Monte Carlo (IMC), a Monte Carlo technique for simulating nonlinear radiative transfer. Specifically, we examine the approximation of absorption and reemission by a spatially continuous artificialscattering process and either a piecewiseconstant or piecewiselinear emission source within each spatial cell. We consider three asymptotic scalings representing (i) a time step that resolves the meanfree time, (ii) a Courant limit on the timestep size, and (iii) a fixed time step that does not depend on any asymptotic scaling. For the piecewiseconstant approximation, we show that onlymore »