Traduit de Anglais vers Français - www.onlinedoctranslator.

com

ETUDE DES CARACTERISTIQUES DE TRANSFERT DE CHALEUR ET DE MASSE DANS LE

PROCESSUS DE SECHAGE D'UN SECHOIR A TAMBOUR BASEE SUR LA METHODE COUPLEE DEM-

CFD Jiawen Hou1, Dongdong WANG1, Fan LI1, Hongxia ZHAO1, Jiangtao ZHANG2, Lijie JIN2*,

Chu-anZHANG3, Dening XIANG3*,Ya-nan CHEN4, Xuehong WU4*

1 Usine de cigarettes Henan Zhongyan Zhumadian，Zhumadian, 463002，Henan, Chine

2Usine de cigarettes de Baoding, China Tobacco Hebei Industrial Co., Ltd., Baoding 071000, Chine ;
3Usine de cigarettes de Wuhan, China Tobacco Hubei Industrial Co., Ltd., Wuhan 430051, Chine ;
4École d'ingénierie énergétique et électrique, Université de l'industrie légère de Zhengzhou, Zhengzhou,
Henan, Chine, 450 000
[email protected]@[email protected]
Afin d'analyser les caractéristiques de transfert de chaleur et de masse du tabac coupé

dans un séchoir à tambour horizontal, la méthode des éléments discrets (DEM) a été utilis

ée pour calculer le modèle de collision de particules et le modèle de transfert de chaleur

couplé DEM-CFD. Avec l'augmentation de la vitesse de rotation, le degré de mélange

augmente progressivement, et lorsque la vitesse de rotation est de 16 tr/min, le degré de

mélange est plus élevé. Lorsque le flux de chaleur du couplage EDEM-CFD est de 0,2 W-1

W, la température des particules augmente progressivement et tend vers une valeur fixe.

Le temps de séjour du matériau a une grande influence sur la qualité de séchage des

particules de matériau. Si le temps de séjour est trop court, la chaleur du matériau sera in

égale, et si le temps de séjour du matériau est trop long, les dommages aux particules de

matériau augmenteront, augmentant le taux d'écrasement.

Mots clés :Couplage DEM-CFD ; transfert de chaleur et de masse ; tabac coupé, séchoir à tambour

1. Préface

Les sécheurs à tambour sont largement utilisés dans le séchage et le traitement des matériaux dans des industries telles

que la chimie, la pharmacie et l'alimentation. En raison de son amélioration significative de l'efficacité de séchage et de la qualité

des matériaux, il a été largement favorisé[1-4]. Quel que soit le type de tambour utilisé, il est directement lié aux vitesses de

rotation du tambour, et la vitesse appropriée du tambour est la base du mélange des particules ou du traitement ultérieur[5, 6].

Le séchoir à tambour peut remuer et cuire le tabac, de sorte que l'agitation et la cuisson efficaces peuvent
favoriser l'échange de chaleur entre les particules de tabac et l'uniformité du mélange. À une certaine température
de paroi, la conduction thermique entre la paroi du tambour et le tabac dépend principalement de la vitesse
d'agitation, qui affecte également l'uniformité des particules[7]. Le mouvement du tabac dans le séchoir est que le
tabac est ramassé et déposé en continu par la plaque intérieure, et le tabac est transporté vers l'avant p
ériodiquement. En même temps, le flux d'air interagit avec le tabac qui tombe, ce qui est facile à transporter et ne
présente pas le phénomène du tabac coincé dans le coin mort, etc., le séchage est plus uniforme[8]. Le
mouvement des particules de tabac dans le séchoir à tambour est très

1
complexe, le mécanisme de mouvement des particules dépend d'une série de facteurs d'état et des caract
éristiques du tabac lui-même.
Avec le développement de la technologie informatique, la simulation numérique est devenue une méthode

efficace pour étudier le mouvement des particules. À l'heure actuelle, le modèle d'écoulement multiphasé d'Euler et le

modèle d'éléments discrets (DEM) sont les principales méthodes de recherche[9]. Le DEM discrétise l'espace de solution

en un réseau d'éléments discrets et relie les éléments adjacents à des éléments de connexion raisonnables en fonction de

problèmes réels ; calcule la force résultante et le moment de force entre une unité et d'autres unités dans diverses

directions, ainsi que les forces externes causées par d'autres champs physiques agissant sur l'unité. Selon la deuxième loi

du mouvement de Newton, l'accélération de l'unité peut être obtenue et intégrée dans le temps pour obtenir la vitesse et

le déplacement de l'élément, et pour obtenir des quantités physiques telles que la vitesse, l'accélération, la vitesse

angulaire, le déplacement linéaire et l'angle de rotation de toutes les unités à tout moment. Lors du calcul du transfert de

chaleur des particules de tabac, de l'air chaud est passé dans le tambour et les paramètres du gaz sont définis par

FLUENT. L'opération est relativement simple et la phase solide, c'est-à-dire la phase particulaire, est définie par le logiciel

EDEM. Français Les deux logiciels réalisent le calcul de couplage grâce au modèle DDPM dans le logiciel FLUENT[10-14].

Scherer et al.[15] ont étudié le séchage par convection de copeaux de bois dans un séchoir rotatif de laboratoire à

chicanes en utilisant des simulations DEM-CFD couplées. Wu et al. [16] ont utilisé le modèle de transfert de chaleur coupl

é DEM-CFD pour étudier la simulation tridimensionnelle de l'écoulement gaz-solide dans un lit fluidisé à particules

zonales flexible. Le comportement dynamique des particules en ruban flexible est étudié de manière vivante. Darius

Markauskas et al.[17] ont utilisé la méthode DEM-CFD pour analyser le processus de transport des particules fibreuses

flexibles.

Le séchage du tabac coupé a un effet important sur la qualité du tabac coupé, mais le processus de séchage est un

processus physique complexe, qui implique trois processus de transfert de chaleur, le processus de transfert de chaleur et de

masse du tabac coupé dans le séchoir à tambour est la clé pour contrôler le processus de séchage du tabac coupé. Les particules

de tabac sont relativement fines et ont un certain degré de flexibilité, et une seule méthode de calcul ne peut pas obtenir

efficacement le mouvement des particules et le processus de transfert de chaleur et de masse. Par conséquent, cet article basé

sur la méthode DEM-CFD pour analyser les caractéristiques de transfert de chaleur et de masse du processus de séchage du

tabac, les résultats de la recherche fournissent des orientations théoriques pour le séchage du tabac coupé.

2. Modèle numérique

2.1. Modèle de collision de particules

L'objectif de cet article est les particules de tabac de la machine de séchage à rouleaux horizontaux à angle droit. La

méthode couplée CFD-DEM est appliquée dans le présent document, les particules sont modélisées au moyen de la th

éorie DEM et traitées comme de véritables particules discrètes. La force d'interaction des particules doit être ajoutée sur

les sphères dans le cadre de la deuxième loi de Newton lorsque des simulations couplées DEM-CFD sont effectuées. Étant

donné que les particules de tabac sont similaires au modèle de balle molle, le modèle de balle molle est généralement le

modèle de contact antidérapant Hertz-Mindlin, qui analyse la collision entre les particules à travers la vibration

d'amortissement de l'oscillateur à ressort simplifié. Les forces normales et tangentielles du modèle sont déterminées

respectivement par la formule (1) et (2)[18] :

(1)
{ ， } (2)

2
 Où: , - coefficients d'élasticité normaux et tangentiels des particules ; , - la quantité de
chevauchement entre la normale et la tangentielle ; , - coefficients d'amortissement normaux et tangentiels ; , -
vitesses normales et tangentielles relatives entre les particules.
Le mouvement des particules est décrit à l'aide des équations de la deuxième loi du mouvement de Newton telles que (3)
et (4) [9, 11] :

∑ ( ) (3)
∑ ( ) (4)
La masse de - i particules dans la formule ; Vitesse de - i particules ; Vitesse angulaire de - je
particules;

- le moment d'inertie de la particule i ; - le diamètre de la particule i ; N - le nombre de collisions

particules; , , - force normale, force tangentielle et moment de frottement de roulement entre les particules i et j.

La méthode des différences centrales est utilisée pour intégrer numériquement les équations de vitesse et de d

éplacement pour (3) et (4), obtenir de manière itérative la nouvelle valeur de déplacement de la particule, puis substituer

la nouvelle valeur de déplacement dans la relation force-déplacement pour obtenir la nouvelle force, de sorte que le

calcul est répété pour suivre le mouvement de la particule à tout moment.

A cela s'ajoute le frottement de roulement, qui peut être résolu par le moment de la surface de contact :

(5)
Où: est le facteur de frottement de roulement ; est la distance entre le centre de masse et le contact

point; est le vecteur vitesse angulaire unitaire de l'objet au point de contact.

2.2. Modèle de transfert de chaleur couplé DEM-CFD

Le modèle de contact Hertz-Mindlin avec conduction thermique a été choisi en raison de l'échange thermique entre

les particules de tabac et le séchoir à tambour. Le processus de transfert thermique dans le cylindre comprend

principalement le transfert thermique entre l'air chaud et les particules sur la paroi intérieure du cylindre de l'échangeur

thermique et les particules de matériau. Le logiciel EDEM ne peut calculer que le transfert thermique entre la paroi int

érieure du cylindre et les particules de matériau, tandis que le transfert thermique par convection entre l'air et les

particules nécessite le couplage d'EDEM et de FLEUNT, et les données de flux thermique sont importées dans l'EDEM

pour chauffer les particules en parallèle.

2.2.1 Modèle de transfert de chaleur par convection

Le transfert de chaleur par convection se produit principalement entre la phase gazeuse et la phase particulaire. En raison

de la différence de température entre la phase gazeuse et la phase particulaire, les deux phases effectuent un transfert de

chaleur, dont l'expression est illustrée par l'équation (6)[19].

（6）
Où, p est la phase particulaire ; f est la phase liquide ; est le flux de transfert de chaleur ; est-ce la chaleur

coefficient de transfert des deux phases. est la surface de la phase particulaire; diff est le
érence de température entre les deux phases.
Le coefficient de transfert de chaleur est calculé par l'équation (7).

（7）
Où, est la conductivité thermique de la phase gazeuse ; est le diamètre superficiel des particules ; est
Numéro de Nussel.

3
2.2.2 Modèle de transfert de chaleur par rayonnement

Pour le modèle de transfert de chaleur radiatif, l'émissivité de la surface des particules doit être définie et la
mise à jour de la température doit être définie dans EDEM pour garantir que la température des particules est
mise à jour. Selon la loi de Stefan-Boltzmann, le transfert de chaleur entre les particules et l'environnement
environnant peut être décrit par l'équation (8) [19].
（8）
Où σ est la constante de rayonnement du corps noir ; est l'émissivité d'une particule sphérique ; est
la température moyenne des particules dans une cellule de grille fermée.

2.2.3 Modèle de conduction thermique de Hertz Mindlin

Pour le modèle de conduction thermique, la simulation EDEM de l'échange de chaleur par flux de particules dans un

équipement rotatif fournit une méthode d'échange de chaleur entre particules, et son expression est présentée dans l'équation

(9) et (10)[20].
（9）

√ （10）
Où, est la conductivité thermique des particules ; est la composante normale du contact
force. est le rayon effectif de la partie en contact ; est le module d'élasticité effectif de la particule.

2.2.4 Transfert thermique transversal entre la géométrie et les particules

Le processus de transfert entre la surface de la paroi du cylindre et les particules est décrit par l'équation
(11)[21] :
（11）
Où est la chaleur transférée par la géométrie à la particule ; T est la température de la
géométrie ou particule; est la conductivité thermique entre la géométrie et les particules.

2.3. Modèle de calcul

Les interactions de données entre CFD et DEM sont mises en œuvre en calculant les forces. La modélisation du

tabac coupé par DEM se fait au niveau de la particule individuelle, tandis que le transfert de chaleur par CFD se fait au

niveau du maillage de calcul. Tout d'abord, Fluent est utilisé pour calculer les informations du champ d'écoulement à un

pas de temps, puis EDEM est lancé pour la même itération temporelle. L'interface de couplage est utilisée pour transférer

la position, le mouvement, le volume, la température et d'autres informations des particules à Fluent, et l'interaction

entre les particules et le fluide est calculée. L'effet du fluide sur les particules sera transmis à EDEM via l'interface sous la

forme d'une force volumique de particule affectant le mouvement des particules, tandis que l'effet sur le fluide agira sur

le fluide via la phase source d'impulsion. Dans le modèle CFD-DEM, le tabac coupé est considéré comme une phase discr

ète, et la deuxième loi de Newton est utilisée pour établir les équations de mouvement de chaque particule.

Le tabac coupé est supposé être composé de sept particules de 5 mm de diamètre avec une longueur en ligne

droite d'environ 23 mm, comme le montre la figure 1. Le modèle géométrique utilisé dans cet article est un séchoir à

tambour horizontal. Le diamètre du tambour est de 1960 mm, la longueur est de 9000 mm, l'épaisseur de la paroi du

tambour est d'environ 50 mm et l'angle entre le tambour et le plan horizontal est de 2°. Comme le montre la figure 2, le s

échoir à tambour est principalement composé d'une chambre avant, d'un cylindre et d'une chambre arrière. À l'entrée de

la chambre avant, une usine de particules est installée pour remplir le cylindre de tabac coupé

4
particules à une vitesse initiale de 2 m/s. les particules définies peuvent générer 10 000 particules sous la forme du

nombre de particules, et la vitesse de formation des particules est de 5 000 particules/s. La vitesse de rotation du cylindre

est réglée à 6 tr/min pendant la formation des particules, et après la formation des particules, la vitesse de rotation du

cylindre est réglée à 9 tr/min, 11 tr/min, 13 tr/min, 15 tr/min, 16 tr/min respectivement. En fait, la quantité de remplissage

des particules est de 50 %, mais en raison de la taille plus grande du tambour, la quantité de calcul est trop importante

pour effectuer un calcul stable. Par conséquent, dans cet article, la quantité de remplissage des particules est fixée à 2 %,

ce qui améliore l'efficacité du calcul, et l'influence de la vitesse du tambour sur l'état de mouvement et la température des

particules peut être obtenue. Les paramètres utilisés dans la simulation sont présentés dans le tableau 1.

Fig. 1 Taille du modèle de particules

Fig. 2 Modèle géométrique du séchoir à tambour (a) la chambre avant (b) le corps du cylindre (c) l'arrière
chambre du tambour respectivement.

Tableau 1 Paramètres de simulation

Propriété de paramètre argument Valeur numérique
Propriétés des particules Densité des particules kg·m-3(densité d'emballage) 125
Module de cisaillement /Pa 1e+07
Coefficient de Poisson 0,25
Propriété du cylindre Densité kg·m-3 7800
Module de cisaillement /Pa 7e+10
Coefficient de Poisson 0,3
Interaction Coefficient de frottement statique entre particule et cylindre 0,5
Coefficient de frottement dynamique entre particule et cylindre 0,01
Coefficient de récupération dynamique de la particule et du cylindre 0,2
Coefficient de frottement statique particule sur particule 0,408
Coefficient de frottement dynamique particule sur particule 0,355
Coefficient de récupération après collision particule sur particule 0,27
Autres paramètres Accélération de la gravité g/(m/s)2) Angle de roulis 9.81
2°
Pas de temps /s 1.8e-05

3. Résultats et analyse

3.1. Degré de mélange des particules

Les particules de tabac dans le séchoir à tambour sont continuellement frites en se déplaçant le long de l'axe du

tambour. Le processus de mélange a une grande influence sur le transfert de chaleur et de masse et sur la

5
L'uniformité du processus de séchage est telle qu'il est nécessaire d'étudier le processus de mélange. Grâce au post-

traitement EDEM, le mélange des particules de tabac peut être observé directement et le nombre de contacts entre les

particules de tabac et les particules de tabac à différents moments peut être enregistré. Un groupe de bacs est établi

dans le post-traitement EDEM pour enregistrer le nombre de contacts entre les particules de tabac et la paroi et la

plaque du tambour. Selon la méthode proposée par Gupta pour décrire le degré de mélange des particules par le

nombre de contacts des particules[22], le rapport entre le nombre de contacts entre les particules de tabac et les

particules de tabac et le nombre total de contacts entre les particules de tabac a été utilisé pour caractériser la dispersion

des particules de tabac au niveau de la paroi et de la plaque du cylindre, la valeur est calculée par le rapport de s

éparation Q entre les particules de tabac coupé et la paroi et la plaque du tambour :

(12)
Dans la formule, est le nombre de contacts entre les particules de tabac et les particules de tabac ; est
le nombre total de contacts (le nombre total de contacts entre les particules de tabac et la paroi du tambour et
l'aube directrice du tambour, et le nombre total de contacts entre les particules de tabac).
Le degré de dispersion des particules de tabac peut être analysé en comparant la taille du rapport de séparation q

entre les particules de tabac et la paroi du tambour et la plaque à aubes directrices. Plus la valeur est grande, plus les

particules de tabac sont proches et plus la surface de contact avec le tambour et la plaque à aubes directrices est petite,

ce qui n'est pas propice à l'uniformité du fil de cuisson. Par analyse statistique du nombre de contacts des particules de

tabac dans le processus de séchage au tambour à air chaud, la courbe du rapport de séparation q entre les particules de

tabac, la paroi du tambour et la plaque à aubes directrices telle que représentée sur la figure 3 change avec le temps

peut être obtenue. La figure 3f montre le rapport de séparation de la machine de séchage à tambour à différentes

vitesses dans un délai de 16 s. On peut voir que le rapport de séparation diminue progressivement avec l'augmentation

de la vitesse de rotation. Lorsque la vitesse de rotation est de 16 tr/min, le rapport de séparation est proche de 0,5 et le

nombre de contacts entre les particules de tabac et la paroi du tambour et la plaque à aubes directrices du tambour est

similaire à celui entre les particules de tabac. C'est-à-dire que le degré de mélange est élevé.
0,8
（un） 0,7（b）
0,7
0,6
0,6 y = Intercepter + B1*x^1 + B2*x^2 + y = Intercepter + B1*x^1 + B2*x^2
équation B3*x^3 + B4*x^4 + B5*x^5 + B6*x^6 + + B3*x^3 + B4*x^4 + B5*x^5 +
poids
B7*x^7
0,5 B6*x^6 + B7*x^7
0,5 dessin 9 tours
Taux de séparation/%

intercepter 11 tours
Taux de séparation/%

poids Non pondéré B1 Non pondéré

0,4
intercepter - 0,07506 ± 0,01284
0,46072 ± 0,03214
0,4 B2 - 0,08317 ± 0,01883
0,46945 ± 0,0448
B1 B3
B2 - 0,092 ± 0,02549 B4 - 0,07597 ± 0,03371

0,3 B3 0,00446 ± 0,00887

8,40341E-4 ± 0,00157
0,3 B5 - 0,00757 ± 0,0111
0,00357 ± 0,00186
B4 B6
- 1,30623E-4 ± 1,48004E-4 B7 - 4,13451E-4 ± 1,64855E-4
B5

0,2 B6 6,83887E-6 ± 7,04955E-6

- 1,27245E-7 ± 1,33555E-7
0,2 Somme des carrés des r
ésidus R-carré (COD)
2,0701E-5 ± 7,403E-6
- 3,88684E-7 ± 1,32183E-7
B7
R-carré ajusté 0,04629
Somme des carrés des résidus 0,03206
0,1
0,97304
0,1 0,988 R au carré (COD)
R-carré (COD)
R-carré ajusté 0,97042
R-carré ajusté 0,98709
11 tours
0,0 9 tours 0,0 Ajustement polynomial à 11 tr/min
Ajustement polynomial à 9 tr/min
- 0,1 - 0,1
0 2 4 6 8 10 12 14 16 0 2 4 6 8 10 12 14 16
Fois Fois

0,7
0,7 (c)
(d)
0,6
0,6
y = Intercepter + B1*x^1 + B2*x^2 + B3* x^3 + 0,5 y = Intercepter + B1*x^1 + B2*x̂ 2 +
B3*x^3 + B4*x̂ 4 + B5*x̂ 5 + B6*x̂
0,5
équation
équation B4*x^4 + B5*x^5 + B6*x^6 + B7*x 6 + B7*x^7
7̂
dessin 15 tours
dessin
0,4
13 tours
poids Non pondéré

0,4
Taux de séparation/%

poids Non pondéré

- 0,0958 ± 0,01807
Taux de séparation/%

intercepter
intercepter - 0,06776 ± 0,01354
B1 0,54607 ± 0,04119
B1 0,46295 ± 0,03267
B2 - 0,15725 ± 0,02974
- 0,11746 ± 0,02498
0,3
B2
B3 0,01524 ± 0,00839 0,3 B3 0,02089 ± 0,00943
- 0,00127 ± 0,00152
- 0,00112 ± 0,00143 B4
B4
B5 2,00501E-5 ± 1,30411E-4
B5 5,00875E-5 ± 1,30249E-4
0,2 B6 - 1,43731E-6 ± 5,98553E-6
0,2
B6
B7
1,14966E-6 ± 5,65738E-6
- 3,84875E-8 ± 9,76195E-8
B7 2,26119E-8 ± 1,09408E-7
Somme des carrés des r 0,06321
Somme des carrés des r 0,03477
0,1 ésidus R-carré (COD) 0,96002
ésidus R-carré (COD) 0,98305
R-carré ajusté 0,98175
0,1 R-carré ajusté 0,95698

13 tours 15 tours
0,0 Ajustement polynomial à 13 tr/min
0,0 Ajustement polynomial à 15 tr/min

- 0,1
0 2 4 6 8 10 12 14 16 0 2 4 6 8 10 12 14 16
Fois Fois

6
 0,7 0,8
(f) (f)
0,6 0,7

0,6
0,5
y = Intercepter + B1*x 1̂ + B2*x^2 + B3*x^3
équation + B4*x^4 + B5*x^5 + B6*x̂ 0,5

Taux de séparation/%
0,4 6 + B7*x^7

Taux de séparation/%
dessin 16 tours

poids Non pondéré 0,4

- 0,06797 ± 0,01874
0,3 intercepter
B1 0,48203 ± 0,04457
9 tours
B2 - 0,12643 ± 0,03353 0,3
B3 0,01268 ± 0,01104 11 tours
0,2 B4 2,0459E-5 ± 0,00185
B5 - 9,64988E-5 ± 1,64003E-4 0,2 13 tours
B6 6,60986E-6 ± 7,36474E-6

0,1 B7
Somme des carrés des r
- 1,40686E-7 ± 1,315E-7
0,04581 0,1
15 tours

ésidus R-carré (COD) 0,96044 16 tours

0,0 R-carré ajusté 0,9566

16 tours 0,0
Ajustement polynomial à 16 tr/min
- 0,1 - 0,1
0 2 4 6 8 10 12 14 16 0 2 4 6 8 10 12 14 16
Fois Fois

Fig. 3 Courbe du taux de séparation du tabac en fonction de la vitesse du rouleau

3.2. Vitesse, temps et distance des particules de tabac

Français Étant donné que toutes les particules ont un coefficient de frottement statique et dynamique égal, le temps mis

par la particule pour se déplacer le long de l'axe dans le même déplacement et le changement de la vitesse de la particule pour

se déplacer le long de l'axe ont été étudiés dans le cas de particules avec seulement la gravité et le frottement. Selon l'analyse du

processus de séchage du séchoir à tambour, lorsque l'angle du tambour est de 2°, la relation entre la vitesse du tambour et le

temps de séchage est illustrée à la Figure 4. On peut conclure du diagramme que le temps de séchage des particules de tabac

coupé dans le tambour diminue linéairement avec l'augmentation de la vitesse du tambour, le temps de séchage était

respectivement de 205,18 s, 175 s, 171,33 s, 150 s et 140,05 s lorsque la vitesse du tambour était de 9 tr/min, 11 tr/min, 13 tr/min,

15 tr/min et 16 tr/min. La courbe de changement dynamique de la vitesse axiale du mouvement des particules est illustrée sur la

Figure 5, et l'image de nuage de changement dynamique de la vitesse radiale du mouvement des particules est illustrée sur la

Figure 6, lorsque le matériau est soulevé par la plaque, la vitesse des particules augmente lorsque le matériau tombe de la

plaque, de sorte que la vitesse a une certaine périodicité, et la courbe de vitesse moyenne n'a pas de tendance évidente à la

hausse ou à la baisse, par conséquent, elle montre également que le mouvement des matériaux a atteint un état relativement

stable. À 16 tr/min, la vitesse axiale moyenne des particules est d'environ 1,75 m/s. La vitesse du tambour affecte le temps de r

étention du tabac coupé dans le tambour, l'augmentation de la vitesse du tambour entraînera la réduction du temps de rétention

du tabac coupé dans le tambour, au contraire, le temps de rétention du tabac à faible vitesse du tambour, par conséquent, la

vitesse du tambour affectera directement le temps de séchage du tabac coupé, affectant ainsi l'effet de séchage. En augmentant

la vitesse du tambour pour atteindre l'objectif d'augmentation de l'efficacité du séchage du tabac coupé, ou en contrôlant la

vitesse de séchage de la qualité du tabac coupé.

250 3.5

205.18s 3.0
200 9 tours
175s 2.5 13 tours

16 tours
Années 150171,33 s
150 2.0
Vitesse/(m/s)

140,05 s
1,5
Fois

100
1.0
9 tours
50 11 tours 0,5
13 tours
15 tours 0,0
0 16 tours

- 0,5
0 2000 4000 6000 8000 0 50 100 150 200
Distance/mm Fois

Fig. 4 Temps de séchage à différentes vitesses de rotation Fig.5 Vitesses axiales à différents régimes de rotation
vitesses

7
 Fig. 6 Distribution de la vitesse des particules dans la direction radiale du tambour

3.3. C
M
engourdi

voir
direct
à travers

rapport

est un pari

（un）9 tours （b）11 tours （c）13 tours

8
 （d）15 tours （et）16 tours

Fig.7 Plage de rideau de matériaux formés par des particules à différentes vitesses de rotation

3.4. Changement de température des particules de tabac

Le modèle de transfert de chaleur couplé DEM-CFD est construit en utilisant l'interface de couplage DDPM, et la vitesse de

rotation du cylindre est fixée à 16 tr/min. Le modèle d'Euler de couplage particules-fluide a été sélectionné dans FLUENT, et l'air

chaud de 113°Ca été introduit dans l'entrée de matériau. Le modèle de transfert de chaleur Hertz-Mindlin et la mise à jour de

température sont sélectionnés dans EDEM. Le modèle de transfert de chaleur de l'EDEM lui-même est comparé au couplage DEM-

CFD, et son flux de chaleur et sa température sont présentés dans les figures 8 et 9. Dans la figure 8, EDEM -0 est le flux de

chaleur de l'EDEM lorsque l'EDEM est calculé séparément, le flux de chaleur de l'EDEM lorsque l'EDEM-CFD-0 est couplé et le flux

de chaleur de CFD lorsque l'EDEM-CFD-1 est couplé. La figure 10 montre le changement de température des particules lorsque le

couplage EDEM-CFD et l'EDEM seul sont calculés. Étant donné que le flux de chaleur généré par EDEM seul est négligeable

comme le montre la figure 8(a), la température des particules change à peine. Comme le montre la figure 9(b), le flux de chaleur

généré par le couplage EDEM-CFD est compris entre 0,2 W et 1 W, de sorte que la température des particules va augmenter

progressivement, mais le taux de changement de température va diminuer progressivement et la température des particules

aura tendance à rester constante. Une image en nuage des changements de température des particules sur 80 s est présentée

dans la figure 10.

1.2
0,0010
(un) (b)
1.0
EDEM-CFD-1
0,0008
0,8
Flux de chaleur/W

Flux de chaleur/W

0,0006 0,6
EDEM-0
EDEM-CFD-0
0,4
0,0004

0,2
0,0002

0,0
0,0000
0 10 20 30 40 50 60 70 80 0 10 20 30 40 50 60 70 80
Fois Fois

Fig. 8 Flux de chaleur (a) généré par EDEM et flux de chaleur (b) généré par fluent à 16 tr/min.

9
 100

90
EDEM-CFD
80
ÉDEM

70

Température/°C
60

50

40

30

20

0 20 40 60 80
Fois

Fig.9 Changement de la température des particules de tabac en fonction du temps

Fig. 10 Image en nuage du changement de température des particules de tabac

3.5. Distribution du temps de séjour des particules

Fig.11 Temps de séjour des particules

10
 Comme le montre la figure 11, la distribution du temps de séjour (RTD) à 16 tr/min a une grande influence sur la

qualité du séchage des particules. Si le temps de séjour est trop court, il entraînera un chauffage inégal des matériaux, si

le temps de séjour du matériau est trop long, il entraînera une augmentation du degré d'endommagement des particules

de matériau et le rendement correspondant sera réduit. La figure 11 montre la distribution du temps de séjour des

particules au cours des 5 premières secondes, le temps de distribution le plus long étant principalement concentré près

de la paroi, ce qui est plus cohérent avec la réalité du séchage du tabac.

4. Conclusions

Français Dans cet article, la méthode de couplage EDEM-FLUENT a été utilisée pour calculer les caractéristiques de

transfert de chaleur et de masse du processus de séchage du tabac dans un séchoir à tambour, les caractéristiques de

transfert de chaleur et de masse du tabac coupé à différentes vitesses de rotation ont été étudiées. Les résultats

montrent que : (1) avec l'augmentation de la vitesse de rotation du séchoir à tambour, le degré de mélange augmente

progressivement. Lorsque la vitesse de rotation est de 16 tr/min, le nombre de contacts entre les particules de tabac et la

paroi du tambour et la feuille du tambour est égal au nombre de contacts entre les particules de tabac, c'est-à-dire un

degré de mélange plus élevé. Le temps de séchage des particules de tabac dans le tambour diminue linéairement avec

l'augmentation de la vitesse de rotation du tambour. Lorsque la vitesse de rotation est de 16 tr/min, la vitesse axiale

moyenne des particules est d'environ 1,75 m/s, l'efficacité de séchage peut être améliorée en augmentant la vitesse de

rotation du tambour. (2) Lorsque le flux de chaleur du couplage EDEM-CFD est compris entre 0,2 W et 1 W, la temp

érature du grain de tabac augmentera progressivement, mais le taux de changement de température diminuera

progressivement et la température du grain de tabac aura tendance à être constante. (3) Le temps de séjour du matériau

a une grande influence sur la qualité de séchage des particules de matériau. Si le temps de rétention du tabac coupé est

trop court, cela entraînera un chauffage inégal du tabac coupé, si le temps de rétention du tabac coupé est trop long, cela

entraînera une augmentation de la sécheresse du tabac coupé.

Remerciements

Les auteurs tiennent à remercier le projet des talents d'innovation scientifique et technologique de

Zhongyuan (234200510011) et le programme d'innovation collaborative de Zhengzhou (majeur et

Programme clé de ZZULI) (2021ZDPY0107)，Projet clé de science et de technologie de China Tobacco

Hebei Co., LTD.(HBZY2023A056) et projet scientifique et technologique clé de China Tobacco

Henan Co., LTD.( AN202130).

Référence

[1] Rezaei, H., S. Sokhansanj, Une étude sur la détermination du temps de séjour des particules solides dans les sécheurs à

tambour rotatif, Technologie de séchage,39. (2021), 11, pp. 1762-1772

[2] Balakrishnan, M., et al.Conception, développement et évaluation d'un séchoir à tambour rotatif pour rhizomes de

curcuma (Curcuma longa L.)Journal de l'ingénierie des procédés alimentaires,45. (2022), 6, pp. 1-16

[3] Mikulionok, IO, Classification des sécheurs à tambour à convection (survol des brevets), Ingénierie
chimique et pétrolière,56. (2020), 7, pp. 588-596
[4] Zhu, W. K., et al.Effet de la déshydratation en deux étapes sur la rétention des composants aromatiques caractéristiques du

tabac séché à l'air chaud dans un séchoir rotatif, Technologie de séchage,34. (2016), 13, pp. 1621-1629

11
[5] Lee, J., et al.Effet des caractéristiques du tissu et des vitesses de rotation du tambour sur les mouvements et les performances

de séchage des vêtements dans un sèche-linge à tambour, Fibres et polymères,23. (2022), 6, pp. 1741- 1752

[6] Namwong, K., et al.Analyse thermique du séchage à l'air chaud des granulés de fumier de poulet dans un séchoir

rotatif horizontal portable modifiéTransfert de chaleur et de masse,58. (2022), 11, pp. 1997-2007

[7] Gui, N., J. Fan, Étude numérique de la conduction thermique des particules granulaires dans des tambours ondulés rotatifs

Journal international des transferts de chaleur et de masse,84. (2015), pp. 740-751

[8] Yang, G., et al.Etude numérique sur la distribution uniforme du champ d'écoulement d'un sécheur à flux d'air

, Hélios,10. (2024), 8, p. e29439

[9] Nascimento, SM, et al.Comparaison entre les approches eulérienne (CFD) et lagrangienne (DEM)
dans la simulation d'un tambour rotatif à volutes, Mécanique numérique des particules,9. (2022), 2,
pp. 251-263
[10] Zhong, W., et al., Modélisation DEM/CFD-DEM des systèmes particulaires non sphériques : d
éveloppements théoriques et applications, Technologie des poudres,302. (2016), pp. 108-152
[11] Ma, H., et al., Une revue des développements récents pour les études CFD-DEM des particules
non sphériques, Technologie des poudres,412. (2022), p. 117972
[12] Quan, X., et al.((2024)Modélisation CFD-DEM de l'écoulement d'hydrogène cryogénique et du transfert de chaleur

dans un lit garni. Durabilité 16

[13] Chen, H., Z. Zheng, Empilement binaire de particules sphériques avec des rapports de taille modérés dans

un fluide visqueux : une étude CFD-DEMJournal brésilien de physique,54. (2024), 4, p. 110

[14] Zhang, X., P. Newell, Simulation du transport de microcapsules dans des milieux fracturés à l'aide de la

CFD-DEM coupléeRecherche sur les ressources en eau,60. (2024), 6, p. e2023WR035728

[15] Scherer, V., et al.Simulation couplée DEM‒CFD du séchage des copeaux de bois dans un tambour rotatif ‒ Conception de d

éflecteurs et réduction du modèle, Carburant,184. (2016), pp. 896-904

[16] Wu, K., et al.Simulation tridimensionnelle de l'écoulement gaz-solide dans un lit fluidisé avec des particules en

ruban flexibleJournal international des écoulements multiphasiques,124. (2020), p. 103181

[17] Markauskas, D., et al.Etude numérique comparative du transport pneumatique de particules allongées
flexibles à travers un coude de tuyau par DEM-CFD, Technologie des poudres,399. (2022), p. 117170

[18] Chen, L., et al.Étude de la collision entre une particule grossière unique et une paroi dans des milieux visqueux à

l'aide de la CFD‒DEM, Journal Européen de Mécanique - B/Fluides,97. (2023), pp. 199-207

[19] Tao, W.,transfert de chaleur (cinquième édition)Presses de l'enseignement supérieur, Pékin, 2019.

[20] Hou, ZJ, Étude sur les caractéristiques de transfert de chaleur du four rotatif électrothermique à alumine,

Université du Shandong, 2021.

[21] Chaudhuri, B., et al., Modélisation du transfert de chaleur en écoulement granulaire dans des cuves rotatives,

Sciences du génie chimique,61. (2006), 19, pp. 6348-6360

[22] Gupta, A., et al. Étude des éléments discrets Mélange dans un mélangeur de taille industrielle，Congrès
mondial des technologies des particules 2010, 'Procédure' (INM GmbH), 2010, p. 1-4

Soumis le 19.4.2024.
Révisé 25.6.2024
Accepté 27.6.2024.

12
13