Page 1
Simulation de réacteurs multicanaux compacts : utilisation de zones
poreuses équivalentes
Ghislain GENIN
Equipe PRISM (dirigée par Jean Marc COMMENGE) Procédés Photoniques, Réactifs, Intensifiés et Systèmes Microstructurés
Laboratoire Réactions et Génie des Procédés
Journées Promotion Procédés Produits NANCY – ENSIC – mardi 25 novembre 2014
Utilisation de la CFD en génie des procédés Gestion de la complexité des systèmes
Page 2
Intensification, miniaturisation des procédés
Réacteurs – Echangeurs multicanaux
Co-courant
Contre-courant
Courants croisés
Difficultés des simulations
Nombre de canaux élevé
Facteurs de formes des canaux importants (L/D)
Nombre de faces / volumes (dessin + maillage)
Besoins en simulation fine 3D
Cartes de température / pression / concentration
Résistance mécanique et chimique
Limiter le nombre d’essais sacrificiels (coûts et délais de fabrication)
Visualiser les maldistributions hydrodynamique / thermique
Equivalence
hydraulique Introduction
Equivalence
thermique
Equivalence
réactionnelle
Conclusion
Perspectives
Page 3
Réduction géométrique Seule une couche élémentaire est considérée
(non visualisation de la maldistribution entre couches)
Nouvelle simplification :
utilisation de milieux poreux équivalents
Série de canaux et parois assimilées à un milieu poreux
… … … …
Diminution du nombre de faces, de volumes, de mailles
Comment définir ce poreux pour qu’il soit équivalent ? Perte de charge
Transfert thermique
Vitesse globale de réaction
Equivalence
hydraulique Introduction
Equivalence
thermique
Equivalence
réactionnelle
Conclusion
Perspectives
Page 4
Résistance équivalente du milieu poreux
Simulation sur un canal Fluide = air
z x
y
Entrée umoy= 1 m/s (Re = 68,5)
Sortie
plans xOz yOz = symétries 1 m
m
Profil de vitesse dans le canal Profil de vitesse (interstitielle) dans le poreux (e = 0,5)
Entrée umoy= 0,5 m/s
e = 0,5
Equivalence
hydraulique Introduction
Equivalence
thermique
Equivalence
réactionnelle
Conclusion
Perspectives
Profil de vitesse établi
Pression nulle
Longueur de 2 cm Mailles hexaédrique de 0,05 mm 160 000 mailles
Page 5
Valeur théorique de la résistance : loi de Poiseuille (dans le cas d’un profil de vitesse uniforme)
eDmR
H
nc 89,0 10.848,2 27
1
eDmR
H
nc
2
5,1 10.2,1 27
2
Résistance du milieu poreux
Résistance supplémentaire
imposée au milieu poreux pour
modéliser la présence des
canaux
21 R
RRP
e 0
2
4
6
8
10
12
0 0,005 0,01 0,015 0,02
longeur z (m)
DP (Pa)
Canal
Poreuxéquivalent ?
DP surestimée : influence du profil de vitesse…
e = 1 m
m
e = 1 mm
uc
e =
1 m
m
up
Résistance équivalente du milieu poreux
Equivalence
hydraulique Introduction
Equivalence
thermique
Equivalence
réactionnelle
Conclusion
Perspectives
uDL
P
H
nc
32
2
D
R
Page 6
Réglage fin de la résistance du milieu poreux
L’augmentation de la
résistance déforme le
profil de vitesse
RP obtenue par itérations pour obtenir DPP = DPC
Résistances sur les 3 axes du milieu poreux
RZ = RP RY = 0
RX = ? RP
Valeur RX élevée pour
modéliser les parois des
canaux
Résistance équivalente du milieu poreux
Equivalence
hydraulique Introduction
Equivalence
thermique
Equivalence
réactionnelle
Conclusion
Perspectives
Page 7
Un exemple : échangeur air-air à courants croisés
10 canaux de section carrée de côté 1 mm séparés par des parois de 1 mm
z
y
x
multicanaux poreux
Maillage structuré hexaèdre de 0,1 mm de côté => 882 000 mailles
Equivalence
hydraulique Introduction
Equivalence
thermique
Equivalence
réactionnelle
Conclusion
Perspectives
e = 10/21 RX = RP
e = 10/21 RZ = RP
Page 8
Echangeur air-air à courants croisés : lignes de courant
Equivalence
hydraulique Introduction
Equivalence
thermique
Equivalence
réactionnelle
Conclusion
Perspectives
RX = 100 RZ
trop faible (dispersion transversale)
RX =105 RZ
« mauvais » distributeur
Plan de symétrie xOz (plan milieu des canaux froids)
canaux poreux
Page 9
Echangeur air-air à courants croisés : pression statique
Equivalence
hydraulique Introduction
Equivalence
thermique
Equivalence
réactionnelle
Conclusion
Perspectives
Cartes de pression statique semblables
maldistribution équivalente ?
Plan de symétrie xOz (plan milieu des canaux froids)
canaux poreux
Page 10
Echangeur air-air à courants croisés : débits massiques
Sorties 1 à 10
Entrées 1 à 10
« Entrées poreux » 1 à 10
« Sorties poreux » 1 à 10
Equivalence
hydraulique Introduction
Equivalence
thermique
Equivalence
réactionnelle
Conclusion
Perspectives
« sorties » 1 et 10 surdimensionnées (+25%)
1 2
1 2
Page 11
Echangeur air-air à courants croisés : débits massiques
Sorties 1 à 10
Entrées 1 à 10
« Entrées poreux » 1 à 10
« Sorties poreux » 1 à 10
Equivalence
hydraulique Introduction
Equivalence
thermique
Equivalence
réactionnelle
Conclusion
Perspectives
1 2
1 2
Symétrie non respectée ? profil de vitesse à l’entrée à affiner Maldistribution correctement modélisée (écarts < 10%)
Page 12
Simulation sur un canal
z x
y
plans xOz yOz = symétries
T 500 K
1 m
m
Fluide = air f = 0,0242 Wm-1K-1 Solide = aluminium s = 202,4 Wm-1K-1
Evolution de Tmoy (pondérée par le débit massique) suivant Oz
Recherche de la conductibilité thermique équivalente (P) permettant d’obtenir
un profil de température identique (minimisation de la somme des carrés des
écarts)
T 300 K
Conductibilité thermique équivalente du milieu poreux
Equivalence
hydraulique Introduction
Equivalence
thermique
Equivalence
réactionnelle
Conclusion
Perspectives
Page 13
Profil de température dans le canal Profil de température dans le poreux
Zone poreuse Fluent
( milieu // )
p = e f + (1-e) s
Solide poreux ≠ aluminium
Conductibilité thermique équivalente du milieu poreux
Equivalence
hydraulique Introduction
Equivalence
thermique
Equivalence
réactionnelle
Conclusion
Perspectives
Page 14
p = 0,847 f
p < p < p
Conductibilité thermique équivalente du milieu
poreux
Equivalence
hydraulique Introduction
Equivalence
réactionnelle
Conclusion
Perspectives
Profil de température dans le poreux
Equivalence
thermique
Page 15
z x
y
Poreux orthotrope ?
Profil de température dans le poreux
PY = 0,870 f
PZ = 5,035 f
Recherche des 3 i complexe risque de rendre le modèle simplifié trop compliqué à mettre en œuvre…
Conductibilité thermique équivalente du milieu
poreux
Equivalence
hydraulique Introduction
Equivalence
réactionnelle
Conclusion
Perspectives
Equivalence
thermique
Page 16
Echangeur air-air à courants croisés : profils de température
Equivalence
hydraulique Introduction
Equivalence
réactionnelle
Conclusion
Perspectives
Equivalence
thermique
canaux
poreux
Plan de symétrie xOz (plan milieu des canaux froids)
Page 17
canaux
poreux
Echangeur air-air à courants croisés : profils de température
Equivalence
hydraulique Introduction
Equivalence
réactionnelle
Conclusion
Perspectives
Equivalence
thermique
Plan de symétrie xOz (plan milieu des canaux chauds)
Cartes de température
semblables
maldistribution thermique
équivalente ?
Page 18
Equivalence
hydraulique Introduction
Equivalence
réactionnelle
Conclusion
Perspectives
Equivalence
thermique
Echangeur air-air à courants croisés : température moyenne en entrées / sorties
Plan de symétrie xOz (plan milieu des canaux
froids)
E1
S1
E10
S10
E1
S1 S10
E10
Cohérence avec les débits et le
fonctionnement contre-courant.
Maldistribution en sortie bien modélisée.
Page 19
Equivalence
hydraulique Introduction
Equivalence
réactionnelle
Conclusion
Perspectives
Equivalence
thermique
Echangeur air-air à courants croisés : température d’entrées des canaux froids ?
S2 S5
S2 S5
Poreux : température le long des « canaux » 2 et 5
Canaux : température le long des canaux 2 et 5
Effets de parois (cf. mono canal)
Page 20
Poreux : profil de vitesse selon Z à l’entrée des « canaux »
Equivalence
hydraulique Introduction
Equivalence
réactionnelle
Conclusion
Perspectives
Equivalence
thermique
Echangeur air-air à courants croisés : température d’entrées des canaux froids ?
Canaux : profil de vitesse selon Z à l’entrée des canaux
S5
S5
Effet de retro mélange plus
marqué à l’entrée du
milieu poreux ( RX trop
élevée)
Effet dû au mauvais distributeur.
Page 21
canaux
Echangeur air-air à courants croisés : profils de température
Equivalence
hydraulique Introduction
Equivalence
réactionnelle
Conclusion
Perspectives
Equivalence
thermique
Plan xOz au milieu de la paroi métallique
poreux
Cartes de température
semblables
DT surestimée (20%)
Page 22
… … Correction des
constantes cinétiques
suivant le rapport des
surfaces
… …
Réaction hétérogène en paroi
Réaction homogène
… … Correction des
constantes cinétiques
suivant la valeur de e
… …
Correction surfacique ou volumique
Equivalence
hydraulique Introduction
Equivalence
réactionnelle
Conclusion
Perspectives
Equivalence
thermique
Page 23
Estimation du gain en nombre de mailles et temps de calcul (vers une simulation d’un module complet ?)
Développement d’un outil de calcul pour estimer Rp et eq en fonction de : • la géométrie des canaux • la nature des matériaux • l’épaisseur des parois
Canaux obstrués : modélisation par un profil géométrique de RP
Méthode utilisée sur un réacteur catalytique de conversion de méthane en hydrogène
Méthode facilitant l’étude de l’effet de la forme et du nombre de canaux : géométrie « poreuse » identique / Rp et eq différents
Méthode applicable pour des études hydrodynamiques sur des procédés avec des forts changements d’échelles (réseaux multitubulaires)
Equivalence
hydraulique Introduction
Equivalence
réactionnelle
Conclusion
Perspectives
Equivalence
thermique
Page 24
Simulation de réacteurs multicanaux compacts : utilisation de zones
poreuses équivalentes
Journées Promotion Procédés Produits NANCY – ENSIC – mardi 25 novembre 2014
Utilisation de la CFD en génie des procédés Gestion de la complexité des systèmes
Merci.
[email protected]