Research Papers: Radiative Heat Transfer

Implementation of High-Order Spherical Harmonics Methods for Radiative Heat Transfer on openfoam

[+] Author and Article Information
Wenjun Ge

School of Engineering,
University of California,
Merced, CA 95343
e-mail: wge@ucmerced.edu

Ricardo Marquez

School of Engineering,
University of California,
Merced, CA 95343
e-mail: rmarquez3@ucmerced.edu

Michael F. Modest

Fellow ASME
School of Engineering,
University of California,
Merced, CA 95343
e-mail: mmodest@ucmerced.edu

Somesh P. Roy

School of Engineering,
University of California,
Merced, CA 95343
e-mail: sroy3@ucmerced.edu

1Corresponding author.

Contributed by the Heat Transfer Division of ASME for publication in the JOURNAL OF HEAT TRANSFER. Manuscript received August 19, 2014; final manuscript received December 29, 2014; published online February 10, 2015. Assoc. Editor: Zhuomin Zhang.

J. Heat Transfer 137(5), 052701 (May 01, 2015) (9 pages) Paper No: HT-14-1548; doi: 10.1115/1.4029546 History: Received August 19, 2014; Revised December 29, 2014; Online February 10, 2015

A general formulation of the spherical harmonics (PN) methods was developed recently to expand the method to high orders of PN. The set of N(N + 1)/2 three-dimensional second-order elliptic PDEs formulation and their Marshak boundary conditions for arbitrary geometries are implemented in the openfoam finite volume based cfd software. The results are verified for four cases, including a 1D slab, a 2D square enclosure, a 3D cylindrical enclosure, and an axisymmetric flame. All cases have strongly varying radiative properties, and the results are compared with exact solutions and solutions from the photon Monte Carlo method (PMC).

Copyright © 2015 by ASME
Your Session has timed out. Please sign back in to continue.


Chandrasekhar, S., 1960, Radiative Transfer, Dover, New York.
Lee, C. E., 1962, “The Discrete Sn Approximation to Transport Theory,” Lawrence Livermore Laboratory, Technical Information Series Report No. LA2595.
Chai, J. C., Lee, H. S., and Patankar, S. V., 1994, “Finite Volume Method for Radiation Heat Transfer,” J. Thermophys. Heat Transfer, 8(3), pp. 419–425. [CrossRef]
Howell, J. R., 1968, “Application of Monte Carlo to Heat Transfer Problems,” Advances in Heat Transfer, Vol. 5, J. P.Hartnett and T. F.Irvine, eds., Academic, New York. [CrossRef]
Jeans, J. H., 1917, “The Equations of Radiative Transfer of Energy,” Mon. Not. R. Astron. Soc., 78(1), pp. 28–36. [CrossRef]
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, 24(4), pp. 373–389. [CrossRef]
Modest, M. F., 2013, Radiative Heat Transfer, 3rd ed., Academic, New York.
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. Transfer, 109(9), pp. 1641–1666. [CrossRef]
Yang, J., and Modest, M. F., 2007, “High-Order P-N Approximation for Radiative Transfer in Arbitrary Geometries,” J. Quant. Spectrosc. Radiat. Transfer, 104(2), pp. 217–227. [CrossRef]
Modest, M. F., 2012, “Further Developments of the Elliptic PN-Approximation Formulation and Its Marshak Boundary Conditions,” Numer. Heat Transfer, Part B, 62(2–3), pp. 181–202. [CrossRef]
Jasak, H., Jemcov, A., and Tukovic, Z., 2007, “OpenFOAM: A C++ Library for Complex Physics Simulations,” International Workshop on Coupled Methods in Numerical Dynamics, IUC, pp. 1–20.
Marquez, R., and Modest, M. F., 2013, “Implementation of the PN-Approximation for Radiative Heat Transfer on OpenFOAM,” ASME Paper No. HT2013-17556. [CrossRef]
Marshak, R. E., 1947, “Note on the Spherical Harmonics Method as Applied to the Milne Problem for a Sphere,” Phys. Rev., 71(7), pp. 443–446. [CrossRef]
Varshalovich, D. A., Moskalev, A. N., and Khersonskii, V. K., 1981, Quantum Theory of Angular Momentum, World Scientific, Singapore.
Ravishankar, M., Mazumder, S., and Kumar, A., 2010, “Finite-Volume Formulation and Solution of the P3 Equations of Radiative Transfer on Unstructured Meshes,” ASME J. Heat Transfer, 132(2), p. 023402. [CrossRef]
Nocedal, J., and Wright, S. J., 2006, Numerical Optimization, 2nd ed., Springer, Berlin.
Press, W. H., 2007, Numerical Recipes: The Art of Scientific Computing, 3rd ed., Cambridge University, Cambridge, UK.
Barlow, R. S., “International Workshop on Measurement and Computation of Turbulent Nonpremixed Flames (TNF),” http://www.sandia.gov/TNF/abstract.html
Modest, M. F., 2003, “Narrow-Band and Full-Spectrum k-Distributions for Radiative Heat Transfer—Correlated-k vs. Scaling Approximation,” J. Quant. Spectrosc. Radiat. Transfer, 76(1), pp. 69–83. [CrossRef]
Modest, M. F., and Riazzi, R. J., 2005, “Assembly of Full-Spectrum k-Distributions From a Narrow-Band Database; Effects of Mixing Gases, Gases and Nongray Absorbing Particles, and Mixtures With Nongray Scatterers in Nongray Enclosures,” J. Quant. Spectrosc. Radiat. Transfer, 90(2), pp. 169–189. [CrossRef]
Wang, A., and Modest, M. F., 2007, “Spectral Monte Carlo Models for Nongray Radiation Analyses in Inhomogeneous Participating Media,” Int. J. Heat Mass Transfer, 50(19–20), pp. 3877–3889. [CrossRef]
Ferziger, J. H., and Perić, M., 2001, Computational Methods for Fluid Dynamics, 3rd ed., Springer, Berlin.


Grahic Jump Location
Fig. 1

Definition of Euler angles for an arbitrary rotation

Grahic Jump Location
Fig. 2

Rotations of 1D slab at angles φ=0,45 deg,-45 deg in the x-y plane. The medium properties increase with low values at the lower wall to higher values at the upper wall.

Grahic Jump Location
Fig. 3

Comparing incident radiation and radiative heat source with analytical solutions of PN and exact solution of 1D slab. (a) Incident radiation G and (b) radiative heat source ∇·q.

Grahic Jump Location
Fig. 4

Incident radiation and radiative source for a square enclosure with various PN approximations. (a) Ck = 1, (b) Ck = 0.1, and (c) Ck = 0.01.

Grahic Jump Location
Fig. 5

Incident radiation for a cylinder enclosure at two axial locations with various PN approximations. (a) z = 0.71 and (b) z = 1.20.

Grahic Jump Location
Fig. 6

Cylinder mesh of scaled Flame D in the analysis

Grahic Jump Location
Fig. 7

Mean temperature and mass fraction fields for scaled Sandia Flame D

Grahic Jump Location
Fig. 8

Radiative source for Sandia Flame D × 4 at two axial locations as calculated with various PN approximations. (a) z = 1.00 m and (b) z = 1.43 m.



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In