A Spectrally Accurate Two-Dimensional Axisymmetric, Tightly Coupled Photon Monte Carlo Radiative Transfer Equation Solver for Hypersonic Entry Flows

2012 ◽  
Vol 134 (12) ◽  
Author(s):  
A. M. Feldick ◽  
M. F. Modest

A two-dimensional axisymmetric ray tracing photon Monte Carlo radiative transfer solver is developed. Like all ray tracing Monte Carlo codes, the ray tracing is performed in 3D, however, arrangements are made to take advantage of the 2D nature of the problem, to minimize computational time. The solver is designed to be tightly integrated into finite volume hypersonic flow solvers and is able to resolve the complex spectral properties of such flows to line-by-line (LBL) accuracy. The solver is then directly integrated into data parallel line relaxation (DPLR), a hypersonic flow solver, and closely coupled calculations are performed.

Author(s):  
Andrew Feldick ◽  
Josh Giegel ◽  
Michael F. Modest

A two-dimensional axisymmetric ray tracing photon Monte Carlo radiative transfer solver is developed. Like all ray tracing Monte Carlo codes, the ray tracing is performed in 3-D, however, arrangements are made to take advantage of the 2-D nature of the problem, to minimize computational time. The solver is designed to be integrated into finite volume hypersonic flow solvers, and is able to resolve the complex spectral properties of such flows to line-by-line accuracy. The solver is then directly integrated into DPLR, a hypersonic flow solver, and closely coupled calculations are performed.


2006 ◽  
Vol 128 (9) ◽  
pp. 945-952 ◽  
Author(s):  
Sandip Mazumder

Two different algorithms to accelerate ray tracing in surface-to-surface radiation Monte Carlo calculations are investigated. The first algorithm is the well-known binary spatial partitioning (BSP) algorithm, which recursively bisects the computational domain into a set of hierarchically linked boxes that are then made use of to narrow down the number of ray-surface intersection calculations. The second algorithm is the volume-by-volume advancement (VVA) algorithm. This algorithm is new and employs the volumetric mesh to advance the ray through the computational domain until a legitimate intersection point is found. The algorithms are tested for two classical problems, namely an open box, and a box in a box, in both two-dimensional (2D) and three-dimensional (3D) geometries with various mesh sizes. Both algorithms are found to result in orders of magnitude gains in computational efficiency over direct calculations that do not employ any acceleration strategy. For three-dimensional geometries, the VVA algorithm is found to be clearly superior to BSP, particularly for cases with obstructions within the computational domain. For two-dimensional geometries, the VVA algorithm is found to be superior to the BSP algorithm only when obstructions are present and are densely packed.


2020 ◽  
Vol 2020 (0) ◽  
pp. 0012
Author(s):  
Tadanobu Sunago ◽  
Michimasa Morita ◽  
Takuma Hori ◽  
Makoto Kashiwagi ◽  
Takuma Shiga

Author(s):  
Xiaodong Lu ◽  
Pei-Feng Hsu

The Monte Carlo (MC) method has been widely used to solve radiative transfer problems due to its flexibility and simplicity in simulating the energy transport process in arbitrary geometries with complex boundary conditions. However, the major drawback of the conventional (or forward) Monte Carlo method is the long computational time for converged solution. Reverse or backward Monte Carlo (RMC) is considered as an alternative approach when solutions are only needed at certain locations and time. The reverse algorithm is similar to the conventional method, except that the energy bundle (photons ensemble) is tracked in a time-reversal manner. Its migration is recorded from the detector into the participating medium, rather than from the source to the detector as in the conventional MC. There is no need to keep track of the bundles that do not reach a particular detector. Thus, RMC method takes up much less computation time than the conventional MC method. On the other hand, RMC will generate less information about the transport process as only the information at the specified locations, e.g., detectors, is obtained. In the situation where detailed information of radiative transport across the media is needed the RMC may not be appropriate. RMC algorithm is most suitable for diagnostic applications where inverse analysis is required, e.g., optical imaging and remote sensing. In this study, the development of a reverse Monte Carlo method for transient radiative transfer is presented. The results of non-emitting, absorbing, and anisotropically scattering media subjected to an ultra short light pulse irradiation are compared with the forward Monte Carlo and discrete ordinates methods results.


Sign in / Sign up

Export Citation Format

Share Document