An acceleration method is proposed particularly for the $P1$ equation. The radiative energy balance is used as a constraint to correct iterative solutions. The method not only accelerates convergence but also preserves the radiative energy balance, the latter being of great importance when radiation calculations are coupled with flow calculations. This acceleration method can be applied to other elliptical problems with boundary conditions of the second and/or the third kind.

1.
Modest, M. F., 1993, Radiative Heat Transfer, McGraw-Hill, New York.
2.
FLUENT, 1995, “Computational Fluid Dynamics Software,” Version 4.
3.
Patankar
,
S. V.
,
1981
, “
A Calculation Procedure for Two-Dimensional Elliptic Situations
,”
Numer. Heat Transfer
,
4
, pp.
409
425
.
4.
Ferziger, J. H., 1998, Numerical Methods for Engineering Application, John Wiley & Sons, New York.
5.
Wachspress, E. L., 1966, Iterative Solution of Elliptic Systems and Applications to the Neutron Diffusion Equations of Reactor Physics, Prentice-Hall, Englewood Cliffs, NJ.
6.
Reed
,
W. H.
,
1971
, “
The Effectiveness of Acceleration Techniques for Iterative Methods in Transport Theory
,”
Nucl. Sci. Eng.
,
45
, No.
3
, pp.
245
254
.
7.
Fiveland
,
W. A.
and
Jessee
,
J. P.
,
1996
, “
Acceleration Schemes for the Discrete Ordinates Method
,”
J. Thermophys. Heat Transfer
,
10
, No.
3
, pp.
445
451
.