Les équations de Navier–Stokes prédisent la perte de charge en reliant le mouvement du fluide, les gradients de pression, la viscosité et la résistance des parois ou des équipements. Dans une installation pilote de génie chimique, les ingénieurs les utilisent pour calculer les profils de vitesse, les contraintes de cisaillement et les pertes de pression dans les conduites, raccords, réacteurs, échangeurs de chaleur, lits garnis et autres équipements d’écoulement. Pour les géométries simples, les équations peuvent être résolues analytiquement ; pour les équipements complexes, elles sont généralement simplifiées sous forme de corrélations d’ingénierie ou résolues numériquement à l’aide de la mécanique des fluides numérique (CFD).

À retenir : Les équations de Navier–Stokes fournissent la physique gouvernante, mais les prédictions pratiques pour les installations pilotes nécessitent également les propriétés du fluide, les conditions aux limites, l’analyse du régime d’écoulement, la géométrie ainsi que des corrélations de friction ou d’équipement validées expérimentalement.

Ce que représentent les équations

Le bilan de quantité de mouvement

Pour un fluide newtonien incompressible à propriétés constantes, l’équation de quantité de mouvement s’écrit généralement comme suit :

[ \rho \left(\frac{\partial \mathbf{u}}{\partial t}

\mathbf{u}\cdot\nabla\mathbf{u}\right) = -\nabla P +\mu\nabla^2\mathbf{u} +\rho\mathbf{g} ]

Ici, (\rho) représente la masse volumique, (\mathbf{u}) la vitesse, (P) la pression, (\mu) la viscosité dynamique et (\mathbf{g}) la gravité.

Le membre de gauche décrit l’accélération du fluide. Les termes du membre de droite représentent les forces de pression, les forces visqueuses et les forces volumiques, telles que la gravité.

Conservation de la masse

L’équation de quantité de mouvement doit être couplée à la conservation de la masse. Pour un fluide incompressible :

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

Cette condition garantit que le champ de vitesse calculé ne crée ni ne détruit de volume de fluide.

Pourquoi la viscosité est importante

La viscosité crée une contrainte de cisaillement entre les couches adjacentes de fluide ainsi qu’entre le fluide et une paroi solide. Selon la loi de Newton sur la viscosité, cette contrainte est liée au gradient de vitesse.

Cette relation permet aux équations de prédire l’intensité avec laquelle la paroi s’oppose à l’écoulement et la manière dont cette résistance se traduit par une diminution de pression le long de l’équipement.

Comment apparaît la perte de charge

Écoulement pleinement développé dans une conduite

Dans une conduite longue, droite et de diamètre constant, le profil de vitesse devient pleinement développé. Le gradient de pression équilibre la résistance visqueuse à la paroi.

Pour un écoulement laminaire permanent dans une conduite circulaire, la solution conduit à la relation de Hagen–Poiseuille :

[ \Delta P

\frac{128\mu LQ}{\pi D^4} ]

où (L) est la longueur de la conduite, (Q) le débit volumique et (D) le diamètre de la conduite.

Cela montre que la perte de charge en régime laminaire est directement proportionnelle à la viscosité, à la longueur et au débit. La forte dépendance en (D^4) explique également pourquoi de petites variations du diamètre peuvent avoir un effet majeur sur les besoins de pompage.

Forme d’ingénierie pour les pertes dans les conduites

Pour la conception pratique, la perte de pression est souvent exprimée à l’aide d’un facteur de friction de Darcy :

[ \Delta P_f

f_D\frac{L}{D}\frac{\rho u^2}{2} ]

Ici, (u) est la vitesse moyenne et (f_D) le facteur de friction de Darcy.

Cette forme est cohérente avec le bilan de quantité de mouvement de Navier–Stokes tout en permettant aux ingénieurs de représenter les effets de paroi au moyen d’un facteur de friction obtenu par la théorie, des corrélations, des abaques ou des mesures.

Pertes singulières et pertes des équipements

Les vannes, coudes, contractions, expansions, zones d’entrée et équipements de procédé créent des pertes supplémentaires :

[ \Delta P_{\text{minor}}

K\frac{\rho u^2}{2} ]

La perte de charge totale est donc supérieure à la contribution de la conduite droite :

[ \Delta P_{\text{total}}

\Delta P_f + \sum \Delta P_{\text{minor}} + \Delta P_{\text{equipment}} ]

Dans une installation pilote, ce profil de pression complet sert à déterminer si les pompes peuvent maintenir le débit requis et si un composant devient un goulot d’étranglement hydraulique.

Détermination du régime d’écoulement

Nombre de Reynolds

Le nombre de Reynolds compare les effets d’inertie et les effets visqueux :

[ Re=\frac{\rho uD}{\mu} ]

Pour un écoulement interne en conduite, un nombre (Re) inférieur à environ 2 000 est généralement considéré comme laminaire, tandis que les valeurs supérieures à environ 3 000 sont généralement considérées comme turbulentes. La zone située entre ces valeurs est transitoire, de sorte que les prédictions dans cette zone nécessitent une attention particulière.

Comportement laminaire

En régime laminaire, les forces visqueuses dominent et le profil de vitesse est relativement ordonné. Pour une conduite circulaire, le facteur de friction de Darcy est :

[ f_D=\frac{64}{Re} ]

La rugosité de la paroi a peu d’influence pratique lorsque l’écoulement est pleinement laminaire. La perte de charge varie donc approximativement de manière linéaire avec le débit.

Comportement turbulent

En régime turbulent, les mouvements fluctuants augmentent le transfert de quantité de mouvement et le cisaillement à la paroi. Le facteur de friction dépend à la fois du nombre de Reynolds et de la rugosité relative.

Une représentation courante est la relation de Colebrook :

[ \frac{1}{\sqrt{f_D}}

-2\log_{10} \left( \frac{\varepsilon/D}{3.7} + \frac{2.51}{Re\sqrt{f_D}} \right) ]

où (\varepsilon/D) est la rugosité relative. Les abaques de Moody et les approximations explicites sont couramment utilisés pour évaluer cette relation.

Dans le régime turbulent pleinement rugueux, le facteur de friction devient largement indépendant du nombre de Reynolds et dépend principalement de la rugosité relative.

Application des équations dans une installation pilote

Mesure expérimentale de la pression

Les installations pilotes mesurent couramment la pression en deux points à l’aide de capteurs de pression différentielle ou de manomètres. La différence de pression mesurée peut être comparée à la perte calculée sur une longueur de conduite ou une section d’équipement connue.

Cette comparaison permet de vérifier si les propriétés du fluide, le facteur de friction, la géométrie et le régime de fonctionnement supposés représentent correctement le système réel.

Établissement d’un profil de pression

Les ingénieurs commencent à partir d’une limite de pression connue et soustraient successivement les pertes au niveau de chaque élément. Le profil peut inclure les conduites, vannes, réchauffeurs, réacteurs, échangeurs de chaleur, filtres et équipements de séparation.

À chaque emplacement, la pression calculée est comparée à la pression de fonctionnement requise et à la pression maximale admissible de conception du composant.

Mise à l’échelle des résultats de l’installation pilote

Pour transférer les résultats entre une unité pilote et un système de plus grande taille, les ingénieurs préservent les groupes sans dimension qui contrôlent le comportement pertinent. L’égalité du nombre de Reynolds est importante lorsque les effets visqueux et inertiels gouvernent l’écoulement.

Pour la mise à l’échelle des pertes de pression, le nombre d’Euler est utile :

[ Eu=\frac{\Delta P}{\rho u^2} ]

Si le modèle et le système à pleine échelle présentent une géométrie comparable et des groupes sans dimension directeurs concordants, les pertes de pression mesurées dans l’installation pilote peuvent contribuer aux prédictions pour le système de plus grande taille.

La similitude n’est pas automatique. Les différences de rugosité, de raccords, d’écoulement en développement, d’éléments internes, de propriétés du fluide et de régime de fonctionnement peuvent limiter la précision d’un changement d’échelle direct.

Extension du modèle au-delà des conduites simples

Lits garnis et milieux poreux

L’écoulement à travers une colonne garnie rencontre la résistance de nombreuses particules plutôt que celle d’une seule paroi continue de conduite. À des vitesses suffisamment faibles, les termes d’accélération peuvent être négligeables, ce qui produit un équilibre entre la pression, la résistance visqueuse et les forces volumiques.

Pour la conception pratique, la perte de charge dans un lit garni est généralement représentée par une corrélation de milieu poreux validée expérimentalement, telle que l’équation d’Ergun, plutôt que par la résolution des équations de Navier–Stokes dans le fluide libre autour de chaque particule.

Mélange et réacteurs

Dans les cuves agitées et les réacteurs, les équations peuvent prédire les schémas de circulation, le cisaillement local, les zones de mélange et les champs de pression. Comme la géométrie et la turbulence sont complexes, les ingénieurs utilisent couramment des modèles CFD étayés par des mesures réalisées dans l’installation pilote.

Le modèle doit inclure la géométrie de l’agitateur, sa rotation, les propriétés du fluide, le traitement de la turbulence ainsi que les conditions appropriées à l’entrée, à la sortie et aux parois.

Systèmes gazeux compressibles

Pour les gaz soumis à des variations de pression importantes, la masse volumique peut varier considérablement. La forme incompressible est alors insuffisante et l’équation de quantité de mouvement doit être couplée à l’équation de continuité, à un bilan énergétique et à une équation d’état.

Cela concerne les compresseurs, les gazoducs, les dispositifs de réduction de pression et les systèmes pilotes de séparation des gaz.

Comprendre les compromis

Les solutions analytiques sont limitées

Les solutions exactes ne sont disponibles que pour des cas idéalisés, tels que l’écoulement laminaire permanent et pleinement développé dans une conduite. Les installations pilotes réelles comprennent des coudes, des raccords, des variations de diamètre, des structures internes complexes et des écoulements non uniformes.

Les corrélations d’ingénierie étendent la théorie à ces cas, mais leur précision dépend de la correspondance entre les conditions de fonctionnement et leur domaine de validation.

La turbulence nécessite une modélisation

Le champ de vitesse turbulent dépendant du temps contient de nombreuses échelles spatiales et temporelles qu’il est généralement trop coûteux de résoudre directement. La CFD utilise donc des modèles de turbulence ou des équations moyennées pour estimer leur effet global.

Différents modèles peuvent produire des prédictions différentes des pertes de pression et du mélange, notamment à proximité des zones de décollement, des équipements rotatifs et des structures internes complexes.

La perte de charge n’est pas toujours proportionnelle au carré du débit

La relation approximative

[ \Delta P\propto \dot{m}^{,2} ]

s’applique principalement au frottement turbulent et à de nombreuses pertes locales lorsque les propriétés du fluide et la géométrie restent constantes. En régime laminaire, la perte de charge est approximativement proportionnelle au débit massique ou volumique.

Les écoulements transitoires, les écoulements multiphasiques, les fluides non newtoniens et les variations de masse volumique peuvent produire un comportement sensiblement différent.

Les mesures nécessitent une interprétation

Une pression différentielle mesurée peut inclure des effets d’altitude, des effets d’accélération, des décalages de l’instrument, des effets d’entrée ou des pertes situées en dehors de la section d’essai visée. L’emplacement et l’étalonnage des capteurs sont donc importants.

Des écarts inattendus peuvent indiquer un encrassement, une obstruction, une fuite, des variations de rugosité, une entraînement de gaz ou une hypothèse incorrecte concernant le régime d’écoulement.

Faire le bon choix en fonction de votre objectif

Utilisez le cadre de Navier–Stokes avec le modèle validé le plus simple permettant de représenter la physique dominante.

Si votre objectif principal est la perte de charge dans une conduite droite : Déterminez le nombre de Reynolds, calculez le facteur de friction approprié et incluez les raccords ainsi que les variations d’altitude dans le bilan de pression total.
Si votre objectif principal est la prédiction du passage de l’échelle pilote à l’échelle industrielle : Faites correspondre le nombre de Reynolds et les autres groupes sans dimension directeurs lorsque cela est possible, puis vérifiez la mise à l’échelle à partir des données mesurées de pression et de débit.
Si votre objectif principal est la performance d’un lit garni ou d’un filtre : Utilisez une corrélation validée de perte de charge pour milieu poreux et tenez compte de la structure des particules, de la perméabilité, de l’encrassement et du régime d’écoulement.
Si votre objectif principal est le mélange ou un équipement complexe : Utilisez la CFD ou des modèles d’ordre réduit, avec des données expérimentales de pression, de vitesse ou de mélange pour validation.
Si votre objectif principal est la capacité de l’équipement : Établissez un profil de pression aux débits de référence et aux débits augmentés, puis comparez chaque pression de fonctionnement aux limites de la pompe et à la pression de conception de l’équipement.

Les équations de Navier–Stokes fournissent le fondement physique, tandis que l’analyse sans dimension, les corrélations d’ingénierie, les méthodes numériques et les mesures réalisées dans l’installation pilote transforment ce fondement en décisions de procédé fiables.

Tableau récapitulatif :

Domaine clé	Principe principal	Application dans l’installation pilote
Physique gouvernante	La conservation de la quantité de mouvement et de la masse relie la pression, la vitesse, la viscosité et la gravité.	Prédire les champs d’écoulement, les contraintes de cisaillement et la distribution de pression.
Écoulement en conduite droite	La relation de Hagen–Poiseuille s’applique à l’écoulement laminaire ; Darcy–Weisbach permet les calculs pratiques de pertes.	Estimer la perte de charge dans les conduites et les besoins de pompage.
Régime d’écoulement	Le nombre de Reynolds distingue les comportements laminaires, transitoires et turbulents.	Sélectionner les méthodes et corrélations appropriées pour le facteur de friction.
Pertes des équipements	Les vannes, raccords, échangeurs de chaleur, réacteurs et filtres ajoutent une résistance locale ou liée à l’équipement.	Établir un profil de pression complet et identifier les goulots d’étranglement.
Systèmes complexes	La CFD et les corrélations d’ingénierie validées étendent les équations de Navier–Stokes aux géométries complexes.	Analyser le mélange, les cuves agitées, les réacteurs et les schémas d’écoulement internes.
Mise à l’échelle	Le nombre de Reynolds, le nombre d’Euler, la géométrie, la rugosité et les propriétés du fluide influencent la similitude.	Transférer plus fiablement les données pilotes vers des systèmes de procédé de plus grande taille.
Validation expérimentale	Les mesures de pression différentielle testent les hypothèses du modèle et détectent l’encrassement ou les obstructions.	Améliorer la confiance dans la conception et le dépannage opérationnel.

Besoin de données fiables sur les écoulements de fluides pour votre prochain projet pilote ? 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 pilotes aident les universités, les instituts de recherche et les entreprises à démontrer le comportement des pertes de charge, à évaluer les pompes et les équipements de procédé, à valider les corrélations d’ingénierie, à soutenir les études CFD et à générer des données pratiques de mise à l’échelle. Collaborez avec LABPARK pour développer des équipements de formation et de recherche axés sur les applications, qui améliorent l’apprentissage pratique, le développement des procédés et la prise de décision en ingénierie.