An extension of implicit Monte Carlo diffusion: Multigroup and the difference formulation

• Published in: Journal of computational physics 2010-08, Vol.229 (16), p.5707-5723
• Main Authors: Cleveland, Mathew A., Gentile, Nick A., Palmer, Todd S.
• Format: Article
• Language: English
• Subjects: ALGORITHMS; Benchmarking; BENCHMARKS; CLASSICAL AND QUANTUMM MECHANICS, GENERAL PHYSICS; Computational efficiency; Computational techniques; DIFFUSION; Exact sciences and technology; Fences; FREQUENCY DEPENDENCE; IMD; IMPLEMENTATION; Mathematical analysis; Mathematical methods in physics; Mathematical models; Monte Carlo methods; Multigroup; Non-grey; NUMERICAL SOLUTION; Physics; PROBABILITY; RADIANT HEAT TRANSFER; RADIATIONS; REFINING; The difference formulation; THERMAL EQUILIBRIUM
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2010.04.004

Implicit Monte Carlo (IMC) and Implicit Monte Carlo Diffusion (IMD) are approaches to the numerical solution of the equations of radiative transfer. IMD was previously derived and numerically tested on grey, or frequency-integrated problems [1]. In this research, we extend Implicit Monte Carlo Diffusion (IMD) to account for frequency dependence, and we implement the difference formulation [2] as a source manipulation variance reduction technique. We derive the relevant probability distributions and present the frequency dependent IMD algorithm, with and without the difference formulation. The IMD code with and without the difference formulation was tested using both grey and frequency dependent benchmark problems. The Su and Olson semi-analytic Marshak wave benchmark was used to demonstrate the validity of the code for grey problems [3]. The Su and Olson semi-analytic picket fence benchmark was used for the frequency dependent problems [4]. The frequency dependent IMD algorithm reproduces the results of both Su and Olson benchmark problems. Frequency group refinement studies indicate that the computational cost of refining the group structure is likely less than that of group refinement in deterministic solutions of the radiation diffusion methods. Our results show that applying the difference formulation to the IMD algorithm can result in an overall increase in the figure of merit for frequency dependent problems. However, the creation of negatively weighted particles from the difference formulation can cause significant numerical instabilities in regions of the problem with sharp spatial gradients in the solution. An adaptive implementation of the difference formulation may be necessary to focus its use in regions that are at or near thermal equilibrium.