Les équations fondamentales sont l'équation de continuité pour un fluide incompressible et les équations de quantité de mouvement de Navier-Stokes. Pour un fluide newtonien, isotherme et incompressible, de masse volumique (\rho) et de viscosité dynamique (\mu) constantes, elles sont :

[ \nabla \cdot \mathbf{v}=0 ]

[ \rho\left(\frac{\partial \mathbf{v}}{\partial t} +\mathbf{v}\cdot\nabla\mathbf{v}\right)

-\nabla p+\mu\nabla^2\mathbf{v}+\rho\mathbf{g}. ]

Ensemble, ces équations déterminent les trois composantes de la vitesse et le champ de pression, sous réserve des conditions appropriées de géométrie de l'équipement, d'entrée, de sortie, de paroi et initiales.

Point clé : L'écoulement newtonien, isotherme et incompressible dans les équipements d'usine pilote est modélisé principalement par la conservation de la masse et la conservation de la quantité de mouvement. Les équations relatives à l'énergie, aux espèces, aux réactions ou aux systèmes multiphasiques ne sont ajoutées que lorsque l'opération unitaire inclut des effets qui ne peuvent être représentés par ces deux seules équations.

L'ensemble des équations directrices

Conservation de la masse : l'équation de continuité

L'équation générale de conservation de la masse est :

[ \frac{\partial \rho}{\partial t}+\nabla\cdot(\rho\mathbf{v})=0. ]

Pour un fluide incompressible de masse volumique constante :

[ \boxed{\nabla\cdot\mathbf{v}=0} ]

Cette condition stipule qu'un élément de fluide ne se dilate pas et ne se contracte pas lors de son déplacement. En coordonnées cartésiennes :

[ \frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} + \frac{\partial w}{\partial z} =0, ]

où (u), (v) et (w) sont les composantes de la vitesse dans les directions (x), (y) et (z).

Conservation de la quantité de mouvement : les équations de Navier–Stokes

Le bilan de quantité de mouvement linéaire est :

[ \boxed{ \rho\frac{D\mathbf{v}}{Dt}

-\nabla p+\mu\nabla^2\mathbf{v}+\rho\mathbf{g} } ]

où la dérivée particulaire est :

[ \frac{D\mathbf{v}}{Dt}

\frac{\partial\mathbf{v}}{\partial t} + (\mathbf{v}\cdot\nabla)\mathbf{v}. ]

Ainsi, l'équation de quantité de mouvement représente un équilibre entre :

L'accélération du fluide : (\rho D\mathbf{v}/Dt)
Les forces de pression : (-\nabla p)
Les forces visqueuses : (\mu\nabla^2\mathbf{v})
Les forces massiques : (\rho\mathbf{g}), généralement la gravité

Le terme convectif ((\mathbf{v}\cdot\nabla)\mathbf{v}) est important dans les coudes, les contractions, les expansions, les mélangeurs et autres zones où la vitesse varie spatialement.

Forme des composantes cartésiennes

Pour une vitesse (\mathbf{v}=(u,v,w)), les équations de quantité de mouvement sont :

[ \rho\left( \frac{\partial u}{\partial t} +u\frac{\partial u}{\partial x} +v\frac{\partial u}{\partial y} +w\frac{\partial u}{\partial z} \right)

-\frac{\partial p}{\partial x} +\mu\nabla^2u+\rho g_x, ]

[ \rho\left( \frac{\partial v}{\partial t} +u\frac{\partial v}{\partial x} +v\frac{\partial v}{\partial y} +w\frac{\partial v}{\partial z} \right)

-\frac{\partial p}{\partial y} +\mu\nabla^2v+\rho g_y, ]

[ \rho\left( \frac{\partial w}{\partial t} +u\frac{\partial w}{\partial x} +v\frac{\partial w}{\partial y} +w\frac{\partial w}{\partial z} \right)

-\frac{\partial p}{\partial z} +\mu\nabla^2w+\rho g_z. ]

Ces trois équations, combinées à la continuité, fournissent quatre équations pour quatre champs principaux : (u), (v), (w) et (p).

Pourquoi le comportement newtonien est important

Loi de viscosité de Newton

Un fluide newtonien présente des contraintes visqueuses qui varient linéairement avec le taux de déformation. Sous forme tensorielle :

[ \boldsymbol{\tau}

2\mu\mathbf{D}, ]

pour un écoulement incompressible, où

[ \mathbf{D}

\frac{1}{2} \left[ \nabla\mathbf{v} + (\nabla\mathbf{v})^{T} \right] ]

est le tenseur du taux de déformation.

Pour une viscosité constante et un écoulement incompressible, la divergence de la contrainte visqueuse se réduit à :

[ \nabla\cdot\boldsymbol{\tau}

\mu\nabla^2\mathbf{v}. ]

Cette réduction est ce qui produit le terme visqueux dans l'équation standard de Navier-Stokes.

Interprétation physique

Le terme visqueux décrit la diffusion de la quantité de mouvement. Il rend compte du transfert de quantité de mouvement causé par les gradients de vitesse, tels que ceux existant entre une paroi de tuyau stationnaire et le fluide en mouvement au centre du tuyau.

Conditions nécessaires pour fermer le modèle

Condition de non-glissement à la paroi

Au niveau d'une paroi solide stationnaire :

[ \mathbf{v}=0. ]

Cette condition est essentielle pour prédire la contrainte de cisaillement à la paroi, la perte de charge et les profils de vitesse dans les tuyaux, les canaux, les échangeurs de chaleur et les équipements garnis.

Pour une paroi mobile, la vitesse du fluide à la paroi est égale à la vitesse de la paroi.

Conditions d'entrée

Une entrée peut être spécifiée par un profil de vitesse,

[ \mathbf{v}=\mathbf{v}_{\text{in}}, ]

ou par un débit massique,

[ \dot{m}

\int_A \rho,\mathbf{v}\cdot\mathbf{n},dA. ]

Un profil d'entrée uniforme est souvent utilisé pour des études simplifiées, tandis que des profils développés ou mesurés sont plus appropriés pour une analyse détaillée en usine pilote.

Conditions de sortie

Les spécifications de sortie typiques incluent une pression connue,

[ p=p_{\text{out}}, ]

ou une condition de sortie appropriée lorsque l'écoulement en aval est suffisamment développé. La condition de sortie sélectionnée peut affecter de manière significative les champs de pression prédits et la recirculation.

Condition initiale

Pour un fonctionnement en régime transitoire, le champ de vitesse initial doit être défini :

[ \mathbf{v}(\mathbf{x},0)=\mathbf{v}_0(\mathbf{x}). ]

Une référence de pression est également requise car la pression dans un fluide incompressible n'est généralement déterminée qu'à une constante arbitraire près.

Simplifications courantes dans les modules d'usine pilote

Écoulement permanent

Si les conditions de fonctionnement ne varient pas de manière significative avec le temps :

[ \frac{\partial\mathbf{v}}{\partial t}=0. ]

L'équation de quantité de mouvement devient :

[ \rho(\mathbf{v}\cdot\nabla)\mathbf{v}

-\nabla p+\mu\nabla^2\mathbf{v}+\rho\mathbf{g}. ]

Cette approximation est couramment utilisée pour les analyses de tuyaux, de canaux et d'échangeurs de chaleur en régime permanent.

Écoulement interne pleinement développé

Pour un écoulement pleinement développé dans un canal ou un tuyau droit, le profil de vitesse ne change pas dans la direction axiale. Les équations directrices peuvent alors se réduire à un équilibre entre le gradient de pression et la résistance visqueuse.

Par exemple, pour un écoulement unidimensionnel entre des plaques parallèles :

[ 0=-\frac{dp}{dx}+\mu\frac{d^2u}{dy^2}. ]

Le profil de vitesse résultant est parabolique lorsque le gradient de pression est constant.

Inertie négligeable

À un nombre de Reynolds suffisamment faible, les effets inertiels sont faibles par rapport aux effets visqueux. L'équation de quantité de mouvement se réduit à la forme d'écoulement de Stokes :

[ 0=-\nabla p+\mu\nabla^2\mathbf{v}+\rho\mathbf{g}. ]

Cette approximation est utile pour les fluides très visqueux, les canaux étroits et les modules microfluidiques ou de laboratoire à faible débit.

Gravité et pression hydrostatique

Si la vitesse est négligeable, l'équation de quantité de mouvement se réduit à :

[ \nabla p=\rho\mathbf{g}. ]

Cela décrit la variation de la pression hydrostatique. Dans les systèmes en écoulement, la gravité peut rester importante dans les tuyauteries verticales, les colonnes garnies, les réservoirs et les configurations liées à la convection naturelle.

Comment les équations s'appliquent aux opérations unitaires

Réseaux de tuyauteries et pompes

Les équations décrivent les distributions de vitesse et les gradients de pression dans toute la tuyauterie. Des formes intégrées ou réduites sont ensuite utilisées pour estimer les pertes de charge et déterminer les besoins en pompage.

Pour une usine pilote, les équations de champ complètes peuvent être résolues par calcul, tandis que les sections de tuyaux plus simples sont souvent représentées par des corrélations de facteur de friction dérivées de ces bilans.

Échangeurs de chaleur

Pour une modélisation strictement isotherme de l'écoulement des fluides, la continuité et la quantité de mouvement déterminent l'hydrodynamique, y compris la distribution de la vitesse, le cisaillement à la paroi et la perte de charge.

Si l'objectif inclut un changement de température ou une puissance thermique, une équation d'énergie doit également être ajoutée. Un bilan courant au niveau du flux est :

[ Q=\dot{m}C_p\Delta T, ]

tandis que la performance de l'échangeur est souvent liée par :

[ Q=UA\Delta T_{\mathrm{lm}}F. ]

Ce sont des relations thermiques supplémentaires, et non des remplacements aux équations d'écoulement.

Réacteurs et équipements de mélange

Dans les réservoirs agités et les réacteurs à écoulement, les équations de Navier-Stokes prédisent la circulation, la recirculation et les champs de vitesse liés au mélange.

La réaction, le transport d'espèces et le comportement du temps de séjour nécessitent des équations de conservation supplémentaires. Les équations d'écoulement des fluides seules ne prédisent pas la conversion ou la concentration.

Lits garnis et colonnes

Dans les équipements garnis, les équations de Navier-Stokes à l'échelle des pores peuvent être utilisées pour une modélisation détaillée, mais les modèles à l'échelle de l'ingénierie les remplacent couramment par des relations de perte de charge en milieu poreux.

Les principes sous-jacents restent la conservation de la masse et de la quantité de mouvement, le garnissage étant représenté par des termes de résistance supplémentaires.

Modules à deux fluides ou multiphasiques

Pour deux fluides immiscibles, chaque phase nécessite ses propres équations de conservation, ainsi que des conditions d'interface. À une interface liquide-liquide, les conditions typiques incluent :

[ \mathbf{v}_1=\mathbf{v}_2 ]

pour la continuité de la vitesse tangentielle et

[ \boldsymbol{\tau}_1\cdot\mathbf{n}

\boldsymbol{\tau}_2\cdot\mathbf{n} ]

pour la continuité de la contrainte visqueuse, lorsque les effets interfaciaux sont traités de manière appropriée.

Les équations pour un seul fluide incompressible sont donc insuffisantes lorsque les frontières de phase, le transfert de quantité de mouvement interphasique ou le changement de phase sont importants.

Comprendre les compromis

Isotherme ne signifie pas que l'énergie est toujours hors de propos

Un modèle isotherme suppose que la température est fixe ou que les effets thermiques n'influencent pas la solution d'écoulement. Il ne fournit pas de distributions de température, de coefficients de transfert thermique ou de puissance thermique.

Pour les échangeurs de chaleur, les réacteurs chauffés ou les systèmes présentant une forte variation de viscosité, une équation d'énergie peut être nécessaire même si l'analyse initiale de l'écoulement utilise une approximation isotherme.

L'incompressibilité a une plage de validité définie

Les liquides sont souvent modélisés comme incompressibles car les changements de masse volumique sont faibles. Les gaz peuvent également être traités comme incompressibles lorsque la vitesse est suffisamment inférieure à la vitesse du son et que les changements de densité induits par la pression sont négligeables.

L'approximation doit être reconsidérée pour les écoulements de gaz à haute vitesse, les grands changements de pression, les fortes variations de température ou les systèmes impliquant des effets liés à la compressibilité.

Une viscosité constante peut être une approximation médiocre

L'équation standard suppose une viscosité (\mu) constante. C'est raisonnable pour de nombreux liquides sur des plages de température et de taux de cisaillement modérées, mais pas pour les fluides fortement dépendants de la température ou les matériaux non newtoniens.

Si la viscosité varie spatialement, le terme visqueux doit généralement être conservé sous forme de divergence de contrainte plutôt qu'écrit simplement comme (\mu\nabla^2\mathbf{v}).

Les hypothèses d'écoulement laminaire ne sont pas requises par Navier-Stokes

Les équations de Navier-Stokes s'appliquent à l'écoulement newtonien, qu'il soit laminaire ou turbulent. Cependant, la résolution directe des échelles turbulentes est exigeante en termes de calcul, de sorte que des modèles de turbulence ou des corrélations validées expérimentalement sont souvent nécessaires pour les équipements à l'échelle pilote.

La « dissipation visqueuse » appartient au bilan énergétique

Les contraintes visqueuses affectent directement l'équation de quantité de mouvement. La chaleur générée par la dissipation visqueuse apparaît dans l'équation d'énergie et est omise d'une formulation strictement isotherme.

Cette distinction évite l'erreur courante consistant à traiter la dissipation visqueuse comme un terme supplémentaire dans l'équation de quantité de mouvement.

Faire le bon choix pour votre objectif

Le modèle directeur correct dépend de ce que l'analyse en usine pilote doit prédire.

Si votre objectif principal est les profils de vitesse et la perte de charge : Utilisez la continuité incompressible plus les équations de quantité de mouvement de Navier-Stokes, avec une masse volumique et une viscosité appropriées aux conditions de fonctionnement.
Si votre objectif principal est la performance du transfert thermique : Ajoutez une équation d'énergie et utilisez la solution d'écoulement pour déterminer la convection, la perte de charge et le comportement du transfert thermique.
Si votre objectif principal est la conversion du réacteur ou la distribution des espèces : Ajoutez des bilans d'espèces ou de réaction aux équations hydrodynamiques.
Si votre objectif principal est l'opération liquide-liquide ou gaz-liquide : Utilisez des équations de conservation spécifiques à la phase et incluez les conditions de vitesse interfaciale, de contrainte et de transfert de masse.
Si votre objectif principal est la conception rapide à l'échelle de l'équipement : Utilisez des modèles réduits validés, tels que des corrélations de facteur de friction, de milieu poreux ou de temps de séjour, tout en reconnaissant leurs limites de fonctionnement.

Avec les hypothèses clairement énoncées et les conditions aux limites correctement spécifiées, ces équations de conservation fournissent une base fiable pour analyser le comportement de l'écoulement des fluides en usine pilote.

Tableau récapitulatif :

Élément du modèle	Relation directrice	Objectif principal
Conservation de la masse	∇ · v = 0	Assure la continuité de l'écoulement incompressible
Conservation de la quantité de mouvement	ρ(∂v/∂t + v · ∇v) = −∇p + μ∇²v + ρg	Prédit les effets de vitesse, de pression, visqueux et de force massique
Viscosité newtonienne	τ = 2μD	Relie la contrainte visqueuse au taux de déformation
Conditions aux limites	Conditions d'entrée, de sortie, de paroi et initiales	Ferme le modèle mathématique
Physique supplémentaire	Équations d'énergie, d'espèces, de réaction ou multiphasiques	Étend le modèle pour le transfert thermique, les réactions ou les phases multiples

Besoin d'équipements d'usine pilote fiables pour l'enseignement de l'écoulement des fluides, le développement de procédés ou la mise à l'échelle ? Contactez LABPARK pour discuter de vos besoins. LABPARK fournit des usines pilotes d'opérations unitaires éducatives et professionnelles pour le génie chimique, les bioprocédés et la biotechnologie, ainsi que pour les applications environnementales et de traitement de l'eau. Nos systèmes modulaires aident les universités, les instituts de recherche et les entreprises à étudier la perte de charge, les profils de vitesse, le mélange, le transfert thermique, la séparation, la réaction et le contrôle des procédés avec des équipements pratiques et axés sur l'application. Indiquez-nous votre procédé cible, vos conditions de fonctionnement, vos objectifs d'enseignement ou de recherche, et l'échelle requise afin que notre équipe puisse recommander une configuration adaptée et soutenir la planification de votre usine pilote.