Modèle D'équations (Réseau Thermique) 2016
Modèle D'équations (Réseau Thermique) 2016
Modèle D'équations (Réseau Thermique) 2016
Une description des équations de masse, de quantité de mouvement et d'énergie qui ont été
utilisées dans les deux études est donnée dans le tableau 2. Ces équations décrivent
l'écoulement incompressible et visqueux dans les cheminées solaires et sont dérivées pour un
volume de contrôle élémentaire dans le collecteur et dans la cheminée D'une centrale
électrique de cheminée solaire. La hauteur du toit du collecteur est progressivement soulevée
vers la tour et l'écoulement est supposé radial à une dimension entre l'entrée du collecteur et
l'entrée de la turbine. Dans la tour, un volume de commande axial arbitraire est choisi entre la
sortie de la turbine et la sortie de la tour. Bernardes énergie Eq. (6) pour le collecteur
comprend la chaleur échangée avec la terre et le toit. La température, la vitesse et par
conséquent la distribution de pression sont obtenues par la résolution des équations de masse,
de moment et d'énergie.[158]
Fig.01
Equation de quantité de mouvement
L'application de l'équation de quantité de mouvement à travers les différents sections du
collecteur est obtenu à partir de la Fig. 2:
l'équation de continuité:
La pression statique p2 est influencée par la composante connectée et sert donc de condition
de couplage pour l'équation d'impulsion pour le collecteur et la cheminée. Avec toute la force
de résistance de l'entrefer du collecteur,
Corrélations utilisées pour la contrainte de cisaillement dans le collecteur (plaque plate) [42].
1
τ =c w ρu2 Δx
2
On obtient pour le coefficient de résistance
Avec
1. 292
cw=
√ Rex
O. O 592
cw=
Au laminaire, lisse, Re x <5×10
5
√Re x
5
0 . 074 1742
cw= −
5
A la turbulence, lisse, 5×10 < Rex < 10
7
√5 Re Rex
x
5
c
A la transition, lisse, on utilise w pour un écoulement turbulent avec Re x =5×10
où τ est la contrainte de cisaillement agissant sur l'air en contact avec la surface du
collecteur.
(X)
Dans l'hypothèse où la chaleur utile, qui est transmise à l'air, le long d'une longueur de la
longueur de l'élément L est constant, la température à la sortie peut être trouvée en intégrant
l'équation ci-dessus
Réseau thermique
Pour tenir compte de l'équilibre énergétique de la couverture, il faut tenir compte du transport
de chaleur par rayonnement infrarouge de grande longueur d'onde. Les coefficients pour les
transitions thermiques suivantes doivent également être déterminés:
• entre la couverture et l'air ambiant
• entre le couvercle et l'écoulement d'air, dans le collecteur
• entre l'absorbeur et l'air en circulation
• entre les deux couvercles (dans le cas d'une double couverture)
• entre le couvercle et le réservoir d'eau.
Un film absorbant (par exemple un film noir) sur la surface de la terre sert d'absorbeur
du collecteur. Lors de l'utilisation d'un film absorbant, la surface de la terre est
considérée comme un absorbeur.
On a utilisé une procédure analogue à / Ong 1995 / pour calculer la température dans le
collecteur.
Le bilan thermique des points du réseau thermique pour un collecteur à simple couvercle avec
un stockage d'eau est déterminé comme suit:
En supposant que la température de l'air varie de façon linéaire le long de chaque section de
collecteur, la température moyenne est alors égal à la moyenne arithmétique
( T f , i−T f , o )
Tf=
2
La chaleur utile transférée au flux d'air en mouvement est écrit en termes de la moyenne du
fluide et la température a l'entrée comme:[3]
q=Γ ( T f −T f ,i ) (5-24
ou
∘
m cp
Γ=
πrL
La figure xx montre le bilan énergétique dans un volume de contrôle du fluide de
Capteur solaire plat. En tenant compte de la variation de l'énergie totale avec le temps et de la
chaleur totale transférée dans le volume de contrôle du fluide, le bilan énergétique sous les
propriétés transitoires du fluide de travail peut s'écrire:
figure xx
collecteur a divisé à 5 parts égales de distance. On pourrait supposer que la terre est un corps
semi-infinitif qui peut stocker l'énergie. Dans chaque section avec l'existence de la vitesse et
le nombre de Reynolds, le coefficient de transfert de chaleur par convection peut être calculé.
La figure 5 montre le schéma du modèle de résistance thermique pour le toit simple vitrage du
du sol. Les équations suivantes du bilan thermique sont obtenues à partir du réseau thermique
aux points en considérant la résistance de contact thermique pour le modèle à simple vitrage:
[47]
Tf1 :
h1 ( T 1 −T f 1 ) +h2 ( T 2−T f 1 )=q}air} } { ¿¿ ¿
T 1 h1 −T f 1 ( h1 +h2 ) +h2 T 2=q}air} } { ¿¿ ¿ (2)
S2 +h32 ( T 3 −T 2 ) +U w ( T 3 −T 2 ) =h2 ( T 2 −T f 1 ) +h r 21 ( T 2 −T 1 ) +h3 ( T 2 −T f 2 )
T2:
−T 1 hr 21−T f 1 h2 +T 2 ( h2 +hr 21+h3 +h 32+U w ) −T f 2 h3 −T 3 ( h 32+U w ) =S2
(3)
h3 ( T 2−T f 2 ) +h 4 ( T 3 −T f 2 ) =0
Tf2
h3 T 2 −T f 2 ( h 3 +h 4 ) +T 3 h4 =0
(4)
h 4 ( T 3 −T f 2 ) +U w ( T 3 −T 2 ) +h32 ( T 3 −T 2 ) +h 5 ( T 3−T f 3 ) =S 3
T3
−T 2 ( U w +h 32) −h 4 T f 2 −T 3 ( h4 +U w +h 32+h 5 ) +T f 3 h5 =S 3
(5)
h5 ( T 3 −T f 3 ) +h6 ( T 4 −T f 3 ) =0
Tf3
T 3 h5 −T f 3 ( h5 +h 6 ) +T 4 h6 =0
(6)
h6 ( T 4 −T f 3 ) +U b ( T 4 −T 4 , 0 ) =0
T4
−T f 3 h 6 +T 4 ( h6 +U b )=U b T 4 , 0
(7)
[ ]
0 0 0 0
0 0 0 0 ¿ ¿ ¿ 0 ¿ 0 ¿ 0 ¿ 0 ¿ ¿ ¿ ¿ ¿ 0 ¿ 0 ¿ 0 ¿ 0 ¿ 0 ¿ ¿ ¿ ¿ 0 ¿ 0 ¿ 0 ¿ 0 ¿ 0 ¿ ¿ ¿ [ ¿ ]=¿
0 ¿ 0 ¿
est le réseau thermique du collecteur avec une simple couverture sans réservoir d'eau
représenté sur la fig 1
Fig. 1
S1 + hr ( T 3 −T 1 ) =hrs ( T 1−T s ) +h w ( T 1−T ∞ ) +h 1 ( T 1−T f 1)
T1 : 21
Tf1 :
h1 ( T 1 −T f 1 ) +h2 ( T 3 −T f 1 )=q}air} } { ¿¿ ¿
T 1 h1 −T f 1 ( h1 +h2 ) +h2 T 3=q}air} } { ¿¿ ¿ (b)
h2 ( T 3 −T f 1 ) +hr 21 ( T 3 −T 1 ) +h5 ( T 3 −T f 3 ) =S 2
T3
−T 1 hr 21−h2 T f 1 −T 3 ( h 2 +hr 21 +h5 )−T f 3 h5 =S 2
(c)
h5 ( T 3 −T f 3 ) +h6 ( T 4 −T f 3 ) =0
Tf3
T 3 h5 −T f 3 ( h5 +h 6 ) +T 4 h6 =0
(d)
h6 ( T 4 −T f 3 ) +U b ( T 4 −T 4 , 0 ) =0
T4
−T f 3 h 6 +T 4 ( h6 +U b )=U b T 4 , 0
(e)
[ ]
0 0
0 0 ¿ ¿ ¿ 0 ¿ 0 ¿ ¿ ¿ ¿ [ ¿ ] =[ ¿ ]
¿0 ¿
En général, les matrices ci-dessus peuvent être écrites comme
[ A ] [ T ] =[ B ]
Le vecteur de température moyenne peut être déterminé par la matrice inversée
[ T ] =[ A ]−1 [ B ]
Coefficients de transfert de chaleur
Le coefficient global de perte de chaleur maximale peut être obtenu par:
U t =( hw +hrs )
Avec
k
h w = Nu
L
et
σε ( T 1 +T s ) ( T 21 +T 2s ) ( T 1 −T s )
h rs=
( T 1 −T ∞ )
Référence [153,154]
La température du ciel propre Ts obtenue par Berdahl et Martin (1984) est donné par
1
[
T s =T ∞ 0 .711+ 0. 0056 ( T ¿ .15 )+ 0. 000073 ( T dp−273. 15 ) 2 +0 . 013 cos ( 15 t ) ] 4
Où t est l'heure de minuit. Le coefficient de transfert de chaleur au sol est donné par
2b
Ub
√ πt
Avec
b=√ kρc p
(
1 1
+ −1
ε1 ε2 ) Référence [153,154]
σ ( T 22 + T 23 ) ( T 2 +T 3 )
h r32=
( 1 1
+ −1
ε2 ε3 )
Les flux de chaleur du rayonnement solaire absorbés par les surfaces sont:
S1 =α 1 H
S2 =τ 1 α 2 H
S3 =τ 2 α 3 H
S 4 =τ 3 α 4 H
Où S1, S2, S3 représente le rayonnement solaire absorbé par le second couvercle, par le
premier couvercle, par le film plastique transparent, par l'absorbeur, respectivement
La transmission et la capacité d'absorption d'une seule
τ ≃τ a τ r (33)
α≃1−τ a (34)
Pour résoudre le problème de l'état de chaleur non stationnaire dans le système de stockage
d'eau, nous considérons un problème de valeur limite homogène pour une plaque plane
infiniment large avec une température limite prescrite. Les distributions de la température et
du flux de chaleur sont déterminées analytiquement. Le coefficient de transfert de chaleur du
système de stockage d'eau peut être calculé par:
[ (
]
)
t1
−δ2k α
q01 1 ∞ 2 Lw sin ( δ k )2 L2w
U w=
Δϑ
= ∑
t k=0 αδ k [ δ k +sin ( δ k ) cos ( δ k ) ]
e
Avec:
δ 1=π /2 , δ 2 =3 π /2 , δ 3 =5 π /2 ,..., δ k =( k −1/2 ) π . [03]
[ T ] =[ A ]−1 [ B ] Le vecteur de température peut être résolu par inversion matricielle. Cependant,
il peut être résolu dans la présente étude par le programme MATLAB. Pour trouver les
corrélations h1, h2 et Nu utilisées pour la convection forcée (plaque plate, température
constante) [03]:
1 pr
Nu x=
√π
√ Rex⋅ 1
(1+ 1. 7 pr )
1
6
4
+21 . 36 pr and
Num , lam=2 Nux
5
Au laminaire Re x <5×10
0 . 037 Re pr
x0 . 8
Num , tur =
( pr )
2
3
1+2. 44 Re −0 . 1 −1
x
5 7
A la turbulence, 5×10 < Rex < 10 , 0 . 6< pr <2000
Num= Nu √ m , lam2
+Nu 2
m , tur
Les corrélations utilisées pour la convection naturelle (plaque plate, température constante)
[12] sont:
1
4
Num=0 . 54 Ra
10 4 < Ra<107 Surface horizontale supérieure chauffée
1
3
Num=0 . 14 Ra
107 < Ra<1011 Surface horizontale supérieure chauffée
1
4
Num=0 . 27 Ra
105 < Ra< 1010 Surface horizontale chauffée inférieure
Cheminée
La cheminée convertit l'énergie thermique produite par le capteur solaire en énergie cinétique.
La différence de densité créée par l'élévation de température dans le collecteur fonctionne
comme force motrice. Le transfert de chaleur qui a lieu à travers la surface de la section de
cheminée est négligeable. Appliquer l'équation de quantité de mouvement à travers une
section différentielle des rendements de la cheminée
√ 2
[∫ ]
Hk
w=
ρ 0 ( ρ0−ρ ) gdz−Δpfrottement
En dehors de la cheminée, la variation de pression, de température et de densité de l'air est
calculée compte tenu de l'atmosphère standard
(
T ∞ ( z ) =T ∞ ( 0 ) 1−
k−1 z
k H0 )
( )
k / ( k−1 )
k −1 z
p∞ ( z )= p∞ ( 0 ) 1+
k H0
( )
1/ ( k−1 )
k−1 z
ρ∞ ( z )=ρ∞ ( 0 ) 1+
k H0
Avec
Rl T∞ (0 )
H 0=
g
etk =1. 235 (standard atmosphère)
La variation de température, de pression et de densité de l'air à l'intérieur de la cheminée est
calculée en tenant compte d'un processus d'expansion adiabatique. En remplaçant k = 1.235
par k = 1.4005 et l'indice ∞ par t, on peut les estimer à partir des équations (26) à tout Z et t.
La vitesse à n'importe quelle section à travers la cheminée est calculée en appliquant
l'équation de continuité dans l'état stationnaire. De plus, les corrélations employées pour le
frottement dans la cheminée (en supposant un tuyau lisse) pour le flux laminaire et turbulent
sont considérées comme éditées dans. [157 pour mes réf].
Turbine et générateur
Le flux de chaleur produit par le collecteur est converti en énergie cinétique (courant de
convection) et en énergie potentielle (chute de pression à la turbine) à travers la cheminée.
Ainsi, la différence de densité de l'air provoquée par l'élévation de température dans le
collecteur fonctionne comme une force motrice. La colonne d'air plus légère dans la cheminée
est reliée à l'atmosphère environnante à la base (à l'intérieur du collecteur) et au sommet de la
cheminée, et elle acquiert alors un ascenseur. Entre la base de la cheminée (sortie du
collecteur) et les environs, une différence de pression
Δp tot est produite. La chute de pression
à travers la turbine peut être exprimée en fonction de la différence de pression totale
1
Δp turb= Δp tot − ρw 2
2
Avec
Hk
Δp tot =∫0 ( ρ0 −ρ ) gdz
La vitesse à la sortie de la cheminée peut être trouvée en utilisant
w=wtot √1−x
Où x est le facteur de chute de pression à la turbine et
w tot est la vitesse obtenue en
négligeant les pertes par frottement.
La puissance théorique utilisable absorbée par la turbine est
P= Δp tot Awtot ηt x √ 1−x
Dans Réf. [100], la chute de pression de la turbine à la proportion optimale de la différence de
pression totale a été calculée à 2/3 (soit 12,5% de plus que pour une éolienne autonome, cette
( )
Δp tur
=16/ 27
Δptot
proportion étant donnée par le facteur de Betz Comme Ainsi, cette hypothèse
de perte de charge constante a conduit à surestimer la taille des passages d'écoulement dans
l'usine et quelque part la puissance de sortie. Selon Réf. [10], le facteur maximal de chute de
( )
Δp tur
Δptot max
pression est calculé à chaque pas de temps comme suit:
( )
Δp tur n−m
=
Δptot max n+1
Considérant que n qui est la puissance du débit volumique dans la cheminée qui fait la
pression de perte due au frottement, est proportionnelle à elle, qui a choisi comme 1,75 selon
la relation du débit volumique avec la perte de pression au régime de turbulence. De plus, m
est la puissance négative du débit volumique qui rend la pression potentielle proportionnelle à
celle-ci, calculée par
η coll
m=−
α absor
La puissance maximale disponible peut être calculée par l'équation suivante.
( )
Δp tur
Pth max = ⋅Δp tot⋅V⋅A
Δp tot max
Et la puissance de sortie obtenue selon le rendement de la turbine par
Pout =ηtur⋅Pth,max
Le modèle théorique suppose que pour un collecteur court, les températures des «limites»
entourant les courants d'air sont uniformes et les températures des courants d'air varient
linéairement le long du collecteur. Un collecteur long peut être supposé être divisé également
en un nombre fini de petits collecteurs ou sections. Les températures de paroi et d'air moyen
de la première section sont égales à la température ambiante. Les coefficients de transfert de
chaleur sont évalués en fonction des valeurs initialement supposées. Un processus itératif est
ensuite créé et les températures moyennes pour la section sont calculées en utilisant les
équations dérivées en employant une inversion de matrice de matrice standard. La
température de l'air d'entrée de la section suivante est égale à la température de sortie de la
dernière section. Le processus itératif est répété jusqu'à ce que toutes les températures
moyennes consécutives diffèrent par moins d'une valeur souhaitée.
http://fvirtman.free.fr/recueil/index.html(codes en C++)
L'analyse technique
Les cotes d'énergie pour les matériaux de couverture, qui ne sont pas dans la gamme
infrarouge longue onde, sont opaques, nécessitent une étude du transfert de chaleur par
rayonnement infrarouge à ondes longues. En annexe, le modèle utilisé pour déterminer les
propriétés optiques des films selon / Whillier, A. 1977 /.
Le débit d'air dans le système est le résultat de la différence de densité entre l'air ambiant et
l'air plus chaud à l'intérieur. Il peut donc être appelé flux de convection libre. En même
temps, cependant, ce flux présente également des caractéristiques d'un écoulement forcé.
Pour cette raison, un transfert de chaleur doit être assumé par la convection libre et forcée.
Les coefficients de transfert de chaleur h nécessaires pour le transfert de chaleur des feuilles
vers l'air s'écoulant sont déterminés en fonction des numéros d'écrou Nu et du diamètre
hydraulique D h
La conductivité thermique k dépendante de la température k de l'air peut être trouvée dans les
tableaux pertinents (par exemple / Wärmeatlas 1997 /) et le diamètre hydraulique est connu,
ainsi les nombres d'éléments nutritifs restent à déterminer. Selon Churchill 1977, le nombre
d'écrous en convection mixte peut être déterminé à partir du nombre d'écrous en convection
forcée et du nombre d'écrous en convection libre:
La valeur 3 est recommandée pour l'exposant. Pour déterminer le nombre d'écrous en
convection forcée, l'état d'écoulement (laminaire-turbulent) dans le collecteur doit d'abord
être estimé en utilisant le nombre de Reynolds. Selon Merker 1987 /, des observations
expérimentales montrent que la couche limite de la plaque est complètement turbulente
seulement d'environ Rex = 3 à 5 x 105. Après avoir défini le nombre de Reynolds local pour
le débit de colonne
Les vitesses d'air correspondantes dans la zone de passage de l'état d'écoulement laminaire à
turbulent peuvent être estimées en connaissant la viscosité cinématique de l'air en circulation.
Le diamètre hydraulique peut être déterminé comme suit:
Zone de passage
Zone humide
La vitesse est
A la frontière de l'écoulement laminaire (Re = 2300) est obtenu à partir de (5-91) avec Dh =
4 m de longueur d'entrée hydrodynamique de Lhyd = 101 m.
/ Na et Lu 1973 ont réalisé une analyse intégrale du problème hydrodynamique de la surface
de passage pour les plaques planaires lisses et ont constaté que la longueur d'entrée
hydrodynamique, en supposant que les couches limites hydrodynamiques se développent à
partir du début des deux plaques est égale à L hyd = 13,75Dh à Re = 2,21 x 105. Les
résultats expérimentaux de / Byrne et al. 1969-70 / sur le paramètre hydrodynamique de
l'impulsion et de l'épaisseur de déplacement sont en accord avec les prédictions de / Na et Lu
1973 /. Selon cette étude, une longueur d'entrée hydrodynamique de L hyd = 55 m pour Re =
2,21 × 105 est obtenue avec Dh = 4 m.
Dans le flux le long d'une plaque plane (fig. 5-9) formant le nombre de Reynolds Rex avec la
longueur de la plaque x, ie. De la plaque commence à un certain PlatLtenlänge ou
certain nombre de Reynolds local, le débit est complètement laminaire. Après passage dans
un soi-disant. La zone de transition, le débit est alors complètement turbulent.
Efficacité du collecteur;
L'annexe contient l'équation détaillée pour déterminer les propriétés physiques de l'air
humide
calcul de simulation:
Afin de pouvoir utiliser le modèle physique développé pour calculer les variables pertinentes,
un programme de calcul a été créé. Un cycle de programme est de fournir les valeurs de la
température, l'humidité relative et le débit massique de l'air dans le système pour une journée
entière. Le langage de programmation utilisé Delphi est généralement bien adapté à la
formulation des sciences technico-scientifiques. Les équations complexes ne permettent pas
une solution fermée de sorte que les méthodes itératives de solution doivent être utilisées. Afin
d'utiliser la capacité de l'ordinateur à effectuer un grand nombre d'opérations informatiques
similaires en un temps très court, l'heure du jour, la longueur du collecteur et la hauteur de la
cheminée ont été divisées respectivement en pas de temps et de longueur. Cette approche
discontinue aboutit à un modèle quasi-stationnaire, tel que décrit dans la section 5.1.2, avec
seulement les paramètres de calcul changeant dans le grand nombre d'opérations identiques.
Le nombre d'étapes de temps et de longueur ne doit pas être trop grand car cela augmente le
temps de calcul. Cependant, comme pour la modélisation des processus thermodynamiques
dans la plante, les temps doivent être suffisamment petits pour que les variations temporelles
des variables pendant un intervalle de temps puissent être négligées. Le raisonnement suivant
est utilisé pour justifier le choix de la longueur des étapes de temps. Le point de départ est
considéré comme étant les variables qui affectent les variations temporelles des variables
d'état pertinentes de l'air. Ce sont les données environnementales décrivant les conditions
météorologiques, ainsi que la température de la surface du sol sous le système. On suppose
que dans une période où les variations temporelles des variables environnementales peuvent
être négligées, la modification résultante de l'état de l'air dans l'installation peut également
être négligée. Cela signifie que le choix de la longueur d'un pas de temps est basé sur
l'enregistrement sensible de la météo. L'enregistrement des données météorologiques en
temps de 600 secondes, comme dans le cas des plantes de Manzanares (Haaf et al., 1986 /),
est courant. Selon ce raisonnement, le modèle de calcul a été utilisé pour diviser un jour en
144 périodes de temps de 600 secondes chacune.
la structure du programme:
Cette section fournit des informations sur la structure du programme de calcul (voir figure 5-
11) et un bref aperçu de la fonction des différents blocs de programme.
Données
L'utilisateur doit saisir les valeurs des paramètres non dépendants du temps. Votre entrée est
contrôlée par menu avec des valeurs par défaut et l'option de modification. Les valeurs des
paramètres suivants doivent être saisies:
paramètres de collection
- Mesures dimensionnelles
• diamètre [m]
• hauteur de la canopée à l'entrée [m]
• Hauteur de Canopy à la sortie [m]
- Les propriétés physiques des couvertures
• transmittance
• réflectance [-]
• émissivité [-]
• coefficient d'extinction [1 / m]
• Épaisseur [m]
• Paramètres de Cheminée
- Mesures dimensionnelles
• le diamètre à l'entrée (inférieure) [m]
• diamètre à la sortie (en haut) [m]
• Hauteur [m]
• épaisseur de la paroi [m]
- Propriétés thermiques
• conductivité thermique de la paroi [W / (m · K)]
• paramètres d'absorption
- Propriétés optiques
• Coefficient d'absorption [-]
• émissivité [-]
- mesures dimensionnelles
• Profondeur [m] (pour visualiser la température au fond)
• épaisseur du réservoir d'eau [m] (si disponible)
• l'absorbance, l'émissivité et la transmissivité de la partie supérieure du réservoir d'eau [-]
• distance entre l'absorbeur et le fond [m] (pour le contact résistance thermique)
- Propriétés thermiques
• La conductivité thermique du sol [W / (m · K)]
• la conductivité thermique du sol [m / s]
• ou le coefficient de pénétration de la chaleur (effusivité thermique) du sol [W s1 / 2 / km²]
• environnement
- Propriétés thermiques et hydrodynamiques
• Vitesse du vent [m / s]
• Température ambiante [K]
• température du point de rosée [K]
• pression ambiante [Pa]
- Propriétés optiques
• visibilité horizontale [m]
• L'épaisseur optique d'extinction à travers les nuages cirrus [m]
• emplacement
• Latitude [degrés]
• longitude [degrés]
• Date
Jour • Jour Nombre de la première à simuler
• Nombre de jours à simuler
• Temps de longitude zone [°]
• albédo de surface [-]
• paramètres arithmétiques
• Erreur maximale [-]
• Nombre d'étapes de longueur [-]
• incrément de temps [s]
• turbine
• facteur de prise de pression [-]
• l'efficacité de l'ensemble du générateur de turbo [-]
Irradiation Subroutine:
Le sous-programme calcule le rayonnement solaire diffus et direct selon les paramètres (A-1)
à (A-19) en fonction des paramètres d'emplacement
optique subroutine:
Pour couvrir les facteurs d'absorption et de transmission en diffus et irradiation directe par
les équations (5-140) sont déterminées à (5-146).
Propriétés subroutine
Système de solution :
Le sous-programme comprend la boucle de calcul sur la longueur de la cheminée et le
collecteur. Les paramètres thermiques et hydrodynamiques présentés dans le tableau 5-2 et le
tableau 5-1 sont déterminés ici.
Présentation des résultats:
Ce programme offre les options suivantes pour afficher les résultats:
• Coefficients de température et de transfert de chaleur de l'air, de la couverture et du
plancher, des vitesses et de l'humidité dans chaque élément de longueur du collecteur à un
moment donné.
• Coefficients de température et de transfert de chaleur de l'air, des vitesses et de l'humidité à
chaque élément de longueur du foyer à un moment donné.
• Débit massique, irradiation, puissance, efficacité, degré de transmission, degré
d'absorption, température ambiante en fonction de l'altitude sur la journée.
Les résultats peuvent être écrits dans un fichier. De plus, des données sont générées pour la
sortie graphique, les représentations suivantes sont possibles:
• Cours de la température de l'air, du sol et de la couverture, de la vitesse et des courants de
chaleur sur la longueur du capteur à tout moment
• La température de l'air, la vitesse et le flux de chaleur sur la longueur de la cheminée à tout
moment
Validation du modèle:
Dans cette section, les résultats des mesures expérimentales de Manzanares et les post-
calculs du modèle appliqué sont présentés et comparés. Le comportement énergétique des
plantes examinées doit être présenté sur la base d'un jour représentatif. À cette fin, les jours 7
et 8 juin 1987 ont été sélectionnés. Les paramètres suivants doivent être considérés comme
des données essentielles:
• Densité de flux de radiation globale
• Température de l'air selon l'altitude
• Humidité relative
• Vitesse du vent
• Facteur d'échantillonnage de pression
Avant le début du calcul de simulation, les valeurs des paramètres constants, qui sont listées
dans le Tableau 5-3, ont été définies.
Valeurs d'entrée pour le programme
valeur Unité
Hauteur de la voilure à l'entrée du collecteur 2m
Hauteur de la voilure à la sortie du collecteur 6m
Diamètre du collecteur: 250 m
Couverture facile Oui -
Transmittance constante Oui,0,90 -
Diamètre de la came à l'entrée 10 m
Diamètre de came à la sortie 10 m
Hauteur du foyer 200 m
Perte de chaleur dans le foyer No -
Absorbance de l'absorbeur 0,93 -
Niveau d'émission de l'absorbeur 0,90 -
Hauteur de rugosité du sol 0,05 m
Conductivité thermique du sol 0,6 W / m · K
Conductivité en température du sol 2,91 x 10-7 m² / s
Location Manzanares -
Latitude 39,03 (nord) °
Longitude: 3,14 (Est) °
Date 07 et 08.06.1987 -
Erreur maximale à chaque étape de temps 0.1%
Nombre de longueurs 200 -
Étapes de temps 600 s
Efficacité de l'ensemble turbogénérateur 0,42
Figue. 5-12 et Fig. 5-13 montrent le cours de la densité de flux de rayonnement global et la
vitesse du vent externe en fonction de l'heure du jour les 7 et 8 juin 1987. Les immersions
significatives dans le cours du rayonnement sont dues à l'effet d'ombrage des transitoires des
nuages. Dans le cas d'un ombrage du dispositif de mesure de rayonnement, l'ensemble du
collecteur est considéré comme ombragé; L'ombrage partiel ne peut donc pas être considéré.
Selon ces chiffres, seuls des vents légers ont eu lieu le jour en question. Figue. 5-14 et Fig. 5-
15 représentent l'évolution de l'humidité relative de l'air extérieur et de la température les 7
et 8 juin 1987 respectivement. Figue. 5-16 et Fig. 5-17 montrent le cours de la température
dans le Collecteur les 7 et 8 juin 1987 respectivement
Référence [publik1]
Continuity:
(1)
(2)
Energy Equation:
(3)
(4)
(5)
The model constants for the turbulence model can be found in the literature Error: Reference
source not found.