Fonction H de Chandrasekhar

La fonction H de Chandrasekhar est utilisée pour la résolution du problème de transfert radiatif unidimensionnel dans un milieu absorbant et diffusant. Elle est définie par une équation intégrale établie par Viktor Ambartsumian et Subrahmanyan Chandrasekhar^[1].

Définition

Fonction H pour diverses valeurs de l'albédo.

La fonction introduite par Subrahmanyan Chandrasekhar est généralement définie par l'équation intégrale établie par Viktor Ambartsumian

H(\mu )=1+\mu H(\mu )\int _{0}^{1}{\frac {\Psi (\mu ')}{\mu +\mu '}}H(\mu ')\,d\mu '

où $\Psi (\mu )$ est une fonction caractéristique décrivant la diffusion dans le milieu. C'est un polynôme pair satisfaisant

\int _{0}^{1}\Psi (\mu )\,d\mu \leq {\frac {1}{2}}

Le cas correspondant à la limite haute est dit conservatif (il conserve la densité de flux d'énergie).

L'isotropie correspond à

2\Psi =\omega _{0}

où $0\leq \omega _{0}\leq 1$ est l'albédo, constant. $\omega _{0}=1$ correspond au cas de la diffusion pure.

Une définition équivalente plus utilisée pour l'évaluation numérique s'écrit

{\frac {1}{H(\mu )}}=\left[1-2\int _{0}^{1}\Psi (\mu )\,d\mu \right]^{1/2}+\int _{0}^{1}{\frac {\mu '\Psi (\mu ')}{\mu +\mu '}}H(\mu ')\,d\mu '

Dans le cas conservatif le premier terme de l'équation ci-dessus s'annule.

Propriétés

$\int _{0}^{1}H(\mu )\Psi (\mu )\,d\mu =1-\left[1-2\int _{0}^{1}\Psi (\mu )\,d\mu \right]^{1/2}$

Dans le cas conservatif cette équation se réduit à

\int _{0}^{1}\Psi (\mu )d\mu ={\frac {1}{2}}

$\left[1-2\int _{0}^{1}\Psi (\mu )\,d\mu \right]^{1/2}\int _{0}^{1}H(\mu )\Psi (\mu )\mu ^{2}\,d\mu +{\frac {1}{2}}\left[\int _{0}^{1}H(\mu )\Psi (\mu )\mu \,d\mu \right]^{2}=\int _{0}^{1}\Psi (\mu )\mu ^{2}\,d\mu$

Dans le cas conservatif cette équation se réduit à

\int _{0}^{1}H(\mu )\Psi (\mu )\mu d\mu =\left[2\int _{0}^{1}\Psi (\mu )\mu ^{2}d\mu \right]^{1/2}

.

Pour une fonction caractéristique correspondante à la diffusion Thomson ou Rayleigh $\Psi (\mu )=a+b\mu ^{2}$ où $a,b$ sont deux constantes satisfaisant $a+b/3\leq 1/2$ et si on définit le moment d'ordre $n$ par $M_{n}=\int _{0}^{1}H(\mu )\mu ^{n}\,d\mu ,\ n\geq 0$ alors

M_{0}=1+{\frac {1}{2}}(aM_{0}^{2}+bM_{1}^{2})

et

(a+b\mu ^{2})\int _{0}^{1}{\frac {H(\mu ')}{\mu +\mu '}}\,d\mu '={\frac {H(\mu )-1}{\mu H(\mu )}}-b(M_{1}-\mu M_{0})

Solution dans le plan complexe

En utilisant la variable complexe $z$ l'équation de définition de H s'écrit

H(z)=1-\int _{0}^{1}{\frac {z}{z+\mu }}H(\mu )\Psi (\mu )\,d\mu ,\quad \int _{0}^{1}|\Psi (\mu )|\,d\mu \leq {\frac {1}{2}},\quad \lim \limits _{\delta \to 0}\int _{0}^{\delta }|\Psi (\mu )|\,d\mu =0

Dans le plan $\Re (z)>0$ la solution est

\ln H(z)={\frac {1}{2\pi i}}\int _{-i\infty }^{+i\infty }\ln T(w){\frac {z}{w^{2}-z^{2}}}\,dw

où la partie imaginaire de $T(z)$ s'annule si $z^{2}$ est réel, c'est-à-dire si $z^{2}=u+iv\equiv u$ . On a alors

T(z)=1-2\int _{0}^{1}\Psi (\mu )\,d\mu -2\int _{0}^{1}{\frac {\mu ^{2}\Psi (\mu )}{u-\mu ^{2}}}\,d\mu

Dans le cas conservatif $0\leq z\leq 1$ la solution est unique. Dans le cas contraire $T(z)=0$ admet les racines $\pm {\frac {1}{k}}$ . Il existe donc une solution donnée par

H_{1}(z)=H(z){\frac {1+kz}{1-kz}}

Approximation

Le développement suivant particulièrement connu car il est à la base de la méthode SN

H(\mu )={\frac {1}{\mu _{1}\cdots \mu _{n}}}{\frac {\prod _{i=0}^{n}(\mu +\mu _{i})}{\prod _{\alpha }(1+k_{\alpha }\mu )}}

où les $\mu _{i}$ sont les racines des polynômes de Legendre $P_{2n}$ et les $k_{\alpha }$ les solutions strictement positives de l'équation caractéristique

2\sum _{j=1}^{n}{\frac {a_{j}\Psi (\mu _{j})}{1-k^{2}\mu _{j}^{2}}}

Les $a_{j}$ sont les poids de la quadrature donnés par

a_{j}={\frac {1}{P_{2n}'(\mu _{j})}}\int _{-1}^{1}{\frac {P_{2n}(\mu _{j})}{\mu -\mu _{j}}}\,d\mu _{j}

D'une façon générale il existe diverses méthodes pour le calcul numérique des fonctions H^[2]^,^[3].

Références

↑ (en) Subrahmanyan Chandrasekhar, Radiative transfer, Dover Publications, 1960 (ISBN 0486-6059-06, lire en ligne)
↑ Rabindra Nath Das et Rasajit Kumar Bera, « Numerical evaluation of Chandrasekhar’s H-function, its first and second differential coefficients, its pole and moments from the new form for plane parallel scattering atmosphere in radiative transfer », sur ArXiv
↑ (en) P. B. Bosma et W. A. de Rooij, « Efficient Methods to Calculate Chandrasekhar's H-Functions », Astronomy and Astrophysics, vol. 126,‎ 1983, p. 283-292 (lire en ligne)

(en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Chandrasekhar's H-function » (voir la liste des auteurs).

Voir aussi

Liens externes

Code de calcul Matlab accessible sur le site suivant : Mike Simcock, « Chandrasekhar's H function », sur MathWorks

[Chandrasekhar-1] (en) Subrahmanyan Chandrasekhar, Radiative transfer, Dover Publications, 1960 (ISBN 0486-6059-06, lire en ligne)

[2] Rabindra Nath Das et Rasajit Kumar Bera, « Numerical evaluation of Chandrasekhar’s H-function, its first and second differential coefficients, its pole and moments from the new form for plane parallel scattering atmosphere in radiative transfer », sur ArXiv

[3] (en) P. B. Bosma et W. A. de Rooij, « Efficient Methods to Calculate Chandrasekhar's H-Functions », Astronomy and Astrophysics, vol. 126,‎ 1983, p. 283-292 (lire en ligne)

[1]

[2]

[3]