Comparison of spherical harmonics method and discrete ordinates method for radiative transfer in a turbulent jet flame
Journal Article
·
· Journal of Quantitative Spectroscopy and Radiative Transfer
- Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)
- Marquette University, Milwaukee, WI (United States)
- University of California, Merced, CA (United States)
Here, in this study, we systematically compared the accuracy and computational cost of two popular solution methods for the radiative transfer equation (RTE): the spherical harmonics method (PN) and the discrete ordinates method (DOM). We first investigated convergence characteristics of different orders of PN and DOM in a series of 1D homogeneous configurations with varying optical thicknesses. Both solvers perform better for optically thicker cases. The accuracy of PN methods increases with its order, , but the gain in accuracy reduces with the increase in , i.e., improvement of P7 over P5 is less than that of P3 over P1. This decreasing trend becomes more prominent as the optical thickness decreases. On the other hand, DOM’s accuracy increases almost linearly with the increase in the number of ordinates (or polar angles in this study) in all cases. While comparing the directional profile of radiative intensity, both solvers perform better when the radiative intensity is more isotropic. These solvers were then connected with a full spectrum k-distribution (FSK) spectral model and used to perform radiation-coupled simulations of a turbulent jet flame in an axi-symmetric cylindrical domain. Results obtained from P1 to P7 approximations for PN, and 2 x 4, 4 x 4, 4 x 8, 8 x 8 finite angles for DOM are compared with that from an optically thin model, and a reference solution from line-by-line (LBL) photon Monte Carlo (PMC) method. The choice of radiation solver shows a noticeable impact on the temperature distribution of the flame. The PN solvers lead to slightly higher radiant fractions and the DOM solvers lead to slightly lower radiant fractions than the PMC benchmark solution. Finally, the computational costs of each of these solvers are also reported and an intermittent evaluation / time blending scheme to improve the computational efficiency of radiation solvers in radiation-coupled simulations are also demonstrated.
- Research Organization:
- Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)
- Sponsoring Organization:
- Exascale Computing Project; National Science Foundation (NSF); USDOE National Nuclear Security Administration (NNSA)
- Grant/Contract Number:
- AC05-00OR22725
- OSTI ID:
- 1906591
- Alternate ID(s):
- OSTI ID: 1908857
- Journal Information:
- Journal of Quantitative Spectroscopy and Radiative Transfer, Journal Name: Journal of Quantitative Spectroscopy and Radiative Transfer Journal Issue: / Vol. 296; ISSN 0022-4073
- Publisher:
- ElsevierCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Similar Records
Comparison of Radiation Models for a Turbulent Piloted Methane/Air Jet Flame: A Frozen-Field Study
Acceleration schemes for the discrete ordinates method
Assessment of the finite-volume method and the discrete ordinate method for radiative heat transfer in a three-dimensional rectangular enclosure
Conference
·
Tue Jun 01 00:00:00 EDT 2021
·
OSTI ID:1807265
Acceleration schemes for the discrete ordinates method
Conference
·
Sat Dec 30 23:00:00 EST 1995
·
OSTI ID:455357
Assessment of the finite-volume method and the discrete ordinate method for radiative heat transfer in a three-dimensional rectangular enclosure
Journal Article
·
Thu Dec 31 23:00:00 EST 1998
· Numerical Heat Transfer. Part B, Fundamentals
·
OSTI ID:328379