La formulation fonction de courant–vorticité simplifie les problèmes d’écoulement rampant en remplaçant une équation du quatrième ordre par deux équations couplées du second ordre. Pour un écoulement de Stokes stationnaire, incompressible et plan, on définit la fonction de courant de manière à satisfaire automatiquement la continuité, puis on la relie à la vorticité par (\nabla^2\psi=-\omega_z). Le problème gouvernant devient une équation de Poisson pour (\psi), couplée à une équation de Laplace pour (\omega_z), ce qui est généralement plus facile à résoudre analytiquement et numériquement.
À retenir : Dans un écoulement rampant, utilisez (\nabla^2\psi=-\omega_z) et (\nabla^2\omega_z=0) au lieu de résoudre directement l’équation biharmonique (\nabla^4\psi=0). Cela réduit l’ordre différentiel, impose naturellement la conservation de la masse et fournit un cadre pratique pour les géométries d’opérations unitaires impliquant des parois, des particules, des canaux et des milieux poreux.
Pourquoi l’écoulement rampant est particulièrement adapté
Les forces visqueuses dominent
Un écoulement rampant se produit lorsque le nombre de Reynolds est très inférieur à un :
[ \mathrm{Re}=\frac{\rho U L}{\mu}\ll 1 ]
Dans ces conditions, les termes inertiels sont négligeables devant les forces visqueuses et de pression. L’équation de la quantité de mouvement se réduit à l’équation de Stokes :
[ -\nabla p+\mu\nabla^2\mathbf{v}+\mathbf{f}=0 ]
où (\mathbf{f}) représente les forces volumiques lorsqu’elles sont pertinentes.
Les équations gouvernantes deviennent linéaires
Comme l’accélération convective non linéaire est absente, la vitesse et la pression dépendent linéairement du mouvement des parois, des gradients de pression et des forces volumiques. Cette linéarité permet la superposition pour les solutions analytiques et améliore la robustesse des solveurs numériques.
L’écoulement réagit sans retard inertiel
Dans l’approximation idéale de l’écoulement de Stokes, le champ de vitesse s’ajuste immédiatement aux variations de pression ou au mouvement des parois. Ce comportement pseudo-stationnaire est utile pour analyser les équipements entraînés lentement, les fluides très visqueux et les écoulements à l’échelle microscopique.
Comment la fonction de courant élimine une contrainte
La continuité est satisfaite automatiquement
Pour un écoulement bidimensionnel incompressible, on définit la fonction de courant (\psi) par
[ u=\frac{\partial\psi}{\partial y}, \qquad v=-\frac{\partial\psi}{\partial x} ]
L’équation de continuité est
[ \frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0 ]
La substitution des définitions de la fonction de courant satisfait cette équation identiquement. Il n’est donc pas nécessaire de résoudre séparément la continuité à chaque étape.
Les lignes de courant sont directement obtenues à partir de (\psi)
Les lignes de niveau constant de (\psi) sont des lignes de courant. La différence entre les valeurs de la fonction de courant de deux lignes de courant représente le débit volumique par unité de profondeur entre celles-ci, selon la convention de signe choisie.
Cette formulation possède donc une valeur à la fois mathématique et physique : la variable numérique décrit directement la structure de l’écoulement et le débit de transport.
La vorticité relie la rotation à la courbure de la vitesse
Pour un écoulement plan, la vorticité hors du plan est
[ \omega_z=\frac{\partial v}{\partial x} -\frac{\partial u}{\partial y} ]
En utilisant les définitions de la fonction de courant, on obtient
[ \omega_z=-\nabla^2\psi ]
ou, de manière équivalente,
[ \nabla^2\psi=-\omega_z ]
Le signe change si l’on utilise une autre convention pour la fonction de courant : la convention doit donc rester cohérente tout au long de la dérivation et du code.
Comment l’équation du quatrième ordre est découplée
Formulation directe en fonction de courant
Le calcul du rotationnel de l’équation de quantité de mouvement de Stokes élimine la pression. Pour une viscosité constante et en l’absence de contribution de forces volumiques rotationnelles, l’équation résultante pour la fonction de courant est
[ \nabla^4\psi=0 ]
Il s’agit de l’équation biharmonique. En coordonnées cylindriques ou sphériques, l’opérateur correspondant est modifié par la géométrie et s’écrit souvent de manière schématique
[ E^4\psi=0 ]
Formulation couplée du second ordre
On introduit la vorticité comme champ indépendant. Le problème du quatrième ordre devient
[ \nabla^2\psi=-\omega_z ]
et
[ \nabla^2\omega_z=0 ]
La première équation est une équation de Poisson pour la fonction de courant. La seconde est une équation de Laplace pour la vorticité.
C’est la simplification centrale : deux équations du second ordre sont généralement plus faciles à discrétiser, à résoudre et à soumettre aux conditions aux limites qu’une seule équation du quatrième ordre.
Pourquoi un ordre inférieur est avantageux
Les discrétisations du second ordre sont largement prises en charge par les méthodes des différences finies, des volumes finis et des éléments finis. Elles nécessitent également des informations aux limites sous une forme qui correspond davantage au problème physique, comme la vitesse de la paroi, l’imperméabilité ou le comportement imposé de la vorticité.
Une équation du quatrième ordre exige des conditions aux limites supplémentaires sur les dérivées, qui peuvent être difficiles à spécifier avec précision au niveau des parois solides, des angles, des interfaces et des surfaces des particules.
Avantages numériques pour les opérations unitaires
La pression est éliminée du calcul principal
La formulation fonction de courant–vorticité permet d’obtenir le champ de vitesse sans résoudre la pression comme inconnue principale. La pression peut être reconstruite ensuite si la perte de charge, la force ou la consommation de puissance sont nécessaires.
Cela est utile dans les problèmes où la structure de l’écoulement est initialement plus importante que le champ de pression détaillé, comme la circulation d’une suspension ou la recirculation autour de particules.
La formulation convient bien aux géométries structurées
Pour les canaux plans, les interstices, les cavités, les films et les passages poreux idéalisés, les deux équations peuvent être discrétisées sur des maillages structurés. Les systèmes linéaires obtenus peuvent être résolus à l’aide de méthodes standard pour les équations elliptiques.
Les étapes numériques typiques sont les suivantes :
- Appliquer les conditions aux limites pour (\psi) et la vitesse.
- Résoudre l’équation de la vorticité.
- Résoudre l’équation de Poisson pour (\psi).
- Calculer (u) et (v) à partir des dérivées de (\psi).
- Itérer jusqu’à convergence des champs couplés.
Pour un écoulement de Stokes pur, les équations sont linéaires. L’itération est principalement liée au couplage, au traitement de la vorticité aux limites, à la géométrie ou aux extensions impliquant des propriétés variables et une physique supplémentaire.
La conservation de la masse est intégrée à la reconstruction de la vitesse
Les vitesses reconstruites à partir de (\psi) satisfont automatiquement l’incompressibilité, à l’exception des erreurs de discrétisation. Cela réduit le risque d’obtenir un champ de vitesse qui viole la continuité, une préoccupation courante dans les formulations pression–vitesse.
Les conditions aux limites deviennent le principal défi numérique
La condition de non-glissement impose les deux composantes de la vitesse sur une paroi solide stationnaire. Sous forme de fonction de courant, l’imperméabilité détermine généralement une valeur constante de (\psi) le long de la paroi, tandis que la vitesse tangentielle fournit des informations sur les dérivées.
La vorticité à la paroi n’est souvent pas imposée directement. Elle doit être déduite de l’équation de la fonction de courant et de la condition de vitesse près de la paroi, à l’aide d’une discrétisation appropriée aux limites.
Avantages analytiques dans les problèmes de génie chimique
La superposition devient directe
Comme l’écoulement de Stokes est linéaire, les solutions peuvent être additionnées. Un mouvement complexe peut être décomposé en composantes plus simples, telles que la translation, la rotation, l’écoulement induit par la pression et le mouvement des parois.
Par exemple, l’écoulement autour d’une particule en mouvement peut être construit en combinant un problème d’écoulement uniforme avec un champ de perturbation satisfaisant les conditions aux limites de la particule.
La symétrie réduit le domaine
Les équations de la fonction de courant et de la vorticité préservent la symétrie géométrique. Si un canal, une disposition de particules ou une cavité est symétrique, seule une partie du domaine peut devoir être analysée.
Cela réduit la complexité algébrique et facilite la vérification des cas limites.
La séparation des variables reste utile
Dans les géométries cartésiennes, cylindriques et sphériques simples, on peut supposer des formes telles que
[ \psi(x,y)=X(x)Y(y) ]
ou utiliser des développements de Fourier, en fonctions propres ou multipolaires. Les équations du second ordre produisent des équations différentielles ordinaires plus faciles à traiter que la manipulation directe de l’équation de la fonction de courant du quatrième ordre.
Les résultats classiques des opérations unitaires sont naturellement retrouvés
Plusieurs résultats familiers peuvent être interprétés dans le même cadre :
- Un écoulement pleinement développé dans un canal possède un profil simple de fonction de courant et une vorticité déterminée par le gradient de pression.
- L’écoulement dans des interstices étroits est principalement gouverné par la résistance visqueuse et peut être approximé à l’aide d’hypothèses de lubrification.
- L’écoulement rampant autour de particules fournit les forces de traînée et les champs de perturbation utilisés dans l’analyse de la sédimentation et des suspensions.
- Les formulations en coordonnées cylindriques décrivent les écoulements idéalisés dans les capillaires, les espaces annulaires et les canaux poreux.
La forme du laplacien et de la relation de vorticité doit être modifiée de manière appropriée pour les géométries cylindriques ou sphériques ; la formule cartésienne ne peut pas être réutilisée telle quelle sans modification.
Applications aux opérations unitaires
Transformation des polymères et fluides très visqueux
Les polymères fondus et les liquides concentrés présentent souvent une viscosité élevée et des nombres de Reynolds caractéristiques faibles. La formulation aide à calculer les profils de vitesse, le cisaillement aux parois, la recirculation et le comportement du temps de séjour dans les passages entraînés lentement.
Pour les fluides non newtoniens, cependant, l’équation harmonique simple de la vorticité ne s’applique généralement plus sous sa forme originale. La viscosité dépend du taux de déformation, de sorte que les équations de quantité de mouvement et de vorticité deviennent couplées par l’intermédiaire du modèle constitutif.
Microcanaux et passages étroits
Les petites échelles de longueur réduisent le nombre de Reynolds, ce qui rend les hypothèses d’écoulement rampant courantes dans les systèmes microfluidiques et de transport à l’échelle microscopique. Les méthodes fondées sur la fonction de courant peuvent résoudre l’écoulement autour d’obstacles, dans les contractions, les expansions et à proximité de frontières mobiles.
Elles sont particulièrement utiles lorsque l’objectif principal est le profil local de l’écoulement, la structure de vitesse indépendante de la pression ou le transport de particules.
Suspensions et transport de particules
L’écoulement autour des colloïdes, des gouttelettes, des bulles et des particules en suspension de plus grande taille détermine la traînée hydrodynamique et la migration. La résolution de (\psi) et de (\omega) fournit le champ de perturbation nécessaire à l’estimation des forces et des interactions entre particules.
La méthode met également en évidence la réversibilité de l’écoulement : l’inversion du mouvement imposé à la particule ou à la paroi inverse le champ de vitesse et les lignes de courant dans la limite idéale d’un écoulement de Stokes linéaire.
Idéalisations des milieux poreux
À l’échelle des pores, les solutions d’écoulement rampant aident à expliquer la résistance visqueuse et le transport induit par la pression dans des passages simplifiés. Elles peuvent également servir de briques élémentaires locales pour des modèles de perméabilité à plus grande échelle.
Un milieu poreux réel nécessite généralement une méthode numérique résolvant la géométrie ou un modèle macroscopique tel que l’écoulement de Darcy. La formulation fonction de courant est directement applicable principalement aux géométries de pores bidimensionnelles ou spécialement réduites.
Comprendre les compromis
Les conditions aux limites de la vorticité ne sont pas toujours simples
L’introduction de la vorticité réduit l’ordre différentiel, mais ajoute un champ dont les valeurs aux limites peuvent être indirectes. De mauvaises approximations de la vorticité à la paroi peuvent produire des contraintes de cisaillement inexactes ou un comportement parasite près de la paroi.
Ce compromis est particulièrement important lorsque la contrainte de cisaillement à la paroi, la traînée ou la perte de charge constitue la grandeur d’intérêt, plutôt que le seul champ de vitesse global.
Les angles et les singularités exigent de la prudence
Les angles vifs, les expansions brusques et les zones de quasi-contact entre particules et parois peuvent générer de forts gradients de vitesse ou des singularités mathématiques. Un maillage grossier ou un traitement des limites d’ordre faible peut produire une vorticité locale trompeuse, même lorsque l’écoulement global semble raisonnable.
Un raffinement local du maillage et une comparaison attentive avec des solutions analytiques limites sont nécessaires.
La pression n’est plus immédiatement disponible
La pression disparaît lorsque l’équation de quantité de mouvement est transformée sous forme de vorticité. Si la perte de charge, la puissance de pompage ou l’équilibre des forces sont nécessaires, la pression doit être reconstruite à partir des équations originales de Stokes ou d’une procédure équivalente de post-traitement.
La formulation simplifie donc la résolution de la vitesse, mais n’élimine pas la pression de l’analyse globale d’ingénierie.
Les hypothèses peuvent ne plus être valables
Les équations de base supposent un comportement incompressible et newtonien, une viscosité constante et une inertie négligeable. Elles peuvent nécessiter des modifications pour :
- Les polymères fondus non newtoniens.
- Les fluides à viscosité variable.
- Les écoulements compressibles ou fortement dépendants de la température.
- Les géométries tridimensionnelles.
- Les effets inertiels lorsque le nombre de Reynolds n’est pas suffisamment faible.
- Les forces volumiques ou les effets multiphysiques qui contribuent à la génération de vorticité.
Les méthodes fondées sur la fonction de courant sont principalement bidimensionnelles
Une fonction de courant scalaire unique est particulièrement pratique pour un écoulement plan incompressible. Les écoulements pleinement tridimensionnels nécessitent généralement des potentiels vectoriels, des formulations vitesse–pression ou d’autres méthodes spécialisées.
La méthode est donc puissante pour les géométries réduites, mais ne doit pas être imposée à des opérations unitaires intrinsèquement tridimensionnelles.
Faire le bon choix selon votre objectif
La formulation est particulièrement efficace lorsque la géométrie et les hypothèses constitutives permettent un modèle bidimensionnel, incompressible et à faible nombre de Reynolds.
- Si votre objectif principal concerne les profils de vitesse ou les lignes de courant : résolvez (\nabla^2\psi=-\omega_z) avec (\nabla^2\omega_z=0), car la continuité est imposée automatiquement et la structure de l’écoulement découle directement de (\psi).
- Si votre objectif principal concerne la robustesse numérique : utilisez des discrétisations couplées du second ordre, avec un traitement soigneux de la vorticité aux parois et des régions d’angle.
- Si votre objectif principal concerne la dérivation analytique : exploitez la linéarité, la symétrie, la séparation des variables et la superposition avant de tenter une solution biharmonique générale.
- Si votre objectif principal concerne la perte de charge ou la puissance de pompage : utilisez la solution fonction de courant–vorticité pour la vitesse, puis reconstruisez la pression à partir de l’équation de quantité de mouvement de Stokes.
- Si votre objectif principal concerne la transformation des polymères ou la rhéologie des suspensions : conservez cette formulation comme point de départ conceptuel, mais introduisez la relation constitutive non newtonienne appropriée avant toute application quantitative.
- Si votre objectif principal concerne une opération unitaire tridimensionnelle : privilégiez une méthode tridimensionnelle vitesse–pression ou fondée sur un potentiel vectoriel, sauf si la géométrie permet une réduction bidimensionnelle justifiée.
La relation entre fonction de courant et vorticité transforme l’analyse des écoulements rampants en un problème d’ordre inférieur, cohérent avec la conservation, tout en préservant les connaissances physiques nécessaires à la conception et à l’interprétation des opérations unitaires du génie chimique.
Tableau récapitulatif :
| Aspect clé | Enseignement pratique |
|---|---|
| Équations fondamentales | Résoudre (\nabla^2\psi=-\omega_z) et (\nabla^2\omega_z=0) au lieu de l’équation biharmonique du quatrième ordre. |
| Avantage principal | Deux équations couplées du second ordre sont plus faciles à discrétiser et à soumettre aux conditions aux limites. |
| Conservation | La définition de la vitesse par l’intermédiaire de la fonction de courant satisfait automatiquement la continuité incompressible. |
| Signification physique | Les contours de la fonction de courant montrent les lignes de courant et le débit ; la vorticité représente la rotation locale du fluide. |
| Utilisations numériques | Adaptée aux canaux, aux cavités, aux interstices étroits, aux particules, aux microécoulements et aux passages poreux idéalisés. |
| Outils analytiques | La linéarité permet la superposition, la réduction par symétrie, la séparation des variables et l’utilisation de solutions classiques. |
| Principaux défis | Le traitement de la vorticité aux parois, les angles vifs, les singularités, la reconstruction de la pression et le raffinement du maillage exigent de la prudence. |
| Limites d’applicabilité | Particulièrement adaptée aux écoulements bidimensionnels, incompressibles, newtoniens, à viscosité constante et à faible nombre de Reynolds. |
| Applications étendues | Les écoulements non newtoniens, à viscosité variable, tridimensionnels ou inertiels nécessitent des formulations modifiées. |
| Résultats d’ingénierie | Utiliser la solution de vitesse pour évaluer les profils, les lignes de courant, le cisaillement, la traînée, le temps de séjour et le transport. |
Besoin d’un équipement fiable pour l’enseignement, la recherche ou le développement de procédés ? Contactez LABPARK pour discuter de vos besoins. LABPARK fournit des installations pilotes d’opérations unitaires éducatives et professionnelles 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 rampants, à étudier le transport induit par la pression, à analyser le comportement des particules et des suspensions et à relier les modèles théoriques aux données pratiques d’exploitation. Nous pouvons vous aider à sélectionner une configuration d’installation pilote adaptée aux objectifs de votre cours, à la capacité de votre laboratoire, à vos besoins de recherche et à vos objectifs de développement de procédés.