Moyennes spectrales
Un problème majeur du transfert radiatif est le nombre extrêmement élevé de raies spectrales d'un gaz, pouvant atteindre ou dépasser plusieurs centaines de milliers, en sus du fond continu. Si l'on ne s'intéresse qu'au seul transfert d'énergie les moyennes spectrales sont utilisées pour réduire notablement le calcul. Il existe diverses moyennes adaptées au cas étudié.
Définitions
modifierOn s'intéresse à un milieu sans diffusion, décrit par l'équation de transfert radiatif stationnaire portant sur la luminance spectrale (ou spectrique[N 1]) dans la direction donnée par la variable Ω[1],[2]
où est le coefficient d'absorption à la fréquence ν et la distribution de Planck écrite en fonction de l'énergie réduite
La luminance peut être réduite à la connaissance de dans le cas du problème unidimensionnel de la propagation dans une direction choisie, s étant la coordonnée mesurant cette direction
Dans le cas d'un milieu homogène (T constant) la solution avec le terme source est
- Le premier terme à droite correspond à l'absorption du terme source (loi de Beer-Lambert). Il tend à disparaître lorsque s devient grand. La fraction de L0 transmise est la transmittance spectrale
- Si le milieu est constitué de n composants les coefficients d'absorption s'additionnent et les transmittances se multiplient
- le second terme correspond à l'émission du milieu. Lorsque ce terme tend vers zéro et la solution du problème vers . Il fait apparaître l'absorbance spectrale
- On a
Moyenne de Rosseland
modifierOn suppose ici que le milieu est proche de l'équilibre thermodynamique, condition réalisée pour un milieu optiquement épais tel que le libre parcours moyen soit petit devant toute longueur caractéristique des gradients présents dans le milieu
On écrit la solution de l'équation de transfert comme le développement au premier ordre
avec
Cette expression s'annule au maximum de la distribution de Planck donnée par
où LambertW est la fonction W de Lambert.
Si on intègre cette expression sur tout le spectre il vient
On réécrit le dernier terme en faisant apparaître la moyenne de Rosseland
On notera que est négatif au-delà du maximum de la fonction de Planck, ce qui privilégie les fortes fréquences dans la moyenne de Rosseland. Corrélativement, augmenter l'absorption aux faibles fréquences diminue cette moyenne.
Un exemple pour lequel le calcul est analytique est le rayonnement continu de freinage (bremẞtrahlung) dans un plasma pour lequel on a, en négligeant l'émission induite
où T est ici la température de translation des électrons. Ceci donne par intégration
Moyennes pour l'absorption et l'émission
modifierEn intégrant l'équation de transfert sur tout le spectre on obtient
où L est la luminance totale
Le premier terme correspondant à l'émission ne pose pas de problème : il suffit d'introduire la moyenne de Planck
où a est la constante radiative.
Le second terme n'est pas connu a priori puisqu'il contient la solution cherchée. On va donc chercher des approximations possibles sous la forme , ce qui constitue une opération délicate puisque luminance et coefficient d'absorption sont fortement corrélés. En fait on utilise un découpage en bandes de largeur variable suivant la méthode adoptée.
Moyenne par bande
modifierOn définit un coefficient d'absorption moyen pris sur un intervalle spectral [3]. Généralement il s'agit de la moyenne de Planck
On suppose ensuite que la luminance moyenne solution du problème dans l'intervalle obéit à une équation analogue à la précédente
où
On voit que la solution tend vers la solution exacte lorsque le nombre d'intervalles tend vers l'infini. En pratique il en faut un grand nombre pour obtenir une bonne précision[4].
Modèles de « bandes étroites »
modifierCe type de modèle a été développé pour traiter de manière aussi précise que possible les raies spectrales[5],[6].
Définissons l'absorbance totale sur un intervalle de fréquence
Connaissant le profil d'une raie il est possible de calculer l'absorbance de cette raie. On peut également connaître les caractéristiques d'un ensemble de raies si on choisit bien son intervalle, lequel est forcément « étroit ». Le choix, pour un gaz donné, du type de raie, de leur espacement, de leur répartition en hauteur, de leur recouvrement (ou non) est le choix du physicien, les options étant nombreuses[2],[5].
Pour le cas où plusieurs composants coexistent on montre aisément que
Cette méthode est limitée par la connaissance précise des caractéristiques a priori du spectre dans chaque intervalle de fréquence et conduit donc à un nombre élevé de bandes.
Méthode de la fonction de distribution
modifierCette méthode, connue en anglais sous le nom de opacity distribution function (ODF) ou cumulative k-distribution (CK method) consiste à rassembler toutes les contributions dans une bande au centre de celle-ci. Dès lors le spectre est caractérisé par la seule distribution statistique des coefficients d'absorption compris entre κ et κ + dκ, ceci pour chaque bande.
La résolution de l'équation de transfert radiatif étant effectué par une méthode quelconque pour chaque échantillon de coefficient (issu de la discrétisation de la fonction de distribution) dans chaque bande il ne reste plus qu'à pondérer le résultat par la densité de probabilité associée.
Une telle méthode permet une précision raisonnable avec un nombre de bandes faible et une discrétisation raisonnable de la fonction de distribution (typiquement 10 x 10)[4].
Notes
modifier- Le vocabulaire est défini par l'ISO 80000-7[1].
- Le terme « opacité » est une traduction littérale de l'anglais « opacity ». Il est couramment utilisé en physique des plasmas et en astrophysique où il désigne ou , ρ étant la masse volumique du milieu.
Références
modifier- (en) Subrahmanyan Chandrasekhar, Radiative transfer, Dover Publications, (ISBN 0486-6059-06, lire en ligne)
- (en) John R. Howell, R. Siegel et M. Pinar Mengüç, Thermal Radiation Heat Transfer, CRC Press, (ISBN 978-1-4398-9455-2)
- (en) Dimitri Mihalas et Barbara Weibel Mihalas, Foundations of Radiation Hydrodynamics, Oxford University Press, (ISBN 0-19-503437-6, lire en ligne)
- (en) A. A. Wray, J.-F. Ripoll et D. Prabhu, « Investigation of the opacity binning approach for solving the shock-generated radiation of the Apollo AS-501 re-entry », Center for Turbulence Research, Proceedings of the Summer Program 2006, (lire en ligne)
- (en) Richard M. Goody et Yuk Ling Yung, Atmospheric Radiation. Theoretical Basis, Oxford University Press, (ISBN 0-19-510291-6)
- (en) W. F. Huebner et W. D. Barfield, Opacity, Springer, (ISBN 978-1-4614-8796-8)