🔢 mathematics

Navier-Stokes-Cahn-Hilliard system in a $3$D perforated domain with free slip and source term: Existence and homogenization

Cet article établit l'existence de solutions faibles pour un système Navier-Stokes-Cahn-Hilliard en dimension trois dans un domaine perforé avec des conditions aux limites de glissement libre et une source non conservative, puis dérive par homogénéisation périodique deux modèles macroscopiques distincts selon l'échelle de la force capillaire, tout en démontrant la convergence de l'énergie libre vers une fonctionnelle homogénéisée.

Auteurs originaux : Amartya Chakrabortty, Haradhan Dutta, Hari Shankar Mahato

Publié 2026-03-17

📖 5 min de lecture🧠 Analyse approfondie

CC BY 4.0

Auteurs originaux : Amartya Chakrabortty, Haradhan Dutta, Hari Shankar Mahato

Article original sous licence CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Ceci est une explication générée par l'IA de l'article ci-dessous. Elle n'a pas été rédigée ni approuvée par les auteurs. Pour une précision technique, consultez l'article original. Lire la clause de non-responsabilité complète

Imaginez que vous essayez de comprendre comment l'huile et l'eau se mélangent (ou se séparent) lorsqu'elles traversent une éponge très fine, comme celle d'une éponge de cuisine, mais à une échelle microscopique où les trous sont minuscules et répétés à l'infini. C'est exactement ce que cette recherche tente de modéliser.

Voici une explication simple de ce travail scientifique, imagée pour tout le monde :

1. Le Défi : L'Éponge et les Deux Fluides

Les auteurs étudient un mélange de deux liquides (comme de l'huile et de l'eau) qui ne veulent pas vraiment se mélanger. Ils utilisent un modèle mathématique appelé Navier-Stokes-Cahn-Hilliard.

Navier-Stokes : C'est la loi qui régit comment les fluides bougent (comme l'eau qui coule dans une rivière).
Cahn-Hilliard : C'est la règle qui explique comment les deux liquides décident de se séparer ou de former des gouttes (comme quand on secoue une vinaigrette).

Le problème se pose dans un milieu poreux (une éponge géante). Les fluides doivent se faufiler entre des obstacles solides (les trous de l'éponge). De plus, les parois de ces obstacles sont "glissantes" : le fluide peut glisser le long des parois sans frottement, comme un patineur sur une glace parfaite, au lieu de s'arrêter net (ce qui est le cas habituel).

2. La Magie des Mathématiques : Deux Scénarios Possibles

Le point le plus intéressant de cette étude est qu'ils ont découvert que le comportement global du mélange dépend d'un seul "bouton de réglage" : la force de tension de surface (la force qui pousse l'huile et l'eau à rester séparées).

Imaginez que cette force est un bouton de volume sur une chaîne stéréo. Les chercheurs ont regardé ce qui se passe quand on tourne ce bouton à deux endroits différents :

Scénario A : Le Bouton est à Zéro (Le Régime "Calme")

Si la force de séparation est très faible (presque nulle), le système se comporte comme un écoulement très lent et tranquille (type "Stokes").

L'analogie : Imaginez une foule de gens marchant très lentement dans un couloir bondé. Personne ne court, personne ne bouscule les autres. Les mouvements sont lents, prévisibles et dominés par la friction.
Résultat : À grande échelle, on ne voit pas de tourbillons ni de courants rapides. Le mélange suit simplement les trous de l'éponge sans créer de mouvement propre. C'est un modèle simplifié, idéal pour des écoulements très lents.

Scénario B : Le Bouton est à un Niveau Moyen (Le Régime "Équilibré")

Si la force de séparation est juste assez forte (ni nulle, ni infinie), la magie opère : le système devient un vrai moteur dynamique (type "Navier-Stokes").

L'analogie : Imaginez maintenant que les gens dans le couloir commencent à courir et à se pousser. La tension entre les deux groupes (huile vs eau) crée de l'énergie. Les gouttes d'huile qui se forment et se déplacent créent des tourbillons et des courants qui entraînent le reste du fluide.
Résultat : À grande échelle, on obtient un modèle complexe où le fluide a sa propre inertie. Il y a des courants, des tourbillons et une interaction dynamique entre la séparation des phases et le mouvement du fluide. C'est comme si l'éponge elle-même commençait à "respirer" ou à créer des courants à cause de la façon dont les fluides se séparent.

3. La Méthode : Le "Zoom" Mathématique

Comment ont-ils fait pour passer de l'échelle microscopique (les trous de l'éponge) à l'échelle macroscopique (l'éponge entière) ?
Ils ont utilisé une technique appelée homogénéisation.

L'analogie : Imaginez que vous regardez une image en très haute définition (les pixels individuels). C'est trop complexe pour comprendre l'image globale. Les chercheurs ont utilisé un outil mathématique (l'opérateur de "dépliement") qui agit comme un zoom arrière magique. Au lieu de compter chaque pixel, ils calculent la "moyenne" de ce qui se passe dans un petit carré, puis ils appliquent cette moyenne à toute l'image.
Cela leur permet de créer de nouvelles équations qui décrivent le comportement global de l'éponge sans avoir à simuler chaque petit trou individuellement.

4. Pourquoi c'est Important ?

Cette recherche est cruciale pour des applications réelles :

Récupération du pétrole : Comprendre comment l'eau injectée pousse le pétrole hors des pores de la roche.
Biologie : Modéliser comment les nutriments circulent dans les tissus biologiques ou les hydrogels.
Matériaux : Concevoir de meilleurs filtres ou matériaux composites.

En Résumé

Cette équipe a prouvé mathématiquement qu'on peut prédire avec certitude comment un mélange de fluides se comportera dans une éponge complexe. Leurs résultats montrent que selon la "force de séparation" des fluides, le système peut soit rester calme et lent, soit devenir un moteur dynamique complexe. Ils ont fourni les équations exactes pour décrire ces deux mondes, offrant ainsi une boussole précieuse pour les ingénieurs et les scientifiques qui travaillent sur les fluides dans les milieux poreux.

1. Problématique et Contexte

L'article étudie un modèle d'interface diffuse pour un mélange binaire incompressible évoluant dans un milieu poreux périodiquement perforé en dimension 3 ( $\Omega^\varepsilon_p \subset \mathbb{R}^3$ ). Le système couplé est de type Navier-Stokes-Cahn-Hilliard (NSCH) et inclut des caractéristiques physiques et géométriques complexes :

Viscosité variable : Le tenseur de viscosité $A^\varepsilon$ dépend de la géométrie microscopique rapide.
Terme source non conservatif : L'équation de Cahn-Hilliard contient un terme source $G(\phi^\varepsilon)$ (bi-Lipschitz), ce qui brise la conservation de la masse (la moyenne spatiale de la concentration n'est ni nulle ni constante dans le temps).
Conditions aux limites mixtes :
- Sur la frontière extérieure $\partial\Omega$ : Condition de non-glissement (no-slip).
- Sur les inclusions solides $\Gamma^\varepsilon_s$ : Condition de glissement libre (free-slip / Navier slip) avec contrainte tangentielle nulle ( $u^\varepsilon \cdot n^\varepsilon = 0$ et $[A^\varepsilon D(u^\varepsilon)n^\varepsilon] \cdot \tau^\varepsilon = 0$ ).
Forçage capillaire : La force de tension de surface est modulée par un paramètre $\lambda_\varepsilon > 0$ qui dépend de l'échelle microscopique $\varepsilon$ .

L'objectif est double : prouver l'existence de solutions faibles pour chaque $\varepsilon$ fixé, puis effectuer l'homogénéisation lorsque $\varepsilon \to 0$ pour obtenir des modèles macroscopiques effectifs.

2. Méthodologie

L'analyse se déroule en deux grandes étapes :

A. Existence et Estimations A Priori (Fixe $\varepsilon$ )

Pour un $\varepsilon > 0$ fixé, les auteurs établissent l'existence de solutions faibles sur un intervalle de temps fini $(0, T)$ .

Schéma de Galerkin : La preuve d'existence utilise un schéma de Galerkin basé sur des fonctions propres de type Stokes adaptées aux conditions aux limites de glissement libre.
Estimations uniformes : Des estimations a priori sont dérivées pour les champs de vitesse $u^\varepsilon$ $u^{ε}$ , de phase $\phi^\varepsilon$ $ϕ^{ε}$ , de potentiel chimique $\mu^\varepsilon$ $μ^{ε}$ et de pression $p^\varepsilon$ $p^{ε}$ .
- Une difficulté majeure réside dans le terme source $G(\phi^\varepsilon)$ qui empêche l'utilisation directe de l'inégalité de Poincaré-Wirtinger classique (car la moyenne de $\phi^\varepsilon$ n'est pas nulle). Les auteurs contournent cela en estimant explicitement la moyenne spatiale $\bar{\phi}^\varepsilon$ via l'inégalité de Gronwall.
- Les estimations sont uniformes par rapport à $\varepsilon$ et à $\lambda_\varepsilon$ , en exploitant la loi de dissipation d'énergie du système.
Outils techniques : Utilisation d'opérateurs de Bogovskiĭ (pour inverser la divergence sur les domaines perforés), d'opérateurs de restriction et d'extension (pour passer du domaine perforé au domaine global), et d'inégalités de Korn et Poincaré avec des constantes indépendantes de $\varepsilon$ .

B. Homogénéisation (Limite $\varepsilon \to 0$ )

L'homogénéisation est réalisée en utilisant la méthode de dépliement périodique (periodic unfolding operator) $T^*_\varepsilon$ , adaptée aux domaines perforés.

Échelle du forçage : Le paramètre de capillarité $\lambda_\varepsilon$ converge vers une limite $\lambda \in [0, +\infty)$ . La force externe $g^\varepsilon$ est supposée être de l'ordre de $\sqrt{\lambda_\varepsilon}$ pour maintenir un équilibre non trivial dans l'équation de la quantité de mouvement.
Deux régimes distincts : L'analyse distingue deux cas selon la valeur de la limite $\lambda$ $λ$ :
1. Régime de capillarité nulle ( $\lambda = 0$ ) : La force capillaire devient négligeable à l'échelle macroscopique.
2. Régime équilibré ( $\lambda \in (0, +\infty)$ ) : La force capillaire et les effets convectifs sont du même ordre de grandeur.

3. Résultats Clés

A. Résultats d'Existence

Existence d'une solution faible $(u^\varepsilon, \phi^\varepsilon, \mu^\varepsilon)$ pour le système NSCH microscopique avec terme source et conditions de glissement libre.
Démonstration de la loi de dissipation d'énergie et de l'existence d'une pression associée $p^\varepsilon$ .

B. Modèles Homogénéisés Effectifs

Selon la limite $\lambda$ , deux modèles macroscopiques distincts émergent :

1. Cas $\lambda = 0$ (Régime de Stokes-Cahn-Hilliard) :

Le terme convectif $(u \cdot \nabla)u$ et le terme d'advection $u \cdot \nabla \phi$ disparaissent à l'échelle macroscopique.
Le système limite est de type Stokes-Cahn-Hilliard (SCH) :
- Équation de quantité de mouvement : Stokes (linéaire) avec un forçage capillaire homogénéisé.
- Équation de phase : Diffusion pure avec source, sans advection macroscopique.
Cela signifie que dans ce régime, les effets d'interface sont trop faibles pour générer un transport macroscopique significatif par convection.

2. Cas $\lambda \in (0, +\infty)$ (Régime Navier-Stokes-Cahn-Hilliard) :

Le système limite conserve la structure non-linéaire complète.
Le système est de type Navier-Stokes-Cahn-Hilliard (NSCH) :
- L'équation de quantité de mouvement contient un terme convectif non linéaire $(u \cdot \nabla)u$ .
- L'équation de Cahn-Hilliard contient un terme d'advection $u \cdot \nabla \phi$ .
L'intensité de ces termes non linéaires est proportionnelle à $\sqrt{\lambda}$ . Cela indique que la capillarité compète avec la viscosité et que le mouvement de l'interface peut induire des écoulements macroscopiques appréciables.

C. Tenseurs Effectifs

Dans les deux cas, les coefficients microscopiques oscillants sont remplacés par des tenseurs effectifs homogénéisés ( $A^{hom}$ , $B^{hom}$ , $C^{hom}$ ) définis via des problèmes cellulaires (cell problems) sur la cellule de référence $Y_p$ . Ces tenseurs capturent l'effet de la géométrie du pore et des conditions aux limites de glissement.

D. Convergence de l'Énergie

Les auteurs prouvent la convergence de l'énergie totale microscopique (normalisée) vers l'énergie du système macroscopique, satisfaisant une loi de dissipation analogue.

4. Contributions et Signification

Nouveauté Mathématique : C'est l'une des premières études rigoureuses combinant l'existence de solutions pour un système NSCH avec un terme source non conservatif et des conditions de glissement libre sur un domaine perforé, suivie d'une procédure d'homogénéisation complète.
Gestion du Terme Source : La méthode développée pour traiter la non-conservation de la masse (estimation de la moyenne spatiale) est une contribution technique importante, car elle permet d'éviter les hypothèses restrictives usuelles sur la moyenne nulle de la concentration.
Impact Physique : Le travail clarifie comment les conditions aux limites de glissement (friction nulle sur les inclusions) et la force capillaire influencent le comportement macroscopique. Il montre une transition de phase entre un régime de Stokes (dominé par la diffusion) et un régime de Navier-Stokes (dominé par la convection) selon l'intensité relative de la capillarité.
Applications : Ces résultats sont pertinents pour la modélisation de processus complexes comme la récupération assistée du pétrole, la dépollution des eaux souterraines, ou la dynamique des hydrogels et des tissus biologiques, où les interfaces diffuses et les milieux poreux jouent un rôle central.

En résumé, cet article fournit un cadre analytique rigoureux pour dériver des modèles macroscopiques fiables à partir de descriptions microscopiques complexes de mélanges biphasiques en milieux poreux, en tenant compte de la dynamique des interfaces et des conditions de glissement réalistes.