See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/301627573

Modélisation du transfert de chaleur à travers un panneau sandwich destiné à

la construction: cas d’application par Comsol multiphysics

Conference Paper · January 2014

CITATIONS READS

0 2,858

1 author:

Abdellatif Zerizer
University M'Hamed Bougara of Boumerdes
59 PUBLICATIONS   188 CITATIONS   

SEE PROFILE

Some of the authors of this publication are also working on these related projects:

Traceability in wood production View project

characterization of algerians oak wood by x-ray tomographic scanner and free software image-j View project

All content following this page was uploaded by Abdellatif Zerizer on 04 October 2017.

The user has requested enhancement of the downloaded file.

 Septième édition
du COlloque FRancophone en Energie, Environnement, Economie et Thermodynamique -
COFRET'14
Paris, CNAM - 23 - 24 - 25 avril 2014

COFRET'14-OF1-147

MODÉLISATION DU TRANSFERT DE CHALEUR À TRAVERS UN

PANNEAU SANDWICH DESTINÉ À LA CONSTRUCTION: CAS
D’APPLICATION PAR COMSOL MULTIPHYSICS

LIMAM, Amel, Abdel, ZERIZER a)* b)

a)b)
Université M’Hamed Bougara de Boumerdès UR-MPE Avenue de l’Indépendance
35000 Boumerdès, Algérie.
*limammaha@yahoo.fr

Résumé
Les activités énergétiques de l’humanité ont atteint un niveau tel, qu’elles perturbent
significativement l’environnement, notamment dans le secteur de la construction qui
représente en Algérie près de 45,50 % de la consommation nationale d’énergie (produits
pétroliers, électricité, gaz naturel…) [1]. Face à cette situation, il semble urgent de proposer
des solutions alternatives aux méthodes de constructions actuelles où le béton et l’acier sont
les principaux matériaux de base, ce qui permettra de mieux envisager un réel développement
durable.
Les constructions sont des structures assez
complexes présentant beaucoup de contraintes en
raison des défauts de géométrie s’alternant avec
différentes conditions limites ; c'est pourquoi, la
simulation spécifique de transfert de chaleur et de
masse dans des constructions est devenue une des
priorités de développement [2].
Dans le présent travail et moyennant l’outil
logiciel Comsol Multiphysics, une modélisation du
transfert de chaleur à travers un panneau sandwich
a été proposée pour le cas d’application d’un
composite constitué de deux peaux en bois du Pin d’Alep collées à une âme en liège
agglomérée noir. Une caractérisation thermique des deux composants de panneau Pin d’Alep
et liège a précédé la modélisation à partir de laquelle on a déterminé les facteurs demandés
tels que la conductivité thermique, la chaleur spécifique…pour les insérer sur l’interface
Comsol.
Il faut souligner que les deux ressources forestières, le chêne-liège (Quercus Suber) et le pin
d’Alep (Pinus halepensis) occupent une bonne place des zones boisées en Algérie [3].
Les résultats de la modélisation obtenue ont montré une bonne isolation thermique du
composite dont les résultats révèlent un gradient de température tout au long de panneau
sandwich durant le passage de flux de la chaleur thermique.
La figure ci-dessus présente le modèle du sandwich proposé pour une âme de 46cm
d’épaisseur.
Mots-clés : Modélisation, panneau sandwich, Comsol multiphysics, Consommation
d’énergie, secteur de construction

1
 Septième édition
du COlloque FRancophone en Energie, Environnement, Economie et Thermodynamique -
COFRET'14
Paris, CNAM - 23 - 24 - 25 avril 2014

Introduction
Durant cette décennie, l’urgence climatique (énergie, production, consommation) et la
nécessité d’un développement durable occupent une place majeure dans tous les débats
mondiaux.
C’est probablement la prise de conscience, finalement récente, que nous vivons sur une petite
planète aux ressources limitées et à l’équilibre fragile, qui a conduit les citoyens du monde à
progressivement placer en priorité le souci environnemental. [4] Le défi actuel consiste à
offrir un réel développement durable à l’humanité, Ces questions primordiales entraînent les
décideurs vers des solutions alternatives aux matériaux habituellement utilisés en construction
(béton et acier).
Selon nos sources bibliographiques [2], l’utilisation du bois dans la construction remonte à
près de 20 000 ans. A cette époque, les constructions étaient à base de matériaux bio-sources
(bois, paille, chanvre, roseaux) lesquelles reviennent en force ces dernières années. Le
transfert de chaleur et de masse est l’un des phénomènes physiques les plus étudiés.
Quoique les constructions à base de bois, isolées en liège soient des structures assez
complexes qui présentent beaucoup de contraintes en raison de défauts géométriques
s’alternant avec différentes conditions limites, la simulation spécifique de transfert de chaleur
et de masse devient une des priorités de calcul [5].
Dans le présent travail, une modélisation du transfert de chaleur à travers un panneau
sandwich en bois de Pin d’Alep et de liège est proposée pour le cas d’application d’un mur
constructif de 1m de longueur soumis aux conditions réelles (rayonnement solaire, convection
à l’intérieure et à l’extérieur). Les résultats obtenus nous ont permis de voir l’influence des
différents paramètres sur la température intérieure et de voir ainsi la meilleure configuration
pour une bonne isolation thermique et donc une économie d’énergie.

1 – Conception et modélisation sous COMSOL MULTIPHYSICS

La modélisation numérique effectuée dans ce travail a été réalisée à l’aide du logiciel

COMSOL MULTIPHYSICS 0.4, basé sur la méthode des éléments finis. Ce logiciel permet
de résoudre différents problèmes physiques ainsi que le choix des conditions aux limites à
travers ces module contenant les équations prédéfinies pour la modélisation du transfert de
chaleur, de la matière (mécanique de fluide, système diphasiques, diffusion, …), les équations
aux dérivées partielles quelconques. L’avantage majeur de COMSOL et la possibilité de
réaliser le couplage entre plusieurs phénomènes physiques de nature différente. [6]
COMSOL multiphysique fournit deux sortes de modes d'opération, le style d'interface
utilisateur graphique est le style de commande en créant des scénarios. Les deux modes
fournissent la commodité pour des utilisateurs surtout. [7]
En quelques étapes, on peut créer un modèle numérique COMSOL, il se fait comme suit :
- L’étape préalable consiste en la création de la géométrie et le choix du maillage.
- La deuxième étape concerne le choix des modules de travail, l’introduction des
constantes et des relations usuelles, la définition des propriétés physiques des
matériaux et le choix des conditions initiales et aux limites. Dans notre travail on a
utilisé le module de transfert de chaleur à travers un panneau sandwich en bois et liège
et un panneau entièrement en bois afin d’évaluer leurs comportements thermiques.
- L’étape de résolution consiste en un choix du schéma de résolution, par exemple les
dimensions et la forme des mailles dont les calculs seront effectués. Dans notre cas, on

2
 Septième édition
du COlloque FRancophone en Energie, Environnement, Economie et Thermodynamique -
COFRET'14
Paris, CNAM - 23 - 24 - 25 avril 2014

a utilisé un maillage standard, le choix du solveur et la tolérance d’approximation du

calcul et de l’exécution du calcul.
- L’étape de post-traitement consiste en l'extraction des résultats de modélisation sous
forme d’images ou de graphiques et leurs validations. Il permet d’exploiter les
résultats (visualisation des champs de vitesse, température, …).

1.1 Détermination des paramètres technologiques à insérer sur l’interface Comsol

On considère le panneau sandwich constitué de trois couches, la couche centrale est en liège
aggloméré noir, d’épaisseur e2 couverte des deux côtés par des peaux en bois de Pin d’Alep
interne et externe d’épaisseur identique e1 et e3. Cette dernière est soumise à une condition de
flux solaire. Les échanges par convection et rayonnement avec les milieux environnants sont
pris en compte sur les deux surfaces interne et externe.
Différents paramètres mécaniques, physiques et thermiques doivent être insérés sur l’interface
de l’utilisateur Comsol. La forme de la structure sandwich fabriquée est présentée sur la figure
01.

Figure 01 : Forme de la structure sandwich fabriquée.

1.1.1 Paramètres propres au liège aggloméré noir

1.1.1.1 La conductivité thermique (λ) et la chaleur spécifique(Cp)

Suivant la norme ISO 8894-1, la mesure de la conductivité thermique et de la chaleur

spécifique a été effectuée en utilisant l’appareil Conductivité mètre (CTmétre) (figure 02)
en régime transitoire grâce à l'utilisation des sondes à choc thermique. L'échantillon de
dimensions (200 * 100 *20) mm est excité par un signal d'entrée obtenu par un élément
chauffant, la réponse du matériau sera recueillie en utilisant un couple thermoélectrique
affichée sur l’écran de système d’acquisition des données.

Figure 02 : Appareil de mesure de la conductivité thermique et la chaleur spécifique (CT mètre).

3
 Septième édition
du COlloque FRancophone en Energie, Environnement, Economie et Thermodynamique -
COFRET'14
Paris, CNAM - 23 - 24 - 25 avril 2014

La valeur de la conductivité thermique de liège aggloméré noir obtenue est de 0.04 W/m.K et
sa Chaleur spécifique est de 1670 J/Kg. °C.

1.1.2.2 La masse volumique (Mv)

Suivant la norme NF B 57-052, les éprouvettes sont découpées en forme de parallélépipède
rectangle, de dimensions (30* 30) cm ayant des faces planes et parallèles. A la température et
à l’humidité ambiante, on a mesuré les dimensions des éprouvettes à l’aide d’une règle,
ensuite on les a pesées en utilisant une balance électronique d’une précision de 10-2. Suivant
cette formule on a obtenu une moyenne de 120 Kg/m3 avec un écart type de 0.04:
ρ = (m / v )*103
m : la masse de l’éprouvette en (g).
V : le volume de l’éprouvette en (cm3).

1.1.2 Paramètres propres aux peaux du Pin d’ Alep

1.1.2.1 Le module de Young (E)

En utilisant la méthode non destructive ‘BING’ (Figure 01), Son principe s’appuie sur
l’évolution des vibrations sur une pièce de bois normalisée ayant les dimensions suivantes
(20, 20, 360) mm. On a identifié les fréquences propres de ces éprouvettes à partir de leurs
réponses à une excitation impulsionnelle, cette excitation avait été appliquée à une des
extrémités de façon à solliciter simultanément tous les modes propres de vibration. C’est par
une analyse du signal acoustique enregistré qu'on a déterminé une ou plusieurs fréquences de
vibration du système afin d’en déduire les propriétés mécaniques de l'échantillon testé.

Figure 01 : Dispositif pour la détermination du module d’élasticité par la méthode ‘BING’.

La valeur moyenne du module d’élasticité obtenue après un teste sur 18 éprouvettes est de
12690 MPa avec un écart type de 1630.

1.1.2.2 La masse volumique (Mv)

En suivant la norme NF B 51-005, les éprouvettes en bois de Pin d’Alep sont coupées à des
dimensions de (20*20*20)mm ensuite pesées en utilisant une balance électronique d’une
précision de 10-2, puis mesurées dans les trois directions par un pied à coulisse numérique à

4
 Septième édition
du COlloque FRancophone en Energie, Environnement, Economie et Thermodynamique -
COFRET'14
Paris, CNAM - 23 - 24 - 25 avril 2014

12% d’humidité. Suivant cette formule on a obtenue une moyenne de 580 Kg/m3 avec un
écart type de 0.033: MV = M / V

1.1.2.3 La conductivité thermique (λ) et la chaleur spécifique(Cp)

Suivant la norme ISO 8894-1, à l’aide d’un CT mètre, la valeur de la conductivité thermique
de bois de pin d’Alep obtenue est de 0.19 W/m.K et sa Chaleur spécifique est de 2500 J/Kg
°C.

1.1.2.4 Coefficients d’échanges convectifs internes et externes (h)

Par définition, un bilan thermique s’écrit de la manière suivante :

φ reçu par la plaque = φ perdu par la plaque (2)

Figure 03 : Schéma du bilan thermique réalisé sur le panneau sandwich (bois-liège).

 φ reçu par le panneau sandwich: elle reçoit donc le flux solaire φsol .
 φ perdu par le panneau sandwich:
Par convection : Afin de déterminer ces deux paramètres, on a caractérisé le flux de chaleur
φconv interne et externe qui est extrait par l’air ambiant (température T f2) du la peau externe à
la température Tp par :
φconv = h.(Tp - Tf2) (2)
La valeur de h dépend du type de convection (forcée ou naturelle). Le vent étant supposé nul,
la convection sera naturelle. Dans ce cas, une corrélation empirique permet d’écrire que :
hext = 1,42. ((Tp – Tf2)/L ) 0,25
hint = 1,42. ((Tp – Tint)/L ) 0,25
Avec L =dimension linéaire caractéristique de la surface d’échange, dans ce cas c’est la
longueur de sandwich qui est égale à 1m.
Au finale :
φconv = 1,42.(Tp - Tf2)1,25 (3)

Par rayonnement : C’est l’échange de la face extérieure en bois avec le ciel (ou
l’atmosphère, considérée comme un corps noir). Le flux que la peau à la température Tp
échange par rayonnement avec le « ciel » à la température T∞ vaut :
φray = v..(Tp4 - T4) (4)

5
 Septième édition
du COlloque FRancophone en Energie, Environnement, Economie et Thermodynamique -
COFRET'14
Paris, CNAM - 23 - 24 - 25 avril 2014

εv = 0.9 émissivité du bois1 ;

T∞ = 0,0552 .Tf21,5
 Φ à travers le panneau: il concerne le flux conductif qui passe à travers la plaque de
bois il est obtenu par la formule suivante :
cond= (Tp-Tint) / ((e1 / λ1)+(e2/ λ2)+(e3/ λ1)) (5)
En remplaçant les termes de l’égalité (1), et les formules (3), (4) et (5) on obtient comme
bilan thermique :
sol+cond = 2.conv + ray
sol+(Tp-Tint)/((e1/λ1)+(e2/λ2)+(e3/λ1))= 2.[1,42.(Tp - Tf2)1,25] + v..(Tp4 - T4) (6)
De cette relation, on calcul la température de la paroi Tp et l’on remplace
sa valeur pour calculer les coefficients d’échanges convectifs ;

Tableau1: Tableau récapitulatif des paramètres technologiques et géométriques à incérer sur

l’interface Comsol
Paramètres géométriques Symboles Valeurs
Longueur de panneau (m) L 01
Largeur de panneau (m) l 0.08
Epaisseur des deux peaux (m) e 0.02
Epaisseur de liège (m) e1 0.046
Epaisseur de panneau (m) ep 0.086
Paramètres technologiques Symboles Valeurs
Module de Young bois (MPa) E 12690
Masse volumique bois (Kg/m3) Mv 580
Conductivité thermique bois (w/m.K) λ 0.19
Chaleur spécifique bois (j/Kg. °C) Cp 2500
Masse volumique liège (Kg/m3) Mv 120
Conductivité thermique liège (w/m.K) λ 0.04
Chaleur spécifique liège (j/Kg. °C) Cp 1670
Coefficient d’échange convectif externe (W/m2.K) hext 15
Coefficient d’échange convectif interne (W/m2.K) hint 20
Température extérieure estivale de base (°C) Tamb 32

2 – Résultats et discussion
On a réalisé deux simulations, les deux configurations sont les suivantes :
- Un panneau sandwich ayant une épaisseur d’isolant de 46mm ;
- Un autre panneau sans isolant entièrement en bois ;
Les maillages obtenus des deux configurations sont présentés sur la figure 04 et les résultats
des simulations (figure 05) pour les mêmes conditions aux limites :
- On tient on compte les phénomènes de convection, rayonnement et conduction ;
- Quantité de chaleur reçue de 200 W/m2 ;
- Température extérieure de base de 32°C.

1
http://fr.wikipedia.org/wiki/%C3%89missivit%C3%A9

6
 Septième édition
du COlloque FRancophone en Energie, Environnement, Economie et Thermodynamique -
COFRET'14
Paris, CNAM - 23 - 24 - 25 avril 2014

Panneau sandwich bois-liège

Panneau bois
Figure 04 : Maillages des deux structures sous comsol.

Panneau sandwich bois-liège Panneau bois

Figure 05 : Model des transferts de chaleur à travers un panneau de bois et un panneau sandwich bois-
liège.
Après une supposition d’isolation des faces hautes et basses des deux panneaux et d’après la
figure 05, On constate que pour une température extérieure de base qui atteint les 32°C, on
remarque un gradient de température décroissant d’une manière progressive de l’extérieur
vers l’intérieur
 pour le panneau sandwich bois-liège, la valeur de la température de la peau extérieure
est de 39°C pour une quantité de chaleur reçue de 200 W/m2, la température intérieure
obtenue est de 33°C, on remarque sur le graphe de la figure 06 qui présente la
variation de la température en fonction de la position de panneau sur l’abscisse (x)
deux comportements décroissant de profil de la température de l’extérieur à l’intérieur,

7
 Septième édition
du COlloque FRancophone en Energie, Environnement, Economie et Thermodynamique -
COFRET'14
Paris, CNAM - 23 - 24 - 25 avril 2014

comportement durant le passage en bois et en liège. Ce gradient considérable (6°C

entre l’intérieur et l’extérieur) est dû à l’efficacité de l’isolant naturel et écologique
utilisé « le liège », ce dernier composé d’une structure cellulaire alvéole possède des
boules d’air statique, considérés comme un bon isolant thermique dans cette structure.

 Pour le panneau en bois, la température intérieure obtenue est de 36°C, on remarque

sur le graphe de la figure06 qui présente la variation de la température en fonction de
la position de panneau toujours sur l’abscisse(x) un comportement décroissant de
profil de la température durant son passage à travers le panneau bois, la différence de
température entre l’intérieur et l’extérieur et de 3°C, cette valeur justifie la nécessité
de l’utilisation des isolants dans les structures en bois destinées à la construction
surtout dans les climats chauds que possède l’Algérie.

Panneau sandwich
Panneau de bois
Figure 06: Profil de la température à travers les deux structures.

Conclusion

Dans cet article, nous avons présenté la modélisation de transfert de la chaleur à travers les
panneaux sandwich destinés à la construction et la comparés avec le transfert de la chaleur à
travers un panneau purement en bois.
La détermination des paramètres technologiques (mécaniques, thermiques et physiques) et
géométriques des deux structures a précédé la modélisation aux éléments finis sous comsol,
L’utilisation de liège (isolant naturel) dans cette structure nous à permis d’avoir un gradient
de température considérable entre l’intérieur et l’extérieur, ce résultat favorise l’utilisation de
liège aggloméré comme isolant dans les constructions.

Perspectives

Les résultats obtenus encourageants pour la suite du développement du système initié à

travers lequel on envisage de simuler et suivre le comportement en saison estivale et
hivernale toute une construction en panneaux sandwich bois-liège;

8
 Septième édition
du COlloque FRancophone en Energie, Environnement, Economie et Thermodynamique -
COFRET'14
Paris, CNAM - 23 - 24 - 25 avril 2014

Références bibliographiques

[1] APRUE, Consommation Energétique Finale de l’Algérie-Chiffres clés Année 2005,

Edition 2007.
[2] V. Gerlich, K. Sulovská and M. Zálešák, COMSOL Multiphysics validation as simulation
software for heat transfer calculation in buildings: Building simulation software validation,
journal of Measurement, Czech Republic, measurement 46 (2013) 2003–2012.
[3] A. Djema, M. Messaoudene, The Algerian Forest: Current Situation and Prospects, 1st
Annual Meeting of the Modelling, Valuing and Managing Mediterranean Forest Ecosystems
for Non-Timber Goods and Services, Palencia-SPAIN, OCT 26-27, 2007.
[4] B. Multon, Y. Thiaux, H. Ben ahmed, Consommation d'énergie, ressources énergétiques et
place de l'électricité, dossier Techniques de l’Ingénieur, 2011.
[5] Suarez, V., et al., Study of the heat transfer in solids using infrared photothermal
radiometry and simulation by COMSOL Multiphysics. Appl. Radiat. Isotopes (2013),
http://dx.doi.org/10.1016/j.apradiso.2013.04.010i.
[6] I.TOMASHCHUK, Assemblage hétérogène cuivre-inox et TA6V-inox par faisceaux de
haute énergie : compréhension et modélisation des phénomènes physico-chimiques, Une thèse
de doctorat. L’université de Bourgogn, 07 octobre 2010.
[7] Xiaolong Wang, Hui Yue, Guangliang Liu, Zhao Zhao, The Application of COMSOL
Multiphysics in Direct Current Method Forward Modeling, 2011 Xi an International
Conference on Fine Geological Exploration and Groundwater & Gas Hazards Control in Coal
Mines, Procedia Earth and Planetary Science 3 ( 2011 ) 266 – 272.

View publication stats