Acceleration schemes for the discrete ordinates method

Discrete ordinates (DO) methods have been developed to solve the multidimensional radiative transport equation (RTE) for applications including combustion processes and other combined-mode heat transfer cases. The convergence of DO methods is known to degrade for optical thicknesses greater than unity, which occur for example in a flame. Acceleration schemes have been developed for use in neutron transport applications, but little work has been done to accelerate convergence of the RTE for radiative heat transfer applications. This article presents several acceleration schemes for the RTE, including successive overtaxation, synthetic acceleration, and mesh rebalance methods. Solution convergence is discussed and demonstrated using two- and three-dimensional examples. Although all methods improve convergence, the mesh rebalance method improves the RTE convergence best. For some conditions, the rebalance method improves convergence dramatically, reducing RTE iterations by an order of magnitude. However, the mesh rebalance method fails to produce convergence of the RTE for large optical thicknesses and fine mesh discretizations. Examples are used to demonstrate that unproved convergence can be obtained by solving the rebalance equation on a coarser grid, which is determined by regrouping the base RTE grid, until an optical thickness of near unity is obtained on the coarse grid. Based on these findings, a general solution strategy is discussed. / /

[1]  R. E. Alcouffe,et al.  Diffusion synthetic acceleration methods for the diamond-differenced discrete-ordinates equations , 1977 .

[2]  Tae-Kuk Kim,et al.  Effect of anisotropic scattering on radiative heat transfer in two-dimensional rectangular enclosures , 1988 .

[3]  J. P. Jessee,et al.  Comparison of discrete ordinates formulations for radiative heat transfer in multidimensional geometries , 1995 .

[4]  S. Patankar,et al.  Finite volume method for radiation heat transfer , 1994 .

[5]  Bart W. Stuck,et al.  A Computer and Communication Network Performance Analysis Primer (Prentice Hall, Englewood Cliffs, NJ, 1985; revised, 1987) , 1987, Int. CMG Conference.

[6]  Philip J. Smith,et al.  Predicting Radiative Transfer in Axisymmetric Cylindrical Enclosures Using the Discrete Ordinates Method , 1988 .

[7]  IMPLICIT SOLUTION SCHEME TO IMPROVE CONVERGENCE RATE IN RADIATIVE TRANSFER PROBLEMS , 1992 .

[8]  W. H. Reed.  The Effectiveness of Acceleration Techniques for Iterative Methods in Transport Theory , 1971 .

[9]  W. Fiveland.  Discrete-Ordinates Solutions of the Radiative Transport Equation for Rectangular Enclosures , 1984 .

[10]  D. Brandt,et al.  Multi-level adaptive solutions to boundary-value problems math comptr , 1977 .

[11]  William L. Briggs,et al.  A multigrid tutorial , 1987 .

[12]  O. Widlund,et al.  Iterative solution of elliptic systems : and applications to the neutron diffusion equations of reactor physics , 1967 .

[13]  J. P. Jessee,et al.  Bounded, High-Resolution Differencing Schemes Applied to the Discrete Ordinates Method , 1997 .

[14]  J. W. Ruge,et al.  4. Algebraic Multigrid , 1987 .

[15]  J. Truelove,et al.  Discrete-Ordinate Solutions of the Radiation Transport Equation , 1987 .

[16]  E. Lewis,et al.  Computational Methods of Neutron Transport , 1993 .

[17]  M. Modest.  Radiative heat transfer , 1993 .

[18]  G. D. Raithby,et al.  A Finite-Volume Method for Predicting a Radiant Heat Transfer in Enclosures With Participating Media , 1990 .

[19]  W. A. Fiveland,et al.  Three-dimensional radiative heat-transfer solutions by the discrete-ordinates method , 1988 .

[20]  J. R. Howell,et al.  Two-Dimensional Radiation in Absorbing-Emitting Media Using the P-N Approximation , 1983 .

[21]  J. Gillis,et al.  Matrix Iterative Analysis , 1961 .

[22]  M. Pinar Mengüç,et al.  Radiation heat transfer in combustion systems , 1987 .

[23]  K. D. Lathrop,et al.  DISCRETE ORDINATES ANGULAR QUADRATURE OF THE NEUTRON TRANSPORT EQUATION , 1964 .