An energy-decreasing algorithm for the finite element approximation of ferronematic equilibrium states

Résumé technique : Un algorithme à énergie décroissante pour l'approximation par éléments finis des états d'équilibre ferronématiques

Formulation du problème
Ce travail traite de l'approximation numérique des états d'équilibre bidimensionnels des ferronématiques. Les ferronématiques sont des systèmes composites constitués de nanoparticules magnétiques dispersées dans un cristal liquide nématique. Les états d'équilibre sont modélisés comme des minimiseurs d'une densité d'énergie libre. Dans la limite asymptotique spécifique considérée ici, où le rapport entre l'échelle de longueur dépendante du matériau et l'échelle de longueur du domaine de calcul est très faible, le problème se réduit à la minimisation de l'énergie harmonique de deux champs vectoriels couplés : le paramètre d'ordre nématique $\mathbf{Q}$ et le vecteur d'aimantation $\mathbf{M}$.

Les deux champs sont soumis à des contraintes de longueur prescrites et à une relation d'orientation non linéaire. Plus précisément, le vecteur d'état $\Psi = (\mathbf{Q}, \mathbf{M})^T$ est contraint de se situer dans une sous-variété lisse, compacte et connexe $\mathcal{N}$ de dimension un de $\mathbb{R}^4$. Le problème mathématique est formulé comme la minimisation du fonctionnel d'énergie harmonique :
$$E(\Psi) := \frac{1}{2} \int_{\Omega} \|\nabla \Psi\|^2_{\mathbb{R}^4 \times 2} \, dx$$
sous la contrainte $\Psi(\mathbf{x}) \in \mathcal{N}$ presque partout dans le domaine $\Omega$, avec des conditions aux limites de Dirichlet prescrites compatibles avec la contrainte.

Méthodologie
Les auteurs proposent une approximation par éléments finis utilisant des éléments affines par morceaux et continus sur une triangulation faiblement acutangle. Un aspect clé de la méthodologie réside dans la gestion de la variété de contrainte non convexe $\mathcal{N}$.

1. Algorithme au niveau continu : Les auteurs conçoivent d'abord un algorithme itératif à énergie décroissante au niveau continu. Cela repose sur une représentation de la variété $\mathcal{N}$ via un vecteur unitaire $\mathbf{n}$ issu du cercle unité $\mathbb{S}^1$. L'algorithme consiste à projeter une mise à jour linéaire sur la sphère afin de garantir que la contrainte soit satisfaite. L'analyse théorique démontre que cette séquence de mises à jour diminue strictement le fonctionnel d'énergie.

2. Algorithme au niveau discret : L'extension de ceci au niveau discret n'est pas triviale car l'identité continue reliant les gradients de $\Psi$ et $\mathbf{n}$ ne vaut pas pour les champs discrets. Les auteurs introduisent deux innovations clés :

   • Propriété de maillage faiblement acutangle : Ils exploitent la propriété de faiblement acutangle de la triangulation pour exprimer l'énergie harmonique discrète comme une somme de contributions locales positives indexées par des paires de nœuds du maillage.
   • Lemme géométrique : Une propriété géométrique cruciale de la variété $\mathcal{N}$ (Lemme 3.2) est utilisée pour relier les mises à jour des deux vecteurs unitaires associés aux nœuds. Cela permet de construire une itération à énergie décroissante discrète.

3. Réalisation computationnelle : Pour résoudre efficacement le problème de minimisation sous contrainte résultant, les auteurs emploient un cadre de décomposition-coordination basé sur une approche de Lagrangien augmenté.

   • La contrainte non linéaire est découplée en introduisant une variable auxiliaire et un multiplicateur de Lagrange.
   • Une méthode itérative de type Uzawa (spécifiquement, une méthode de direction alternée des multiplicateurs ou ADMM) est utilisée pour trouver le point selle du Lagrangien augmenté.
   • Pour améliorer l'efficacité computationnelle, la fonction de couplage non linéaire dans l'étape de mise à jour primale est linéarisée à chaque itération, transformant le sous-problème en un système linéaire résoluble par des méthodes directes (décomposition QR).

Contributions clés

• Développement d'algorithmes : L'article présente un nouvel algorithme itératif à énergie décroissante spécifiquement adapté au problème d'équilibre ferronématique sous contraintes de longueur.
• Garanties théoriques : Les auteurs fournissent une preuve rigoureuse que l'algorithme discret proposé garantit une séquence de valeurs d'énergie non croissante, à condition que la triangulation soit faiblement acutangle. Ceci est réalisé grâce à un lemme géométrique spécifique concernant la variété de contrainte.
• Approximation non conforme : La méthode utilise une approximation non conforme où la contrainte est imposée strictement aux nœuds du maillage, tandis que les fonctions discrètes ne prennent leurs valeurs dans la variété qu'à ces nœuds.
• Résolveur efficace : L'intégration de la stratégie de décomposition-coordination avec un schéma Uzawa/ADMM linéarisé offre une approche computationnellement viable pour résoudre le système non linéaire couplé.

Résultats numériques
L'algorithme est testé sur deux exemples impliquant des conditions aux limites de Dirichlet lisses sur un domaine carré, avec des solutions analytiques connues.

• Convergence : Les expériences numériques démontrent que l'erreur d'énergie diminue lorsque le maillage est affiné. Pour le premier cas test, une convergence d'ordre un est observée à la fois dans la semi-norme $H^1$ et la norme $L^2$. Pour le second cas test, plus complexe, une convergence d'ordre un est observée dans la semi-norme $H^1$, tandis qu'une convergence d'ordre deux est atteinte dans la norme $L^2$.
• Sensibilité aux paramètres : L'étude analyse l'influence du paramètre d'augmentation ($\zeta$), du paramètre de mise à jour ($\rho$) et de la tolérance Uzawa ($\varepsilon_{pri}$). Les résultats indiquent que, bien que des valeurs plus élevées de $\zeta$ augmentent considérablement le nombre d'itérations, elles n'améliorent pas nécessairement le rapport précision-coût. Des tolérances plus strictes sont requises sur des maillages plus fins pour prévenir les oscillations dans la décroissance de l'énergie de la boucle externe.
• Performance : La méthode converge avec succès vers des états d'équilibre avec des erreurs d'énergie contrôlées, validant ainsi la propriété théorique de décroissance de l'énergie.

Signification
L'article prétend fournir la première analyse numérique et cadre algorithmique spécifiquement pour le problème d'équilibre ferronématique dans la limite asymptotique décrite. Bien que la minimisation de l'énergie harmonique sous contraintes de longueur soit bien explorée pour des modèles plus simples (tels que le modèle d'Oseen–Frank ou le modèle d'Eriksen), ce travail étend ces concepts au système ferronématique couplé. La contribution principale est le développement d'un algorithme robuste par éléments finis à énergie décroissante qui gère le couplage non linéaire spécifique et les contraintes géométriques des ferronématiques, offrant un outil pratique pour la simulation de ces matériaux dans les applications d'optique, de capteurs et de télécommunications.