Library Subscription: Guest
Home Archives Officers Future meetings American Society of Thermal and Fluids Engineering

AN EFFICIENT MONTE CARLO-BASED SOLVER FOR THERMAL RADIATION IN PARTICIPATING MEDIA

Joseph A. Farmer
Department of Mechanical Engineering, Marquette University, Milwaukee, WI, USA

Somesh P. Roy
Department of Mechanical Engineering, Marquette University, Milwaukee, WI, USA.

DOI: 10.1615/TFEC2019.rad.027584
pages 1565-1573


KEY WORDS: Radiation, Monte Carlo, HPC, Convergence, Low-discrepancy sequence

Abstract

Monte Carlo-based solvers, while well-suited for accurate calculation of complex thermal radiation transport problems in participating media, are often deemed computationally unattractive for use in the solution of real-world problems. The main disadvantage of Monte Carlo (MC) solvers is their slow convergence rate and relatively high computational cost. This work presents a novel approach based on a low-discrepancy sequence (LDS) and is proposed for reducing the error bound of a Monte Carlo-based radiation solver. Sobols sequence – an LDS generated with a bit-by-bit exclusive-or operator – is used to develop a quasi-Monte Carlo (QMC) solver for thermal radiation in this work. Preliminary results for simple radiation problems in participating media show that the QMC-based solver has a lower error than the conventional MC-based solver. At the same time, QMC does not add any significant computational overhead. This essentially leads to a lower computational cost to achieve similar error levels from the QMC-based solver than the MC-based solver for thermal radiation.

Purchase $20.00 Check subscription Publication Ethics and Malpractice Recommend to my Librarian Bookmark this Page