These Mo Rille
These Mo Rille
présentée par :
MORILLE Benjamin
Composition du jury :
Je remercie l’ensemble des personnes qui m’ont permis d’en être à écrire ces quelques
mots aujourd’hui :
Je remercie mes parents qui ont été présents pour me soutenir et m’accompagner tout
au long des étapes importantes de ma vie, tout au long de mes études.
Je remercie celle qui partage ma vie et qui sait être, de jour en jour, de plus en plus
précieuse. Merci pour ton aide, ton soutien, ta présence au quotidien.
Je remercie mes compagnons de covoiturage qui ont égayé mes trajets pendant mes
deux années de classe préparatoire et ont eux aussi contribué à leur manière à mon parcours.
François, merci pour ces discussions toujours aussi posées et intéressantes et cette amitié
sincère qui dure malgré la distance. Julien, merci pour ton humour et ta répartie dont je ne sais
malheureusement m’inspirer pour écrire ces quelques mots.
Je remercie mes amis basketteurs qui me permettent de m’épanouir dans ma vie
extraprofessionnelle, gage de stabilité. Merci à mes coéquipiers qui savent supporter mon
caractère sur le terrain, merci aux filles qui supportent mes consignes depuis maintenant trois
saisons.
Je remercie l’ensemble de l’équipe du laboratoire, pour la bonne humeur qui y règne et
pour les bons mets qui ont émoustillé nos papilles durant ces trois années.
Je remercie l’ADEME et la région Pays de la Loire pour la confiance qu’ils ont placée
en moi, pour le soutien qui m’a permis de mener ces travaux de recherche.
Je remercie Cédric Garnier, ingénieur ADEME, qui a suivi et encadré mes travaux ;
Thierry Boulard pour sa grande expérience dont il a su me faire profiter dès qu’il y en a eu
besoin ; et l’ensemble des personnes qui ont apporté un avis critique sur mes travaux tout au
long de l’avancement de ma thèse.
Je remercie Christophe Migeon, pour son encadrement, ses remarques avisées, et sa
disponibilité remarquable durant ces trois années.
Je remercie enfin Pierre Emmanuel Bournet, mon directeur de thèse, qui a toujours été
présent pour répondre à mes questions, orienter et coordonner mes travaux, et sans qui cette
thèse ne serait pas ce qu’elle est aujourd’hui.
3
Résumé
L’objectif de cette thèse est de développer un modèle numérique instationnaire
permettant de simuler l’évolution du climat distribué sous serre pour des conditions de nuit en
présence de chauffage et pour des conditions diurnes où l’influence du rayonnement solaire
est prépondérante. La thèse s’appuie sur des approches expérimentales et numériques qui ont
pour ambition d’être appliquées à deux situations, hivernales et pré-estivales correspondant
aux problématiques auxquelles les producteurs doivent faire face : gestion de la température et
de l'humidité avec une préoccupation de réduction des consommations énergétiques.
Des campagnes de mesures ont été menées au printemps 2010 et à l’hiver 2011 pour,
d'une part, améliorer la compréhension des interactions entre les paramètres climatiques, et
d'autre part, fournir des données d’entrée et de validation aux modèles numériques.
Les campagnes de mesures de nuit (hiver 2011) en présence de chauffage ont mis en
évidence les besoins plus importants en chauffage et les températures de toiture plus basses
lors des nuits claires. Les campagnes de mesure menées au printemps 2010 ont permis mettre
en évidence l’influence prépondérante du rayonnement solaire sur la température et la
transpiration du végétal ainsi que le phénomène de stockage-restitution de chaleur dans le sol.
Enfin, la forte sensibilité du taux d’humidité ainsi que la moindre dépendance de la
température vis-à-vis de la « qualité » de la ventilation a été démontrée. Les données
expérimentales sont exploitées afin de mettre en œuvre les modèles de transpiration adaptés
aux conditions de culture sous serre. Les bilans hydriques réalisés sur le végétal divisé en
deux couches fournissent des informations précieuses sur la distribution de la transpiration et
des contributions de chaque couche.
Les simulations numériques sont réalisées en 2D instationnaire, en application sur une
situation nocturne (nuit claire et nuit nuageuse) et une situation diurne (journée parfaitement
ensoleillée). Elles intègrent deux éléments novateurs : la résolution de l'ETR dans le couvert
végétal et l'utilisation d'une routine permettant de simuler la condensation. L'évolution des
paramètres climatiques simulés est validée sur la base des données expérimentales. Il apparait
cependant que les écarts de température et d'humidité entre l'air dans et au dessus du végétal
sont sous estimés. Enfin, les modèles numériques sont exploités pour réaliser des bilans
hydriques et énergétiques.
Mots clés : Climat distribué, expérimental, numérique, CFD, instationnaire, condensation,
modèle de culture, équation de transfert radiatif, bilans hydriques, bilans énergétiques
4
Abstract
The objectives of this thesis are to develop a numerical tool which would be able to provide
exhaustive information on the evolution of the distributed climate in a greenhouse. The thesis
is based on experimental and numerical approaches and aims for giving information on the
main concerns of the growers: the temperature and the humidity management with a special
attention giving to the reduction of the energy consumptions. This is the reason why the study
considers two situations in which these main concerns are the most pronounced.
Experimental investigations on the climatic parameters have been led at the end of
spring 2010 and at winter 2011 in order to increase the comprehension of the interaction
between the physical mechanisms governing the climate and to provide data carry on the
numerical part of this study.
The higher heating needs and lower roof wall temperature during clear nights were
evidence as well as most important role played by the solar radiation during days especially
on temperature, the crop transpiration and the phenomenon of storage and restitution of heat
in the ground. Then, the sensible ventilation efficiency impact on humidity and the lower on
temperature are evidenced. Last, experimental data were also exploited through the
application of transpiration models adapted to greenhouse conditions. Water balance on a crop
divided into two layers provides supply information on the transpiration distribution.
The unsteady 2D simulations consider night (clear, cloudy) and daytime (under clear
sky conditions) situations. Innovative improvements are implemented into existing models :
the RTE is solved in the canopy volume and provide the exact absorbed radiation in each crop
layer all along the day and a condensation sub-routine is activated to simulate humidity rate
evolution during the night. Models are validated against a large number of experimental
climatic parameters. Simulations show good agreements with experimental data even if slight
underestimations of humidity and temperature differences between above and inside crop air
evidence the need of improvements on this point. Last, water and energy balances are
obtained thanks to these models.
5
Nomenclature
Nomenclature
Lettres latines
I : Luminance cd.m-2
Jm : flux de diffusion de l’espèce m kg.m-2.s-1
K : coefficient d’extinction -
K : énergie cinétique turbulente m2.s-3
K0 : perméabilité du milieu poreux m2
L : longueur du chemin optique m
7
Nomenclature
& ′′′
m : Variation de la masse d’air humide kg.m-3.s-1
& ′h′′2o
m : Flux massique de vapeur d’eau kg.m-3.s-1
mm : masse de l’espèce m kg
Nu : nombre de Nusselt -
Nφ : discrétisation de l’angle solide φ -
Nθ : discrétisation de l’angle solide θ -
p : Pression Pa
Patm : pression atmosphérique Pa
2 -3
Pr : production d’énergie cinétique turbulente m .s
Pr : nombre de Prandtl -
Prt : nombre de Prandtl turbulent -
Pva : pression vapeur de l’air Pa
Pvf : pression vapeur de la feuille Pa
Pvsat : pression vapeur saturante Pa
Pvsat-air : pression vapeur saturante de l’air Pa
Pvsat-feuille : pression vapeur saturante de l’air à la température de la feuille Pa
q• : puissance volumique thermique échangée W
Q : chaleur échangée J
Q• : puissance thermique échangée W
Ql : flux de chaleur latente W.m-3
Qs : flux de chaleur sensible W.m-3
R0 : rayonnement à la surface du végétal W.m-2
ra : résistance aérodynamique s.m-1
Re : nombre de Reynolds -
RGLO : rayonnement grande longueur d’onde W.m-2
Rn : rayonnement net W.m-2
rs : résistance stomatique s.m-1
rs,min : résistance stomatique minimale s.m-1
Sct : nombre de Schmidt -
Sh : terme source de l’équation de l’énergie W.m-3
8
Nomenclature
Caractères grecs
α0 : amplitude de la sollicitation du sol K
β : coefficient de dilatation à pression constante K-1
γ : constante psychrométrique Pa.K-1
Γs : rayonnement absorbé par le sol W.m-2
∆ : gradient de pression de vapeur saturante entre l’air et la feuille Pa.K-1
9
Nomenclature
Vecteurs et tenseurs
r
δFext : forces extérieures infinitésimales exercées sur le fluide N.m-3
r
δFsurfacique : forces surfaciques infinitésimales exercées sur le volume de fluide N.m-2
r
δFvolumique : forces volumiques infinitésimales exercées sur le volume de fluide N.m-3
r
q : vecteur flux radiatif W.m-2
τ : tenseur des contraintes visqueuses Kg.m-1.s-2
10
Table des matières
11
Table des matières
12
Table des matières
13
Table des matières
14
Table des matières
15
Chapitre 1 : Introduction
Chapitre 1. Introduction
I. Contexte
1. Contexte de la thèse
La thèse s’est déroulée au sein de l’unité de recherche EPHor (Environnement Physique de la
plante Horticole) d'Agrocampus Ouest à Angers. Elle s’inscrit dans la continuité des travaux
réalisés par cette équipe et notamment de deux précédentes thèses : celle de Sid-Ali Ould
Khaoua portant sur la modélisation de l'aération naturelle et du microclimat des serres en
verre et celle d'Abderzak Kichah concernant plus particulièrement le climat au niveau du
végétal avec toujours comme objectif la modélisation du climat sous serre.
La présente thèse a été cofinancée par la région Pays de la Loire et par l'ADEME (Agence de
l’Environnement et de la Maitrise de l’Energie) et s'inscrit dans le programme prioritaire de
l'ADEME relatif aux "bâtiments économes". L'objectif principal est d'améliorer les
connaissances liées au fonctionnement climatique des serres et ce dans l'optique de proposer
des solutions permettant de minimiser les consommations énergétiques et hydriques.
2. Contexte général
Les serres de production sont des abris destinés à la production de plantes ornementales,
légumières ou fruitières et sont utilisées pour soustraire ces cultures aux contraintes
climatiques extérieures afin de gérer au mieux les besoins des plantes, accélérer leur
croissance et étendre les périodes de production à l'ensemble de l’année. Les serres
constituent donc un enjeu économique majeur pour les horticulteurs et pour les maraîchers.
En France, selon une étude menée par l’ADEME (2007), les surfaces destinées à la
production sous serre représentent 2900ha de cultures maraîchères (tomate et concombre
principalement) et 2500ha de cultures ornementales (fleurs coupées et plantes en pots) dont au
total 2600ha sont des serres chauffées consommant 6TWh par an. La consommation
énergétique annuelle moyenne des serres destinées au maraîchage est de 320kWh/m² et pèse
pour 22% dans les charges financières des exploitants. En ce qui concerne le domaine
ornemental, la consommation énergétique moyenne est de 160kWh/m² représentant 11% des
charges des exploitations. L'énergie utilisée pour le chauffage des serres est essentiellement
fossile, énergie dont le coût n'a cessé d'augmenter lors des dernières années. La problématique
de la maîtrise des consommations énergétiques est donc à présent une des priorités de la
17
Chapitre 1 : Introduction
profession, l'objectif étant de continuer de procurer aux plantes le climat dont elles ont besoin
tout en réduisant les dépenses énergétiques.
La compréhension et la maîtrise de l'ensemble des mécanismes physiques et
écophysiologiques qui gouvernent le climat sous serre a été et est l'objet de nombreuses
études afin de permettre aux concepteurs de serres de proposer des solutions adaptées aux
besoins des producteurs. Alors que les études étaient, il y a encore peu de temps,
essentiellement réalisés dans le but de réduire les coûts d'investissements (prix des serres),
elles sont à présent principalement menées dans l'optique de diminuer les coûts de
fonctionnement (coûts énergétiques, consommation en eau, autres intrants...). Outre les
serristes et les exploitants, les consommateurs sont eux aussi concernés par les enjeux liés à la
maîtrise de l'ensemble des phénomènes régissant le climat des serres et ce pour trois raisons :
• La première est liée à leur envie de consommer des produits en dehors des saisons
naturelles de production.
• La seconde est relative au fait qu'ils représentent le dernier maillon de la chaîne de
production sous serre, et devraient donc également bénéficier de la réduction générale
des coûts de production.
• Enfin, la sensibilité grandissante du consommateur à la qualité des produits constitue
la troisième raison.
Ces enjeux rejoignent les principaux axes du développement durable qui reconsidère la
croissance économique en prenant en compte les aspects environnementaux via une
production qualitative plutôt que quantitative.
II. Problématique
La serre, en plus de protéger les plantes des conditions climatiques extérieures, doit permettre
de créer un climat favorable au développement des cultures. Etre capable de contrôler et
d'optimiser ce climat est le principal souhait des exploitants.
1. Gestion de la température
La capacité d'une serre à créer un climat favorable aux cultures s'explique à la fois par le fait
qu'elle constitue un espace abrité du vent mais aussi et surtout par le fait que le volume de la
serre est réchauffé par ce qu'on appelle justement l'effet de serre. En effet, le moteur principal
du fonctionnement d'une serre est le rayonnement solaire qui, "emprisonné" par les parois de
18
Chapitre 1 : Introduction
2. Gestion de l'humidité
Outre la régulation de la température, la gestion de l'humidité est aussi un enjeu important de
la climatique des serres. Toute culture rejette, dans une proportion plus ou moins grande selon
le type de culture, de l'eau dans l'air sous forme de vapeur : c'est la transpiration. Ce
phénomène tend à augmenter l'humidité de l'air. Un niveau trop élevé de l'humidité de l'air est
un frein au bon développement des plantes et peut contribuer à l'apparition de condensation
(point de rosée) favorisant alors le développement de maladies fongiques et de moisissures.
Afin de réguler cette humidité, la ventilation est privilégiée mais cela impose d'extraire l'air à
l'extérieur de la serre (remplacé par de l'air plus sec mais frais) ce qui pose des problèmes en
hiver. En effet, les producteurs se trouvent devant une situation paradoxale : ils doivent
19
Chapitre 1 : Introduction
chauffer le volume de la serre mais également le ventiler si bien qu'une part de l'énergie
dépensée s'avère l'être inutilement. Des systèmes de déshumidification par condensation
(pompe à chaleur à déshumidification), des systèmes d'absorption de l'humidité ou bien des
échangeurs de chaleur sur l'air extrait sont parfois utilisés pour limiter ces pertes. Ces
aménagements permettent des économies substantielles mais sont encore loin d'être
généralisés.
3. Gestion des intrants
Parmi les possibilités offertes pour réguler la température et l'humidité du volume de la serre,
l'arbitrage est opéré selon des critères économiques, c’est-à-dire en recherchant une moindre
consommation d'intrants (énergie et eau en particulier).
Le coût des énergies utilisées pour l'exploitation des serres représente actuellement 22% du
prix de revient pour les productions maraîchères et environ 11% pour l'horticulture
ornementale (ADEME, 2007). La part importante de ces coûts ainsi que l'augmentation du
prix des énergies conduit à privilégier de plus en plus les solutions économes, en particulier
préférer la ventilation naturelle à la ventilation forcée, optimiser l'utilisation des apports
énergétiques solaires et limiter l’usage des systèmes de refroidissement et de
déshumidification (ou pour le moins opter pour des systèmes de refroidissement et de
déshumidification économes en énergie).
Cependant, la ventilation naturelle est très dépendante des conditions climatiques extérieures
et le rayonnement, principal moteur du climat dans la serre, est fonction de la saison et de
l'ensoleillement. Ainsi pour pouvoir choisir les solutions les moins coûteuses en énergie, il
importe de maîtriser les facteurs exerçant une incidence sur le climat sous serre, de bien
comprendre les phénomènes physiques le régissant et de maîtriser leurs interactions. La mise
au point d’un outil de prédiction de la distribution du climat sous serre en fonction des
conditions climatiques extérieures devrait permettre une meilleure gestion et un meilleur
contrôle du climat pour répondre aux besoins des cultures.
III. Questionnements
Le développement d’un outil permettant de prédire le climat sous serre en fonction des
conditions climatiques extérieures données requièrt une connaissance exhaustive des
mécanismes physiques régissant le climat. Il est alors essentiel de savoir quels sont les
20
Chapitre 1 : Introduction
facteurs affectant le plus le climat et comment le climat sous serre évolue en fonction de la
variation des paramètres climatiques extérieurs au cours d’une journée.
Cette interrogation se décline selon les questionnements suivants :
quel est l’impact du rayonnement sur l’évolution du climat dans la serre ?
dans quelles proportions la ventilation influence-t-elle la température et le taux
d’humidité ?
comment le végétal réagit-il à des conditions climatiques variables et dans quelle
mesure modifie-t-il le climat ?
Ces questionnements, non exhaustifs, peuvent trouver diverses réponses selon le moment de
la journée. Il s’agit donc de bien prendre en compte les mécanismes physiques qui affectent le
climat et son évolution temporelle. Ainsi, il pourrait être apporté des éléments de réponses
aux questions telles que :
à quel moment et à quels endroits les risques de condensation sont-ils les plus
importants ?
comment les besoins en chauffage évoluent-ils à l’échelle de la journée ?
quand l’irrigation des plantes est-elle nécessaire ?
Les questions concernant l’évolution du climat peuvent être nombreuses, il n’en reste pas
moins qu’elles ne peuvent pas être décorrélées d’interrogations concernant la distribution
spatiale des paramètres climatiques :
la ventilation permet-elle de « réguler » l’humidité et la température en tout point de la
serre ?
à quel endroit la condensation est-elle le plus susceptible d’apparaître ?
est-ce que le chauffage se répartit uniformément et sert-il à répondre aux besoins des
cultures ?
quelle est la localisation idéale d’un capteur permettant de piloter le climat sous serre et
d’obtenir des conditions climatiques optimisées pour le développement du végétal ?
L’ensemble de ces questionnements met donc en exergue les besoins de disposer à la fois
d’information sur la distribution spatiale et l’évolution temporelle du climat sous serre.
IV. Objectifs
Afin de réduire au maximum les consommations d'intrants dans les serres (d'énergie et d'eau
principalement), prédire au mieux la distribution du climat dans les serres et son évolution est
21
Chapitre 1 : Introduction
donc un atout majeur. De nombreux travaux de recherche ont été entrepris depuis les années
1980 pour développer des modèles prédictifs et ces modèles ont été améliorés
progressivement. Ainsi, l'homogénéité spatiale des paramètres physiques a constitué
l'hypothèse de base des premiers modèles. S'en sont suivies des approches plus complexes
permettant de rendre compte des hétérogénéités du climat dans la serre. Ces approches ont pu
être mises en œuvre grâce au développement des logiciels de CFD (Computational Fluid
Dynamics) et à l'accroissement de la puissance de calcul des ordinateurs. Alors
qu'initialement, l'attention était essentiellement portée à l'aéraulique et donc à l’efficacité de la
ventilation des serres, les modèles ont intégré au fur et à mesure l'ensemble des phénomènes
physiques régissant le climat de la serre : les échanges convectifs, le rayonnement, l'impact du
végétal... Cependant l’essentiel des travaux concernant le climat distribué sous serre porte sur
des études réalisées en régime stationnaire.
Le but des travaux de recherche de cette thèse est de progresser sur l'élaboration d'un
modèle instationnaire. Il sera donc nécessaire d'intégrer dans les modèles, les mécanismes
dont l'évolution temporelle a un impact sur le climat. La prise en compte de la course du
soleil ainsi que celle de l'inertie du sol seront les points principaux des travaux menés.
Ces améliorations seront entreprises par étapes en investiguant différents mécanismes de
manière découplée afin de pouvoir les utiliser dans un modèle considérant la serre dans son
intégralité. Deux situations particulières seront envisagées : la première se focalisera sur une
configuration estivale, période lors de laquelle l'impact de la course du soleil est le plus
sensible ; la seconde se concentrera sur l'étude des conditions nocturnes en hiver, au moment
où la demande en chauffage est la plus importante. Dans les deux situations, le
développement d'un outil de simulation permettra, après validation, de réaliser des bilans
hydriques et énergétiques.
V. Démarche et Méthodes
22
Chapitre 1 : Introduction
3ème étape :
Echelle de la serre
Intégration des couplages réalisés aux
étapes 1 et 2
23
Chapitre 1 : Introduction
- correspondant aux grandeurs physiques calculées par le logiciel de simulation qui seront
utilisées dans les étapes de validation du modèle numérique.
24
Projet de thèse
Etude
Campagnes de Données d'entrée et de bibliographique :
mesure validation climatique des
serres, transferts
radiatifs, interaction
avec le végétal
Modèles
Caractérisation Evaluation de Echange de chaleur au
Evapo- l'extinction du niveau du végétal et du
rayonnement dans le Simulations
Transpiration sol
Numériques
couvert végétal
(a) (b)
27
Chapitre 2 : Le système serre
28
Chapitre 2 : Le système serre
29
Chapitre 2 : Le système serre
Figure 4 : Dispositifs de serre existants : cultures sur tablettes (a) et sur géotextile (b),
ouvrants latéraux (c) et de toiture (d), écrans d'ombrage (e et f), chauffage localisé (g).
30
Chapitre 2 : Le système serre
La serre servant de support des campagnes de mesures pour la thèse ainsi que l’ensemble des
éléments qui la composent sont décrits dans le chapitre 4.
31
Chapitre 2 : Le système serre
Figure 5 : Description des phénomènes physiques mis en jeu dans une serre.
1. Principes généraux
1.1. Aéraulique des serres
L’aéraulique regroupe l’ensemble des phénomènes physiques ayant trait au déplacement de
l’air. Ce déplacement de l’air peut résulter des forces d’inertie (i.e. de la vitesse du
déplacement de l’air) lorsque la vitesse du vent est importante ou des forces de flottabilité
lorsque la vitesse de l’air est faible. Dans le premier cas, la différence de pression entre deux
points de la serre est la cause de la mise en mouvement de l’air. L’air se déplace du lieu à
haute pression vers celui à basse pression. Dans le second cas, les différences de température
à l’intérieur de l’air sont le phénomène moteur. En effet, du fait de la gravité, l'air chaud, de
masse volumique plus faible, a tendance à monter tandis que l'air froid, plus lourd, descend.
Ce mécanisme induit donc un mouvement des masses d’air. Ces effets résultants des forces de
flottabilité sont en quelque sorte équivalents à ceux de la poussée d'Archimède et doivent leur
existence à la gravité terrestre.
Pour traiter de l’aéraulique des serres, il convient de considérer d’une part l’intérieur de la
serre et d’autre part l’extérieur de la serre. A l’extérieur la vitesse du vent est de l’ordre de
quelques mètres par seconde alors qu’à l’intérieur elle dépasse rarement les 0.5m.s-1 (Day et
Bailey, 1999). Quelle que soit la serre, on constate que la vitesse de l'air à l'intérieur de la
32
Chapitre 2 : Le système serre
serre ne dépasse que rarement 20% de la vitesse à l'extérieur de la serre (Bartzanas et al., 2004
; Mistriotis et al., 1997b ; Wang et al., 1999a)
L'aération des serres consiste en un échange d'air entre l'intérieur et l'extérieur de la serre.
Pendant la période estivale, la température peut atteindre des valeurs élevées susceptibles
d'endommager les cultures. De même, l’humidité dépasse aussi parfois les valeurs acceptables
pour les plantes. Une humidité trop faible peut amener les cultures à un état de stress ; une
humidité trop importante réduit la croissance des plantes, altère le processus de reproduction
et favorise le développement de maladies (Baptista et al., 2008). L’aération de la serre par
l’intermédiaire des ouvrants (ou un autre système d’extraction ou d’évacuation de l’humidité)
est par conséquent nécessaire afin de ménager un climat optimal à la croissance des végétaux
tout en augmentant les échanges de dioxygène et de gaz carbonique (Nebbali, 2008).
Cette aération est étroitement liée aux conditions climatiques externes notamment à la vitesse
du vent et à la différence de pression entre l’intérieur et l’extérieur de la serre générée par le
champ extérieur des vitesses (Boulard et Baille, 1995).
Même dans le cas de serres en verre fermées, des échanges d’air ont lieu entre l’intérieur et
l’extérieur de la serre. Ils sont principalement la conséquence de la perméabilité et des fuites
plus ou moins importantes de la serre (étanchéité non parfaite). Les avis divergent sur
l'influence du vent sur ces fuites : Wang et Deltour (1999) ont observé qu'il n'y avait pas
d'influence alors que Ruther (1985) a corrélé le taux de renouvellement de l'air dû aux fuites à
la vitesse du vent.
Que ce soit pour les serres ouvertes ou fermées, des échanges d’air ont donc lieu entre
l’intérieur et l’extérieur de la serre. Ces échanges occasionnent des transferts thermiques et
hydriques qu’il est important de contrôler afin de maîtriser les productions sous serre.
Bien que en général le vent extérieur constitue le moteur de l’aération de la serre, la différence
de température entre l'intérieur et l'extérieur de la serre peut également jouer un rôle important
notamment quand le vent extérieur est faible (Baptista et al., 1999) et/ou lorsque la serre est
fermée. En effet, l’hétérogénéité des températures induit une hétérogénéité de la distribution
des densités de l'air. En fait, il existe des écarts de température entre les différentes surfaces à
l’intérieur de la serre (paroi, sol, plantes,...) et l'air. Ces différences induisent un
réchauffement (ou refroidissement) de l’air et donc des gradients de masse volumique entre
33
Chapitre 2 : Le système serre
l’air à proximité de ces surfaces et l’air du reste de la serre. Ce phénomène est aussi connu
sous le nom "d'effet cheminée" (Roy et al., 2002).
Dans une serre fermée le vent n'influe que sur les fuites. La vitesse de l'air est alors
relativement faible (<0.1m.s-1, Wang et al., 1999a). Les forces de flottabilité deviennent le
principal moteur des mouvements internes de l’air. Plus généralement, les forces de
flottabilité jouent un rôle important dès que la vitesse de l'air à l’intérieur de la serre est
inférieure à 0.5 m.s-1 (Mistriotis et al., 1997a).
1.2. Echanges thermiques
Les trois modes de transfert de chaleur coexistent dans une serre : la convection, la
conduction et le rayonnement. Ils interagissent entre eux et de ce fait le climat de la serre est
en partie la résultante de ces interactions.
1.2.a) Convection
Les mouvements de l'air induisent des échanges de chaleur au niveau des interfaces avec les
éléments solides dans la serre. Ces échanges de chaleur associés au déplacement d'un fluide
correspondent au mécanisme de convection. Il s’agit donc de transfert de chaleur avec
déplacement de matière.
La convection joue un rôle essentiel dans les transferts de chaleur qui surviennent à l’échelle
de la serre. A l’intérieur de la serre, la convection est responsable d’un transfert de chaleur
entre, d'une part, le sol, les tablettes de culture, la végétation et les parois de la serre et d’autre
part, l’air environnant. A l’extérieur de la serre, le vent dissipe l’énergie emmagasinée par le
toit et les parois de la serre évacuant ainsi une partie de la chaleur accumulée dans la serre.
Selon l’intensité du vent (pour l'extérieur) et de la vitesse d'air (pour l'intérieur), la convection
peut être libre (transfert essentiellement associé aux gradients de température), forcée
(transfert principalement dû à l’aération de la serre/au vent) ou mixte (intensité des deux
phénomènes du même ordre de grandeur). L’intensité des échanges thermiques est
proportionnelle à la différence de température et dépend fortement du régime de convection
(Roy et al., 2002), cependant, les échanges de chaleur sont généralement moins importants en
convection libre, qu'en convection forcée. En effet, le mouvement du fluide plus important en
convection forcée permet d'évacuer la chaleur et de maintenir un gradient thermique
important entre l'air et les surfaces en contact avec l'air.
34
Chapitre 2 : Le système serre
L’estimation avec précision des coefficients d’échanges est donc cruciale dans l’évaluation de
l’intensité des échanges thermiques et des pertes énergétiques.
1.2.b) Conduction
35
Chapitre 2 : Le système serre
bien qu’il est nécessaire de ventiler la serre (cf. II.1.1.a) pour dissiper le surplus d’énergie.
Inversement la nuit, quand les températures sont trop faibles, le recours à divers systèmes de
chauffage est alors nécessaire. Cette situation assez paradoxale met en lumière les économies
d’énergie qui pourraient être réalisées grâce à une meilleure gestion de l’énergie stockée le
jour dans le sol.
1.2.c) Rayonnement
36
Chapitre 2 : Le système serre
Tout corps peut absorber une onde électromagnétique (Figure 7). Il peut également la
réfléchir ou la transmettre (être traversé). La capacité d’un corps à réfléchir, absorber et
transmettre une onde dépend des propriétés de ce corps appelées respectivement coefficient
de réflexion, d’absorption et de transmission. Ces propriétés sont aussi dépendantes de la
longueur d’onde de l’onde incidente. Ainsi le verre est transparent vis-à-vis des courtes
longueurs d’onde et transmet donc une grande part du rayonnement de courte longueur d’onde
reçu alors qu’il absorbe la quasi-totalité du rayonnement de grande longueur d’onde
(Tableau 1).
Rayonnement réfléchi
Rayonnement absorbé
Rayonnement transmis
37
Chapitre 2 : Le système serre
Figure 9 : Longueur d'onde absorbée par le végétal pour réaliser la photosynthèse [3].
Les éléments de la serre émettent un rayonnement de grande longueur d'onde dans toutes les
directions. Ce rayonnement n'est pas transmis par le verre mais est en grande partie absorbé
et/ou réfléchi : il contribue à chauffer les parois de la serre et donc, par convection, l’air de la
serre.
38
Chapitre 2 : Le système serre
La quantité d’eau contenue dans l’air sous forme de vapeur d’eau a une influence sur la
transpiration des cultures (cf II.2.2.b) ainsi que sur le risque de condensation.
La transpiration des plantes est le phénomène qui permet la circulation de la sève. Plus l’écart
de pression de vapeur d’eau (appelé déficit de pression de vapeur d’eau) entre l’air près du
végétal et l'air à l'intérieur des feuilles (chambre sous stomatique) est élevé, plus la
transpiration est favorisée. Cette transpiration est alors une source de vapeur d’eau dans l’air
au voisinage des cultures. Il résulte de ce phénomène une humidité plus élevée à proximité du
végétal que dans le reste de la serre. Pour une serre ouverte, les vitesses de l'air sont plus
élevées que pour une serre fermée et cette humidité est convectée. L'aération naturelle permet
donc de renouveler l'air à proximité du végétal, ce qui maintient un gradient de pression de
vapeur d'eau important entre le végétal et l'air. La transpiration étant d'autant plus importante
que ce gradient est élevé, ce transfert de vapeur d'eau grâce au mouvement de l’air favorise
donc les pertes hydriques du végétal. Ce transport de vapeur d'eau peut néanmoins avoir lieu
lorsque la vitesse de l'air est faible du simple fait de gradient de concentration de vapeur d’eau
dans l’air. En effet, en présence de différence de concentration de vapeur d'eau, des
phénomènes de diffusion tendent à réduire les hétérogénéités. Il s'agit en fait d'un transfert au
niveau microscopique des composantes de l'air sans mouvement macroscopique de l'air. Ce
transfert a lieu depuis les zones les plus humides (proches du végétal) vers les zones
présentant les concentrations de vapeur d'eau les plus faibles ; en particulier, vers les parois
lorsqu'il y a condensation, puisque le phénomène de condensation "soutire" de la vapeur d'eau
à l'air.
39
Chapitre 2 : Le système serre
Dans une serre, il est préférable d'éviter la condensation pour deux raisons principales:
- l'apparition d'eau libre sur les feuilles des plantes favorise le développement de certaines
maladies (fongiques comme le botrytis par exemple),
- la formation de gouttes sur les parois de la serre réduit sensiblement la transmission du
rayonnement.
La condensation de la vapeur d'eau contenue dans l'air se produit lorsque de l'air humide est
refroidi jusqu'à atteindre le point de rosée (Figure 10). En effet, l'eau se condense quand la
pression de vapeur d'eau dans l'air atteint la pression de vapeur saturante (courbe de
saturation). Cette pression de vapeur saturante augmente avec la température. Ainsi, un air
chaud peut présenter une humidité absolue importante sans qu'il y ait de condensation. Par
contre, à humidité absolue constante, si cet air est refroidi, alors son humidité relative (i.e le
rapport entre la pression de vapeur et la pression de vapeur saturante, usuellement exprimé en
pourcentage) augmente jusqu'à atteindre 100% c’est-à-dire le point de rosée : il y a alors
condensation.
40
Chapitre 2 : Le système serre
1.3.b) CO2
Les tablettes de culture et le substrat des pots, ont un impact moindre que le sol sur le climat.
Néanmoins, il convient d'avoir conscience des effets qu'ils engendrent. Si ces éléments jouent
un rôle moindre que le sol c'est tout d'abord parce que leur inertie thermique est largement
plus faible. Par ailleurs, leurs dimensions sont assez faibles si bien que la conduction dans ces
éléments n'entre que peu en jeu. En revanche, le rayonnement qu'ils reçoivent ainsi que celui
qu'ils émettent doit être pris en compte de même que l'effet de leur température sur la
convection de l'air.
Le sol, quant à lui, a une grande inertie. Le flux de chaleur conductif dans le sol par
conduction est donc le point déterminant. De jour le sol est chauffé par le rayonnement
solaire, la chaleur est conduite dans le sol où elle est progressivement emmagasinée. La nuit,
en l’absence de rayonnement solaire, l'air de la serre devient plus froid que le sol ; la chaleur
est alors restituée.
L'évolution de la température à la surface du sol peut être assimilée à un signal périodique
(Guyot, 1999, Kimball, 1982).
41
Chapitre 2 : Le système serre
Les cultures interfèrent avec le mouvement de l’air dans la serre. En effet, les cultures
constituent un obstacle pour l’air si bien que lorsque l’air pénètre dans le couvert végétal,
l’écoulement subit des pertes de charges. Ces pertes de charges engendrées lors du passage de
42
Chapitre 2 : Le système serre
l’écoulement dans le couvert végétal sont d’autant plus importantes que la densité du couvert
est élevée (Molina Aiz et al., 2006 ; Sase et al., 2012).
La distribution des vitesses et/ou la distribution des pressions dans le couvert végétal sont
donc modifiées lorsque la vitesse de l'air est importante. En revanche pour une serre fermée,
les vitesses d'air sont faibles si bien qu'aucune modification des distributions de vitesse entre
l’air dans et au dessus du végétal n'a été constatée (Wang et al., 1999a)
En plus d'avoir une influence sur le champ de vitesse, le passage de l'air dans le végétal est la
cause du développement d'une couche limite au voisinage des feuilles des cultures. Cette
couche limite est un frein aux différents transferts (dont les transferts hydriques) entre le
végétal et l'air. Il s'agit de la résistance aérodynamique. Cette résistance aérodynamique est un
facteur limitant la transpiration des cultures.
L’évaporation au niveau du substrat et des tablettes se produit notamment au moment de
l’arrosage car l'humidité du substrat des plantes augmente et un déséquilibre avec l'humidité
de l'air apparaît. A l'interface air/substrat, l'air plus sec a donc tendance à se charger en eau.
L'évaporation du substrat est proportionnelle à cet écart d'humidité.
La transpiration est un phénomène entièrement lié à l’activité des cultures. Les plantes
transpirent principalement à travers les stomates des feuilles en même temps qu’elles
absorbent du CO2 permettant ainsi la photosynthèse. La transpiration de la plante modifie
donc la teneur en eau de l'air de la serre.
La transpiration des cultures est liée au déficit de pression de vapeur d’eau (différence entre la
pression de vapeur saturante à la température du végétal et la pression de vapeur de l'air
ambiant). Les résistances aérodynamique et stomatique modulent la transpiration du végétal.
La première est principalement fonction des caractéristiques de l'écoulement de l'air alors que
la seconde représente la résistance aux transferts hydriques exercée par les stomates : elle
dépend, outre du type de plante, et des conditions climatiques de l’air à proximité de la feuille.
La transpiration des plantes est un phénomène au cours duquel l'eau passe de l'état liquide à
l'état gazeux. Ce changement d'état nécessite de l'énergie (chaleur latente) ce qui concourt
ainsi au refroidissement des plantes et de l’air environnant.
43
Chapitre 2 : Le système serre
La condensation peut avoir lieu dans les conditions décrites précédemment (cf. I.1.3.a) au
niveau des végétaux et au niveau des parois (lorsqu'elles sont froides). Elle entraîne des
échanges thermiques inverses de ceux de la transpiration : l’eau passe de l’état gazeux à l’état
liquide en libérant de la chaleur.
L'impact des végétaux sur le climat se traduit par le au refroidissement du système serre en
période estivale mais cause une augmentation de l'humidité qui, si elle n'est pas contrôlée,
peut entraîner de la condensation lorsque la température diminue. Certes la condensation
réchauffe l'air puisqu'elle libère de la chaleur, mais les risques sanitaires sont de loin plus
importants si bien que les exploitants essaient au maximum de réduire les risques de
condensation.
44
Chapitre 3 : Formalisation mathématique
La mécanique des fluides est l’étude des fluides (liquides et gaz) et des forces internes
associées qui permet entre autre de décrire le comportement d’un fluide (son déplacement,
l’évolution de ses paramètres physiques, …). Les équations décrivant le mouvement d'un
fluide découlent à la fois de la mécanique classique (principe fondamental de la dynamique)
45
Chapitre 3 : Formalisation mathématique
Ces équations de bilan sont appliquées sur chaque élément de volume de l'espace d'étude
discrétisé. Dans le cas d'un écoulement turbulent (cas fréquemment rencontré dans les serres),
des équations supplémentaires doivent être résolues. Les équations de bilans "classiques"
doivent alors être adaptées pour intégrer ce caractère turbulent de l’écoulement. Ces
adaptations imposent d'introduire de nouvelles inconnues. Il faut alors faire des choix de
modèles de fermeture (introduction de nouvelles équations de transport pour les quantités
turbulentes) permettant la résolution de l’ensemble des équations.
46
Chapitre 3 : Formalisation mathématique
∂ρ r
+ ∇ ⋅ ( ρv ) = S m
∂t 1
{
424
b
3 {
c
(1)
a
avec :
a : variation locale de la masse en fonction du temps ;
r
b : variation de masse due au champ de vitesse v = f (u , v, w) : la particule fluide de masse
r
volumique ρ est convectée à la vitesse v .
c : termes sources ou puits correspondant à la création ou à la disparition de matière (par ex,
dans une serre, l’air perd de sa masse lorsque la vapeur d’eau se condense sur une paroi).
Cette équation ne suffit pas à elle seule pour décrire l’aéraulique de la serre. La description de
la dynamique de l’écoulement (champ de vitesse) s’appuie en plus sur la résolution de
l’équation fondamentale de la dynamique (équation de conservation de la quantité de
mouvement).
2. Equation de conservation de la quantité de mouvement
A partir de la relation fondamentale de la dynamique qui traduit l’équilibre des forces
exercées sur une particule fluide, on peut écrire :
D(ρv ) ∂ (ρv )
r r
+ ρv ⋅ grad (v ) = ∑ δF ext
r r
= (2)
Dt ∂t
Les forces extérieures sont de deux natures :
- forces de volume : forces dues à des champs de forces externes (ex : forces dues au champ
de gravité, poussée d’Archimède,…) ;
- forces de surfaces : forces de pression, forces de frottement ;
47
Chapitre 3 : Formalisation mathématique
∂(u i ) ∂u
= − ∂p + µ ∂ ²u i − ρgδ i 3 + ρFi
ρ + ρu j ⋅ i (8)
∂t ∂x j ∂xi ∂x ²
j
Ces 3 équations ainsi que celle de conservation de la masse permettent de traiter les
problèmes d’aéraulique pure tant que l’on est en régime de convection forcée (Papadakis et
al., 1992), c'est-à-dire tant que les forces de flottabilité sont négligeables devant les forces
d’inertie. Dans le cas de problèmes en régime de convection libre ou mixte (Lamrani et al.,
2001) qui correspondent par exemple aux cas des serres fermées ou ouvertes mais avec de
faibles vitesses d’air à l’extérieur, les forces de flottabilité ne peuvent pas être négligées. Ces
forces de flottabilité sont induites par les gradients de densité, eux-mêmes associés aux
gradients de température et d'humidité à l'intérieur du fluide. Dès lors, il faut disposer d’une
équation décrivant au moins la variation de la température de l’air : il s’agit de l’équation de
l’énergie.
48
Chapitre 3 : Formalisation mathématique
∂
(U int + Ec ) = − ∫∫ ρ U int + v² .vr.nrdS + W ⋅ + Q ⋅ (11)
∂t S 2
Il faut également exprimer la puissance W ⋅ des forces qui s'exercent sur le volume de la
particule fluide :
- forces de gravité
- forces de pression
- force de frottement visqueux
rr rr rr
W ⋅ = ∫∫∫ ρg.vdV + ∫∫ − p.n.v dS + ∫∫ − σ .n.v dS
1V
4243 1S
4243 S
1 4243 (12)
Puissance du poids Puissance des forces de pression Puissance des forces de frottement visqueux
On a aussi :
Q ⋅ = ∫∫∫ q ⋅ dV (13)
V
gaz carbonique) comme un gaz parfait. Au final, on obtient (en notation tensorielle) :
49
Chapitre 3 : Formalisation mathématique
∂T ∂T ∂ ∂T ∂u i
ρC p + u j ρC p = λ − p +Φ
{ + Sh
1
424 ∂3
t ∂x j ∂x j ∂x j ∂xi e { (14)
a
14243 14243 1 424 3 f
b c d
∂(u i ) ∂u ∂p ∂ ²u i
ρ + ρu j ⋅ i = − + µ − ρgδ i 3 + ρF (15)
∂x ²
∂t ∂x j ∂ x i j
ρC ∂T + u ρC ∂T = ∂ λ ∂T − p ∂u i + Φ + S
p ∂t
∂x j ∂x j ∂x j ∂xi
j p h
On dispose donc de cinq équations : une pour la conservation de la masse, trois pour la
conservation de la quantité de mouvement et une pour la conservation l’énergie. Il y a 6
inconnues : la masse volumique, les trois composantes de la vitesse, la pression et l’énergie
interne.
Il est donc nécessaire de disposer d'une autre équation afin de pouvoir résoudre le système. Il
s'agira d'une équation d'état pour l'air.
5. Equation de Boussinesq
L'équation d'état utilisée découle de l'hypothèse d'incompressibilité du fluide, hypothèse
v
valable dès lors que le nombre de Mach, M = est inférieur à 0.3 ; v étant la vitesse de
c
l'écoulement du fluide et c la vitesse du son. Le nombre de Mach est largement inférieur à 0.3
puisque dans les problèmes liés aux serres, les vitesses de l’air ne dépassent que rarement
0.5 m.s-1 à l’intérieur et quelques mètres par seconde à l’extérieur. Néanmoins,
l’approximation de fluide incompressible ne permet pas de prendre en compte les variations
50
Chapitre 3 : Formalisation mathématique
de densité dues aux variations de température et donc de rendre compte des mouvements du
fluide résultant de gradients de température (i.e la convection naturelle).
En présence de gravité, des différences de température induisent des hétérogénéités de densité
qui peuvent, sous l’effet de la poussée d’Archimède, mettre le fluide en mouvement.
L’approximation de Boussinesq permet de prendre en compte ces variations de densité tout en
ne considérant pas les effets de compression liés à la vitesse des ondes sonores.
Cette approximation a été proposée par Oberbeck (1879) et Boussinesq (1903). Elle peut être
expliquée comme suit :
Il existe des écoulements dans lesquels la température, et par conséquent la masse volumique,
varient peu. Cependant la force d’Archimède est la cause du mouvement du fluide. Donc la
variation de la masse volumique et des autres propriétés thermophysiques du fluide est
négligeable partout sauf dans la force d’Archimède.
La masse volumique ρ en un point de l’écoulement peut être exprimée par un développement
en série de Taylor autour d’un état de référence ρ ref , Tref et en se limitant à un développement
du premier ordre :
ρ = ρ ref [1 − β (T − Tref )] (16)
1 ∂ρ
avec : β = − , le coefficient de dilatation volumique du fluide (air) à pression
ρ ref ∂T p
1
constante. Pour un gaz parfait β = ≈ 3.10 −3 K −1 pour la gamme de températures considérée
T
dans notre étude.
Le terme de la force d’Archimède g (ρ − ρ ref ) est du même ordre de grandeur que les termes
51
Chapitre 3 : Formalisation mathématique
∂T ∂T ∂ ∂T
ρC p + u j ρC p = λ + Φ + Sh (18)
∂t ∂x j ∂x j ∂x j
Le terme de dissipation visqueuse Φ étant négligeable devant les transferts de chaleur par
diffusion, l’équation d’énergie devient :
∂T ∂T ∂ ∂T S h
+uj = a ⋅ + (19)
∂t ∂x j ∂x j ∂x j ρc p
λ
avec a = , la diffusivité thermique du fluide en m 2 s −1 .
ρc p
mm
où y m est la fraction massique d'une espèce m définie par : y m = avec mm la masse de
∑ mi
l'espèce m, en kg.
A : variation temporelle,
B : variation convective,
C : effet des gradients de concentration sur la diffusion de l'espèce m,
D : termes source et puits de l'espèce m.
Le terme J m s'exprime avec la loi de Fick :
52
Chapitre 3 : Formalisation mathématique
J m = − ρ ⋅ Dm ⋅ ∇y m (21)
où Dm est le coefficient de diffusion de l'espèce m dans l'air.
Le coefficient de diffusion de la vapeur d'eau dans l'air est de 0.251.104cm².s-1.
On peut écrire autant d'équations qu'il y a d'espèces considérées mais pour n espèces, n-1
équations suffisent puisqu'on a en plus la relation ∑y i =1 .
∂
(ρ ⋅ y m ) + ∇(ρ ⋅ y m v ) = −∇J m + S m
r
∂t
On dispose donc de 6 équations avec 6 inconnues : les 3 composantes de la vitesse, la
température, la pression et la fraction massique de la vapeur d'eau le système est donc fermé.
Cependant, ces équations bien que valables pour des écoulements turbulents ne permettent pas
de rendre compte des fluctuations turbulentes puisque qu'elles sont d'échelle trop petite en
comparaison des méthodes de résolution, des choix de maillage et des capacités de calcul des
instruments utilisés. Des modifications doivent donc être apportées afin de prendre en compte
le caractère turbulent de l’écoulement via la modélisation du comportement moyen des
quantités turbulentes.
53
Chapitre 3 : Formalisation mathématique
La turbulence désigne l'état d'un fluide, liquide ou gaz, dans lequel la vitesse présente en tout
point des fluctuations aléatoires entraînant la création de tourbillons dont la taille et la
localisation varient sans arrêt. Les écoulements turbulents se caractérisent donc par une
apparence très désordonnée, un comportement difficilement prévisible (i.e. chaotique) et
l'existence de nombreuses échelles spatiales et temporelles interférant entre elles.
Le comportement complexe des écoulements turbulents est la plupart du temps abordé par la
voie statistique. Ould Khaoua (2006) rapporte que le système d’équations de Navier-Stokes
décrit précédemment reste valable pour les écoulements turbulents rencontrés sous serre
(Tennekes et Lumley 1972, Lesieur, 1997). Cependant les moyens de calcul actuels
permettent difficilement une simulation directe de la turbulence sur des volumes aussi
importants que ceux considérés dans les serres. Comme il est en effet impossible de simuler
des tourbillons dont la taille est inférieure à celle des mailles considérées dans l’approche
numérique, on a alors recours à une approche statistique pour inclure les mécanismes
turbulents dans l’approche numérique.
symbole ′ ( ) , les fluctuations ou écarts par rapport à ces moyennes. Si on utilise les propriétés
usuelles de la moyenne :
(u + v ) = u + v
(λ ⋅ u ) = λ ⋅ u
(u ⋅ v ) = u ⋅ v + u ′ ⋅ v ′ (24)
∂u ∂u
∂t = ∂t
Et les relations suivantes :
54
Chapitre 3 : Formalisation mathématique
u i ⋅ u j = u i ⋅ u j + u i′ ⋅ u ′j
u i ⋅ T = u i ⋅ T + u i′ ⋅ T ′ (25)
( ) ( )
ui′ xi , t = p ′ xi , t = 0
En injectant l’équation (23) dans les équations du système (22), on obtient les équations de
Reynolds :
∂u i
∂x = 0
i
∂u ∂u 1 ∂p ∂ ∂u i
i +uj i = − + ν ⋅ − u i′ ⋅ u ′j + f i
∂t ∂x j ρ ∂xi ∂x j ∂x j
(26)
∂T + u ∂T = ∂ a ⋅ ∂T − u ′ ⋅ T ′ + S h
∂t ∂x j ∂x j ∂x j ρC p
j j
∂ ∂ ym u j
( )
ρ ⋅ ym + ρ ⋅ =
∂
∂x j
ρ
∂y
⋅ Dm ⋅ m − ρ u 'j y m' + S m
∂t ∂x j ∂x j
Le tenseur symétrique - de second ordre - des fluctuations des vitesses u i′.u ′j est appelé
tenseur de Reynolds.
∂ p
N.B Par souci de simplification, les termes − + gz et − gβ (T0 − T )δ i 3 sont remplacés
∂xi ρ
1 ∂p
respectivement par : − et + f i respectivement dans le système d’équations (26).
ρ ∂xi
La décomposition de Reynolds introduit 12 inconnues supplémentaires, les six termes de
corrélations turbulentes u i′.u ′j plus trois termes u i′.T ' et trois termes u i' y m' . Il est donc
nécessaire d'avoir recours à des modèles de turbulence pour déterminer les contraintes de
Reynolds et fermer le système d’équations.
55
Chapitre 3 : Formalisation mathématique
∂u j ∂u i 2
+ k + ν t ∂u k
u i′.u ′j = −ν t +
.δ ij
∂xi ∂x j 3 ∂x k
∂T ν ∂T
u i′ ⋅ T ′ = − a t ⋅ =− t ⋅ (27)
∂xi Prt ∂xi
∂y ν ∂y
u 'j y m' = − Dw t ⋅ w = − t ⋅ w
∂xi Sct ∂xi
1
k étant l’énergie cinétique turbulente : k = u i′u i′
2
La viscosité turbulente ν t n’est pas une propriété intrinsèque du fluide, elle dépend de
2 ∂u k
l’écoulement lui-même. Le terme δ ij ν t tient compte de la dilatation volumique du
3 ∂x k
fluide, il est nul pour un écoulement incompressible. at est la diffusivité turbulente. En règle
νt
générale, on admet que le nombre de Prandtl turbulent Prt = de l’énergie vaut Prt = 1 . De
at
νt
même pour Sc t = (nombre de Schmidt), on admet qu'il est égal à 1. Le problème se
Dw t
∂k ∂k ∂ νt ∂k
+uj = ν + + Pr + ε + f ki (28)
∂t ∂x j ∂x j σk ∂x j
Avec :
Pr : production d’énergie cinétique turbulente par action du mouvement moyen,
ε : dissipation visqueuse,
f ki : production due à l’action des forces de gravité,
56
Chapitre 3 : Formalisation mathématique
∂ε ∂ε ∂ ν t ∂ε ε
∂t
+uj = ν +
∂x j ∂x j
+
σ ε ∂x j k
( )
⋅ Cε 1 Pk + Cε 3 ⋅ f ki − Cε 2 ⋅ ε
3 123 (31)
144 42444
144424443 II III
I
57
Chapitre 3 : Formalisation mathématique
I : transport par diffusion turbulente et visqueuse ; le transport par pression fluctuante est
également englobé dans ce terme. Le nombre de Prandtl-Schmidt turbulent associé à ε est
noté σ ε ,
II : le premier terme traduit la production de la dissipation visqueuse par l’action du
mouvement moyen. Quant au deuxième terme, il exprime la production ou la destruction de la
dissipation due aux forces volumiques, qui peut être liée au champ gravitationnel par
exemple. Dans des situations de turbulence homogène, Lumley (1972) a montré que les
interactions avec le mouvement moyen sont négligeables par rapport aux autres contributions
(Ha-Minh, 1991),
III : terme de dissipation de ε par action de la viscosité.
Ce système composé des équations (27), (28), (30) et (31), forme le modèle k − ε dit standard
avec les valeurs des constantes (Tableau 2) évaluées par Launder et Spalding (1974) à partir
de résultats d‘expériences de décroissance d’une turbulence homogène à l’aval d’une grille.
C ′µ = C µ ⋅ C d Cε 1 Cε 2 Cε 3 σk σε
0.09 1.44 1.92 0.7 1.0 1.3
Tableau 2 : Constantes du modèle k - ε.
D'autres modèles dérivés du modèle k − ε existent : le modèle k − ε realizable et le modèle
k − ε RNG (ReNormalization Group). Le premier modèle inclut un terme supplémentaire
dans l'équation du taux de dissipation afin de mieux prendre en compte le caractère
excessivement diffusif du modèle k − ε standard. Le second modèle dérivé du modèle k − ε
standard améliore le calcul des contraintes et prend mieux en compte le caractère
tourbillonnaire de l'écoulement. Mais le modèle k − ε standard demeure actuellement l’un des
plus utilisés, notamment dans les codes de calcul industriels et dans les applications liées aux
serres. Afin de modeliser les phénomènes turbulents rencontrés sous serre, Nebbali (2008) a
montré que le modèle k − ε standard présente un bon compromis entre complexité, précision
et temps de calcul (John et Whittle, 1992). Reichrath et Davies (2002) en revanche ont montré
que le modèle k − ε RNG donnait dans un certain nombre de cas de figure des résultats plus
précis. Pour autant, la complexité supplémentaire du modèle k − ε RNG n'apporte pas une
amélioration significative des résultats. Il ne sera par conséquent pas retenu dans la suite de ce
travail.
58
Chapitre 3 : Formalisation mathématique
Il est important de noter que puisque le modèle k − ε ne peut pas être appliqué à proximité
des parois, il est le plus souvent associé à une loi de paroi (Chassaing, 2000).
Les transferts thermiques dans les couches du sol d'une serre doivent être pris en compte dès
lors qu'on envisage une approche instationnaire. En effet, le sol a une inertie importante si
bien qu'il joue un rôle de stockage (le jour) et de restitution de chaleur (la nuit). A notre
connaissance, Tong et al. (2009) et Nebbali et al. (2012) sont les premiers à avoir engagé des
travaux de modélisation du climat distribué sous serre considérant l'aspect dynamique du sol
en régime instationnaire.
Le sol n'est le lieu d'aucun transport de matière, le seul mode de transfert est la conduction.
Donc, pour le cas d'un sol homogène, un bilan thermique sur un élément de volume permet
d'établir l'équation de la chaleur pour un transfert unidimensionnel :
∂T ∂ 2T ∂T ∂ 2T
ρc p = λ 2 ou =a 2 (32)
∂t ∂z ∂t ∂z
qui peut se généraliser aux cas de transferts tridimensionnels à :
∂T ∂T
ρc p = λ ⋅ ∆T ou = a ⋅ ∆T (33)
∂t ∂t
λ
avec a = la diffusivité thermique du sol.
ρ ⋅ cp
Il existe pour des cas simples des solutions analytiques à cette équation (Taine et Petit, 2003).
Pour le cas d'un sol à une température initiale homogène Ti et dont la surface supérieure (z=0)
est soumise, à t=0, à la température T0 constante, la surface inférieure étant maintenue à T0, la
solution est :
T ( z , t ) = Ti + (T0 − Ti ) × erf z
(34)
2⋅ a ⋅t
Pour le cas d'un sol à une température initiale homogène T0 et dont la surface (z=0) est
soumise, à t=0, à une température variant de manière périodique Te = T0 ⋅ (1 + cos(ω 0 ⋅ t )) , la
solution est :
ω 0 ω o
T ( z , t ) = Te + (T0 − Te ) × erf
z
+ T0 exp − ⋅ z ⋅ cos ω 0 ⋅ t − ⋅z (35)
2⋅ a ⋅t 2 a 2 a
59
Chapitre 3 : Formalisation mathématique
Enfin, pour le cas d'un sol à une température initiale homogène T0 et dont la surface est
soumise, à t=0, à un flux variant de manière périodique ϕ e = ϕ 0 ⋅ (1 + cos(ω 0 ⋅ t )) , la solution
est :
∞
2 ⋅ϕ0 t
T ( z , t ) = T0 + ∫ erfc(u ) ⋅ du
λ s ρ s C ps z
2⋅ a ⋅t
(36)
ϕ0 a ω 0 ω0 π
+ exp − ⋅ z ⋅ cos ω 0 ⋅ t − ⋅ z −
λ ω0 2a 2a 4
Ces expressions analytiques seront utilisées, dans une phase de mise au point du logiciel de
mécanique des fluides pour valider le modèle de conduction dans le sol en régime
instationnaire et déterminer la conductivité thermique du sol ainsi que la condition initiale de
la distribution de la température dans le sol.
60
Chapitre 3 : Formalisation mathématique
où
s : position
r
Ω : vecteur direction
61
Chapitre 3 : Formalisation mathématique
r
Ω ' : vecteur direction de diffusion
-1
κ : coefficient d’absorption (m )
62
Chapitre 3 : Formalisation mathématique
r ∞ r r
q =∫
0 4π
∫ λ ).Ω.dΩ.dλ
I (Ω (39)
∂T ∂T ∂ ∂T S h
+uj = a ⋅ + + div(q ) (40)
∂t ∂x j ∂x j ∂x j ρC p
On notera que la divergence de q s'écrit comme la différence d'un terme relatif à l’émission
propre de l'élément de volume à la température T et d'un terme relatif à l’absorption dans le
milieu.
∞ ∞ r
r
div ( q ) = 4π ∫ κ λ I λ (T ( s )). d λ − ∫ κ λ ∫ I λ ( s , Ω ).d Ω .d λ
0
(41)
0 4π
14 442 4 4 4 3 0 14 442 4 4 4 3
A B
2hc ²
où I λ (T ( s )) =
0
63
Chapitre 3 : Formalisation mathématique
3. Modèle numérique
Le verre possédant des propriétés optiques différentes selon la longueur d’onde, il est donc
nécessaire de distinguer deux gammes de longueur d’onde : les courtes longueurs d'onde
(CLO) [0,1-3µm] et les grandes longueurs d'onde (GLO) [3-100µm] dans l’approche de
modélisation. Ce modèle bi-bande a été choisi pour sa capacité à séparer les influences des
rayonnements atmosphériques (grandes longueurs d’onde) et solaires (courtes longueurs
d’onde) (Figure 14). Il permet de calculer le flux de chaleur par rayonnement dans toutes les
directions et pour toutes les longueurs d’ondes en utilisant la méthode numérique dite des
ordonnées discrètes (DO).
64
Chapitre 3 : Formalisation mathématique
milieu poreux (prise en compte de l'interaction mécanique du végétal avec l'air de la serre),
sera le lieu de création (termes sources) et disparition (termes puits) d'énergie et de vapeur
d’eau estimées via la résolution de bilans thermique et hydrique.
1. Milieu poreux
Le couvert végétal constitue un obstacle à l’écoulement de l’air, induisant ainsi des pertes de
charge. Le modèle traditionnellement utilisé pour intégrer l’influence aérodynamique du
couvert végétal sur l’écoulement d’air est celui assimilant les plantes à un milieu poreux
(Darcy, 1856). Un milieu poreux est un matériau constitué d’une matrice solide, supposée
rigide (induisant une perte de pression), parcourue par un réseau de pores interconnectés.
Le modèle du milieu poreux est donc employé pour modéliser l’influence aérodynamique du
couvert végétal et se formalise par l’addition d’un terme source dans l’équation de
conservation de la quantité de mouvement. Cet effort de traînée peut être exprimé comme suit
(Boulard et Wang, 2002) :
Sϕ = − LAIv ⋅ C d v 2 (42)
où LAIv (m²/m3) est la densité de surface de feuille c'est-à-dire la surface foliaire par unité de
volume, Cd (sans dimension) est le coefficient de traînée de la totalité du couvert végétal et v
(m s-1) est la vitesse de l’air. Ce coefficient de traînée a été évalué par Haxaire (1999) pour un
m2
LAIv de 14.55 3 pour un couvert de tomates et est égal à 0.32. Molina-Aiz et al. (2004)
m
pour une culture de melon ainsi que Boulard et Wang (2002) pour des laitues ont utilisé cette
valeur dans leurs études. Molina-Aiz et al. (2006) et Teitel et al. (2008) ont également mis en
évidence la faible variabilité de ce coefficient pour des cultures de tomates, de poivrons,
d'aubergines et de haricots (valeurs comprises entre 0.25 et 0.36 pour des LAIv compris entre
3.6 et 12.6).
Pour les faibles vitesses d'air, les milieux poreux peuvent être traités par l’équation de Darcy-
Forchheimer qui correspond au cas de figure où on n’est plus dans le cadre d’application de la
loi de Darcy, i.e. pour des nombres de Reynolds>5. Cette équation exprime la force de traînée
dans le milieu poreux comme la somme d’un terme proportionnel à la vitesse (pertes
visqueuses dues aux frottements, Norton et al., 2007) et d’un terme quadratique de la vitesse
(résistance due aux effets de l’inertie, Fluent User’s Guide, 2006).
65
Chapitre 3 : Formalisation mathématique
µ Cf
S ϕ = − v + ρ 0.5 v 2 (43)
K0 K0
où µ (kg.m-1.s-1) est la viscosité dynamique de l’air, K (m²), la perméabilité du milieu poreux
et Cf (sans dimension), un facteur de résistance inertielle. Pour la gamme de vitesses d’air
typiquement observées dans le couvert végétal, le terme du second ordre de l’équation est
prépondérant et les pertes visqueuses (terme du premier ordre) peuvent être négligées
(Boulard et Wang, 2002). La valeur de Cf est déduite des deux équations précédentes :
Cf
= − LAIv ⋅ C d (44)
K 0.5
Fluent permet également de modéliser le terme de perte de charge S ϕ par une loi de puissance
Cf
Avec C 0 = ρ 0.5
(kg.m-2) et C1=2. Il suffit donc de connaître le coefficient Cd pour
K
modéliser le terme source S ϕ .
66
Chapitre 3 : Formalisation mathématique
(T − Ta )
Qs = 2 LAIv ⋅ ρ a c p a = 2 LAIv ⋅ hs (T f − Ta )
f
(47)
ra
Avec :
ra : la résistance aérodynamique de feuille (s.m-1),
ρa : la masse volumique de l'air (kg.m-3),
Cpa : la capacité calorifique de l’air (J.K-1.kg-1)
Tf : la température de feuille (K),
Ta : la température de l'air au voisinage de la feuille (K),
hs : le coefficient d'échange convectif de chaleur sensible (W.K-1.m-3),
LAIv : la densité foliaire (m² de feuille/m3 de couvert végétal).
Le flux de chaleur latente Ql (transpiration) échangé au niveau du couvert végétal peut être
exprimé par une relation similaire au flux de chaleur sensible. Pour le cas d’une espèce dite
hypostomatique dont les feuilles ne transpirent que d’un seul côté (face inférieure dans le cas
de l’Impatiens) son expression est :
(ω − ωa )
Ql = Lv ⋅ LAIv ⋅ ρ a
f
(48)
ra + rs
Lv est la chaleur latente de vaporisation de l'eau (kJ.kg-1),
ωf , l'humidité absolue de l'air à l'intérieur de la feuille considéré comme saturé (kg.kg-1),
ωa, l'humidité absolue de l'air au voisinage de la feuille (kg.kg-1),
rs, la résistance stomatique de feuille (s.m-1).
67
Chapitre 3 : Formalisation mathématique
des ambitions de la thèse est de résoudre complètement l’équation de transferts radiatifs afin
de pouvoir connaitre le rayonnement absorbé pour chaque couche du végétal. Dans tous les
cas, le rayonnement net Rn est connu.
La température de l'air est également connue via la résolution de l'équation de l'énergie.
La fraction massique de l’eau dans l’air ωa est calculée via la résolution de l’équation de
conservation de la vapeur d’eau.
Quant au terme ωf, l'humidité absolue de l’air contenu dans la feuille (chambre sous
stomatique) il est calculé en supposant que l’air à ce niveau est saturé en eau. Il s’agit donc de
l’humidité absolue calculée à partir de la pression de vapeur saturante de l’air contenu dans la
feuille, l'air étant considéré comme étant à la température de la feuille :
Pvsat (T f )
ω f = 0.622 ⋅
Patm − Pvsat (T )
f
(49)
3802.7 442.7
2
ra (ω − ω a )
Rn − Lv ⋅ LAIv ⋅ ρ a f
T f = Ta + (50)
2 LAIv ⋅ ρc p ra + rs
Tf dépend donc de l’humidité absolue qui elle-même dépend de Tf. Un processus itératif est
donc nécessaire pour déterminer Tf, en supposant une température initiale de feuille Tf, on
peut estimer l'humidité absolue ωf et donc recalculer la valeur de la température. Si cette
température est suffisamment proche de la température initialement choisie, on la conserve
sinon on réitère les calculs en prenant la nouvelle valeur de Tf pour évaluer l'humidité
absolue. On recommence autant de fois que nécessaire jusqu’à ce que la température calculée
soit suffisamment proche de la température précédemment calculée (Figure 16).
68
Chapitre 3 : Formalisation mathématique
Tf initial
Pvsat(Tf )
ωf = 0.622⋅
Patm − Pvsat(Tf )
ra (ω − ω a )
R n − Lv ⋅ LAIv ⋅ ρ a f
T f = Ta +
2 LAIv ⋅ ρc p ra + rs
Tf final
69
Chapitre 3 : Formalisation mathématique
appliqué à la feuille assimilée à une plaque plane de largeur lg et de longueur Lg, sont
disponibles dans la littérature (Tableau 3). Les différentes expressions du nombre de Nusselt
moyen pour les différentes configurations envisageables et les différents régimes de
convection et d'écoulement sont consignées dans le tableau 3.
hs ⋅ d
Nu = ⋅ (54)
λa
Les expressions des nombres adimensionnels sont les suivantes :
g ⋅ β ⋅ ∆T ⋅ d 3 ρ a
2
ρa ⋅ v ⋅ d µ a ⋅ c pa
Gr = ; Re = ; Pr = (55)
µa 2
µ λa
Où g est l'accélération de la pesanteur (m.s-2)
1
β est le coefficient de dilatation, dont une approximation est (K-1).
Ta
Dans la littérature traitant des applications liées aux serres, on trouve de nombreuses
expressions pour évaluer la résistance aérodynamique d'une feuille. Elles sont exprimées soit
directement en fonction de la vitesse et/ou de la température, soit via l'expression du
coefficient d'échange ou celle du nombre de Nusselt. Afin de permettre la comparaison, ces
expressions ont été harmonisées et exprimées en fonction du nombre de Nusselt pour les
différents modes de convection (Tableaux 4, 5 et 6) et des dimensions caractéristiques de la
feuille (Lgfeuille et lgfeuille sont la longueur et la largeur maximales moyennes de la feuille)
En convection libre (Tableau 4), le nombre de Nusselt s'exprime en fonction de la différence
de température (nombre de Grashof) élevée à une puissance dépendant du caractère de
l'écoulement. On remarque que le coefficient de proportionnalité utilisé par les auteurs
présente une variabilité relativement importante même s'il reste à chaque fois du même ordre
de grandeur.
70
Chapitre 3 : Formalisation mathématique
En convection mixte (Tableau 5), les expressions utilisées par les différents auteurs
proviennent de la même source (les travaux de Stanghellini, 1987), et sont donc toutes
semblables. Elles résultent d'une combinaison des expressions obtenues en convection libre et
en convection forcée mais sont de nature un peu différente de celle de la plaque plane.
Enfin en convection forcée (Tableau 6), le nombre de Nusselt dépend principalement de la
vitesse de l'écoulement (nombre de Reynolds) élevée à une puissance dépendant du caractère
de l'écoulement. La variabilité du coefficient de proportionnalité est faible d'un auteur à
l'autre.
On remarquera cependant qu'il demeure un certain flou quant au choix fait par les auteurs sur
la dimension caractéristique utilisée. Le plus souvent, on ne sait pas quel choix est fait ni
quelle est la valeur utilisée. Les expressions du nombre de Nusselt pour les plaques planes
utilisent la longueur de la plaque plane (dans le sens de l'écoulement) comme dimension
caractéristique. La difficulté pour l'application de ces formules réside aussi dans le manque
d'information de la direction de l'écoulement par rapport à la feuille.
Afin de prendre en compte le fait qu'une feuille n'a pas la géométrie d'une plaque plane
rectangulaire, la dimension caractéristique considérée est d = γ 0 ⋅ L g feuille où L g feuille est la
longueur de la feuille et, pour une géométrie correspondant à celle de la feuille de l’Impatiens
de Nouvelle Guinée, γ 0 = 0.581 si l'écoulement est perpendiculaire à l'axe principal de la
feuille ou γ = 0.643 si l'écoulement est parallèle à l'axe principal de la feuille (Parkhurst et al.,
1968, cité par Wiebelt et Henderson, 1978).
71
Type de Régime de
Orientation de la plaque plane
convection l’écoulement
Re 2 << Gr
Turbulent
Nu = 0.54 ⋅ Gr
1
0.04 ⋅ α 1
Nu = 0.09 + ⋅ Gr 3
5 1 1
Nu = 0.13 ⋅ Gr 3
Nu = 0.12 ⋅ Gr 3
π
10 9 < Gr < 1012
1
3
3
Laminaire 1 1
3 2 0.17 3 + 3⋅c (α )
Nu = 0.68 ⋅ Re 2 + 0.57 ⋅ Gr Nu = 0.68 ⋅ Re 2 + Gr
3 3 3 3 3 3
5 4
Nu = 0.68 ⋅ Re + ⋅ Gr 5
Mixte 10 < Gr < 10
3 9
1 + ⋅c(α )
5
Re 2 ≈ Gr
Turbulent 1
Nu = 0.03 ⋅ Re 5 + 12.1 ⋅ Gr
12 3
Re 2 >> Gr Turbulent 4
Nu = 0.03 ⋅ Re 5
Re > 5 ⋅ 10 5
L'écoulement dans le couvert végétal étant laminaire (Re<105 et 104<Gr<109), les expressions
retenues seront donc les suivantes :
Convection libre :
1
Nu = 0.40 ⋅ Gr 4
Convection mixte :
(
Nu = 0.40 ⋅ Gr + 6.92 ⋅ Re 2 ) 1
4
Convection forcée :
1
Nu = 0.65 ⋅ Re 2
73
Expression réexprimée Type de végétal Dimension caractéristique Expression utilisée par
1
Nu = 0.72 ⋅ Gr 4 Chrysanthème Non précisée Yang (1995)
1
Nu = 0.54 ⋅ Gr 4 Laitue Non précisée Boulard et Wang, 2002
1
Nu = 0.50 ⋅ Gr 4 Rose Non précisée Katsoulas et al., 2007
Nu = 0.25 ⋅ Gr 0.3 Rose Non précisée Katsoulas et al., 2007
2
1 d=
Nu = 0.37 ⋅ Gr 4 Geranium 1 + 1 Montero et al., 2000
l g feuille Lg feuille
1
Nu = 0.24 ⋅ Gr 4 Ficus benjamina L. Non précisée Denden et Lemeur, 1999
1
Nu = 0.49 ⋅ Gr 4 Tomate Non précisée Wang et Boulard, 1999
1
Nu = 0.12 ⋅ Gr 3
Tomate Non précisée Wang et Boulard, 1999
(écoulement turbulent)
Tableau 4 : Expressions du nombre de Nusselt en convection libre disponibles dans la litterature relative aux cultures sous serres.
Nu = 0.63
∆T
d
+ 0.83 ⋅ v( d) 1
2
Bananier d = lg feuille + Lg feuille Demrati et al.,2007
Tableau 5 : Expressions du nombre de Nusselt en convection mixte disponibles dans la litterature relative aux cultures sous serres.
Expression réexprimée Type de végétal Dimension caractéristique Expression utilisée par
1
Gazon, pommier, café, coton, maïs,
Nu = 0.61 ⋅ Re 4 d=0.07m Sentelhas et al., 2006
melon, soja, raisin, tomate
1
Nu = 0.69 ⋅ Re 4 Poivron Diamètre moyen de la feuille Möller et al., 2004
4
Nu = 0.036 ⋅ Re 5 Laitue Non précisée Boulard et Wang, 2002
1
Nu = 0.60 ⋅ Re 4 Rose Non précisée Katsoulas et al., 2007
2
1 d=
Nu = 0.60 ⋅ Re 4 Geranium 1 + 1 Montero et al., 2000
lg feuille Lg feuille
1
Nu = 0.69 ⋅ Re 4 Tomate Non précisée Boulard et al., 2004
1
Nu = 0.69 ⋅ Re 4 Tomate Non précisée Wang et Boulard, 1999
4
Nu = 0.32 ⋅ Re 5
Tomate Non précisée Wang et Boulard, 1999
(écoulement turbulent)
Tableau 6 : Expressions du nombre de Nusselt en convection forcée disponibles dans la litterature relative aux cultures sous serres.
Chapitre 3 : Formalisation mathématique
76
Chapitre 3 : Formalisation mathématique
Différentes fonctions ont été utilisées pour modéliser l'influence du rayonnement sur la
résistance stomatique (Tableau 7).
Expression Type de végétal Expression utilisée par
f 1 (G ) = 1 + [exp(a1 ⋅ (G − b1 ))]
−1
Tabac Avissar et al, 1985 (57)
f 1 (G ) = 1 +
a1
Tomates Stanghellini, 1987
(G − b1 ) ⋅ c1 (58 )
1 + a1 ⋅ G
f 1 (G ) = non précisé Farquhard, 1978 (59)
1 + b1 ⋅ G
a1 + ⋅G 9 espèces
f 1 (G ) = Baille et al., 1994a (60)
b1 + ⋅G ornementales
f 1 (G ) = 1 + [exp(a1 ⋅ (G − b1 ))]
−1
Bananier Demrati et al., 2007 (61)
Tableau 7 : Expression de la fonction exprimant l'influence du rayonnement sur la
résistance stomatique (a1, b1, c1 : constantes expérimentales).
Boulard et al. (1991) pour une culture de tomates a trouvé que l'expression (57) minimisait les
écarts avec les observations réalisées. Baille et al. (1994a) pour son étude sur la résistance
stomatique de neuf plantes ornementales (dont l'Impatiens) et Demrati et al. (2007) pour une
culture de bananier ont également utilisé cette formulation.
Pour modéliser l'influence de la température sur la résistance stomatique d'une feuille de
f 3 (T ) = 1 + a2 ⋅ (T − T0 ) , T0 étant la température
2
tomates, Stanghellini (1987) a choisi
correspondant à la résistance stomatique minimale, alors que Boulard et al. (1991) ont trouvé
(
que la fonction f 2 (T ) = 1 + exp a 2 ⋅ (T − b2 ) ) permettait d'obtenir un écart plus faible avec les
mesures.
Stanghellini (1987) a modélisé l'influence du déficit de saturation avec la
fonction f 3 (D ) = 1 + a3 ⋅ (D − D0 ) . D'après Boulard et al. (1991), il n'est pas nécessaire
2
d'utiliser une fonction pour modéliser l'impact du déficit de pression de vapeur car le
rayonnement, la température et le déficit de pression de vapeur sont fortement corrélés. En
résumé, deux facteurs (dont obligatoirement le rayonnement) suffisent à décrire l'impact des
paramètres climatiques sur la résistance stomatique. Ceci explique pourquoi Baille et
77
Chapitre 3 : Formalisation mathématique
al. (1994a) ont identifié un impact du déficit de pression de vapeur sans impact de la
température sur la résistance stomatique alors que Boulard et al. (1991) avaient dressé le
constat inverse.
Baille et al. (1994) étant à notre connaissance les seuls à avoir abordé l'impact des conditions
climatiques sur la résistance stomatique de l'Impatiens, la formulation retenue sera celle
utilisée dans leur étude :
rs = rs , min f 1 (G ) f 2 (D )
a1 + ⋅G
Avec f 1 (G ) = et f 2 (D ) = 1 + a 2 ⋅ (D − D0 )
2
(62)
b1 + ⋅G
Les valeurs obtenues par Baille et al. (1994a) sont : rs ,min = 150s.m −1 , a1 = 421W .m −2 ,
78
Chapitre 3 : Formalisation mathématique
79
Chapitre 4 : Dispositifs expérimentaux
La serre utilisée pour les campagnes de mesure est un compartiment de 100 m² (10 m x 10 m)
faisant partie d’une serre de type Venlo de 3000 m² de superficie. La serre est orientée
quasiment suivant l’axe nord-sud si bien que le faîtage est sensiblement perpendiculaire à la
course du soleil. La serre en verre est équipée (Figure 17):
- de tablettes sur lesquelles sont cultivées les plantes,
81
Chapitre 4 : Dispositifs expérimentaux
Ouest Est
Lampes
Chauffage Tablettes
82
Chapitre 4 : Dispositifs expérimentaux
et leur surface totale est estimée à l’aide d’un système optique. Ces plantes sont ensuite
remplacées par d’autres pour maintenir un nombre de plantes constant sur les tablettes.
Les plantes sont irriguées avec une solution fertilisante par subirrigation pendant 15 min, 2 à 3
fois par jour selon les campagnes de mesure.
Comme évoqué précédemment, les campagnes de mesure ont pour objectif d’explorer
l’évolution temporelle de la distribution des paramètres climatiques. Avant de décrire les
campagnes de mesures menées, les capteurs utilisés pour mesurer la température de l’air, les
températures de paroi, les températures de feuille, l’humidité de l’air et les types de
rayonnements (courte longueur d’onde, grande longueur d’onde, net, PAR,…) sont décrits. La
méthode utilisée pour évaluer la transpiration du végétal est également détaillée.
1. Les capteurs
1.1. Mesures des températures de surfaces
Les températures des parois de la serre, des tablettes et de la surface du sol de la serre sont
mesurées à l’aide de sondes platine plates à résistance PT100 (TC ONLINE, France)
(Figure 18).
Les températures des feuilles sont mesurées à l’aide de thermocouple Cu-Cs (Cuivre-
Constantan) collés sur la face inférieure des feuilles. Les températures dans le sol bétonné de
la serre sont également acquises avec ce type de thermocouple.
1.2. Mesures de la température et de l'humidité de l'air
Les caractéristiques de l’air sont mesurées à l’aide de sonde de température et d’humidité
Vaisala HMP45C (Figure 19). Ces sondes sont positionnées à l’intérieur d’abris ventilés de
façon à ce que les mesures ne soient pas affectées par le rayonnement.
83
Chapitre 4 : Dispositifs expérimentaux
a) b) c)
Figure 20 : Capteur de rayonnement global (a), photosynthétiquement actif (b) et de
courte et grande longueur d'onde (c).
1.4. Mesure de la transpiration
La transpiration des plantes est mesurée à l’aide d’une balance Tableau 8 placée en dessous
des tablettes.
4 à 5 plantes (en fonction de leur taille) sont disposées sur un plateau relié à la balance
(Figure 21) par des tiges qui traversent les tablettes. La variation de poids mesurée est
uniquement due à la transpiration des plantes, l’évaporation du substrat étant négligeable (cf.
Chapitre 5.IV.3.1.c).
Les données récoltées pendant les périodes d’irrigation sont inappropriées et sont par
conséquent retirées des fichiers de données. De façon à perturber le moins possible les
expériences, les périodes d'irrigation sont fixées, dans la mesure du possible, en dehors de la
période d'étude, à savoir en début et fin de nuit pour les campagnes du printemps 2010 et en
84
Chapitre 4 : Dispositifs expérimentaux
cours de journée pour les campagnes de l'hiver 2011. Néanmoins des arrosages ont du être
réalisés en milieu de journées lors des campagnes du printemps 2010 afin d’éviter le stress
hydrique des plantes.
85
Chapitre 4 : Dispositifs expérimentaux
a) b) c) d)
Figure 22 : Capteurs de vitesse de l'air boule chaude (a), anémomètre sonique 2D (b),
anémomètre à coupelle (c), anémomètre sonique 3D (d).
2. Centrales d'acquisition
L’ensemble des capteurs utilisés lors des campagnes de mesure est relié à deux centrales
d’acquisition Campbell (CR7 et CR5000). Les capteurs sont interrogés par les centrales toutes
les 3 secondes. Puis les valeurs sont moyennées sur 10 minutes avant d’être stockées. Seules
les balances qui permettent de quantifier l’évapotranspiration n’utilisent pas ces centrales.
Chacune d’elles est reliée à un ordinateur qui les interroge toutes les 10 minutes afin de
récupérer la valeur qu’elles affichent.
3. Incertitude de mesures des capteurs
L'interprétation des résultats numériques et des écarts entre valeurs expérimentales et
numériques requiert en premier lieu la connaissance de la précision de mesure des capteurs.
Quand la mesure d'un capteur sert aux données d'entrée, une incertitude forte peut inciter à
réaliser un test de sensibilité. Sur les données de validation, la connaissance de l'incertitude de
la mesure peut justifier en partie l’écart entre mesure et simulation et fixe, dans tous les cas, la
valeur minimale de l'écart en dessous duquel les écarts ne seront pas significatifs. La plage de
mesure ainsi que la précision de l'ensemble des capteurs utilisés est précisée dans le
tableau 8 :
86
Chapitre 4 : Dispositifs expérimentaux
Nom et
Capteur Plage de mesure Incertitude
Constructeur/Distributeur
Vaisala HMP45C,
Humidité de l’air 0.8% à 100% 2% d’humidité relative
Cambell Scientific Ltd
Vaisala HMP45C,
Température de l’air -39.2°C à 60°C 0.2°C
Cambell Scientific Ltd
5% de la valeur PAR/LE,
Cellule silicium PAR
mesurée Solems
NR-lite 2,
Pyrradiomètre +/-30W/m²
Kipp and Zonen
CSAT3, Cambell
Anémomètre 3D 0 à 30m.s-1 0.08 m.s-1
Scientific Ltd
0 à 32kg 1g
Balances Melter-Toledo
0 à 16kg 0.1g
Tableau 8 : Précision des appareils de mesures et plage d'utilisation.
V. Protocoles expérimentaux :
1. Situation hivernale
La campagne de mesure en condition nocturne pendant la période hivernale a été menée de
février à mars 2011. L'objectif de cette campagne de mesure est de caractériser le climat sous
serre en présence de chauffage afin d'améliorer la compréhension des phénomènes physiques
87
Chapitre 4 : Dispositifs expérimentaux
mis en jeu. La campagne de mesure permet aussi de fournir des données d'entrée et de
validation pour les simulations numériques envisagées.
3m
4 tablettes entrée
6m
balance
Compartiment
a)
b)
Figure 23 : Schéma du dispositif expérimental des campagnes de mesures hivernales :
vue de dessus (a) et vue de côté (b).
Durant la nuit, les ouvrants sont maintenus fermés. Le chauffage au sol est en fonctionnement
de manière à maintenir dans la serre une température de consigne égale à 16°C. Si jamais le
88
Chapitre 4 : Dispositifs expérimentaux
flux de chauffage au sol n'est pas suffisant pour atteindre cet objectif, le chauffage par tube
aérien est utilisé. Les plantes sont irriguées tous les matins à 11h. La disposition de l'ensemble
des capteurs utilisés est présentée sur la Figure 23 :
Durant la nuit, seuls les gradients verticaux de température (T) et d’humidité relative (HR)
sont mesurés (gradients horizontaux faibles). Trois capteurs sont utilisés :
- un capteur en dessous des tablettes (HR1 et T1, HR2 et T2)
- un capteur au milieu de la végétation (HR3 et T3, HR4 et T4)
- un capteur au dessus du couvert végétal (HR5 et T5, HR6 et T6)
Le rayonnement net est mesuré à la fois au dessus et en dessous du végétal. La différence
entre les deux mesures permet d’accéder au rayonnement absorbé par le végétal.
La transpiration et la température des feuilles à 3 hauteurs dans le couvert à 5 positions
différentes sont mesurées pour quantifier l'activité du végétal.
La mesure du flux de chaleur au niveau du sol permet de connaître l'intensité du chauffage qui
est variable comme précisé précédemment. Les thermocouples placés en profondeur dans le
sol (0 cm, 2.5 cm et 5 cm) permettront d'estimer la valeur de la conductivité thermique du sol.
Enfin une attention particulière est donnée à la mesure des températures des parois de toiture
qui sont fortement affectées par les pertes radiatives.
2. Situation estivale
Les campagnes de mesures en période estivale, réalisées de mai à juin 2010, se sont déroulées
en deux étapes :
une campagne préliminaire qui avait pour but, d'une part, de caractériser la
décroissance du rayonnement dans le végétal afin d’évaluer le coefficient d’extinction
du rayonnement et de l’intégrer dans le sous modèle de culture du logiciel de
simulation ; d’autre part, de quantifier la part de l’évaporation du substrat dans le
processus d’évapotranspiration
une seconde campagne dont le but était de fournir des données d’entrée et de validation
aux simulations numériques.
89
Chapitre 4 : Dispositifs expérimentaux
Est
Ouest
a)
b)
Figure 24 : Disposition des tablettes dans le compartiment d'étude (campagne
préliminaire) vue de dessus (a) et de côté (b).
Pour la campagne préliminaire, 210 plantes sont disposées sur les 18 m² de tablette. Leur
hauteur moyenne est de 26.5 cm et l'indice foliaire de la culture est 3.5 m².m-². Deux balances,
placées à proximité l’une de l’autre, sont utilisées (Figure 24). Dans une première phase de
cette campagne, 5 plantes sont positionnées sur chacune des deux balances qui mesurent la
variation du poids des plantes pour donner accès à l'évapotranspiration. Cette première étape
nous a permis de vérifier que les deux lots de plantes étaient similaires, c’est-à-dire qu’ils
possédaient la même évapotranspiration.
Dans une seconde phase, les pots de l'une des deux balances ont été recouverts de film
plastique afin de limiter au maximum l’évaporation du substrat et ainsi mesurer uniquement la
90
Chapitre 4 : Dispositifs expérimentaux
transpiration. La balance 1 (Figure 24), sur laquelle les pots sont équipés de film plastique,
mesure donc uniquement la transpiration des plantes. La balance 2, sur laquelle les pots ne
sont pas équipés de film plastique, mesure toujours l'évapotranspiration des plantes
(évaporation + transpiration). Par différence, il est possible de quantifier la part de
l'évaporation dans l'évapotranspiration.
Parallèlement, des capteurs de rayonnement PAR ont été installés au dessus du végétal
(15 cm) et dans le couvert végétal (tous les 5.5 cm) (Figure 25) afin d'estimer la valeur du
coefficient de décroissance du rayonnement de courte longueur d'onde dans le couvert (cf.
Chapitre 3.VI.2.3). Un CNR1, placé au dessus du couvert végétal et à proximité des balances,
mesure séparément le rayonnement de courte et de grande longueur d'onde à la fois incident et
réfléchi. Cette mesure a permis de vérifier que le rayonnement PAR au dessus du végétal
représente 44% du rayonnement de courte longueur d'onde incident ce qui est conforme aux
données de la littérature (Britton et Dodd, 1976).
91
Chapitre 4 : Dispositifs expérimentaux
a)
P1 P2 P3 P4 P5
1.4m
~ 1.2m
~ 1.1m
~ 1.0m
0m
- 0.01m
- 0.03m
- 0.07m
b)
Figure 26 : Disposition des tablettes dans le compartiment d'étude (campagne
"complète") vue de dessus (a) et de côté (b).
Pour cette deuxième partie des campagnes de mesure, les plantes ayant grossi, elles ont été
espacées. Il y a ainsi 48 plantes sur chaque tablette et 4 sur chacune des balances. Leur
hauteur moyenne est de 32.5 cm et l'indice foliaire du couvert végétal est de 4.2 m².m-². Les
ouvrants de toiture sont maintenus complètement ouverts.
92
Chapitre 4 : Dispositifs expérimentaux
Les balances sont placées dans l’axe de la course du soleil (Figure 26a). Les températures de
feuille (Tf1 à Tf15) sont également mesurées pour rendre compte de l'activité du végétal à
trois hauteurs (feuille du bas, du milieu et du haut de la plante). Une ligne de capteur
mesurant le rayonnement (PAR1 à PAR5), l'humidité (HR1 à HR5) et la température (T1 à
T5) de l'air est également installée sur l'axe de la course du soleil (Figure 26b). La distribution
verticale de la température et de l'humidité est explorée à l'aide des capteurs T2b, T2c, HR2b
et HR2c. Enfin, les températures dans le sol (Ts1 à Ts12) sont mesurées afin d'estimer les flux
de chaleur mis en jeu dans le phénomène de stockage-restitution d'énergie dans le sol. Ces
mesures permettent également de déterminer la valeur de la conductivité du sol et de valider
la capacité du logiciel de simulation à simuler les transferts thermiques dans le sol. Elles
permettent également de déterminer la condition initiale du champ de température dans le sol
pour les simulations numériques.
Le CNR1 est placé au dessus du couvert végétal (15 cm) et à proximité de la balance 2. Le
rayonnement net est mesuré sous le végétal afin de permettre de disposer de plus amples
informations concernant le rayonnement de grande longueur d’onde. En effet, le rayonnement
de courte longueur d’onde peut être déduit des mesures réalisées avec les cellules silicium
(rayonnement PAR) ; par différence, on a donc accès au rayonnement net de grande longueur
d’onde.
93
Analyse de l’évolution spatio-
Chapitre 5 :
temporelle du climat sous serre
I. Objectifs
Les objectifs de cette partie sont de présenter l'exploration des données expérimentales issues
des deux campagnes de mesures. La distribution des paramètres climatiques et leur évolution
est présentée dans un premier temps de nuit, en hiver, puis, dans un second temps, la même
analyse réalisée pour la situation de jour en période estivale est exposée. Enfin, les résultats
expérimentaux permettant de fournir les paramètres nécessaires aux simulations numériques
et de valider ou infirmer les hypothèses émises.
Les mesures collectées pour deux nuits contrastées sont présentées : une nuit avec un ciel
dégagé (3 mars 2011) et une nuit avec un ciel nuageux (24 février 2011).
14
400 Nuit claire
Nuit nuageuse
350 12
Rayonnement GLO(W/m²)
Temperature (°C)
300 10
250 8
200 6
150
4
100 Nuit claire
Nuit nuageuse 2
50
0
0 20h 22h 0h 2h 4h 6h 8h 10h
20h 22h 0h 2h 4h 6h 8h 10h Heure
Heure
Figure 27 : Evolution du rayonnement de Figure 28 : Evolution de la température à
grande longueur d’onde incident à l’extérieur de la serre pour une nuit
l’extérieur de la serre pour une nuit nuageuse et une nuit claire.
nuageuse et une nuit claire.
Les nuits avec un ciel clair sont caractérisées par un rayonnement reçu de grande longueur
d'onde plus faible (>250W.m-²) que pour les nuits nuageuses (>350W.m-²) (Figure 27). Les
valeurs du rayonnement de grande longueur d’onde sont constantes durant la nuit, ce qui
illustre des conditions de nébulosité qui se maintiennent au cours de la nuit.
95
Analyse de l’évolution spatio-
Chapitre 5 :
temporelle du climat sous serre
Ce rayonnement de grande longueur d’onde reçu, plus faible pendant les nuits claires,
implique des échanges radiatifs plus importants et donc une température de l’air extérieur plus
froide pendant les nuits claires. On observe que la température diminue durant toute la nuit
lors des nuits claires alors que durant les nuits nuageuses elle diminue au début de la nuit, puis
se stabilise jusqu’au lever du jour (Figure 28).
Les conditions climatiques de ces deux nuits sont donc contrastées tant en terme de
rayonnement qu’en terme de température extérieure. Ces conditions externes affectent
l’évolution du comportement climatique de la serre.
2. Evolution du climat dans la serre
2.1. Température de l'air
20
19
18
17
Temperature (°C)
16
15
14
T6
13 T5
12 T4
11 T3
10 T2
9 T1
8
20h 22h 0h 2h 4h 6h 8h 10h
Heure
Figure 29 : Evolution de la température de l’air à différentes positions à l’intérieur de la
serre (cf. figure 23 pour la position des capteurs).
La Figure 29 présente l’évolution temporelle de la température à l’intérieur de la serre pour
les différentes positions décrites dans le chapitre précédent (Figure 23). Seule l’évolution des
températures pour la nuit claire est présentée puisque l’évolution des distributions des
températures durant les nuits nuageuses est semblable.
Pour l’ensemble des sondes, on observe que la température mesurée varie peu ce qui n’est pas
surprenant puisque, comme précisé précédemment, le chauffage fonctionne afin de maintenir
une température de 16°C dans le compartiment durant la nuit. Les faibles variations de
température observées (vers 0h et 7h) sont liées aux variations d’intensité du chauffage.
Les informations importantes à retenir concernent donc les différences de température à
l’intérieur de la serre. On remarque que le gradient de température sous les tablettes (T1 et
96
Analyse de l’évolution spatio-
Chapitre 5 :
temporelle du climat sous serre
T2) est relativement faible (écart de température de 0.25°C en moyenne) tout comme celui à
l’intérieur du végétal (T3 et T4 : 0.35°C en moyenne). Le gradient horizontal de température
de l’air au dessus du végétal est quant à lui très faible (différence de 0.1°C en moyenne, T5 et
T6).
En revanche, l’écart de température entre l’air sous la tablette et celui au dessus du végétal est
plus sensible ; il atteint 1°C en moyenne et dépasse même ponctuellement 1.3°C. L’écart entre
la température au dessus du végétal et celle dans le végétal est également sensible (0.7°C). Le
gradient de température le plus important se situe donc au niveau de la tablette : l’écart moyen
entre la température de l’air sous tablette et celle dans le végétal est de 1.7°C. On pourra donc
considérer que :
- la température sous les tablettes ne présente pas d’importants gradients verticaux,
- la température de l’air dans le végétal est homogène,
- la température au dessus du végétal l’est également.
La situation suivante est donc observée : l’air sous les tablettes est plus chaud que celui au
dessus des tablettes ce qui doit amener à une situation où, du fait des forces de flottabilité, un
brassage de l’air s’instaure. On peut donc s’attendre à ce que s’établisse un écoulement de
l’air sous forme d’une cellule rotative. Cependant, l’écart entre les mesures de température T5
et T6 n’étant pas significatif, il ne permet pas de présager du sens de rotation de l’écoulement
d’air. A priori le sens de rotation sera déterminé par le lieu où les gradients de température
sont les plus forts. Il est donc intéressant d’observer pour cela les différences de température
entre les parois.
97
Analyse de l’évolution spatio-
Chapitre 5 :
temporelle du climat sous serre
20
20
18
18
16
16
14
Temperature (°C)
14
Temperature (°C)
12
Tp1 12
10 Tp2 10
8 Tp3
Tp4 8
6
6 Tp1
4
4 Tp2
2 Tp3
2
0 Tp4
0
20h 22h 0h 2h 4h 6h 8h 10h
20h 22h 0h 2h 4h 6h 8h 10h
Heure
Heure
Pour les deux types de nuits, la température Tp1 de la paroi Ouest de la serre est supérieure à
la température Tp4 de la paroi Est et les deux températures de toiture (Tp2 et Tp3) présentent
une évolution semblable. Le gradient le plus important se situe donc entre la paroi latérale
Ouest (Tp1) et la toiture (Tp2). Ainsi, les forces de flottabilité sont plus importantes de ce
côté de la serre, entraînant un mouvement de l’air du bas vers le haut plus important. Le sens
de rotation de l’écoulement d’air est donc supposé être dans le sens horaire.
2.3. Flux de chaleur associé au chauffage
L’intensité du chauffage au sol est régulée afin de maintenir une température de 16°C dans la
serre durant la nuit. Si l’écart entre la température mesurée et la température de consigne est
trop important alors on a recours au système de chauffage aérien.
Les Figures 32 et 33 présentent respectivement les flux de chauffage au sol (mesure des
fluxmètres) et aérien pour une nuit claire et pour une nuit nuageuse.
98
Analyse de l’évolution spatio-
Chapitre 5 :
temporelle du climat sous serre
70
Flux Nuit Claire
60 Flux Nuit Nuageuse
Flux Chauffage (W/m²)
50
40
30
20
10
0
20h 22h 0h 2h 4h 6h 8h
Heure
Le flux de chauffage aérien a été évalué en estimant le coefficient d’échange moyen entre l’air
et les parois du tube de chauffage à partir de la connaissance de la température moyenne de
l’eau chaude qui y circule. En mesurant la longueur de l’ensemble des tuyaux de chauffage,
on estime la puissance cédée à l’air. La valeur obtenue est ramenée à la surface du
compartiment de serre afin de pouvoir comparer cet apport à celui du chauffage au sol. Le
flux de chauffage au sol est plus important pour les nuits claires que pour les nuits nuageuses,
ce qui était prévisible compte tenu des pertes radiatives plus importantes. Néanmoins la
différence de chauffage entre les deux nuits est plus marquée en début et en fin de nuit. A
minuit, les flux sont identiques. Ensuite, les pertes continuent de croître pour la nuit claire (la
température extérieure continue à diminuer) alors qu’elles restent constantes pour la nuit
nuageuse (la température extérieure est stable).
99
Analyse de l’évolution spatio-
Chapitre 5 :
temporelle du climat sous serre
d'humidité. L'humidité sous les tablettes (HR1 et HR2, Figure 23) ne présente pas de fort
gradient et est quasiment égale à l'humidité mesurée par la sonde HR6. En revanche, un écart
sensible est mesuré entre le taux d'humidité côté ouest de la serre (HR5) et côté est (HR6).
Cette observation laisse présager du sens de l'écoulement d'air. En effet, l'humidité plus
importante du côté Ouest résulte d'une accumulation de la vapeur d'eau dans l'air au contact
du végétal. Un écoulement d'air tournant dans le sens horaire pourrait donc être supposé. Ceci
correspond bien à l’hypothèse réalisée à partir de l’analyse des températures mesurées.
Enfin, les différences les plus importantes (20 points d’humidité relative environ) se situent
entre l'air au dessus du végétal (HR6) et l'air dans le végétal (HR3). Contrairement à la
température, un gradient sensible d'humidité (7 points d’humidité relative environ) est observé
dans le végétal (HR3 et HR4).
HR2
80 HR1 80
75
70 70
65 HR6
60 60 HR5
HR4
55
HR3
50 50
HR2
45 HR1
40 40
20h 22h 0h 2h 4h 6h 8h 10h 20h 22h 0h 2h 4h 6h 8h 10h
Heure Heure
100
Analyse de l’évolution spatio-
Chapitre 5 :
temporelle du climat sous serre
0.015
0.015 HS6
HS6
Humidité spécifique (kg eau/kg air)
0.014 HS5
18 18
température (°C)
température (°C)
16 16
14 14 Tf moyen
Tf haut Ta interieur
Tf milieu Ta au dessus
12 12
Tf bas T tablette
10 10
8 8
20h 22h 0h 2h 4h 6h 8h 10h 20h 22h 0h 2h 4h 6h 8h 10h
Heure Heure
Figure 38 : Evolution temporelle de la Figure 39 : Evolution temporelle de la
température de feuille à 3 niveaux dans le température moyenne de feuille, de l’air
végétal. dans le végétal, de l’air au dessus du
végétal et de la tablette.
Les températures de feuilles mesurées à trois hauteurs différentes (Figure 38) ne présentent
pas d'importants écarts. Les écarts maximums observés (0.3°C) ne sont pas significatifs.
Le végétal peut donc être considéré comme ayant, de nuit, une température homogène égale à
la température moyenne des températures mesurées. En revanche, les écarts de température
entre le végétal et l'air ou la tablette sont sensibles (Figure 39). Le végétal est à la fois plus
froid que l'air en son sein de 0.6°C et que l'air au dessus du couvert de 1.2°C en moyenne. Un
101
Analyse de l’évolution spatio-
Chapitre 5 :
temporelle du climat sous serre
écart encore plus sensible est observé entre la température du végétal et celle de la tablette
(2.2°C en moyenne).
1000
J1
900
J2
Rayonnement solaire (W/m²)
800 J3
700
600
500
400
300
200
100
0
0:00 4:00 8:00 12:00 16:00 20:00
Heure (h)
Jour J1 J2 J3
Text (°C) 20.7 21.2 24.1
Vext (m.s-1) 5.45 3.1 2.5
Tableau 9: Température et vitesse de l’air moyennes lors des 3 journées étudiées.
102
Analyse de l’évolution spatio-
Chapitre 5 :
temporelle du climat sous serre
34
32 J1
J2
30
Température de l air (°C)
J3
28
26
24
22
20
18
16
14
0:00 4:00 8:00 12:00 16:00 20:00
Heure (h)
103
Analyse de l’évolution spatio-
Chapitre 5 :
temporelle du climat sous serre
faible et les écarts de température s'observent principalement entre l'air au dessus du végétal et
celui dans la culture (Figure 42).
32
T5
30
T6
T7
Température de l air (°C)
28
26
24
22
20
18
16
14
0h 4h 8h 12h 16h 20h 0h
Heure (h)
100
HR5
90 HR6
Humidité relative de l air (%)
HR7
80
70
60
50
40
30
0h 4h 8h 12h 16h 20h 0h
Heure (h)
104
Analyse de l’évolution spatio-
Chapitre 5 :
temporelle du climat sous serre
Néanmoins, ces écarts d'humidité relative étant assujettis aux valeurs de température, il est
plus pertinent de comparer les humidités absolues (Figure 44). La distribution verticale de
l’humidité absolue confirme les observations faites pour l’humidité relative.
0.015
HS7
Humidité absolue de l air (kg/kg) 0.014
HS6
0.013 HS5
0.012
0.011
0.01
0.009
0.008
0.007
0.006
0.005
0:00 4:00 8:00 12:00 16:00 20:00 0:00
Heure (h)
J1
0.014
Humidité absolue de l air (kg/kg)
J2
J3
0.012
0.01
0.008
0.006
105
Analyse de l’évolution spatio-
Chapitre 5 :
temporelle du climat sous serre
Jour J1 J2 J3
Vext (m.s-1) 5.45 3.1 2.5
HAext (kg.kg-1) 7.21 9.03 9.22
Tableau 10 : Vitesse de l’air extérieur et humidité absolue.
Bien que l'évolution de l'humidité ne soit pas identique pour les 3 journées, une même
distribution spatiale de l'humidité absolue s'établit au cours de la journée (Figure 46) ; cette
distribution particulière est constatée pour les 3 journées observées ce qui amène à conclure
que l’intensité de la ventilation n'a pas d'influence prépondérante sur l'hétérogénéité spatiale
de l'humidité absolue.
x 10
-3 6:00 x 10
-3 9:00 x 10
-3 12:00
1 1 1
Différence d humidité
Différence d humidité
Différence d humidité
0 0 0
absolue (kg/kg)
absolue (kg/kg)
absolue (kg/kg)
-1 -1 -1
-2 -2 -2
-3 -3 -3
-4 -4 -4
1 2 3 4 5 1 2 3 4 5 1 2 3 4 5
Position du capteur Position du capteur Position du capteur
x 10
-3 15:00 x 10
-3 18:00 x 10
-3 21:00
1 1 1
Différence d humidité
Différence d humidité
Différence d humidité
0 0 0
absolue (kg/kg)
absolue (kg/kg)
absolue (kg/kg)
-1 -1 -1
-2 -2 -2
J1
-3 -3 -3 J2
J3
-4 -4 -4
1 2 3 4 5 1 2 3 4 5 1 2 3 4 5
Position du capteur Position du capteur Position du capteur
106
Analyse de l’évolution spatio-
Chapitre 5 :
temporelle du climat sous serre
Ces observations ont fait l’objet d’une communication orale lors du congrès Greensys 2011
en Grèce, communication appuyée par la rédaction d’un article (Morille et al., 2012).
2.3. Rayonnement
Le rayonnement reçu par la végétation est présenté sur la Figure 47. Etant donnée la forte
influence des ombres portées de la structure de la serre sur les valeurs mesurées, on a réalisé
un lissage de l'influence des ombres en effectuant une moyenne glissante sur 2h des valeurs
mesurées. Il est plus aisé ensuite de comparer les évolutions temporelles du rayonnement reçu
par la végétation pour les cinq positions (Figure 26) des capteurs de rayonnement situés au
dessus des plantes (Figure 48).
700 700
650 650 P1
600 600 P2
550 P3
550
500 P4
500 P5
Rg-in (W/m²)
450
Rg-in (W/m²)
400 450
350 400
300 350
250 300
200
250
150
100 200
50 150
0 100
0:00 4:00 8:00 12:00 16:00 20:00 8:00 10:00 12:00 14:00 16:00
Time (h) Time (h)
107
Analyse de l’évolution spatio-
Chapitre 5 :
temporelle du climat sous serre
- l'impact des ombres portées sur le végétal peut différer en raison de l’évolution de la
course du soleil au cours des 3 jours,
- le comportement des plantes peut présenter une variabilité.
350
J1
300 J2
Evapotranspiration (W/m²)
J3
250
200
150
100
50
0
0:00 4:00 8:00 12:00 16:00 20:00
Heure (h)
30
Tf7
28 Tf8
Température de feuille (°C)
Tf9
26
T6
24
22
20
18
16
14
0h 4h 8h 12h 16h 20h 0h
Heure (h)
Figure 50 : Evolution temporelle de la température de feuille à 3 niveaux du végétal.
(Temps universel). (cf. Figure 26)
Comme évoqué précédemment, les feuilles des couches supérieures sont plus chaudes que
l’air dans le végétal (T6) alors que les feuilles des couches inférieures sont plus froides. Il en
résulte donc un flux de chaleur sensible positif pour l’air (l’air gagne de la chaleur) au niveau
108
Analyse de l’évolution spatio-
Chapitre 5 :
temporelle du climat sous serre
des feuilles de la couche supérieure et un flux de chaleur négatif pour l’air (l’air cède de la
chaleur aux feuilles) au niveau des couches inférieures du végétal.
2.5. Flux de chaleur dans le sol
0 0 0
Le flux thermique horizontal est positif entre les positions 1 et 2 (Figure 52), nul entre les
positions 2 et 3 (Figure 53) et négatif entre les positions 3 et 4 (Figure 54) ce qui signifie que
la chaleur du sol se propage depuis les côtés de la serre vers son centre. Ceci était prévisible
puisque les positions 2 et 3 sont sous les tablettes et ne sont donc pas exposées au
109
Analyse de l’évolution spatio-
Chapitre 5 :
temporelle du climat sous serre
Variation du flux thermique dans le sol (point1) Variation du flux thermique dans le sol (point2)
Flux7-8
Flux Thermique (W/m²)
Flux8-9 Flux11-12
200 200
100 100
0 0
-100 -100
0h 4h 8h 12h 16h 20h 0h 0h 4h 8h 12h 16h 20h 0h
Heure Heure
Figure 55 : Variation du flux thermique vertical dans le sol pour les 4 positions. (Temps
universel).
L’étude du flux thermique vertical met en évidence les phénomènes de stockage de chaleur
dans le sol durant la journée et de restitution pendant la nuit. En effet, le flux est positif entre
10h et 22h environ (le sol emmagasine donc de la chaleur) et négatif entre 22h et 10h (le sol
restitue alors l’énergie précédemment emmagasinée). Le flux proche de la surface (entre 1 et
3 cm de profondeur) est systématiquement supérieur au flux mesuré entre 3 et 7 cm de
profondeur confirmant donc l’atténuation du flux thermique avec la profondeur et l’existence
d’une profondeur au-delà de laquelle les variations journalières des conditions climatiques
extérieures n’ont aucune influence sur la température.
Les valeurs des flux présentées ici ont été obtenues en considérant une valeur de la
conductivité thermique du béton de 2.39 W.m-1.K-1, valeur issue de résultats présentés dans le
Chapitre 7.
110
Analyse de l’évolution spatio-
Chapitre 5 :
temporelle du climat sous serre
1. Décroissance du rayonnement
Pour évaluer le coefficient d'extinction du rayonnement dans la culture d'Impatiens de
Nouvelle-Guinée qui a été utilisée comme support de notre expérimentation, des capteurs de
rayonnement PAR équi-répartis ont été disposés dans le couvert végétal (Figure 56).
Figure 56 : Disposition des capteurs de rayonnement dans le couvert végétal pour une
plante de 32.5 cm.
Pour mesurer la décroissance du rayonnement dans le couvert végétal, deux journées bien
ensoleillées ont été retenues, caractérisées par un profil de rayonnement solaire en forme de
cloche (22 et 24 avril 2010). Le coefficient d’extinction est estimé au midi solaire avec
différentes conditions expérimentales afin de s’assurer que ce coefficient de décroissance, qui
est une propriété optique du couvert végétal, ne varie pas selon les conditions expérimentales.
Les résultats obtenus en utilisant un écran d’ombrage déployé (24 avril) sont en premier lieu
présentés puis ils sont suivis des résultats sans écran d’ombrage (22 avril) où le rayonnement
direct sera maximum.
La Figure 57 présente la décroissance du rayonnement (cf Chapitre 3.VI.2.3) mesurée dans le
couvert végétal au midi solaire dans une serre équipée d'un écran d'ombrage. Le LAI cumulé
z
représentant l’axe des abscisses correspond au terme LAI . et sous-entend l’hypothèse
z0
d’une distribution régulière de la densité de feuille selon la verticale. Le rayonnement est
normalisé en divisant les valeurs de rayonnement mesurées par le rayonnement incident. On
trace donc :
G (z )
z
− K ⋅ LAI ⋅
=e z0
(65)
G0
111
Analyse de l’évolution spatio-
Chapitre 5 :
temporelle du climat sous serre
G (z )
z
−⋅1.5307⋅ LAI ⋅
=e z0
G0
Selon Wang (2001), la loi de Beer sous estime la valeur du rayonnement résiduel et donc
surestime toujours la fraction de rayonnement absorbée. L'erreur peut parfois atteindre 50%
de la valeur réelle. Elle peut même être parfois plus importante si les rayonnements direct et
indirect ne sont pas considérés séparément. Pour des feuillages denses (LAI compris entre 4,5
et 5,7), Mercado et al. (2006) ont constaté des erreurs de 100% de la valeur mesurée pour les
valeurs de LAI cumulé inférieures à 2.
Ces observations concordent avec les résultats obtenus. L'écart entre la modélisation par la loi
de Beer Lambert (10% en moyenne) et les mesures demeure néanmoins senblable aux valeurs
rapportées par Wang et Mercado.
Dans le cas d'une serre sans écran d'ombrage, le coefficient de décroissance n'est pas
identique (Figure 58) :
112
Analyse de l’évolution spatio-
Chapitre 5 :
temporelle du climat sous serre
G(z )
z
−⋅1.599⋅ LAI ⋅
=e z0
G0
Le coefficient de décroissance obtenu dans le cas d'une serre avec écran d'ombrage (1,53) est
assez peu différent de celui dans le cas d'une serre sans écran d'ombrage (1,60). En revanche,
l'erreur commise en modélisant l'extinction du rayonnement par la loi de Beer Lambert est
plus importante (15%) ; elle reste néanmoins inférieure à ce qu'ont pu observer Wang et
Mercado.
On notera que la décroissance du rayonnement a été mesurée avec des capteurs de type PAR
(longueur d'onde comprise entre 400 et 700 µm). Or, les valeurs du coefficient d'extinction
qui nous intéressent concernent la plage couverte par le rayonnement solaire (de 350 à
2500µm). En analysant le spectre solaire, Monteith (1993) a montré qu’il fallait diviser le
coefficient d’extinction kPAR par 1.34 pour obtenir le coefficient K. Nos valeurs obtenues
corrigées pour le coefficient d'extinction sont donc respectivement 1,15 et 1,17 avec et sans
écran. On obtient donc des valeurs semblables qui mènent à la conclusion de l’indépendance
de la valeur du coefficient d’extinction par rapport à l’utilisation d’un écran.
Ces valeurs semblent supérieures à celles que l'on peut trouver dans la littérature
(cf. Chapitre 3.VI.2.3). Ceci peut s'expliquer par la non validité d'une hypothèse, à savoir
l'uniformité de l’indice foliaire LAI sur toute la hauteur du couvert végétal. Il s'avère qu'en
effet la moitié supérieure du végétal représente 70% de la surface des feuilles (et donc la
moitié inférieure 30%).
113
Analyse de l’évolution spatio-
Chapitre 5 :
temporelle du climat sous serre
A B
Couvert
végétal
sin (θ ) =
L
(1)
L'
on en déduit pour le LAI "traversé" par le rayonnement direct:
LAI
LAI ' =
sin (θ )
(67)
Les conditions expérimentales sont telles que les mesures de rayonnement ne sont disponibles
que sur la plage horaire 10h-17h (8h-15h heure solaire) car en dehors, les ombres des autres
compartiments de la serre interfèrent avec les mesures du fait que le rayonnement solaire
n'entre plus directement dans le compartiment de serre. Sur cette plage horaire, d'après la
Figure 60, la hauteur du soleil varie entre 30° et 48° en avril et entre 38° et 62° en juin.
114
Analyse de l’évolution spatio-
Chapitre 5 :
temporelle du climat sous serre
115
Analyse de l’évolution spatio-
Chapitre 5 :
temporelle du climat sous serre
Lors des campagnes préliminaires, les deux balances ont été placées l'une à côté de l'autre
dans des conditions climatiques quasi-identiques. Cinq plantes de taille équivalente ont été
disposées sur chaque plateau. La mesure de la variation de masse des plateaux permettait donc
d’accéder à l’évapotranspiration des plantes.
Dans un premier temps on s’est assuré que les plantes transpiraient de la même manière sur
les deux balances :
116
Analyse de l’évolution spatio-
Chapitre 5 :
temporelle du climat sous serre
-100
10
-150
5
-200
0
-250
-5
-300
-10
-350
-15
-400
12h 13h 14h 15h 16h 17h -20
Heure (14/04/10) 12h 13h 14h 15h 16h 17h
Heure (14/04/10)
15
-100
10
Evapotranspiration (g)
Evaporation (g)
5
-200
0
-300 -5
-10
-400
-15
-500 -20
20h 0h 4h 8h 12h 20h 0h 4h 8h 12h
Heure (17-18/04/10) Heure (17-18/04/10)
Sur une durée de 18h entre deux irrigations, comprenant une nuit complète,
l'évapotranspiration de la balance 1 atteint 471g (Figure 67) tandis que la transpiration de la
balance 2 est de 462g (Figure 68). L'écart représente moins de 2% avec une incertitude de
l'ordre de grandeur de 2% du fait de la différence de gabarit (et donc d'évapotranspiration) des
plantes sur les deux balances.
117
Analyse de l’évolution spatio-
Chapitre 5 :
temporelle du climat sous serre
Pour les mesures effectuées, les écarts entre les deux balances avant ou après la mise en place
des films plastiques vont dans le même sens.
Ces éléments nous permettent d'affirmer que dans les conditions de nos campagnes de
mesures, l'évaporation représente une part négligeable de l'évapotranspiration. Par
conséquent, elle sera négligée dans le travail de modélisation.
118
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
I. La transpiration
La transpiration est un des indicateurs de l’activité du végétal. En effet c’est par le biais de ce
mécanisme que la sève circule et monte jusqu’au niveau des feuilles pour leur apporter les
éléments nécessaires à la production de matière organique. Ce phénomène est initié par
l’ouverture des stomates des feuilles qui favorisent les échanges hydriques avec l’air. L’eau
est alors libérée dans l’air sous forme gazeuse. L’eau contenue dans la plante passe donc de
l’état liquide à l’état gazeux. Ce changement d’état s’accompagne d’une consommation
d’énergie qui permet aux plantes de se prémunir d’une élévation trop importante de leur
température. La transpiration est donc un phénomène consommateur d’énergie. La plante tire
cette énergie du rayonnement absorbé, ce qui contribue au final, à une moindre augmentation
de la température du végétal en comparaison de ce qu’elle serait si les feuilles étaient un
matériau inerte. Par conséquent le flux de chaleur sensible échangé entre le végétal et l’air est
plus faible que si le végétal était un matériau inerte ce qui conduit à une moindre
augmentation de la température de l’air dans la serre.
On peut ainsi considérer que la transpiration est un élément contribuant au refroidissement de
la serre. Plus elle est importante, moins la température de la serre augmente du fait du
rayonnement.
119
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
Il existe deux principales méthodes utilisées pour évaluer la transpiration en fonction des
paramètres climatiques. La première appelée "méthode directe" est celle que nous avons
retenue dans le modèle de CFD permettant de simuler l'interaction du végétal avec le climat.
La seconde, très utilisée dans les travaux expérimentaux, est la méthode de Penman Monteith.
1. Méthode directe
La méthode directe relie proportionnellement la transpiration au déficit de pression de vapeur
de la feuille, c'est-à-dire à la différence entre la pression de vapeur saturante de l’air à la
température de la feuille et la pression de vapeur de l’air. Par analogique électrique, la
transpiration est inversement proportionnelle à la résistance au transfert hydrique, résistance
composée de la résistance stomatique et de la résistance aérodynamique.
120
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
L'impatiens de Nouvelle Guinée est une plante qui ne transpire que par la surface inférieure
des feuilles. L'expression de la transpiration évaluée selon la méthode directe est différente
selon que l’on considère une feuille ou un couvert végétal (Figure 69):
- pour un m² de feuille :
(ω − ωa ) ρ ⋅ C pa (Pvsat − feuille − Pva )
E=ρ qui s'écrit aussi E =
f
(en W.m-² de feuille)
ra + rs γ ra + rs
- pour un couvert végétal :
(ω − ωa ) ρ ⋅ C pa (Pvsat − feuille − Pva )
E = LAI ⋅ ρ qui s'écrit aussi E = LAI
f
(en W.m-² de surface
ra + rs γ ra + rs
au sol)
rs
RS=rs/LAI
ra RA=ra/LAI
Figure 69 : Résistances aux transferts hydriques pour une feuille et pour un couvert
végétal.
121
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
virtuelle possédant les propriétés moyennes des feuilles du couvert végétal. Ce concept a
permis l'élaboration d'un grand nombre de modèles de transpiration pour une large gamme de
cultures sous serres : rosier (Baille et al., 1990), tomate (Boulard et al., 1991), concombre
(Yang, 1995), laitue (Boulard et Wang, 2002) ou bien encore 9 plantes ornementales dont
l'impatiens (Baille et al, 1994a).
La formule de Penman Monteith, découle du bilan énergétique sur le couvert végétal et
s'exprime pour un végétal amphistomatique (stomates sur les deux faces) :
ρ a ⋅ LAI ⋅ c pa ⋅ Da
∆ ⋅ (Rn − Γs ) +
ra
Lv E m = (71)
r
∆ + γ 1 + s
ra
Et pour un végétal hypostomatique (stomates sur une face):
2 ⋅ ρ a ⋅ LAI ⋅ c pa ⋅ Da
∆ ⋅ (Rn − Γs ) +
ra
Lv E m = (72)
r
∆ + 2 ⋅ γ 1 + s
ra
avec
Em : flux d'évaporation massique (kg.m-2.s-1)
Lv : chaleur latente de vaporisation de l'eau (=2450J.kg-1)
LAI : indice foliaire (m² de feuille .m-² de sol)
Rn : rayonnement net total (W.m-2)
Γs : rayonnement absorbé par le sol (W.m-2), considéré négligeable pour les couverts denses,
ρa : masse volumique de l'air (kg.m-3)
cpa : chaleur spécifique de l'air (J.kg-1.K-1)
Da : déficit de pression de vapeur (hPa)
γ : constante psychométrique (hPa.K-1)
∆ : gradient de pression de vapeur saturante entre l'air et la feuille (hPa.K-1)
ra : résistance aérodynamique (s.m-1)
rs : résistance stomatique (s.m-1)
122
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
D f − Da (P
vsat − feuille − Pva ) − (Pvsat − air − Pva ) Pvsat − feuille − Pvsat − air
∆= = = (73)
T f − Ta T f − Ta T f − Ta
Les modèles de transpiration ont été utilisés pour évaluer la transpiration des cultures
d’Impatiens utilisées durant les campagnes expérimentales. Ces travaux ont fait l’objet d’un
article publié dans la revue Scientia Horticulturae (Vol 152, 2013, pp 80–91) qui l’objet de ce
chapitre :
This implies better knowledge of the transpiration process in order to adapt water inputs to
plant needs. Several simple models have been developed to estimate transpiration. The
Penman Monteith (PM) model and the Direct Method (DM) are two of the most widespread
transpiration models. Yet the PM model was primarily developed for crops grown in open
field conditions and its validity in the case of sheltered crops is only seldom questioned.
123
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
Moreover, the PM model is rarely compared to the DM model, although this could be a
The aim of the study is thus, to analyse the ability of the classical PM model to assess
the transpiration rate inside a greenhouse with a New Guinea Impatiens crop in order to
the PM model with measurements revealed severe discrepancies. A detailed assessment of the
assumptions underlying the PM model showed that the location of the climatic parameters
used to implement the model was a key factor affecting the results. In particular, it was
demonstrated that the temperature and humidity considered in the PM model should be taken
inside the crop and not above the crop, leading to a so-called PM-Modified (PMM) model.
The PMM and the DM models were then compared, leading to similar results in agreement
Secondly, the capacity of the PMM and the DM models to investigate the vertical
distribution of the transpiration inside the crop was discussed. After dividing the canopy into
two layers, the PMM and DM models were applied and validated against experimental data.
Even if such a procedure does not really improve the quality of the results compared to the
Highlights:
- Assumptions are explored, improvements are proposed and PM is compared to Direct Model
- Multilayer model does not provide more accurate results than one layer model
124
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
Nomenclature:
A shelf area, m²
125
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
Rn-abs1 net radiation flux absorbed by the upper mid-layer of the crop, W m-2
Rn-abs2 net radiation flux absorbed by the lower mid-layer of the crop, W m-2
ti time, s
Ta air temperature, K
Tl leaf temperature, K
PM Penman-Monteith
DM Direct Model
PMM Penman-Monteith-Modified
126
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
1. Introduction
The scarcity of water resources has led growers to improve their irrigation strategies in
order to provide crops with their exact water requirements. This requires predicting the plant
water uptake that depends both on climatic conditions and plant physiology. Several models
make it possible to predict plant transpiration. Among them, two models are primarily used in
greenhouse horticulture: the Penman-Monteith model (Monteith and Unsworth, 1990) and the
so-called Direct Model which directly ensues from the formalization of the convective
The Direct Method (DM) links the transpiration to the leaf-to-air vapour pressure
deficit. It is used with climatic data recorded either inside the crop (Table 1) (Boulard et al.,
2004; Fatnassi et al., 2004; Yang, 1995) or above the crop (Kittas et al., 2000; Leonardi et al.,
2000). This model is generally used for evaluating transpiration in Computational Fluid
Dynamics (CFD) studies (Table 1) but could be applied in experimental studies as well. The
Penman-Monteith (PM) equation is derived from both the Direct Method equation and the
The PM equation was primarily developed for open field conditions. This method is
based on a big-leaf and perfectly-mixed-tank approach that assumes homogeneity of both the
thermodynamic conditions within the canopy and the air above the plants. It expresses crop
transpiration as a function of the average air-to-air vapour pressure deficit and the net
radiative flux that are measured at only one point in the greenhouse, normally just above the
crop (Table 1). Authors, however, do not always provide information about the exact location
of the sensors. At best, they mention that they are placed above the crop, without any further
details. Other models developed by Penman (1948), Stanghellini (1987) and Fynn (1993) are
127
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
also available to evaluate transpiration. These models are very similar to the Penman-
Monteith model. The main difference lies in the complexity with which the resistances to
water vapour transfers are considered. The historical evolution and the different
improvements of these models are well detailed in the papers of Prenger et al. (2002), Kirnak
Up until now, few studies have been undertaken to compare the performances of
transpiration models under greenhouse conditions. Prenger et al. (2002) compared the
Penman, PM, Stanghellini and Fynn models, while Villarreal-Guerrero et al. (2012) compared
the PM and Stanghellini models with the so-called heat balance based model (Table 2).
Prenger et al. (2002) obtained the best performance with the Stanghellini model for evaluating
transpiration rate of a Red Sunset Maple nursery. In the case of a tomato crop, Lopez-Cruz et
al. (2008) also obtained the best results with the Stanghellini model. Lastly, for bell pepper
and tomato crops Villarreal-Guerrero et al. (2012) obtained the best agreement with both the
knowledge, however, a comparison with the Direct Method has not yet been done. Since PM
and Stanghellini models seem to have similar performances according to the above studies, it
appears relevant to compare only one of them with the DM model. To do this, the PM model
will be used in the present work since it is more widely used than the Stanghellini model,
which may be acceptable in open field studies. However, in greenhouses and especially in the
vicinity of the crop, several studies have highlighted the heterogeneity of these parameters
128
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
(Demrati et al., 2007; Yang, 1995), casting doubts on the validity of the PM model in its
classical form in greenhouses. For that reason, the aim of the present study is to address the
following points:
In order to answer these questions, the transpiration of a New Guinea Impatiens crop
grown inside a glasshouse compartment was considered. The PM and DM models were first
described, focusing on their application conditions, in particular. The PM model was then
implemented with the commonly adopted choices for climatic parameters and compared to
find the cause of the observed discrepancies between evaluated and measured transpiration
rates. The PM model was then tested with climatic parameters placed in different locations in
order to adapt this model to indoor greenhouse climatic conditions. Using this adaptation, the
Finally, the PMM and DM models were applied to a crop that was divided into two horizontal
layers, tested, compared and then validated against measurements. The vertical distribution of
the transpiration and the contribution of each layer to the total transpiration were also
investigated.
According to the DM approach, for a hypostomatic plant, the latent heat flux due to
129
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
where E is the canopy transpiration rate (kg m-2 s-1), λ is the latent heat of water evaporation
(J kg-1), ρ a is the air density (kg m-3), c pa is the specific heat of air (J kg-1 K-1), γ is the
psychrometric constant ( γ = 66 Pa K-1), ra and rs are the leaf aerodynamic and stomatal
resistances (s m-1), respectively, LAI is the Leaf Area Index (m² m-²), Pvsat-l is the saturated
vapour pressure of the air (assumed at the leaf temperature) in the stomata chamber (Pa), Pv-a
is the water vapour pressure of the air (Pa), and Dl (=Pvsat-l -Pv-a) is the leaf-to-air vapour
One of the difficulties of the DM application is that it requires knowledge of the leaf
temperature. In most previous experimental studies (Table 1), even if vertical leaf-temperature
gradients were observed inside the canopy (Demrati et al., 2007; Yang, 1995), only one value
was used in the DM model to characterise the whole canopy temperature. Classically, this
within the canopy (Kittas et al., 2000; Yang, 1995) or placed at pre-determined levels
(Boulard et al., 1991). Concerning the air temperature, the studies listed in Table 1 show that
the above-canopy air is often considered as the reference temperature to be used to express
parameters in Eq. 0 (Bartzanas et al., 2004; Boulard et al., 1991; Katsoulas et al., 2001a;
Kittas et al., 2000). This choice was largely discussed by Yang (1995) who clearly showed
that a model based on the within-canopy air temperature provided the best results. Moreover,
Yang’s conclusion is emphasized by the fact that the aerodynamic resistance characterises the
resistance to water vapour transfer of the leaf boundary layer and not of the canopy boundary
layer.
130
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
The Penman-Monteith model for hypostomatic plants is derived from the combination
of the energy balance equation and the latent/sensible heat flux expressions:
2 ⋅ ρ a ⋅ LAI ⋅ c pa ⋅ Da
∆ ⋅ Rn +
ra
λE =
r
∆ + 2 ⋅ γ 1 + s
ra
where R n is the net radiation flux absorbed by the crop (W m-2), Da = Pvsat −a − Pv −a is the air-
to-air vapour pressure deficit (Pa) and Ta is the air temperature (K).
The ∆ -parameter is the approximation of the slope of the saturated water vapour pressure
Pvsat −l − Pvsat −a
∆= (Pa K-1)
Tl − Ta
where Tl is the leaf temperature (K). In several studies, ∆ is assumed to be equal to the slope
of the saturated water vapour pressure curve, which can be evaluated from the air temperature
only:
The advantage of this formulation is that the leaf temperature is no longer required. This
formulation is all the more reliable since Ta is close to Tl. Other formulations of ∆ are
available in the literature. For example, Kirnak et al. (2002), Prenger et al. (2002) and
Baille (1994b) chose constant values: ∆ = 1.45 ⋅ 10 2 Pa K-1 for the night and
131
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
Concerning the net radiation flux Rn absorbed by the crop, in the studies available, it is
systematically assimilated to the net radiation measured above the crop Rn-a, thus neglecting
the radiation exchanges below the crop (between the crop and the ground or shelf). Baille
(1994a) even assumed that the net radiation above the crop could be assimilated to the
applicable method to estimate plant transpiration since, contrary to the DM model, it does not
use the leaf temperature (which still remains difficult to measure) as a parameter. Moreover,
given the fact that the PM model was built from the so-called “Big-Leaf” approach, meaning
that the canopy is represented by a virtual big leaf that exchanges heat with the above-canopy
air, all the required climatic parameters (air temperature and humidity, radiation and air
speed) are measured above the crop (Table 1). For that reason, the PM model has been
extensively applied for a wide range of crops under greenhouse conditions with varying
Both PM and DM require the knowledge of parameters such as the Leaf Area Index, the
leaf aerodynamic resistance ra and the leaf stomatal resistance rs. According to Villarreal-
Guerrero et al. (2012) these parameters are sometimes difficult to measure or may depend on
several climatic parameters such as the radiation reaching the crop, vapour pressure deficit,
For evaluating ra, the most physical method consists of considering the leaf as a flat
plate and deducing the value of the resistance from the determination of the heat exchange
coefficient induced by the airflow using dimensionless numbers (Appendix A). However, on
the basis of the literature, it can be seen that all the climatic parameters needed to evaluate the
132
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
value of the aerodynamic resistance are not always recorded at the same locations (Table 1),
casting doubts on the relevance of the method. The choice of the right location of the climatic
parameters to be used in the modelling approach was discussed by Yang (1995). According to
this author, and contrary to what is generally done, the method of evaluating resistances to
vapour transfers based on dimensionless numbers requires the use of inside crop climatic
parameters. Using the climatic parameters measured above the crop involves one major
drawback: it means that both water vapour transfer resistances between the crop and the
within-canopy air, and between the latter and the above-canopy air are taken into account,
The leaf aerodynamic resistance value can also be obtained from the inversion of the
transpiration model. This inversion technique may lack genericity, given that the dataset used
for validating the transpiration rates is the same as the one used for estimating the parameters
required for building the model. Moreover, a transpiration evaluation under different climatic
conditions could lead to substantial differences with measurements, which again highlights
the lack of genericity of this method. Yang (1995) claimed that the inversion method is only
useful for calibrating constants of models with climatic parameters measured at a particular
location. Consequently, the values of resistances obtained from the inversion technique
depend on the location of the sensors, thus once again casting doubts on the genericity of
these models.
To simplify the implementation of the models, some authors considered the leaf
(Bartzanas et al., 2004; Demrati et al., 2007; Majdoubi et al., 2009), thus reducing the number
of required climatic parameters and simplifying the experimental setup. This simplification is
133
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
however only acceptable when the convective regime is almost the same during the studied
time period (in particular when the air velocity and the difference between air and leaf
temperature are constant) or when the aerodynamic resistance is negligible compared to the
The leaf stomatal resistance rs is the resistance of water vapour transfer through the
transpiring crop surface. It depends on the physiological behaviour of the crop and especially
on the degree of aperture of the stomata, which is directly linked to the climatic conditions
(Eq. 6). The evaluation of rs using transpiration models is the main objective of several
studies (Baille et al., 1994b; Boulard et al., 1991; Boulard et al., 2004; Kittas et al., 2000;
Rana et al., 1998) or may be one of their objectives (Katsoulas et al., 2001b; Villarreal-
Guerrero et al., 2012; Yang, 1995). To obtain access to the requested parameter, these authors
inverse the transpiration models and deduce the resistance values from measured transpiration
rates. This resistance can be modelled by a multiplicative function (Boulard et al., 1991):
where rs,min is the minimal leaf stomatal resistance, G0 (W m-2) is the global radiation that
reaches the crop, and Da (Pa) is the air-to-air vapour pressure deficit.
For a tomato crop, Boulard (1991) demonstrated that only two of these three functions were
parameters. As a result, the stomatal resistance may be written either as a function of G0 and
Ta or as a function of G0 and Da. Boulard (1991) and Stanghellini (1987) suggested the
a1 + G0
f 1 (G0 ) = (1) and f 2 (Da ) = 1 + a 2 ⋅ (Da − D0 )
2
(2)
b1 + G0
134
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
To our knowledge, Baille et al. (1994b) and Kichah et al. (2012) are the only ones to
have conducted a study to evaluate the stomatal resistance of an Impatiens crop and its
dependence on climatic parameters. This study reveals that Ta is not required as a free
parameter.
The application of the transpiration models to a crop in a two layers approach requires
the value of the radiation at mid-height in order to estimate the net radiation absorbed by each
layer, as well as to calculate the leaf stomatal resistance of the lower layer. The short wave
radiation decay inside the crop due to the absorption of the leaves is provided by Beer’s law:
G ( z ) = G0 ⋅ exp − k ⋅ LAI ⋅
z
z0
where G(z) is the global radiation at a depth z inside the crop, G0 is the global radiation at the
top of the crop, k is the extinction coefficient, z is the depth inside the crop and z0 is the height
of the crop. The net radiation absorbed by each layer may then be evaluated according to the
following expressions:
z /2
Rn − abs1 = Rn − a − Rn − a ⋅ exp − k ⋅ LAI ⋅ 0
z0
z /2
Rn − abs 2 = Rn − a ⋅ exp − k ⋅ LAI ⋅ 0 − Rn −u
z0
temperature (air and leaves) and air humidity, as well as plant characteristics, are needed to
estimate crop transpiration via PM and DM models. In the present section, the experimental
135
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
setup implemented for the purpose of the study and the corresponding recorded climatic
glasshouse oriented north-south and located in Angers (47° 28’ North, 0° 33’ East and 39 m
Scarlet’) were used as plant models and uniformly placed on four shelves, 80 cm above the
ground, to create a dense canopy with a total area of 18 m². The average plant height was 30
cm and the density was 10.6 plants per m². During the field survey, the canopy covered 100%
of the shelves, no shading screen was used and the roof vents were fully open to avoid high
temperatures. Plants were watered twice during the night (at 8 pm and at 4 am) by flooding
the shelves with a complete nutrient solution, and flowers were regularly removed. The
experiments were carried out from May to June 2010. Three clear and sunny days were then
selected for further investigations: 20th May (D1), 3rd June (D2) and 23rd June (D3).
view of the experimental device is depicted in Fig. 1 and described below. The above- and
30-cm-long silicon solar cells (PAR/LE, Solems, Palaiseau, France, +/-10%). The value of the
mid-height radiation inside the crop s deduced from these sensors. According to Monteith
(1993), the extinction coefficient deduced from silicon cells (PAR radiation only) must be
multiplied by 1.34 to obtain the resulting value of k for the whole global radiation spectrum.
In the present study, a value of k = 1.15 was found. Inside the canopy, the cells were equally
distributed every 7.5 cm over the plant height. The above- and under-canopy net radiation
136
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
fluxes (Rn-a and Rn-u, respectively) were measured by two pyrradiometers (NR-lite 2, Kipp and
Zonen, Delft, The Netherlands, +/- 30W/m²). The net radiation Rn absorbed by the crop was
then deduced as the result of these two radiation fluxes. A sonic anemometer (CSAT3,
Campbell Scientific Ltd., Antony, France) also monitored the air velocity (V) 15 cm above the
crop. The temperature (Ta, +/- 0.1°C) and relative humidity (RH, +/- 2% HR) of the air were
measured by shielded and aspirated sensors (Vaisala HMP45C, Campbell Scientific Ltd.,
Antony, France) in the vicinity of the balance at three levels: 15 cm above the crop (subscript
1), at mid-height inside the crop (subscript 2) and at the bottom of the plants, just above the
pots (subscript 3). Finally, the temperatures of leaves of the upper (Tl1), middle (Tl2) and
bottom (Tl3) part of the plant were recorded with copper-constantan (Cu-Cs) thermocouples
All of the above-mentioned parameters were measured every 3 s and averaged online
over 10 min periods with a data logger system (CR3000 and CR7, Campbell Scientific Ltd.,
Antony, France). Data were then analysed using a moving average technique (over 2 h
periods) in order to filter the noise associated with shadows that occasionally covered the
sensors.
Time evolution of the climatic parameters inside the greenhouse as well as the outside
radiation during these experiments can be found in the paper of Morille et al. (2012).
The Leaf Area Index (LAI) was measured by a destructive method using a camera
integrating system. The method was applied four times during the experimental period by
selecting three plants placed on the balance. For the purpose of the study, the plants were
divided into two horizontal layers in order to obtain the LAI vertical distribution.
Measurements showed that the average LAI of the whole plant was 4.2 m2 of leaves by m2 of
137
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
ground covered by the crop, and indicated that the leaves of the upper- and lower-mid layer of
the plant represented 69% and 31% of the total leaf area, respectively. The leaf width and
length were also measured and the values obtained were then averaged to evaluate the
Switzerland, +/- 0.1g) located approximately at the centre of the shelves. Five plants were
placed on the balance. The corresponding latent heat flux was calculated in W m-2 of crop
(m(t i +1 ) − m(t i )) j
λE = λ ⋅
A ⋅ (t i +1 − t i ) n
where λ is the latent heat of vaporisation (J kg-1), m(ti) is the mass (kg) given by the balance at
time ti (s), A (m²) is the area of the shelves, n is the number of plants on the balance tray and j
is the number of plants on the shelves. Data from the balance were directly stored every ten
minutes and then averaged over 2 h periods. Preliminary experiments using plastic films
placed on the pot surfaces showed that the substrate evaporation was negligible (lower than
2% of the total evapotranspiration) compared to the crop transpiration so that no plastic film
3. Results
In this section, as was generally done in previous studies, the PM model was tested
using the climatic parameters measured 15 cm above the crop (Rn-a, RH1, Ta1, V) with the
approximation of ∆ described by Eq. 0 and the net radiation flux measured above the crop
(Rn-a). However, even if authors often choose a constant value for the leaf aerodynamic
138
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
resistance, it was calculated here according to the heat exchange coefficient method described
in Appendix A using d=0.05 m. This method also provides information on the airflow
convection regime and on the temporal variation of the calculated leaf aerodynamic
resistance.
Results indicate that the airflow convection regime is mainly free or mixed inside the
crop, depending on the period of the day. By night, low air-to-leaf temperature gradients are
present, so that the aerodynamic resistance (Fig. 2) is close to 300 s m-1. During the day, the
radiation that reaches the crop induces higher air-to-leaf temperature gradients, the convection
regime becomes punctually free and the aerodynamic resistance varies between 150 and 200 s
m-1. The corresponding mean aerodynamic resistance throughout the day is 220 s m-1. These
values are relatively different from the value of 495 s m-1 obtained by Baille (1994) using the
Leaf stomatal resistance (Fig. 2) is calculated with Eq. 0 using the values of the
parameters determined by Kichah et al. (2012) for a New Guinea Impatiens crop:
Since solar radiation has a strong impact on the stomatal resistance, its estimated value
at night (1100 s m-1) is far from the diurnal one (250 s m-1). The difference with Baille’s
diurnal value (210 s m-1) remains low but is amplified for the nocturnal value (1540 s m-1).
The evolutions over time of the measured and estimated latent heat fluxes for the three
selected days are shown in Fig. 3: (a) D1, (b) D2 and (c) D3. At night, the calculated
transpiration is quite correct, but great overestimations are observed as soon as the sun rises.
The estimated transpiration can be twice as high as the measured one. Similar observations
were reported by Prenger et al. (2002), Lopez-Cruz et al. (2008) and Villarreal-Guerrero et al.
(2012), but with smaller differences (at most 25%). At a daily time scale, the integrated value
139
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
of the evaluated transpiration is 62.7% higher than the measured one. The following sections
of the paper present investigations which were performed to explain these differences: The
validity of simplifications and the choice of the best location of climatic parameters used in
of the two major approximations used in the PM equation was analysed: a modified
expression for ∆ 0 and Rn-a (above canopy net radiation) instead of Rn (net absorbed
radiation).
Concerning the influence of the ∆ -expression, results show that the difference
between the exact value of ∆ 0 and the estimated one (Eq. 4) does not exceed 11%, inducing
differences lower than 1% on the instantaneous estimated transpiration and lower than 0.1%
net radiation Rn-a, results indicate that the relative difference between these two values can
only be very large during the night when the above canopy radiation is low. Thus, even if the
evaluated transpiration differs from 30% to 40% at night with this approximation, its impact
on the daily integrated transpiration is not relevant (less than 4%) because transpiration rates
It can therefore be deduced that these two approximations do not really have an impact
on the model results and can therefore not explain the observed difference between measured
and estimated transpiration rates. Nevertheless, in the following stages of the study, the net
radiation absorbed by the crop instead of the above-canopy net radiation will be used to
140
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
On the basis of the conclusions of Yang (1995), the choice of the location of the
measured climatic parameters was examined since it may explain the differences between the
PM model and experiments. Focusing on the transpiration rate, the slope and the coefficient
of determination R² of the regression line between the estimated and the measured values
were calculated for a set of combinations of leaf temperature and air temperature/humidity.
Results are reported in Table 3 for day D1. They clearly show that the choice of the air
climatic parameter location is more crucial than the leaf temperature to obtain reliable results
with the PM equation. In particular, the use of air humidity and temperature above the crop is
inadequate to estimate crop transpiration (slope close to 1.5). The choice of Ta3 & HR3 or Ta2
& HR2 provides similar and better correlations. For that reason, averaging both the within-
canopy air characteristics and leaf temperatures appears to be a good compromise to serve as
reference parameters to be used for the PM application, which will hereafter be referred to as
the Penman-Monteith-Modified (PMM) model. This result is consistent with New Guinea
Impatiens metabolism. Indeed, since Impatiens stomata lie on the underside of the leaves
(hypostomatic plant) and due to the fact sensible temperature and humidity gradients exist
inside the canopy, it is logical to estimate the water vapour heat exchanges by considering the
The results obtained with the PMM and the DM models are compared in this section.
For both models, the average within-canopy air temperature (Ta2+Ta3)/2 and humidity
Tl -average=(Tl1+Tl2+Tl3)/3.
141
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
Fig. 4 compares the temporal evolution of the estimated latent heat fluxes due to transpiration
with the measured one throughout day D1. Both models are globally in agreement with the
experiments, but with weak overestimations from midday to night. Similar observations were
reported for the two other days selected D2 and D3 (results not shown). This can be explained
by the position of the radiation sensor (15cm above the plants on the tray) which is smaller
than the balance tray and is therefore later impacted by the shadows of the next greenhouse
compartment. A sudden jump is observed at 2 pm: it corresponds to the end of the shadow
impact of the shading screen folded on the plants placed on the balance.
In order to analyse the accuracy of the models, the estimated values of the latent heat
flux were plotted against the measured ones for the three selected days (Fig. 5). The slope of
the linear regression, the correlation coefficient R² and the value of the model
with the PM model are also added in Table 4. The slope of the linear regression confirms the
benefit of using the inside crop climatic data for both the PMM and the DM models since
slopes are close to 1. The values of the integrated simulated latent heat flux throughout the
day and the measured one are in agreement. They reveal the ability of the two models to
evaluate the daily transpiration. The present analysis also reveals that the DM model provides
slightly better results than the PMM model (the slope of the linear regression is closer to 1),
but the R² coefficient reveals more dispersed values. However, the differences remain tenuous
as shown by the ANOVA analysis of the differences between measured and evaluated
transpiration obtained for the three selected days with PMM and DM model respectively
(Table 5). So it would be hasty to draw conclusions about the better performance of one
142
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
4. Discussion
Using the inside air parameters to evaluate transpiration not only leads to more reliable
results but also raises questions about the extension of the transpiration models to a crop
divided into several layers. Such an application could improve our knowledge of the vertical
distribution of crop transpiration and provide a more accurate analysis of the transpiration
process. Moreover, since the PM model was developed according to the Big Leaf theory, the
question of its extension to a crop divided into several layers seems relevant. In this section,
the ability of the PMM and DM model to evaluate the transpiration of a crop divided into two
layers is evaluated, and the contributions of the upper and the lower mid-part of the canopy to
the transpiration are described, as well as its evolution throughout the day. The study
considers only two layers as a first step of an innovative method and also because of
experimental constraints due to the height of the canopy (only 30 cm). The climatic
The latent heat flux due to crop transpiration is estimated for each layer by applying Eq. (1)
and Eq. (2). The daily time evolution of the sum of the contributions of each crop layer is
compared with the measured transpiration (Fig. 6). Results obtained with the DM and PMM
models are in good agreement with measurements, which is confirmed by the comparison of
estimated values integrated over the day with the measured ones (Table 7): the latent heat flux
(integrated value, i.e. the total amount of latent heat exchanged during a day) simulated by the
PMM model is overestimated by less than 8.1% for the three selected days, whereas the DM
model underestimates the measured latent heat flux by 11% at most. Both present weak
dispersion values (R²>0.90, except for day D3 with R²=0.82). Both the DM and PMM models
143
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
therefore show their ability to two layers and applying the transpiration models leads to
The results obtained with the two-layer approach are first compared to those obtained with the
one-layer method. Dividing the crop into qualitatively good results. The two-layer method is
more accurate in terms of the latent heat flux estimation than the one-layer approach (Table 4
and Table 7). In particular, the PMM model with two layers only overestimates transpiration
by 8.5% at most, whereas overestimation is greater than 12.7% with the one-layer approach.
The statistical test on the differences between evaluated and calculated transpirations for the
144 values of each selected day shows that there are no significant differences between the
performances of the two models when considering two layers (Table 8). The same statistical
test also shows that there is also no significant difference between the one and two layer
approaches (Table 9 and 10). That means that it probably would not be relevant to multiply
sensors inside the cover to obtain only a small improvement of the transpiration estimate.
The two-layer approach, however, provides useful information on the contribution of each
part of the crop. Fig. 7 shows the time evolution of the contribution of each layer to the total
transpiration for the DM and PMM models. The results from both models are consistent
throughout the day, and the contribution of the upper part is about 80% of the total
At night, the PMM model reveals a lower predominance of the upper layer (60%), whereas
according to the DM model, its contribution (75%) barely decreases under the diurnal
contribution. In fact, since the temperature and humidity gradients are slight during the night,
144
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
observed the difference is due to the distribution of the leaf area index (69% for the upper
5. Conclusion
to a New Guinea Impatiens crop grown inside a greenhouse. The Penman-Monteith model
was first tested and compared to measurements and to the Direct Method. Both models were
applied with attention to the genericity of the procedure, especially in the evaluation of the
aerodynamic resistance to water transfers. In its usual form, i.e., considering the climatic
conditions just above the crop, the PM model produces highly overestimated transpiration
rates. Questioning the assumptions of the PM approach, the present work demonstrates that
choosing the within canopy climatic parameters instead of above canopy ones, the model
leads to good agreement with measurements with a quality of results comparable to that
obtained by the Direct Method. It is then concluded that the PM model could be successfully
applied under greenhouse conditions, provided that some adjustments are made (thus defining
The study then explores the extension of the PM-Modified and DM models to analyse
the vertical distribution of the transpiration inside the crop. To do this, the crop was divided
into two horizontal layers and the parameters to be used in the model implementation were
adapted. This approach barely improved the accuracy of the transpiration estimate, meaning
that it would not be relevant in practice since it requires quite heavy additional experimental
means to obtain the required parameters. However, it does provide new information on the
vertical distribution of transpiration that could, at least partially, make up for the lack of
145
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
Acknowledgments
This work was supported by ADEME (Environment and Energy Management Agency) and
by the Region Pays de la Loire. We thank D. Lemesle, and A. Travers for their valuable
List of figures
Fig. 3: Evolution over time of the transpiration for each of the three selected days;
Fig. 4: Comparison between the evolutions over time of the transpiration evaluated
by the DM and the PMM models with the measurements for day D1.
Fig. 5: Transpiration rates evaluated by DM (a) and PMM (b) models versus
measurements.
Fig. 6: Comparison between the evolutions over time of the transpiration evaluated
by the two-layer DM and PMM models with the measurements for day D1.
during D1.
List of tables
coefficient of determination (for each study, the best model is given in bold type)
146
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
Tab. 3: Correlation between the PM estimated and measured transpiration rate for
Tab. 4: Correlations between measured and estimated transpiration rates with PM,
Tab. 5: Analysis of variance of the differences between measured and evaluated transpiration
with PMM and DM models (1 layer).
Tab. 6: Synthesis of the climatic parameters used to evaluate transpiration of the crop divided
into two layers.
Tab. 7: Correlations between measured and evaluated transpiration with transpiration models
applied to a crop with a two-layer approach.
Tab. 8: Analysis of variance of the differences between measured and evaluated transpiration
with PMM and DM models (2 layers).
Tab. 9: Analysis of variance of the differences between measured and evaluated transpiration
with 1-layer and 2-layer PMM models.
Tab. 10: Analysis of variance of the differences between measured and evaluated transpiration
with 1-layer and 2-layer DM models.
147
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
1000
900
800
700
600
500
400
300
200
100
0
0 2 4 6 8 10 12 14 16 18 20 22 0
Hour of the day
Fig. 2: Time evolution of aerodynamic and stomatal resistances (Day D1)
D1 D2 D3
550 550 550
PM PM PM
500 500 500 Measurements
Measurements Measurements
450 450 450
Transpiration (W/m²)
Transpiration (W/m²)
Transpiration (W/m²)
148
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
350
PM-Modified
DM
300
Measurements
Transpiration (W/m²)
250
200
150
100
50
0
0 2 4 6 8 10 12 14 16 18 20 22 0
Hour of the day
Fig. 4: Comparison between the evolutions over time of the transpiration evaluated by
the DM and the PMM models with the measurements for day D1.
350 350
Calculated Transpiration (W/m²)
300 300
250 250
200 200
150 150
100 100
D1 D1
D2 D2
50 50
D3 D3
0 0
0 50 100 150 200 250 300 350 0 50 100 150 200 250 300 350
Measurements (W/m²) Measurements (W/m²)
(a) (b)
Fig. 5: Transpiration rates evaluated by DM (a) and PMM (b) models versus
measurements.
350
PMM 2 layers
300 DM 2 layers
Measurements
Transpiration (W/m²)
250
200
150
100
50
0
0 2 4 6 8 10 12 14 16 18 20 22 0
Hour of the day
149
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
Fig. 6: Comparison between the evolutions over time of the transpiration evaluated by
the two-layer DM and PMM models with the measurements for day D1.
Fig. 7: Evolution over time of the transpiration contribution of each layer (proportion of
the total transpiration) evaluated by two-layer DM and PMM models for day D1.
150
Model Exp/Num Tl Ta RHa Radiation Air speed
3 levels (4 thermocouples
Boulard et al., 1991 DM/PM Exp Inside the greenhouse Inside the greenhouse Rg inside greenhouse Inside greenhouse
per level)
Boulard et Jemaa, 1993 PM/SPM Exp Not measured Inside the greenhouse Inside the greenhouse Rg inside greenhouse Not measured
Baille et al., 1994a PM Exp Not measured Above the crop Above the crop Rg above the crop Not measured
Baille et al., 1994b PM/SPM Exp Not measured Above the crop Above the crop Rn Rg above the crop Not measured
Randomly distributed Above and inside the Above and inside the
Yang, 1995 DM Exp Above and inside the crop PAR above the crop
positions and averaged crop crop
Lorenzo et al., 1998 SPM Exp Not measured Inside the greenhouse Inside the greenhouse Rg inside greenhouse not measured
Rana et al., 1998 PM Exp Not measured Above the crop Above the crop Rn above the crop Above the crop
Denden et Lemeur, 1999 PM Exp By radiometers Above the crop Above the crop Rn, Rg above the crop Above the crop
Kittas et al., 2000 DM/PM Exp 8 positions and averaged Inside the greenhouse Inside the greenhouse Rn inside the greenhouse Not measured
Leonardi et al., 2000 DM Exp Not measured Inside the greenhouse Inside the greenhouse Rg inside the greenhouse Not measured
Katsoulas et al., 2001a SPM/ DM Exp 8 positions and averaged Inside the greenhouse Inside the greenhouse Rg inside the greenhouse Not measured
6 positions per each of 3
Montero et al., 2001 SPM/DM Exp Inside the greenhouse Inside greenhouse Rn, Rg above the crop Not measured
levels
Kirnak et al., 2002 PM Exp 3 levels Near the crop Near the crop Rn near the crop Not measured
Inside the greenhouse Inside greenhouse Inside the greenhouse
Boulard et Wang, 2002 DM Exp & Num Not measured Rg (32 positions)
(24 positions) (24 positions) (24 positions)
Yes, without further
Fatnassi et al., 2003 DM Num Above the crop Above the crop Rn above the crop Not measured
information
Bartzanas et al, 2004 DM Num Not measured Inside the greenhouse Inside the greenhouse Rg inside the greenhouse Inside the greenhouse
Above and inside the
Boulard et al., 2004 DM Exp 3 leaves at the same height Above and inside the crop Rg above the crop Above the crop
crop
Above and inside the
Fatnassi et al., 2004 DM Num 3 leaves at the same height Above and inside the crop Rg, Rn above the crop Above the crop
crop
Medrano et al., 2004 SPM Exp Not measured Inside the greenhouse Inside greenhouse Rg inside the greenhouse Not measured
Medrano et al., 2005 SPM Exp Not measured Inside the greenhouse Inside greenhouse Rg above the crop Not measured
Fatnassi et al., 2006 DM Exp Not measured not measured not measured Not measured Not measured
Above and inside the Above and inside the
Demrati et al., 2007 PM Exp 3 levels Above and inside the crop Rn above the crop
crop crop
Gasquez et al., 2008 SPM Exp Not measured Inside the greenhouse Inside the greenhouse Rg inside the greenhouse Not measured
Medrano et al., 2008 SPM Exp Not measured Inside the greenhouse Inside the greenhouse Rg inside the greenhouse Not measured
Colla et al., 2009 SPM Exp Not measured Inside the greenhouse Inside the greenhouse Rg inside the greenhouse Not measured
Yes, without further
Majdoubi et al., 2009 DM Exp & Num Inside the greenhouse Inside the greenhouse Rn above the crop Not measured
information
Above and inside the Above (Rn/Rg) and inside
Kichah et al., 2012 DM Num 3 levels Above and inside the crop Above the crop
crop (PAR) the crop
Crop solarimeter above the Above and inside the
Villareal-Guerrero et al., 2012 DM/PM Exp Above and inside the crop Rg above the crop Above the crop
crop crop
Tab. 1: Location of the climatic parameters used in studies dealing with transpiration evaluation in greenhouses.
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
152
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
Model
Slope of linear
Model Day R² underestimation/overestimation
regression
(integrated values)
D1 1.68 0.95 +62.7%
PM D2 1.69 0.90 +73.3%
D3 1.79 0.92 +78.4%
D1 1.05 0.91 +7.8%
DM D2 0.94 0.90 -4.8%
D3 1.06 0.78 +11.1%
D1 1.14 0.93 +12.7%
PMM D2 1.12 0.92 +15.0%
D3 1.13 0.91 +16.5%
Tab. 4: Correlations between measured and estimated transpiration rates with PM, DM
and PMM for the three selected days.
153
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
Appendix A
The leaf aerodynamic resistance is the resistance to the water vapour and heat transfers
induced by the boundary layer that develops along the leaf due to the air movement:
ρa ⋅ c p
ra = ⋅
hs
where hs is the heat exchange coefficient expressed as a function of the Nusselt number Nu
according to the flat plate theory (Fig. ). This method requires the characteristic dimension d
(m) of the leaf, which is evaluated from the expressions suggested by Montero (2001):
2
d=
1 1
+
l w
where l and w are the length (m) and the width (m) of the leaf, respectively.
154
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
Determination of heat
exchange coefficient hs
hs ⋅ d ρa ⋅ V ⋅ d g ⋅ β ⋅ ∆T ⋅ d 3 ρ a
2
Nu = ⋅ Re = Gr =
λa µa µa 2
where Nu is the Nusselt number, hs (W m-2 K-1) is the heat exchange coefficient, d (m) is the
characteristic dimension of the flat plate, λa (W m-1 K-1) is the air thermal conductivity, Re is
the Reynolds number, ρ a (kg m-3) is the air density, V (m s-1) is the air speed, µa (Pa s) is the
air dynamic viscosity, g (m s-2) is the acceleration of gravity, β (K-1) is the volumetric thermal
expansion coefficient and ∆T (K) is the temperature difference between the flat plate and the
air.
155
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
Laminar 1
Nu = 0.54 ⋅ Gr 4
10 < Gr < 10
4 9
Free
Re 2 << Gr
Turbulent 1
Nu = 0.12 ⋅ Gr 3
10 < Gr < 10
9 12
Laminar 1
Nu = 0.68 ⋅ Re 2 + Gr 4
3 3 3
Re 2 ≈ Gr Turbulent 1
Nu = 0.03 ⋅ Re 5 + 12.1 ⋅ Gr
12 3
10 < Gr < 10
9 12
Laminar 1
Nu = 0.56 ⋅ Re 2
Forced Re < 3 ⋅ 10 5
Re 2 >> Gr Turbulent 4
Nu = 0.03 ⋅ Re 5
Re > 5 ⋅ 10 5
Tab. A.1: Expressions of the Nusselt number for the flat plate
Bibliography
Baille, M., Baille, A., Delmon, D., 1994a. Microclimate and transpiration of greenhouse rose
crops. Agricultural and Forest Meteorology 71, 83-97.
Baille, M., Baille, A., Laury, J.C., 1994b. A simplified model for predicting
evapotranspiration rate of 9 ornamental species vs climate factor and leaf area. Sci. Hortic.
Amsterdam 59, 217-232.
Bartzanas, T., Boulard, T., Kittas, C., 2004. Effect of vent arrangement on windward
ventilation of a tunnel greenhouse. Biosystems Engineering 88, 479-490.
Boulard, T., Baille, A., Mermier, M., Villette, F., 1991. Mesures et modélisation de la
résistance stomatique foliaire et de la transpiration d'un couvert de tomates de serre.
Agronomie 11, 259-274.
156
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
Boulard, T., Fatnassi, H., Roy, J.C., Lagier, J., Fargues, J., Smits, N., Rougier, M.,
Jeannequin, B., 2004. Effect of greenhouse ventilation on humidity of inside air and in
leaf boundary-layer. Agricultural and Forest Meteorology 125, 225-239.
Boulard, T., Jemaa, R., 1993. Greenhouse tomato crop transpiration model application to
irrigation control. Acta Horticulturae 335, 381-388.
Boulard, T., Wang, S., 2002. Radiactive and convective heterogeneity in a plastic tunnel:
consequences on crop transpiration. Plasticulture 121, 23-35.
Colla, G., Rouphael, Y., Cardarelli, M., Rinaldi, S., Rea, E., 2009. Water use efficiency and
tranpiration of greenhouse geranium crop. Acta Horticulturae 807, 271-276.
Demrati, H., Boulard, T., Fatnassi, H., Bekkaoui, A., Majdoubi, H., Elattir, H., Bouirden, L.,
2007. Microclimate and transpiration of a greenhouse banana crop. Biosystems
Engineering 98, 66-78.
Fatnassi, H., 2001. Modélisation et caractérisation du microclimat et de l'hétérogénéité
climatique dans une serre de grande surface équipée de filets anti-insectes. Université IBN
ZOHR. Maroc, Maroc.
Fatnassi, H., Boulard, T., Bouirden, L., 2003. Simulation of climatic conditions in full-scale
greenhouse fitted with insect-proof screens. Agricultural and Forest Meteorology 118, 97-
111.
Fatnassi, H., Boulard, T., Lagier, J., 2004. Simple indirect estimation of ventilation and crop
transpiration rates in a greenhouse. Biosystems Engineering 88, 467-478.
Fatnassi, H., Boulard, T., Poncet, C., Chave, M., 2006. Optimisation of greenhouse insect
screening with computational fluid dynamics. Biosystems Engineering 93, 301-312.
Fynn, R.P., Al Shooshan, A., Short, T.H., McMahon, R.W., 1993. Evapotranspiration
measurements ans modeling for a potted chrysanthemum crop. Transaction of the ASAE
36, 1907-1913.
Gazquez, J.C., Lopez, J.C., Baeza, E., Perez-Parra, J.J., Fernandez, M.D., Baille, A.,
Gonzalez-Real, M., 2008. Effects of Vapour Pressure Deficit and Radiation on the
Transpiration Rate of a Greenhouse Sweet Pepper Crop. Proceedings of the International
Workshop on Greenhouse Environmental Control and Crop Production in Semi-Arid
Regions, 259-265.
157
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
Katsoulas, N., Baille, A., Kittas, C., 2001a. Effect of misting on transpiration and
conductances of a greenhouse rose canopy. Agricultural and Forest Meteorology 106,
233-247.
Katsoulas, N., Kittas, C., Baille, A., 2001b. Estimating transpiration rate and canopy
resistance of a rose crop in a fan-ventilated greenhouse. Proceedings of the International
Symposium on Growing Media and Hydroponics, 303-309.
Kichah, A., Bournet, P.E., Migeon, C., Boulard, T., 2012. Measurement and CFD simulation
of microclimate characteristics and transpiration of an Impatiens pot plant crop in a
greenhouse. Biosystems Engineering 112, 22-34.
Kirnak, H., Hansen, R., Keener, H., Short, T., 2002. An evaluation of physically based and
empirically determined evapotranspiration models for nursery plants. Turkish journal of
agriculture and forestry 26, 355-362.
Kittas, C., Katsoulas, N., Baille, A., 2000. Influence of misting on the diurnal hysteresis of
canopy transpiration rate and conductance in a rose greenhouse. Proceedings of the
International Conference and British-Israeli Workshop on Greenhouse Techniques
Towards the 3rd Millennium, 155-161.
Leonardi, C., Baille, A., Guichard, S., 2000. Predicting transpiration of shaded and non-
shaded tomato fruits under greenhouse environments. Sci. Hortic. 84, 297-307.
Lopez-Cruz, I.L., Olivera-Lopez, M., Herrera-Ruiz, G., 2008. Simulation of Greenhouse
Tomato Crop Transpiration by Two Theoretical Models. Proceedings of the International
Workshop on Greenhouse Environmental Control and Crop Production in Semi-Arid
Regions, 145-150.
Lorenzo, P., Medrano, E., Sanchez-Guerrero, M.C., 1998. Greenhouse crop transpiration: An
implement to soilless irrigation management, In: Carpena, E.M. (Ed.), International
Symposium on Water Quality and Quantity in Greenhouse Horticulture, pp. 113-119.
Majdoubi, H., Boulard, T., Fatnassi, H., Bouirden, L., 2009. Airflow and microclimate
patterns in a one-hectare Canary type greenhouse: An experimental and CFD assisted
study. Agricultural and Forest Meteorology 149, 1050-1062.
Medrano, E., Alonso, F.J., Cruz Sanchez-Guerrero, M., Lorenzo, P., 2008. Incorporation of a
Model to Predict Crop Transpiration in a Commercial Irrigation Equipment as a Control
Method for Water Supply to Soilless Horticultural Crops. Proceedings of the International
158
Application des modèles de
Chapitre 6 :
transpiration : bilans hydriques
159
Résultats des simulations
Chapitre 7 :
numériques
I. Mise en œuvre
161
Résultats des simulations
Chapitre 7 :
numériques
…
Champ de vitesse, de température,
d’humidité spécifique,…
L’énergie
Flux radiatif absorbé
La masse
Flux de transpiration
Flux de condensation
Termes source /
Termes puits
162
Résultats des simulations
Chapitre 7 :
numériques
Les propriétés physiques et optiques de ces matériaux sont renseignées dans le Tableau 11.
Seules les propriétés optiques des matériaux semi-transparents sont nécessaires, les autres
propriétés optiques des autres matériaux sont prises en compte dans la définition de
l’émissivité dans le choix des conditions limites. Ainsi le coefficient d'absorption et l'indice
de réfraction du béton et de l'aluminium n'ont pas besoin d'être renseignés. On notera que la
chaleur spécifique de la vapeur d'eau est calculée par le logiciel à l'aide d'une fonction
polynomiale de la température :
C p = 1563 + 1.6 ⋅ T − 2.9 ⋅ 10 −3 ⋅ T 2 + 3.2 ⋅ 10 −6 ⋅ T 3 − 1.15 ⋅ 10 −9 ⋅ T 4 (74)
Les propriétés physiques du béton sont relativement variables selon le type de béton. Ne
connaissant pas précisément le type de béton constituant le sol de la serre, un des enjeux des
simulations numériques est l'identification de ces propriétés. Comme expliqué dans la partie
« formalisation mathématique », la propriété qui caractérise les transferts thermiques par
λ
conduction dans le sol est la diffusivité thermique a = . La masse volumique et la
ρ ⋅Cp
capacité calorifique retenue sont donc celles indiquées dans le tableau 11. Quant à la
conductivité thermique, la valeur indiquée est celle déduite des mesures et simulation via une
méthode détaillée dans ce chapitre.
Coefficient Coefficient
Capacité Conductivité Indice de
Masse d'absorption d'absorption
calorifique : thermique : réfraction
Volumique : (Courte (Grande
Cp λ (sans
ρ (kg.m-3) -1 -1 -1 -1
longueur longueur
(J.K .kg ) (W.m .K ) unité) :
d'onde) : d'onde) :
N
κ (m-1) κ (m-1)
Verre 2500 800 1.15 30 10000 1.4
Béton 2300 880 3.76 - - -
Aluminium 2719 871 202.4 - - -
Air 1.22 1006 0.0242 0 0 1
Vapeur f(T)
0.554 0.0261 0.54 0.54 1
d'eau (cf. Eq 74)
Tableau 11 : Propriétés physiques et optiques des matériaux.
163
Résultats des simulations
Chapitre 7 :
numériques
Les propriétés physiques du mélange air-vapeur d'eau sont calculées par le logiciel Fluent en
fonction de sa composition. En ce qui concerne les propriétés optiques, le coefficient
d'absorption est égal à la somme des coefficients d'absorption de chaque constituant du
mélange pondéré par des facteurs fonction de la constitution du mélange (Fluent 2012). L'air
représentant une très grande proportion du mélange et le coefficient d'absorption de la vapeur
d'eau étant faible, le mélange air-vapeur d'eau sera presque totalement transparent.
Les propriétés optiques du mélange air-vapeur d'eau seront modifiées pour le volume
correspondant au couvert végétal afin de prendre en compte la décroissance du rayonnement
dans le végétal. L'identification du coefficient d'absorption à considérer fait également l'objet
d'un développement de ce chapitre.
Equation Critère
Continuité 10-3
Vitesse 10-4
Energie cinétique turbulente k 10-5
Dissipation visqueuse Epsilon 10-5
DO 10-6
Energie 10-7
H2O 10-5
Tableau 13 : Critère de convergence sur les résidus.
4. Etapes de simulation
Les simulations numériques présentées dans la suite de ce chapitre sont organisées de la
manière suivante :
164
Résultats des simulations
Chapitre 7 :
numériques
• Dans une seconde partie, des simulations appliquées à des cas de nuit en période
hivernale (et donc en présence de chauffage) sont présentées. L'étude de ce type de situation a
fait jusqu’à présent l'objet de relativement peu de travaux numériques (Montero et al., 2005,
Piscia et al., 2012). Un élément novateur de la présente étude est la prise en compte la
modélisation de la condensation via un sous modèle adapté.
165
Résultats des simulations
Chapitre 7 :
numériques
Les deux parties pré-citées : l'application à un cas de nuit en période hivernale et l'application
à un cas de jour au printemps feront l'objet de bilans thermiques et hydriques dans une
dernière partie de ce chapitre.
1. Objectifs
Modéliser les transferts de chaleur dans le sol n’est pas indispensable si on travaille sur de
courtes durées (de l’ordre de l’heure). En revanche pour des travaux réalisés à l’échelle de la
journée, il est indispensable de considérer l’effet stockage-restitution de chaleur. Cette
propension à emmagasiner de la chaleur et à la restituer dépend de la capacité calorifique du
sol : plus elle sera élevée plus la capacité de stockage de chaleur sera importante.
Il est fréquent qu’en demi-saison (début de printemps et automne), l’énergie captée par la
serre en journée soit trop importante si bien qu’il est nécessaire de ventiler la serre pour
dissiper le surplus d’énergie (au risque sinon d’endommager les cultures). Inversement pour
la nuit, quand les températures sont trop faibles, le recours à divers systèmes de chauffage est
alors nécessaire. Cette situation assez paradoxale met en lumière les économies d’énergie qui
pourraient être réalisées grâce à une meilleure gestion de l’énergie stockée le jour dans le sol.
Les études de CFD considérant le couplage entre le sol et le climat de la serre sont peu
nombreuses. Nebbali (2008) a proposé les premiers travaux incluant l’impact du phénomène
de stockage-restitution sous serre.
Une approche instationnaire nécessite aussi la prise en compte de la course du soleil et donc
de la variation de l'aspect directionnel du rayonnement solaire direct au cours de la journée.
Ceci a une incidence sur la part du rayonnement atteignant le sol - et donc sur la quantité
d'énergie stockée dans le sol pendant la journée - mais aussi sur la part de rayonnement
absorbée par le végétal. En effet, l'atténuation du rayonnement dans un milieu semi-
transparent dépend de la longueur du chemin parcouru par un rayon lumineux dans ce milieu.
Cette variation doit être prise en compte.
Comme nous l’avons indiqué dans la partie « démarche et méthode », les différents
mécanismes physiques qui coexistent sous serre sont fortement couplés si bien qu’il est
complexe de modéliser et valider l’impact des transferts thermiques dans le sol ou celui des
transferts radiatifs dans le couvert végétal en les considérant indépendamment de l'impact des
166
Résultats des simulations
Chapitre 7 :
numériques
autres mécanismes physiques. On se propose donc dans une première étape de mener des
simulations dans des configurations simplifiées pour ensuite augmenter le degré de
complexité et à terme intégrer dans le modèle l’ensemble des mécanismes physiques.
On se propose donc dans cette partie préliminaire de modéliser les transferts thermiques dans
un sol en béton en régime instationnaire et de comparer les résultats aux solutions analytiques
ainsi que de simuler la décroissance du rayonnement dans un couvert végétal.
167
Résultats des simulations
Chapitre 7 :
numériques
Des conditions limites adiabatiques sont imposées sur les deux limites latérales et sur la limite
inférieure du domaine de calcul. A la limite supérieure, puisqu'on peut considérer que le sol
est soumis à différentes sollicitations thermiques (flux convectif ou température imposée et
flux radiatif), trois types de conditions aux limites sont envisagés :
- une température constante,
- une température variant de manière sinusoïdale en fonction du temps,
- un flux de chaleur variant de manière sinusoïdale en fonction du temps.
Pour chacune des situations, trois points ont été choisis respectivement à 3.5, 12.5 et 25 cm de
profondeur et l’évolution temporelle de la température à ces points a été comparée aux
solutions analytiques. Dans tous les cas, le sol est initialement à une température de 270K sur
toute la profondeur.
Les caractéristiques thermiques retenues pour le béton sont les suivantes :
- masse volumique : ρ s = 2300kg .m −3
168
Résultats des simulations
Chapitre 7 :
numériques
λs
La diffusivité thermique est donc a = = 1.87 ⋅ 10 − 6 m².s −1 = 6.68 ⋅ 10 −3 m².h −1 .
ρsC p s
Le pas de temps utilisé pour les simulations numériques est de 6 minutes et la simulation est
menée sur 24h, soit 240 pas de temps.
T ( z , t ) = Ti + (T0 − Ti )× erf z
(75)
2 ⋅ a ⋅ t
169
Résultats des simulations
Chapitre 7 :
numériques
Figure 73 : Evolution temporelle de la température d’un sol dont la surface est soumise
à une température constante.
Avec T0 = 290 K et α 0 = 10 K
La solution analytique pour cette configuration est :
ω 0 ω 0
T ( z , t ) = Ti + (T0 − Ti ) × erf
z
+ α 0 exp − ⋅ z ⋅ cos ω 0 ⋅ t − ⋅z (76)
2⋅ a ⋅t 2 a 2 a
Figure 74 : Evolution temporelle de la température d’un sol dont la surface est soumise
à une température sinusoïdale.
170
Résultats des simulations
Chapitre 7 :
numériques
∞
2 ⋅ϕ0 t
T ( z , t ) = T0 + ∫ erfc(u ) ⋅ du
λρc p z
2⋅ a⋅t
(77)
ϕ0 a ω 0 ω0 π
+ exp − ⋅ z ⋅ cos ω 0 ⋅ t −
⋅z−
λ ω0 2a 2a 4
Figure 75 : Evolution temporelle de la température d’un sol dont la surface est soumise
à un flux sinusoïdal.
Pour les 3 configurations, les simulations (Figures 73 à 75) restituent très fidèlement les
variations de température prédites par les solutions analytiques.
Il s'agit à présent d'utiliser les données expérimentales pour identifier la conductivité
thermique du sol.
171
Résultats des simulations
Chapitre 7 :
numériques
- La seconde consiste à comparer la valeur du flux mesurée par les fluxmètres à la surface
du sol avec la valeur estimée de ce flux à l'aide des mesures de température en profondeur
(campagne hiver 2011).
La température mesurée à la surface du sol lors des campagnes de mesures du printemps 2010
est utilisée pour modéliser une condition à la limite supérieure du domaine sous la forme
Te = T0 + α 0 ⋅ sin (ω 0 ⋅ t ) . La difficulté est que la température mesurée à la surface du sol sous
les tablettes ne suit pas la même évolution que celle mesurée en dehors des tablettes.
En supposant que le sol est soumis à une condition de température imposée à sa limite
supérieure, la distribution de température, en régime permanent, suit la fonction :
ω 0 ω 0
T ( z , t ) = To + α 0 exp − ⋅ z ⋅ cos ω 0 ⋅ t − ⋅z (78)
2 a 2 a
En utilisant les mesures de température sous les tablettes, endroit où le sol n'est jamais soumis
au rayonnement solaire direct et donc où l'hypothèse de température imposée à la limite
ω 0
supérieure est a priori valide, on évalue les valeurs de To , de α 0 exp − ⋅ z et du
2 a
ω0
déphasage − ⋅ z (Tableau 14) :
2a
1 cm 3 cm 7 cm
To 295.6 295.69 295.68
ω0
− ⋅z 0 -0.06 -0.26
2a
ω0
α 0 exp − ⋅ z -1.75 -1.51 -1.33
2a
172
Résultats des simulations
Chapitre 7 :
numériques
ω0
On en déduit deux valeurs de : 4.37 et 4.43. La diffusivité thermique est déduite à partir
2a
2π
de la valeur moyenne et en utilisant ω 0 = avec Tsignal , la période du signal égale à un
Tsignal
La comparaison entre la mesure du flux conductif par les fluxmètres utilisés lors de la
campagne de mesure de l'hiver 2011 et une estimation de ce flux à partir des températures
∂T
mesurées dans le sol ( λ ) donne accès à une estimation de la conductivité thermique. Afin
∂z
d'être certain que la comparaison n'est pas biaisée par le rayonnement solaire atteignant le sol,
la comparaison est restreinte aux mesures nocturnes. En utilisant les mesures réalisées dans le
sol sous les tablettes, la valeur moyenne de la conductivité thermique estimée est de 2.25
W.m-1.K-1. En utilisant les mesures réalisées dans le sol en dehors des tablettes, la valeur de la
conductivité thermique estimée est de 2.54 W.m-1.K-1. La valeur moyenne de ces deux valeurs
sera retenue à savoir 2.39 W.m-1.K-1.
Cette valeur correspond à une diffusivité thermique a = 1.183 ⋅ 10 −6 m².s-1. Cette valeur est
différente de la valeur obtenue à partir des mesures des campagnes du printemps 2010. On
teste donc l'impact du choix de cette valeur sur les températures simulées dans le sol.
La Figure 76 présente l'évolution de la température dans le sol à une profondeur de 5 cm.
173
Résultats des simulations
Chapitre 7 :
numériques
La différence de température induite par le choix de la diffusivité entre les deux valeurs
estimée précédemment est inférieure à 0.1K, ce qui est inférieur à l'incertitude de la mesure.
L'impact du choix est donc très faible et la première valeur estimée est retenue.
2.4. Détermination de la condition initiale
La conduite des simulations numériques en régime instationnaire requiert la connaissance de
la distribution des températures du sol au temps t = 0. En particulier, il est nécessaire de
connaître la distribution de la température en profondeur dans le sol car elle peut avoir un
impact sensible sur l'évolution du climat au cours de la journée. Cette partie a donc pour
objectif de déterminer quelle est cette distribution de température au temps t = 0. Pour cela, en
considérant comme précisé précédemment un régime permanent, l'évolution spatio-temporelle
de la température à la surface du sol est modélisée en utilisant les données expérimentales.
Puis, l'évolution de la température dans le sol est simulée sur plusieurs périodes et validée
avec les températures en profondeur dans le sol. Ainsi, l'expression analytique de la
distribution de température au temps t = 0 est retenue comme condition initiale.
On propose donc d'en tenir compte en considérant une condition limite de la forme :
Te ( x, t ) = T0 ( x ) + α 0 ( x ) ⋅ sin (ω 0 ⋅ t ) (79)
En prenant l'origine du repère en haut à gauche du domaine de calcul (x = 0 à la limite gauche
du domaine de calcul sur la Figure 71), on obtient :
- en x = 0.95 m : Te (0.9, t ) = 299.3 − 5 ⋅ sin (ω 0 ⋅ t )
En faisant l'approximation d'une symétrie par rapport à l'axe x = 4.9 (milieu de la géométrie),
on peut retenir une fonction parabolique pour les fonctions T0 (x ) et α 0 ( x ) . On obtient ainsi la
formulation suivante pour la condition à la limite supérieure :
T0 ( x ) = 0.278 ⋅ ( x − 4.9) + 294.9
2
(80)
174
Résultats des simulations
Chapitre 7 :
numériques
175
Résultats des simulations
Chapitre 7 :
numériques
176
Résultats des simulations
Chapitre 7 :
numériques
La condition initiale dans le sol est choisie en considérant que la distribution de température
découle de la variation de la température en régime permanent et que la variation à l'échelle
de la journée n'est pas influencée par la variation annuelle de la température dans le sol.
Les résultats des étapes précédentes sont donc utilisés à t = 0 pour initialiser la distribution de
la température dans le sol :
ω 0
Te ( x, z ) = T0 ( x ) + α 0 ( x ) ⋅ sin − z (82)
2 a
177
Résultats des simulations
Chapitre 7 :
numériques
178
Résultats des simulations
Chapitre 7 :
numériques
Puisque l'étude est menée pour des journées claires, on considère que l’atténuation du
rayonnement atmosphérique est négligeable devant celle du rayonnement solaire, sa
décroissance dans le couvert végétal est donc négligée. Pour cette étape préliminaire, le
modèle bi-bande est activé mais aucun rayonnement de grande longueur d'onde n'est imposé
et seule la décroissance du rayonnement de courte longueur d'onde est explorée.
La largeur du faisceau lumineux correspondant à la source du rayonnement doit être définie.
Le rayon apparent du soleil, c'est à dire l'angle solide sous lequel le soleil est vu depuis la
Terre est égal à 0.5°
3.3. Coefficient d'absorption du couvert végétal
Le coefficient d'absorption influence la décroissance du rayonnement global dans le végétal.
Théoriquement cette décroissance suit une loi de Beer Lambert : G ( z ) = G0 ⋅ e −α ⋅ z où z est le
179
Résultats des simulations
Chapitre 7 :
numériques
180
Résultats des simulations
Chapitre 7 :
numériques
181
Résultats des simulations
Chapitre 7 :
numériques
182
Résultats des simulations
Chapitre 7 :
numériques
Cette variation a une forte influence sur le rayonnement absorbé dans la couche (maille)
verticale du végétal. La Figure 87 présente le rayonnement absorbé par quelques mailles
composant le couvert végétal. Le rayonnement absorbé calculé par le logiciel de simulation
est comparé au rayonnement absorbé en utilisant l'approximation selon laquelle le coefficient
d'absorption ne varie pas durant la journée, ce qui correspond à l'approximation réalisée
jusqu'à présent dans la routine simulant l'interaction du végétal avec le climat. Les mailles
sont numérotées de haut en bas (la maille 1 est donc située dans la plus haute couche du
végétal et la maille 10 est située dans la plus basse couche du végétal). La résolution de l'ETR
et les résultats qui en découlent montrent que l'approximation sous estimait largement le
rayonnement absorbé dans les plus hautes couches du végétal (20% de la hauteur du végétal)
et le sur estimait dans les autres couches. La résolution de l'ETR entraînera donc une
distribution différente des paramètres climatiques témoin de l'activité du végétal (température
et transpiration du végétal).
L’objectif de l’étude d’un nuit en période hivernale est de simuler l’évolution du climat dans
le serre en présence de chauffage afin d’être capable de caractériser les pertes énergétiques.
Les phénomènes de condensation sont pris en compte ce qui permettra d’évaluer le moment et
le lieu où apparaît la condensation durant la nuit.
183
Résultats des simulations
Chapitre 7 :
numériques
liées aux choix des conditions limites n’influent pas sur le réalisme et les résultats obtenus. Le
domaine de calcul est donc large de 19.6 m et haut de 10 m (soit 4.1 m au dessus du faitage).
Les ouvrants sont complètement fermés. Le maillage est structuré quadratique et comporte
36960 mailles de dimension moyenne 8x6.5 cm.
Vent
Ouest Est
Compartiment 1 Compartiment 2
Zone de condensation
Paroi 2 de la toiture
Tuyau de chauffage
184
Résultats des simulations
Chapitre 7 :
numériques
près de la toiture, lieu où la condensation est susceptible d’intervenir, et au niveau des tubes
de chauffage aérien (Figure 89). Une étude de sensibilité sur la densité du maillage et son
impact sur les résultats obtenus a été menée en considérant 5 maillages différents (de 9200
mailles à 47700 mailles). Les résultats ne varient pas entre le maillage de 36960 mailles et
celui de 47700 mailles, d’où le choix retenu.
Le volume correspondant au couvert végétal est discrétisé de la même manière pour les deux
compartiments mais seule l’interaction du végétal dans le compartiment 1 (à l’ouest) est
simulée. De même, le chauffage aérien n’est considéré que dans le compartiment 1.
185
Résultats des simulations
Chapitre 7 :
numériques
a) b)
Figure 91 : Variations des conditions aux limites pour les nuits claire (a) et nuageuse (b)
(cf. Figure23 pour la position des capteurs).
3. Conditions initiales
Les conditions initiales à 0h sont fixées de la manière suivante : un cas stationnaire a été
considéré ; les conditions limites utilisées sont les mesures relevées à 23h. En faisant
converger ce cas stationnaire, on dispose alors d’un champ de vitesse d’air établi à 0h qui est
conservé pour constituer la situation initiale de l’écoulement d’air dans le compartiment.
Ensuite, certains paramètres sont modifiés avec les valeurs mesurées à 0h :
186
Résultats des simulations
Chapitre 7 :
numériques
l’humidité absolue et la température de l’air dans le volume de la serre sont fixées aux
valeurs déduites des mesures au dessus du couvert végétal ;
l’humidité absolue et la température de l’air dans le couvert végétal sont fixées aux
valeurs déduites des mesures dans le couvert végétal ;
la température de feuille est la température moyenne des 15 mesures de température de
feuille ;
les températures de toiture sont également fixées à partir de valeurs mesurées
correspondantes.
Toutes ces valeurs sont synthétisées dans le Tableau 16.
d’eau dans l’équation de transport d’espèce au niveau des mailles adjacentes aux parois où la
condensation survient (Piscia et al., 2012) :
ω h2o ∂ω h 2o A cell wall
& ′h′′2o =
m ρD
(ω h2o − 1) ∂n Vcell
Autrement dit, lorsqu’il y a condensation, une certaine quantité de vapeur d’eau est
retirée à l’air. Les effets thermiques liés au changement de phase de l’eau de l’état
vapeur à l’état liquide sont églement pris en compte en ajoutant un terme source dans
l’équation de l’énergie qui est proportionnel à la masse d’eau condensée. Ceci est
réalisé à partir d’une UDF (Annexe F) qui :- récupère la température de paroi pour la
paroi définie comme lieu potentiel de condensation et évalue, à l’aide de la pression de
187
Résultats des simulations
Chapitre 7 :
numériques
- ajoute une quantité d’énergie dans l’équation de conservation de l’energie pour les
mailles adjacentes à la paroi sur laquelle il y a de la condensation.
E cond = m
& 'h'' 2O × Lv
Tout ceci necessite le calcul des propriétés physique de l’eau (« H2O property routines »)
5. Performance du modèle
Le modèle numérique développé en prenant en compte l’ensemble des phénomènes physiques
décrits dans le chapitre 2 et formalisés dans le chapitre 3 ainsi que la condensation a été validé
sur la base de mesures réalisées lors de la nuit claire. L’évolution de l’ensemble des
paramètres climatiques simulés est présentée et comparée aux valeurs expérimentales. Le
modèle est par la suite utilisé et validé avec les conditions de la nuit nuageuse.
5.1. Température dans la serre
188
Résultats des simulations
Chapitre 7 :
numériques
189
Résultats des simulations
Chapitre 7 :
numériques
Le champ de température simulé met aussi en évidence que les températures les plus faibles
apparaissent au niveau de la toiture de la serre, justifiant le fait que la condensation apparait
essentiellement sur ces parois.
5.2. Température dans le couvert végétal
190
Résultats des simulations
Chapitre 7 :
numériques
191
Résultats des simulations
Chapitre 7 :
numériques
de vapeur d’eau dans l’air du coté est de la serre. Cette distribution est liée à l’écoulement de
l’air dans la serre et indique que ce paramètre dépend fortement du mouvement de l’air
Figure 97 : Champ d'humidité relative (%) dans la serre à 5h et zoom sur la toiture.
a)
b)
Figure 98 : Champ d'humidité absolue (kg.kg-1) dans la serre à 1h (a) et à 5h (b).
192
Résultats des simulations
Chapitre 7 :
numériques
Figure 99 : Evolution des températures de toiture. (cf Figure 23 pour la position des
capteurs)
L'écart le plus important constaté entre les simulations et les mesures concerne les
températures de toiture (entre 3.4°C et 4.3°C d'écart selon le moment de la nuit). Cependant
l'évolution au cours de la nuit (diminution de la température de 2°C environ) est la même pour
les simulations que pour la mesure (Figure 99).
Des études de sensibilité sur la vitesse de l'air à l'entrée du domaine de calcul et sur la
température de ciel ont été menées afin de quantifier leur impact sur la température de toiture.
Même en utilisant une vitesse 10 fois supérieure à celle utilisée pour les simulations (1m.s-1),
ou en abaissant la température de ciel de 50°C, la température de paroi diminue de moins d'un
degré. Ces deux paramètres n’expliquent donc pas les écarts constatés pour les températures
de parois simulées. L'hypothèse selon laquelle l'impact de la température intérieure sur la
température de toiture serait trop important est donc évoquée. Puisque les températures dans
la serre sont correctement simulées, cette hypothèse se traduirait par une surestimation du
coefficient d'échange air/paroi, ce qui recoupe à nouveau l'hypothèse d'un brassage d'air et
d'une vitesse de l'air trop importante dans la serre.
5.5. Vitesse de l'air
Avant de discuter l'hypothèse d'un brassage d'air trop important dans la serre évoqué
précédemment, la distribution des vecteurs vitesse est présentée (Figure 100). Les mesures
expérimentales faisaient penser à un écoulement sous forme d’une cellule rotative dans le sens
horaire (Chapitre 5.II.2.2). Les simulations indiquent également un écoulement dans le sens
horaire au niveau des parois avec cependant deux cellules contra rotatives qui se situent au
centre de la serre au lieu d’une seule. L’établissement de ces deux cellules s’explique par la
193
Résultats des simulations
Chapitre 7 :
numériques
présence des tubes de chauffage. Cette distribution des vitesses explique la distribution de
l’humidité absolue décrite précédemment : la principale cellule rotative contribue à la
convection de l’humidité du coté Est de la serre.
194
Résultats des simulations
Chapitre 7 :
numériques
qu’il ne l’est dans la réalité et ce d’autant plus que les tablettes de cultures n’occupent que 3m
sur les 10m de large de la serre (dans le plan perpendiculaire à celui retenu pour les
simulations)
5.6. Application au cas d’une nuit nuageuse
Le modèle développé pour une nuit claire et validé sur la base de mesures expérimentales est
à présent appliqué à une nuit nuageuse afin de vérifier qu’il peut tout aussi bien prédire le
climat sous serre avec des conditions radiatives et des conditions de chauffage très différentes.
Dans un souci de concision, seule la précision du modèle appliqué à une nuit nuageuse est
comparée à celle du modèle appliqué à une nuit claire, ce qui permet également de synthètiser
les résultats pour la nuit claire. Les évolutions des paramètres climatiques simulés et
comparés aux données expérimentales sont présentées dans l’annexe A.
Pour rendre compte de la précision du modèle, un indicateur de type RMS a été retenu :
( )
2
1 n *
RMS = ∑ yi − yi
n i =1
avec
n : le nombre de mesure
yi* : la valeur issue de la simulation
yi : la valeur issue de l’expérimentation
Les résultats obtenus pour les deux types de nuit sont consignés dans le tableau 17 :
Nuit Claire Nuit Nuageuse
T1 (°C) 2,07 1,39
T3 (°C) 0,06 0,31
T5 (°C) 1,40 1,17
Tp2 (°C) 3,88 2,92
Tf (°C) 0,29 0,39
HR5 (%) 2,95 5,81
Tableau 17 : Comparaison des RMS de la nuit claire et de la nuit nuageuse.
La précision du modèle appliqué à une nuit nuageuse est meilleure que pour la nuit claire pour
la température sous les tablettes et la température au dessus du végétal, identique pour la
température de feuille et moins bonne pour la température dans le végétal mais elle reste
cependant excellente (RMS = 0.3°C). Une meilleure corrélation simulation-mesure est
également obtenue pour les températures de toiture mais l’écart reste tout de même important
195
Résultats des simulations
Chapitre 7 :
numériques
(RMS = 2.92). Enfin, la précision est également moins élevée pour la simulation de
l’humidité relative.
5.7. Bilans hydriques et énergétiques
Ces modèles de simulations développés permettent de disposer de nombreuses informations.
Ils permettent en particulier de réaliser des bilans hydriques et énergétiques. Les premiers
fournissent des facteurs explicatifs sur l’évolution de l’humidité dans la serre, ainsi que des
ordres de grandeurs sur les flux de condensation et de transpiration. Les seconds permettent
principalement de quantifier et caractériser les pertes énergétiques de la serre au cours de la
nuit.
196
Résultats des simulations
Chapitre 7 :
numériques
Figure 103 : Evolution de l'humidité absolue moyenne dans la serre au cours de la nuit
claire.
Pour clore, cette partie bilan hydrique, il peut être intéressant de présenter la distribution de la
condensation le long de la toiture (Figure 104 et 105, l'origine de l'abscisse de cette figure
correspond à la limite ouest du compartiment). En effet, il s'avère que la condensation
apparaît essentiellement sur la toiture ouest du compartiment et que la distribution de flux
condensé sur cette paroi est différente pour la nuit claire et la nuit nuageuse. La condensation
influence naturellement la distribution de l'humidité absolue au voisinage de la paroi et
explique pourquoi elle présente les valeurs les plus faibles à cet endroit.
197
Résultats des simulations
Chapitre 7 :
numériques
Ces valeurs doivent néanmoins être relativisées du fait de l’erreur observée entre la mesure et
la simulation en ce qui concerne les températures de paroi qui est un paramètre prépondérant
dans la régulation du flux de condensation.
Les bilans énergétiques réalisés sur la serre concernent principalement la caractérisation des
déperditions thermiques au niveau de la toiture. En effet, les apports énergétiques dans la
serre sont connus et sont ceux mesurés lors des campagnes expérimentales. Ils sont rappelés
dans le tableau 18 :
Nuit Nuageuse Nuit Claire
Chauffage Chauffage Chauffage Chauffage
Total Total
Sol Aérien Sol Aérien
0h 32 1 33 28 88 116
1h 40 11 50 26 52 78
2h 41 21 62 29 55 83
3h 43 2 45 36 57 93
4h 46 0 46 36 53 89
5h 49 1 50 34 52 86
6h 49 1 50 35 56 90
Tableau 18 : Récapitulatifs des flux de chauffage (en W.m-2 de sol).
Puisque la température de la serre est maintenue à 16°C par régulation des flux de chauffage,
il n'y a pas de phénomène de stockage de chaleur dans l'air. Le flux de chauffage est ajusté
pour compenser les pertes de la serre qui ont lieu exclusivement au niveau la toiture. Comme
Ces pertes ont deux causes principales : les échanges radiatifs et les échanges convectifs avec
198
Résultats des simulations
Chapitre 7 :
numériques
l'extérieur. La contribution de chacune de ces pertes aux déperditions totales est évaluée et
consigné dans le tableau 19 :
Déperditions nuit claire (W/m²) Déperditions nuit nuageuse (W/m²)
totales Radiatives Convectives totales radiatives convectives
1h 1114 790 324 759 515 244
2h 1098 812 286 804 534 270
3h 1153 817 336 853 562 291
4h 1174 819 355 868 569 299
5h 1166 816 350 880 573 307
6h 1183 820 363 855 547 308
7h 1183 817 366 855 547 308
Tableau 19 : Quantification des pertes radiatives et convectives au niveau de la toiture.
Les pertes convectives sont du même ordre de grandeur pour les deux types de nuit et ne
dépendent que de la différence de température entre la toiture et l'air extérieur puisque les
conditions extérieures sur la vitesse de l'air sont identiques. Puisque les températures de
toiture sont surestimées dans les simulations en comparaison des mesures, cet écart de
température est en réalité certainement plus important et donc les pertes convectives données
ici peuvent être sous estimées.
Toujours est-il que les déperditions de chaleur au niveau de la toiture sont plus importantes
pour la nuit claire du fait des pertes radiatives plus importantes que pendant la nuit nuageuse
(+48% en moyenne).
Le tableau 20 présente la proportion des pertes radiatives et convectives. Aucune variation
significative n'est constatée au cours de la nuit. La principale différence a trait à la
comparaison des deux types de nuit. La part des pertes radiatives représente en moyenne 71%
des pertes durant une nuit claire et 66% durant une nuit nuageuse.
199
Résultats des simulations
Chapitre 7 :
numériques
Les simulations menées pendant le jour ont pour but principal de prendre en compte la course
du soleil. Elles intègrent donc les développements et validations exposés dans les étapes de
simulation préliminaire :
l’impact du sol : la conduction de la chaleur dans le sol et donc les phénomènes de
stockage/restitution d’énergie selon le moment de la journée.
la prise en compte de la décroissance du rayonnement à travers le modèle DO appliqué
à la végétation et donc sa variation au cours de la journée selon l’angle sous lequel le
rayonnement solaire traverse le couvert végétal.
Une journée parfaitement ensoleillée a été retenue (20 mai 2010) afin de mener ces
simulations numériques. Les ouvrants de la serre sont complètement ouverts tout au long de la
journée et les écrans d’ombrages complètement repliés. Les conditions de culture des
compartiments adjacents ne sont pas contrôlées, en particulier la position des écrans
d’ombrage, si bien que du fait d’ombres portées par ceux-ci sur le compartiment d’étude, il a
fallu restreindre la période d’étude à la plage horaire 8h-16h (heure solaire).
1. Domaine de calcul et maillage
Le domaine de calcul (Figure 106) comprend le compartiment de serre étudié et le volume
d’air au dessus de celui-ci. Du fait de contraintes techniques qui empêchent de définir des
conditions limites radiatives variables au niveau de l’entrée et de la sortie d’air (cf. IV.2), le
domaine est étendu au compartiment en amont et à celui en aval afin que le rayonnement
solaire incident sur le compartiment de serre étudié soit bien représentatif de la réalité durant
la période d’étude (8h - 16h). Le domaine de calcul est donc large de 29.4 m et haut de 10 m.
Vent
Ouest Est
200
Résultats des simulations
Chapitre 7 :
numériques
Le maillage retenu (Figure 107) est semblable au maillage utilisé pour les conditions de nuit,
structuré quadratique et comporte 51810 mailles de dimension moyenne 8.9x6.3 cm.
Comme pour les conditions de nuit, le volume correspondant au couvert végétal est discrétisé
pour les trois compartiments mais seule l’interaction avec le végétal est considérée dans le
compartiment central (compartiment 2).
201
Résultats des simulations
Chapitre 7 :
numériques
202
Résultats des simulations
Chapitre 7 :
numériques
r
S compensati on = κ λ ∫ λ
I ( s , Ω ).d Ω .d λ (83)
4π
Cette correction est réalisée via une UDF spécifique (annexe F).
4. Conditions initiales
Les conditions initiales sont pour une part fixées à partir des mesures réalisées dans la serre à
minuit (Tableau 21) et pour une autre part à partir du travail réalisé dans les étapes
préliminaires pour les températures dans le sol dont il est primordial de prendre en compte
l'inertie.
Tf (K) 287.6
Humidité absolue de l’air
0.0053
dans la serre (kg.kg-1)
Humidité absolue dans le
0.0069
végétal (kg.kg-1)
Température de toiture 275.3
Tableau 21 : Conditions initiales des simulations de jour.
5. Performance du modèle
5.1. Vitesse de l'air
Les ouvrants sont, pour cette configuration, complètement ouverts. Les échanges de masse
d'air entre l'extérieur et l'intérieur de la serre sont donc importants et les paramètres
climatiques de l'air extérieur (vitesse, température, humidité) ont une influence sensible sur le
climat dans la serre comme l'exploitation des campagnes de mesure l'a mis en exergue
(cf. Chapitre 5.III). Il est donc important de commencer par présenter la distribution des
vecteurs vitesse dans la serre (Figure 110). Conformément à ce à quoi on pouvait s'attendre
(Ould Khoua, 2006), une cellule rotative s'instaure dans le compartiment de serre. Les vitesses
(Figure 110) sont sensiblement plus élevées que pour les cas de nuit et atteignent 1 m.s-1 au
maximum dans le compartiment (au contact des parois).
203
Résultats des simulations
Chapitre 7 :
numériques
Figure 110 : Champ des vecteurs vitesse (m.s-1) de l'air simulée à 12h.
204
Résultats des simulations
Chapitre 7 :
numériques
Plage horaire
d'étude
205
Résultats des simulations
Chapitre 7 :
numériques
Plage horaire
d'étude
Figure 114 : Evolution de l’humidité de l’air simulée et mesurée dans le couvert végétal.
5.4. Température du sol
La température à la surface du sol est comparée aux mesures réalisées (Figure 115). La
température simulée est une moyenne sur toute la surface du sol du compartiment 2 et la
température mesurée à laquelle on compare les simulations est une moyenne des températures
du sol sous et en dehors des tablettes pondérée par les surfaces qu'elles représentent. La
moyenne de l'écart absolue entre mesure et simulations pour les valeurs de la plage d'étude est
de 0.3°C avec un écart maximum de 0.7°C à 10h. Le modèle restitue donc correctement
l'interaction du sol avec le climat aérien de la serre. Il reste à comparer les flux de chaleur
stockés et restitués par le sol ce qui fera l'objet d'une partie des bilans énergétiques (cf. V).
Plage horaire
d'étude
206
Résultats des simulations
Chapitre 7 :
numériques
Plage horaire
d'étude
Figure 117 : Champ de température (K) de l'air simulée dans le végétal à 12h.
207
Résultats des simulations
Chapitre 7 :
numériques
La température des feuilles les plus hautes (40 cm dans le couvert végétal) est en effet 6°C à
7°C supérieure aux feuilles les plus basses (10 cm) ce qui correspond à l'ordre de grandeur des
mesures expérimentales. La comparaison des températures moyennes de feuille (mesurées et
simulées) présentées sur la Figure 121 ne montre en effet que de faibles écarts (1°C en
moyenne) compte tenu de l'impact de la position du capteur de température sur la température
moyenne expérimentale, notamment pour les couches les plus hautes. En effet, entre 30cm et
40cm, des écarts de température simulée atteignent jusqu'à 4°C. La mesure des températures
dans cette couche, la plus haute du végétal, comptant pour 1/3 dans la température moyenne,
la position exacte du capteur peut avoir une influence importante sur la température moyenne
du végétal.
Figure 119 : Evolution des températures de feuille simulées pour différentes hauteurs
dans le végétal.
208
Résultats des simulations
Chapitre 7 :
numériques
Plage horaire
d'étude
Plage horaire
d'étude
RMS
Tair (°C) 0,28
HR (%) 3,59
Tair Végétal (°C) 1,54
HR Végétal (%) 13,95
Tf (°C) 1,20
Ts (°C) 0,40
Transpiration (W/m²) 35,36
Tableau 22 : RMS pour les simulations de jour.
Il faut retenir que le modèle utilisé restitue bien la dynamique des températures et humidités
dans la serre ainsi que pour les températures de feuille et la transpiration du couvert.
209
Résultats des simulations
Chapitre 7 :
numériques
Cependant, la température et l’humidité dans le couvert végétal ne sont pas bien simulées du
fait de l’existence de gradients trop faibles entre l'air dans et au dessus du végétal.
1. Bilans hydriques
Le modèle de climat distribué permet d’aborder finement les échanges au niveau du végétal,
en particulier de connaitre précisément la distribution spatiale de la transpiration et son
évolution au cours de la journée. L’application de modèles de transpiration à un végétal divisé
en deux couches horizontales sur la base de données expérimentales (Chapitre 6) a fourni des
résultats qui donnent des éléments de validation pour les simulations numérique. En
particulier il a été montré que pendant la nuit que la contribution de la couche supérieure du
végétal à la transpiration totale est comprise entre 60% et 70% selon le modèle de
transpiration utilisé. Cependant, compte tenu de l’existence de gradients d’humidité et de
température quasiment nuls et de la distribution de la densité de feuille, une contribution
proche de 70% était plutôt attendue. De jour, la contribution de la couche supérieure est
nettement prépondérante (entre 80% et 85%).
La Figure 122, issue des simulations numériques, présente les contributions de la moitié
supérieure et inférieure du végétal à la transpiration totale. Durant le jour, la couche
supérieure du végétal contribue pour 80% à la transpiration totale du végétal, ce qui recoupe,
malgré une légère sous estimation, les résultats obtenus à partir des modèles de transpiration
utilisant les données expérimentales. Durant la nuit, les deux couches contribuent dans la
même proportion (52%-48%). L’écart avec les résultats expérimentaux se justifie par le fait
que le végétal est considéré comme étant homogène en termes de distribution de feuille, ce
qui ne reflète pas exactement les conditions expérimentales. Cela explique également la légère
sous estimation de la contribution de la couche supérieure à la transpiration totale durant la
journée et met aussi en évidence, la nécessité d’aller vers plus de précision dans la
discrétisation du végétal.
210
Résultats des simulations
Chapitre 7 :
numériques
211
Résultats des simulations
Chapitre 7 :
numériques
Figure 123 détaille l’évolution de toutes ces composantes au cours de la journée. Ainsi on
constate qu’au midi solaire, 65% du flux solaire atteint la surface du végétal. 42% de ce
rayonnement absorbé contribue au phénomène de transpiration alors que les 58% restant
représentent la chaleur sensible échangée avec l’air pour l’ensemble du végétal (la distribution
verticale de cette chaleur sensible échangée avec l’air étant relativement hétérogène) qui
contribue au réchauffement de l’air.
La nuit, en revanche, les échanges de chaleur sensible sont inversés si bien que c’est l’air qui
cède de la chaleur au végétal.
212
Résultats des simulations
Chapitre 7 :
numériques
autres grandeurs mises en jeu. Cela n’est pas le cas de l’énergie stockée dans le sol qui
représente un volume bien plus important et induit des échanges convectifs qui
contribuent à l’augmentation ou à la diminution de la température de l’air de la serre.
Enfin l’air de la serre peut aussi faire l’objet d’un bilan énergétique : il échange de
l’énergie avec les éléments de la serre par convection et contribue aussi à faire varier sa
température. Une part sera convectée et évacuée à l’extérieur de la serre par transfert de
masse d’air, l’autre sera « conservée » à l’intérieur de la serre sous forme de capacité
calorifique. L’air ayant une masse volumique très faible, cette seconde part peut,
compte tenu des variations de température, être négligée (le flux maximal échangé
représente 7 W.m-2). La part d’énergie évacuée via la ventilation est estimée en utilisant
les valeurs de température de l’air en sortie et en entrée afin de réaliser un bilan (le débit
d’air entrant dans le volume étant égal au débit d’air sortant).
L’ « utilisation » du rayonnement solaire au niveau de la serre est synthétisée sur la
Figure 124.
Rayonnement
solaire
213
Résultats des simulations
Chapitre 7 :
numériques
214
Conclusions et perspectives
Conclusions et perspectives
Des campagnes de mesures ont été menées au printemps 2010 et à l’hiver 2011 pour
répondre aux besoins concernant les deux situations étudiées. Ces campagnes de mesures ont
pour objectif d’explorer, dans un premier temps, la distribution des paramètres climatiques
dans la serre, en particulier dans l’environnement proche du végétal afin de pouvoir mieux
comprendre l’influence des uns sur les autres. Dans un second temps, ces campagnes de
mesures permettent de fournir des données d’entrée et de validation pour les modèles
numériques développés.
Les campagnes de mesures de nuit en présence de chauffage ont mis en évidence les
besoins plus importants en chauffage et les températures de toiture plus basses lors des nuits
215
Conclusions et perspectives
claires. La distribution des paramètres climatiques dans la serre sont en revanche semblables
pour les deux types de nuit puisque du fait de la régulation en température de la serre, le
climat établi est comparable.
Les campagnes de mesure menées au printemps 2010 ont permis de mettre en
évidence l’influence prépondérante du rayonnement solaire sur la température et la
transpiration du végétal. Les gradients de température et d’humidité de l’air qui s’amplifient
dès lors que le rayonnement solaire influence le climat dans la serre ont été quantifiés, en
particulier dans l’environnement proche du végétal. Le phénomène de stockage-restitution de
chaleur dans le sol a été mis en évidence et analysé. Enfin, la sensibilité du taux d’humidité
ainsi que la moindre dépendance de la température vis-à-vis de la « qualité » de la ventilation
a été démontrée.
Les données expérimentales ont été exploitées afin de mettre en œuvre les modèles de
transpiration. Les modalités d’application de ces modèles prévus initialement pour des
cultures en condition extérieure ont été discutées et des solutions sont proposées pour les
appliquer à des végétaux cultivés sous serre. En guise de prolongement, l’application des
modèles de transpiration à un végétal divisé en 2 couches a été validée. Cette approche permet
de fournir de précieuses informations sur la distribution de la transpiration du végétal,
améliore la compréhension des mécanismes intervenant au niveau du couvert et permet
l’identification des principaux paramètres qui affectent la transpiration.
216
Conclusions et perspectives
d’un sous modèle de culture interagissant avec le climat est une approche novatrice. Les
résultats (température et humidité) obtenus avec les simulations sont en bon accord avec les
relevés expérimentaux (écarts moyen de 1°C environ) effectués in-situ. Seules les
températures de toiture présentent une surestimation importante (3°C environ). Les gradients
de température simulés plus faibles que ceux constatés expérimentalement s’expliquent par un
brassage de l’air trop important à l’intérieur de la serre. Ce brassage excessif explique
également pourquoi les températures de toiture sont surestimées puisqu’il favorise le
réchauffement de la toiture par l’air intérieur. L’exploitation de cette situation a permis
d’identifier le moment où apparait la condensation, de localiser les endroits où les risques de
condensations sont les plus importants et de mettre en évidence et de quantifier la contribution
prépondérante des pertes radiatives aux pertes thermiques de la serre, contribution d’autant
plus prépondérante que les nuits sont claires.
Les simulations appliquées à la situation diurne intègrent les développements
préliminaires suivants :
- la prise en compte des transferts de chaleur dans le sol permet de mettre en évidence et de
quantifier le phénomène de stockage-restitution de la chaleur par le sol ;
- la résolution de l’équation de transfert radiatif permet de prendre parfaitement en compte
l’influence du caractère directionnel du rayonnement solaire.
Les valeurs simulées avec le modèle surestiment en moyenne d’1°C environ les
températures. L’humidité au dessus du végétal est correctement simulée (écart de 4 points
d’humidité). En revanche les paramètres climatiques simulés dans le couvert végétal
présentent de plus grands écarts avec les mesures. Bien que les gradients de température de
feuille soient bien restitués, ceux de température et d’humidité de l’air dans et au dessus du
végétal sont trop faibles.
Les bilans énergétiques réalisés à partir de ce modèle permettent d’identifier la
manière dont est « utilisé » le rayonnement solaire pénétrant dans la serre. Ils permettent
également de quantifier l’impact du végétal sur le climat de la serre.
Les écarts entre mesures et simulations les plus importants ont été identifiés entre l’air
dans le végétal et au dessus du végétal. Les gradients trop faibles mettent en évidence le fait
que la prise en compte du végétal peut être améliorée. En particulier, il faudrait dépasser
217
Conclusions et perspectives
Une fois de telles évolutions intégrées, le modèle pourrait être utilisé et faciliter le
développement de logiciels d’aide à la décision, à la conception et/ou à la gestion/régulation
des paramètres climatiques.
218
Références bibliographiques
Références bibliographiques
Bailey, B.J., Day, W., 1999. The use of models in greenhouse environmental control. Acta
Hort. (ISHS) 491, 93-100.
Baille, M., Baille, A., Delmon, D., 1994a. Microclimate and transpiration of greenhouse rose
crops. Agricultural and Forest Meteorology 71, 83-97.
Baille, M., Baille, A., Laury, J.C., 1994b. A simplified model for predicting
evapotranspiration rate of 9 ornamental species vs climate factor and leaf area. Sci.
Hortic. Amsterdam 59, 217-232.
Baille, M., Baille, A., Tchamitchian, M., 1990. A simple model for the estimation of
greenhouse transmission: influence of structures and internal equipment. Acta
Horticulturae 281, 35-46.
Baptista, F.J., Bailey, B.J., Meneses, J.F., 2008. Comparison of humidity conditions in
unheated tomato greenhouses with different natural ventilation management and
implications for climate and botrytis cinerea control. Acta Horticulturae 801:1013-
1020
Baptista, F.J., Bailey, B.J., Randall, J.M., Meneses, J.F., 1999. Greenhouse ventilation rate:
theory and measurement with tracer gas techniques. J. Agric. Eng. Res. 72, 363-374.
Bartzanas, T., Boulard, T., Kittas, C., 2004. Effect of vent arrangement on windward
ventilation of a tunnel greenhouse. Biosystems Engineering 88, 479-490.
Bezari, S., Bouhdjar, A., Aït Messaoudène, N., 2007. Etude du microclimat d'une serre tunnel
équipée d'un dispositif de stockage thermique dans l'eau. Energies renouvelables 07,
307-313.
Bibi Triki, N., Benyarou, F., Benyoucef, B., El Sayed, M.A.R., 1999. Bilan thermique et
méthode d'estimation des besoins énergetiques de la serre agricole tunnel à double
paroi constrituée dans sa face Nord d'un mur en panneau sandwich. Revue des
Energies Renouvelables, 77-82.
Boulard, T., 1985. Etude du gisement de chaleur récupérable sur l'air d'une serre :
comparaison avec les besoins en chauffage. Station de Bioclimatologie, Montfavet.
Boulard, T., Baille, A., Mermier, M., Villette, F., 1991. Measurement of stomatal resistance
and transpiration in a greenhouse tomato canopy - Comparison between a 1 layer and
a multilayer model. Agronomie 11, 259-274.
219
Références bibliographiques
Boulard, T., Fatnassi, H., Roy, J.C., Lagier, J., Fargues, J., Smits, N., Rougier, M.,
Jeannequin, B., 2004. Effect of greenhouse ventilation on humidity of inside air and in
leaf boundary-layer. Agricultural and Forest Meteorology 125, 225-239.
Boulard, T., Wang, S., 2002. Experimental and numerical studies on the heterogeneity of crop
transpiration in a plastic tunnel. Computers and Electronics in agriculture 34, 173-190.
Bournet, P.E., Ould Khaoua, S.A., Boulard, T., 2007. Predicted effects of vents arrangements
on the ventilation and energy transfers in a multi-span glasshouse using a bi-band
radiation model. Biosystems Engineering. Volume 98, Issue 2, p. 224–234
Boussinesq, J., 1903. Théorie analytique de la chaleur mise en harmonie avec la
thermodynamique et avec la théorie mécanique de la lumière, volume 2, Gauthier-
Villars.
Britton, C.M., Dodd, J.D., 1976. Relationships of photosynthetically active radiation and
shortwave irradiance. Agricultural Meteorology 17, 1-7.
Brown, R.H., Blaser, R.E., 1968. Leaf area index ipi pasture growth. Herb. Abstr. 38, 1-9.
Chassaing, P., 2000. Turbulence en mécanique des Fluides. Analyse du phénomène en vue de
sa modélisation à l'usage de l'ingénieur, In: Cépaudès-éditions (Ed.), Institut National
Polytechnique de Toulouse.
Darcy, H., 1956. Les fontaines publiques de la ville de Dijon, Paris.
Demrati, H., Boulard, T., Fatnassi, H., Bekkaoui, A., Majdoubi, H., Elattir, H., Bouirden, L.,
2007. Microclimate and transpiration of a greenhouse banana crop. Biosystems
Engineering 98, 66-78.
Fatnassi, H., Boulard, T., Poncet, C., Chave, M., 2006. Optimisation of greenhouse insect
screening with computational fluid dynamics. Biosystems Engineering 93, 301-312.
Fleury, B., 1982. Contribution à la modélisation d'un couvert végétal associatif. Modèle
dynamique d'un couvert de trèfle violet (Trifolium pratense L.) en culture pure.
Mémoire de fin d'études E.N.I.T.A., Dijon.
Fluent, 1998. Fluent5 User's Guide. Fluent Incorporated, Lebanon, USA.
Fluent, 2012. Fluent13.1 User's Guide. Fluent Incorporated, Lebanon, USA.
Guyot, G., 1999. Climatologie de l'environnement. De la plante aux écosystèmes.Masson,
Paris (F.).
Ha-Minh, H., 1991. Physique et modélisation de la turbulence en écoulement de fluides.
Ecole de printemps de mécanique des fluides numérique.
220
Références bibliographiques
Haxaire, R., 1999. Caractérisation et modélisation des écoulements d'air dans une serre.
Université de Nice Sophia Antipolis, Nice, pp. 1-148.
Jolliet, O., 1988. Modélisation du comportement thermique d'une serre horticole. Ecole
Polytechnique Fédérale de Lausanne.
Jones, P.J., Whittle, G.E., 1992. Computational Fluid Dynamics for Building Air Flow
Prediction - Current Status and Capabilities. Building and Environment 27, 321―338.
Kichah, A., Bournet, P.E., Migeon, C., Boulard, T., 2012. Measurement and CFD simulation
of microclimate characteristics and transpiration of an Impatiens pot plant crop in a
greenhouse. Biosystems Engineering 112, 22-34.
Kittas, C., 1994. Détermination du coefficient global de transmission de chaleur à travers la
paroi d'une serre. Agricultural and Forest Meteorology 69, 205-221.
Kittas, C., Draoui, B., Boulard, T., 1995. Quantification du taux d'aération d'une serre à
ouvrant continu en toiture. Agriculture and Forest Meteorology 77, 95-111.
Lamrani, M.A., Boulard, T., Roy, J.C., Jaffrin, A., 2001. Airflows and temperature patterns
induced in a confined greenhouse. J. Agric. Eng. Res. 78, 75-88.
Launder, B.E., Spalding, D.B., 1974. The numerical computation of turbulent flows. Comput
Method. Appl. M., 269-289.
Lesieur, M., Compte, P., Lamballais, E., Métais, O., Silvestrini, J., 1997. Large-Eddy
Simulations of shear flows. Journal of Engineering Mathematics 32, 195-215. .
Makhlouf, S., 1988. Expérimentation et modélisation d’une serre solaire à air avec stockage
par chaleur latente assisté par pompe à chaleur de déshumidification. Laboratoire
Ecothermique, CNRS - Valbonne, Université de Nice.
Mercado, L.M., Huntingford, C., Gash, J.H.C., Cox, P.M., Jogireddy, V., 2006. Improving the
representation of radiation interception and photosynthesis for climate model
applications. Tellus 59B, 553-565.
Mesmoudi, K., Soudani, A., Serir, L., 2008. Modèle de bilan énergétique d'une serre sans
couvert végétal. Energies Renouvelables 11, 51-64.
Mistriotis, A., Arcidiacono, C., Picuno, P., Bot, G.P.A., Scarascia-Mugnozza, G., 1997a.
Computational analysis of ventilation in greenhouses at zero-and low-wind-speeds.
Agricultural and Forest Meteorology 88, 121-135.
221
Références bibliographiques
Mistriotis, A., Bot, G.P.A., Picuno, P., Scarascia-Mugnozza, G., 1997b. Analysis of the
efficiency of greenhouse ventilation using computational fluid dynamics. J. Agric.
Eng. Res. 85, 217-228.
Molina-Aiz, F.D., Valera, D.L., Gil, J.A., Pena, A.A., 2006. Optimisation of Almeria-type
greenhouse ventilation performance with computational fluid dynamics. International
sysposium on greenhouse cooling : methods technologies and plant response, Almeria,
Spain.
Monteil, C., Amouroux, M., 1997. Analyse du comportement thermique du sol d'une serre
agricole par simulation dynamique. Archives ouvertes.
Monteil, C., Issanchou, G., Amouroux, M., 1991. Modèle énergétique de la serre agricole.
Journal of physics 3, 429-455.
Monteith, J.L., 1993. Using tube solarimeters to measure radiation intercepted by the crop
canopies ans to analyse stand groth. Delta-T Devices
Montero, J.I., Munoz, P., Anton, A., Iglesias, N., 2005. Computational fluid dynamic
modelling of night-time energy fluxes in unheated greenhouses, In: VanStraten,
G.B.G.P.A.V.W.T.M.M.L.M.F. (Ed.), Proceedings of the International Conference on
Sustainable Greenhouse Systems, Vols 1 and 2, pp. 403-409.
Morille, B., Migeon , C., Bournet, P.E., 2012. Effect of the heterogeneity of the radiation
distribution on the crop activity. Acta Hort. (ISHS) 952, 755-762.
Nebbali, R., 2008. Modélisation de la dynamique du climat interne distribué dans une serre de
culture. Thèse de doctorat de l'Université de Franche-Comté
Norton, T., Sun, D.W., Grant, J., Fallon, R., Dodd, V., 2007. Applications of computational
fluid dynamics (CFD) in the modelling and design of ventilation systems in
agricultural industry : a review. Bioresource technology 98, 2386-2414.
Papadakis, G., Frangoudakis, A., Kyristis, S., 1992. Mixed, forced ans free convection heat
transfer at a greenhouse cover. J. Agric. Eng. Res. 51, 191-205.
Piscia, D., Montero, J.I., Baeza, E., Bailey, B.J., 2012. A CFD greenhouse night-time
condensation model. Biosystems Engineering 111, 141-154.
Prandtl, L., 1925. Z. angew. Math. Mech 5, 136-139.
Reichrath, S., Davies, T.W., 2002. Computational fluid dynamics simulations and ventilation
of the pressure distribution on the roof of a commercial multi-span Venlo-type
glasshouse. J. Wind Eng. Ind. Aerod. 90, 139-149.
222
Références bibliographiques
Roy, J.C., Boulard, T., Kittas, C., Wang, S., 2002. Convective and ventilation transfers in
greenhouses, Part 1: the greenhouse considered as a perfectly stirred tank. Biosystems
Engineering 83, 1-20.
Ruther, M., 1985. Natural ventilation rates of closed greenhouses. Acta Horticulturae 170,
185-191.
Seginer, I., 1984. On the night transpiration of greenhouse roses under glass or plastic cover.
Agriculture and forest Meteorology 30, 203-216.
Siegel, R., Howell, J., 2002. Thermal Radiation Heat Transfer, Taylor & Francis, New York,
NY.
Stanghellini, C., 1987. Transpiration of greenhouse crops. An aid to climate management.
Agricultural University, Wageningen.
Taine, J., Petit, J.P., 2003. Transferts thermiques: introduction aux sciences des transferts.
Dunod.
Tennekes, H., Lumley, J.L., 1972. A First Course in Turbulence. MIT Press, Cambridge, New
York.
Tong, G., Christopher, D.M., Li, B., 2009. Numerical modellinf of temperature variations in a
chinese solar greenhouse. Computers and Electronics in Agriculture 68, 129-139.
Von Elsner, B., Briassoulis, D., Waaijenberg, D., Mistriotis, A., Von Zabeltitz, C., Grataud,
J., Russo, G., Suay-Cortes, R., 2000a. Review of structural and functional
characteristics of greenhouses in European Union Countries, Part II : typical designs.
J. Agric. Eng. Res. 75, 111-126.
Von Elsner, B., Briassoulis, D., Waaijenberg, D., Mistriotis, A., Von Zabeltitz, C., Grataud,
J., Russo, G., Suay-Cortes, R., 2000b. Review of structural and functional
characteristics of greenhouses in European Union Countries: Part I, design
requirements. J. Agric. Eng. Res. 75, 1-16.
Wang, K., 2001. Simulation model for population dynamics of Frankliniella occidentalis
(Thysanoptera : Thripidae) on greenhouse cucumber. Environmental entomology 30,
1073 -1081
Wang, S., Boulard, T., Haxaire, R., 1999a. Air speed profiles in a naturally-ventilated
greenhouse with a tomato crop. Agricultural and forest Meteorology 96, 181-188.
223
Références bibliographiques
Wang, S., Boulard, T., Haxaire, R., 1999b. Experimental and numerical studies of airflow,
temperature and humidity distributions in a greenhouse Tunnel, Proceedings of the
99th International Conference on Agricultural Engineering, III, pp. 34-38.
Wang, S., Deltour, J., 1999. Lee-Side ventilation-induced air movement in a large-scale
multi-span greenhouse. J. Agric. Eng. Res. 74, 103-110.
Wang, X., Zhang, Y., Riskowski, G.L., Ellis, M., 2002. Measurement and analysis of dust
spatial distribution in a mechanically ventilated pig building. Biosystems Engineering
81, 225-236.
Wiebelt, J.A., Henderson, J.B., 1978. Theoretical thermal modeling of leaf with experimental
verification. Agricultural Meteorology 19, 101-111.
Yang, X., 1995. Comments on 'thermal and aerodynamic conditions in greenhouses in relation
to estimation of heat flux and evapotranspiration. Agricultural and Forest Meteorology
77, 131-136.
224
Sitographie
Sitographie
[1] http://fr.wikipedia.org/wiki/Fichier:Transfert_radiation.png
[2] http://www.cnrs.fr/cw/dossiers/dosclim1/motscles/Images/spectreRS.html
[3] http://www.enzyme.wikibis.com/chlorophylle.php
[4] http://www.piaf-archives.org/espace-formation/file.php/10/m8section3_web_web/co/
m8section3_8.html
[5] http://www.sceptiques.qc.ca/activites/conferences/fevrier-2008
[6] http://www-energie2.arch.ucl.ac.be/donn%C3%A9es%20climatiques/1.3.3.3.htm
225
Liste des figures
Figure 4 : Dispositifs de serre existants : cultures sur tablettes (a) et sur géotextile (b),
ouvrants latéraux (c) et de toiture (d), écrans d'ombrage (e et f), chauffage localisé (g)......... 30
Figure 5 : Description des phénomènes physiques mis en jeu dans une serre. ........................ 32
Figure 9 : Longueur d'onde absorbée par le végétal pour réaliser la photosynthèse [3]. ......... 38
227
Liste des figures
Figure 22 : Capteurs de vitesse de l'air boule chaude (a), anémomètre sonique 2D (b),
Figure 25 : Répartition des capteurs de rayonnement PAR dans le couvert végétal. .............. 91
228
Liste des figures
Figure 36 : Evolution temporelle l’humidité absolue pour la nuit claire. .............................. 101
Figure 37 : Evolution temporelle de l’humidité absolue pour la nuit nuageuse. ................... 101
229
Liste des figures
Figure 52 : Variation du flux thermique horizontal dans le sol entre les positions 1 et 2...... 109
Figure 53 : Variation du flux thermique horizontal dans le sol entre les positions 2 et 3...... 109
Figure 54 : Variation du flux thermique horizontal dans le sol entre les positions 3 et 4...... 109
Figure 55 : Variation du flux thermique vertical dans le sol pour les 4 positions. (Temps
Figure 56 : Disposition des capteurs de rayonnement dans le couvert végétal pour une
Figure 57 : Décroissance du rayonnement dans le couvert végétal avec écran d’ombrage. .. 112
Figure 58 : Décroissance du rayonnement dans le couvert végétal sans écran d’ombrage. .. 112
230
Liste des figures
Figure 66 : Différence d’évapotranspiration des plantes entre les deux balances (balance
Figure 67 : Evapotranspiration des plantes sur la balance 1 (sans film). ............................... 117
Figure 69 : Résistances aux transferts hydriques pour une feuille et pour un couvert
Figure 73 : Evolution temporelle de la température d’un sol dont la surface est soumise à
Figure 74 : Evolution temporelle de la température d’un sol dont la surface est soumise à
Figure 75 : Evolution temporelle de la température d’un sol dont la surface est soumise à
profondeurs............................................................................................................................. 175
profondeurs............................................................................................................................. 175
231
Liste des figures
Figure 91 : Variations des conditions aux limites pour les nuits claire (a) et nuageuse (b)
feuille : simulations et mesures. (cf Figure 23 pour la position des capteurs) ....................... 190
Figure 97 : Champ d'humidité relative (%) dans la serre à 5h et zoom sur la toiture. ........... 192
Figure 98 : Champ d'humidité absolue (kg.kg-1) dans la serre à 1h (a) et à 5h (b). ............... 192
232
Liste des figures
Figure 99 : Evolution des températures de toiture. (cf Figure 23 pour la position des
Figure 102 : Evolution du flux de transpiration et de la condensation au cours de la nuit. ... 197
Figure 103 : Evolution de l'humidité absolue moyenne dans la serre au cours de la nuit
Figure 104 : Distribution de la condensation le long de la toiture pour une nuit claire ......... 198
Figure 105 : Distribution de la condensation le long de la toiture pour une nuit nuageuse ... 198
Figure 110 : Champ des vecteurs vitesse (m.s-1) de l'air simulée à 12h................................. 204
Figure 112 : Evolution de la température de l’air simulée et mesurée.au dessus du végétal . 205
Figure 114 : Evolution de l’humidité de l’air simulée et mesurée dans le couvert végétal. .. 206
Figure 115 : Evolution de la température de surface du sol simulée et mesurée. .................. 206
Figure 117 : Champ de température (K) de l'air simulée dans le végétal à 12h..................... 207
Figure 118 : Champ de température (K) de feuille simulée à 12h. ........................................ 208
233
Liste des figures
Figure 119 : Evolution des températures de feuille simulées pour différentes hauteurs
Figure 125 : Evolution temporelle du bilan d’énergie réalisée sur la serre............................ 214
234
Liste des tableaux
Tableau 3: Expressions du nombre de Nusselt moyen pour une plaque plane. ....................... 72
Tableau 9: Température et vitesse de l’air moyennes lors des 3 journées étudiées. .............. 102
Tableau 17 : Comparaison des RMS de la nuit claire et de la nuit nuageuse. ....................... 195
Tableau 18 : Récapitulatifs des flux de chauffage (en W.m-2 de sol). ................................... 198
Tableau 19 : Quantification des pertes radiatives et convectives au niveau de la toiture. ..... 199
235
Liste des tableaux
Tableau 20 : Proportion des pertes radiatives et convectives au niveau de la toiture. ........... 199
236
Annexe A
Annexe A
Figure A1 : Evolution des températures de l’air simulées et mesurées pour une nuit
nuageuse.
237
Annexe A
238
Annexe B
Annexe B
Keywords: CFD, climate, greenhouse, airflow, radiation, heat transfer, transpiration, crop
Abstract
In order to improve crop quality and reduce consumption of water, energy and
fertilizer in greenhouses, it is crucial to have a precise knowledge of the internal
distributed climate. During the day, because of the path of the sun, there are important
variations in the amount of solar radiation reaching the crop, which strongly impacts
climate and therefore crop activity. Up to now, this mechanism was scarcely considered
in CFD studies. In this work, a 2D pseudo-unsteady CFD approach is implemented to
compute the radiation, temperature and humidity variations at plant and greenhouse
levels during a solar day in summer. The studied greenhouse is a Venlo glasshouse
containing an Impatiens crop under oceanic temperate conditions. The Navier-Stokes
equations were solved using the standard k-ε turbulence model and a bi-band radiative
submodel distinguishing short and long wavelength radiations was employed. A user-
defined crop submodel including the aerodynamic drag of the crop considered as a
porous medium as well as the sensible and latent heat exchanges occurring at plant level
were included. Boundary conditions were set each hour, with varying radiation
intensity, diffuse fraction and direction, soil heat flux and inlet temperature. The model
makes it possible to assess the climate evolutions of the air and leaf temperature
distributions inside the greenhouse.
INTRODUCTION
Heterogeneity of solar radiation, temperature and humidity in greenhouses strongly
affect crop activity which impacts production in terms of uniformity and quality (Bartzanas et
al., 2004). Among the latter variables, the effect of solar radiation reaching the crop is crucial
since it governs plant growth through photosynthesis, induces a rise of leaf temperature, and
directly intervenes in sensible and latent heat exchanges between the crop and the ambient air.
In order to better control plant development and provide growers with useful resource
management tools (e.g. water supply control), it is hence essential to get a quantitative
understanding of the climate heterogeneity and particularly of the solar radiation distribution.
To reach this goal, the Computational Fluid Dynamics (CFD) tool is now commonly
employed, yielding fast and reliable results (Norton et al., 2007). The first simulations used to
model flow and temperature distribution in greenhouses (e.g. Mistriotis et al., 1997) paved the
239
Annexe B
way to more realistic ones implementing a canopy submodel to take account of the influence
of the crop on the microclimate (Boulard and Wang, 2002; Majdoubi et al., 2009) and
modelling the radiative mechanisms by solving the Radiative Transfer Equation (RTE) for the
luminance (Bournet et al., 2007; Baxevanou et al., 2008).
However, solar radiation intensity, diffuse fraction and orientation vary during the day.
For a more realistic approach, it is therefore important to study the climate on a diurnal time
scale and include the effect of the path of the sun in the simulations. In this prospect, Boulard
and Wang (2002) used a submodel to calculate the distribution of solar radiation within the
greenhouse but did not take account of complex mechanisms such as light scattering due to
cladding material or multiple reflections. Considering an arc type tunnel plastic greenhouse,
Fidaros et al. (2007) and Baxevanou et al. (2008) used a similar approach to calculate the path
of the sun but dealt with the latter mechanisms by solving the RTE. Nevertheless, they
neglected the crop transpiration in their model.
In the present study, a commercial CFD code is used to investigate the effect of solar
and atmospheric radiations combined with the interactions of the crop with the local climate
during a solar day inside a Venlo glasshouse. Particular attention is paid to the heterogeneity
of radiation, temperature and latent and sensible heat flux distributions at canopy level and to
the influence of solar radiation distribution on the energy balance of the crop as a function of
time.
Numerical Model
Two-dimensional simulations were performed using the commercial code Fluent 6.3.
The CFD approach consists in discretizing space according to the finite volume method and
solving the transport Navier-Stokes equations over the discretized space for the quantities of
interest. The Boussinesq model was activated and the standard k-ε (turbulent kinetic energy
and dissipation rate) closure (Launder and Spalding, 1974) was chosen as turbulence model,
because it offers a good compromise between turbulence description accuracy and calculation
cost (Bartzanas et al., 2007).
The code makes it possible to include source terms in the transport equation to take
account of specific processes. The radiative flux q rad which appears in the energy
conservation equation was determined from a radiation submodel by integrating the
monochromatic radiative flux over the whole wavelength spectrum. The Discrete Ordinates
(DO) model was used to solve the Radiative Transfer Equation (RTE) for a finite number of
discrete solid angles, each associated with a fixed vector direction. This model permits the
modelling of radiation at semi-transparent walls and is suitable for media with an absorption
coefficient that varies across the spectrum. In this study, a bi-band radiation model taking into
account both short wavelength [0.1-3µm[ solar radiation and long wavelength [3-100µm]
atmospheric radiation was activated (Bournet et al., 2007). For each spectral band, the
refractive index n and the absorption coefficient κ λ were defined for semi-transparent media
(glass and polyethylene) whereas the emissivity was set for opaque surfaces (ground).
Crop model
The crop is considered as a porous medium which exerts a drag effect on the airflow.
The resulting pressure loss is represented by a source term in the momentum equations
calculated from the formula provided by Wilson (1985). The canopy gets a certain amount of
solar radiation G (W m-2), and exchanges both a sensible heat flux Qs (W m-2) with the
240
Annexe B
ambient air and a latent heat flux Ql (W m-2) by transpiration. Hence, the energy balance on
the leaves may be written as:
Gabs − Ql − Qs = 0 (1)
with Qs = 2 LAI v ρ a C p ( Tc − Ta ) / ra (2)
and Ql = λ LAI v ρ a ( w f − wa ) /( rs + ra ) (3)
Gabs is the absorbed radiation in each cell of the crop and is approximated by the Beer’s law
along the vertical axis of the canopy (extinction coefficient = 0.44). The incident radiation
reaching the top of the canopy is directly accessed from the resolution of the RTE.
LAI v is the leaf area density defined as the ratio of the leaf area index (2.91 for the
considered Impatiens crop) by the plant height (0.2 m), ρ a is the air density (Kg m-3), C p is
the specific heat of air (J Kg-1 K-1), Ta and Tc are the temperature of the air and leaves
respectively and ra is the aerodynamic resistance (in s m-1). λ is the latent heat of water
vaporisation in J kgw-1, and wf and wa are the absolute humidity of the leaves at saturation and
of the greenhouse air (deduced from the resolution of the transport equation for the water
vapour mass fraction) in kgwater kg-1. The aerodynamic resistance ra depends on the local
computed air speed and climatic conditions (Bartzanas et al., 2004), on the characteristic
length of the leaves and on the air velocity. The stomatal resistance rs is obtained from an
empirical formula as a function of the vapour pressure deficit and global radiation (Baille et
al., 1994).
The temperature of the leaves Tc can then be obtained by combining equations (1) to (3):
ra
Tc = (Gabs − Ql ) + Ta (4)
2 LAI v ρ a C p
A specific user defined routine was developed to link both the sink of momentum and the
sensible and latent heat source terms with the corresponding conservation equations in the
CFD code.
241
Annexe B
the sun onto the plane of the greenhouse (Fig. 2). The following formulae were imposed as
boundary conditions and were used to calculate the direct and diffuse solar radiation intensity
during the day:
− LF
0.9 + 9.4 sin( H S )
Gdirect = 1370 e (8)
Gdiffuse = 54.8 sin( H S )( TA − 0.5 − sin( H S ) ) (9)
where Hs (°) is the height of the sun (which depends on latitude, time and day of the year and
LF the Linke atmospheric factor defined as:
LF = 2.4 + 14.6 β + 0.4( 1 + 2 β ) ln( Pw ) (10)
where β =0.005 and Pw is the water vapor pressure.
The chosen approach was pseudo-unsteady at the diurnal scale with a one hour time step
meaning that after reaching convergence with the desired tolerance, the boundary conditions
corresponding to the following hour were imposed and the calculation was run again from the
previous converged field, which allowed faster convergence.
242
Annexe B
incident solar radiation at the top of the canopy are consistent with the radiation intensity set
as boundary condition (given a 0.8 transmissivity of the glass) and coherent with results
obtained in other studies (Bournet et al., 2009).
Figure 5 shows the simulated transpiration rate along the shelves at different hours in
the upper 2 cm layer of the crop (i.e. 25.5 cm above the shelves). For the transpiration rate,
two sets of three peaks can be observed, in phase with the solar radiation peaks (but of lower
magnitude), highlighting the crucial role played by radiation on latent heat exchanges.
Transpiration is hence increasing when close to midday with a maximum of around 250 W m-
². However, at midday, transpiration is lower at the left ends of the shelves, which could be
explained by the air flow distribution. At 9 am and 3 pm, the same heterogeneity as that one
described previously for the incoming radiation from one side of the shelf to another can be
seen. It is manifestly caused by the position of the sun in the sky and the location of the
openings. Transpiration is however higher at 3 pm than at 9 am because of the higher
temperature, which impacts the values of the stomatal resistance as well as the absolute
humidity of the crop and the air through water vapour pressure variations.
A comparison of the leaf temperature calculated by the crop sub-model with the air
temperature (i.e. the temperature of the porous medium) in the upper 2 cm layer of the culture
at different time is shown in Fig. 6. Both temperatures reach their maximum at around
midday. The difference between air temperature and leaf temperature varies between -0.1 and
2 degrees during the day and is higher at noon, i.e. when the incoming radiation is maximal.
At 6 and 9 am, the leaf temperature almost equals the air temperature. Simulated air
temperature is higher than leaf temperature at 12 am, 3 and 6 pm, stressing the fact that the
plants effectively lose heat by transpiration. This process is enhanced by a high level of solar
radiation reaching the crop. Furthermore, leaf temperature is higher in zones where air
velocity is lower, i.e. at the left of both shelves, because of reduced air cooling. The air
temperature varies to up to 2 degrees inside the canopy along the shelves. These spotted
heterogeneities may impact the activity and growth of the crop, consequently resulting in a
reduction of the homogeneity and quality of the production.
CONCLUSION
Modelling of the distributed climate inside a Venlo glasshouse was performed using a
CFD tool taking into account the combined effect of incident radiation and interactions at
crop level. Variations in radiation, temperature and humidity distribution over the day were
monitored by implementing evolutive boundary conditions based on experimental
measurements corresponding to the path of the sun on July 1st. The Radiative Transfer
Equation was solved with the help of the Discrete Ordinate method. A user-defined routine
was used to add specific source terms in the model in order to take account of the sink of
momentum due to the friction forces on the leaves and of volumetric sensible and latent heat
source fluxes inside the crop.
The simulations offer valuable information about the evolution of radiation
distribution reaching the crop during the day, which is responsible for important
heterogeneities of temperature and transpiration along the shelves. This knowledge may be
used to develop solutions aiming at reducing these heterogeneities, e.g. shading and thermal
screens, moving shelves or entirely diffuse glass.
243
Annexe B
Literature Cited
Baille, M., Baille, A. and Laury, J.C. 1994. A simplified model for predicting
evapotranspiration rate of nine ornamental species. Climate factors and leaf area. Scientia
Horticulturae, 59:217-232.
Bartzanas, T., Boulard, T. and Kittas, C. 2004. Effect of vent arrangement on windward
ventilation of a tunnel greenhouse. Biosystems Engineering 88(4):479-490.
Bartzanas, T., Kittas, C., Sapounas, A., Nikita-Martzopoulou, C. 2007. Analysis of airflow
through experimental rural buildings: Sensitivity to turbulence models. Biosystems
Engineering 97: 229-239.
Baxevanou, C., Bartzanas, T., Fidaros, D. and Kittas, C. 2008. Solar Radiation distribution in
a tunnel greenhouse. Acta Hort. 801:855-862.
Boulard, T. and Wang, S. 2002. Experimental and numerical studies on the heterogeneity of
crop transpiration in a plastic tunnel. Computers and Electronics in Agriculture 34:173-
190.
Bournet, P.E., Ould Khaoua S.A. and Boulard, T. 2007. Numerical prediction of the effect of
vents arrangements on the ventilation and energy transfers in a multi-span glasshouse using
a bi-band radiation model. Biosystems Engineering 98:224-234.
Bournet, P.E., Kichah, A., and Chassériaux, G. 2009. Analysis of radiation and
transpiration heterogeneity inside a greenhouse Impatiens crop using Computational Fluid
Dynamics. Acta Hort., in press.
Fidaros, D., Baxevanou, C., Bartzanas, T., and Kittas, C. 2007. Thermal Behaviour of a
Tunnel Arc Greenhouse during a Solar Day. Acta Hort. 801(2):893-900.
Launder, B.E. and Spalding, D.B. 1974. The numerical computation of turbulent flows.
Computational Method for Applied Mechanical Engineering 3:269-289.
Majdoubi, H., Boulard, T., Fatnassi, H. and Bouirden, L. 2009. Airflow and microclimate
patterns in a one-hectare Canary type greenhouse: an experimental and CFD assisted study.
Agricultural and Forest Meteorology, 149(6-7): 1050-1062.
Mistriotis, A., Bot, G.P.A., Picuno, P., and Scarascia-Mugnozza, G. 1997. Analysis of the
efficiency of greenhouse ventilation using computational fluid dynamics. Journal of
Agricultural Engineering Research 85:217-228.
Norton, T., Sun, D., Grant, J., Fallon, R., and Dodd, V. 2007. Applications of computational
fluid dynamics (CFD) in the modelling and design of ventilation systems in the agricultural
industry: A review. Bioresource Technology 98: 2386–2414.
Richards, P.J. and Hoxey, R.P. 1993. Appropriate boundary conditions for computational
wind engineering models using the k - ε turbulence model. Journal of Wind Engineering
and Industrial Aerodynamics 46-47:145-153.
Wilson, J.D 1985. Numerical studies of the flow through a windbreak. J. Wind Eng. Ind.
Aerodynamics 21:119-154.
244
Annexe B
Tables
Table 1. Values used in the model
Figures
West East
35 m
σ Tsky4 σ Tsky4
20 cm Canopy
Uniform outside air temperature 7.5 cm Container
1000
Radiation intensity (W/m²)
800
600
400
200
0
0 12 24
Fig. 2. Solar global ( ), direct ( ) and diffuse ( ) radiation during the day (July 1st) in the
region of Angers under clear skies.
245
Annexe B
West East
6 am
12 am
6 pm
Fig. 6. Distribution of the computed leaf (filled symbols) and air (hollow symbols)
temperature at 1.005 m height.
246
Annexe C
Annexe C
Abstract
The crop activity inside a greenhouse is strongly affected by the local
microclimatic conditions (radiation, temperature and humidity of air in particular). In
order to reduce the consumption of energy and water, the control of these parameters is
of prime interest. A better knowledge of the local climate may also help to formalize
processes in CFD studies. The objectives of the present study are: (i) to analyse the
evolution of the microclimate in the vicinity of the canopy all day long and (ii) to assess
the spatial heterogeneity of the radiation distribution, air humidity and temperature
above the crop, and crop evapotranspiration. In this prospect, a 100 m2 greenhouse
compartment located in Angers (47° 28’ N, 0° 33’ E) was equipped with a set of sensors
including temperature and humidity probes as well as solar silicium cells and
pyranometers. Crop evapotranspiration was measured at two locations with balances.
Under rather clear sky conditions, the horizontal distribution of the climatic parameters
was analysed as well as their temporal evolutions which were compared to each other.
At a short time scale (2 hours), strong radiation and humidity spatial-heterogeneities are
observed without impact on temperature distribution: different evolutions of their
characteristics all over the day, mainly due to shadows, are reported, but averaging the
values of the climatic parameters over time (from 8 am to 4 pm) erases these
heterogeneities. The investigations of the temporal evolution of the climatic parameters
evidence the role played by the radiation. Knowing both the spatial distribution and
temporal evolution of the crop activity could help to select hypotheses for CFD
simulations as well as to identify the plant requirements and to adapt a strategy for
heating and irrigation.
INTRODUCTION
During the last three decades, numerous studies on the climate inside a greenhouse
have been conducted considering the climatic parameters as homogenous inside the
greenhouse. In these global approaches, the inside climate was described by the values of
physical parameters at one point (Monteil et al., 1991; Roy et al., 2002). Nevertheless, within
the context of environmental crisis accompanied with a high increase of the cost of energy,
247
Annexe C
growers need to optimize the climate management, especially in the vicinity of the crop in
order to reduce the water and energy consumption. This requires a precise knowledge of the
climate distribution next to the crop. Several studies assessing the heterogeneity of the climate
inside a greenhouse have been published in recent years (Bartzanas et al., 2002; Boulard and
Wang, 2002a; Fidaros et al., 2007; Suay et al., 2008). However, most of them focus on the
temperature distribution and only few of them deal with the analysis of the horizontal
heterogeneity of other climatic parameters such as radiation, humidity or transpiration.
Boulard and Wang (2002a) for instance studied the plant transpiration heterogeneity in a
plastic greenhouse and showed that the openings were the main cause of the climate
heterogeneity. A better knowledge of the climate heterogeneity (radiation distribution in
particular) could help growers to estimate the water consumption and to improve the quality
of the crop.
The objectives of this experimental study are to assess and characterize the
heterogeneity of the climatic parameter distributions (downward global radiation, air
temperature and humidity) just above an ornamental potted-plant crop. The temporal
evolutions of these climatic parameters are also investigated all day long.
Experimental device
The experimental greenhouse was located in Angers (47° 28’ North, 0° 33’ Est and 39
m altitude). It consists of a 100 m² (10 m x 10 m) compartment of a Venlo type greenhouse.
The main axis of the greenhouse was orientated North South. Potted Impatiens (Novae-
guinea, cv. ‘Sonic Scarlet’) were used as plant model and uniformly placed on four shelves,
80 cm above the ground, to create a dense canopy with a total area of 18 m². The plant height
was comprised between 28 cm and 32 cm and its density was 10.6 plants per m². The canopy
covered 100% of the shelves and the Leaf Area Index (LAI) was 4.4 m2 of leaves per m2 of
ground covered by the crop.
During the experiments, no shading screen was used and the roof vents were fully
opened to avoid high temperatures. Plants were watered twice during the night (at 8 pm and at
4 am) by flooding the shelves with a complete nutrient solution and flowers were regularly
removed.
248
Annexe C
time ti (s), A (m²) is the shelf area, n is the number of plants on the tray of the balance and m
is the number of plants on the shelves,
- two sonic anemometers give access to the air speed value above the crop (Vin),
- temperature of upper leaves (Tl) was recorded with thermocouples glued to the underside of
leaves.
Several sensors were also placed outside the greenhouse in order to measure air
temperature (Ta-out), absolute humidity (AHout), air speed (Vout) and solar radiation (Rg-out).
Data were recorded every 3 seconds with a data logger system (Campbell Scientific
Micro logger, CR3000 and CR7) and then averaged over ten minute periods. Results were
then analysed using a moving average technique (over 2 hour periods) in order to partly
smooth the noise associated to shadows occasionally covering the sensors. Data from the
balances were directly stored every ten minutes and then averaged over two hour periods.
Experiments were conducted from the beginning of May to the end of June 2010 during a 35
day period. In the present study, three days were specifically analyzed: May 21st (D1), June
3rd (D2) and June 24th (D3) with clear and sunny weather conditions and a global radiation
reaching about 900 W m-2 (Rg-out) at midday solar time (Fig. 2).
Solar radiation
Global radiation not only acts on the photosynthetic activity of the plants, but also
appears to be the main factor controlling the crop transpiration in greenhouses. For each
selected sunny day, the outside solar radiation (Rg-out) follows a bell shape evolution. As
expected, the inside global radiation (Rg-in) strictly follows the evolution of Rg-out with a lower
intensity. The corresponding transmission coefficient is comprised between 31% and 67%,
depending on the moment of the day. The temporal evolutions of this transmission coefficient
(Fig. 3) at three positions above the crop (P2, P3 and P4) are in agreement with the
observations of Soriano et al. (2004) who obtained a similar evolution (within the range 50%-
75% in summer) by using greenhouse scale models in Spain. The lower values obtained in the
present study are mainly due to greenhouse structure shadows.
The differences of the transmission coefficient between two spatial locations are often
of the order of 20% except at the end of the day when the western sensor falls into a
continuous shadows. These results are in agreement with those of Papadakis et al. (1998).
Consequently, the spatial distribution of the Rg-in downward radiation above the crop is not
homogenous at a given time (Fig. 4a). The ratio of radiation measured at two points above the
canopy may exceed 1.5. Different reasons can explain this spatial heterogeneity:
- the greenhouse-structure shadows do not reach each sensor at the same time and during the
same duration,
- the incidence angle of the sun rays is not the same when they reach the vents, the roof or the
other walls (Suay et al., 2008). Moreover the solar rays may sometimes cross twice the glass,
through the roof opening and then through the roof so that the radiation is locally lower.
Averaging the radiation over a 8-hour period (from 8 am to 4 pm) along the day (Fig.
4b) smoothes the spatial heterogeneity of the Rg-in radiation. Excluding point 5 which falls
inside a shadow area after 2 pm so that the average is reduced, the maximum difference
(between point 1 and 3) is not really significant (around 10%). Since D3 corresponds to the
summer solstice, the averaged value of Rg-in is slightly higher.
249
Annexe C
Air humidity
To investigate the time evolution of air humidity and temperature just above the crop,
the values obtained at the central point (P3) were taken as reference (Fig. 5a and b). The
temporal evolution of air humidity above the canopy at P3 varies from one day to another
(Fig.5a). Indeed, as the roof vents were opened, the outside air humidity, the outside air speed
(Table 1) as well as the amount of crop evapotranspiration governs the local air humidity just
above the plants. The evapotranspiration contributes to the increase of air humidity whereas
the ventilation flux may enhance or reduce its value depending on the outside air water
vapour content and speed. The humidity balance inside the greenhouse therefore depends on
the relative weight of these two phenomena.
The temporal evolutions of humidity above the plants during D2 and D3 are similar
since the outside climatic conditions were quite the same for these two days (Table 1); similar
plant evapotranspirations are also reported (see hereafter). The rise of humidity at the
beginning of the day may be attributed to the rise of the evapotranspiration induced by solar
radiation. The lower absolute humidity on sunny day D1 is due to the lower outside air
absolute humidity combined with higher ventilation induced by a higher outside air speed.
The differences of air absolute humidity between each spatial point and the reference
point were plotted every 3 hours over 24 hours (Fig. 6). At night, the spatial distribution of air
absolute humidity is almost homogeneous: humidity difference between locations is lower
than the measurement incertitude of the sensor (5.10-4 gwater g-1dry air). During the day, the
humidity distribution substantially discloses a symmetry shape, the maximum occurring at P3.
This symmetry may be explained by three hypotheses:
- a symmetric distribution of plant evapotranspiration,
- an air flow that conveys the humidity from the edge to the middle of the shelves,
- a higher water vapour concentration in the middle of the shelves whereas the humidity is
better evacuated at the ends of the shelves.
Since global radiation is the main factor influencing crop activity and since its spatial
evolution is not symmetric at a given time all day long as mentioned above, the distribution of
the evapotranspiration is not perfectly symmetric.
The investigation of the air flow direction above the crop remains delicate because the
recorded air speed is often below the threshold provided by the sensor manufacturer.
However, the Richardson number (Ri) may be assessed from a rough estimate of the air speed
and from the leaf-to-air temperature difference (Table 2). Since the value of Ri is about 6, it
can be concluded that the convection mode is locally free or mixed so that an upward flow
may occur near the plant. This upward flow would be supplied by air coming from the edge
and passing through the canopy. These possibilities added to the phenomenon of better
dissipation at the ends of the shelves may explain the observed humidity distribution along the
shelves.
Air temperature
The temporal evolution of the air temperature above the crop is strongly linked to the
evolution of the global radiation (Fig. 6b) with a time offset (maximum reached between 2:30
pm and 3:30 pm). The difference of air temperature for similar sunny days may be explained
by differences of the outside air temperature (Table 3) and by the heat storage phenomenon of
the greenhouse. The higher the outside air temperature (Ta-out), the higher the inside air
temperature (Ta-in) is. However, the air-to-air temperature difference (Ta-in-Ta-out) is about half
the values obtained by Mesmoudi (2010) for the same kind of greenhouse but under semi arid
climate.
250
Annexe C
All over the day, the air temperature difference between the different points (P1 to P5)
remains low. The standard deviation does not exceed 0.5°C which is of the order of
magnitude of the sensor measurement incertitude.
Evapotranspiration
The temporal evolution of the evapotranspiration measured at P4 (Fig. 7a) is
qualitatively similar to that of global radiation with a maximum reached at roughly the same
time. At a fixed location, for sunny days, similar temporal evolution of the evapotranspiration
is reported with small differences attributed to a difference in water vapour content of the air.
Differences can also be explained by the influence of shadows on the radiation reaching the
crop. Yet, integrating evapotranspiration over 8-hour periods (from 8 am to 4 pm) smoothes
these differences which remain lower than 10% (Fig. 7b) and may be considered as
insignificant.
Comparisons between the temporal evolution of the evapotranspiration at positions P2
and P4 for a fixed day show differences (Fig. 8a) that can be mainly attributed to shadows
which impact radiation reaching the crop. However, integrating the evapotranspiration over 8-
hour periods shows that differences between the two balances do not exceed 20% (Fig. 8b).
CONCLUSION
Investigations of the temporal evolution and distribution of the climate parameters in
the vicinity of an Impatiens crop have been conducted. The global radiation received by the
crop discloses spatial heterogeneities (approximately 20% difference from one point to
another) at a given time. At a daily time scale, these heterogeneities are smoothed (lower than
10%). The temporal evolution of the air temperature above the crop is strongly linked to the
outside radiation but presents a time offset. However, the spatial heterogeneities observed on
the radiation distribution have no significant influence on the temperature distribution of the
air above the crop. In contrast, the temporal evolution of the absolute humidity above the crop
results of a balance between the quantity of water vapour provided by the crop and the
characteristics of the outside air (humidity, air speed). Anyway, for each studied day, a spatial
heterogeneity of the absolute humidity distribution was observed with higher air absolute
humidity in the middle of the shelves than at their edge. The evapotranspiration of the crop is
in phase with the evolution of the outside radiation and presents spatial heterogeneity at a
given time not strong enough however to be attributed to anything else than the shadows.
This study established the absence of heterogeneity of the temperature distribution and
a symmetric distribution of the humidity so it is suggested that no other reasons than solar
radiation can explain the evapotranspiration spatial differences.
Acknowledgements
Thanks to D. Lemesle and A. Travers for their valuable technical assistance
throughout the measurement.
Literature cited
Bartzanas, T., Boulard, T. and Kittas, C. 2002. Numerical simulation of the airflow and
temperature distribution in a tunnel greenhouse equipped with insect-proof screen in the
openings. Comput. Electron. Agri. 34:207-221.
Boulard, T and Wang, S. 2002a. Experimental and numerical studies on the heterogeneity of
crop transpiration in a plastic tunnel. Comput. Electron. Agri. 34:173-190.
Fidaros, D.K., Baxevanou, C.A., Bartzanas, T. and Kittas, C. 2008. Thermal behaviour of a
251
Annexe C
Tables
Table 1: Average outside air speed and absolute humidity for the three selected days
Day of measurement D1 D2 D3
Vout (m s-1) 5.45 3.1 2.5
AHout (10-3) (kg kg-1) 7.21 9.03 9.22
Table 2: Characterisation of the convection mode (spatial and temporal average from 8 am to
4 pm)
Day of measurement D1 D2 D3
Vin (m s-1) 0.23 0.22 0.22
Tl-Ta-in (K) 0.62 0.50 0.53
Ri 7.7 5.71 6.26
Table 3: Maximum inside, maximum outside, difference between maximum outside and
inside and average outside air temperature
Day of measurement D1 D2 D3
Ta-in max (°C) 28.2 31.4 33.2
Ta-out max (°C) 24.6 25.7 28.1
Ta-in-Ta-out max (°C) 3.6 5.7 5.1
Ta-out average (°C) 20.7 21.2 24.1
252
Annexe C
1000
Figures D1
900
D2
800 D3
700
Rg-out (W/m²)
600
500
400
300
200
100
0
0:00 4:00 8:00 12:00 16:00 20:00
Time (h)
Fig. 1. Sketch of the experimental setup Fig. 2. Temporal evolution of the
measured outside solar radiation.
1
Transmission coefficient (Rg-in/Rg-out)
P2
0.9
P3
0.8 P4
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
8:00 10:00 12:00 14:00 16:00 18:00
Time (h)
Fig. 3. Spatio-temporal evolution of solar-radiation transmission-coefficient (D1).
700
a) P1
b)
650
600 P2
P3
550
P4
500 P5
Rg-in (W/m²)
450
400
350
300
250
200
150
100
8:00 10:00 12:00 14:00 16:00
Time (h)
Fig. 4. Temporal variation for D1 (a) and averaged spatial values for the three selected
days (b) of the global radiation Rg-in reaching the crop.
a) 0.016 b) 34
D1
0.015 32 D1
D2
D2
Air absolute humidity (g/g)
0.014 D3 30
D3
Air temperature (°C)
28
0.013
26
0.012
24
0.011
22
0.01
20
0.009 18
0.008 16
0.007 14
0:00 4:00 8:00 12:00 16:00 20:00 0:00 4:00 8:00 12:00 16:00 20:00
Time (h) Time (h)
Fig. 5. Temporal evolution of the air absolute humidity (a) and temperature (b) at position
P3 above the crop for the three selected days.
253
Annexe C
x 10
-3 6:00 am x 10
-3 9:00 am x 10
-3 12:00 am
1 1 1
Absolute humidity
Absolute humidity
Absolute humidity
difference (g/g)
difference (g/g)
difference (g/g)
0 0 0
-1 -1 -1
-2 -2 -2
-3 -3 -3
-4 -4 -4
1 2 3 4 5 1 2 3 4 5 1 2 3 4 5
Position of the sensor Position of the sensor Position of the sensor
x 10
-3 3:00 pm x 10
-3 6:00 pm x 10
-3 9:00 pm
1 1 1
Absolute humidity
Absolute humidity
Absolute humidity
difference (g/g)
difference (g/g)
difference (g/g)
0 0 0
-1 -1 -1
-2 -2 -2 D1
-3 -3 -3 D2
D3
-4 -4 -4
1 2 3 4 5 1 2 3 4 5 1 2 3 4 5
Position of the sensor Position of the sensor Position of the sensor
a) 350 b)
D1
300 D2
Evapotranspiration (W/m²)
D3
250
200
150
100
50
0
0:00 4:00 8:00 12:00 16:00 20:00
Time (h)
Fig. 7. Temporal evolution of the evapotranspiration at P4 (a) and average
evapotranspiration compared with the average global radiation (b) for the three
selected days.
a) 300 b)
P2
250 P4
Evapotranspiration (W/m²)
200
150
100
50
0
0:00 4:00 8:00 12:00 16:00 20:00
Time (h)
Fig. 8. Comparisons of the temporal evolution of the evapotranspiration (a) and average
evapotranspiration (b) between P2 and P4 for three selected days.
254
Annexe D
Annexe D
Abstract
Following the rise of energy costs and fuel scarcity, energy becomes today one of
the main greenhouse crop production expenses. In this context, growers have to adapt
their production system and climate management strategy in order to remain
competitive, particularly at night when heat supply is often required and risk of
condensation is enhanced. Numerical tools such as Computational Fluid Dynamics
(CFD) can help predict the climate inside a greenhouse. CFD however is still scarcely
applied in a dynamic way (Nebbali et al., 2011) and under night conditions (Montero et
al., 2005). The aim of this study is to implement an unsteady 2D CFD model to predict
the time-evolution of greenhouse temperature and humidity distributions all night long
in winter, taking account of both radiative transfers and crop interactions with the
inside climate. Beyond the classical conservation equations, a radiative submodel and a
crop-submodel were activated. The latter takes account of both the heat and
transpiration fluxes induced by the plants and their associated mechanical resistances
which depend on the transfer mechanisms from the plant to the air. The boundary
conditions were modified at each time step. Simulations were performed for a clear
night on the basis of data collected in winter 2011 inside a 100 m2 Venlo glasshouse with
Impatiens New Guinea crop grown on shelves. Numerical results were validated against
data recorded inside the greenhouse: roof temperature, inside air temperature,
transpiration flux and air humidity. They highlight the influence of the dynamic
boundary conditions on the evolution of the microclimate inside the greenhouse. This
study demonstrates the ability of the CFD code to simulate the greenhouse climate
evolution with realism. It also suggests a potential exploitation for designing the heating
devices in order to optimize the inside climate and energy consumption for various
outside conditions.
255
Annexe D
INTRODUCTION
During winter at night, it is common to keep the greenhouse vents closed as long as
possible in order to limit heat losses and avoid excessive heating. But following the decrease
of the inside air temperature combined with plant transpiration, the relative humidity may
strongly increase, which may enhance risks of condensation and development of fungal
diseases. The most widespread method to cope with this problem is to evacuate excess
humidity by opening the vents. However, with this method, while humidity is evacuated
through the vents, cool air enters the greenhouse so that heating is required to restore a
temperature adapted to the crop needs. Such a procedure is however not satisfactory due to
the fact it causes an important energy waste.
In the context of natural resource scarcity, optimizing energy consumption becomes
one of the main concerns of growers who take care to minimize their production cost to
remain competitive. They pay a particular attention to better control their greenhouse system
management in order to reduce greenhouse inputs (energy, water, CO2…). This implies a
better knowledge of the physical mechanisms driving the climate inside greenhouses.
To reach this goal, Computational Fluid Dynamics is an efficient tool which makes it
possible to access the entire distribution of the climatic parameters inside the greenhouse. The
device also enables to test different contrasted situations without expensive costs. Up to now
however, only few CFD studies were devoted to night time conditions in greenhouses when
the heating and dehumidification needs are nevertheless the highest. Montero et al. (2005)
studied the night-time greenhouse climate using a 2D steady state model including a radiative
submodel which takes account of the radiative flux exchanged between the greenhouse and
the sky. No heating system was used in their study and the impact of the crop on the
greenhouse climate (and therefore the humidity rate) was not considered. Piscia et al. (2012)
validated a transient 3D model which simulates condensation phenomena that may occur on
walls. Yet, in this study, the crop was only considered as a source of water vapor. No
interaction with the climate was considered so that the crop behavior was independent of the
climatic conditions, which remains questionable.
The objective of the present study is to validate an unsteady 2D model including mass
and radiative transfers, as well as the interaction between crop and inside greenhouse climate
in order to predict the time-evolution of both temperature and humidity distributions all night
long.
256
Annexe D
Campbell Scientific Ltd., Antony, France) monitored the air speed (V) 15 cm above the crop.
The temperature of the leaves of the upper (Tl1), middle (Tl2) and bottom (Tl3) part of the crop
were recorded with copper-constantan (Cu-Cs) thermocouples glued to the underside of plant
leaves and then averaged. A balance was used to estimate the crop evapotranspiration and two
fluxmeters (Captec, France) measured the greenhouse ground surface heat fluxes. Outside the
greenhouse, an aspirated Vaisala HMP45C sensor measured the air temperature (Tout) while
the downward and upward long wave radiations were recorded by a CNR1 pyrradiometer
(Kipp & Zonen, The Nederlands). All the above-mentioned parameters were measured every
3 s and averaged online over 10-min periods with a data logger system (CR3000 and CR7,
Campbell Scientific Ltd., Antony, France). The experiments were carried out from February
to March 2011. For the purpose of the study, a clear night (March 3rd) was analysed.
Numerical Method
The CFD simulations were carried out with the commercially available CFD package,
AnsysFluent 13. This numerical tool solves the 2D convection diffusion equations with the
finite volume technique. Assuming that the fluid is uncompressible, the transport equations
may be written in the general form:
∂φ ∂Uφ ∂Vφ
+ + = Γ∆φ + S φ (1)
∂t ∂x ∂y
where Φ represents the concentration of the non-dimensional transported quantity, namely
momentum, mass (air and water vapour) and energy, U and V are the components of the
velocity vector, Γ is the diffusion coefficient, and SΦ is the source term. The density depends
not only on the temperature, but also on the water vapour content in the air and the buoyancy
force is added as a source term in the momentum equations. The standard k-ε turbulence
model (Launder and Spalding, 1974) was adopted. Radiative transfers were simulated through
the resolution of the radiative transfer equation using the Discrete Ordinate (DO) method
(Bournet et al., 2007). From the mechanical point of view, the crop may be considered as a
porous medium obstructing the inside airflow and inducing a pressure loss estimated from the
Darcy-Forchheimer’s law. The crop activity is modelled by a user defined function which
calculates the energy balance over the leaves Rn=E+H (where Rn is the absorbed net
radiation flux, and H is the sensible heat flux) in order to estimate the latent heat flux (E) due
to plant transpiration, which is a source term of the water vapour conservation equation.
The calculation domain was divided into 36960 cells restricted to the experimental
greenhouse compartment and the outside air just above. In order to reduce the interference of
the outlet condition on the calculated flow, the adjacent downward compartment was also
included in the calculation domain. The initial conditions were obtained by calculating under
steady state conditions the distribution of the climatic parameters performed by the CFD tools
using the boundary conditions at 11 pm. The initial humidity rate was deduced from
measurements at the same time. The boundary conditions (Fig. 2) were directly inferred from
the measurements and modified at each time step (Fig. 3). Along the lateral walls of the
greenhouse, fixed temperatures were set while a heat flux was imposed along the ground
surface in order to take account of the heating. The heat flux intensity corresponded to the
average values of the two ground fluxmeters. An inlet airflow with a fixed temperature was
considered at the upwind vertical boundary of the calculation domain whereas an outflow
condition with zero gradients except for the pressure was chosen at the exit. Lastly, a radiation
sky temperature (Tsky) was introduced at the upper limit of the calculation domain and along
257
Annexe D
the inlet and outlet boundaries. Its value was deduced from the long wave radiation
measurements according to the following expression:
Wsky −radiation = σTsky
4
(2)
where σ is the Stefan Blotzmann constant (σ=5.67.10 W m K ). Calculations were carried
-8 -2 -4
out under unsteady conditions considering a 10 minute time step from 11 pm to 6 am.
CONCLUSION
An unsteady CFD model including convective and radiative transfers as well as the
interaction between the crop and the local microclimate was developed to study the night
evolution of the greenhouse climate. The present study highlights the ability of the numerical
tool to simulate the temperature inside the greenhouse even if the model appears to
258
Annexe D
overestimate mixing inside the shelter. Overestimations of the inside air speeds also enhance
heat exchanges along the roofs, which results in overestimated roof temperatures. The main
weakness of the model lies in its disability to correctly simulate processes driving the
humidity. The predicted humidity rate quickly rises probably because condensation and
leakages for instance are not included although they strongly regulate the water vapour
content of the greenhouse. The integration of these mechanisms would probably improve the
performance of the model and the prediction of the humidity rate evolution inside the
greenhouse in particular.
Literature Cited
Bournet, P.E., Ould Khaoua, S.A. and Boulard, T. 2007. Numerical prediction of the effect of
vent arrangements on the ventilation and energy transfers in a multi-span glasshouse using
a bi-band radiation model. Biosystems Engineering 98:224-234.
Montero, J.I., Munoz, P., Anton, A. and Iglesias, N. 2005. Computational fluid dynamic
modelling of night-time energy fluxes in unheated greenhouses, Acta Hort. 691(1): 403-
410.
Piscia, D., Montero, J.I., Baeza, E. and Bailey, B.J. 2012. A CFD greenhouse night-time
condensation model. Biosystems Engineering 111:141-154.
Launder, B. E. and Spalding, D. B. 1974. The numerical computation of turbulent flows.
Computer Methods in Applied Mechanics and Engineering 3(2): 269-289.
259
Annexe D
Fig 4. Time evolution of measured (dots) and simulated (lines) inside air temperature.
260
Annexe D
Fig. 7. Time evolution of inside crop air and average leaf temperature.
1 am 2 am
Fig. 8. Relative humidity distribution at 1 am and 2 am.
261
Annexe E
Annexe E
Abstract
Humidity control is a key factor in greenhouse management as it impacts the
growth of most greenhouse crops. Excessive humidity and the presence of droplets may
enhance the development of fungal diseases or affect optical properties of the cover
material. For that reason, growers have to adapt their climatic strategy in order to
warranty an optimal environment to the plants. The humidity level inside the
greenhouse results from a water vapour balance between the evaporation of the ground,
the transpiration of the plants, the ventilation and the condensation. Numerical tools
such as CFD are today widely used to predict the climate inside greenhouses. Up until
now however both transient simulations and condensation phenomena have been hardly
addressed in CFD simulations. The aim of this work is to analyse into details the heat
and mass transfers inside the building. In this prospect, 2D transient CFD simulations
considering a monospan glasshouse compartment were conducted, focusing on night-
time cases in winter. The model solves the classical RANS equations in combination with
the k-ε closure for turbulence. Radiative exchanges between the sky and the components
of the greenhouse are included through a specific submodel. A user defined function is
activated to take account of both heat and mass transfers induced by the crop. A water
vapour sink along the walls is also defined to retrieve water from the ambient air when
condensation occurs. Simulations were performed for a clear night and validated against
data recorded inside a 100 m2 Venlo glasshouse with an Impatiens New Guinea crop
grown on shelves. They enhance the ability of the model to predict both air and wall
temperatures of the greenhouse. Condensation appears to strongly impact the resulting
relative air humidity inside the greenhouse and may thus not be omitted from the CFD
approach.
INTRODUCTION
Under low outside temperature conditions, relative humidity inside greenhouses may
strongly increase until the dew point is reached and condensation starts. Condensation mainly
263
Annexe E
occurs on the cover which is the coldest surface inside the greenhouse, particularly at night
under clear sky conditions (Montero et al., 2005). High humidity combined with free water
appearance may lead to a loss of crop quality. As mentioned by Katsoulas et al. (2007) it
creates favorable conditions for the development of fungal diseases and enhances leaf
necrosis, calcium deficiencies and soft and thin leaves. Consequently, the crop production
may be affected due to anatomical changes or delays in the maturity of the plants. The most
commonly adopted strategy to cope with excessive humidity consists in combining heating
and ventilation. But this results in high energy consumption and inside climate heterogeneity
(Campen et al., 2009). Beside ventilation, other techniques have been developed such as
dehumidification by hygroscopic absorption or by condensation (Campen, 2009).
Humidity therefore deserves a particular attention given its impact on the crop. Among
numerical tools available to analyze its characteristics, CFD appears to be powerful. But only
few CFD studies of greenhouse climate include humidity, or condensation processes (Piscia et
al., 2012). According to Pieters and Deltour (1997) however, condensation on the cover is
shown to be the most important sink of water from the inside air at night and may contribute
to 85% to 89% of the removed water vapor. From a dynamic climate model of the
greenhouse, they also showed that fossil heating requirements were underestimated by 15%
for a glasshouse if condensation was neglected. In their study, Piscia et al. (2012) considered
the crop as a constant source of water vapor thus neglecting the interaction of the crop with
the ambient climatic conditions. Moreover, they activated a condensation model determining
the water uptake from the air, but omitting the associated heat flux along the walls.
The present work focuses on the evolution of the air water vapor content of a
glasshouse compartment, taking account of radiative transfers, heating from the ground and
from heat pipes, sensible and latent heat contributions of the crop to the ambient air, and
condensation. Another objective is to assess the impact of condensation on the inside air and
wall temperature, on the inside relative humidity and on leaf temperatures.
264
Annexe E
underside of plant leaves. Outside the greenhouse, an aspirated Vaisala HMP45C sensor
measured the air temperature (Tout) while the downward and upward long wave radiations
were recorded by a CNR1 pyrradiometer (Kipp & Zonen, The Nederlands). All the above-
mentioned parameters were measured every 3 s and averaged online over 10-min periods with
a data logger system (CR3000 and CR7, Campbell Scientific Ltd., Antony, France).
Numerical Method
CFD is based on the resolution of the transient distribution of physical quantities F
such as mass (air and water vapour), momentum and energy which, in 2D form, is governed
by the following convection diffusion transport equation:
∂φ ∂Uφ ∂Vφ
+ + = Γ∆φ + S φ (1)
∂t ∂x ∂y
where U and V are the components of the velocity vector, G is the diffusion coefficient, and
SF is the source term. Turbulence was modelled using the standard k-ε turbulence closure
which is widely used and was validated in several greenhouse CFD studies (Bournet et al.,
2010). Buoyancy effects were included through a linearized dependence of the air density on
temperature and water vapour mixture.
A radiative submodel was activated to take account of long wavelength radiative
transfers occurring both between the components of the greenhouse and between the
atmosphere and the roof. The resolution of the radiative transfer equation (ETR) makes it
possible to calculate the radiative flux which is added to the energy equation as a source term.
The Discrete Ordinate (DO) method (Bournet et al., 2007) was chosen to calculate the
radiation component as it is particularly suitable for semi-transparent materials. The optical
properties of the glass within the far infra-red wavelength range [3-100µm] are an
absoptivity=0.9; a reflectivity=0.1 and a 0 transmissivity (power absorption coefficient= 104
m-1). The refractive index for the glass was set to 1.52.
A specific User Defined Function was activated to take account of the interaction of
the crop with the ambient environment. The plants were assimilated to an equivalent
homogenous porous medium which induces a drag force calculated from the Darcy
Forchheimer low (Boulard and Wang, 2002). Moreover, the crop was divided into elementary
volumes in which the heat and water vapour mass transfers were modelled by sink or source
terms. These terms consist in a sensible heat term which appears in the energy equation, and a
water vapour mass term which affects the water vapour transport equation (Kichah et al.,
2012). They are calculated from an energy balance over the crop.
In the present case, the ground of the greenhouse is made of concrete, meaning that no
evaporation occurred along the floor. Preliminary field surveys were also launched to estimate
the relative weight of the evaporation from the pots compared to the transpiration of the
plants. They showed that evaporation represented less than 2% of the total evapotranspiration
and could therefore be neglected in the simulations.
A specific submodel was also activated to take account of the condensation which
occurs along the inner roof surfaces (no condensation is expected along the vertical walls in
the present case, due to the fact they remain warmer than the saturation temperature). The
water uptake from the ambient air, which happens when the local temperature goes beyond
the dew point, is expressed as a mass flux sink term in the water vapour transfer equation. The
corresponding source term for the energy equation can be obtained by multiplying the mass
flux by the latent heat of vaporization. These source terms are applied to the inside domain
265
Annexe E
cells in contact with the roof. The total rate of mass condensation flux is calculated with Eq. 2
(Bird, et al.,1960; Bell, 2003):
ω h2o ∂ω h 2o A cell wall
m& ′h′′2o = ρD (2)
(ω h2o − 1) ∂n Vcell
: volumetric mass source of water vapor in a computational cell (kg m-3s-1).
: mass fraction of water vapor.
: mixture density (kg m-3).
D: diffusion coefficient (m2s-1).
n: direction normal to the wall and assumed to be the y direction in the formulation.
: face area of a computational cell at wall (m²).
: volume of a computational cell at wall (m3).
The calculation domain was divided into 36900 cells restricted to the experimental
greenhouse compartment plus the adjacent one leeward (to avoid interference with the outlet),
and to the outside air just above. Preliminary simulations were undertaken to get grid
independence results. Calculations were carried out under unsteady conditions considering a
10 minute time step from midnight to 6 am. Initial conditions (velocity, pressure and
temperature) were inferred from a steady state simulation performed with the boundary
conditions recorded at 11 pm. The initial relative humidity was reset from the value recorded
at midnight. The boundary conditions (Fig. 2) were established from measurements and
updated at each time step (Fig. 3). Fixed temperatures were set on the lateral walls while a
heat flux was imposed both along the ground surface and locally inside the greenhouse to
account for overhead heating. A fixed temperature was set at the upwind vertical boundary
domain and an outflow condition at the exit of the calculation domain. The measured infrared
radiation was imposed along the inlet, outlet and upper boundaries of the calculation domain,
through a sky temperature.
RESULTS AND DISCUSSION
The CFD simulations were carried out with the commercially available CFD package,
AnsysFluent 13. For the purpose of the study, a clear night (March 3rd) was analysed.
Simulations were run for 36 ten-min time step consisting each of 2000 iterations. The purpose
was to analyse the time evolution of the climate inside the greenhouse and to evidence the
condensation process occurrence. The prescribed boundary conditions and heat flux source
terms are plotted in Fig. 3. During the night, the air, lateral wall and sky temperature remained
relatively stable with a decrease of less than 3 K over a 6 hour period. The ground heating
started rapidly and was completed by overhead heating once ground heating was insufficient
to maintain the set point temperature. It reached 60 W m-2 at all night.
The time evolution of the inside air at three locations is also shown in Fig. 4. Both
measured and predicted temperatures follow the same trend with a slight temperature drop
over the night (less than 1 K). As reported by measurements, predicted roof temperatures (Fig.
5) remain lower than the greenhouse air temperature, the maximum difference between these
temperatures being 10 K. The model also shows its ability to predict lower temperatures
inside the crop than above or below the canopy. This may be explained by the fact that the
local inside air being warmer than the leaves, the air inside the canopy gives up heat to the
plants (Fig. 6).
The time evolution of the measured and predicted relative humidity is plotted in Fig. 7.
To estimate the model performance, the RMS (Root Mean Square) was calculated for each
variable:
266
Annexe E
where and are respectively the measured and simulated values. The following results
were obtained: RMS= 3.54 K for roof temperature, RMS=1.55 K for greenhouse air
temperature, RMS= 0.11 K for leaf temperatures and RMS= 1.48 %RH for relative humidity
measured above a canopy HR3 (Fig.1). Predicted results are therefore very close to the
measurements except for the roof temperatures where they are higher than measured ones.
(Fig. 4 and 5). These higher simulated temperatures of the roof may be due to measurement
errors of long wave radiations used as boundary condition in the sky. Indeed, the sensitivity
tests which were conducted showed that the temperatures of the roof are much influenced by
the long wave radiations of the sky.
The simulated transpiration fluxes are also very close to measured ones with and order
of magnitude of 10-5 kg.s-1(Fig.8). They disclose a decrease from 0 to 6 am, which may be
explained by the increase in the relative humidity inside the greenhouse (which reduces the
gap between the humidity of air inside greenhouse and leaves; the relative humidity of leaves
being assumed to be equal to 100% RH).
The flux of condensation (Fig.9) increased with time all over the night until a balance with the
transpiration flux was reached.
CONCLUSION
A comprehensive study of the air, roof and leaf temperature and air humidity inside a
glasshouse during the night was undertaken. Previous simulations undertaken for a similar
case, but without taking condensation into account showed that the model was able to predict
air temperature evolution with a relatively good accuracy, but failed in simulating leaf
temperature and air relative humidity correctly (Morille et al., 2012). The air quickly reached
water vapour saturation, which was not the case during the experiments. The present study
combines a radiative submodel, together with a crop and a condensation submodel. The
condensation process occurred mainly along the roof which was cooler than the outside and
ambient air, due to important radiative losses induced by low sky temperatures (clear night).
The roof clearly appears as a key factor regulating the inside relative humidity and thus inside
water content as it participates to water uptake from the ambient air.
The CFD model showed its ability to correctly predict the evolution of the main
characteristics of both the inside air, the greenhouse cover and the canopy: temperature and
humidity in particular. It provides promising information on the greenhouse behaviour and
could be applied to explore the influence of several climatic or greenhouse parameters such as
the sky temperature, the plant characteristics, or the optical properties of the cover.
Literature Cited
Bell, B. (2003). Application Brief: Film Condensation of Water Vapor. Lebanon, New
Hampshire: Fluent, Inc.
Bird, R.B., W.E. Stewart, and E.N. Lightfoot. (1960). Transport Phenomena, New York: John
Wiley & Sons.
Boulard, T. and Wang, S. (2002). Experimental and numerical studies on the heterogeneity of
crop transpiration in a plastic tunnel. Computers and Electronics in Agriculture, 34, 173-
190.
Bournet, P.-E., Boulard, T., 2010. Effect of ventilator configuration on the distributed climate
of greenhouses: a review of experimental and CFD studies. Comput. Electron. Agric.
74, 195–217.
267
Annexe E
Bournet, P.E., Ould Khaoua, S.A. and Boulard, T. 2007. Numerical prediction of the effect of
vent arrangements on the ventilation and energy transfers in a multi-span glasshouse
using a bi-band radiation model. Biosystems Engineering 98:224-234.
Campen, J., Kempkes, F. and Bot, G. 2009. Mechanically controlled moisture removal from
greenhouses. Biosystems Engineering, 102(4):424-432.
Campen, J. 2009. Dehumidification of Greenhouses, PhD Thesis, Wageningen University,
Wageningen, NL, 117 pages.
Katsoulas, N., Manoloaraki, C., Kittas, C. and Bartzanas, T. 2007. Effects of Anti-Drip Cover
Materials on Microclimate and Production of a Hydroponic Cucumber Crop. Acta Hort.
801 (2):267_274.
Kichah A., Bournet P.E., Migeon C. 2012. Measurement and CFD simulations of the
transpiration and leaf temperatures of an Impatiens pot plant crop in a greenhouse.
Biosystems Engineering 112, 22-34.
Launder, B. E. and Spalding, D. B. 1974. The numerical computation of turbulent flows.
Computer Methods in Applied Mechanics and Engineering 3(2): 269-289.
Montero, J.I., Munoz, P., Anton, A. and Iglesias, N. 2005. Computational fluid dynamic
modelling of night-time energy fluxes in unheated greenhouses. Acta Hort. 691(1): 403-
410.
Morille, B., Genez, R., Migeon, C., Bournet, P.E. and Bouhoun Ali, H. 2012. CFD
Simulations of the Distributed Climate Time-Evolution inside a Glasshouse at Night, ,
International conference of agricultural engineering CIGR-Ageng, July 8-12, Valencia
Spain, 6p.
Piscia, D., Montero, J.I., Baeza, E. and Bailey, B.J. 2012. A CFD greenhouse night-time
condensation model. Biosystems Engineering 111:141-154.
List of Figures
268
Annexe E
Fig. 4. Time evolution of the measured and simulated inside air temperatures.
269
Annexe E
Fig.6. Time evolution of the measured and simulated inside crop air and average leaf
temperatures.
Fig. 7. Time evolution of the measured and simulated inside air relative humidities
above the crop.
Fig. 8. Time evolution of the simulated Fig. 9. Time evolution of the simulated
transpiration flux. condensation flux.
270
Annexe F
#include "udf.h"
DEFINE_GRAY_BAND_ABS_COEFF(user_gray_band_abs,c,t,nb)
{
real abs_coeff = 0;
switch (nb)
{
case 0 :
if (THREAD_VAR(t).fluid.porous)
abs_coeff = 9;
else
abs_coeff = 0;
break;
case 1 : abs_coeff = 0; break;
}
return abs_coeff;
}
#include "udf.h"
#include "sg.h"
#include "mem.h"
/******************************************************************************/
/* **********************Initialisation du domaine *********************************/
/******************************************************************************/
/******************************************************************************/
/************************* Caractérisation de la végétation **************************/
/******************************************************************************/
/******************************************************************************/
/******************Definition des coefficients bifacteur pour calculer rs *****************/
/******************************************************************************/
# define A 421
# define B 41
# define C 0.02
# define Dmax 2
# define rsmini 150 /* Résistance stomatique minimale */
# define lambda 0.0262 /* conductivité thermique de l'air */
# define Visc 1.8e-5 /* viscosité dynamqiue de l'air */
271
Annexe F
/******************************************************************************/
/********** Source de température végétation (Flux chaleur sensible) ************/
/******************************************************************************/
/* on récupère la valeur de la température de la feuille stockée dans user defined memory 0*/
Tinterne=C_UDMI(cell,thread,0); /* C_UDMI(cell,thread,0) variable Tinterne*/
/* on récupère la valeur de l'humidité absolue de la feuille stockée dans user defined memory 1*/
Humfeuille=C_UDMI(cell,thread,1); /* C_UDMI(cell,thread,1) variable
Humidite*/
/******************************************************************************/
/*************** calcul de la résistance aérodynamique *****************************/
/******************************************************************************/
Vitesse_x=C_U(cell,thread);
Vitesse_Y=C_V(cell,thread);
norme=sqrt(pow(Vitesse_x,2)+pow(Vitesse_Y,2));
delta_t=fabs(Tinterne-Tambiant);
Gr=9.81*1/293*delta_t*pow(L,3)*pow(1.22,2)/pow(Visc,2);
Re=1.22*norme*L/Visc;
if (Re*Re<0.1*Gr)
{Nusselt=0.40*pow(Gr,0.25);}
if ((Re*Re>0.1*Gr)&&(Re*Re<10*Gr))
{Nusselt=0.40*pow(Gr+6.92*Re*Re,0.25);}
if (Re*Re>10*Gr)
{Nusselt=0.60*pow(Re,0.5);}
RA=1.22*Chp/(lambda*Nusselt/L);
C_UDMI(cell,thread,2)=RA;
/******************************************************************************/
/********************* Calcul de la résistance stomatique ***************************/
/******************************************************************************/
272
Annexe F
G=C_STORAGE_R_XV(cell,thread,SV_DO_IRRAD,0);
C_UDMI(cell,thread,3)=G;
D=6.11*exp((17.25*(Tambiant-zero))/(237.8+(Tambiant-zero)));
DPV= (D-((1013*Humair)/(0.622+Humair))); /* Déficit de saturation de l'air en mbar
*/
C_UDMI(cell,thread,4)= DPV;
f1=(A+G)/(B+G);
if ((0.1*DPV)>Dmax)
f2=1+(C*pow((0.1*DPV-Dmax),2));
else
f2=1;
RS=rsmini*f1*f2;
C_UDMI(cell,thread,5)=RS;
/******************************************************************************/
/*********************** Calcul de la resistance totale *****************************/
/******************************************************************************/
RT = RA+RS;
/******************************************************************************/
/********************** calcul du coefficient phi1 ********************************/
/******************************************************************************/
phi1=C_R(cell,thread)*Chp*2*LAIV/RA;
/******************************************************************************/
/********************* Calcul de la temperature interne ****************************/
/******************************************************************************/
RayAbs=0;
/*Rayonnement absorbé dans la maille*/
for (nb=0;nb<sg_bee_gees;nb++)
{
real abs_coeff;
real refractive_index;
/* coefficient d'absorption en multibande */
if (nb==0)
{
abs_coeff = 9.0;
}
if (nb==1)
{
abs_coeff = 0.0;
}
RayAbs+=abs_coeff*C_DO_IRRAD(cell,thread,0);
}
C_UDMI(cell,thread,6)=RayAbs;
273
Annexe F
newtemp=Tambiant+RA*RayAbs/(2*C_R(cell,thread)*Chp*LAIV)-(RA*Lv/(2*Chp)*(Humfeuille-
(C_YI(cell,thread,0))))/RT;
/* relaxation*/
Tinterne=(1-rx)*Tinterne+rx* newtemp;
C_UDMI(cell,thread,0)= Tinterne;
phi2=Tinterne*phi1;
/* source term */
source=phi2-phi1*C_T(cell,thread);
C_UDMI(cell,thread,7)=source;
/******************************************************************************/
/****************** Source de vapeur d'eau végetation ********************/
/******************************************************************************/
/* on récupère la valeur de la température de la feuille stockée dans user defined memory 0*/
Tinterne=C_UDMI(cell,thread,0); /* C_UDMI(cell,thread,0) variable Tinterne*/
/* on récupère la valeur de l'humidité absolue de la feuille stockée dans user defined memory 1*/
Humfeuille=C_UDMI(cell,thread,1); /* C_UDMI(cell,thread,1) variable
Humidite*/
/******************************************************************************/
/*************** calcul de la résistance aérodynamique *****************************/
/******************************************************************************/
Vitesse_x=C_U(cell,thread);
Vitesse_Y=C_V(cell,thread);
norme=sqrt(pow(Vitesse_x,2)+pow(Vitesse_Y,2));
delta_t=fabs(Tinterne-Tambiant);
Gr=9.81*1/293*delta_t*pow(L,3)*pow(1.22,2)/pow(Visc,2);
Re=1.22*norme*L/Visc;
Nusselt=0.40*pow(Gr+6.92*Re*Re,0.25);
RA=1.22*Chp/(lambda*Nusselt/L);
274
Annexe F
C_UDMI(cell,thread,2)=RA;
/******************************************************************************/
/***************** Calcul de la résistance stomatique *******************************/
/******************************************************************************/
G=C_STORAGE_R_XV(cell,thread,SV_DO_IRRAD,0);
C_UDMI(cell,thread,3)=G;
D=6.11*exp((17.25*(Tambiant-zero))/(237.8+(Tambiant-zero)));
DPV= (D-((1013*Humair)/(0.622+Humair))); /* Déficit de saturation de l'air en mbar
*/
C_UDMI(cell,thread,4)= DPV;
f1=(A+G)/(B+G);
if ((0.1*DPV)>Dmax)
f2=1+(C*pow((0.1*DPV-Dmax),2));
else
f2=1;
RS=rsmini*f1*f2;
C_UDMI(cell,thread,5)=RS;
/******************************************************************************/
/********************* Calcul de la resistance totale *******************************/
/******************************************************************************/
RT = RA+RS;
/******************************************************************************/
/*********************** calcul du coefficient phi1 ********************************/
/******************************************************************************/
phi1=C_R(cell,thread)*LAIV/RT;
/******************************************************************************/
/************************* Calcul de l'humidité **********************************/
/******************************************************************************/
Pvsat=6.11*exp((17.25*(Tinterne-zero))/(237.8+(Tinterne-zero)));
Newhum=0.622*Pvsat/(1013-Pvsat);
/* relaxation*/
Humfeuille=(1-rx)*Humfeuille+rx*Newhum;
C_UDMI(cell,thread,1)= Humfeuille;
phi2=Humfeuille*phi1;
/* terme source */
275
Annexe F
source=phi2-phi1*C_YI(cell,thread,0);/* Kg / m3 s */
C_UDMI(cell,thread,11)=source;
/******************************************************************************/
/******************************************************************************/
/******************** Initialisation de domaine **********************/
/******************************************************************************/
DEFINE_INIT(INIT_Solution, domain)
{
cell_t c;
Thread *thread;
thread_loop_c (thread, domain)
{
begin_c_loop (c, thread)
{
C_UDMI(c,thread,0)= 288.9;
C_UDMI(c,thread,1)= 0.011;
}
end_c_loop (c, thread)
}
}
Compensation du rayonnement
/******************************************************************************/
/***************** Source compensation de rayonnement absorbé ****************/
/******************************************************************************/
RayAbs1=0;
for (nb1=0;nb1<sg_bee_gees;nb1++)
{
/* coefficient d'absorption en multibande */
if (nb1==0)
{
abs_coeff1 = 9.0;
}
if (nb1==1)
{
abs_coeff1 = 0.0;
}
RayAbs1+=abs_coeff1*C_DO_IRRAD(cell,thread,0);
276
Annexe F
}
C_UDMI(cell,thread,13)=RayAbs1;
source=-RayAbs1;
return source;
}
#include "udf.h"
*/
/* CONSTANTS */
/* */
#define ID 3 /* ID is the zone label for the condensing */
/* surface in the Fluent case file */
#define MW_H2O 18.01534 /* molecular weight of h2o */
#define MW_AIR 28.966 /* Molecular weight of air */
277
Annexe F
#define F_YI_FLUX(f,t,i)F_STORAGE_R(f,t,SV_SPECIES_FLUX_0+i)
/* FUNCTION DECLARATIONS */
double psat_h2o(double); /* computes sat. pressure of h2o */
double steam_hfg(double); /* computes latent heat of h2o */
double t_c(double t_sat); /* used to compute expression used in steam_hfg */
double molefrac(double,double,double);/* computes mole fraction of species */
/* in a binary mixture given mass */
/* fraction and molecular weights */
double massfrac(double,double,double);/* computes mass fraction of species */
/* in a binary mixture given mole */
/* fraction and molecular weights */
/* GLOBAL VARIABLES */
static double sr_ur_alpha = 0.01; /* allows under-relaxation factor in */
/* surface reaction calculation to be */
/* changed in the middle of a run */
enum{
MVOL,
HFG,
QFLUX,
MFSAT,
MFLUX,
};
278
Annexe F
DEFINE_INIT(init_function,domain)
{
/* The initialization function initalizes the under-relaxation rate */
/* for the condensation mass flux */
int n,k;
cell_t c;
Thread *t,*tw,*tf;
/* initialize sr_ur_alpha to SURFACE_RATE_ALHPA */
sr_ur_alpha = RP_Get_Real("wallcond/underrelax");
wall_mf_ur = RP_Get_Real("wallcond/wall_mf_ur");
DEFINE_ADJUST(update,domain)
{
sr_ur_alpha = RP_Get_Real("wallcond/underrelax");
wall_mf_ur = RP_Get_Real("wallcond/wall_mf_ur");
}
DEFINE_ON_DEMAND(post_process)
{
/* This macro takes value of heat flux at wall and stores it in */
/* user-defined memory */
279
Annexe F
280
Annexe F
c0 = F_C0(f,thread);
t0 = THREAD_T0(thread);
new_mf = C_YI(c0,t0,0);
old_mf = mfface;
F_PROFILE(f,thread,position) = old_mf+wall_mf_ur*(new_mf - old_mf);
}
}
end_f_loop(f,t)
}
DEFINE_SOURCE(h2o_condensation,c,t,dS,eqn)
{
/*********************************************************************/
/* */
/* 9/21/01 Reworking old surface reaction macro */
/* */
/* FACE VARIABLES (from solver): */
/* */
/* tface temperature (K) */
/* s1 h2o mass fraction */
/* xf[] coordinates of face centroid (m) */
/* */
281
Annexe F
/* */
/* CELL VARIABLES(from solver): */
/* */
/* tcell temperature (K) */
/* mfcell h2o mass fraction */
/* pcell static pressure (absolute) (Pa) */
/* dcell density (kg/m3) */
/* mdiff mass diffusivity (m2/s) */
/* xc[] coordinates of cell centroid (m) */
/* */
/* OTHER VARIABLES */
/* yf1 h2o mole fraction at wall */
/* yfsat mole fraction of saturated h2o at wall */
/* mffsat mass fraction of saturated h2o at wall */
/* ph2o_f h2o partial pressure at wall (Pa) */
/* yc1 h2o mole fraction */
/* ph2o_c h2o partial pressure in cell */
/* tsat_f saturation temperature of h2o at wall */
/* tsat_c saturation temperature of h2o at cell */
/* grad_mf magnitude of gradient of h2o mass fraction normal to wall */
/* mw1 molecular weight of h2o */
/* mw2 molecular weight of air */
/* sat =1 if tcell <= tsat; =0 otherwise */
/* mflux condensation mass flux (kg/m2-s) */
/* dn dist. from cell centroid to face centroid in normal dir. */
/* At face area */
/* s0 vector from face centroid to cell centroid */
/* e_xi face unit normal vector */
/* A vector in face normal dir. with length equal to face area */
/* */
/* FUNCTIONS */
/* double psat_h2o(tsat); computes sat. pressure of h2o */
/* double molefrac(massfrac1,m1,m2) computes mole fraction */
int sat,n;
double tface,pface,yface;
double tcell,mfcell,pcell,dcell,mdiff;
double yf1,yc1,ph2o_f,ph2o_c,tsat_f,tsat_c,grad_mf,dn,mflux;
double yfsat,mffsat;
double At,sum_At;
double s1,mw1,mw2;
double xc[ND_ND],xf[ND_ND]; /* cell,face centroids */
double s0[ND_ND],A[ND_ND],p0[ND_ND];/* other vectors */
double delta_p;
double old_rate,new_rate;
double source;
real rat1;
real op_pres = RP_Get_Real("operating-pressure");
Domain *d = Get_Domain(1);
Thread *tf;
face_t f;
source = 0;
mflux = 0.0; /* initialize mflux out of loop because it is */
/* possible that a cell will have more than
one */
/* face on condensing surface
*/
282
Annexe F
sum_At = 0.0;
new_rate = 0.0;
old_rate = 0.0;
mw1 = MW_H2O;
mw2 = MW_AIR;
rat1 = 0.0;
c_face_loop(c,t,n)
{
tf = C_FACE_THREAD(c,t,n);
if(THREAD_TYPE(tf) == THREAD_F_WALL && THREAD_ID(tf) == ID )
{
f = C_FACE(c,t,n);
old_rate = C_UDMI(c,t,MVOL);
/* Initialize sat,mflux */
sat = 0;
mflux = 0.0;
/* Face values */
s1 = F_YI(f,tf,0); /* h2o mass fraction */
pface = F_P(f,tf)+op_pres; /* F_P is a gauge pressure */
tface = F_T(f,tf); /* face temperature (K) */
yface = molefrac(s1,MW_H2O,MW_AIR); /* yface is h2o mole fraction */
ph2o_f = yface*pface; /* h2o partial pressure (Pa) */
/* at face */
delta_p = ph2o_f - psat_h2o(tface); /* Difference between face and */
/*
saturation pressure */
/*****************************************************************************/
/******** COMPUTE CONDENSATION RATE (kg/m2-s) IF T <= TSAT ******************/
/*****************************************************************************/
if(sat == 1)
{
/* The next line is should be valid for laminar and turbulent flow. */
/* F_FLUX has units of kg/s */
mflux = F_YI(f,tf,0)*F_YI_FLUX(f,tf,0)/MAX(SMALL, 1.0-F_YI(f,tf,0));
283
Annexe F
*/
temp = F_T(f,ft);
h2o_enth = Enthalpy(sp,temp,0.,op_pres,NULL,NULL);
source = h2o_enth*C_UDMI(c,ct,MVOL)/MAX(SMALL,F_YI(f,ft,0));
}
}
dS[eqn] = 0.0;
return source;
}
/*******************************************************************/
/* H2O PROPERTY ROUTINES */
/*******************************************************************/
284
Annexe F
var1 = PSAT_A*(tsat-PSAT_TP);
ans1 = sum1*(H2O_TC/tsat-1.0);
/* compute exponential to determine result */
/* psat has units of Pascals */
psat = H2O_PC*exp(ans1);
return psat;
}
/* compute tc */
tc = t_c(t_sat);
/* sum e's */
n = 1; /* initialize n */
sum_e = 0.; /* initialize sum_e */
while(n<8)
{
sum_e += e[n]*pow(tc,(double)n);
++n;
}
/* compute y */
y = L_A + L_B*pow(tc,1./3.) + L_C*pow(tc,5./6.) + L_D*pow(tc,7./8.) + sum_e;
return y * L_P;
}
285
Annexe F
denom = m2+m1*((1./mass_frac)-1.);
x = num/denom;
}
return x;
}
286