Abstract

The method of spherical harmonics (or PN) is a popular method for approximate solution of the radiative transfer equation (RTE) in participating media. A rigorous conservative finite-volume (FV) procedure is presented for discretization of the P3 equations of radiative transfer in two-dimensional geometry—a set of four coupled, second-order partial differential equations. The FV procedure presented here is applicable to any arbitrary unstructured mesh topology. The resulting coupled set of discrete algebraic equations are solved implicitly using a coupled solver that involves decomposition of the computational domain into groups of geometrically contiguous cells using the binary spatial partitioning algorithm, followed by fully implicit coupled solution within each cell group using a preconditioned generalized minimum residual solver. The RTE solver is first verified by comparing predicted results with published Monte Carlo (MC) results for two benchmark problems. For completeness, results using the P1 approximation are also presented. As expected, results agree well with MC results for large/intermediate optical thicknesses, and the discrepancy between MC and P3 results increase as the optical thickness is decreased. The P3 approximation is found to be more accurate than the P1 approximation for optically thick cases. Finally, the new RTE solver is coupled to a reacting flow code and demonstrated for a laminar flame calculation using an unstructured mesh. It is found that the solution of the four P3 equations requires 14.5% additional CPU time, while the solution of the single P1 equation requires 9.3% additional CPU time over the case without radiation.

1.
Modest
,
M. F.
, 2003,
Radiative Heat Transfer
, 2nd ed.,
Academic
,
New York
.
2.
Siegel
,
R.
, and
Howell
,
J.
, 2001,
Thermal Radiation Heat Transfer
, 4th ed.,
Taylor and Francis
,
London
.
3.
Olfe
,
D. B.
, 1967, “
A Modification of the Differential Approximation for Radiative Transfer
,”
AIAA J.
0001-1452,
5
(
4
), pp.
638
643
.
4.
Modest
,
M. F.
, 1989, “
The Modified Differential Approximation for Radiative Transfer in General Three-Dimensional Media
,”
J. Thermophys. Heat Transfer
0887-8722,
3
(
3
), pp.
283
288
.
5.
Modest
,
M. F.
, 1990, “
The Improved Differential Approximation for Radiative Heat Transfer in Multi-Dimensional Media
,”
ASME J. Heat Transfer
0022-1481,
112
, pp.
819
821
.
6.
Fiveland
,
W.
, and
Jamaluddin
,
A.
, 1991, “
Three-dimensional Spectral Radiative Heat Transfer Solutions by the Discrete Ordinates Method
,”
J. Thermophys. Heat Transfer
0887-8722,
5
(
3
), pp.
335
339
.
7.
Chai
,
J. C.
,
Lee
,
H. S.
, and
Patankar
,
S. V.
, 1994, “
Finite-Volume Method for Radiative Heat Transfer
,”
J. Thermophys. Heat Transfer
0887-8722,
8
, pp.
419
425
.
8.
Yuen
,
W.
, and
Takara
,
E.
, 1994, “
Development of a General Zonal Method for Analysis of Radiative Transfer in Absorbing and Anisotropically Scattering Media
,”
Numer. Heat Transfer, Part B
1040-7790,
25
, pp.
75
96
.
9.
Shen
,
Z.
,
Smith
,
T.
, and
Hix
,
P.
, 1983, “
Linearization of Thermal Radiation Terms for Improved Convergence by Use of the Zone Method
,”
Numer. Heat Transfer, Part B
1040-7790,
6
, pp.
377
382
.
10.
Smith
,
T.
, and
Severin
,
S.
, 1985, “
Development of the Finite Analytic With Radiation Method
,”
ASME J. Heat Transfer
0022-1481,
107
, pp.
735
737
.
11.
Chai
,
J. C.
,
Lee
,
H. S.
, and
Patankar
,
S. V.
, 1993, “
Ray Effect and False Scattering in the Discrete Ordinates Method
,”
Numer. Heat Transfer, Part B
1040-7790,
24
, pp.
373
389
.
12.
Mathur
,
S. R.
, and
Murthy
,
J. Y.
, 2001, “
Acceleration of Anisotropic Scattering Computations Using Coupled Ordinates Method (COMET)
,”
ASME J. Heat Transfer
0022-1481,
123
(
3
), pp.
607
612
.
13.
Mazumder
,
S.
, 2005, “
A New Numerical Procedure for Coupling Radiation in Participating Media With Other Modes of Heat Transfer
,”
ASME J. Heat Transfer
0022-1481,
127
(
9
), pp.
1037
1045
.
14.
Yang
,
J.
, and
Modest
,
M. F.
, 2007, “
Elliptic PDE Formulation of General, Three-Dimensional High-Order PN-Approximations for Radiative Transfer
,”
J. Quant. Spectrosc. Radiat. Transf.
0022-4073,
104
(
2
), pp.
217
227
.
15.
Modest
,
M. F.
, and
Yang
,
J.
, 2008, “
Elliptic PDE Formulation and Boundary Conditions of the Spherical Harmonics Method of Arbitrary Order for General Three-Dimensional Geometries
,”
J. Quant. Spectrosc. Radiat. Transf.
0022-4073,
109
, pp.
1641
1666
.
16.
Bayazitoglu
,
Y.
, and
Higenyi
,
J.
, 1979, “
The Higher-Order Differential Equations of Radiative Transfer: P3 Approximation
,”
AIAA J.
0001-1452,
17
, pp.
424
431
.
17.
Ratzel
,
A. C.
, and
Howell
,
J. R.
, 1983, “
Two-Dimensional Radiation in Absorbing-Emitting-Scattering Media Using the P-N Approximation
,”
ASME J. Heat Transfer
0022-1481,
105
, pp.
333
340
.
18.
Menguc
,
M. P.
, and
Viskanta
,
R.
, 1985, “
Radiative Transfer in Three-Dimensional Rectangular Enclosures Containing Inhomogeneous, Anisotropically Scattering Media
,”
J. Quant. Spectrosc. Radiat. Transf.
0022-4073,
33
(
6
), pp.
533
549
.
19.
Ou
,
S. C.
, and
Liou
,
K. N.
, 1982, “
Generalization of the spherical Harmonic Method to Radiative Transfer in Multi-Dimensional Space
,”
J. Quant. Spectrosc. Radiat. Transf.
0022-4073,
28
(
4
), pp.
271
288
.
20.
McClarren
,
R. G.
,
Evans
,
T. M.
,
Lowrie
,
R. B.
, and
Densmore
,
J. D.
, 2008, “
Semi-Implicit Time Integration for PN Thermal Radiative Transfer
,”
J. Comput. Phys.
0021-9991,
227
, pp.
7561
7586
.
21.
McClarren
,
R. G.
, 2007, “
Spherical Harmonics Methods for Thermal Radiation Transport
,” Ph.D. thesis, University of Michigan, Ann Arbor, MI.
22.
Cheney
,
W.
, and
Kincaid
,
W.
, 1985,
Numerical Mathematics and Computing
, 2nd ed.,
Brooks-Cole
,
Belmont, MA
.
23.
Ferziger
,
J.
, and
Peric
,
M.
, 1999,
Computational Methods for Fluid Dynamics
, 2nd ed.,
Springer-Verlag
,
Berlin
.
24.
Kumar
,
A.
, and
Mazumder
,
S.
, “
Coupled Solution of the Species Conservation Equations Using Unstructured Finite-Volume Method
,”
Int. J. Numer. Methods Fluids
0271-2091, to be published.
25.
Kumar
,
A.
, and
Mazumder
,
S.
, 2007, “
A Low-Memory Block-Implicit Solver for Coupled Solution of the Species Conservation Equations on an Unstructured Mesh
,”
Proceedings of the International Mechanical Engineering Congress and Exposition IMECE2007
, Paper No. 2007-42517.
26.
Sung
,
K.
, and
Shirley
,
P.
, 1992, “
Ray Tracing With a BSP Tree
,”
Computer Graphics Gems III
,
D.
Kirk
, ed.,
AP Professional
,
San Diego, CA
, pp.
271
274
.
27.
Karypis
,
G.
, and
Kumar
,
V.
, 1998, “
A Fast and High Quality Multilevel Scheme for Partitioning Irregular Graphs
,”
SIAM J. Sci. Comput. (USA)
1064-8275,
20
, pp.
359
392
.
28.
Saad
,
Y.
, 2003,
Iterative Methods for Sparse Linear Systems
, 2nd ed.,
SIAM
,
Philadelphia, PA
.
29.
Selle
,
L.
,
Lartigue
,
G.
,
Poinsot
,
T.
,
Koch
,
R.
,
Schildmacher
,
K. U.
,
Krebs
,
W.
,
Prade
,
B.
,
Kaufmann
,
P.
, and
Veynante
,
D.
, 2004, “
Compressible Large Eddy Simulation of Turbulent Combustion in Complex Geometry on Unstructured Meshes
,”
Combust. Flame
0010-2180,
137
, pp.
489
505
.
You do not currently have access to this content.