Ventillation Naturelle 2

République Algérienne Démocratique et Populaire
Ministère de l’Enseignement Supérieure et de la Recherche scientifique

Université Mouloud MAMMERI de Tizi-Ouzou Faculté de Génie de la
construction Département de Génie mécanique
MEMOIRE DE FIN D’ETUDE

En vue de l’obtention du diplôme de Master en Génie mécanique
Option : Froid Chauffage Et Climatisation
Thème :
Etude et simulation numérique de la ventilation naturelle dans
une maison individuelle.
Encadré par : M. FOUCHANE SIDALI Réalisé par :

M. REDAOUI SAMIR
M. TAMACHE MOULOUD
Devant le jury composé de :
Président : M. FERROUK Mohammed
Examinateur : M. BELTAGY Hani
Promotion 2019/2020
DEDICACE
Qu’ALLAH nous pardonne.
Qu’ALLAH nous guide dans le bon chemin.
Tous les mots ne sauraient exprimer la gratitude, l’amour, le respect, la reconnaissance,
c’est tout simplement que : Je dédie cette thèse de master à :
A ma très chère mère : Tu représentes pour moi la source de Tendresse et
l’exemple de dévouement qui n’a pas cessé de m’encourager. Tu as fait plus qu’une mère
puisse faire pour que ses enfants suivent le bon chemin dans leur vie et leurs études.
A Mon très cher Père : Aucune dédicace ne saurait exprimer l’amour, l’estime, le
dévouement et le respect que j’ai toujours pour toi papa. Rien au monde ne vaut les
efforts fournis jour et nuit pour mon éducation et mon bien être. Ce travail est le fruit de
tes sacrifices que tu as consentis pour mon éducation et ma formation le long de ces
années.
A ma très chère fiancée :
A mes très chères sœurs :
Vous avez toujours été présentes pour moi. Je vous souhaite un avenir plein de
joie, de bonheur, de réussite et de sérénité. Je vous exprime à travers ce travail mes
sentiments de fraternité et d’amour. Que dieu vous protège et consolide les liens sacrés
qui nous unissent.
A mes très chers amis :
A mon binôme qui est avant tout un ami TAMACHE MOULOUD avec qui j’ai
passé une année formidable où l’on a partagé des moments de joie et de malheur.
Je ne peux trouver les mots justes pour vous exprimer mon affection. En
témoignage de l’amitié qui nous uni et des souvenirs de tous les moments que nous avons
passés ensemble.
A tous mes camarades et mes amis qui m’ont beaucoup aidé, sans oublié les gens
de la promotion de structures.
Mes enseignants de génie Mécanique :
Qui m’ont éclairé sur ce chemin du savoir et qui doivent voir dans ce travail la fierté
d'un savoir bien acquis.
A tous ceux qui m’aiment et que j’aime.
REDAOUI SAMIR
Dédicace
J’ai toujours pensé faire ou offrir quelque chose à mes parents en signe
de reconnaissance pour tout ce qu’ils ont consenti des efforts rien que pour
me voir réussir, Et voilà, l’occasion est venue.
A ceux qui m’ont donné la vie, symbole de beauté, et de fierté, de
sagesse et de patience.
A ceux qui sont la source de mon inspiration et de mon courage, à
qui je dois de l’amour et la reconnaissance.
A mes parents.
A tous mes frères et sœurs.
A mon binôme REDAOUI Samir et sa famille.A tous mes Amis.
A tous les étudiants de la faculté génie de la construction ainsi que les
étudiants de l’université de Mouloud Mammeri.
REMERCIEMENTS
Nous remercions avant tout Allah, de nous avoir donner la force
d’accomplir ce travail.
Nous tenons à remercier tous les membres de nos familles pour leur soutien
moral et leurs encouragements durant toutes nos années d'études.
Nous tenons à remercier notre encadreur M. FOUCHANE SIDALI
d’avoir accepté d’encadrer notre travail et pour son assistance et sa contribution,
Nous tenons à remercier les membres de jury qui ont bien voulu accepter
d’examiner ce travail.
Nous tenons à remercier tous nos enseignants et ceux qui ont participé à
notre formation.
Nous tenons à remercier enfin nos collègues pour leur aide et leur
encouragement.
NOMENCLATURE
Symbole Signification Unité
K Coefficient de transmission surfacique W/m².

°C
Kmoy. Coefficient de transmission surfacique moyen W/m²
Kg Coefficient de transmission surfacique global W/m².

°C
Kw Coefficient de transmission thermique de la fenêtre, de la porte ou W/m².
de la porte-fenêtre °C
Kg Coefficient surfacique en partie centrale du vitrage W/m².
°C
Kf Coefficient surfacique moyen de la menuiserie W/m².
°C
Kp Coefficient surfacique en partie centrale du panneau opaque W/m².
°C
Kmr,tot Coefficient surfacique moyen de la façade rideau W/m².
°C
Kmri Coefficient surfacique moyen du module i W/m².
°C
Kwn Coefficient de transmission surfacique du vitrage nu W/m².
°C
kl Coefficient de transmission linéique W/m.°
C
S Surface intérieure de la paroi m²
he Coefficient extérieur d’échanges superficiels W/m².°C
tbi Température intérieure de base °C
tbe Température extérieure de base °C
R Résistance thermique m².°C/W
c Capacité thermique massique J/kg.°C
z Différence de niveau entre le niveau de la face supérieure du m

plancher et le niveaudu sol
Tau Coefficient de réduction de température −
Vh Volume habitable du local ou hors œuvre m3
Qs Débit supplémentaire dû au vent m3/h
Qvréf Débit de ventilation extrait de référence m3/h
Qvmin Débit de ventilation extrait minimal de référence m3/h
Qvmax Débit de ventilation extrait maximal de référence m3/h
P pression atmosphérique Pa
Po Perméabilité à l’air d’un ouvrant sous une différence de pression de m3/h.
1 Pa m²
PV Perméabilité à l’air d’un volume sous une différence de pression de m3/h
1 Pa
PL Perméabilité à l’air d’un logement sous une différence de pression m3/h
de 1 Pa
ev Coefficient d’exposition au vent −
Notion des apports
A Apports calorifiques W
AE Apports calorifiques effectifs W
AI Apports calorifiques internes W
APO Apports calorifiques par une paroi opaque W
AREN Apports calorifiques dus au renouvellement d’air W
AINF Apports calorifiques dus aux infiltrations d’air extérieur W
AT Apports calorifiques totaux W
AV Apports calorifiques à travers les parois vitrées W
AVE Apports par rayonnement solaire à travers les parois vitrées W
AVT Apports par transmission à travers les parois vitrées W
Aabs Apports absorbé W
Atra Apports traversant W
BF Facteur de by-pass -
C Coefficient de correction (selon les cas) -
CS Coefficient de simultanéité -
E Ecart diurne °C
DH Degré hygrométrique -
ENT Enthalpie kJ/kgas
EAT Ecart annuel de température °C
FS Facteur solaire -
HR Humidité relative %
HS Humidité spécifique d’air gvap./kga
s
I Rayonnement W/m²
K Coefficient surfacique de transmission calorifique W/m².°C
L Latitude degré
M Masse rapportée à la surface du plancher kg/m²
N Coefficient d’amortissement -
PC Pouvoir calorifique J/m3
PP Pression partielle de vapeur d’eau Pa
TH Température humide de l’air °C
TO Temps officiel h
TR Température de rosée °C
TS Température sèche de l’air °C
TSV Temps Solaire Vrai h
V Volume m3
VS Volume spécifique m3/kgas
Indices

a Coefficient réglementaire ................................... W/m².°C
alt Altitude.................................................................... m
c Coefficient réglementaire ........................................ W/m².°C
e épaisseur ................................................................... m
e’ Coefficient réglementaire................... W/m².°C
h Coefficient d’échange superficiel ............................ W/m².°C
k Coefficient linéique d’émission calorifique………. W/m.°C
m Masse ...................................................................... kg
qm Débit massique ....................................................... kg/s
qv Débit volumique...................................................... m3/s
qvo Débit volumique d’infiltration pour un ouvrant…... m3/h.m²
r Résistance thermique d’échanges superficiels....... m².°C/W

z Différence de niveau (paroi en contact avec le sol) m
..........
α Facteur d’absorption -
λ Conductivité thermique W/m.°C
symboles Signification Unité

ΔTS Différence de températures sèches…………………………… °C
∆te Différence de température équivalente………………………. °C
∆tes Différence de température équivalente non corrigée (conditions
°C
standard) pour une paroi ensoleillée …………………………….
∆tem Différence de température équivalente non corrigée (conditions
°C
standard) pour une paroi l’ombre
Symbole Signification
DTR Document thermique règlementaire
PCM Matériau a changement de phase
AI Apports internes
PH Paroi horizontale
PV Paroi verticale
PVI Paroi vitrée
an air neuf
inf d’infiltration
l latent
lnc local non climatisé
max maximal
min minimal
réf. s référence sensible
sup supérieur
Lettres Gracques
β: Coefficient de la dilatation
ε: taux de dissipation de l’énergie cinétique
μ: viscosité dynamique [Kg/m. s]
Г: Coefficient de diffusion
Φ: Grandeur quelconque
Re : Nombre de Reynolds.
Pr : Nombre de Prandtl.
Nu : Nombre de Nusselt.
Ra : Nombre de Rayleigh.
Gr : Nombre de Grashoff.
ABREVIATION
EDP : équation à la dérivée partielle.
MDF : Méthode des différences finies.
MEF : Méthodes des éléments finies.
MVF : Méthodes des volumes finies.
CFD: Computational Fluid Dynamics.
Gambit: Geometry and Mesh building intelligent Toolkit.
Quick: Quadratique Upwind interpolation for convective Kinetics.
PRESTO : pressure staggering option.
Résumé
Le thème du présent travail vise à étudier le comportement thermique d’une maison

individuelle rafraîchis naturellement (dans notre cas se sont des systèmes de ventilation
naturelle), et à préciser l’influence des systèmes sur la distribution des températures et des
vitesses de l’air dans la maison. Pour cela, nous nous sommes attachés à visualiser les
écoulements d’air et à connaître les champs de vitesse et de température à l’aide d’un code de
calcul Fluent –Gambit (Computational fluid dynamic, CFD). Nous avons choisi pour atteindre ce
but, la modélisation à l’aide du code de champs Fluent, dans l’étude de quatre cas.
Summary
The topic of the present job aims at studying the thermal behavior of a single-family home
cooled of course (in our case, these are systems of natural aeration), and to specify the influence
of systems on the distribution of temperatures and speed of air in the home. For it, we set out to
show the flows of air and to know the fields of speed and temperature with the aid of a code of
counting Fluent – gambit (dynamic Computational fluid, CFD). We chose to attain this purpose,
modeling with the aid of the code of fields Fluent, in the study of four cases.
‫ملخص‬
‫يهدف موضوع هذا العمل إلى دراسة السلوك الحراري للمنزل المنفصل المبرد بشكل طبيعي (في حالتنا أنظمة التهوية‬
.‫ وتوضيح تأثير األنظمة على توزيع درجات الحرارة وسرعات الهواء في المنزل‬،)‫الطبيعية‬
Fluent ‫ ركزنا على تصور تدفقات الهواء ومعرفة مجاالت السرعة ودرجة الحرارة باستخدام كود حساب‬،‫للقيام بذلك‬
‫ في دراسة‬،Fluent ‫ اخترنا النمذجة باستخدام رمز المجال‬،‫( لتحقيق هذا الهدف‬CFD،‫–ديناميكي الموائع الحسابية‬Gambit)
.‫أربع حاالت‬
Introduction générale
La ventilation est un secteur clé du bâtiment, dont le rôle est d’assurer un air sain et
confortable toute l’année, tout en minimisant les consommations énergétiques. Notre travail
consiste à apporter des éléments de réponses et des pistes d’amélioration pour l’élaboration de
stratégies de ventilation, à partir de travaux expérimentaux et numériques.
Une mauvaise maîtrise de la ventilation d'un local aura des conséquences néfastes sur la
qualité de l'air intérieur, la conservation du local et les consommations d'énergie. Dans la vie de
ce dernier, à cause du simple vieillissement ou dans le cadre de réhabilitations ou de
restructurations, les installations de ventilation doivent faire l'objet de modifications.
La ventilation naturelle est un élément clé d’une architecture traditionnelle et aussi un

facteur essentiel pour l’architecture climatique. Elle permet également d’évacuer les odeurs et les
polluants qui s’y accumulent, d’éliminer l’excès d’humidité et de fournir aux appareils à
combustion l’oxygène dont ils ont besoin pour fonctionner sans danger pour notre santé. La
ventilation naturelle est à nouveau au goût du jour, à la fois dans les habitations individuelles et
dans les bâtiments publics. Renouveler l’air n’est pas le seul objectif de la ventilation. Elle est
également utilisée dans les régions arides pour modifier la température intérieure et pour créer
une sensation de fraîcheur dans les bâtiments.
Il existe trois approches permettant d’étudier le problème du cheminement de l'air dans
un local :
L’expérience, la simulation numérique et l'utilisation de formules simplifiées permettant de

déterminer les vitesses et températures d'air dans le local, sachant que ces deux paramètres sont
déterminants pour l'évaluation du confort.
Notre travail penche vers la simulation numérique de l’écoulement et le transfert de chaleur
réalisé à partir des équations de transfert. La méthode numérique est celle des volumes finis.
Objectif du travail
L’objectif du présent travail est premièrement calculer le bilan thermique de la maison,

ensuite la simulation numérique de l’écoulement laminaire et turbulent d’un mouvement de l’air
dans un local.
La maitrise du code fluent précédera tout autre objectif, nous essayerons de nous
familiariser avec le préprocesseur gambit pour la construction d’une géométrie avec génération de
2
Introduction générale
maillage et incorporation des conditions aux limites. Ensuite nous procéderons au code de
simulation proprement dit.
1. Contenu du mémoire
Ce travail est composé de quatre chapitres.

Le premier chapitre contient les généralités sur le confort et les différents types de
ventilation.Le deuxième chapitre comportera le calcul d’un bilan thermique d’un local.
Le troisième chapitre comportera une formulation numérique et présentera le générateur de
maillage GAMBIT et le code de calcul FLUENT.
Dans le dernier chapitre nous évaluerons et discuterons les résultats des simulations trouvé
àl’aide de FLUENT.
Le mémoire se termine par une conclusion générale qui retrace le bilan de cette étude et les
recommandations proposées pour la continuité de cette étude.
3
Chapitre I : Généralités sur le confort
thermique
I.1 Introduction :
Assurer une sensation de chaleur en hiver et se préserver des fortes chaleurs en été est
depuis longtemps un souci majeur pour les concepteurs des bâtiments. D’ailleurs, un des
objectifs de l’architecture réside dans la satisfaction des occupants par le bien être thermique.
Le confort peut être défini comme le degré de désagrément ou de bien-être produit par
les caractéristiques de l'environnement intérieur d’un bâtiment. Une telle définition considère
une interaction entre l'individu et l'espace qui l'entoure, c'est-à-dire, entre des conditions
ambiantes physiquement mesurables et certaines conditions individuelles qui affectent notre
perception.
La qualité de vie à l'intérieur de l'espace a été souvent rapprochée à une appréciation

thermique en premier lieu.
Si un logement n'est pas ventilé, l'oxygène viendrait à manquer, les microbes se
développeraient, l'humidité favoriserait l'apparition de moisissures. Ainsi, une bonne
ventilation est nécessaire pour renouveler l'air mais aussi pour évacuer les polluants toxiques
et autres allergènes. Ce qui va nuire au confort.
Pour garantir une atmosphère saine dans toutes les pièces du bâtiment et préserver le
confortet la santé de ses habitants, deux solutions :
La première est la ventilation naturelle. Zoom sur un procédé économique, écologique
et efficace, qui n'exige aucun élément mécanique consommant de l'électricité à outrance ; et la
ventilation mécanique.
La seconde est la ventilation mécanique (ou artificielle) est destinée à suppléer de façon
totale ou partielle une ventilation absente ou inefficace.
Une ventilation a un rôle qui est l’obtention d'un confort d'été en favorisant les échanges
thermiques convectifs et évaporatifs. L'augmentation du renouvellement d'air permet
d’accroître les échanges avec l'air extérieur et de refroidir le bâtiment lorsque la température
de l’air extérieur est inférieure à celle de l'air intérieur. Le renouvellement d'air doit être limité
quand les températures s’inversent.
I.2 Confort thermique :
I.2.1 Définition
Le confort thermique est une sensation de bien-être lorsqu'on est exposé à une ambiance
intérieure. Le confort thermique ne dépend pas exclusivement de la température, mais
4
thermique
également des conditions d'humidité de l'air intérieur, des éventuels courants d'air, du
niveau de respirabilité de l'air ou de qualité d'air intérieur. La température à elle seule dépend
d'une température résultante sèche, sorte de moyenne des températures intérieures et
rayonnées par les différents corps et parois.
A titre d'exemple, le confort thermique ne peut être atteint à 22°C de température sèche
intérieure avec des murs froids alors qu'à 18°C soit 4 degrés en moins, il est atteint et de plus
de manière plus satisfaisante, avec un rayonnement homogène des corps en présence, y
compris l'individu qui ressent les différents rayonnements. Autre sensation analogue, en
montagne alors que l'air est à +10°C, avec un bon rayonnement du soleil, le confort thermique
peut être atteint facilement.
Le confort c’est tout ce qui contribue au bien-être, a la commodité de la vie matérielle

Bien-être : sensation agréable, procurée par la satisfaction Des besoins physiques
l’absence de contraintes psychologique
Commode : qui se prête aisément et d’une manière appropriée à l’usage qu’on fait. Le
confort thermique a deux composantes le confort (global) et le confort (local) est associe aux
non uniformités proche du corps humain. La difficulté majeure est de caractériser une
sensation humaine par une valeur numérique sachant que le confort thermique est une notion
très complexe faisant intervenir physiologique et psychologique
I.2.2 Les paramètres affectant le confort thermique :

Les paramètres physiques d’ambiance, au nombre de quatre, sont la température de
l’air, latempérature moyenne radiante, la vitesse de l’air, et l’humidité relative de l’air
Les paramètres liés à l’individu, ils sont multiples, on recense notamment deux
paramètres principaux qui sont l’activité et la vêture de l’individu ;
Les Paramètres liés aux gains thermiques internes, gains générés dans l’espace par
des sources internes autres que le système de chauffage. (Éclairages, appareils électriques,
postesinformatiques…etc).
I.2.3 Paramètres liés à l’ambiance extérieure
I.2.3.a La température de l’air ambiant :

La température de l’air, ou température ambiante (Ta), est un paramètre essentiel du
confort thermique. Elle intervient dans l’évaluation du bilan thermique de l’individu au
niveau des échanges convectifs, conductifs et respiratoires. Dans un local, la température de
5
thermique
l’air n’est pas uniforme, des différences de températures d’air se présentent également en plan
à proximité des surfaces froides et des corps de chauffe [1].
I.2.3.b La vitesse de l’air :

La vitesse de l’air joue un grand rôle dans les échanges convectifs et évaporatifs elle
intervient dans la sensation de confort thermique de l’habitant dès qu’elle est supérieure à
0,2m/s [2]. A l’intérieur des bâtiments, on considère généralement que l’impact sur le confort
des occupants est négligeable tant que la vitesse de l’aire ne dépassant pas cette vitesse, sauf
en cas de mauvaise conception du bâtiment ou du système d’aération. Elle peut, en revanche,
être tenue pour responsable de l’apparition d’inconforts locaux, liés à la présence de courants
d’air froids ou chaudslocalisés
I.2.3.c L’humidité relative de l’air :

L’humidité relative est la mesure de la vapeur d’eau contenue dans l’air à une
température donnée. Elle est exprimée en pourcentage de la quantité maximale de vapeur
d’eau que l’air pourrait contenir à la même température. Elle est due essentiellement à
l’occupation humaine, l’activité des habitants et l’humidité extérieure. Une ventilation
adaptée permet de résoudre ce problème. L’inconfort n’apparait que pour une humidité
relative inférieure à 30% ou supérieure à70%. Alors pour un confort optimal et pour une
température de l’air aux environs de22 °C, on peut dès lors recommander que l’humidité
relative soit gardée entre 40 et 65 %.
Fig (I.1): Polygone de confort hygrothermique, température et humidité[28].

1 : Zone à éviter vis-à-vis des problèmes
1 Et 3 : Zones à éviter vis-à-vis des développements de bactéries et de
microchampignons.
2 : Zone à éviter vis-à-vis des développements d'acariens.
4 : Zone de confort thermique. 6
thermique
I.2.4 Paramètres liés à l’individu :

I.2.4.a La vêture :
Le port de vêtements influence les échanges de chaleur avec l’environnement. Les
vêtements permettent de créer un microclimat sous-vestimentaire, à travers leurs résistances
thermiques. Leur rôle essentiel est de maintenir le corps dans des conditions thermiques
acceptables, été comme hiver. La vêture a un rôle primordial d’isolant thermique, notamment
en période hivernale et dans toutes les ambiances froides, ce rôle est pris en compte à travers
la définition d’un indice de vêture, exprimé en caractérisant la résistance thermique d’un
vêtement (figure I-2).
00 0,1 0,3 0,5 0,8 1,0 1,5 3
Fig. (I.2) : valeurs exprimées en Clo des tenues vestimentaires[29].
I.2.4.b L’activité humaine :
L’activité humaine est un paramètre essentiel pour la sensation thermique de l’individu,

définissant directement le métabolisme de l’individu, c’est à dire la quantité de chaleur
produite par le corps humain. Dans le cas d’une très forte activité, elle peut être responsable
de sensations d’inconfort chaud, même en présence de conditions météorologiques très
favorables. Il est à noter toutefois que, dans le cas d’une activité classique de bureau, les
plages de variation du métabolisme demeurent limitées.
I.2.5 Paramètres liés aux gains thermiques internes :
Les apports internes comprennent toute quantité de chaleur générée dans l’espace par 7
thermique
des ressources internes autres que le système de chauffage.
Avec l’essor de la technologie et des besoins électrique (éclairage, électroménager,), les

apports de chaleur internes ont fortement augmenté. Les appareils électriques transforment en
effet quasiment toute l’énergie qu’ils consomment en chaleur, Les postes informatiques sont
également de vraies sources de chaleur et les occupants constituent eux aussi une autre source
d’apports internes par leur métabolisme. Ces gains de chaleur dépendent du type du bâtiment,
du nombre des utilisateurs et de son usage.
D’après Hugues Boivin [4] le confort de l’espace est directement influencé par le taux
de ces gains internes (figure I-3), on peut dire que ces apports sont inévitables dès lors que les
locaux sont habités. Il faut noter cependant que ces apports sont variables selon le
comportement des occupants, et qu’ils constituent donc un facteur d’aggravation de
l’inconfort chaud. Seuls, une bonne ventilation et un comportement adéquat de l’occupant
peuvent réduire ces apports ou leur influence sur la température intérieure [5].
Fig. (I.3) : Gains thermiques internes d’un espace[29].

I.2.6 Facteurs d’inconfort thermique :
Malgré la réalisation d’un confort thermique global, des zones d’inconforts sont
susceptibles d’être observées dans les bâtiments. Un environnement thermique inégal peut
être la source d’inconfort pour certaines parties du corps. L'insatisfaction thermique peut être
causée par un inconfort, causée par un refroidissement ou un réchauffement non désiré d'une
partie du corps (tête, pieds, ou mains), par exemple, un courant d’air.
8
thermique
Un inconfort local peut également être dû à des différences de températures

anormalement élevées entre la tête et la cheville, avec un sol trop chaud ou trop froid, ou à
une asymétrie de rayonnement thermique. Ainsi, le confort thermique peut être affecté par
plusieurs facteurs, citons [6] :
Le courant d’air local, l’asymétrie de la température de rayonnement, la différence

verticale de la température de l’air et la température des planchers.
I.2.6.a Effet des courants d'air :
La perception d'un courant d'air localisé notamment au niveau de la nuque ou du visage

est un élément d'inconfort. Cette perception du courant d'air dépend de la vitesse de l'air, de la
température de l'air, de la zone du corps concernée. Les courants d’air provoquent aussi une
sensation de froid due à une convection assez importante entre la peau et l’air ambiant. Des
courants d’air excessifs en hiver conduisent souvent les occupants à augmenter la température
intérieure pour contrecarrer la sensation de froid. La norme recommande une vitesse d'air
moyenne inférieure à 0,15 m/sen hiveret à 0,25 m/s en été lors d'un travail sédentaire [7].
I.2.6.b Effet de l'asymétrie d'un rayonnement thermique :

Les asymétries du rayonnement sont dues, à la présence d’une paroi chaude ou froide
telle qu’un plafond ou un plancher chauffant, un vitrage chaud ou froid. En outre, le gradient
vertical de température est aussi une source d’inconfort. S’il est suffisamment élevé, il peut
apparaître une sensation de chaud au niveau de la tête ou de froid au niveau des pieds, même
si le corps est en état d’équilibre thermique conclut que l’asymétrie de température radiante
doit être inférieure à 10°C pour une paroi verticale froide (baie vitrée en hiver), et 5°C pour
un plafond chaud (plafond chauffant).
I.2.6.c Effet de gradient thermique vertical de l'air :

En général, les températures sont plus élevées en hauteur donc au niveau de la tête, La
norme admet une différence de température d'air maximum de 3°C entre 0,1 m du sol (niveau
deschevilles) et 1,1 m du sol (niveau de la tête chez une personne assise) [8].
9
thermique
I.2.6.d Effet de la température du sol :
Hoffman J B [9]. Précise qu’une température de plancher trop élevée ou trop basse
entraîne un inconfort au niveau des pieds. Plusieurs auteurs ont effectué des recherches sur ce
sujet et selon Olsen BW [10], les températures optimales de sol pour les personnes chaussées
et à la neutralité thermique sont de 23 °c pour les personnes debout et de 25°c pour les
personnes assises, avec un minimum de 6% d’insatisfaits.
D'autres éléments influençant le confort thermique ont été étudiés : l'âge, le sexe, la
nourriture, la localisation géographique, la couleur des murs, le bruit, la lumière. Toutefois, il
ne peut exister de règles de confort universelles du fait de grandes variations géographiques et
intra- individuelles [11], sauf que Les médecins de l’habitat sain proposent les valeurs
suivantes pour chaque facteur du confort [12] :
* Température des murs : 22 ±2°C
* Humidité relative entre 30 et70%
* Température du sol : 19 à24°C
* Vitesse de l’air : inférieure à 0,20 m/s ;
* Différence de température entre deux murs d’une même pièce doit être inférieure à
10°C
* Différence de température entre le sol et le plafond doit être inférieure à 5°C.
I.2.7 L’évaluation de confort thermique :

L’évaluation du confort thermique dans les espaces est un paramètre capital dans toute
conception [13]. Les premières recherches se sont basées sur les enquêtes de terrain avec des
questionnaires en classifiant la sensation thermique (très chaud, neutre et très froid) ainsi que
sur les essais de laboratoires sous des conditions climatiques artificielles [14].
Cette évaluation a conduit plusieurs chercheurs à développer et à élaborer des

méthodesde prédiction des niveaux de confort à l’intérieur des bâtiments.
La méthode de Steven Szokolay :
La méthode de S. Szokolay consiste à définir les zones de contrôle potentiel, projetées

parétapes sur le diagramme psychrométrique.
La zone de confort : 10
thermique
Définir par un point neutre, obtenu par l’intersection de la température en entre et

l’humidité relative portée à 50% d’humidité relative.et délimité par quatre points de
températures humides.
La zone de contrôle potentiel pour les périodes froides : C’est la zone de chauffage
passif, délimitée par un seul point de la température humide porte à la courbe de saturation
HB=100% où :
La zone de contrôle potentiel pour les périodes chaudes :
Zone d’effet de masse : La détermination de la variation moyenne de la température

pour le mois le plus chaud permet d’obtenir le premier point de cette zone.
Zone d’effet de masse avec ventilation nocturne :
Zone de ventilation naturelle-mouvement d’air sur la surface de la peau :
Zone de refroidissement évaporatif direct : Les limites supérieures et inférieures de

cette zone sont tracées à partir des lignes de températures humides du deuxième et troisième
point de la zone déconfort.
Zone de refroidissement évaporatifs indirect : Délimitée par un seul point ou la

limite supérieure de cette zone correspond à la ligne 14g/Kg d’humidité absolue.
Représentation des conditions climatiques qui prévalent :
Chaque mois est représenté par une ligne dessinée entre deux points :
 La température moyenne maximale mensuelle avec l’humidité relative

minimale.
 La température moyenne minimale mensuelle avec l’humidité relative

maximale.
11
thermique
Fig (I.4) : diagramme psychrométrique [29].
12
thermique
I.3 La ventilation :
I.3.1 Définition :
La ventilation est l'action qui consiste à créer un renouvellement de l'air, par
déplacement dans un lieu clos. Elle est mise en œuvre dans les lieux où l'oxygène risque de
manquer, ou bien où des polluants et autres substances indésirables (humidité, par exemple)
risqueraient de s'accumuler en son absence : logements, bureaux, magasins, salles de
spectacles, d'enseignement, ouvrages souterrains, tunnel routier, atelier industriel, mine….etc.
La ventilation a pour vocation d’évacuer l’air vicié des logements en le renouvelant

par del’air frais. La ventilation mécanique désigne tous les dispositifs motorisés d’évacuation
ou d’insufflationd’air frais.
I.3.2 Principaux types de ventilation d’une maison :
Dans des maisons de mieux en mieux isolées, choisir un système de ventilation

performant est indispensable, en particulier dans les pièces humides. Ventilation mécanique
contrôlée (VMC) simple flux ou double flux, système thermodynamique, hygroréglable ou
non, il n’est pas toujours facile de s’y retrouver parmi les différents types de ventilations
proposés.
I.3.2.1 La ventilation naturelle :
La ventilation naturelle était utilisée avant la généralisation des systèmes de ventilation

mécanique contrôlé VMC. Elle utilise le principe du tirage naturel ou les dépressions d’air
créées par le vent. Avec un système de ventilation naturelle, l’air extérieur entre dans la
maison à travers des grilles d’aération situées en bas des logements, dans les pièces
principales de l’habitation.
L’air circule dans le logement et monte au fur et à mesure qu’il se réchauffe,
pour êtreévacué par des grilles d’aération placées en hauteur.
Si un logement n'est pas ventilé, l'oxygène viendrait à manquer, les microbes se
développeraient, l'humidité favoriserait l'apparition de moisissures. Ainsi, une bonne
ventilation est nécessaire pour renouveler l'air mais aussi pour évacuer les polluants toxiques
et autres allergènes. Pour garantir une atmosphère saine dans toutes les pièces du bâtiment et
préserver la santé de ses habitants, une solution existe : la ventilation naturelle, un procédé
économique, écologique et efficace, qui n'exige aucun élément mécanique consommant de
l'électricité à outrance.
13
thermique
I.3.2.1.a Principe de fonctionnement de la ventilation naturelle :
La ventilation naturelle se sert des différences de pression provoquées par le vent afin
d'assurer la circulation de l'air d'une pièce à l'autre en continu. Sachant que l'air chaud est plus
léger que lorsqu'il est froid, il monte et génère un balayage permanent dans le bâtiment. L'air
extérieur pénètre tout simplement par des grilles d'arrivée d'air installées dans les pièces de
vie, puis passe dans tout le logement grâce à des grilles de transfert. Les pièces humides et
plus chaudes - la cuisine, les WC et la salle de bains - sont équipées de bouches d'extraction,
de sorte que l'air vicié puisse retourner dehors. Il passe alors par des conduits à tirage naturel
coiffés d'extracteurs statiques, qui ressortent par le toit. Aucun ventilateur n'intervient, ce qui
représente un atoutconsidérable en termes de consommation électrique.
Fig. I.5 la ventilation dans une maison individuelle [30].
I.3.2.1.b Les différents types de ventilation naturelle :
b.1 Le tirage thermique (effet de cheminée) :
L’effet de cheminée est expliqué par la poussée d’Archimède qui fait monter dans la
cheminée l’air chaud, dont la densité (et la pression) est plus petite que celle de l’air froid. Le
tirage de la cheminée est d’autant plus fort que la différence thermique est elevée.il permet
14
thermique
d'évacuer l'air chaud selon le principe de la convection. En effet, celui-ci remonte

naturellement pour s'échapper en partie haute. Il crée simultanément une dépression en partie
basse de manière proportionnelle à la différence de température entre l'intérieur du logement
et l’extérieur. Dans les pièces humides, il convient donc d'installer des grilles d'aération en
hauteur pour permettre l'évacuation de l'air intérieur et en bas pour favoriser l'arrivée de l'air
frais de l'extérieur. Combiné à un extracteur, l'effet de tirage du conduit vertical s'en trouve
renforcé car, en tournant avec la force du vent, il produit une dépression supplémentaire
[15].
Fig. (I.6) : Tirage thermique [31].

b.2 La ventilation traversant :
Ce type de ventilation consiste à tirer une partie du vent. Pour en disposer, il convient de
ménager des entrées et des sorties d’air dans le logement. Faire appel à un professionnel qui
vous conseillera sur le choix et la disposition de la ventilation est la meilleure manière de tirer
parti de cette ventilation naturelle. Pour ce faire :
– Il faut installer des entrées d’air face au vent dominant.
– Il faut installer des sorties d’air à l’opposé.
De cette manière, le vent crée une pression qui permet de faire entrer l’air extérieur
d’un côté, tandis qu’il crée une dépression qui chasse l’air chaud et pollué de l’autre.
La ventilation traversant engendre ainsi un courant d’air qui permet le renouvellement
de l’air intérieur. Cette méthode nécessite cependant un agencement spécifique requérant,
idéalement, des pièces de vie situées du côté du vent dominant.
Fig. (I.7) : Ventilation Traversant[31].

15
thermique
b.3 Ventilation naturelle mono-orientée :

S’il n’y a qu’une seule ouverture dans le local, un mouvement d’air se crée sous le
double effet de la pression extérieure du vent et de la différence de température entre
l’intérieur et l’extérieur. Ceci crée une surpression dans la partie basse de la fenêtre et une
dépression dans la partie haute. Pour une même vitesse de vent, ce moteur est cinq à six
fois plus faible que celuid’une ventilation traversant.
Fig (I.8) : Ventilation naturelle mono-orientée[31].
b.4 Ventilation naturelle hybride :
Ce type de ventilation permet d'assurer une aération optimale du domicile. Pour cela,
une assistance mécanique à basse pression est ajoutée. Le système prend le relais dès l'instant
que les conditions de tirage naturel ne sont plus réunies. Ce type de ventilation est donc
surtout utile durant les périodes de froid ou de canicule, car le débit d'air est stabilisé et, de ce
fait, les déperditions énergétiques grandement limitées [16].
Fig. (I.9) : ventilation naturelle hybride[32].
16
thermique
I.3.2.1.c Avantages et inconvénients de la ventilation naturelle :
La ventilation naturelle possède de multiples atouts :

• Nul besoin de moteur pour fonctionner ! Elle n'entraîne donc pas de
consommationénergétique ni de nuisance sonore ;
• La pression du vent et les différences de température constitue le seul moteur, ce qui

réduitvotre impact sur l'environnement ;
 Il n'y a presque pas d'entretien à faire ;
 Il s'agit d'une solution économique. Les grilles d'aération et les bouches

d'extraction affichent des prix très abordables dans les magasins de bricolage ou encore sur
Internet ;
En rénovation, la ventilation naturelle est souvent privilégiée puisqu'elle évite de

s'encombreravec des gaines de pulsion et une unité double flux indépendante
Mais ce concept révèle aussi quelques limites :
 La ventilation naturelle demeure tributaire de la météo. En conséquence, son

débit varie en permanence. Pour pouvoir ainsi profiter du tirage thermique, la température
extérieure doit être inférieure à celle au sein du logement ce qui est rarement le cas durant la
saison estivale. Dans ce contexte, l'aération est insuffisante, ce qui risque de causer de
l'humidité et de la condensation ;
 La ventilation naturelle ne contrôle pas les quantités d'air qui transitent au sein
du bâtiment. Ceci s'explique par les différences de températures entre dedans et dehors, les
mouvements du vent et la pression atmosphérique. En hiver, le risque de sur-ventilation est
grand et, donc, de déperditionthermique ;
 Certes le tirage naturel se réalise aisément en hiver puisque l'intérieur est plus
chaud. Cependant, en été, il peut littéralement s'inverser, l'air pouvant se déplacer des pièces
humides vers les espaces de vie ;
 Pour une construction neuve, elle peut ne pas suffire. En effet, l'isolation est plus
performante que dans l'ancien. Résultat : la ventilation naturelle risque de donner lieu à des
pertesde chaleur en hiver et de la condensation en été ;
 La qualité de l'air intérieure n'est pas forcément garantie. Le vent peut
effectivement compromettre le renouvellement de l'air selon le sens où il souffle, des fenêtres
ouvertes peuvent également perturber le bon fonctionnement de la ventilation, et ce n'est sans 17
thermique
parler des ouvertures qui laisser passer les bruits ainsi que la pollution. Dans un
environnement urbain, il est donc indispensable de prévoir un traitement acoustique des
grilles ;
 Les grilles d'ouverture peuvent être source d'inconfort, principalement en hiver, à
moins qu'elles soient placées suffisamment en hauteur pour ne pas ressentir l'air froid ;
 Les ouvertures dans les façades ne sont pas forcément très esthétiques.
I.3.2.2La ventilation mécanique
Ventilation assurée par une installation mécanique qui permet l'évacuation ou

l'insufflation forcéed'air frais.
I.3.2.1.d Les différents systèmes de ventilation mécanique :
a.1 Ventilation mécanique contrôlé (VMC) :

Ensemble des dispositifs mécaniques assurant le renouvellement de l’air dans les pièces,
notamment les pièces humides (cuisine, salle de bain, toilettes) par une ventilation minimum
et permanente. Intégrée au bâtiment, la VMC fonctionne avec une centrale de ventilation
forçant l'extraction de l'air et assurant ainsi la qualité de l'air intérieur.
Fig. (I.10) : Ventilation mécanique contrôlé[32].
a.2 Ventilation à simple flux :

Le système de ventilation à simple flux est le plus utilisé actuellement. Il comporte des
bouches d’extraction auto réglables reliées à un moteur. Ce dernier assure l’extraction de l’air
chargé d’humidité provenant de l’intérieur de l’habitation. L’alimentation en air extérieur
se fait la plupart du temps grâce à des entrées d’air situées sur les fenêtres. Peu onéreuse et
facile à installer, la VMC simple flux fonctionne en continu et peut provoquer des courants
d’air. Elle est parfois un peu bruyante et occasionne des déperditions de chauffage.
18
thermique
Fig. (I.11) : Ventilation à simple flux[32].

La VMC hygroréglable :
La VMC hygroréglable règle son débit d’air en fonction de l’hygrométrie de son

environnement. Les déperditions d’énergie sont ainsi moins importantes qu’en
utilisant une VMC simple flux classique, car le flux d’air est ajusté en permanence.
Il existe 2 types de VMC hygroréglables :
a.2.1 La VMC hygroréglable de type A : seules ses bouches d’aération

s’ajustent au taux d’humidité. Ses entrées d’air sont semblables à celles d’une VMC
classique à simple flux. L’apport d’air est donc constant même si l’extraction d’air est
régulée.
a.2.2 La VMC hygroréglable de type B : à la différence de celle de type A, ses
entrées d’air sont également hygroréglables permettant une régulation plus performante
de son débit. Lesdéprédations thermiques sont aussi plus faibles.
Fig. (I.12) : ventilation mécanique contrôle par hygroréglables[32].
19
thermique
a.3 La ventilation à double flux :

La ventilation à double flux est constituée de deux réseaux de circulation d’air. Elle
ne comporte pas plusieurs entrées d’air disposées dans les pièces principales, mais une
seule, installée généralement dans les combles. Son grand avantage est de permettre
de récupérer les calories de l’air vicié au moyen d’un échangeur thermique, sans que les
deux circuits ne se mélangent, afin de réchauffer l’air frais arrivant de l’extérieur. Ainsi,
elle permet d’éviter tout courant d’air et de limiter les déperditions d’énergie dans le
logement. La VMC double flux est un système plus complexe à installer et plus
encombrant. Son débit d’air peut être ajusté pour chaque pièce, afin d’assurer une
efficacité maximale. Dotée d’un système inversé, elle peut aussi rafraîchir un logement
en été. Elle est très souvent hygroréglable.
Fig. (I.13) : ventilation mécanique double flux[33].

Actuellement l’un des systèmes les plus intéressants sur le marché, la VMC
thermodynamique est l’association d’une VMC double flux avec une pompe à chaleur
géothermique. Elle permet de bénéficier d’une aération efficace, ainsi que d’un un
système de chauffage, car elle restitue les calories récupérées dans le sol pour chauffer la
maison. Elle est réversible et peut aussi rafraîchir le logement en été. Cette solution
utilisant une énergie renouvelable et gratuite suffit dans certains cas à chauffer jusqu’à
90 % une maison passive.
20
thermique
Fig. (I.14) : VMC double flux thermodynamique[34].
a.4 La ventilation mécanique par Insufflation (VMI) :

La ventilation mécanique par insufflation (VMI) inverse le principe de
fonctionnement dela VMC. Plutôt que d’aspirer l’air frais, elle le projette vers l’intérieur
de l’habitation, tout en le filtrant et en le préchauffant. Elle ne nécessite pas de
passage de gaines, ce qui est un avantage en rénovation lorsque de gros travaux ne
sont pas envisageables. Elle convient particulièrement aux personnes sensibles aux
allergies, car elle assèche l’air trop humide en le préchauffant et elle contient des filtres
très performants. En revanche elle est un peu plus chère que les autres systèmes et
encore peu répandue.
Fig. (1.15) : Ventilation simple flux par insufflation[34].

21
thermique
a.5 La VMC gaz :

Ce type de ventilation existe surtout dans les logements collectifs, chauffés avec un
ancien modèle de chaudière à gaz. Le fonctionnement de la VMC gaz ressemble à celui
de la VMC simple flux. Cependant, elle évacue également les gaz de combustion issus d’un
chauffe-eau au gaz ou d’une chaudière. Son utilisation est soumise à une réglementation
très stricte et, raccordée à la chaudière, elle peut stopper son activité en cas de problème au
niveau de l’extracteur d’air.
Fig. (I.16) : ventilation mécanique contrôle au gaz[34].
a.6 La ventilation mécanique répartie :

La ventilation mécanique répartie fonctionne comme une VMC : elle comporte des
entrées d’air disposées en haut des ouvertures du logement. Mais à l’inverse des systèmes de
VMC, elle extrait l’air à partir d’aérateurs indépendants, situés dans chaque pièce. Ce
dispositif est à réserver à des logements anciens, dans lequel passer des gaines demanderait
des travaux trop importants. Elle ne respecte pas les préconisations actuelles en matière
d’économies d’énergie.
Fig. (I.17) : ventilation mécanique repartie[34].

22
thermique
I.3.3 Principe d’écoulement :

L’écoulement est provoqué par la différence de masse volumique entre l’air intérieur
(chaud) et l’air extérieur (froid). L’air chaud possède une masse volumique plus faible que
l’airfroid et tend à monter. L’air froid est plus lourd et tend à rester au sol. L’air chaud et
humide s’élève vers les cheminées et est remplacé par de l’air froid et sec. Les grandes
surfaces des entrées de l’air assurent de faibles vitesses qui ne causent pas de stress aux
occupants.
Système de Ventilation Ventilation mixte Ventilation mixte Ventilation

ventilation naturelle Entrée naturelle Entrée mécanique mécanique
Caractéristiques Sortie mécanique Sortie naturelle
Principales
Utilisation des Utilisation courante Utilisation des Utilisation
forces convectives pour des bâtiments forcesconvectives générale
Domaine d’application existantes ateliers relativement bas existantes
hautset étroits
Possibilité de contrôle dela Non Non oui oui
distribution spatiale de l’air
introduit
Possibilité de contrôle dela Non Non oui oui
qualité (température humidité,
pureté) de l’airIntroduit
Possibilité de contrôle de la Non Non Effet thermique : Oui Oui
pression à l’intérieur du (dépression) (dépression) Sans effet thermique : (surpression
bâtiment (surpression où Non ou
Dépression dépression)
Possibilité de récupérationde Non Oui Non Oui
chaleur sur l’air extrait
Indépendance vis-à-vis du
vent
- Des entrées d’air Non Non Oui Oui
- Des sorties d’air Non Oui Non Oui
Problèmes particuliers Existence Existence de Solution
de courants d’air courants d’air onéreuse
pourles gros
débits
d’air
Tableau (I.1) : Comparaison des différents systèmes de ventilation générale. [17]
I.3.3.1 Critères d’évaluation de l’efficacité de ventilation :

Les normes et réglementations en vigueur proposent des débits recommandés pour
satisfaire une bonne qualité d’air, mais ne qualifie pas réellement la performance de la
ventilation. L’efficacité de la ventilation peut être exprimée à l’aide d’indices représentatifs de
23
thermique
sa capacité à renouveler l’air et à éliminer les polluants.
I.3.3.2 Efficacité du renouvellement d’air :
L’efficacité de renouvellement de l’air intérieur d’un système de ventilation caractérise

la rapidité à laquelle l’air d’une pièce est remplacé en comparaison avec le renouvellement
d’air le plus rapide. On peut déterminer ce terme expérimentalement à l’aide de méthodes
de gaz traceurs, mais il reste difficile à évaluer en cas de ventilation naturelle ou de bâtiments
trèsperméables, du fait de la forte variabilité du renouvellement d’air.
Les différents indices d’efficacité du renouvellement d’air sont construits à partir des
temps caractéristiques de l’air, que sont l’âge de l’air, le temps résiduel et le temps de
résidence :
 L’âge de l’air :
C’est le temps écoulé depuis l’entrée dans un local d’une particule de l’air. Il caractérise
la « fraîcheur de l’air » et s’exprime différemment suivant la méthode de gaz traceur
utilisée et selon si on s’intéresse localement à des zones d’intérêt particulières ou de manière
globale dans l’enceinte.
 Le temps résiduel :
Il s’agit du temps mis par une molécule d’air quittant un point P donné du domaine
étudié pour atteindre la sortie d’air.
 Le temps de résidence et le temps nominal :

Le temps de résidence est le temps que met une molécule d’air venant du soufflage (ou
de l’extérieur) et passant un point P donné du domaine étudié pour atteindre la bouche
d’extraction de ce domaine. Il correspond à la somme de l’âge de l’air et du temps résiduel de
l’air.
Un temps de résidence prolongé peut entraîner un risque de mauvaise qualité de l’air

intérieur. On définit le temps de résidence minimum de l’air ou temps nominal, comme étant
l’inverse du taux de renouvellement d’air.
 Indice d’efficacité de renouvellement d’air :

Selon le type de ventilation, cet indice sera compris entre 0 et 1, il correspondant au
minimum de Le rendement de ventilation est considéré efficace s’il est compris entre 0,5 et
1,0.
24
thermique
Les meilleures efficacités de renouvellement d’air sont obtenues pour une ventilation par
piston où l’air vicié est déplacé par l’air frais, généralement de bas en haut, l’air soufflé ne se
mélangeant ainsi pas à l’air de la pièce (Figure 1.18). Dans le cas d’une ventilation par
mélange consistant à assurer une homogénéité parfaite de l’air frais avec l’air vicié
l’efficacité n’est que de 50% car le temps de résidence de l’air dans le domaine est égal au
double du temps nominal de renouvellement d’air.
Enfin, si le système de ventilation génère des mouvements d’air avec des zones mortes
provoquées par des phénomènes de recirculation ou des court-circuit entre les bouches de
soufflage et les grilles d’extraction, le rendement peut descendre en dessous de 50%. [18]
Fig I.18 : Modes de ventilation, (a) par piston, (b) par mélange, et (c) en court-circuit [34].
I.3.3.3 Taux de renouvellement d’air :

En ventilation naturelle, on définit le taux de changement d’air (ACH), qui indique le
nombre de fois en une heure où la totalité du volume d’air d’un espace est renouvelé. Ce
dernier est exprimé par le rapport de débit volumique et le volume de la pièce considérer, ce
rapport est l’ACH. Cet indice donné par ASHARE comme [19]
60 × 𝑄
𝐴𝐶𝑃𝐻 =
𝑉𝑜𝑙
Où :
ACPH = nombre de changements d'air par heure ; des valeurs plus élevées
correspondent à une meilleure ventilation.
Q = débit volumétrique d'air en pieds cubes par minute (cfm).
Vol = volume d’espace L × l × H, en pieds cubes. 25
thermique
3.6 × 𝑄
𝐴𝐶𝐻 =
𝑉𝑜𝑙
Où :
ACPH = nombre de changements d'air par heure ; des valeurs plus élevées
correspondent à une meilleure ventilation.
Q = débit volumétrique d'air en litres par seconde (L / S).
Vol = volume d’espace L × W × H, en mètre cube.
Conclusion :
Exposé à un environnement thermique, le corps humain réagit par une interaction

dynamique mobilisant un ensemble de réactions rétroactives, volontaires et involontaires, qui
permettent de contrôler les échanges thermiques avec cet environnement afin de satisfaire son
confort thermique.
L’environnent thermique est caractérisé par quatre grandeur physiques (la

température de l’air, la température due au rayonnement, l’humidité et la vitesse de l’air). Ces
variables réagissent avec l’activité et la vêture du corps humain pour établir son état
thermique et constituent ensemble les six paramètres de base de l’échange thermique entre
l’homme et son environnement. Mais au-delà de ces variables, la perception thermique d’un
environnement peut être influencée par des variables physiologiques, psychologiques et
sociologiques.
L'étude et l'évaluation du confort thermique permettent de classe quatre (4) zones :

 Inconfort : expression de l'insatisfaction ou du désagrément de façon
continue ou répétitive, ne permettant pas à la nuisance de se "faire oublier » ;
 Confort dégradé : résulte de l'apparition épisodique (temporelle) ou locale

(spatiale) d'un certain inconfort non persistant ou ressenti localement mais de façon peu
intense (par exemple, avoir légèrement froid aux pieds) ;
 Confort : peut résulter de 2 possibilités :

Se déduit de l'absence exprimée de désagrément ou d'insatisfaction (échelle affective);
Se déduite l'incapacité qu'ont les gens à s'exprimer sur leur état thermique (Échelle
perceptive) ; cette notion est proche de l'indifférence.
 Confort optimal : résulte de l'expression de l'agrément du bien être exprimé
par l'individu par rapport au climat perçu (l'exemple le plus convaincant est celui de la 26
thermique
douche :on ajuste la température de l'eau pour une sensation d'agrémenter cherchée).
Le système de ventilation a une importance capitale sur l’amélioration de la qualité de
l’air intérieur.
L’aération, qu’elle soit réalisée de façon naturelle ou mécanique, est indispensable pour
obtenir un air ambiant de bonne qualité. Elle a pour priorité de satisfaire les besoins d’hygiène
et de confort et vise d’une part à apporter de l’air neuf, afin d’éviter les situations de
confinement et de permettre le bon fonctionnement des appareils à combustion ; d’autre part à
évacuer l’air vicié par différents polluants, et à lutter contre l’humidité et les condensations.
27
CHAPITRE II : Calcul des apports thermiques
Ⅱ.1 Introduction
Le (DTR) est un document technique réglementaire qui a pour objet de fixer les
méthodes de détermination des apports calorifique et des déperditions thermique, vérification
de la conformité à la règlementation thermique d’été et d’hiver des bâtiments.
Les méthodes de détermination des apports calorifique et des déperditions thermiques

du (DTR) s’appliquent aux locaux à usage : d’habitation, d’hébergement, de bureaux,
d’accueil, de réunion, de restauration.
Ⅱ.2 Les apports caloriques
Ⅱ.2.1 Définition
Les apports calorifiques (appelés aussi gains) d’un local sont égaux à la somme des
apports de chaleur sensible et latente, provenant d’une source intérieure ou extérieur du local,
pour des conditions intérieures et extérieures déterminées, et ne tenant pas compte des apports
dus à l’installation.
Ⅱ.2.2 Calcul des apports calorifiques
Les apports calorifiques doivent être déterminés selon les étapes suivantes :
 Définition des zones thermiques,
 Détermination de l’intervalle de temps critique,
 Calcul des apports calorifiques pour chaque volume thermique et pour toutes les
heures situées dans l’intervalle du temps critique,
 Calcul des apports calorifiques effectifs,
 Calcul de la puissance frigorifique,
 Adoption d’un système de conditionnement d’air.
Ⅱ.2.3 Les apports calorifiques sensibles As et latents Al (W) :
Ils sont donnés par :
- As = APO+ AV +AIs +AINFs (II.1)
- Al = AIl + AINFl (II.2)
28
- APO (en W) représente les apports par les parois opaques,
- AV (en W) représente les apports à travers les parois vitrées,
- AIs et AI1 (en W) représentent les parties sensibles et latentes des apports
internes,
- AINFs et AINF1 (en W) représentent les parties sensibles et latentes des

apports dus aux infiltrations d’air,
Ⅱ.2.4 Les apports calorifiques effectifs sensibles en (W) AEs et latents AEl :
ils sont donnés par :
AEs = (CΔas × As) + (BF × ARENs) (II.3)
AEl = (CΔal × Al) + (BF × ARENl) (II.4)
- As(en W) représente les gains sensibles,
- Al(en W) représente les gains latents,
- ARENset ARENi(en W) représentent les parties sensibles et latentes des apports

dus à la ventilation des locaux,
- BF ou facteur de by-pass exprime la partie de l’air extérieur (air neuf) non traité
par l’installation de climatisation (imperfection de l’appareil de traitement), et qui parvient au
local sans modifications (BF=0.40)
- CΔas est un coefficient de majoration des gains sensibles qui prend en compte les
gains supplémentaires (Échauffement du ventilateur, réseau de conduits d’air traversant des
locaux non conditionnés)
- CΔallest un coefficient de majoration des gains latents qui prend en compte les
gains supplémentaires (tels que ceux dus aux fuites d’air éventuelles dans les réseaux de
conduits d’air)
Ⅱ.2.5 Les apports calorifiques totaux sensibles ATs et latents ATi :
Ils sont donnés par :

ATS = (CΔAS ×AS) + ARENS (II.5)
29
ATI = (CΔAL ×AI) + ARENL (II.6)
Ⅱ.2.6 Vérification réglementaire

La somme des apports calorifiques par les parois vitrées et les parois opaques aériennes
doit vérifier au mois de Juillet à 15 h TSV, pour une température sèche intérieure de 27°C, la
relation ci-après :
APO (15 h) + AV (15 h) ≤ 1.05 × AREF (15 h) (II.7)
- APO (en W) désigne les apports calorifiques à travers les parois opaques
aériennes,
- AV (en W) désigne les apports calorifiques à travers les parois vitrées,
- AREF (en W) désigne les apports calorifiques de référence, Les apports
calorifiques de référence AREF sont donnés par :
AREF = AREF, PH+ AREF, PV + AREF, PVI (II.8)
- AREF, PH (en W) désigne les apports calorifiques de référence à travers les parois
opaques horizontales,
- AREF, PV (en W) désigne les apports calorifiques de référence à travers les parois
opaques verticales,
- AREF,PVI (en W) désigne les apports calorifiques de référence à travers les parois
vitrées,
Ⅱ.2.7 Apports de référence des parois opaques horizontales

Les apports calorifiques de référence des parois opaques horizontales AREF,PH sont
calculés par la formule suivante :
AREF,PH = Σ (a × SINT × ΔTSREF,PH) (II.9)
- a (en W/m².°C) est un coefficient lié à la nature de la construction et fonction de

la zone climatique;
- Sint (en m²) désigne la surface de la paroi horizontale comptée de l’intérieur ;
- ΔTSREF, PH (en °C) est la différence de température de référence pour les parois
horizontales.
30
La valeur du coefficient a :
La valeur de a est égale à :
 1,90 W/m². °C pour les planchers hauts (plafonds) en contact avec un local non
conditionné (quelle que soit la zone climatique) ;
 2,70 W/m². °C pour les planchers bas en contact avec un local non conditionné (quelle
que soit la zone climatique) ;
 Donnée dans le tableau 1 pour les planchers hauts en contact avec l’extérieur
Tableau II.1 : valeurs de a
Valeurs de a (en W/m².°C)
Zone Logement Logement en immeuble

Climatique Individuel collectif, bureau
A 1.10 1.10
B 1.10 0.9
B’ 1.10 0.9
C 1.10 0.85
Ⅱ.2.8 Apports de référence des parois opaques verticales

Les apports calorifiques de référence des parois opaques verticales Aréf, PV sont
calculés par la formule suivante :
AREF,PV = Σ (c × Sint × ΔTSREF,PV) (II.10)
- c (en W/m².°C) est un coefficient lié à la nature de la construction et fonction de la

zone climatique,
- Sint (en m²) désigne la surface de la paroi verticale comptée de l’intérieur,
- ΔTS REF, PVI (en °C) est la différence de température de référence pour la paroi opaque
verticale
Ⅱ.2.9 Apports de référence des parois vitrées

Les apports de référence pour les parois vitrées sont donnés par la somme des apports
de référence de chaque paroi vitrée.
31
 Les apports de référence AREF, PVI pour une paroi vitrée en contact avec l’extérieur sont
donnés par :
AREF, PVI = AVEREF + AVTREF (II.11)
- AVE REF (en W) désigne les apports de référence dus à l’ensoleillement,
- AVT REF (en W) désigne les apports de référence dus au gradient de température. Les
apports de référence dus à l’ensoleillement à travers une paroi vitrée AVEREF sont donnés par
:
AVE REF = [SVENS× It + (SV – SVENS) × Id] FS REF× NPVI REF (II.12)
- SV (en m²) est la surface totale vitrée,
- SVENS (en m²) est la surface vitrée ensoleillée à 15 h TSV
- It (en W/m²) est le rayonnement total maximal réel pour l’orientation et la

latitude considérées,
- Id (en W/m²) est le rayonnement diffus maximal réel,
- FSREF est un facteur solaire de référence qui égal à :
- 0,15 pour les locaux à usage d’habitation et d’hébergement,
- 0,38 pour les locaux à usage de bureaux.
- NPVIREF est le coefficient d’amortissement relatif aux gains de référence des

parois vitrées pour l’orientation considérée ; les valeurs de NPVIREF
Ⅱ.2.10 Apports à travers les parois opaques

Ⅱ.2.10.a Parois extérieures
Les parois opaques extérieures sont celles qui sont en contact direct avec l’air extérieur
(Verticales ou horizontales).
Les apports de chaleurs à travers une paroi opaque à un instant t, APO(t), sont donnés
par la formule suivante :
APO(t) = 1,2×Kété× SINT ×Δte(t) (II.13)
- 1,2 (sans dimension) est un coefficient de majoration tenant compte des apports
latéraux linéiques (à travers les ponts thermiques),
32
- Kété (en W/m². °C) est le coefficient de transmission en partie courante de la
paroi considérée pour l’été,
- SINT (en m²) est la surface intérieure totale de la paroi considérée ; pour les
toitures en pente, on prendra la projection horizontale de la surface ;
- Δte(t) (en °C) est la différence équivalente de température à l’heure t ;
Ⅱ.2.10.b Le coefficient Kété des parois opaques est donné par la formule suivante :
1
- Kété = ∑ (II.14)
𝑅
- ∑R (en m². °C/W) représente la somme des résistances thermiques des

différentes couches de matériaux constituant la paroi.
Ⅱ.2.10.c La différence de température équivalente Δte(t) en °C est donnée par : Paroi

ensoleillé :
Δte(t) = Δtes(t) + CΔte + [Δtem(t) - Δtes(t)] × (II.15)

.
Ⅱ.2.10.d Paroi à l’ombre 24 h / 24
Δte(t) = Δtes(t) + CΔte (II.16)
- Δ tes(t) (en °C) est la différence de température équivalente à l’heure t en considérant

que la paroi est à l’ombre ;
- CΔ te (en °C) est un facteur de correction dû, d’une part, à la différence maximale
ΔTSmax de la température sèche de l’air entre l’extérieur et l’intérieur pour le mois
Considéré, et d’autre part à l’écart diurne de base Eb pour le mois considéré ;
- Δ tem(t) (en °C) est la différence équivalente de température à l’heure t pour
L’orientation de la paroi considérée ;
- It,b (en W/m²) est le rayonnement total de base pour le mois,

- It,b(40) (en W/m²) est le rayonnement total de base pour le mois de Juillet,
Le facteur d’absorption α des parois extérieures est donné ci-après selon que la paroi est
peinte ou non.
33
Le facteur d’absorption α d’une paroi extérieure peinte est donné en fonction de la
couleur de la face extérieure de la paroi :
- Face extérieure de couleur sombre (bleu foncé, rouge foncé, brun foncé), α = 0,90 ;
- Face extérieure de couleur moyenne (vert clair, bleu clair, gris clair), α = 0,70 ;
- Face extérieure de couleur claire (blanc, crème), α = 0,50 ;
- Face extérieure de couleur noire mat, α = 1.
- Le facteur d’absorption α d’une paroi extérieure non peinte est donné par le tableau
ci- dessous en fonction du matériau constituant la face extérieure de la paroi.
Tableau II.2 : facteurs d’absorption

Nature du matériau Facteur d’absorption Nature du matériau Facteur
d’absorption
Asphalte 0.90 pierre calcaire blanche 0.60
feutre bitumé 0.90 granit poli 0.45

SABLE 0.70 acier poli 0.06
ARDOISE 0.93 Cuivre poli 0.04
TUILE 0.80 Aluminium poli 0.04
carrelage claire 0.30 Tôle galvanisée 0.70
marbre poli 0.60 Amiante ciment 0.96
tôle oxydée 0.80 Peinture d’aluminium 0.60
tôle rugueuse 0.94 Fonte brute 0.80
béton, brique 0.75 Gravier 0.80
Ⅱ.2.10.e Parois intérieures

Les apports de chaleur à un instant t, APO(t), traversant une paroi opaque en contact
avec deux locaux conditionnés sont donnés par la formule suivante :
APO(t) = Kété × SINT × [TSA – TSB ,I] (II.17)
- Kété (en W/m².°C) est le coefficient de transmission global de la paroi considérée,
- SINT (en m²) est la surface intérieure de la paroi considérée,
- TSA (en °C) est la température sèche de l’air intérieur du local adjacent,
34
- TSB,I (en °C) est la température sèche de l’air intérieur du local considéré.
Ⅱ.2.11 Parois en contact avec le sol :

Les apports par transmission à travers les parois en contact avec le sol sont calculés en
fonction de la différence de niveau z. La différence de niveau z est : pour un plancher, la
différence comptée positivement entre le niveau de la face supérieure du plancher et le niveau
du sol, pour un mur enterré, la différence comptée positivement entre le niveau du sol et le
niveau considéré du mur.
Les apports par transmission pour une paroi en contact avec le sol, plancher ou mur
enterré, sont considérés comme nuls au-delà de z = zlim (z > zlim) dont les valeurs sont données
dans le tableau ci-dessous :
Tableau I.3 : valeurs de z lim (en m)
Valeur de z Lim
Température du local Zone climatique
conditionnée A B, C B’ C1 C2 C3
altitude<500 m
24 3,80 5,40 5,60 7,80 7,40 5,80
25 3,40 5,00 5,20 7,40 7,00 5,40
26 3,00 4,60 5 ,80 7,00 3,60 5,00
27 2,60 4,20 5,40 5,60 3,20 4,60
Formule de calcul
APO = Kete× SAPP (TSM – TsI,B) (II.18)

II.2.12 Apports à travers les parois vitrées :
Les apports à travers les parois vitrées intérieures AVT(t) sont dus uniquement à la
différence de température de part et d’autre de la paroi. Ils sont donnés par :
AVT(t) = KETE × SOUV × [(TSe(t) - Clnc) – TSB,I] (II.19)
- SOUV (m²) est la surface de l’ouverture dans la paroi opaque,
- TSe(t) (°C) est la température extérieure sèche à l’heure t,
- TSb,i (°C) est la température intérieure sèche de base pour le local considéré
- Clnc est un coefficient correcteur,
Les gains dus au rayonnement solaire à travers les parois vitrées extérieures AVE(t)
sont donnés par :
AVE(t) = [SVENS × It + (SV – SVENS) × Id] × FS × NPVI(t) (II.20)
35
- SV (m²) est la surface totale vitrée,
- SVENS (m²) est la surface vitrée ensoleillée,
- It (W/m²) est le rayonnement total maximal réel,
- Id (W/m²) est le rayonnement diffus maximal réel,
- NPVI(t) représente le coefficient d’amortissement relatif aux gains par

ensoleillement à travers les parois vitrées à l’heure t considérée,
- FS est le facteur solaire du vitrage,
Ⅱ.2.13 Détermination du coefficient d’amortissement

Le coefficient d’amortissement NPVI(t) est donné dans les tableaux ci-après en fonction
de la durée de fonctionnement des installations, de la présence ou non de protections (stores,
volets, rideaux, etc.), de la masse M rapportée à la surface du plancher, de l’orientation de la
paroi vitrée et de l’heure solaire vraie.
En l’absence d’informations, on pourra prendre :

- Locaux à usage d’habitation, d’hébergement et de vente :
- Zones A, B, B’ et C : 12 heures de fonctionnement,
- Zones D1, D2 et D3 : 16 heures de fonctionnement ;
- Locaux à usage d’enseignement, de bureaux, d’accueil et de réunion : 12 heures de
fonctionnement ;
- Locaux de restauration et à usage artisanal : 16 heures de fonctionnement.
Pour l’utilisation des tableaux I.4, toute paroi vitrée sans protection intérieure est
considérée comme une paroi vitrée sans protection (même si elle comporte une protection
extérieure, ou si elle est protégée de l’ensoleillement direct par des saillies).
La masse M rapportée à la surface du plancher est donnée par :
M= ∑mext+1/2∑msep+∑mi (II.21)
∑Spla
- mext (kg) est la masse des parois séparant la zone thermique considérée et
l’environnement extérieur (murs extérieurs, toiture, plancher éventuel en contact avec le sol) ;
- msep (kg) est la masse des parois séparant la zone thermique considérée et les
autres zones thermiques (cloisons, planchers intermédiaires, plafonds) ;
36
- mi (kg) désigne toute masse fixe (cloisons, murets, cage d’escalier, etc.) se
trouvant au sein de la zone thermique considérée.
- Spl (en m²) est la surface de plancher pour la zone thermique considérée ; si une
zone thermique a été définie sur plusieurs niveaux, Spl est la somme des surfaces de plancher.
Pour le calcul de M, on ne tiendra pas compte de la masse des parois vitrées.

Tableau I.4 Valeur de npv
Ⅱ.2.14 facteurs solaire

Le facteur solaire FS est égal au rapport du flux total transmis par le vitrage antisolaire
(FTas) au flux transmis à travers un vitrage ordinaire (FT vo) dans les conditions de référence.
En d’autres termes, FS = 1 pour un vitrage ordinaire, FS < 1 pour un vitrage antisolaire.
37
Pour réduire le coût de la climatisation, on est donc amené à installer des vitrages
antisolaires. Les valeurs de FS sont données dans le tableau ci-dessus :
Tableau II.5 : facteur solitaire des vitrages spéciaux en simple épaisseur
Apport Nature Coefficient de simultanéité Cs

des Locaux
apports
Bureaux Logement Vente Artisanal enseignements
Réunions hébergement Restauration industriel
Accueils
Occupant Sensible 0.80 0.50 0.80 0.90 1.00
Latent
Moteurs Sensible 0.60 0.50 0.85 0.85 0
Électriques
Appareils Sensible 0.60 0.50 0.85 0.85 0

électrique Latent
Eclairages Sensible 0.70 0.35 0.90 0.85 0.90

fluorescents non
encastrés
Eclairage Sensible 0.70 0.35 0.90 0.85 0.90
fluorescent
encastres,
éclairage
incandescent
non encastré
Eclairage Sensible 0.70 0.35 0.90 0.85 0.90
encastré dans un
faux-plafond
avec plénum de
reprise
Appareils à gaz Sensible 0 0.50 0.85 0.85 0
Latent
Réservoir Sensible 0 1.00 1.00 1.00 0
Latent
Evaporation libre Sensible 0 1.00 1.00 1.00 0
Latent
Vapeur vive Sensible 0 0 1.00 1.00 0

Latent
Tuyauterie, 1.00 1.00 1.00 1.00 1.00
conduits d’airs Sensible
38
Ⅱ.2.14 Le rayonnement total maximal réel It

Le rayonnement total maximal réel It et le rayonnement diffus maximal réel Id [W/m²]
sont donnés par :
It = [CCADRE× CLIMP× CALT× CROSE] It,b (II.22)
Id = [CCADRE× CLIMP× CALT× CROSE] Id,b (II.23)
Où :
- CCADRE est un coefficient tenant compte de la nature du cadre ;
- CLIMP est un coefficient tenant compte de la limpidité de l’atmosphère ;
- CALT est un coefficient tenant compte de l’altitude ;
- CROSE est un coefficient tenant compte de la valeur du point de rosée
- It,b et Id,b (en W/m²) sont les rayonnements total et diffus de base
 Le coefficient de correction CCADRE est donné ci-après :
- CCADRE= 1,17 pour un cadre métallique,
- CCADRE = 1 pour un cadre en bois, ou en PVC.
 Le coefficient de correction CLIMP est donné ci-après :
- CLIMP = 1 pour une atmosphère limpide ; on pourra adopter une atmosphère limpide
en zone rurale, hors agglomération (oasis, hameaux, etc.) ;
- CLIMP= 0,92 pour une atmosphère peu limpide ; on pourra adopter une atmosphère peu
limpide pour les zones suburbaines, les zones résidentielles des villes ;
- CLIMP = 0,87 pour une atmosphère obscure ; on pourra adopter une atmosphère
obscure pour le centre des grandes villes, les zones industrielles, à proximité d’industries
polluantes.
 Le coefficient de correction CALT est obtenu en augmentant sa valeur de 0,7 % par

300m en prenant CALT= 1 pour l’altitude 0 (par exemple, si l’altitude est égale à 349 m,
(CALT= 1,007)
 Le coefficient de correction Cros est donné dans le tableau II.6 en fonction de la zone
climatique.
39
Tableau II.6 valeur de Cros en fonction des zones thermiques
Zone Cros
A alt<500 0,99
500≤alt<1000 1,02
alt>1000 1,02
B alt<500 1,03
500≤alt<1000 1,05
alt>1000 1,07
B1
toutes 1,10
B2
Ⅱ.3 APPORTS DE CHALEUR INTERNES

On désigne par apports de chaleur internes, ou gains internes, les quantités de chaleur
dégagées sous forme latente ou sensible à l’intérieur des locaux conditionnés.
Les gains internes prévus par ce le DTR 2012 sont ceux ayant pour origine les
occupants, les machines entraînées par des moteurs électriques, les appareils électriques,
l’éclairage, les appareils à gaz, les réservoirs, l’évaporation libre, l’introduction de vapeur
vive, les tuyauteries et les conduits d’air intérieurs.
Ⅱ.3.1 Principes de calcul

Les apports internes, lorsqu’ils existent, sont calculés en les considérants constants
pendant la durée de fonctionnement de l’installation de climatisation.
Les apports internes sont calculés en considérant un amortissement des gains sensibles.
On utilise pour les besoins du calcul un coefficient dit d’amortissement noté NAI(t).
Le coefficient d’amortissement NAI(t) est à affecter séparément pour chaque apport

interne.
Les apports internes sont déterminés en multipliant, pour chaque apport, les gains
maximaux par un coefficient de simultanéité (noté CS) spécifique à chaque apport.
Ⅱ.3.2 Les apports de chaleur internes AL(t) a l’instant t sont données par la
formule suivante :
(II.24)
- AISJ (W) représente la partie sensible de l’apport interne j ;
40
Tableau II.7 : coefficients de simultanéité de gains internes
Apports Nature des Coefficient de simultanéité
apports Locaux
Bureaux logement restauration artisanal enseignement
Occupants Sensible 0,8 0,5 0,8 0,9 1
Eclairages Sensible 0,7 0,35 0,9 0,85 0,9
Ⅱ.3.3 Gains dus aux occupants

Les gains dus aux occupants sont une source de chaleur sensible et latente, pour les
calculs on doit tenir compte de nombre des occupants dans chaque pièce du bâtiment.
Si une zone thermique est définie, le nombre des occupants et de 4 occupants pour tout
le bâtiment pour chaque pièce la moyenne des occupants et d’un occupant.
Les gains internes pour un occupant sont donnés dans le tableau ci-dessus :
Tableau II.8 : gains dus aux occupants
Degré Exemple Température sèche du local (°C)
d’activité d’applicati 28 27 26 24 21
on Apports sensible (sen) et latent (lat)
en w/pers
Sen lat Sen lat Sen lat Sen lat Sen lat
Assis,travail Logement, 52 79 58 73 63 69 71 60 83 94
très léger
bureaux 52 94 58 88 64 83 74 72 85 62
Ⅱ.3.4 Gains dus aux machines électriques :

Les appareils peuvent constituer à la fois une source de chaleur sensible et latente, ou
seulement une source de gains sensibles.
Une machine électrique utilisée dans un local dégage intégralement l’équivalent

calorifique de la puissance électrique moyenne absorbée
Pour la détermination des apports dus aux machines électriques, on utilisera soit :
- Les indications données par le fabricant ;
- Les formules du tableau Ⅱ.9pour le calcul de la puissance absorbée, cette puissance

constituant les gains sensibles ; pour certains appareils (de cuisines, de restaurants, etc.), il y a
lieu de tenir compte aussi des gains latents ; on pourra alors utiliser les formules relatives aux
gains par évaporation
- À défaut, les valeurs des tableaux ci-après :
41
Tableau II.9 : Apports calorifiques dus aux appareils ménagers
gains à admettre (en w)

appareils Sens Lat total
Appareils ménagers
four électrique 3000 610 3610
four micro-onde 1400 - 1400
cuisinière électrique 3000 1500 4500
Télévision 150 - 150
réfrigérateur(2001) 175 - 175
machine à laver (3kg) 3000 1500 4500
micro-ordinateur 450 - 450
moulin à café 500 140 640
Mixer 1000 277 1277
fer à repasser 500 280 780
Aspirateur 200 - 200
Ⅱ.3.5 Gains dus à l‘éclairage :

AI = ∑(WnCmeCcr) (II.25)
- Wn ( W) est la puissance nominale de l’ampoule ou du tube fluorescent ;
- Cmeest un coefficient de majoration ; il est égal à :
• 1,2 pour les lampes à incandescence,
• 1,25 pour les tubes fluorescents ;
- Ccr est le pourcentage de chaleur résiduelle correspondant à la part d’énergie restante

dans la salle ; Ccr est :
• Égal à 1 pour les installations dont l’éclairage n’est pas raccordé à un système
d’extraction d’air,
• Est donné dans le tableau 10 quand l’éclairage est raccordé à un système d’extraction
d’air.
Tableau II.10: valeurs de Ccr
Débit d’air pour 100w de puissance de 20 30 50 100

lampe ou de tube (m3 /h)
Ccr 0.6 0.5 0.4 0.35
42
Si l’on ne connait pas la puissance installé pour l’éclairage, on calculera les gains dus à
L’éclairage en utilisant le tableau ci-dessus :
Tableau II.11 : gains dus à l’éclairage

Désignation du local ou Intensité Puissance dégagée des lampes
ou
genre d’activité lumineuse
des tubes (W/m2)
recommandé
(lux) lampes à lampes
incandescence fluorescentes
Entrepôt, pièces 120 25 8
d’habitation, théâtre
Bureaux, salle de cours 250 55 16
Salles de lecture, 500 110 32
laboratoires, magasins
Montages fins, ateliers, 750 170 50
supermarchés
Dessin industriel, grands 1000 - 65
bureaux,
Locaux industriels 1500 - 100
montages très fins
43
Introduction :
Le local étudié est un projet à un seul niveau, destiné à usage d’habitation. Le local
étudié estsans voisinage, c’est notre volume thermique étudié.
Présentation générale
Le travail consiste à évaluer, en premier lieu, le bilan thermique été du local, le calcul
desapports de référence et en fin procéder à la vérification réglementaire.
Conformément au DTR, la vérification réglementaire des locaux à usage d’habitation, de

bureaux et d’hébergement, doit s’effectuer selon les étapes suivantes :
Calcul pour l’ensemble du local (supposé conditionner) à 15h TSV, pour le mois de
Juillet :
-des apports par les parois opaques aériennes.
-des apports par les parois vitrées.
-calcul des apports calorifiques de référence.
-vérification de la conformité à la réglementation thermique d’été des locaux. A titre de

rappel, la somme des apports calorifiques par les parois vitrées et les parois opaques doit
vérifier au mois de Juillet à 15 h TSV, pour une température sèche intérieure de 24°C, la
relation ci-après :
APO (15 h) AV (15 h) 1,05AREF (15 h) [w]
Données relatives au site :
Le projet est un logement d’une superficie de 69.26 m², implanté à TIGZIRT sur mer.
Lesparamètres géographiques sont :
Altitude : 150 m.
Latitude : 36. °EST
Longitude : 04° Nord.
Orientation du local : Nord
Les paramètres thermiques de base sont :
Zone climatique « A » Alt < 500 m.
44
Plan de local étudié :
45
Les conditions de base extérieures :
- La température de base extérieure sèche : Tsbe =34°C.
- L’écart annuel de température EAT =31°C.
- L’écart diurne =9°C.
- L’humidité relative extérieure : Hre =60%.
- Humidité spécifique : g = 20,2 g de vapeur d’eau /kg d’air sec ;(D’après le

diagrammepsychrométrique).
- Les conditions de base intérieures :
- La température intérieure : Ti =24°C.
- L’humidité relative intérieure : Hri =45%Humidité spécifique : g = 11 g de

vapeurd’eau/ kg d’air sec (D’après le diagramme Psychrométrique).
 Le calcul consiste à déterminer les paramètres suivants :
 Apports à travers les parois opaques extérieures (APO(t)).
 Apports à travers les parois vitrées (AV(t)).
 Calcul des charges internes et externes du local.
Coefficient global de transmission « K »

1/K=∑R +1/hi +1/heen w/m2.°C
Calcul des résistances thermiques :
 mortier e = 0,02 m et 𝜆 = 1,4 W/m. °C R = e/ 𝜆 = 0,02 / 1,4 = 0,014m².°C/W
 Brique e = 0,1 m et 𝜆 = 0,65 W/m.°C R = e/ 𝜆 = 0,1 / 0,65 = 0,15m². °C/W
 Lame d’air e = 0,05 m 𝜆 = 0,26W/m.°C
R = e/ 𝜆 = 0,05 / 0,26 = 0,19 m². °C/W
 Brique e = 0,1 m et 𝜆 = 0,65 W/m.°C R = e/ 𝜆 = 0,1 / 0,65 = 0,15m². °C/W
 Mortier e = 0,01 m et 𝜆 = 1,4 W/m.°C R = e/ 𝜆 = 0,01 /1,4 = 0,007m². °C/W
 Plâtre e = 0.01m 𝜆 = 0,35 W/m.°C R = e/ 𝜆 = 0,01 / 0,35 = 0,03m². °C/W
Résistance superficielle :
1/he + 1/hi = 0.14 m². °C/W
46
1/K= [(0,014 + 0, 15 +0, 19 + 0, 15 +0,007 +0.03) +0, 14)]
1/K = 0.69 m². °C/W.
K=1.45W/m²°c.
Tableau II.12: coefficient de transmission

Désignation Compositio e (m) Λ Ρ (Kg /m3) R 1/hi+1/he 1/K K
n (W/m°C) (m²°C/W)
(m²°C/W) (m²°C/W) (W/m
²°C)
Mur Mortier 0.02 1.4 0.014 0.14 0.69 1.45
Brique 0.1 0.65 0.15
Lame d’air 0.05 0.26 0.19
Brique 0.1 0.65 0.15
Mortier 0.01 1.4 0.007
Plâtre 0.01 0.35 0.03
Plancher Béton armé 0.15 1.75 0.086 0.14 0.226 4.65
Haut
Plancher bas Béton armé 0.15 1.75 0.086 0.14 0.226 4.65
Porte 5.8
Fenêtre 3.8
II. Calcul des apports et déperditions thermiquesCalcul des apports

Chambre I : on commence par la chambre une comme indiqué dans le plan, elle est
constituée de trois façades exposées à l’extérieur.
Les apports thermiques Du séjour sont indiqués dans le tableau récapitulatif suivant :
Tableau II.2 : tableau récapitulatif des apports thermiques
Pièce mur APO AV Ais AIL Ainfs AifL As Ats Aref Aref•(APO+AV)
séjour sud 97.4 156.79 426.17 90 77,95 67.95 757.32 1035.18 253.21 254.19
nord 57.48 0 520.41 90 208.32 181.6 692.17 967 82.85 57.48
ouest 210 694.75 426.17 90 396.7 345.8 1727.62 2054 1548.11 904.75
total 364.88 851.54 1372.75 270 682.97 595.35 3177.11 4047.18 1884.17 1216.42
Mur Sud
Apport a travert les parois opaques :
APO(t)=1.2K étéSint DTe(t)

DTe(t) = Δtes(t) + cΔte+ 𝛼[ΔTem(t) − Δtes(t)] 𝐼𝑡 ; 𝑏
9 𝐼𝑡 ; 𝑏 (40)
DTe(t) = 13.9 −0.1 + 0.6[−13.9)] ×1.13
47
DTe(t) =4.37C Sint=12.82 m2. Kété = 1.45 W/°m2 Eb =9 °C CΔte = − 0,1°C

Δtem (15h) =0°C
Δtes (15h)=13.9°C. It,b=217W/m2 I t,b (40) =217W/m2
D’où
APO (15h) = 1.2 * 1.45* 12.82* 4.37
APO (15h) = 97.4 W/m²

Apports à travers les parois vitrées (AV(t))
AV (t) = AVT (t) + AVE (t)
AV (t) = 76.36+80.43=156.79 W/m²

Apport due au gradient de température (par transmission) (AVT(t))
AVT (t) = 1,2 * Kété * Souv * [Tse (t) - Tsb,i]
Pour : Kété = 3.8W/m2 ; Souv = 1.68m2; Tse=34°C; Tsbi = 24°C.
AVT(15h) = 1.2 * 3.8 * 1.68*[34 –24]

AVT(15h) = 76.60 W/m²
Apport par ensoleillement (AVE(t))
AVE (t) = [Svens * It + (SV - SVens) * Id] * Fs * Npvi (t)
FS= 0. 45
Npvi (15h) = 0.68 Id,b =48,5 W/m2 Id =55.11

It,b= 246.6 W/m2, d’ou: It = 225.4 W/m2 Avec:
Ccadre =1 Climp =0,92 Calt =1,0035 Cros =0.99
AVE (15h) = [1 * 225.4 + (1.68-1) * 55.11] * 0.45* 0.68
AVE (15h) =80.43 W

Calcul des apports de chaleur par introduction d’air extérieur
ARENS (t) = 0,320 * Qvan * (Tse(t) - Tsb,i)
ARENL (t) = 0,797 *Qvan*Max [(Hsb,e- Hsb,i);0]
Pour un nombre d’occupants : N= 3 personnes,
Qvan min = 25 m3 /h
48
ARENS (15h) = 0,320 * (3*25) * (34 - 24)
ARENS (15h) = 240w
ARENL (t) = 0,797 *75*Max [(14.5- 11);0]
ARENL (t) = 209.21 w
Tableau II.13: apports de chaleurs par introduction d’air extérieur
LOCAL Nbr S en (m2) qvan (m3 /h) Tse-Tsbi ARENS(w) Hbe-Hbi ARENL(w)
occupants
Logement 3 1.68 75 10 240 3.5 209.21
Apports dus aux infiltrations d’air
AINFS (t) = 0,320 * QVinf * ( Tse(t) - Tsb,i )
AINFL (t) = 0,797 * QVinf* Max [(Hsb,e-Hsb,i) ; 0]
Avec: Qv inf =Σ(Qvo inf * Souv) (m3/ h)

Qvo inf = 14.5 m3/ h.
Souv = 1.68 m2.
Qv inf = 24.36 m3/ h
AINFS (t) = 0,320 * 24.36* (34 - 24)
AINFS (t) =77.95 w

AINFL (15h) = 0,797 * 24.36* (14,5 - 11)
AINFL (t) = 67.95 W
Tableau II.14: Apports des infiltrations d'air

LOCAL Nbr S en (m2) qvan (m3 Tse-Tsbi AINFs (w) Hbe-Hbi AINFL (w)
occupants /h)
Logement 3 1.68 75 10 77.95 3.5 67.95
Apports internes
AI(t) = Σj (Csj * AIsj * NAIj) + Σj (Csj * AILj)
Tableau II.15: apports de chaleurs interne

Chambre I Asoc Asm Ase AIs AIL
Mur sud 101.17 285 40 426.17 90
49
Apport de chaleur par les occupants
Asoc = Cso * AIso * NAIoALoc = Cso * AILO

Avec: Ais=71(W/personne)
AIL= 60 (W/personne) Cso=0,5 NAIo=0,95
Asoc = 3*0,5*71*0,95
Asoc = 101.17 W
Aloc = 3*60*0,5Aloc = 90W.
Tableau II.16: Apports par les occupants

LOCAL Nbr Aiso Cso NAIo Asoc (w) AILOc Aloc (w) AOc totale
Occupants (W/personne) totale(w)
Logement 3 71 0.5 0.95 101.17 60 90 191.17
Apports dus aux machines et appareillages Asm=Csm * Alsm * NAIm Alm = Csm *
ALlmMachine :
Ordinateur (Aism = 450W/ordinateur et Ailm = 0W/ordinateur).
Télévision (Aism = 150W/télévision et Ailm = 0W/ télévision)
AlL = 0 w
AlS = 285 w
Tableau II.17: Apports dus aux machines et appareillages
LOCAL Type Nbr As/app Cs Nai Asm (w) Alm(w) Aism T (w)
Logement Television 1 150 0.5 0.95 71.25 0 285
Micro-ordinateur 1 450 0.5 0.95 213.8 0
Apports dus aux éclairages
Pour 4 Lampes fluorescentes : Wn = 8 W/m²,

AIse = Σ ( Wn * Cme * Ccr)
Lampes fluorescentes : Wn = 8 W/m², Ccr = 1, Cme = 1,25,
AIse = 4*(1,25 * 1* 8)
AIse = 40 W/m²
Tableau II.18: Apports dus aux Eclairages

Local Wn (w) Cme Cr Aise (w)
Logement 8 1.25 1 40
50
Les apports sensibles

As = APO + AV + AIs + AINFs
As =97.4+156.7+426.17+77.95
As = 757.32 w /m²
Tableau II.19: Apports sensibles
Local APO AV AIs AINFs As

Logement 97.4 156.8 426.17 77.95 757.32
Les apports latents
Al = 90 w
Les apports effectifs

AEs = (CDas× As) + (BF × ARENs)
AEs = (1.05×757.32) + (0.4 × 240)

AEs= 891.18 w
Tableau II.20: Apports effectifs sensibles
Local CPas As BF ARENs AEs

Logement 1.05 757.32 0.4 240 891.18
AEl = (CDal× Al) + (BF × ARENl)

AEL = (1.05× 90) + (0.4 ×209.21)
AEl= 178.18 w
Tableau II.21: Apports effectifs latents
Local CPas Al BF ARENl AEl
Logement 1.05 90 0.4 278.95 178.18
Les apports totaux
ATs = (CDas× As) + ARENs
ATs =1035.18w
Vérification réglementaire :
APO (15 h) + AV (15 h) ≤ 1,05 × Aréf (15 h) (W)
Aréf = Aréf,PH+ Aréf,PV + Aréf,PVI (W)
51
Calcul des apports calorifiques de référence
Aréfréf,PV= Σ (1.4 × 12.82×6.88)

Aréfréf,PV= 123.48w
Apports de référence des parois vitrées
Aréf, PVI = AVEréf +AVTréf
AVE = [SVens* It + (SV- SVens) Id] *FS* NPVI,réf
AVEréf (15h) = [1 * 225.4+ (0.68) * 55.11]* 0.45* 0.68
AVEréf (15h) = 80.44w

AVTréf = e × Souv× ΔTSréf, PVI
AVTréf (15h)= 4.2*1.68*7
AVTréf (15h) =49.39w

D’où :
Aréf = Aréf, PH + Aréf, PV + Aréf, PVI
Aréf =(123.38+129.83)
Aréf = 253.21
1,05 × Aréf (15 h) =265.87
APO (15h) +AV (15h) = 254.2w/m²
Aréf (15 h) × 1,05 ≥ APO (15 h) + AV (15 h)
Donc l’exigence réglementaire est vérifiée
Mur nord:
Apport à travers les parois opaques

𝛼
DTe(t) = Δtes(t) + cΔte+ [ΔTem(t) − Δtes(t)] 𝐼𝑡;𝑏
9 𝐼𝑡;𝑏(40)
DTe(t) = 0 −0.1 + 0.6[4.4] ×1.03

DTe(t) =2.62°C
52
Sint =5.99m2.Kété =1.45 W/°m2 Eb = 9°CCΔte = − 0,1°C
Δtem (15h) =4.4 °C Δtes (15h)=0°C. It,b=48.5W/m2It,b (40) =47W/m2
D’ou
APO (15h) = 1,2 * 1.45 *12.61 * 2.62 APO (15h) = 57.48 W

AV (t) = AVT (t) + AVE (t)
AV (t) = 0 w
Apport à travers la porte
APO (15h) = 1,2 * 5.8 *1.89 * 2.62APO (15h) =34.46 W

AVT(15h) = 0W
Apport par ensoleillement (AVE (t))
AVE (t) = 0 w
Calcul des apports de chaleur par introduction d’air extérieur
ARENS (t) = 0,320 * Qvan * ( Tse(t) - Tsb,i )
Qvan min = 25 m3 /h
ARENS (15h) = 0,320 * (3*25) * (34 - 24)
ARENS (15h) = 240w
ARENL (t) = 0,797 *3*25*Max [(14.5- 11);0]

ARENL (t) = 209.21 w
Tableau II.22: Apports par introduction d’air extérieur
occupants
Logement 3 14.5 75 10 240 3.5 209.21

Souv = 1.89m2.
53
Qv inf = 21.7m3/ h
AINFS (t) =208.32 w
AINFL (t) = 181.6 W
Tableau II.23 : Apports par infiltration
occupants /h)
Logement 3 1.89 65.1 10 208.32 3.5 181.6
Apports internes
Tableau II.24 : Apports interne
Sejour Asoc Asm Ase AIs AIL
Mur nord 111.36 285 124.05 520.41 90

Asoc = Cso * AIso * NAIo ALoc = Cso * AILO
Avec: Ais=71(W/personne) AIL=60(W/personne)
Cso=0,5 NAIo=0,95
Asoc = 3*0.5*71*0.95Asoc = 101.17 w
Aloc = 3*60*0.5Aloc = 90 w
Tableau II.25 : Apports par occupants

occupants (W/personne) totale (w)
Logement 3 71 0.5 0.95 101.17 60 90 191.17
Apports dus aux machines et appareillages
Asm=Csm * Alsm * NAIm
Alm = Csm * ALlmMachines :

Télévision (Aism = 150W/télévision et Ailm = 0W/ télévision).
AlL = 0w
AlS = 285 w
54
Logement Télévision 1 150 0.5 0.95 71.25 0 285
Micro-ordinateur 1 450 0.5 0.95 213.8 0

AIse = Σ ( Wn * Cme * Ccr)
Lampes fluorescentes : Wn = 8 W/m², Ccr = 1, Cme = 1,25,
AIse = 4*(1,25 * 1* 8)AIse = 40 W



As = APO + AV + AIs + AINFsAs = 692.17w
Tableau II.28: Apports sensibles

Logement 57.48 0 426.17 208.32 692.17
Les apports latents

Al = 90 w

AEs = (CDas× As) + (BF × ARENs)AEs = (1.05× 692.17) + (0.4 ×240)
AEs= 823 w
Tableau III.28: Apports effectifs sensibles
AEl = (CDal × Al) + (BF × ARENl)
AEL = (1.05× 90) + (0.4 ×209.21)

AEl= 178.2 w
55
Chapitre III : Simulation numérique de la ventilation naturelle
Tableau II.29: Apports effectifs latents

Logement 1.05 90 0.4 209.21 178.2
Les apports totaux

ATs = (CDas× As) + ARENsATs =967w
APO (15 h) + AV (15 h) ≤ 1,05 × Aréf (15 h) (W)
Calcul des apports calorifiques de référence
Aréfréf,PV= Σ (c × Sint ×ΔTSréf,PV)Aréfréf,PV= Σ (1.4 × 12.61×4.66)

Apports de référence des parois vitrées
Aréf, PVI = AVEréf +AVTréf wAVE = 0
AVT = e × Souv× ΔTSréf, PVI AVT(15h)=0 w

D’oùAréf = Aréf, PH + Aréf, PV + Aréf, PVIAréf = 82.85 w
1,05 × Aréf (15 h) =86.37 APO (15h) +AV (15h) = 57.48
D’ou 1,05 × Aréf (15 h) ≥ APO (15 h) + AV (15 h)
Donc l’exigence réglementaire est vérifiée.Mur ouest:
Apport a travers les parois opaques

𝛼
DTe(t) = Δtes(t) + cΔte+ [ΔTem(t) − Δtes(t)] 𝐼𝑡;𝑏
9 𝐼𝑡;𝑏(40)
DTe(t) = 10.6 −0.1 + 0.6[−10.6]×1DTe(t) =4.14°C
Sint =29.61 m2. Kété = 1.45 W/°m2°C Eb = 9 °C CΔte = − 0,1

Δtem (15h) = 0 °C Δtes (15h) = 10.6°C.
It,b= 516 W/m2 It,b (40)=516 W/m2
56
D’où
APO (15h) = 1,2 * 1.45 * 29.16* 4.14APO (15h) = 210 w

AV (t) = AVT (t) + AVE (t) AV (t) = 1002.55 w
AVT (t) = 1,2 * Kété * Souv * [Tse (t) - Tsb,i] Pour:
Kété = 6.8W/m2 ; Souv = 8.55m2; Tse=34°C ; sbi = 24°C.
AVT(15h) = 1.2 * 6.8 * 8.55*[34 –24]
AVT(15h) =697.68 W
Apport par ensoleillement (AVE(t))

AVE (t) = [Svens * It + (SV - SVens) * Id] * Fs * Npvi (t)
Fs= 0. 45 Npvi (15h) = 0.68 Id,b =48,5 W/m2Id =50.37
It,b= 155.5 W/m2, d’ou: It = 161.5 W/m2 Avec:
Ccadre =1 Climp =0,92 Calt =1,0035 Cros =0.99
AVE (15h) = [5.09 * 161.5 + (8.55-5.09) * 50.37] * 0.45* 0.6
AVE (15h) =304.87W
Calcul des apports de chaleur par introduction d’airextérieur
Pour un nombre d’occupants : N= 3 personnes,
Qvan min = 25 m3 /h
ARENS (15h) = 0,320 * (3*25) * (34 - 24)
ARENS (15h) = 240w
ARENL (t) = 0,797 *75*Max [(14.5- 11);0]
ARENL (t) = 209.21 w
Tableau II.30: apports de chaleurs par introduction d’air extérieur
occupants
Logement 3 1.68 75 10 240 3.5 209.21
57
Avec: Qv inf =Σ(Qvo inf * Souv) (m3/ h)
Qvo inf = 14.5 m3/ h. Souv = 8.55 m2
Qv inf = 123.97m3/ h
AINFS (t) = 0,320 * 123.97* (34 - 24)
AINFS (t) =396.7 w
AINFL (15h) = 0,797 * 123.97* (14,5 - 11)
AINFL (t) = 345.8W
Tableau II.31: Apports des infiltrations d'air
occupants /h)
Logement 3 8.55 14.5 10 396.7 3.5 345.8
Apports internes
Tableau II.32: apports de chaleurs interne
Chambre I Asoc Asm Ase AIs AIL
Mur sud 101.17 285 40 426.17 90
Asoc = Cso * AIso * NAIo
ALoc = Cso * AILO
Avec: Ais=71(W/personne)
AIL= 60 (W/personne)
Cso=0,5
NAIo=0,95
Asoc = 3*0,5*71*0,95
58
Asoc = 101.17 W
Aloc = 3*60*0,5
Aloc = 90W.
Tableau II.33: Apports par les occupants

occupants (W/personne) totale(w)
Logement 3 71 0.5 0.95 101.17 60 90 191.17
Apports dus aux machines et appareillages
Asm=Csm * Alsm * NAIm
Alm = Csm * ALlmMachine :

Télévision (Aism = 150W/télévision et Ailm = 0W/ télévision).
AlL = 0 w
AlS = 285 w

Ise = Σ ( Wn * Cme * Ccr)Lampes fluorescentes :
Wn = 8 W/m², Ccr = 1, Cme = 1,25,AIse = 4*(1,25 * 1* 8)AIse = 40 W.

As = APO + AV + AIs + AINFs
As =210+1002.57+426.17+396.7=2035.42 w
Tableau II.36 : Apports sensibles
Logement 210 1002.55 426.17 396.7 2035.42
Les apports latents

Al = 90 w
AEs = (CDas× As) + (BF × ARENs)
59
AEs = (1.05×2035.42) + (0.4 ×204)
AEs= 2218.8 w.
Tableau II.37 : apports effectifs sensibles
Local CPas As BF ARENs AEs
Logement 1.05 2035.42 0.4 204 2218.8
AEl = (CDal× Al) + (BF × ARENl)
AEL = (1.05× 90) + (0.4 ×209.21)AEl=178.18 w
Tableau II.38 : apports effectifs sensibles

Logement 1.05 90 0.4 209.21 178.18
II.1.3 Les apports totaux

ATs = (CDas× As) + ARENsATs = 2341.2 w
Vérification réglementaire
APO (15 h) + AV (15 h) ≤ 1,05 × Aréf (15 h) (W)
Aréfréf,PV= Σ (c × Sint × ΔTSréf,PV)
Aréfréf,PV= Σ (1.4 × 29.61× 9.86)
Aréf, PVI = AVEréf +AVTréf [W]
AVE = [SVens* It + (SV- SVens) Id] *FS* NPVI,réf
AVE (15h) = [5.09 * 535.92+ (3.46) * 50.37]* 0.45* 0.68
AVE (15h) = 888.04 w
AVT = e × Souv× ΔTSréf, PVI
AVT(15h)= 4.2*8.55*7
AVT(15h) =251.37w
D’où: Aréf = Aréf, PH + Aréf, PV + Aréf, PVI
Aréf =408.7+251.37+888.04 =1548.11

1,05 × Aréf (15 h) = 1625.51w
APO (15h) +AV (15h) = 1212.55 w
D’où APO (15 h) + AV (15 h) ≤ 1,05 × Aréf (15 h)
Donc l’exigence réglementaire est vérifiée.
60
Conclusion
D’après le bilan thermique effectués, le flux de chaleurs joue un rôle très important sur
la ventilation pendant la période estivale.
Le choix du type du matériau de construction c’est un révélé fondamental dans la

détermination du besoin énergétique total à cause de la différence dans la conductivité
thermique entre les matériaux ; chaque matériau à une performance énergétique différente.
61
Ⅲ.1. Introduction :
Dans ce chapitre ont étudié la discrétisation des équations aux dérivées partielles a pour
le but de transformer le problème continu en un problème discret, de façon à pouvoir assurer
sa résolution par des méthodes numériques telle que la méthode des éléments finis, la
méthode des différences finis et la méthode des volumes finis. Cette dernière est la plus
utilisée actuellement dans les codes de simulation numérique CFD (Computationnel Fluid
Dynamics).
L’étude de notre problème, implique la résolution d’un système d’équations

mathématiques, aux dérivées partielles. Cela nécessite une approche par une simulation
numérique.
Nous avons opté dans cette étude pour l’utilisation du code de calcul FLUENT version
6.2.16 et comme solveur, et le GAMBIT version6.2.16 pour la description géométrique.
Ⅲ.2Calcul numérique
Le modèle mathématique constitué par un système d’équations aux dérivées partielles

(EDP) et transformé à l’aide d’une méthode de discrétisation dans un système d’équations
algébriques.
Le principe de calcul numérique est illustré par la figure (Ⅲ.1). [21]
Phénomène physique
Modèle mathématique
Système d’EDP non linéaires sur un domaine

continu
Modèle discrétisation + algorithme
Système des équations algébriques
Linéaire sur un domaine décret
Méthode de résolution
Solution sur un domaine décret
Figure (Ⅲ.1) : principe de calcul numérique.
62
Ⅲ.3 Les différentes méthodes numériques de discrétisation
En vue du passage d’un problème exact (continu) au problème approché (discret), on

dispose de plusieurs techniques concurrentes et complémentaires, les plus connues sont :
 Méthode des différences finis(MDF).

 Méthode des éléments finis(MEF).
 Méthode de volumes finis(MVF).
Ⅲ.3.a Méthode de différences finies

La méthode des différences finies est l’origine des premières simulations dans la
mécanique des fluides et transfert de chaleur (1960 à 1965). Son principe à discrétiser les
équations aux Nœuds de la grille en calculant chaque dérivée partielle à l’aide du
développement de Taylor, pour obtenir les équations linéaires qui lient les variables inconnues
en un nœud aux valeurs de ces mêmes variables aux nœuds voisins. La précision du schéma
dépend du nombre de nœud de l’ensemble retenu pour approximer une dérivée donnée et du
maillage effectué (en général utilisation de maillage régulier et structuré). [22]
Ⅲ.3.b la Méthode des éléments finis

La technique des éléments finis discrétise l’espace à l’aide d’éléments géométriques
simples (triangle ou quadrangles en général). Comme elle permet de modéliser des géométries
très complexes, elle est parfois préférée à la méthode des volumes finis. Ensuite, la forme
forte des équations est remplacée par la forme faible dans laquelle les inconnues sont
approximées par une combinaison linéaire de fonction de base dont le support est un élément.
La méthode des éléments finis a vu les premières utilisations pour résoudre les équations de
Navier Stockes dans les années soixante-dix. Dans la méthode des éléments finis les variables
inconnues sont calculées par une combinaison linéaire de fonctions.
L’approche éléments finis se décompose en deux cas :
 Approximation nodale des fonctions inconnues par sous-domaines (fonctions

d’interpolation).
 Utilisation de la méthode des résidus pondérés qui permet de construire une
formulation intégrale à partir des équations aux dérivées partielles (la méthode Galerkin est la
plus fréquemment utilisée). C’est une méthode conservative sur un domaine.
63
Ⅲ.3.c Méthode des volumes finis
Cette méthode consiste à subdiviser le domaine physique de l’écoulement en éléments

de volume plus ou moins régulier dans lesquels les équations, décrivant l’écoulement, écrites
sous une forme conservative, sont intégrées.
Le but de cette méthode est de convertir l’équation différentielle générale en un système

d’équation algébrique en mettant en relation les valeurs de la variable considérée faux points
adjacents d’un volume de contrôle typique. Cela est obtenu par l’intégration de l’équation
différentielle gouvernante dans ce volume de contrôle avec l’approximation des différents
termes (convection, diffusion et source).
Ⅲ.4 La méthode des volumes finis, elle consiste en:
 La division du domaine de calcul en des volumes de contrôle discret par un

maillage de calcul.
 L’intégration des équations gouvernantes sur les volumes de contrôle
individuels, afin de construire les équations algébriques pour les variables dépendantes
discrètes(les inconnus), tels que les vitesses, pressions, températures et les autres scalaires
conservés.
 La linéarisation des équations discrétisées et la solution du système
d’équations linéaires résultants pour pouvoir adapter les valeurs des variables dépendantes.
Avantages et inconvénients :
 De nombreux codes de simulation numérique en mécanique des fluides
reposent sur cette méthode : Fluent, Star CD, CFX, Fine turbo. Cette méthode est bien
adaptée aux lois de conservation.
 Permet des maillages non-structurés (géométries complexes).
 Difficulté à monter en ordre.
Ⅲ.5 Principe de la méthode des volumes finis

Les étapes des résolutions de la méthode des volumes finis sont les suivantes :
Ⅲ.5.a. Maillage
C’est la subdivision du domaine d’étude en grilles longitudinales et transversales dont
l’intersection représente un nœud, où on trouve les composantes u et v de la vitesse qui se
trouve au milieu de segments relient deux nœuds adjacents. La discrétisation du domaine est
obtenue par un maillage constitué d’un réseau de points.
Les grandeurs scalaires pression, température, sont stockées dans le nœud P du
64
maillage, tandis que les grandeurs vectorielles u et v sont stockés aux milieux des segments
reliant les nœuds. L’équation du transport est intégrée sur le volume de contrôle associe aux
variables scalaires est les équations de quantités de mouvement sont intégrée sur le volume de
contrôle associe aux composants de vitesse.
Le volume de contrôle de la composante longitudinale u est décalé suivant la direction x

par rapport au volume de contrôle principale, celui de la composante transversale v est décalé
suivant la direction de y.
Ce type de maillage est dit maillage décalé ; il permet une bonne approximation des flux
convectifs est une meilleure évaluation de des gradients de pression ainsi une stabilisation
numérique de la solution. [23]
Ⅲ.5.b. Discrétisation des équations gouvernantes
L’étape principale de la méthode des volumes finis est l’intégration des équations
régissantes pour cheque volume de contrôle. Les équations algébriques déduites de cette
intégration rendent la résolution des équations de transport plus simple.
Chaque nœud est entouré par un ensemble de surfaces qui comporte un élément de
volume.
Toutes les variables du problème et les propriétés du fluide sont stockées aux nœuds de
cet élément.
Les équations régissant l’écoulement, sont présentées sous les formes moyennées dans
un repère cartésien (x, y, z):
Ⅲ.5.c. Les équations gouvernantes
Le problème considéré est gouverné par les équations couplées de Navier Stokes pour
un fluide incompressible, tenant compte de l’hypothèse de Boussinesq, et de l’énergie dans
laquelle le terme de la dissipation visqueuse a été négligé. La projection des équations (Ⅲ.1),
(Ⅲ.2) et (Ⅲ.3) dans un plan de coordonnés cartésienne (x, y) donne respectivement :
∂P ∂
+ ∂x (Puj ) = 0 (Ⅲ.1)
∂t i
Notons que toutes équations régissant l’écoulement interne en turbomachine peuvent

être écrites sous la forme générale suivante (équation de transport de la variable φ):
65
∂ y ∂ ∂ϕ
(ρ∅ ) + ⋅ (ρuϕj ) = (Γeff (∂x )) + S∅ (Ⅲ.2)
∂t ∂xj ∂JCj j
𝑈j : est la composante de la vitesse suivant la direction xj.
𝜌 : La masse volumique.
Г : est le coefficient de diffusion.
S : est le terme source.
ɸ : Une des variables suivantes ϕ = 1, 𝑢, 𝑣,w, k,𝗌,e𝑡𝑇.
Ⅲ..5. c.1. Equation de continuité
6𝑢 6𝑣
+ =0 (Ⅲ. 3)
6𝑥 6𝑦
u et v étant les composantes du champ de vitesse V (u,v.) dans la direction x et y
respectivement.
Ⅲ.5.c .2 Equation de quantité de mouvement
Suivant (ox):
𝜕𝑣 𝜕𝑣 1 𝜕𝑃 𝜕2 𝑈 𝜕2 𝑣
𝑢⋅ +𝑣⋅ =− ⋅ +( 2
+ ) (Ⅲ. 4)
𝜕𝑥 𝜕𝑦 𝑃 𝜕𝑥 𝜕𝑥 𝜕𝑦 2
Suivant (oy):
𝑢𝜕𝑣 𝑣𝜕𝑣 1 𝜕𝑃 𝜕2 𝑣 𝜕2 𝑣
𝜕𝑥
+
𝜕𝑦
=−
𝜌 𝜕𝑦
+𝑣(
𝜕𝑥 2
+
𝜕𝑦 2
) − 𝑔𝛽(𝑇 − 𝑇0 ) (Ⅲ.5)
III.5.3.c Equation de l’énergie

𝜕𝑇 𝜕𝑇 𝜕2 𝑇 𝜕𝑇
𝑢⋅ +𝑣 = 𝑎( + ) (Ⅲ. 6)
𝜕𝑥 𝜕𝑦 𝜕𝑥 2 𝜕𝑥 2
L’équation différentielle de la continuité, de la quantité de mouvement et de l’énergie

formant le modèle mathématique de l’écoulement de la convection naturelle laminaire.
Où :
𝑢 : est la composante de la vitesse suivant la direction x.
P : est la composante de la vitesse suivant la direction y.
𝜌 : La masse volumique.
𝑇 : Température du local.
𝑇0 : Température d’entrée du local.
𝑎 : Diffusivité thermique.
66
Ⅲ.6 Approximation Boussinesq
Pour des configurations d’écoulement dans une enceinte, la résultante des forces
externes est limitée aux forces gravitationnelles de sorte que l’influence de la convection
naturelle provoquée par des gradients de la température soit évidente. Quand les forces de
flottabilité se produisent dans l’écoulement, la conception des propriétés physiques fixes
n’assortit plus le comportement de l’écoulement du fluide. L’effet de la convection naturelle
est pris en considération dans l’équation de quantité de mouvement par la variation de la
masse volumique.
Si en prenant un état thermodynamique de référence (masse volumique 𝜌0 et

température 𝑇0
L’équation d'état pour la masse volumique 𝜌 devient en utilisant le développement de

Taylor [16] :
𝑇−𝑇𝑜
𝜌 = 𝜌0 (1 − + ⋯) (Ⅲ.7)
𝑇0
En se limitant à un développement du premier ordre :
𝜌 = 𝜌0(1 – 𝛽 (𝑇 − 𝑇0)) (Ⅲ.8)
Le terme β de l’équation (III.8) est le coefficient de dilatation isobare du fluide,

Équations de conservations.
1 𝜕𝑃
𝛽=− (Ⅲ.9)
𝜌 (𝜕𝑇)𝑃=𝑐𝑠𝑡
Convection naturelle, créant une force motrice volumétrique (Archimède) à laquelle

s’oppose notamment une force de frottement visqueux
L’approximation de Boussinesq consiste à négliger les variations de la masse

volumique.
Ⅲ.7 Le modèle de turbulence

Il existe différents niveaux de modèles, se distinguent par leur degré de complexité,
c'est-à-dire par le nombre d’équations de transport supplémentaires introduites pour les
quantités turbulentes pour fermer le problème : on cite ainsi le modèle de turbulence à zéro,
une ou deux équations. Les modèles de turbulence les plus utilisés à l’heure actuelle dans les
codes CFD sont les modèles à deux équations et plus particulièrement le modèle de
turbulence k-ε standard qui est le plus répandu. [24]
67
Trois modèles de turbulence sont disponibles dans le code Fluent :
 Le modèle k-ε standard.
 Le modèle RSM.
 Le modèle RNG k-ε.
Ⅲ.8 Le modèle k-ε
Dans cette étude on a utilisé le modèle k-ε standard. Il est un modèle semi empirique
basé sur les équations de transport de l’énergie cinétique turbulente k et de sa dissipation ε
Dans le modèle k-ε la viscosité turbulente est exprimée en fonction de k et ε par la formule
suivante : [25]
K²
µ t = CµP (Ⅲ.10)
ε
Où : ε appelé taux de dissipation de l’énergie cinétique turbulente.
𝐶𝜇 : est une constante, dont sa valeur est égale à 0.09.
Les équations gouvernantes dans le régime turbulent sont toujours les équations de
Navier Stocks, l’équation de continuité et l’équation d’énergie, en ajoutant deux équations de
transport dites équations de fermeture du modèle K-ε standard. [26]
Ⅲ.9 Les nombres adimensionnels
Les résultats expérimentaux sont souvent exprimés en termes de nombres

adimensionnels, afin de permettre une application plus générale que celles pour lesquelles les
résultats sont obtenus. Les Nombres adimensionnels permettent donc une validation des
résultats par rapport aux travaux antérieurs.Quelques nombres adimensionnels sont présentés
dans ce qui suit:
Ⅲ.9.a Le nombre de Reynolds Re
Le nombre de Reynolds Re est le rapport des forces d’inertie aux forces de viscosité,
dans un écoulement considéré. Il caractérise la nature de l’écoulement (laminaire,
turbulent…).Il est défini comme suit:
Re=𝜌𝑉𝐿⁄𝜇 (Ⅲ.11)
Pour aller d’un régime laminaire vers un régime turbulent, cela nécessite le passage
68
d’une zone de transition caractérisée par son nombre de Reynolds critique qui a pour valeur :
𝑅e𝑐=2300, pour le cas d’une conduite.
Ⅲ.9.b Le Nombre de Prandtl Pr
Le Nombre de Prandtl Pr est le rapport de la viscosité cinématique 𝜇/𝜌 à la diffusivité

thermique 𝑎 = 𝗁⁄𝜌𝐶𝑃. Il caractérise les propriétés thermiques du fluide, en mettant en rapport
la distribution de vitesses (gouvernée par la viscosité cinématique) et la distribution de
températures (gouvernée par la diffusivité thermique), on le définie par :
μCp
𝑃𝑟 = (Ⅲ. 12)
λ
Cp: Capacité calorifique massique (J/𝑘𝑔.𝐾)
λ : coefficient d’échange conductive (W/m2.K)
μ : viscosité cinématique (Pa.s)
Ⅲ.9.c Le Nombre de Nusselt Nu
Le Nombre de Nusselt Nu est le rapport du flux thermique 𝜑, transmis par convection à

un flux thermique transmis par conduction à travers le fluide, il est donné par la relation :
ℎ𝐷
𝑁𝑈 = (Ⅲ. 13)
λ
Ⅲ.9.d Le nombre de Rayleigh Ra

Le nombre de Rayleigh Ra est le paramètre de contrôle de la convection thermique. Plus
le nombre de Rayleigh est grand, plus la convection est intense.
Le nombre de Rayleigh peut s’interpréter comme une mesure de l’importance du

mécanisme responsable de l’instabilité du fluide (la poussée d’Archimède) par rapport aux
mécanismes de freinage (la diffusivité thermique et la viscosité). Il est également possible de
voir le nombre de Rayleigh comme la différence de température. Le nombre de Rayleigh est
donné par la relation :
(Ⅲ. 14)
Avec :
 Gr - nombre de Grashof
 Pr - nombre de Prandtl
 g - accélération de la pesanteur (m/s2)
69
 Lc - longueur caractéristique (m)
 Ts - température de la paroi (K)
 T∞ - température du fluide loin de la paroi (K)
 ν - viscosité cinématique (m2/s)
 α - Diffusivité thermique (m2/s)
 β - coefficient de dilatation thermique volumétrique (K−1)
Pour Ra < 109 : écoulement est laminaire, Ra > 109 : écoulement est Turbulent.
Ⅲ.9.e Nombre de Grashof
Le nombre de Grashof compare les forces d’Archimède (ou forces de flottabilité) aux
forces visqueuses qui s’exercent sur un élément de fluide dont la température diffère de celle
du milieu environnant.
(Ⅲ. 15)
avec :
 g - accélération de la pesanteur
 β - coefficient de dilatation
 ΔT - différence de température
 Lc - longueur caractéristique
 ρ - Masse volumique du fluide
 μ - viscosité dynamique
Ⅲ.10 Logiciel de simulation gambit
Le Gambit est un logiciel de maillage 2D/3D ; préprocesseur qui permet de mailler des
domaines de géométrie de différentes forme d’un problème de CFD (Computationnel Fluide
Dynamics). Il génère des fichiers (.Msh) pour Fluent. Ce dernier est un logiciel qui permet de
résoudre et simuler des problèmes de mécanique des fluides et de transferts thermiques par la
méthode des volumes finis. Le Gambit regroupe trois fonctions :
 Définition de la géométrie du problème.

 Le maillage et sa vérification.
 La définition des frontières (des conditions aux limites) et définitions des

domaines de calculs. [27]
70
Fig (Ⅲ.2) : Interface du logiciel GAMBIT.

Ⅲ.11 Maillage sous gambit
Pour réaliser un maillage d’une forme géométrique dans le logiciel Gambit on
effectuera les étapes suivantes :
Ⅲ.11.a Démarrage de gambit
Le chemin de l’application de Gambit est le suivant :/Fluent.Inc/ntbin/ntx86/Gambit.exe
Vous pouvez créer un raccourci dans la barre des tâches. S’il y a un problème
d’exécution, supprimez tous les fichiers.Mesh.lok dans le répertoire:/Fluent.Inc/ntbin/ntx86 et
relancez Gambit.exe.
Ⅲ.11.b Construction de la géométrie
Le but de la construction de la géométrie est de définir les domaines de calcul qui seront
des faces dans un problème 2D et des volumes dans un problème3D.
Pour effectuer une simulation par Gambit-Fluent il faut suivre une logique de calcul.
Cette simulation faite à partir des processus avec une approche interactive, étape par
étape.
71
Cette logique est résumée dans l’organigramme suivant :
Géométrie Maillage Simulation
Non Setup Solution Résultats

Structuré (Fluent) (Fluent) (Fluent)
Général Solution Plots

Méthode
Matériel Graphique
Solution
Model Contrôle
Moniteurs
Conditions
Aux limites
Solution
Valeurs
Lancer le
de
Calcul
référence
Fig(Ⅲ.3) :''Processus CFD' 'pour la simulation numérique pour Gambit-Fluent.
72
Fig (Ⅲ.4) :les opérations pour construire une géométrie.
Ⅲ.11.c Maillage
La génération du maillage (2D ou 3D) est une phase très importante dans une analyse
CFD, vu l’influence de ses paramètres sur la solution calculée. Ce menu permet de mailler en
particulier une ligne de la géométrie, à savoir disposer les nœuds avec des conditions
particulières (utilisation d’un ratio pour modifier la pondération du maillage, application de
formes différentes de maillage).
Ⅲ.11.c.1 Choix du type de maillage
Le nombre de cellules nécessaires pour résoudre un problème dépend de la précision

attendue. Il est important que le maillage soit resserré dans toutes les zones où les variables
présentent de forts gradients, en particulier au voisinage des parois. Afin de choisir un
maillage optimal pour obtenir des résultats numériques acceptables. On fait varier le nombre
de cellules et on observe l’évolution des résidus en fonction du nombre d’itérations.
Ⅲ.11.c.2Différents types de maillage

La méthode numérique des volumes finis est employée pour discrétiser le domaine
physique et les différentes équations aux dérivées partielles tridimensionnelles des modèles
mathématiques appliqués à l’écoulement considéré. Elle est basée sur l’intégration de ces
équations sur des volumes finis obtenus par découpage du champ de l’écoulement suivant un
maillage imposé par le logiciel Fluent. Selon la forme géométrique, il faut bien choisir une
maille adaptable à la géométrie, les mailles utilisées par Fluent sont les suivant : [28]
Tétraèdre Hexaèdre Pyramide

prisme
Fig(Ⅲ.5) : Exemple des maillages utilisées en volumes finis par Fluent.
Ⅲ.11.c.3 Maillage structuré (quadra/hexa)

Il est beaucoup plus facile de le générer en utilisant une géométrie à multi bloc, il
73
présente les avantages suivants :
 Economique en nombre d’éléments, présente un nombre inférieur de maille

parrapport à un maillage non structuré équivalent.
 Réduit les risques d’erreurs numériques car l’écoulement est aligné avec le
maillage.Ses inconvénients :
 Difficile à le générer dans le cas d’une géométrie complexe
 Difficile d’obtenir une bonne qualité de maillage pour certaines géométries

complexes
Ⅲ.11.c.4 Le maillage non structuré
Les éléments de ce type de maillage sont générés arbitrairement sans aucune contrainte
quant à leur disposition.
 Ses avantages :
• Peut-être généré sur une géométrie complexe tout en gardant une bonne
qualité deséléments.
• Les algorithmes de génération de ce type de maillage (tri/tétra) sont très
automatisés.
 Ses inconvénients :
• Très gourmand en nombre de mailles comparativement au maillage structuré
• Engendre des erreurs numériques (fausse diffusion) qui peuvent être plus
importantesi l’on compare avec le maillage structuré.
Fig(Ⅲ.6):types des maillages.
Ⅲ.11.c.5 Maillage hybride
Maillage généré par un mélange d’éléments de différents types, triangulaires ou

quadrilatéraux en2D, tétraédriques, prismatiques, ou pyramidaux en 3D.
74
 Ses avantages :
• Combine entre les avantages du maillage structuré et ceux du maillage non

structuré
Ⅲ.11.d Techniques générales de génération du maillage
Pratiquement, il n’existe pas de règle précise pour la création d’un maillage valable,
cependant il existe différentes approches qui permettent d’obtenir une grille acceptable. Nous
pouvons résumer ces règles ainsi :
- Maintenir une bonne Qualité des éléments.
- Assurer une bonne Résolution dans les régions à fort gradient.
- Assurer un bon lissage dans les zones de transition entre les parties à
maillage finet les parties à maillage grossier.
- Minimiser le nombre Total des éléments (temps de calcul raisonnable).
Ⅲ.11.e Qualité d’un maillage
La génération d’une très bonne qualité de maillage est essentielle pour l’obtention d’un
résultat de calcul précis, robuste et signifiant.
- Une bonne qualité de maillage repose sur :
- Une bonne résolution dans les régions présentant un fort gradient (couches
limites, ondes de choc…etc.).
Enfin, la qualité de maillage a un sérieux impact sur la convergence, la précision de la
solution et surtout sur le temps de calcul.
Ⅲ.11.f Lissage
Le changement dans la taille des éléments de maillage d’une zone maillée à un autre
doit être graduel, la variation de la taille des éléments de deux zones adjacentes ne doit pas
dépasser 20%.
Ⅲ.12 Présentation de Fluent
Fluent est un programme informatique conçu pour la simulation des écoulements de
fluide et du transfert de chaleur dans des géométries complexes.
Il présente une grande flexibilité des résultats et une adaptation avec n'importe quel type
de maillage, et il permet le raffinement du maillage en fonction des conditions aux limites, des
dimensions. Cela est très utile dans les régions à gradients importants.
Fluent écrit en langage ‘C’, il emploie toute la flexibilité et la puissance qu'offre ce

langage. Il utilise l'architecturé serveur client nécessaire au calcul parallèle sur plusieurs
75
machines. Fluent dispose d'un outil de graphisme pour l'affichage des résultats et leur
exploitation.
On peut aussi exporter les résultats vers un autre logiciel de graphisme, et l'option (User
Define Fonction) UDF permet de résoudre des équations additionnelles ou des termes sources
additionnels définis par l'utilisateur.
Ainsi, à titre non exhaustif, il a les capacités de modélisation suivantes :
 Ecoulements 2D ou 3D ;
 Régime permanents ou transitoires ;
 Ecoulement incompressibles ou incompressibles, non visqueux, laminaires
Ou turbulents, Newtonien ou non ; 

 Transfer de chaleur, par conduction, par convection ou radiatif ;
 Les écoulements avec changements de phases ;
 Ecoulements en milieu poreux.
Fig (Ⅲ.7): Schéma général qui représente la position du Fluent dans les étapes de calcul
Numérique.
Fluent est un code de calcul pour modéliser les écoulements des fluides et les transferts
thermiques dans ou autour des géométries complexes. Il peut résoudre des problèmes
d’émailles non structurées qui peuvent être produites pour des géométries complexes, avec
une facilité relative. Les types de mailles supportées sont des mailles, en 2D, triangulaires ou
quadrilatérales, ou en 3D tétraédriques, hexaédriques, pyramidales, et les mailles (hybrides)
76
mixtes.
Fluent supporte le langage de programmation ‘C’, et utilise pleinement la flexibilité et la
puissance offertes parce langage (allocation de la mémoire dynamique).
Fluent s’utilise à travers une interface graphique. L’utilisateur avancé peut adapter ou
augmenter selon les besoins, l’interface en écrivant des macros et des fonctions de menu, afin
d’automatiser certaines procédures.
Ⅲ.12.a Interface du code Fluent
On peut démarrer 4 versions de Fluent 2D, 3D, 2DDP, 3DDP ayants la même interface.
Fig (Ⅲ.8) : Interface de code Fluent.

 2D (2Dimensions)
 2DDP (2 Dimensions Double Précision) (utiliser dans notre travail)
 3D (3 Dimensions)
 3DDP (3DimensionsDoublePrécision)
Fig (Ⅲ.9) : Fenêtre de la résolution
77
Chapitre III Simulation numérique de la ventilation naturelle
- File : pour les opérations concernant les fichiers : lecture,
sauvegarde,importation…etc.
- Grid : pour la vérification et la manipulation du maillage et de géométrie.
- Define : pour définir les phénomènes physiques, les matériaux et les

conditionsaux limites.
- Solve : pour choisir les équations à résoudre, les schémas de

discrétisation, lesfacteurs de sous relaxation, les critères de convergence et pour initialiser et
procéder au calcul.
- Adapte : pour l'adaptation du maillage selon plusieurs paramètres.
- Surface : pour créer des points, de lignes et des plans nécessaires à

l'affichage desrésultats.
- Display et plot : pour l'exploitation des résultats.

- Report : pour afficher des rapports contenant les détails du problème traité.
- Parallèle : pour le calcul parallèle.
- Help : pour l'exploitation du contenu du code
Ⅲ.12.b Choix des paramètres de FLUENT
Ⅲ.12.b.1 Simples précisions ou double précision
FLUENT offre deux modes de calcul : le mode (double précision) et le mode (simple
précision). Dans le mode (double précision) : les nombres à virgule flottante sont représentés
en utilisant 64 bits, alors que le mode (simple précision) : utilise une représentation à 32 bits.
L’outre, le mode (double précision) est préconisé, pour les écoulements impliquant des
longueurs d’échelles très disparates, comme dans le cas d’un canal très longet mince.
Ⅲ.12.c Choix de la formulation du solveur
Sous(FLUENT), on peut choisir entre plusieurs formulations du solveur :

La formulation « ségrégâtes », ou isolée « implicite »: Cette formulation résout les
équations de continuité, de quantité de mouvement et quand c’est nécessaire celle de
l’énergie, séquentiellement, C’est-à-dire isolées les unes des autres « implicite par défaut ».
Le solveur isolé est classiquement employé pour les écoulements incompressibles à
modérément compressibles.
78
Ⅲ.12.b.2 La formulation (couple), ou couplée (implicite ou explicite) :
cette option permet aux équations gouvernantes d’être résolues simultanément, c’est- à-
dire couplées les unes avec les autres. Cependant, les autres scalaires, tels que les quantités de
la turbulence, sont traités isolément. Initialement, ce mode a été conçu pour les écoulements
compressibles à grandes vitesses. Ceci lui donne un avantage pour le traitement des
écoulements hautement couplés (forte interdépendance entre la densité, l’énergie et les
moments) avec des forces de volumes. Il faut signaler que le solveur couplé implicite requiert
presque le double de la mémoire qu’utiliserait le solveur isolé, alors que le solveur couplé
explicite vient au milieu, en terme de besoins en ressources, mais converge plus lentement
que la formulation implicite et n’est conseillé que pour les écoulements instationnaires.
Ⅲ.12.c Schémas de discrétisation

Sous « FLUENT », les variables stockées au centre de la cellule doivent êtres
interpolées aux faces du volume de contrôle. Il est possible de choisir entre différents schémas
de discrétisation pour les termes convectifs des équations gouvernantes, alors que les termes
visqueux sont automatiquement discrétisés au second ordre pour plus de précision. Il reste que
la discrétisation au premier ordre procure une meilleure convergence, alors que le « Second
Ordre Upwind Schème » est de rigueur pour les écoulements non alignés au maillage.
Ⅲ.12.d Choix de schéma d’interpolation de la pression
Dans la plupart des cas, le schéma « standard » est acceptable pour des écoulements
spécifiques.
On peut choisir par miles options suivantes :
 Le schéma « PRESTO » (Pressure staggering Option) est approprié pour les

écoulements dans des domaines fortement courbés.
 Le schéma au « second ordre » est à utiliser pour l’écoulement
compressible et pouraméliorer la précision en écoulements incompressibles.
 Le schéma linéaire « linear » est disponibles comme alternative dans le cas
ou les autres options ont des difficultés de convergence ou génèreraient des comportements
non physiques.
79
Ⅲ.12.d.1 Le schéma "PRESTO"
Le schéma "PRESTO" (pressure staggering option): est approprié pour les écoulements
hautement tourbillonnaires, à grande vitesse de rotation ou les écoulements dans des domaines
fortement courbés. Il utilise le bilan de masse dans un volume de contrôle décalé pour la face
ou est calculée la pression décalée. Ce schéma est inspiré de l'idée de la grille décalée valable
avec un maillage structuré proposé par Patankar (1980).
Ⅲ.13 Condition aux limites
Le tableau suivant représente les types et les caractéristiques de différente zone
Tableau (Ⅲ.1) : les paramètres de solveur (conditions aux limites).
Zone Types Caractéristiques

Fenêtre Velocity-inlet V=1(m/s)
ouest
Porte (nord) Outflow -
Fenêtre sud Outflow -
Mur ouest Wall ɸ= 53.82 W/m2
Mur est Wall T=297.16 K
Mur nord Wall ɸ =66.68W/m2
Mur sud Wall ɸ =71.39 W/m2
Paramètres de sous relaxation

Le tableau suivant représente les valeurs des différentes variables de sous relaxation :
Tableau (Ⅲ.2) : les valeurs des différentes variables de sous relaxation.
Variables Valeurs
Pression 0.3
Quantité de mouvement 0.7
Force de volume 1
Densité 1
Energie cinétique 0.8
turbulent
Dissipation turbulent 0.8
Viscosité turbulente 1
Energie 1
80
Algorithme de résolution
Le tableau (Ⅲ.3) présente les différents modèles utilisé sous Fluent.
Tableau (Ⅲ.3) :les différents modèles utilisés sous solveur Fluent.
Variables Least squares cell based

Accouplement pression–vitesse Simple
Pression standard
Quantité de mouvement First ordre upwind
Energie cinétique turbulent First ordre upwind
Dissipation turbulent First ordre upwind
Energie First ordre upwind
Ⅲ.14 Procédure de résolution
Après avoir vérifié le maillage de la géométrie en Gambit, on lance la lecture du fichier

contenant l’extension. msh. La création du fichier de donné (extension.cas) est nécessaire en
Fluent pour la vérification et la manipulation du maillage et de la géométrie ; puis on passe au
type de solveur. Notre cas est tridimensionnel et l'écoulement est stationnaire. Les autres
étapes définissent les choix du fluide, matériaux et les conditions aux limites.
La dernière étape, c'est choisir les équations à résoudre, les schémas de discrétisation,
les facteurs de sous relaxation, les critères de convergence ainsi que l'initialisation du calcul
qui dépend des conditions aux limites. [20]
Ⅲ.15 Contrôle de convergence
Pour s'assurer de la convergence des solutions, on présente la courbe des résidus (Figure
III.10) des différents paramètres de l’écoulement dans le local. Le logiciel fluent déclare la
convergence dés que les itérations de tous les résidus calculés atteignent la valeur du critère
de convergences introduites pour chaque paramètre. Le tableau suivant contient le résidu
compté pour chaque équation.
81
Fig ((Ⅲ.10) : la courbe des résidus (Pour le cas N°4).

Le tableau suivant résume les différentes valeurs des grandeurs représentées dans la
figure.
Tableau (Ⅲ.4) :Les résidus pour les différentes équations
Quantité du Quantité du Equation de Equation de k Epsilon ε

équation
mouvement mouvement continuité L’énergie
Selon x Selon y
Résiduel 10−4 10−4 10−4 10−4 10−4 10−4
Conclusion
A travers ce chapitre, on perçoit mieux les difficultés que l’on peut rencontrer lors de la
modélisation numérique d’un problème en mécanique des fluides même très simple.
En effet, les étapes de la modélisation numérique : la création de la génération du
maillage, l’implémentation du problème sous « Fluent », le contrôle de la convergence et
même le traitement des données de la simulation sont fortement interdépendants.
Par ailleurs, l’analyse des différents paramètres à régler qu’offre le logiciel « FLUENT»
nous a permis d’adopter une méthode de résolution optimale.
La simulation de tout problème formulé via des équations aux dérivés partielle est de
plus en plus répondue au niveau industriel.
IV Introduction :
Dans ce chapitre on a présenté les différentes géométries du local que nous allons
82
Chapitre IV Résultats et discussion
simulés au quatrième chapitre, tous les dimensions du locale sont en mètre, la longueur de
local en tous les cas est 13.16 mètres et la hauteur est 2.9 mètres. Une largeur de 5metres et
une superficie de 65.8 m².
IV.1 But de l’étude :
Le but de la présente mémoire d’une façon générale c’est l’étude de l’influence de la

Ventilation dans un local, car la température de l’air et le débit de chaleur jouant un rôle
très important dans la circulation de l’air dans les locaux.
L’étude est basée sur la variation du flux de chaleur et l’adaptation de la géométrie

du local pour la réalisation de la ventilation naturelle.
Il faut noter ici que pour des difficultés de simulation en 3D dans les logiciels de
CFD on a préféré de négliger la profondeur (troisième dimension) et réaliser seulement une
étude et simulation en 2DDP, car cette dernière donne des résultats très proche à celle
effectué en 3D et satisfaisante pour comprendre la circulation naturelle de l’air dans le
local.
IV .2 Cas de notre étude :
IV .2.1 séjour fermée :
Dans ce cas on a imposé :

83
Densité de Flux ; Mur ouest ɸ=53.82 W/m²

Mur sud
ɸ=71.39 W/m²
Mur nord
ɸ=66.68 W/m²
Fig (IV.1): la géométrie du local étudié

et simulé
Fig (IV.2) : courbe de résidus
84
Fig (IV.3): Contours des températures.
Fig. (IV.4) : diagramme des températures
IV.2.1.a Discussion des résultats
La température dans le local est presque uniforme verticalement à cause de

l’infiltration de flux pour chaque mur ; sur le côté du mur ouest en remarque que la
température et max a cause de flux qui est égale à ∅ = 2054 𝑤 par rapport au mur
nord qui et égale à ∅ = 967 𝑤 et le mur sud qui est égaleà ∅ = 1035.18 𝑤 , la
température varie du mur ouest au mur est. L’augmentation est estimé d’une moyen
de (10 °C/m).
IV .2.2 séjour avec fenêtre ouest (Ouverture mono façade)

85
Dans ce cas on a imposé des conditions aux limites :

La vitesse de l’air dans le régime laminaire V= 1 m/s La vitesse de l’air dans le régime
turbulent V= 2.1 m/sFlux ;
Mur ouest ɸ=45.9 W/m²

Mur sud ɸ=71.39 W/m²
Mur nord ɸ=66.68 W/m²
Fig (IV.5) : la géométrie du local à simuler
a- Régime laminaire b-Régime turbulent
Fig (IV.6): les courbes des residus
86
Fig (IV.7) : contour de la température
a- Régime laminaire
b- Régime turbulent
Fig (IV.8) : Diagramme de température

87
IV.2.2.a Discussions des résultats : la température
On remarque que la distribution de la température du local est homogène au centre est

élevée dansles coins parce que la ventilation ne couvre pas toutes les zones du séjour ;
Dans le régime turbulent la température et plus basse que celle de régime laminaire et ça
s’est due àla vitesse d’infiltration du vent plus la vitesse est grande plus la température va
baisser.
Nous observons que la zone de confort est devenue très grande que le cas
précédent et lachaleur placé sur la surface intérieure du séjour donc la solution est efficace
mais n’est pas optimal.
a-régime laminaire b- Régime turbulent

Fig (IV.9) : Contour de vitesse
Fig (IV.10) : Vecteurs vitesses
88
Fig (IV.11) : Diagrammes des vitesses
IV.2.2.b Discussions des résultats des vitesses
On peut observer la vitesse à l’entrée de la fenêtre est constante alors que dans le séjour la
vitesse est en changement considérable et elle varie selon la position ; l’écoulement de
l’air se sépare au milieu en formant deux tourbillons symétrique ce qui va assurer une
meilleure ventilation du localet ça afin de ressentir le confort voulu. Ceci correspond bel
et bien au comportement attendu.
D’après la figure (IV.11) on peut observer la différence de cette vitesse dans le régime
laminaire onremarque que la vitesse du vent est au max dans les deux tourbillons crée alors
que dans le régime turbulent la vitesse est au maximum à l’entrée de la fenêtre.
89
a- Régime laminaire b- Régime turbulent
Fig (IV.12) : Contour de la pression
Fig (IV.13) : Diagramme des pressions

IV.2.2.c Discussions des résultats de la pression
L’influence de la fenêtre est efficace est suffisant pour notre but et donne les meilleures
performances que celle de local fermée.
90
Nous avons besoin de modification sur le local et nous pensants qu’avoir plus
d’ouverturesest la solution parfaite pour dégager une bonne quantité de chaleur.
IV .2 .3 séjour ouvert par l’ouest (baie vitrée) et par le nord (porte) : Dans ce cas on a
imposédes conditions aux limites :
La vitesse de l’air dans le régime laminaire V=
1 m/s La vitesse de l’air dans le régime
turbulent V= 2.1 m/s
Mur sud ɸ=71.39 W/m²

Fig (IV.14) : Géométrie du local à simulé avec deux ouvertures côté ouest et nord
a-régime laminaire b- régime turbulent
Fig (IV.15) : les Courbes des résidus

91
a-régime laminaire b- régime turbulent
Fig (IV.16) : Contour de température
a-régime laminaire
b-régime turbulent

IV.2.3.a Discussion des résultats de température :
On remarque dans ce cas que l’effet de tirage thermique a une influence considérable
sur la température de séjour, le vent infiltrer par la fenêtre sort dans la direction du nord ce
92
qui donne un meilleur confort sur ce côté.
Dans le régime laminaire la zone de confort est plus importante que celle ressentie
dans le casturbulent vu que la vitesse d’infiltration et plus grande.
a-régime laminaire b-régime turbulent
Fig (IV.18) : Contour de vitesse
Fig (IV.19) : Vecteur vitesse
a-régime laminaire
93
b-régime turbulent
IV .2.3.b Discussion des résultats des vitesses
Dans ces figures on constate que le mouvement de circulation de l’air est suivant la
directiondu nord dans les deux cas.
Dans le régime laminaire la vitesse de l’air est faible ce qui donne une meilleure
propagation de ce dernier ; alors que le cas turbulent la propagation est faible et c’est du a
la grande vitesse d’infiltration.

Fig (IV.21) : profils des vitesses
a-régime laminaire
94
b-régime turbulent
4.3.c) Discussion des résultats
On voit que la pression et au maximum à la sortie de la porte ou elle atteint (p=175

Pascal) ; Dans ce cas l’effet de tirage thermique crée une zone de dépression dans le coin
nord-ouest.
Ici la figure montre que la pression n’est pas uniforme et surtout dans la zone nord-
ouest où il y a création d’une zone de dépression. L’ouverture nord montre son importance
dans la climatisation naturelle. Néanmoins, cette modification n’est pas suffisante pour des
températures très élevées donc on a besoin d’autres modifications pour améliore la zone de
confort et plus précisément pour la circulation, la ventilation et le renouvellement de l’air.
95
IV.2 .4 séjour ouvert par l’ouest (bais vitrée) et par le nord (porte), par le sud
(fenêtre) :Nous avons imposé des conditions aux limites :
La vitesse de l’air dans le régime laminaire V= 1 m/s La vitesse de l’air dans le régime
turbulent V= 2.1 m/s

Mur sud ɸ=68.51W/m²
La modification proposée se base sur la provocation du mouvement de l’air
(circulation d’air)par le principe de l’effet de tirage thermique. On réalise une ouverture sur
le mur nord et une autre dans le mur sud pour laisser l’air chaud cumulée (Pression et
température un peu élevés) à s’évacuée vers l’extérieur et plus la pression est grande plus
la vitesse de sortie de l’air augmente.
Fig (IV.23) : la géométrie simulée
Fig (IV.24) : les courbes des résidus

96
Fig (IV.25) : profils des températures

97
IV.2.4.a Discussions des résultats des températures
On constate que dans ce cas la température a diminué par rapport aux autres cas est ça s’est
due à lameilleure circulation de l’air dans le séjour ;
La densité et petite donc les molécules de fluide s’écartent donc température augmente
on serapprochant de la paroi ouest ; plus la vitesse de vent grandi plus la zone de confort est
large.
Fig (IV.27) : profils des températures
Fig (IV.28) : profil des vecteurs vitesses
98

IV.2.4. b Discussions des résultats des vitesses
D’après les figures (IV.30), nous remarquons que la vitesse est maximale à la sortie de la porte est à
la sortie de la fenêtre sud ; dans le régime turbulent la vitesse est double de celle de régime
laminaire mouvement de devient remarquable et elle est amortie à côté du mur est du local et séparé
en deux l’un suit la direction du nord (porte) l’autre part vers la direction du sud (fenêtre sud). C’est
un point d’amélioration recherché mais elle nécessite encore d’amélioration.

Fig (IV.30) : Contour de pression
99
a-Régime laminaire
b-Régime turbulent
IV.2.4.c Discussions des résultats de la pression
On voit que la pression et au max au voisinage de la sortie de la porte et la fenêtre ;

Dans ce cas l’effet de tirage thermique crée une zone de dépression dans le coin nord-ouest et
sud-ouest duséjour.
Ici les figures montrent que la pression n’est pas uniforme et surtout dans les zones de
dépressions ; Les ouvertures nord est sud montrent leurs importances dans la climatisation naturelle.
Conclusion :
Cette modification n’est pas suffisante pour des températures très élevées, on a
besoin d’autres modifications pour améliore la zone de confort et plus précisément pour la
circulation, la ventilation et le renouvellement de l’air. Néanmoins, c’est la meilleure des
initiatives pour un meilleur confort dans notre étude.
Envisagé des solutions alternatives comme une ventilation mécanique résoudra le problème
de l’inconfort dans tous les coins et recoins de notre séjour.
10
0
Conclusion générale
Conclusion :
Au cours de cette étude, on a d’abord entamé des notions de base de la ventilation et ces
différents types. On distingue trois mécanismes de ventilation, ventilation naturelle,
ventilation mécanique et ventilation mixte (une combinaison entre les deux).
La ventilation dans les locaux d’une façon générale elle permet l’amélioration des
conditions de confort telle que la circulation d’air, le renouvellement d’air et l’hygiène par
l’air neuf et purs.
La ventilation naturelle est une stratégie passive, sans moyen mécanique, de maintenir
un environnement intérieur confortable. Elle joue, également un grand rôle dans la réduction
de la consommation énergétique. Le principe de base de la ventilation est le phénomène de
convection naturelle due aux différences de température qui ont pour effet de provoquer le
mouvement de l’air donc tirage de l’air.
Dans ce présent travail, à l’aide de la simulation numérique basée sur la méthode des
volumes finis par le code Gambit et Fluent, nous avons montré l’importance de la ventilation
naturelle à travers des propositions de modification sur géométrie d’un locale simple pour but
d’améliorer le confort des occupants. Avec ces propositions on a réussi de réaliser les
améliorations suivantes :
Partant de la discrétisation du domaine d’étude (maillage de la géométrie), passant par

la modélisation numérique des écoulements des fluides. Le problème a été abordé selon une
approche numérique, en utilisant le code commercial CFD(FLUENT).
En été, les niveaux de température intérieure atteints dans un local non climatisé
dépendent principalement de la température extérieure, des apports internes de chaleur, des
apports de chaleur dus au soleil, de la ventilation et de l’inertie du local. En fonction de
l’orientation, Il s’agit essentiellement de réduire les apports de chaleur extérieure.
On a déterminé grâce à la simulation numérique les zones de confort et d’inconfort, et

on aillustré virtuellement les mouvements d’air dans notre local.
Enfin. On retiendra de cette étude que le local étudié, présente des insuffisances en
101
Conclusion générale
matière de prise en charge du confort intérieur ; pour remédier à ce problème d’insuffisance

on va opter pour une ventilation mixte ventilation mécanique et ventilation naturelle.
102
BIBLIOGRAPHIE
1- NEUF : « Climat intérieur/ confort, Santé, confort visuel » Revue européenne
d’architecture N°77, novembre - décembre 1978. P 12.
2- Liébard, A.et De Herde, A. « Traité d’architecture et d’urbanisme bioclimatiques ». Ed. Le

Moniteur. Paris 2005 p. 30a
3- Hugues Boivin « la ventilation naturelle Développement d’un outil d’évaluation du

potentiel de la climatisation passive et d’aide à la conception architecturale », Mémoire de
maitrise, université Laval Québec, 2007, p 115
4- Izard-L. « Architectures d'été construire : pour le confort d’été », Edition Edisud, 1994,
p141
5- Maalej, J., « Emetteurs de chaleur dans les bâtiments : comportement thermique et étude
des performances ». Thèse de Doctorat soutenue le 19 septembre 1994, l’Université de
Valenciennes,1994, p 151
6- Morinne, M. « Travail à la chaleur et confort thermique ». Les notes scientifiques et

techniques de l’INRS, NST 184, décembre 1999
7- Vogt J.J. « confort physiologique. Technique de l’ingénieur »...B2180. 1995.p 10
8- Hoffman J.B. « Ambiances climatisées et confort thermique » les actes du COSTIC. P110
9- Olsen B.W. « thermal confort requirement fir floors occuppied by people with bare feet »,
cite in: Mansouri Y. (2003), opcit, p273
10- Vogt J.J. « confort physiologique ». Techniques de l’ingénieur, document N° B 2180. 1995,
p 10
11- Deoux. S, « Le guide de l’habitat sain », Andorra : édition medieco, Avril 2002, p.211
12- Givoni, B. « L’Homme, L’Architecture et le Climat ». Edition, Le Moniteur, paris ,1978,
p39
13- Tixier, N. « De la notion de confort à la notion d’ambiance » in revue du laboratoire
cresson de l’école d’architecture de Grenoble et CNRS Ambiances architecturales et
urbaines, France. 2007
14- J. Koffi, Analyse Multicritère Des Stratégies De Ventilation En Maisons Individuelles »,
Thèse De Doctorat, université de la rochelle, France, (2009).
15- Hugues Boivin, la ventilation naturelle développement d’un outil d’évaluation du potentiel
de la climatisation passive et d’aide à la conception architecturale, maître ès sciences
(M.Sc.), université Laval Québec, (2007).
16- Adrien Dhalluin, Etude de stratégies de ventilation pour améliorer la qualité
environnementale intérieure et le confort des occupants en milieu scolaire, thèse de doctorat,
université de la Rochelle, 2013.
17- R. Bassiouny, N.S.A. Koura, An Analytical and Numerical Study of Solar Chimney Use for
Room Natural Ventilation, Energy and Buildings, vol. 40, pp. 865–873, (2008).
18- Modélisation numérique du transfert thermique, Méthode des volumes finis.
20- DTR C3.2 C3.4
21- Éric Goncalves, Résolution numérique, discrétisation des EDP et EDO, septembre (2005).
22- Jan.Halma, Tony Arts et Jaroslav Fort, numerical solution of unsteady and steady transonic
flow, Elsevier2003.
23- Z. Akchiche, Etude de comportement d’une cheminée solaire en vue de l’isolation
thermique, Mémoire de magistère, Université Kasdi Merbah, Ouargla, (2011).
24 André Lallemand, Dynamique des fluide réels, Techniques de l’Ingénieur, traité Génie
énergétiqueBE8157.
25- Mourad Gahgah, Influence des régimes d’écoulement sur les performances énergétiques
des cheminées solaires, Mémoire de magistère, Ecole Nationale Supérieure Polytechnique,
Alger, (2008).
26- Documentation FLUENT6.1.22 united states of America2003.
27 M. Slack. Cyclonic separator. Fluent Europe Ltd.GTC1-CT99-10030.
27- Documentation Fluent.http://courses.cit.cornell.edu/fluent.(16/05/2010).
Les références des figures
28- Figure I.1

https://energieplus-lesite.be/theories/confort11/le-confort-thermique-d1
29- Figure I.2
THELLIER, Françoise. « L’homme et son environnement thermique – Modélisation ».
Université de Paul Sabatier de Toulouse, p 65 (1999).
30- Figure I.(3-5)
https://www.google.com/search?q=la%20ventilation&tbm=isch&hl=fr&tbs=rimg
31- Figure I.(6-8)
https://www.souchier-boullet.com/prescription/prescription-gestion-energetique/ventilation-
naturelle-intelligente-vni/
32- Figure I.(9-12)
33- Figure (I.13)
http://www.saph-isolation-toiture.fr/VMC-double-flux.html
34- Figure I.(14-18)

Ventillation Naturelle 2

Transféré par

Droits d'auteur :

Formats disponibles

Ventillation Naturelle 2

Transféré par

Informations du document

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Ventillation Naturelle 2

Transféré par

Droits d'auteur :

Formats disponibles

République Algérienne Démocratique et Populaire

Ministère de l’Enseignement Supérieure et de la Recherche scientifique

MEMOIRE DE FIN D’ETUDE

Option : Froid Chauffage Et Climatisation

Encadré par : M. FOUCHANE SIDALI Réalisé par :

Devant le jury composé de :

Président : M. FERROUK Mohammed

Examinateur : M. BELTAGY Hani

Nous remercions avant tout Allah, de nous avoir donner la force

Nous tenons à remercier notre encadreur M. FOUCHANE SIDALI

d’avoir accepté d’encadrer notre travail et pour son assistance et sa contribution,

Symbole Signification Unité

K Coefficient de transmission surfacique W/m².

Kg Coefficient de transmission surfacique global W/m².

he Coefficient extérieur d’échanges superficiels W/m².°C

tbi Température intérieure de base °C

tbe Température extérieure de base °C

R Résistance thermique m².°C/W

c Capacité thermique massique J/kg.°C

z Différence de niveau entre le niveau de la face supérieure du m

Notion des apports

Symbole Signification Unité

Symbole Signification Unité

r Résistance thermique d’échanges superficiels....... m².°C/W

symboles Signification Unité

Le thème du présent travail vise à étudier le comportement thermique d’une maison

La ventilation naturelle est un élément clé d’une architecture traditionnelle et aussi un

L’expérience, la simulation numérique et l'utilisation de formules simplifiées permettant de

L’objectif du présent travail est premièrement calculer le bilan thermique de la maison,

Ce travail est composé de quatre chapitres.

La qualité de vie à l'intérieur de l'espace a été souvent rapprochée à une appréciation

I.2 Confort thermique :

Le confort c’est tout ce qui contribue au bien-être, a la commodité de la vie matérielle

I.2.2 Les paramètres affectant le confort thermique :

I.2.3 Paramètres liés à l’ambiance extérieure

I.2.3.a La température de l’air ambiant :

I.2.3.b La vitesse de l’air :

I.2.3.c L’humidité relative de l’air :

Fig (I.1): Polygone de confort hygrothermique, température et humidité[28].

I.2.4 Paramètres liés à l’individu :

00 0,1 0,3 0,5 0,8 1,0 1,5 3

Fig. (I.2) : valeurs exprimées en Clo des tenues vestimentaires[29].

I.2.4.b L’activité humaine :

L’activité humaine est un paramètre essentiel pour la sensation thermique de l’individu,

I.2.5 Paramètres liés aux gains thermiques internes :

des ressources internes autres que le système de chauffage.

Avec l’essor de la technologie et des besoins électrique (éclairage, électroménager,), les

Fig. (I.3) : Gains thermiques internes d’un espace[29].

Un inconfort local peut également être dû à des différences de températures

Le courant d’air local, l’asymétrie de la température de rayonnement, la différence

I.2.6.a Effet des courants d'air :

La perception d'un courant d'air localisé notamment au niveau de la nuque ou du visage

I.2.6.b Effet de l'asymétrie d'un rayonnement thermique :

I.2.6.c Effet de gradient thermique vertical de l'air :

I.2.6.d Effet de la température du sol :

* Température des murs : 22 ±2°C

* Humidité relative entre 30 et70%

* Température du sol : 19 à24°C

* Vitesse de l’air : inférieure à 0,20 m/s ;

I.2.7 L’évaluation de confort thermique :