Singularities in phase separation models: a spectral element approach for the nonlocal Cahn-Hilliard equation

Titre

Singularités dans les modèles de séparation de phases : une approche par éléments spectraux pour l'équation de Cahn–Hilliard non locale.

1. Problématique et Contexte

L'équation de Cahn–Hilliard (CH) classique est un modèle macroscopique décrivant la séparation de phases dans les alliages binaires. Ce papier se concentre sur sa version non locale, qui intègre des interactions à longue portée via un noyau de convolution singulier.

Le système étudié est défini par :
$$\frac{\partial \rho}{\partial t} = \nabla \cdot \left( m(\rho) \nabla \mu \right), \quad \mu = F'(\rho) - K * \rho$$
où :

• $\rho$ est la densité de phase.
• $m(\rho)$ est la mobilité (constante dans ce travail).
• $F(\rho)$ est un potentiel logarithmique singulier (dérivé de l'énergie libre), qui garantit que la solution reste physiquement dans l'intervalle $(-1, 1)$, évitant ainsi les états purs $\pm 1$.
• $K$ est un noyau d'interaction de type Newtonien ($K(x) \sim \log \|x\|$ en 2D), qui présente une singularité faible à l'origine.

Défis numériques :
L'approximation numérique de ce système est difficile pour deux raisons principales :

1. La nature non locale du terme de convolution $K * \rho$ conduit à des opérateurs denses, coûteux en calcul.
2. La singularité du noyau $K$ à l'origine et la singularité du potentiel $F'$ aux bords de l'intervalle $[-1, 1]$ posent des problèmes de stabilité et de précision pour les schémas de collocation standards.

2. Méthodologie

Les auteurs développent un cadre numérique efficace et contrôlé par l'erreur basé sur une méthode pseudo-spectrale multiforme (multishape).

A. Décomposition de l'intégrale singulière

Pour traiter la singularité du noyau Newtonien, l'intégrale de convolution est décomposée en deux parties pour un rayon de coupure $\varepsilon$ :
$$(K * \rho)(x) = \int_{\Omega \setminus B_\varepsilon(x)} K(x-y)\rho(y)dy + \int_{\Omega \cap B_\varepsilon(x)} K(x-y)\rho(y)dy$$

• Partie régulière ($I_1$) : L'intégrale sur le domaine extérieur à la boule $\varepsilon$ est calculée numériquement.
• Partie singulière ($I_2$) : L'intégrale sur la petite région autour de $x$ est approximée par un développement de Taylor d'ordre zéro (approximation par morceaux constants) :
  $$I_2 \approx \rho(x) \int_{x - (\Omega \cap B_\varepsilon(x))} K(u) du$$
  Les auteurs dérivent des formules fermées exactes pour ces intégrales géométriques, ce qui permet de calculer le terme de correction analytiquement sans erreur de quadrature.

B. Méthode des Éléments Spectraux (SEM) et Multiforme

Pour évaluer la partie régulière $I_1$ avec une haute précision :

• Ils utilisent la bibliothèque logicielle MultiShape, qui implémente une méthode d'éléments spectraux basée sur la méthode de "patching" (assemblage).
• Le domaine de calcul (le complément de la petite boule autour du point de collocation) est discrétisé en un assemblage de quadrilatères (un "multiforme").
• Cela permet d'utiliser des quadratures de haute précision (Clenshaw–Curtis) sur des géométries complexes générées dynamiquement autour de chaque point de collocation.
• Une matrice de convolution précalculée est construite, permettant une évaluation rapide de l'opérateur non local.

C. Analyse Mathématique

Le papier étend la théorie d'existence et d'unicité aux domaines de type Lipschitz (avec coins), ce qui est crucial pour les maillages numériques standards. Ils prouvent également :

• La propriété de séparation stricte instantanée : la solution s'éloigne uniformément des phases pures $\pm 1$ pour tout $t > 0$.
• La convergence vers un équilibre unique avec un taux de décroissance polynomial en temps, grâce à l'inégalité de Łojasiewicz–Simon.

3. Résultats Clés

A. Validation de la méthode spectrale

• Les tests numériques montrent une convergence spectrale pour l'approximation du potentiel Newtonien.
• L'erreur d'approximation de la partie singulière est de l'ordre de $O(\varepsilon^2 \log \varepsilon)$ et est indépendante du maillage.
• La méthode permet d'atteindre des résolutions élevées avec des ressources computationnelles limitées (les matrices de convolution sont précalculées une seule fois).

B. Comportement dynamique du système

Les auteurs ont simulé le système avec différentes conditions initiales et paramètres de mise à l'échelle du noyau ($\eta$) :

1. Diffusion pure : Pour des valeurs positives de $\eta$, le système tend vers un état homogène.
2. Séparation de phases : Pour des valeurs négatives de $\eta$, le système subit une séparation de phases, formant des structures stables.
3. Effet du potentiel régularisé : L'utilisation d'un potentiel régularisé ($C^3$) permet d'accélérer considérablement les simulations à long terme en évitant la rigidité (stiffness) introduite par la singularité logarithmique, tout en conservant une précision proche de la machine pour les états stationnaires.
4. Mélanges de noyaux : En combinant le noyau Newtonien avec un mollifieur, les auteurs ont observé la formation de structures riches et complexes (motifs énergétiques) qui ne sont pas présentes avec le noyau Newtonien seul.

4. Contributions Principales

1. Extension théorique : Preuve de l'existence, de l'unicité et de la convergence vers l'équilibre pour le système CH non local avec potentiel singulier sur des domaines de Lipschitz (y compris les polygones).
2. Nouvelle méthode numérique : Développement d'une méthode d'éléments spectraux adaptée aux noyaux singuliers via une décomposition de coupure et une approximation géométrique exacte, une approche nouvelle dans la littérature pour ce type de problème.
3. Analyse d'erreur rigoureuse : Établissement d'estimations d'erreur analytiques et vérification numérique de la convergence.
4. Résultats numériques à haute résolution : Première présentation de solutions numériques haute résolution pour le système CH non local avec potentiel logarithmique et noyau Newtonien, révélant des dynamiques complexes (séparation de phases, motifs).
5. Logiciel open-source : Mise à disposition du code source pour la reproduction des résultats.

5. Signification et Perspectives

Ce travail comble un vide important dans la littérature numérique concernant les modèles de dynamique de densité fonctionnelle (DDFT) non locaux avec des potentiels singuliers. La méthode proposée surmonte les limitations des méthodes spectrales classiques face aux singularités de noyaux.

Les résultats ouvrent la voie à des applications plus complexes, notamment :

• La résolution de systèmes couplés (ex: Cahn–Hilliard–Navier–Stokes).
• L'étude de problèmes de contrôle optimal pour des modèles biologiques (comme le modèle Hawkins–Daruud pour la croissance tumorale).
• La modélisation de réactions chimiques médiées par des interactions non locales.

En résumé, l'article fournit un cadre robuste, mathématiquement justifié et numériquement efficace pour simuler des phénomènes de séparation de phases complexes régis par des interactions non locales singulières.