P1 Approximation Theory
The P1 approximation is available as a radiation discretization method in The Radiation in Participating Media Interface.
The P1 approximation is the simplest approximation provided by the method of spherical harmonics method (PN-method). This approximation provides additional accuracy compared to a Rosseland approximation even if it remains a very simple method. The P1 method relies on the following hypotheses:
The media is optically thick media: τ >>1, where τ is the optical thickness defined by the integral of absorption coefficient, κ, along a typical optical path:
From a computational point of view this approximation has a limited impact because it introduces only one additional degree of freedom for G, which is a scalar quantity and adds a heat source or sink to the temperature equation to account for radiative heat transfer contributions. This method, however, fails to accurately represent cases where the radiative intensity propagation dominates over its diffusivity or where the scattering effects cannot be described by a linear isotropic phase function.
The P1 approximation accounts for the radiation transfer equation
by solving following equation for (Ref. 20):
(4-95)
where
DP1is the P1 diffusion coefficient, defined as
a1 is the linear Legendre coefficient of the scattering phase function
Qr is the radiative heat source:
(4-96)
When scattering is modeled as isotropic, a1=0 and the P1 diffusion coefficient reduces to
The following boundary condition applies (Ref. 20):
where qrnet is the net radiative heat flux at the boundary.
Radiation in Participating Media
For the Participating Medium (Radiation in Participating Medium Interface) feature node, the equation Equation 4-95 is implemented.
In addition Qr, defined by Equation 4-96, is added as an heat source in the heat transfer equation:
Opaque Surface
The Opaque Surface (Radiation in Participating Medium and Radiation in Absorbing-Scattering Medium Interfaces) boundary condition defines a boundary opaque to radiation and defines the incident intensity on a boundary:
The Opaque Surface feature accounts for the net radiative heat flux, qr, net, in the heat balance.
Two cases are considered, depending on surface emissivity value:
Gray wall: ε is a number between 0 and 1
Gray Wall
The radiative heat flux at the boundary depends on the surface emissivity, ε:
with
Black Wall
The radiative heat flux at the boundary expression simplifies to
with
Incident Intensity
The Incident Intensity (Radiation in Participating Medium and Radiation in Absorbing-Scattering Medium Interfaces) node defines a boundary that receives incident radiative intensity Iext and that is transparent for outgoing intensity. On these boundaries, the relation between G, qrnet (net radiative heat flux) and Iext (incident radiative intensity) is
by defining
there is
which defines the heat radiative heat flux and also contributes to G boundary condition: