Chandrasekhar's H-function

In atmospheric radiation, Chandrasekhar's ''H''-function appears as the solutions of problems involving scattering, introduced by the Indian American astrophysicist Subrahmanyan Chandrasekhar.Sparrow, Ephraim M., and Robert D. Cess. "Radiation heat transfer." Series in Thermal and Fluids Engineering, New York: McGraw-Hill, 1978, Augmented ed. (1978). The Chandrasekhar's ''H''-function $H(\mu)$ defined in the interval $0\leq\mu\leq 1$, satisfies the following nonlinear integral equation

:$H(\mu) = 1+\mu H(\mu)\int_0^1 \fracH(\mu') \, d\mu'$

where the characteristic function $\Psi(\mu)$ is an even polynomial in $\mu$ satisfying the following condition

:$\int_0^1\Psi(\mu) \, d\mu \leq \frac$.

If the equality is satisfied in the above condition, it is called ''conservative case'', otherwise ''non-conservative''. Albedo is given by $\omega_o= 2\Psi(\mu) = \text$. An alternate form which would be more useful in calculating the ''H'' function numerically by iteration was derived by Chandrasekhar as,

: \frac = \left -2\int_0^1\Psi(\mu) \, d\mu\right + \int_0^1 \fracH(\mu') \, d\mu'.

In conservative case, the above equation reduces to

:$\frac= \int_0^1 \fracH(\mu')d\mu'$.

Approximation

The ''H'' function can be approximated up to an order $n$ as

: $H(\mu) = \frac\frac$

where $\mu_i$ are the zeros of Legendre polynomials $P_$ and $k_\alpha$ are the positive, non vanishing roots of the associated characteristic equation

:$1 = 2 \sum_^n \frac$

where $a_j$ are the quadrature weights given by

:$a_j = \frac\int_^1 \frac \, d\mu_j$

Explicit solution in the complex plane

In complex variable $z$ the ''H'' equation is

: $H(z) = 1- \int_0^1 \frac z H(\mu)\Psi(\mu) \, d\mu, \quad \int_0^1 , \Psi(\mu), \, d\mu \leq \frac, \quad \int_0^\delta , \Psi(\mu), \, d\mu \rightarrow 0, \ \delta\rightarrow 0$

then for $\Re (z)>0$, a unique solution is given by

:$\ln H(z) = \frac \int_^ \ln T(w) \frac \, dw$

where the imaginary part of the function $T(z)$ can vanish if $z^2$ is real i.e., $z^2 = u+iv = u\ (v=0)$. Then we have

: $T(z) = 1- 2 \int_0^1 \Psi(\mu) \, d\mu - 2 \int_0^1 \frac \, d\mu$

The above solution is unique and bounded in the interval $0\leq z\leq 1$ for conservative cases. In non-conservative cases, if the equation $T(z)=0$ admits the roots $\pm 1/k$, then there is a further solution given by

:$H_1(z) = H(z) \frac$

Properties

*\int_0^1 H(\mu)\Psi(\mu) \, d\mu = 1-\left -2\int_0^1\Psi(\mu) \, d\mu \right. For conservative case, this reduces to $\int_0^1 \Psi(\mu)d\mu=\frac$.
*\left -2\int_0^1\Psi(\mu) \, d\mu\right \int_0^1 H(\mu) \Psi(\mu) \mu^2 \, d\mu + \frac \left int_0^1 H(\mu)\Psi(\mu)\mu \, d\mu\right2 = \int_0^1 \Psi(\mu)\mu^2 \, d\mu. For conservative case, this reduces to \int_0^1 H(\mu)\Psi(\mu) \mu d\mu = \left \int_0^1 \Psi(\mu)\mu^2d\mu\right.
*If the characteristic function is $\Psi(\mu)=a+b\mu^2$, where $a, b$ are two constants(have to satisfy $a+b/3\leq 1/2$) and if $\alpha_n = \int_0^1 H(\mu)\mu^n \, d\mu, \ n\geq 1$ is the nth moment of the ''H'' function, then we have
:$\alpha_0 = 1 + \frac (a\alpha_0^2 + b \alpha_1^2)$
and
:$(a+b\mu^2) \int_0^1\frac\,d\mu'=\frac-b(\alpha_1-\mu\alpha_0)$

See also

* Chandrasekhar's X- and Y-function

External links

*MATLAB function to calculate the ''H'' function https://www.mathworks.com/matlabcentral/fileexchange/29333-chandrasekhar-s-h-function

References

{{Reflist

Special functions
Integral equations
Scattering
Scattering, absorption and radiative transfer (optics)