Les gradients de pression et les conditions aux limites interfaciales sont évalués en combinant les pertes de charge mesurées avec des modèles d'écoulement basés sur la conservation. Dans une installation pilote liquide-liquide, les ingénieurs mesurent la pression en plusieurs points axiaux, calculent le gradient de pression à partir de la perte de charge sur la distance, et le comparent aux prédictions du bilan de quantité de mouvement. Le modèle applique des conditions de non-glissement aux parois solides, la continuité de la vitesse tangentielle à l'interface fluide-fluide, et la continuité de la contrainte de cisaillement tangentielle ; pour une interface plane et stable, les deux couches liquides partagent le même gradient de pression axial.

Le gradient de pression contrôle les profils de vitesse et la perte de charge, tandis que les conditions aux limites interfaciales déterminent comment la quantité de mouvement et le transfert sont échangés entre les liquides non miscibles. Une évaluation fiable nécessite de traiter les effets de paroi, la géométrie de l'interface, la viscosité, la tension superficielle et le mouvement interfacial possible comme un problème couplé.

Établissement du gradient de pression

Mesure de la perte de charge le long de l'installation pilote

Les gradients de pression sont généralement évalués à l'aide de prises de pression, de transducteurs ou de manomètres positionnés à des emplacements axiaux connus. Pour deux points de mesure,

[ \frac{dP}{dz} \approx \frac{P_2-P_1}{z_2-z_1} ]

où (P_1) et (P_2) sont les pressions mesurées et (z_1) et (z_2) sont leurs positions axiales.

Une valeur négative indique une perte de charge dans le sens de l'écoulement. En pratique, plusieurs points de mesure sont préférables car ils révèlent les effets d'entrée, les régions pleinement développées, les pertes locales et la récupération de pression.

Application du bilan de quantité de mouvement

Pour un écoulement stable et pleinement développé dans un canal droit, l'équation de quantité de mouvement axiale pour chaque couche liquide a la forme générale

[ 0=-\frac{dP}{dz}+\mu_i\frac{d^2u_i}{dy^2}+\rho_i g_z ]

où (u_i), (\mu_i) et (\rho_i) sont la vitesse, la viscosité et la densité de la phase (i).

Si la gravité n'a pas de composante axiale, le gradient de pression entraîne la diffusion visqueuse de la quantité de mouvement. Les différences de viscosité produisent donc des courbures de profil de vitesse différentes, même si le gradient de pression axial est commun aux deux phases.

Distinction entre gradients de pression et effets hydrostatiques

Lorsque le canal est incliné ou vertical, le changement de pression mesuré inclut à la fois des contributions frictionnelles et hydrostatiques. Les données de pression doivent donc être corrigées de l'élévation avant de comparer les mesures avec un modèle d'écoulement pleinement développé.

Pour une phase de densité (\rho), la contribution hydrostatique dépend de (\rho g) et du changement d'élévation. Dans un écoulement liquide-liquide stratifié, cette correction peut différer entre les sections si l'arrangement des phases ou la rétention (holdup) change.

Application des conditions de paroi et d'interface

Condition de non-glissement sur paroi solide

À chaque paroi fixe du canal, la vitesse du liquide est fixée à zéro :

[ u_i=0 ]

Cette condition s'applique indépendamment à la phase en contact avec la paroi. Elle établit le comportement de la couche limite et affecte fortement la perte de charge totale, en particulier dans les canaux étroits à l'échelle pilote.

Si la paroi est mobile, poreuse ou partiellement glissante, la condition de non-glissement doit être remplacée par une condition appropriée de quantité de mouvement ou de perméation à la paroi.

Continuité de la vitesse tangentielle

À une interface liquide-liquide, la vitesse tangentielle est normalement continue :

[ u_{t,1}=u_{t,2}=u_{t,I} ]

Cela signifie que les deux liquides n'ont pas de saut fini de vitesse tangentielle à l'interface sous l'hypothèse habituelle du milieu continu.

La vitesse interfaciale commune peut être nulle pour une interface stationnaire ou non nulle lorsque l'interface est entraînée vers l'aval.

Continuité de la contrainte de cisaillement tangentielle

Pour une interface propre sans contraintes interfaciales supplémentaires, les contraintes de cisaillement tangentielles correspondent :

[ \mu_1\frac{\partial u_{t,1}}{\partial n}

\mu_2\frac{\partial u_{t,2}}{\partial n} ]

Les gradients de vitesse peuvent différer car les viscosités diffèrent. Un liquide plus visqueux nécessite généralement un gradient de vitesse plus faible pour transmettre la même contrainte de cisaillement.

Les tensioactifs, la viscosité interfaciale, l'adsorption ou les contraintes de Marangoni modifient cette condition. Dans ces cas, le saut de contrainte tangentielle est équilibré par la force interfaciale pertinente plutôt que d'être fixé à zéro.

Évaluation du bilan des contraintes normales

Approximation de l'interface plane

Pour une interface plane et stable avec une variation négligeable de la contrainte interfaciale normale, la pression est continue à travers l'interface :

[ P_1=P_2 ]

Par conséquent, dans un écoulement parallèle pleinement développé, le gradient de pression axial doit être le même dans les deux couches liquides :

[ \left(\frac{dP}{dz}\right)_1

\left(\frac{dP}{dz}\right)_2 ]

Cela ne signifie pas que les phases ont des vitesses égales ou des taux de cisaillement égaux. Leurs distributions de vitesse restent dépendantes de la viscosité, de l'épaisseur de la couche et du gradient de pression commun.

Correction de l'interface courbe

Pour une interface courbe, la tension superficielle crée un saut de pression décrit par la relation de Young-Laplace :

[ P_1-P_2=\gamma\kappa ]

Pour une interface sphérique, (\kappa=2/R), ce qui donne

[ P_1-P_2=\frac{2\gamma}{R} ]

où (\gamma) est la tension interfaciale et (R) est le rayon de courbure.

Ainsi, la pression elle-même n'a pas besoin d'être continue à travers une interface de gouttelette ou de bulle courbe. Le saut de pression est un effet de contrainte normale et ne doit pas être confondu avec l'exigence d'un gradient de pression axial commun dans un modèle d'écoulement stratifié plat et pleinement développé.

Inclusion des contraintes normales visqueuses

Lorsque l'interface est fortement déformée ou que l'écoulement est hautement dynamique, le bilan des contraintes normales inclut les contraintes normales visqueuses ainsi que la pression et la tension superficielle. Un modèle pratique doit alors résoudre la courbure de l'interface et les gradients de vitesse locaux plutôt que d'imposer une simple continuité de pression.

Modélisation du transport interfacial

Interfaces stationnaires et mobiles

Pour une interface stationnaire, le transport à travers la frontière est décrit par la différence entre les flux arrivant de chaque phase. Pour une interface mobile, la vitesse de l'interface doit être incluse car le matériau peut être transporté à la fois par le mouvement de la phase et le mouvement de l'interface.

Le bilan ponctuel général est

[ \left[ \left(\mathbf{F}-b\mathbf{v}_I\right)_B

\left(\mathbf{F}-b\mathbf{v}_I\right)_A \right]\cdot\mathbf{n}_I

B_S ]

où (\mathbf{F}) est le flux total, (b) est la quantité extensive transportée par unité de volume, (\mathbf{v}_I) est la vitesse de l'interface et (B_S) est un terme source interfacial.

Contributions convectives et diffusives

Avec le flux total décomposé en parties convective et diffusive,

[ \mathbf{F}=b\mathbf{v}+\mathbf{f} ]

le bilan d'interface devient

[ \left[ \left(\mathbf{f}+b(\mathbf{v}-\mathbf{v}_I)\right)_B

\left(\mathbf{f}+b(\mathbf{v}-\mathbf{v}_I)\right)_A \right]\cdot\mathbf{n}_I

B_S ]

Cette forme est utilisée pour représenter le transfert de masse interfacial, le transfert de chaleur, la réaction chimique, l'accumulation de tensioactifs et le changement de phase.

Relation entre l'interface et la surface de contact

Le taux de transport total dépend non seulement du flux interfacial local mais aussi de la surface interfaciale disponible. Pour les gouttelettes ou bulles dispersées, une estimation couramment utilisée est

[ a=\frac{6(1-\epsilon)}{d_{\mathrm{av}}} ]

où (a) est la surface interfaciale spécifique, (\epsilon) est la porosité ou fraction de vide de la phase continue, et (d_{\mathrm{av}}) est le diamètre moyen de la phase dispersée.

Des gouttelettes ou bulles plus petites augmentent la surface interfaciale par unité de volume, à condition que la dispersion reste stable et que la géométrie supposée soit applicable.

Comprendre les compromis

L'interprétation du gradient de pression peut être trompeuse

Une perte de charge mesurée à travers une installation pilote peut inclure le frottement aux parois, les pertes à l'entrée et à la sortie, les raccords, les changements d'élévation, la redistribution des phases et la déformation interfaciale. Traiter toute la perte mesurée comme un gradient frictionnel pleinement développé peut produire des prédictions incorrectes de vitesse et de rétention.

Les prises de pression doivent donc être placées dans une section géométriquement uniforme, et les pertes locales doivent être mesurées ou modélisées séparément.

Les hypothèses d'interface plane ont des limites

Le résultat du gradient de pression axial partagé s'applique plus directement à un écoulement parallèle, stable et pleinement développé avec une interface suffisamment plane. Il devient moins fiable près des entrées, sorties, coudes, contractions, gouttelettes dispersées, vagues ou changements rapides de rétention de phase.

Dans ces régions, la position et la courbure de l'interface doivent être résolues conjointement avec les champs de vitesse et de pression.

Les concepts de couche limite nécessitent la bonne géométrie

Un gradient de pression favorable accélère généralement l'écoulement et retarde le décollement de la couche limite, tandis qu'un gradient de pression défavorable ralentit l'écoulement et peut favoriser le décollement. Ces concepts sont importants près des coudes, expansions, contractions et corps internes, mais ils ne remplacent pas les conditions interfaciales biphasiques utilisées dans un canal stratifié droit.

Par exemple, un gradient de pression défavorable le long d'une paroi solide peut créer un décollement même lorsque l'interface liquide-liquide elle-même est correctement modélisée.

Les propriétés interfaciales peuvent varier pendant le fonctionnement

La concentration en tensioactifs, la tension interfaciale, la taille des gouttelettes et la composition des phases peuvent changer pendant l'extraction ou la réaction. Ces changements peuvent introduire des contraintes de Marangoni, altérer la mobilité de l'interface et changer la surface interfaciale spécifique.

Un modèle d'interface propre et constante est donc une base de référence utile, mais les données pilotes doivent être vérifiées pour détecter des signes de contamination, de coalescence, de rupture ou de réaction interfaciale.

Faire le bon choix pour votre objectif

La méthode d'évaluation doit correspondre à la question physique et au niveau de détail requis.

Si votre objectif principal est la perte de charge : Mesurez les pressions sur une section axiale uniforme, corrigez l'élévation et les pertes locales, et comparez le gradient résultant avec le bilan de quantité de mouvement biphasique.
Si votre objectif principal est la distribution de vitesse : Appliquez des conditions de non-glissement aux parois et faites correspondre la vitesse tangentielle interfaciale et la contrainte de cisaillement en utilisant le gradient de pression commun mesuré ou modélisé.
Si votre objectif principal est le transfert de masse interfacial : Utilisez le bilan de conservation de l'interface mobile, incluez les flux convectifs et diffusifs, et estimez la surface interfaciale disponible à partir de la taille des gouttelettes ou des bulles.
Si votre objectif principal est la stabilité de l'interface : Incluez le bilan des contraintes normales, la tension superficielle, la courbure, la rétention de phase et tout gradient de pression défavorable pouvant causer une déformation ou un décollement.
Si votre objectif principal est la CFD ou la simulation de processus : Définissez explicitement les conditions de paroi, d'entrée, de sortie et interfaciales, puis validez les prédictions de pression, de vitesse, de distribution de phase et de taux de transfert par rapport aux mesures pilotes.

Une évaluation solide de l'installation pilote traite les gradients de pression et les conditions interfaciales comme des contraintes de conservation couplées plutôt que comme des ajustements empiriques indépendants.

Tableau récapitulatif :

Focus de l'évaluation	Méthode clé	Considérations importantes
Gradient de pression	Mesurer la pression en plusieurs points axiaux et calculer ΔP/Δz	Corriger l'élévation, les pertes locales, les effets d'entrée et la récupération de pression
Conditions de paroi	Appliquer le non-glissement ou des conditions de glissement/perméation appropriées	Les effets de paroi influencent fortement les profils de vitesse et la perte de charge
Vitesse d'interface	Faire correspondre la vitesse tangentielle à travers l'interface liquide-liquide	Tenir compte des interfaces stationnaires ou mobiles
Cisaillement d'interface	Faire correspondre la contrainte de cisaillement tangentielle pour une interface propre	Inclure les effets de tensioactif, de Marangoni ou de viscosité interfaciale le cas échéant
Contrainte normale	Utiliser la continuité de pression pour les interfaces planes ou le saut de pression de Young-Laplace pour les interfaces courbes	Inclure les contraintes normales visqueuses lors d'une forte déformation ou d'un écoulement dynamique
Transport interfacial	Équilibrer les flux convectifs et diffusifs à travers l'interface	Considérer le mouvement de l'interface, la réaction, le changement de phase et la surface interfaciale
Validation du modèle	Comparer les prédictions CFD ou du bilan de quantité de mouvement avec les données pilotes	Valider la pression, la vitesse, la rétention, la distribution de phase et les taux de transfert

Planifiez une installation pilote multiphasique plus fiable

Besoin d'évaluer la perte de charge, la rétention de phase, le transfert de masse interfacial ou la performance de mise à l'échelle ? 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 la biotechnologie, ainsi que les applications de traitement de l'eau et de l'environnement. Nos systèmes aident les universités, les instituts de recherche et les entreprises à mener des formations pratiques, à valider des modèles basés sur la conservation, à étudier le transport liquide-liquide et multiphasique, et à générer des données fiables pour le développement de processus et la mise à l'échelle.