Dans Quelles Conditions L’hypothèse Du Continuum Fluide Est-Elle Valide Pour Modéliser La Dynamique Des Fluides Dans Les Installations Pilotes D’opérations Unitaires ?

L’hypothèse du continuum fluide est valide lorsque le libre parcours moyen moléculaire est très inférieur à la longueur caractéristique de l’écoulement. Cette condition est quantifiée par le nombre de Knudsen, (\mathrm{Kn} = \ell_0/L), où (\ell_0) est le libre parcours moyen et (L) la plus petite dimension pertinente de l’équipement ou du champ d’écoulement. Pour les installations pilotes conventionnelles d’opérations unitaires, (\mathrm{Kn} \ll 1), de sorte que la densité, la vitesse, la température et la pression peuvent être considérées comme des champs variant de manière continue et être modélisées à l’aide des équations de bilan en milieu continu.

À retenir : La modélisation en milieu continu est généralement appropriée lorsque (\mathrm{Kn} < 0.01), et particulièrement fiable lorsque cette valeur est très inférieure à ce seuil. L’hypothèse peut échouer dans les gaz fortement raréfiés, les équipements sous vide poussé, les passages de taille micro- ou nanométrique, ou à proximité de régions dont les dimensions approchent le libre parcours moyen moléculaire.

Ce que requiert l’hypothèse du continuum

L’équipement doit être grand par rapport aux échelles moléculaires

La longueur caractéristique (L) doit être choisie comme la plus petite dimension qui contrôle la physique locale de l’écoulement. Selon le problème, il peut s’agir du diamètre d’une conduite, de l’épaisseur d’une couche limite, de la taille des pores, de la hauteur d’un canal ou de la distance sur laquelle la température et la vitesse varient de manière significative.

L’exigence essentielle est la suivante :

[ \mathrm{Kn}=\frac{\ell_0}{L}\ll 1 ]

Lorsque ce rapport est faible, un volume local représentatif contient suffisamment de molécules pour que les fluctuations moléculaires se compensent par moyennage. Les propriétés macroscopiques peuvent alors être attribuées à chaque point du champ d’écoulement.

Les collisions moléculaires doivent établir un équilibre local

L’approximation du continuum suppose que les molécules interagissent suffisamment fréquemment pour que les propriétés thermodynamiques locales, telles que la densité, la température et la pression, aient une signification. Dans les gaz, cela exige que la dimension de l’écoulement soit très supérieure au libre parcours moyen moléculaire.

Pour les liquides dans des conditions ordinaires, l’espacement moléculaire est extrêmement faible et les interactions intermoléculaires sont fréquentes. Par conséquent, le comportement continu est généralement valide dans les équipements conventionnels des installations pilotes, sauf si le procédé implique des géométries exceptionnellement petites ou des conditions physiques inhabituelles.

La pression de fonctionnement ne doit pas rendre le gaz fortement raréfié

Pour les gaz, le libre parcours moyen augmente lorsque la pression diminue. Une installation pilote qui se trouve nettement dans le régime continu à la pression atmosphérique peut se rapprocher du régime d’écoulement avec glissement ou du régime de transition sous vide poussé.

La pression doit donc être évaluée conjointement avec l’échelle de l’équipement. Le critère pertinent est la valeur locale de (\mathrm{Kn}), en particulier dans les régions à basse pression, les étranglements étroits, les pores et à proximité des interfaces sous vide.

Comment interpréter le nombre de Knudsen

Régime continu

Pour environ

[ \mathrm{Kn}<0.01 ]

les équations classiques du milieu continu et la condition aux limites de non-glissement sont généralement appropriées. Il s’agit du régime normal pour les conduites, les réservoirs, les colonnes, les échangeurs thermiques et de nombreux systèmes de mélange fonctionnant à des pressions ordinaires.

Régime d’écoulement avec glissement

Pour environ

[ 0.01\lesssim\mathrm{Kn}<0.1 ]

le fluide en volume peut encore être modélisé avec les équations du milieu continu, mais les effets de paroi peuvent nécessiter des conditions aux limites de glissement ou de saut de température. La modélisation standard de Navier-Stokes avec condition de non-glissement peut alors produire des erreurs mesurables.

Régimes transitoire et moléculaire libre

Lorsque (\mathrm{Kn}) approche ou dépasse (0.1), les équations du milieu continu deviennent progressivement moins fiables. Pour des valeurs plus élevées, des méthodes moléculaires telles que la méthode de simulation directe de Monte-Carlo ou des approches fondées sur la théorie cinétique peuvent être nécessaires.

Ces régimes sont peu courants dans les opérations unitaires ordinaires des installations pilotes, mais peuvent apparaître lors de procédés sous vide, de fuites de gaz, dans de petits pores et dans des dispositifs à l’échelle microscopique.

Quelles équations deviennent valides

Les cinq équations de bilan macroscopiques

Lorsque l’hypothèse du continuum et les hypothèses constitutives appropriées sont satisfaites, l’écoulement peut être décrit par :

Une équation de conservation de la masse
Trois équations de conservation de la quantité de mouvement
Une équation de conservation de l’énergie

Ces équations déterminent les champs macroscopiques après avoir été complétées par des relations constitutives, des conditions aux limites et des données sur les propriétés thermodynamiques.

Les propriétés de transport doivent être bien définies

Le modèle nécessite des valeurs locales de propriétés telles que :

[ \rho=\rho(T,p), \qquad c_p=c_p(T,p), \qquad \mu=\mu(T,p), \qquad \lambda=\lambda(T,p) ]

Ici, (\rho) est la densité, (c_p) la capacité thermique isobare, (\mu) la viscosité dynamique et (\lambda) la conductivité thermique.

Pour un fluide monophasique proche de l’équilibre thermodynamique local, ces propriétés peuvent généralement être obtenues à partir de corrélations de propriétés à l’équilibre ou d’une équation d’état appropriée.

Les lois constitutives doivent correspondre au fluide

Les équations de Navier-Stokes supposent en outre une relation constitutive entre les contraintes et la déformation. Pour un fluide newtonien isotrope, la contrainte visqueuse est linéairement liée au taux de déformation.

La loi de Fourier est également utilisée pour la conduction thermique :

[ \mathbf{q}=-\lambda\nabla T ]

Ces hypothèses conviennent généralement aux gaz ordinaires et aux liquides simples, mais peuvent nécessiter des modifications pour les fluides non newtoniens, les matériaux anisotropes, les fluides viscoélastiques ou les écoulements présentant de forts effets de non-équilibre.

Conditions au-delà du comportement continu

L’équilibre thermodynamique local doit s’appliquer

Le principe de l’état local suppose que chaque élément de fluide suffisamment petit peut être décrit à l’aide de relations d’équilibre, même si les propriétés varient d’un point à l’autre.

Cette hypothèse devient discutable lorsque les gradients sont extrêmement raides, que la relaxation moléculaire est lente par rapport à l’échelle de temps de l’écoulement ou que le gaz est fortement raréfié.

Le modèle de phase et de composition doit être approprié

Un fluide monophasique, monocomposant ou un mélange bien caractérisé constitue le cadre le plus simple pour la modélisation en milieu continu. Les écoulements multiphasiques, les systèmes réactionnels et les écoulements présentant une diffusion massique importante nécessitent des équations de bilan et des modèles constitutifs supplémentaires.

L’hypothèse du continuum elle-même n’interdit pas la modélisation multiphasique, mais chaque phase doit être représentée de manière appropriée et le comportement interphasique doit être résolu ou modélisé.

La viscosité et la conductivité thermique doivent rester physiquement admissibles

Pour respecter une production d’entropie non négative, la viscosité dynamique et la conductivité thermique doivent être positives. Leurs valeurs doivent également rester appropriées sur la plage locale de température, de pression et de composition.

Les termes de force volumique et de source de chaleur doivent être inclus lorsqu’ils sont pertinents

Négliger la génération volumique de chaleur n’est acceptable que lorsque l’absorption du rayonnement, la réaction chimique, l’échauffement par effet Joule ou d’autres sources sont négligeables. Si ces effets sont importants, ils doivent apparaître explicitement dans l’équation de l’énergie.

De même, la gravité, les forces électromagnétiques et les autres forces volumiques doivent être incluses lorsqu’elles influencent de manière significative l’opération.

Comprendre les compromis

Un faible nombre de Knudsen ne garantit pas une simulation précise

La validité du milieu continu indique uniquement si une description par champs macroscopiques est appropriée. Elle ne garantit pas la précision du modèle de turbulence, de la cinétique réactionnelle, de la fermeture multiphasique, de la condition aux parois, du maillage ou de la corrélation de propriétés.

Un résultat de CFD en milieu continu peut donc être mathématiquement cohérent tout en échouant à prédire correctement la perte de charge, le mélange, le transfert thermique ou la séparation.

La longueur caractéristique peut varier au sein d’une même unité

Utiliser uniquement le diamètre de l’équipement peut masquer une rupture locale du comportement continu. Les petits pores d’un lit garni, les films minces, les étranglements, les couches limites et les chemins de fuite peuvent avoir des dimensions caractéristiques bien plus petites que celles de la cuve globale.

La plus petite échelle locale physiquement pertinente doit être évaluée, en particulier lors d’un fonctionnement à basse pression.

Les conditions de non-glissement peuvent échouer avant le modèle continu

Pour des nombres de Knudsen modérés, le fluide peut encore obéir aux équations continues de la quantité de mouvement dans le volume, tandis que la condition classique de non-glissement à la paroi devient inexacte. C’est pourquoi la décision concernant le caractère continu ou non continu et celle concernant les conditions aux limites doivent être traitées séparément.

La validation à l’échelle pilote reste importante

Les mesures des profils de vitesse, des pertes de charge, des taux de transfert thermique et des zones de séparation peuvent vérifier si les modèles continus et constitutifs choisis représentent correctement l’opération unitaire réelle.

La concordance avec des cas de référence tels que l’écoulement laminaire pleinement développé dans une conduite est utile, mais ne remplace pas la validation pour des géométries ou des conditions de fonctionnement complexes.

Comment appliquer cela à votre installation pilote

Les vérifications suivantes constituent une procédure pratique d’évaluation préliminaire :

Si votre objectif principal concerne les opérations unitaires conventionnelles sur des liquides ou des gaz : calculez le (\mathrm{Kn}) local en utilisant la plus petite dimension d’écoulement pertinente ; lorsque (\mathrm{Kn}<0.01), les équations standard de Navier-Stokes en milieu continu et de l’énergie sont généralement justifiées.
Si votre objectif principal concerne les équipements sous vide, à l’échelle microscopique ou à milieu poreux : recalculez (\mathrm{Kn}) aux pressions les plus faibles et dans les passages les plus étroits, et envisagez des modèles d’écoulement avec glissement ou des modèles cinétiques si (\mathrm{Kn}) n’est pas très inférieur à (0.01).
Si votre objectif principal concerne les transferts de chaleur et de matière : vérifiez l’équilibre local, la conduction de Fourier, la diffusion et les corrélations de propriétés en plus du contrôle de (\mathrm{Kn}).
Si votre objectif principal concerne la prédiction CFD : considérez la validité du milieu continu comme la première étape d’évaluation seulement, puis validez les lois constitutives, les conditions aux limites, la qualité du maillage, les modèles de fermeture et les résultats par rapport aux mesures réalisées sur l’installation pilote.

Pour les installations pilotes ordinaires fonctionnant avec des dimensions supérieures à l’échelle microscopique et à des pressions non raréfiées, un nombre de Knudsen très faible constitue normalement une base fiable pour la modélisation macroscopique de la dynamique des fluides.

Tableau récapitulatif :

Régime d’écoulement	Nombre de Knudsen	Approche recommandée	Applications typiques
Continu	Kn < 0.01	Les équations standard de Navier-Stokes et les frontières sans glissement sont généralement appropriées	Conduites, réservoirs, colonnes, échangeurs thermiques, mélangeurs
Écoulement avec glissement	0.01 ≤ Kn < 0.1	Utiliser les équations du milieu continu avec des conditions aux parois de glissement ou de saut de température	Systèmes à basse pression et passages étroits
Transitoire	0.1 ≤ Kn < 1	Les modèles continus deviennent moins fiables ; envisager des méthodes cinétiques	Équipements sous vide, petits pores, canaux à l’échelle microscopique
Moléculaire libre	Kn ≥ 1	Utiliser une modélisation moléculaire ou cinétique, telle que la méthode DSMC	Gaz fortement raréfiés et systèmes sous vide poussé

Vérification essentielle : évaluez le nombre de Knudsen local en utilisant la plus petite dimension d’écoulement pertinente et la pression de fonctionnement la plus faible.

Vous avez besoin d’équipements fiables pour installations pilotes en dynamique des fluides, transfert de chaleur et de matière ou validation de procédés ? Contactez LABPARK pour discuter de vos besoins. LABPARK fournit des installations pilotes d’opérations unitaires à vocation éducative et professionnelle pour le génie chimique, les bioprocédés et les biotechnologies, ainsi que pour les applications environnementales et le traitement de l’eau. Nos systèmes aident les universités, les instituts de recherche et les entreprises à démontrer les principes des écoulements continus, à former les ingénieurs, à valider les modèles de procédés et à produire des données pratiques d’exploitation grâce à des équipements et à un accompagnement axés sur les applications.