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.

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 )\underset{0}{\overset{1}{\int }}\frac{\mathrm{\Psi }({\mu }^{\prime })}{\mu +{\mu }^{\prime }}H({\mu }^{\prime })d{\mu }^{\prime }$

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

$\underset{0}{\overset{1}{\int }}\mathrm{\Psi }(\mu )d\mu \le \frac{1}{2}$

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

L'isotropie correspond à

$2\mathrm{\Psi }={\omega }_0$

où $0\le {\omega }_0\le 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 )}={[1-2\underset{0}{\overset{1}{\int }}\mathrm{\Psi }(\mu )d\mu ]}^{1/2}+\underset{0}{\overset{1}{\int }}\frac{{\mu }^{\prime }\mathrm{\Psi }({\mu }^{\prime })}{\mu +{\mu }^{\prime }}H({\mu }^{\prime })d{\mu }^{\prime }$

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

Propriétés

• $\underset{0}{\overset{1}{\int }}H(\mu )\mathrm{\Psi }(\mu )d\mu =1-{[1-2\underset{0}{\overset{1}{\int }}\mathrm{\Psi }(\mu )d\mu ]}^{1/2}$

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

$\underset{0}{\overset{1}{\int }}\mathrm{\Psi }(\mu )d\mu =\frac{1}{2}$

• ${[1-2\underset{0}{\overset{1}{\int }}\mathrm{\Psi }(\mu )d\mu ]}^{1/2}\underset{0}{\overset{1}{\int }}H(\mu )\mathrm{\Psi }(\mu ){\mu }^{2}d\mu +\frac{1}{2}{[\underset{0}{\overset{1}{\int }}H(\mu )\mathrm{\Psi }(\mu )\mu d\mu ]}^2=\underset{0}{\overset{1}{\int }}\mathrm{\Psi }(\mu ){\mu }^{2}d\mu$

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

$\underset{0}{\overset{1}{\int }}H(\mu )\mathrm{\Psi }(\mu )\mu d\mu ={[2\underset{0}{\overset{1}{\int }}\mathrm{\Psi }(\mu ){\mu }^{2}d\mu ]}^{1/2}$.

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

$M_0=1+\frac{1}{2}(a{M}_0^2+b{M}_1^2)$

et

$(a+b{\mu }^2)\underset{0}{\overset{1}{\int }}\frac{H({\mu }^{\prime })}{\mu +{\mu }^{\prime }}d{\mu }^{\prime }=\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-\underset{0}{\overset{1}{\int }}\frac{z}{z+\mu }H(\mu )\mathrm{\Psi }(\mu )d\mu ,\underset{0}{\overset{1}{\int }}|\mathrm{\Psi }(\mu )|d\mu \le \frac{1}{2},\lim {\mathrm{\backslash limits}}_{\delta \to 0}\underset{0}{\overset{\delta }{\int }}|\mathrm{\Psi }(\mu )|d\mu =0$

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

$\ln H(z)=\frac{1}{2\pi i}\underset{-i\mathrm{\infty }}{\overset{+i\mathrm{\infty }}{\int }}\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\underset{0}{\overset{1}{\int }}\mathrm{\Psi }(\mu )d\mu -2\underset{0}{\overset{1}{\int }}\frac{{\mu }^2\mathrm{\Psi }(\mu )}{u-{\mu }^2}d\mu$

Dans le cas conservatif $0\le z\le 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\mathrm{\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_{2^{\prime }n}({\mu }_j)}\underset{-1}{\overset{1}{\int }}\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 HModèle:,.

Voir aussi

• Méthode de Wiener-Hopf
• Problème de Milne

Liens externes

• Code de calcul Matlab accessible sur le site suivant : Modèle:Site web