Transfert de rayonnement - Définition

Source: Wikipédia sous licence CC-BY-SA 3.0.
La liste des auteurs de cet article est disponible ici.

Géométrie sphérique

Écrire l'équation de transfert en géométrie sphérique

Les modèles d'étoiles les plus simples font l'hypothèse que ces objets possèdent la symétrie sphérique, c'est-à-dire qu'ils sont constitués d'un empilement de couches sphériques individuellement homogènes. Les quantités physiques, dont l'intensité du rayonnement en un point, ne dépendent alors spatialement parlant que de la variable r, distance du point au centre de l'étoile. Dans ce modèle l'intensité est une fonction du point choisi et de l'angle θ compris entre 0 et π donnant la direction du faisceau lumineux par rapport au rayon. Formellement :

I_\nu\,=\,I_\nu(r,\; \theta)

Le schéma ci-contre indique la géométrie utilisée et permet d'écrire l'équation de transfert.

La dérivée de l'intensité le long du trajet du photon est :

\frac{dI}{ds}= \frac{\partial I}{\partial r} \frac{dr}{ds} + \frac{\partial I}{\partial \theta} \frac{d\theta}{ds}\ .

Comme on le voit sur la figure on a les relations

\begin{cases} dr/ds\,=\cos\theta\\ d\theta/ds\,= -(\sin\theta/r)\ , \end{cases}

de sorte que l'équation de transfert devient

\cos\theta\frac{\partial I_\nu}{\partial r}\,-\,\frac{\sin\theta}{r}\frac{\partial I_\nu}{\partial\theta}= -\kappa_\nu\,I_\nu + \epsilon_\nu \ .

Intégrons cette équation sur l'angle solide

d\omega\,=\,2\pi\,\sin\theta\,d\theta\,.

Le deuxième terme de gauche de l'équation conduit à une intégrale qui peut s'intégrer par parties comme :

2\pi\,\int_0^{\pi} (\sin\theta)^2(\partial I_\nu/\partial \theta)\,d\theta = 2\pi\,\int_0^\pi 2 I_\nu \sin\theta\,\cos\theta\,d\theta \equiv 2 H_\nu\,,

H est le flux défini plus haut.

Le membre de gauche de l'équation de transfert intégré sur tout l'angle solide donne donc

\frac{dH_\nu}{dr} - \frac{2}{r}H_\nu\,\equiv \frac{1}{r^2}\frac{d}{dr}(r^2H_\nu)\,.

L'équation de transfert complète intégrée sur tout l'angle solide s'écrit donc

\frac{1}{r^2}\frac{d}{dr}(r^2 H_\nu) = 4\pi (\epsilon_\nu - \kappa_\nu\,J_\nu)\ ,

J est l'intensité moyenne du rayonnement.

Si l'équilibre radiatif est réalisé cela signifie par définition qu'en intégrant sur toutes les fréquences l'énergie absorbée est égale à l'énergie émise. On a donc

4\pi \int_0^\infty (\epsilon_\nu - \kappa_\nu\,J_\nu)\,=\,0\ .

En considérant alors le membre de gauche on obtient le résultat

H \equiv \int_0^\infty H_\nu\,d\nu\,= \frac{C}{r^2}\,,

dans lequel C est une constante, indiquant que le flux par unité de surface varie comme l'inverse du carré du rayon de la couche traversée. Physiquement cela signifie que le flux total (4π r2 H ) franchissant les couches de l'étoile reste constant, ce qui est normal si ces dernières ne contiennent aucune source d'énergie.

Température effective d'une étoile

En résumé, la température effective d'une étoile est la quantité T telle que l'énergie totale par unité de temps émise par l'astre s'exprime sous la forme

E \,=\, 4\pi R^2 \sigma T_{{\rm eff}}^4\,,

R est le rayon de l'étoile supposée sphérique. Autrement dit, par définition

T_{\rm eff}\,=\,[E/(4\pi R^2\sigma)]^{1/4}\,.

La quantité E se mesure en erg/s dans le système CGS ou en watt en unités SI.

Comme application numérique calculons la température effective du Soleil. Cet astre émet une luminosité totale

E \equiv L_{\odot} =3,826 \times 10^{26}\;{\rm W}

et son rayon est

R_\odot\,=\,6,96\times10^{8}\ {\rm m}\,.

Avec σ = 5,67×10-8 W⋅m-2⋅K -4 on trouve

(T_\odot)_{\rm eff}\,=\,5770\ {\rm K}\,.

La température effective d'une étoile est un paramètre servant à mesurer la quantité totale d'énergie lumineuse émise par cet astre. On la définit comme la température du corps noir qui émettrait dans l'espace la même quantité d'énergie. Comme une étoile ne rayonne pas comme un corps noir, cette quantité n'est pas une « vraie » température. Elle a néanmoins un certain rapport avec la température des couches superficielles de l'étoile.

Considérons un corps noir de température T. Un tel corps noir émet une intensité de lumière I indépendante de la direction égale à la fonction de Planck, soit :

I_\nu\,=\,B_\nu (T)\,.

D'après la définition de l'intensité I le flux à la fréquence ν émergeant d'une ouverture de surface unité faite dans la cavité du corps noir est

F_\nu\,=\, 2\pi\,\int_0^1 I_\nu (\mu)\,\mu d\mu = 2\pi B_\nu(T)\int_0^1\mu d\mu = \pi B_\nu(T)\,.

Le flux total émis est donc

F \,=\,\int_0^\infty F_\nu d\nu = \pi \int_0^\infty B_\nu (T) d\nu = \sigma T^4\,,

où σ est la constante de Stefan-Boltzmann. En sens inverse, à flux donné, la température donnée par cette relation sera par définition la température effective de l'étoile. Pour obtenir le flux total émis par toute l'étoile, on intègrera encore cette quantité sur toute sa surface.

Calculons maintenant directement le flux en intégrant l'équation de transfert. Dans le cadre de l'approximation d'Eddington, la fonction source du rayonnement est donnée en fonction de la profondeur optique par la fonction linéaire

S(\nu, \,\tau_\nu)\,=\,S_0(\nu)\,[1+(3/2)\tau_\nu]\,.

L'intensité émergente à la fréquence ν dans la direction arccosμ est égale à

I_\nu (\mu) \,=\,\int_0^\infty S(\nu, t) e^{-t/\mu}\,dt/\mu = \int_0^\infty S_0(\nu)[1+(3/2)t]\, e^{-t/\mu}\,dt/\mu= S_0(\nu)[1 + (3/2)\mu]\,.

Le flux traversant une surface unité est

F_\nu \,=\, 2\pi\int_0^\infty S_0(\nu) [1 + (3/2)\mu] \,\mu d\mu = 2 \pi S_0 (\nu)\,.
Diagramme de Hertzsprung-Russell créé par Richard Powell et diffusé sur Wikipédia avec sa permission. La luminosité totale d'une étoile est portée en fonction de la température effective de sorte que sur une même verticale l'ordonnée d'un point mesure la surface de l'étoile.

Si on fait de plus l'hypothèse de l'ETL, la fonction source S(ν) à la surface est égale à la fonction de Planck B(T) à la température T de surface (qui maintenant, dans le cadre de l'hypothèse de l'équilibre thermodynamique local, est une « vraie » température), de sorte que l'on obtient en intégrant sur les fréquences

F\,=\, \int_0^\infty F_\nu \,d\nu = 2 \pi \int_0^\infty B_\nu (T_0)\,d\nu = 2 \sigma T_0^4\,.

En comparant les deux expressions du flux à travers une surface unité, on voit que

T_{\rm eff}^4\,=\,2\,T_0^4

et par conséquent

T_0/T_{\rm eff}\,=\,(1/2)^{1/4} = 0,841\,.

Par comparaison la résolution exacte du cas gris donne

T_0/T_{\rm eff}\,=\,(\sqrt{3}/4)^{1/4} = 0,8114

ce qui montre que l'approximation d'Eddington est assez exacte (au moins dans certaines circonstances). On voit que le paramètre de température effective est assez voisin de la température du gaz de surface et possède de ce fait une certaine signification physique (surtout si on se contente d'ordres de grandeur) mais cela n'implique pas que le rayonnement émergent lui-même suive la loi de Planck.

Pour reconnaître une étoile dans le ciel et en déduire ses paramètres physiques, les astronomes utilisent sa température effective et sa luminosité totale. Lorsque l'on porte les étoiles dans un diagramme basé sur ces deux paramètres (le justement célèbre diagramme de Hertzsprung-Russell, ou « diagramme HR »), les points ne se répartissent pas au hasard mais selon des séquences définies (comme la séquence principale). L'analyse de la structure de ce diagramme a permis d'établir la théorie de l'évolution d'une étoile. Les points ayant même abscisse dans ce diagramme ont même température effective et par conséquent, par définition de cette température effective, même énergie lumineuse émise par unité de surface. Il en résulte que l'ordonnée du point représentatif d'une étoile sur une même verticale est une simple mesure du carré du rayon de l'astre. Les étoiles les plus grosses (géantes ou supergéantes) sont situées dans le haut du diagramme et les plus petites (naines), dans la partie basse.

Page générée en 0.112 seconde(s) - site hébergé chez Contabo
Ce site fait l'objet d'une déclaration à la CNIL sous le numéro de dossier 1037632
A propos - Informations légales
Version anglaise | Version allemande | Version espagnole | Version portugaise