Revue Phys. Appl.

23 (1988) 1161-1168 JUIN 1988, 1161

Classification
Physics Abstracts
44.00

Modélisation et simulation du fonctionnement d’un séchoir solaire

agroalimentaire
M. El Hisse et M. Fournier

IMP du CNRS, Université de Perpignan, Avenue de Villeneuve, 66025 Perpignan Cedex, France

(Reçu le 22 décembre 1987, révisé le 11 février 1988, accepté le 8 mars 1988)

Résumé. Un modèle mathématique simulant le fonctionnement d’un séchoir solaire agroalimentaire 2014
2014

modèle SERESOL basculant 2014 est proposé. Les paramètres mal connus ont été identifiés à partir
d’expériences menées en laboratoire dans un tunnel de séchage ventilé à atmosphère contrôlée. Le modèle a
ensuite été validé au cours d’une expérience réelle de séchage in situ dans le séchoir solaire. Des résultats de
simulations d’opérations de séchage sous différents climats sont alors présentés.
Abstract. A mathematic model for describing the behaviour of an agroalimentary solar dryer
2014

rocking 2014

model SERESOL - is proposed. The parameters used in the equations are derived from experimental
laboratory data, and the model has been valided by a real experience in situ. Then some results of simulations
of solar drying operations under different climates are presented.

Nomenclature

Majuscules :
AS : Air sec.
Aw : activité de l’eau (Sd)
CPa : chaleur massique à pression constante de l’air (J. kg-1 °C-1)
Cp f : chaleur massique à pression constante du fruit (J. kg-1 °c- 1)
CH : paramètre de l’équation
Dm : coefficient de diffusivité (m2s-1)
Df : diamètre de la tranche de fruits (m)
FF : facteur de forme (Sd)
Go : densité de débit massique d’air (kg. s-1 m- 2)
Ha : flux d’enthalpie molaire de l’air (W/mole)
HR : humidité relative de l’air (%)
Ma : débit massique de l’air (kg.s-1)
Ms : masse sèche du fruit (kg)
Pf : pression partielle de vapeur d’eau à la surface du produit (Pa)
P va : pression partielle de vapeur d’eau dans l’air (Pa)
P’s : pression partielle de vapeur saturante (Pa)
PT : pression totale (Pa)
PH : paramètre de l’équation
Sel : surface claie (m2)
Sech : surface d’échange des fruits (m2)
si : surface de l’isolant au niveau d’une claie (m2)
Ta : température de l’air (K)
Te : température extérieure (K)

Article published online by EDP Sciences and available at http://dx.doi.org/10.1051/rphysap:019880023060116100

1162

Tf : température du fruit (K)

V : volume (m3)
Vb : volume d’un boîtier (Sel x h) (m3)
Xa : humidité absolue de l’air (kg/kg.AS)

Minuscules :

a, b, c : coefficients à identifier caractéristiques du fruit étudié.

ca : conductance de l’isolant (WK-1)
ep : épaisseur de la couche de fruit (tranche) (m)
ei : épaisseur du caisson (paroi + isolant) (m)
h : distance entre deux claies (m)
r : rayon (coordonnées polaires, sphère ou cylindre) (m)
rv : constante des gaz parfaits (J.kg-1 K-1)
t : temps (s)
x : abscisse (coordonnée cartésienne) (m)
w : teneur en eau du produit base sèche (kglkg.Ms)
w : teneur en eau moyenne du produit base sèche (kg/kg.Ms)

Lettres grecques :

a T : coefficient de transfert de chaleur convectif surfacique air-produit (W.m- 2 °C-1)

03B1v : coefficient de transfert de chaleur convectif volumique air-produit (W.m-3 °C-1)
03B2m : coefficient de transfert de masse produit-air (m.s-1)
À : conductivité thermique (W. m-1 K-1)
àHv : chaleur latente de vaporisation de l’eau (J.kg-1)
/1p : variation de pression (Pa)
/1T : variation de température (K)
Ot : pas de temps (s)
Ox : pas d’espace (m)
p : masse volumique de l’air (kg.m-3)

Introduction. 1. Le séchoir solaire SERESOL : Modèle basculant.

Un prototype de séchoir solaire agroalimentaire a 1.1 DESCRIPTION ET FONCTIONNEMENT. - C’est

été conçu et réalisé par une entreprise régionale : il un séchoir solaire indirect (Fig. 1), le fruit à sécher
permet de sécher entre 10 et 20 kg de produit frais étant soumis à un transfert thermique uniquement
en 2 ou 3 jours. Cet appareil peut fonctionner en convectif. Il bénéficie d’une autonomie énergétique
convection naturelle ou en convection forcée : si le (séchoir passif) par son fonctionnement possible en
ventilateur est alimenté par photopiles, l’autonomie convection naturelle. Un ventilateur alimenté par un
du séchoir est assurée.
Nous proposons ici une modélisation mathémati-
que de ce séchoir qui permet de simuler son fonction-
nement dans les conditions les plus diverses ou bien
d’optimiser sa géométrie. L’originalité du travail
réside essentiellement dans la manière dont le pro-
blème a été appréhendé : le modèle décrit le fonc-
tionnement du séchoir solaire système macrosco-
-

pique - en tenant compte des phénomènes qui

interviennent au sein de la matière échelle micros-
-

copique.
Ce travail de modélisation est complété par une
application à un procédé thermique d’une technique
de l’automatique : l’identification paramétrique qui
Fig.1. -
Schéma du séchoir solaire SERESOL modèle
permet lorsque certains paramètres intervenant dans
basculant.
les équations sont malconnus de les identifier à
partir de mesures expérimentales. [Schema of solar rocking dryer SERESOL.]
 1163

panneau de photopiles lui conserve son autonomie la résolution par une méthode classique de triangula-
mais en régime forcé. risation permet de déduire les valeurs des 6 incon-
nues, en particulier la valeur de la température de
Il est constitué de 3 parties distinctes [1] : l’air de sortie des capteurs en fonction de la tempéra-
1) Une partie fixe constituant le caisson de dessi- ture de l’air à l’entrée et de l’ensoleillement reçu [2].
cation composée de six claies horizontales en acier
déployé inoxydable (surface totale des claies : 2. Modélisation du caisson de dessication.
0,16 m2 x 6 ~1 m2) l’arrière du caisson est isolé par
une feuille de polystyrène revêtue d’un film d’alumi-
2.1 DESCRIPTION DU PHÉNOMÈNE DE SÉCHAGE. -
nium.
Depuis les travaux de Sherwood [3] la théorie du
2) Une partie basculante destinée à opérer aisé- phénomène de séchage distingue la succession de 2
ment le chargement et le déchargement des claies ; phases essentielles, après une brève période de mise
une charnière assure la rotation. Sous la couverture en température.
transparente en plaque alvéolaire de polycarbonate - La phase I est une phase pendant laquelle la
on peut distinguer :
vitesse d’évaporation de l’eau interne du produit (ou
-

le bas (0,72 m2) qui est en fait l’absorbeur du

allure de séchage) est constante. Cette phase corres-
capteur solaire. Réalisé de polyuréthane
en mousse
pond à l’élimination de l’eau d’imprégnation faible-
émulsionnée (20 pores par inches) il est noirci dans ment liée à la matière solide - l’eau est libre et la
la masse. Lors de la fermeture du séchoir, par
surface du produit est saturée - : c’est un phéno-
rotation vers le bas, l’absorbeur poreux vient coiffer mène d’évaporation superficielle sous l’effet de
le dessus du caisson de dessication : des tendeurs l’écart entre la pression de vapeur d’eau dans une
latéraux, en caoutchouc, assurent par pression la couche mince d’air entourant le produit, et la
bonne étanchéité de l’ensemble.
-

le haut (0,50 m2) qui joue le rôle de cheminée

pression de vapeur d’eau dans l’air de séchage. Cette
solaire grâce à la présence d’un absorbeur métalli-
phase est souvent imperceptible pour de nombreux
produits dans les conditions habituelles de séchage.
que, recouvert d’une peinture noir mat à dépôt de -

La phase II commence à la fin de la période

carbone. Lors de la fermeture du séchoir, le bas de
cette cheminée se trouve dans le prolongement de la précédente lorsque la teneur en eau moyenne du
partie haute du caisson de dessication ; la section de produit atteint une valeur critique. A partir de cette
valeur le front de vaporisation de l’eau commence à
passage de l’air est en arrière de l’absorbeur métalli-
que.
régresser vers l’intérieur du produit. Le phénomène
basse du séchoir est fixé un autre
prépondérant est alors la diffusion interne de l’eau.
3) Dans la partie La modification de la surface effective de transfert et
capteur également de type absorbeur poreux, dis- la résistance rencontrée par la vapeur d’eau dans la
posé parallèlement au sol servant de préchauffage de couche externe provoquent alors un ralentissement
l’air (0,58 m2). de l’allure de séchage.
Un ventilateur placé à l’entrée de ce capteur
assure la circulation de l’air dans le séchoir. Les deux phénomènes physiques intervenant dans
chacune des phases (évaporation superficielle et
1.2 MODÉLISATION DES CAPTEURS SOLAIRES DE diffusion interne de l’eau) étant de nature différente,
TYPE ABSORBEUR POREUX. - Nous avons choisi un nous sommes amenés à proposer un jeu d’équations
modèle simple en ne considérant que les températu- par phase. Ces deux systèmes ne diffèrent en réalité
res moyennes. L’inertie thermique des capteurs est que par une équation : celle donnant l’expression de
faible donc un modèle écrit en régime permanent, la vitesse de séchage.
suffit pour représenter ce générateur d’air chaud.
La température de la matrice poreuse est supposée 2.2 ECRITURE DES SYSTÈMES D’ÉQUATIONS DU
uniforme et l’écoulement de l’air y est considéré MODÈLE. - Le caisson de dessication est composé
unidimensionnel. de 6 claies superposées (Fig. 2) que l’on peut sché-
L’épaisseur de cette matrice est telle que le matiser comme une succession de 2 boîtiers :
rayonnement ne la traverse pas (2 cm). -

boîtier plein contenant le produit à sécher,

un

L’écriture des bilans thermiques au niveau : balayé par l’air.

-
un boîtier vide situé dessus, dans lequel s’effec-
-

de la couverture vitrée
-
de l’air à l’entrée du capteur tue le mélange de l’air ascendant, sortant du boîtier
-

de l’air dans la matrice poreuse plein et l’air sec et chaud arrivant latéralement après
-

de la matière poreuse traversée de l’absorbeur poreux.

-

de l’isolant C’est ce mélange d’air qui aborde ensuite le

- de l’air à la sortie du capteur, conduit à un boîtier plein de la claie supérieure suivante et ainsi
système de 6 équations facilement linéarisables dont de suite.
1164

-
entrée de la couche de fruits : Ta’, Ma’,
Xa’
-
sortie de la couche de fruits : Ta2, Ma2,
Xa2
-
entrée latérale au-dessus des fruits : Ta°, Ma°,
Xa°
-
sortie de la claie après le mélange d’air : Ta,
Ma, Xa.

2.2.2.1 Modèle de la phase 1 ou phase évapora-

Fig.2. -

Schéma d’une claie constituée d’un boîtier plein toire.-

C’est la phase de début de séchage pendant
et d’un boîtier vide.
laquelle on enlève l’eau d’imprégnation faiblement
[Schema of a tray compound with full box and empty box.] liée au produit, et qui donc migre facilement sans
rencontrer de résistance à l’intérieur du produit.
2.2.1 Hypothèses. L’évaporation de l’eau, dans ce cas, est un phéno-
mène de surface, la répartition de la teneur en eau
-

L’air de séchage dans le produit reste uniforme. Dans le séchage d’un

-
Au niveau d’une claie, la température, l’humi-
dité et le débit de l’air sont homogènes et uniformes
produit biologique la durée de cette phase est brève
devant la phase II. Un séchage lent facilite sa mise
avant et après le mélange. La distribution du débit
en évidence.
d’air tout le long de la paroi-capteur du caisson est
déterminée en calculant le profil des vitesses d’air à
-
Bilans dans le boîtier plein entre la couche de
fruits et l’air de séchage (partie inférieure d’une
travers l’absorbeur poreux [2].
-
L’air humide et la vapeur d’eau sont considérés claie).
comme des gaz parfaits. Au niveau de l’air de séchage
-
L’inertie thermique de l’air est négligeable Bilan
devant celle des fruits.
thermique :
La variation de la puissance thermique de l’air
-

Les fruits pendant la traversée de la couche de fruits est égale à

-

La surface d’échange entre l’air et les fruits la somme des puissances échangées par convection
reste constante en cours de séchage. par l’air lors du balayage des fruits, et par conduction
- Tous les fruits d’une même claie, sont à la à travers la paroi isolante.
même température et leur teneur en eau initiale est
identique.
-

Les échanges
-
Les pertes par conduction sont calculées à Bilan massique :
travers le fond et les côtés du caisson. -

Conservation de la masse d’air qui traverse la

- Les échanges par rayonnement avec les fruits couche de fruits :
peuvent être négligés. En effet d’une part les échan-
ges avec l’absorbeur peuvent être estimés à moins de
2 % des apports convectifs, d’autre part les échanges -
Conservation de la masse d’eau à la sortie de
avec le fond, dont la température est très proche de cette couche.
celle des fruits, sont encore plus faibles. L’accroissement de la teneur en eau de l’air
-

Les profils internes d’humidité sont supposés provient de l’évaporation de l’eau du fruit :
identiques pour chaque fruit.
2.2.2 Formulation du modèle. - Un système
d’équations est écrit pour chaque phase du séchage,
le passage de la phase évaporatoire (I) à la phase Au niveau de la couche de fruits
diffusionnelle (II) correspond à une activité de l’eau Bilan thermique :
dans le fruit légèrement inférieure à l’unité (par La variation de l’enthalpie des fruits est égale à la
exemple Aw 0,99). =
somme des puissances échangées par convection
Les équations suivantes découlent des bilans ther- avec l’air et des puissances utilisées pour vaporiser
miques et massiques établis au niveau d’une seule l’eau des fruits :
claie, il en sera de même pour l’ensemble des autres
claies.
Nous adoptons les notations suivantes pour décrire
l’état de l’air (température, débit et humidité absolue
de l’air : T, M, X) à différents niveàux :
 1165

Bilan massique ou vitesse de séchage : Ceci permet de réduire de moitié la taille des
La variation d’évaporation est réglée essentielle- calculs. En tout instant la teneur en eau, en tout
ment par l’effet de pompe en surface, dont le terme point du fruit est ainsi déterminée en fonction de
moteur est exprimé par une différence de pression : l’échange convectif en surface. ,

Si on suppose que dans une couche toutes les

tranches de pomme ont une même répartition
d’humidité (tranches identiques et de même épais-
seur) la valeur moyenne de la teneur en eau de la
couche de produit, est obtenue en chaque instant par
-

Equation du mélange d’air dans le boîtier vide intégration sur le volume du fruit de la distribution
(partie supérieure d’une claie) de la teneur en eau interne. Seule cette valeur
-

Conservation des débits massiques d’air sec : moyenne est directement mesurable expérimentale-
ment

Conservation des débits massiques d’eau conte-

nue dans l’air :
2.3 DÉTERMINATION DES COEFFICIENTS APPARAIS-
SANT DANS LE MODÈLE.
-

Bilanenthalpique lors du mélange d’air : 2.3.1 Echange thermique par convection entre l’air et
L’enthalpie de l’air à la sortie de la claie est égale à les fruits.-

Les fruits coupés en tranches sont

la somme des enthalpies des courants d’air ascendant les claies.
disposés sur
et latéral, diminuée des pertes thermiques à travers Pour décrire d’une façon globale l’échange thermi-
la paroi isolante du fond du caisson.
que entre l’air et cette couche de fruits, nous avons
retenu la corrélation de Lof et Hawlay [4], donnant
le coefficient de transfert thermique volumique :
2.2.2.2 Modèle de la phase II ou phase diffusion-
nelle. -

En deçà de la teneur en eau critique du

fruit, wc,i le phénomène prépondérant est donc la
diffusion de l’eau interne. Le coefficient d’échange thermique surfacique a T
Le transfert de matière au sein du fruit peut alors est alors déduit simplement en utilisant la relation
être bien représenté par l’équation diffusionnelle de suivante :
FICK, appliquée à un cylindre.
En effet dans notre cas, une tranche de pomme
peut être facilement assimilée à un cylindre dont le
séchage ne se fait que par les faces planes, la surface 2.3.2 Conductance de l’isolant.
cylindrique (peau) présente une bonne étanchéité au
-

L’échange à tra-
flux d’eau. L’épaisseur de cette tranche est très
vers l’isolant fait intervenir sa conductance ca définie
inférieure à son diamètre. par :
Donc l’équation générale de transfert de masse
peut s’écrire dans le cas où le coefficient Dm ne
dépend que de la température du fruit.
On peut également tenir compte de la convection sur
la face extérieure de l’isolant, et utiliser un coeffi-
cient global de conductance.
avec les conditions initiales :
2.3.3 Echange massique entre l’air et les fruits. -

L’hypothèse de Lewis est souvent utilisée pour

caractériser le rapport du transfert de chaleur à celui
et les conditions aux limites aux bords de la tranche : de masse en introduisant la fonction de Lewis,
F(Le), définie par :

Une autre condition aux limites peut être déduite

de la symétrie de la tranche :
Pour un écoulement de l’air turbulent et des pres-
sions de vapeur négligeables devant la pression
totale ; la fonction F(Le) est voisine de l’unité.
1166

Nous poserons c = F (Le ) et nous identifierons cette Le modèle de séchage utilisé est donc un modèle
valeur. simplifié puisqu’il n’y a pas de mélange d’air latéral.
-

Pression partielle de vapeur saturante : Les valeurs des coefficients obtenus sont alors :
P,S est exprimée par la formule de J. Bertrand [5]
utilisable entre 0 °C et 200 °C.

Ces valeurs sont caractéristiques du produit (pom-

mes Golden taillées en tranche).
avec : Pys en atmosphère
2.5 RÉSOLUTION DU MODÈLE ET VALIDATION. -
Teneur en eau à l’équilibre :
-

Les systèmes d’équations établies pour chaque phase

Nous avons choisi la corrélation proposée par
de séchage sont discrétisés puis résolus par utilisation
Henderson [6] qui semble bien adaptée aux produits
de la méthode des différences finies.
biologiques : Les systèmes d’équations doivent donc être dou-
blement indicés pour tenir compte à la fois du temps
et du numéro de la claie ; l’équation de transfert de
masse à l’intérieur du produit en phase II (équation
CH et PH sont caractéristiques du produit et devront de diffusion) est triplement indicée, car on doit
faire l’objet d’une identification expérimentale.
rajouter la distribution de la teneur en eau à
2.4 DÉTERMINATION l’intérieur du fruit. D’autre part cette équation de
PAR IDENTIFICATION DES
diffusion et ses conditions aux limites sont exprimées
COEFFICIENTS MAL CONNUS. - La méthode d’iden-
en termes de différences finies suivant un schéma
tification paramétrique, méthode propre à l’automa-
est basée sur la minimisation de l’écart explicite ; les dérivées spatiales sont écrites sous
tique quadra- forme de différences centrées, les dérivées temporel-
tique entre les du modèle et celles du
réponses les sous forme de différences régressives.
processus réel, soumis à la même excitation.
La validation du modèle peut alors se faire en
La détermination de cinq paramètres, nécessite
comparant les résultats expérimentaux et ceux obte-
l’emploi d’une technique de descente (Powell) : nus par le modèle pour des conditions identiques.
- trois coefficients apparaissent dans l’expression L’expérience utilisée doit être différente de celle qui
des coefficients d’échange thermique et massique a, a servi à l’identification des paramètres.
b, c. Les résultats présentés sur la figure 3 sont satisfai-
-

deux dans l’expression de la teneur en eau à sants : les écarts décelés restent inférieurs à 2 %.
l’équilibre : CH et PH. La courbe présentée concerne les valeurs de perte
Les résultats expérimentaux permettant l’identifi- de masse au cours du temps. La validation aurait très
cation sont obtenus dans un tunnel de séchage bien pu être faite sur les valeurs des températures de
ventilé à atmosphère contrôlée. l’air ou des fruits.

Fig. 3. -
Validation du modèle sur une expérience de séchage en tunnel de séchage ventilé.

[The model validation on an apple drying experiment in a ventilated dryer-stove.]

 1167

3. Résultats de simulation. une montée trop forte en température du produit

risque de provoquer une détérioration des fruits en

Un programme de simulation est un outil à utilisation
les rendant ainsi inaptes à toute consommation.
multiple. Par exemple il convient parfaitement Sur les courbes de la figure 5, les pertes de masse
pour : obtenues les deux premiers jours d’une opération
-

étudier le comportement d’un système réel menée sous les deux climats les plus extrêmes (sec et
dans un environnement donné et prédire ses perfor-
humide) peuvent être comparées.
mances ; On constate dans cette figure 5 que la perte en eau
-
étudier l’influence d’un paramètre sur l’évolu- durant la nuit étant de toute façon très faible, il est
tion de la réponse du modèle : s’il s’agit d’un donc préférable d’arrêter la ventilation et de fermer
paramètre géométrique le design de l’appareil pourra l’entrée et la sortie du séchoir.
être optimisé à moins de frais que lors d’une Enfin, si la teneur en eau de l’air est élevée on doit
réalisation expérimentale. monter plus haut en température pour commencer à
Comme application nous nous limiterons ici à la sécher. Pendant une journée le produit sèche moins
présentation des résultats de simulation du fonction- longtemps, la perte de masse est plus faible et la
nement du séchoir sous deux environnements clima- durée totale du séchage est donc plus longue.
tiques différents :
-
un climat chaud et sec qui peut être illustré par Conclusion.
la région de Perpignan (France) au printemps ; Le modèle mathématique de simulation du fonction-
-
un climat chaud et humide, beaucoup plus nement du séchoir solaire étudié, a été validé grâce à
défavorable à une opération de séchage naturel mais une expérience de séchage de pommes Golden
qui justifie bien plus l’utilisation d’un séchoir solaire. coupées en tranche. La confrontation des résultats
Ce climat peut illustrer la réalité par exemple
théoriques et expérimentaux est satisfaisante, prou-
d’Abidjan (Côte d’Ivoire). vant ainsi que notre mise en équation des 2 phases
Le séchoir est rempli avec 6 kg de pommes successives de ce processus représente bien les 2
Golden coupées en tranches de 1 cm d’épaisseur. La phénomènes d’évaporation superficielle (Phase I) et
figure 4 simule le déroulement d’une opération de de diffusion interne (Phase II).
séchage à Perpignan sur quatre jours. Alors qu’en La technique d’identification paramétrique mise
début de séchage la température des fruits est très en oeuvre implique des expériences simples, effec-
différente de la température de l’air (jusqu’à 20 °C tuées dans un tunnel de séchage ventilé à atmosphère
d’écart), en fin de séchage on remarque une unifor- contrôlée [7].
misation de ces températures vers la température de Pour chaque produit à sécher un jeu de paramètres
l’air en sortie des capteurs. spécifiques, nécessaires au programme de simula-
D’ailleurs ceci doit être surveillé attentivement car tion, est ainsi déterminé.

Fig. 4. Simulation d’une opération de séchage de pommes à Perpignan de durée 4 jours. - - - - - Température des
-

fruits ; -. -. - Température de l’air à l’entrée de la claie ;..... Température de l’air à la sortie de la couche de fruits.
[Apple drying simulation on 4 days in Perpignan. - - - - - Fruit temperature ; - . - . - inlet air temperature ; - outler air
temperature.]
1168

Fig.5. Comparaison d’opérations de séchage de deux

-

jours simulées à Perpignan et à Abidjan. Perte de masse

cumulée pour un chargement de 6 kg de produit frais.

[Comparison of two drying simulation during 2 days in Perpignan and Abidjan. Cumulated loss mass of the fruits for a
6kg loading.]

Le modèle permet alors de simuler le séchage d’un Le modèle présenté ici, valable à tout débit, est
produit dans des conditions climatiques diverses. donc applicable et permettra le calcul des valeurs
La modélisation de ce séchoir en convection caractéristiques du séchage (température, teneur en
naturelle reste à faire : dans ce type de fonctionne- eau de l’air ou des fruits) avec la valeur du débit d’air
ment seul est modifié le débit d’air qui n’est plus fixe recalculée à chaque instant.
mais dépend de sa température et de son humidité à L’optimisation des paramètres géométriques ou
l’entrée et à la sortie du séchoir. de fonctionnement (débit d’air) peut aussi être un
objectif de ce travail.

Bibliographie

[1] FOURNIER, M., THEMELIN, A., Expérimentation de [4] LOF, G.O.G., HAWLEY, R. W., Unsteady state heat
séchoirs solaires à fruits de faible capacité : les transfert between air and loose solids, Ind. Eng.
modèles SERESOL et M5003, Séchage solaire Chem. 40 (1948) p. 106.
et développement rural, ENSAM-Bordeaux [5] BRUHAT, G., Thermodynamique (Masson, Paris)
(1983) pp. 105-121. 1968.
EL HISSE, M., Modélisation et simulation d’un séchoir
[6] HENDERSON, S. M., A basic concept of equilibrium
[2] moisture, Agricultural Engineering (1952) 33
solaire agroalimentaire, Thèse 3e cycle, Perpi-
p. 29-32.
gnan (1987). [7] PALLET, D., FOURNIER, M., THEMELIN, A., Modélisa-
[3] SHERWOOD, T. K., Drying of solids, Ind. Eng. Chem. tion, identification et simulation d’un séchoir
21 (1929) p. 12-16 ; 21 (1929) p. 976-980 ; 22 solaire à bois, Revue Phys. Appl. 22 (1987)
(1930) p. 132-137 ; 24 (1934) p. 307-310. p. 1399-1409.