Full Content is available to subscribers

Subscribe/Learn More  >

Finite-Volume Solution of the P3 Equations of Radiative Transfer and Coupling to Reacting Flow Calculations

[+] Author Affiliations
Mahesh Ravishankar, Sandip Mazumder, Ankan Kumar

The Ohio State University, Columbus, OH

Paper No. HT2009-88014, pp. 1-11; 11 pages
  • ASME 2009 Heat Transfer Summer Conference collocated with the InterPACK09 and 3rd Energy Sustainability Conferences
  • Volume 1: Heat Transfer in Energy Systems; Thermophysical Properties; Heat Transfer Equipment; Heat Transfer in Electronic Equipment
  • San Francisco, California, USA, July 19–23, 2009
  • Conference Sponsors: Heat Transfer Division
  • ISBN: 978-0-7918-4356-7 | eISBN: 978-0-7918-3851-8
  • Copyright © 2009 by ASME


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 pre-conditioned Generalized Minimum Residual (GMRES) solver. The RTE solver is first verified by comparing predicted results with published Monte Carlo (MC) results for a benchmark problem. 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 4 P3 equations requires 14.5% additional CPU time, while the solution of the single P1 equation requires 9.3% additional CPU time over the 10 equations that are solved for the reacting flow calculations.

Copyright © 2009 by ASME



Interactive Graphics


Country-Specific Mortality and Growth Failure in Infancy and Yound Children and Association With Material Stature

Use interactive graphics and maps to view and sort country-specific infant and early dhildhood mortality and growth failure data and their association with maternal

Citing articles are presented as examples only. In non-demo SCM6 implementation, integration with CrossRef’s "Cited By" API will populate this tab (http://www.crossref.org/citedby.html).

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 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