UNIVERSITE DE YAOUNDE 1 UNIVERSITY OF YAOUNDE 1
FACUTE DES SCIENCES FACULTY OF SCIENCE DEPARTEMENT
DE PHYSIQUE DEPARTEMENT OF PHYSICS
LABORATOIRE : ENERGIE ET ENVIRONNEMENT
NIVEAU : MASTER 2
UE : PHYS 5229/ PROPAGATION DU RAYONNEMENT SOLAIRE DANS
L’ATMOSPHERE
RESPONSABLE DE L’UE : Dr. WANDJI NYAMSI WILLIAM
GROUPE D’EXPOSE : 1
Thème : Résolution de l’ETR -
Méthode des ordonnées discrètes.
ARMEL ZAMBOU KENFACK 17O2884
𝑀𝐵𝑂𝑈𝑂𝑀𝐵𝑂𝑈𝑂 𝐹𝐼𝐹𝐸𝑁 𝑀𝐴𝑆𝐻𝑂𝑈𝐷𝐼 15𝑄2292
𝑁𝑍𝑂𝑈𝑁𝑇𝑂𝑈 𝑁𝐺𝐴𝑀𝑊𝐴 𝑃𝐴𝑈𝐿𝐼𝑁 15𝐻2278
𝑀𝐴𝐺𝑁𝐸 𝑇𝐴𝐿𝐿𝐴 𝑆𝑇𝐸𝑃𝐻𝐴𝑁𝐼𝐸 𝐶𝐻𝐴𝑅𝑀𝐸𝐿𝐿𝐸 16𝐽2339
NOMS ET MARICULE DES ETUIDIANTS :
𝐵𝐴𝑌𝐴𝑂𝐿𝐴 𝑆𝑇𝐸𝑉𝐸 15𝐼2832
𝑇𝑆𝑂𝐹𝐴𝐶𝐾 𝑇𝐼𝑂𝐺𝑂 𝐺𝐴𝐸𝐿 𝐹𝑅𝐴𝑁𝐶𝐼𝑆 14𝑁2180
𝑁𝑍𝑂𝐾𝑂𝑈 𝑁𝐺𝐴𝑇𝑆𝐼𝑁𝐺 𝑅𝑂𝐷𝑅𝐼𝐺𝑈𝐸 15𝐻2934
[ 𝑇𝐴𝐶𝐻𝐼𝐷𝐴 𝑈𝐿𝑅𝐼𝐶𝐻 𝐷𝐸𝑅𝐾𝑆𝐸𝑁 13𝐿2552
ANNEE ACADEMIQUE : 2020/2021
INTRODUCTION
Au fur et à mesure de l’évolution technologique, nous constatons que
l’influence de la thermique est essentielle de par les différents modes de transfert de
chaleur. De plus, le transfert radiatif est important du fait qu’il régit la propagation du
rayonnement dans un milieu. De ce fait, l’équation de transfert radiatif (ETR) qui
caractérise le rayonnement nécessite une résolution afin d’en déterminer les
grandeurs physiques telles que la luminance, l’éclairement, la transmittance et bien
d’autres. La résolution de l’ETR peut se faire de façon analytique mais dans la plupart
des cas notamment en tenant compte du contexte réel, elle se fait de façon
numérique grâce aux codes de transfert radiatif (CTR). Dans ces codes sont
implémentés différentes méthodes ou approximation pour la résolution numérique
de l’ETR parmi lesquelles nous avons la méthode des ordonnées discrètes (MOD). La
méthode des ordonnées discrètes, introduite par Chandrasekhar au début des
années 1960, puis améliorée par Fiveland et Truelove, consiste à discrétiser l’espace
angulaire en un nombre fini de directions pour représenter la dépendance angulaire
de la luminance. La discrétisation angulaire conduit à transformer l’équation de
transfert radiatif en un jeu d’équations, chacune correspondant à une direction
discrète. La méthode des ordonnées discrètes s’apparente donc à une méthode aux
différences finies pour la dépendance angulaire de la luminance. L’intégrale de
diffusion sur toutes les directions de la sphère est remplacée par une somme discrète
sur toutes les directions. La méthode aux ordonnées discrètes consiste donc à réécrire
l’équation de transfert radiatif sous la forme d’un système de n équations
correspondant aux différentes directions.
I. PRINCIPE DE LA MOD :
La MOD introduit une quadrature angulaire, permettant d’obtenir un système
d’équations intégro-différentielle, dont la taille dépend du nombre de directions
discrètes de propagation du rayonnement. L’utilisation d’une quadrature angulaire
implique l’emploi d’un poids directionnel ωk associé à chaque direction discrète
correspondant à l’angle solide de rang k. De plus, la luminance est supposée
constante sur chaque angle solide discret et le calcul de la luminance dans chaque
direction discrète est effectué par chacune de ces équations. Les directions de
propagation sont discrétisées en un nombre fini de directions. Les équations
intégro-differentielles et se mettent alors sous la forme d’un système d’équations
différentielles à coefficients constants. Elle permet de passer de l’équation de
transfert radiatif, qui est de type intégro-différentielle à un système d’équation
RESOLUTION DE L’ETR-MOD
2
algébriques obtenues en remplaçant la représentation directionnelle par un
ensemble discret de directions.
II. DISCRETISATION SPATIALE
On considère une enceinte bidimensionnelle de géométrie quelconque contenant
un milieu semi-transparent. Le système de coordonnées est schématisé sur la figure
1. On intègre I'ETR sur un élément i de surface 𝑆𝑖 (figure 2). On a :
(a)
Figure 1 : Système de coordonnées
L'application du théorème de Gauss dans le cas d'un maillage quadrangulaire
donne :
(b)
Les intensités moyennes et centrales sont données respectivement par :
RESOLUTION DE L’ETR-MOD
3
(c)
Figure 2 : Volume de contrôle
III. DISCRETISATION ANGULAIRE
Dans la méthode des ordonnées discrètes, on remplace l'intégrale par des
quadratures pondérées : (d)
Où 𝝎𝒊 les sont les poids associés aux directions 𝛀𝒊 , Donc I ‘équation de transfert
radiatif se ramène à :
(e)
De la même façon, l'équation des conditions aux limites s'écrit :
(f)
RESOLUTION DE L’ETR-MOD
4
De plus, l'équation de Ir divergence du flux devient :
(g)
IV. METHODE NUMERIQUE ET ALGORITHME
En équilibre radiatif, le champ de température peut être obtenu par la connaissance
𝑖
des luminances centrales 𝐼𝑚 On a alors à résoudre le système constitué des
équations (e) et (g). Toutefois il est nécessaire de connaître les luminences latérales
𝑖
𝐼𝑘𝑚 , sur chaque côté k de l'élément i pour chaque direction Ω𝑚 . Des relations de
fermeture sont alors à introduire. Pour le faire, on résoud l'ETR dans le cas
monodimensionnel. La solution de I'ETR le long d'un cheminement défini par la
direction Ω est donnée par :
(h)
Où 𝒔𝒊 est la position initiale du parcours, 𝒔𝒇 est la position finale et t =𝒔𝒊 − 𝒔𝒇 . On
pose pour un élément i, et pour une direction m :
(i)
L'équation (i) devient alors :
RESOLUTION DE L’ETR-MOD
5
(j)
Dans ce travail on propose un nouveau schéma de différenciation pour résoudre les
équations de transfert radiatif en utilisant des maillages quadrangulaires. La figure 3
présente une maille quadrangulaire déformée dans le cas où 3 cotés (côtes 1, 3,4)
reçoivent le rayonnement provenant d'un seul côté (côte 2). Il s'agit de trouver les
luminances moyennes sur chacun des côtés éclairés par le rayonnement provenant
du côte 2.
Figure 3: 3 côtes reçoivent le rayonnement d'un seul côté.
❖ La luminance moyenne du côte 1 est donnée par :
RESOLUTION DE L’ETR-MOD
6
(K)
Puisque l'éclairement du côté 1 provient du côté 2, on a alors d'après l’équation (j) :
(L)
On intègre cette expression en utilisant la relation (k) pour trouva la luminance
Moyenne 𝑰𝒊𝒍𝒎 :
On obtient les relations (m) et (n) :
L'équation (m) se réduit à :
(o)
❖ Pour Ie côté 4 la luminance moyenne s’écrit :
RESOLUTION DE L’ETR-MOD
7
On obtient :
❖ Le développement pour 𝑰𝒊𝟑𝒎 est similaire et il a été rapporté à l'annexe A. Les
relations obtenues pour𝑰𝒊𝟏𝒎 , 𝑰𝒊𝟒𝒎 et 𝑰𝒊𝟑𝒎 sont appelées équations
caractéristiques. Globalement il y a trois cas possibles :
• 𝟏𝒆𝒓 𝒄𝒂𝒔 : Trois côtés reçoivent le rayonnement d'un seul côté.
• 𝟐𝒆𝒎𝒆 𝒄𝒂𝒔 : Deux côtés reçoivent le rayonnement des 2 autres côtés.
• 𝟑𝒆𝒎𝒆 𝒄𝒂𝒔 : Un côté reçoit le rayonnement des 3 autres côtés.
Nous référons à l'annexe A pour le traitement des 𝟐𝒆𝒎𝒆 𝒄𝒂𝒔 et 𝟑𝒆𝒎𝒆 𝒄𝒂𝒔 .
V. ALGORITHME DE RESOLUTION
RESOLUTION DE L’ETR-MOD
8
L'enceinte à analyser peut-être maillée avec un agencement quelconque de surfaces
quadrangulaires de façon non structurée. A cette fui, différents types de générateurs
de maillage peuvent être utilisés. Dans le cas présent, nous avons utilisé le code
COSMOS auquel nous avons adjoint un tableau de connectivité pour le repérage des
côtés des surfaces quadrangulaires.
L'algorithme développé reconnaît les cas considérés (1° , 2° ou 3° cas) selon
L’orientation des éléments par rapport à la direction de propagation donnée, et il
applique la forme appropriée de l'équation caractéristique.
En équilibre radiatif, l’algorithme de détermination du champ de température
se déroule comme suit :
a) Choix d'un premier champ de température,
b) Initialisation des luminances centrales (moyennes) et latérales de chaque
élément et introduction des conditions aux limites physiques de l'enceinte.
c) Pour une direction donnée de la quadrature, détermination des luminances
latérales moyennes de chaque élément.
d) Pour la même direction de la quadrature, détermination des luminances
centrales en procédant itérativement :
RESOLUTION DE L’ETR-MOD
9
Où n -1 représente l'itération précédente, et 𝑰𝒐𝒊(𝒏−𝟏) est calculé à partir de la
relation :
e) Vérification du critère de convergence en luminance latérale :
Sinon retour à l'étape c).
CONCLUSION
Grosso-modo la méthode des ordonnées discrètes est l’une des méthodes les plus populaires,
en raison de la précision obtenue pour un faible temps de calcul, la rendant intéressante pour traiter
des problèmes instationnaires. Toutefois, pour résoudre l’ETR, l’utilisation d’une discrétisation
angulaire associée à ce type de résolution spatiale comporte deux défauts :
–La diffusion numérique, due aux schémas de fermeture utilisés, peut engendrer des valeurs
de luminance négatives ou supérieures à celles de la solution physique ;
–L’effet de rayon, lié à la discrétisation angulaire, peut générer des discontinuités dans le flux
de chaleur conduisant à des solutions non physiques.
Merci pour votre aimable attention
RESOLUTION DE L’ETR-MOD
10
