Comparison between performances of Monte Carlo method and method of lines solution of discrete ordinates method

2005 ◽  
Vol 93 (1-3) ◽  
pp. 115-124 ◽  
Author(s):  
Gokmen Demirkaya ◽  
Faruk Arinç ◽  
Nevin Selçuk ◽  
Isil Ayranci
Keyword(s):  
Monte Carlo ◽  
Monte Carlo Method ◽  
Method Of Lines ◽  
Discrete Ordinates ◽  
Discrete Ordinates Method

Experimental and Numerical Studies of Short Pulse Propagation in Model Systems

2002 ◽  
Author(s):  
Sunil Kumar ◽  
Zhixiong Guo ◽  
Janice Aber ◽  
Bruce Garetz
Keyword(s):  
Monte Carlo ◽  
Pulse Propagation ◽  
Short Pulse ◽  
Model Systems ◽  
Discrete Ordinates ◽  
Pulsed Lasers ◽  
Discrete Ordinates Method ◽  
Numerical Studies ◽  
Tissue Phantoms ◽  
Good Agreement

In this paper experimental and numerical studies of the propagation of short-pulsed lasers through scattering and absorbing media are presented. Experimental results of a 60 ps pulse laser transmission in tissue phantoms are obtained and compared with Monte Carlo simulations. Good agreement between the Monte Carlo simulation and experimental measurement is found. Three models are developed for the simulation of short pulse transport. Benchmark comparisons among the Monte Carlo (MC), transient discrete ordinates method (TDOM) and transient radiation element method (TREM) are conducted.

Equation Solving DRESOR Method for Radiative Transfer in Three-Dimensional Isotropically Scattering Media

2014 ◽  
Vol 136 (9) ◽  
Author(s):  
Zhifeng Huang ◽  
Huaichun Zhou ◽  
Guihua Wang ◽  
Pei-feng Hsu
Keyword(s):  
Monte Carlo ◽  
Monte Carlo Method ◽  
Computation Time ◽  
Discrete Ordinates ◽  
Distinct Advantage ◽  
Reverse Monte Carlo ◽  
Scattering Media ◽  
Equation Solving ◽  
Radiative Intensity ◽  
The Monte Carlo Method

Distributions of ratios of energy scattered or reflected (DRESOR) method is a very efficient tool used to calculate radiative intensity with high directional resolution, which is very useful for inverse analysis. The method is based on the Monte Carlo (MC) method and it can solve radiative problems of great complexity. Unfortunately, it suffers from the drawbacks of the Monte Carlo method, which are large computation time and unavoidable statistical errors. In this work, an equation solving method is applied to calculate DRESOR values instead of using the Monte Carlo sampling in the DRESOR method. The equation solving method obtains very accurate results in much shorter computation time than when using the Monte Carlo method. Radiative intensity with high directional resolution calculated by these two kinds of DRESOR method is compared with that of the reverse Monte Carlo (RMC) method. The equation solving DRESOR (ES-DRESOR) method has better accuracy and much better time efficiency than the Monte Carlo based DRESOR (original DRESOR) method. The ES-DRESOR method shows a distinct advantage for calculating radiative intensity with high directional resolution compared with the reverse Monte Carlo method and the discrete ordinates method (DOM). Heat flux comparisons are also given and the ES-DRESOR method shows very good accuracy.

Angular False Scattering in Radiative Heat Transfer Analysis Using the Discrete-Ordinates Method With Higher-Order Quadrature Sets

2013 ◽  
Author(s):  
Brian Hunter ◽  
Zhixiong Guo
Keyword(s):  
Monte Carlo ◽  
Radiative Transfer ◽  
Phase Function ◽  
Radiative Heat ◽  
Higher Order ◽  
Heat Transfer Analysis ◽  
Discrete Ordinates ◽  
Equal Weight ◽  
Participating Media ◽  
Discrete Ordinates Method

The SN quadrature set for the discrete-ordinates method is limited in overall discrete direction number in order to avoid physically unrealistic negative directional weight factors. Such a limitation can adversely impact radiative transfer predictions. Directional discretization results in errors due to ray effect, as well as angular false scattering error due to distortion of the scattering phase function. The use higher-order quadrature schemes in the discrete-ordinates method allows for improvement in discretization errors without an overall directional limitation. In this analysis, four higher-order quadrature sets (Legendre-Equal Weight, Legendre-Chebyshev, Triangle Tessellation, and Spherical Ring Approximation) are implemented for determination of radiative transfer in a 3-D cubic enclosure containing participating media. Radiative heat fluxes, calculated at low direction number, are compared to the SN quadrature and Monte Carlo predictions to gauge quadrature accuracy. Additionally, investigation into the reduction of angular false scattering with sufficient increase in direction number using higher-order quadrature, including heat flux accuracy with respect to Monte Carlo and computational efficiency, is presented. While higher-order quadrature sets are found to effectively minimize angular false scattering error, it is found to be much more computationally efficient to implement proper phase function normalization for accurate radiative transfer predictions.

Comparing the Monte Carlo Method to the Discrete Ordinates Methods in Fluent for Calculating Radiation Heat Transfer in a Small Particle Receiver

2016 ◽  
Author(s):  
Eugene Cho ◽  
Fletcher Miller
Keyword(s):  
Monte Carlo ◽  
Solar Radiation ◽  
Boundary Conditions ◽  
Fortran Program ◽  
Working Fluid ◽  
Discrete Ordinates ◽  
Ansys Fluent ◽  
Discrete Ordinates Method ◽  
Directional Distribution ◽  
Solar Receiver

The Combustion and Solar Energy Laboratory (C&SEL) at San Diego State University is developing a Small Particle Heat Exchange Receiver (SPHER) to absorb and transfer heat from concentrated solar radiation to a working fluid for a gas turbine. The SPHER is to be used with a Concentrated Solar Power (CSP) system where a heliostat field highly concentrates solar radiation on the optical aperture of the SPHER. The solar radiation is volumetrically absorbed by a unique carbon nanoparticle gas mixture within the cavity of the SPHER. This research focuses on comparing a Computational Fluid Dynamics (CFD) model using the ANSYS FLUENT Discrete Ordinates (DO) Model and a program developed by the C&SEL which uses a Monte Carlo Ray Trace (MCRT) method to calculate the spatial and directional distribution of radiation for an idealized solar receiver geometry. Previous research at the C&SEL has shown successful implementation of the MCRT method to calculate the spatial and directional distribution of radiation for an idealized solar receiver geometry. The MCRT method is highly accurate and will serve as the benchmark solution for this research. However the MCRT code takes several days to run, is inflexible to geometry changes, and is cumbersome to implement as the MCRT code needs to be rewritten for each new receiver geometry being considered. These factors necessitate the need to find an alternate method that accurately calculates the spatial and directional distribution of radiation for a solar receiver and can be efficiently implemented for various receiver geometries being studied. The Discrete Ordinates method is a new method for solving the Radiative Transport Equations (RTE) using a FORTRAN program, developed by the C&SEL, and the ANSYS FLUENT Discrete Ordinates model for calculating the RTE. The FORTRAN program calculates the proper inlet radiation boundary conditions that ANSYS FLUENT uses for calculating the RTE. The methodology used for determining the correct CFD mesh, radiative boundary conditions, optimal number of DO theta and phi discretization, as well as the optical properties of the working fluid will be presented in this paper. The main focus of this research is to compare two different methods for solving the Radiative Transport Equations within the idealized SPHER. The solution data for several cases using the previous coupled MCRT method and the ANSYS FLUENT Discrete Ordinates method is presented for both a collimated and diffuse gray radiation approximation. The case studies focus on researching how the MCRT method and Discrete Ordinates method differ when comparing critical receiver parameters such as the mean outlet temperature, wall temperature profile, outlet tube temperature profile, and total receiver efficiency while keeping the total inlet radiation flux of 5 MW and inlet mass flow rate of 5 kg/s constant. This research also presents a study on the optimal Discrete Ordinates angular discretization, as well as a study to determine the solution’s dependence on the number of inlet boundary conditions imposed on the window.

Sign in / Sign up

Export Citation Format

Share Document