A splitting scheme for the wave maps equation at low regularity

Résumé technique : Un schéma de fractionnement pour l'équation des cartes d'ondes à faible régularité

Énoncé du problème
Ce travail traite de l'approximation numérique de l'équation des cartes d'ondes tridimensionnelle à valeurs dans la sphère $S^2$. L'équation, donnée par $(\partial_t^2 - \Delta)u = u(|\partial_t u|^2 - |\nabla_x u|^2)$, est une équation d'ondes géométrique fondamentale. Bien que le bien-posé local de cette équation dans le cadre continu soit bien compris pour des données initiales dans $H^s$ avec $s > d/2$ (spécifiquement $s > 3/2$ en 3D) grâce à l'exploitation de la « structure nulle » de la non-linéarité, l'analyse numérique de tels schémas à ce niveau de faible régularité reste un défi.

Les méthodes numériques standard, telles que les approches par éléments finis ou les schémas de fractionnement de Lie, reposent généralement sur des méthodes énergétiques qui exigent une régularité supérieure ($s > d/2 + 1$) pour prouver la convergence. Étendre ces preuves de convergence au régime sous-critique ($s > d/2$) est difficile car :

1. La structure nulle (annulations entre ondes parallèles) est essentielle au bien-posé mais n'est généralement pas préservée par les discrétisations standard.
2. La non-linéarité implique des dérivées temporelles, qui introduisent des sources d'erreur significatives dans les schémas de discrétisation temporelle.
3. Les méthodes existantes qui préservent la contrainte sphérique (par exemple, les méthodes de projection) exigent souvent une grande régularité pour contrôler l'erreur de projection, les rendant incompatibles avec une analyse de faible régularité.

Méthodologie
Les auteurs proposent et analysent un schéma de fractionnement de Lie filtré conçu spécifiquement pour préserver la structure nulle au niveau discret. La méthodologie comprend plusieurs composants clés :

• Reformulation : Avant la discrétisation, la non-linéarité est réécrite en utilisant l'identité $|\partial_t u|^2 - |\nabla u|^2 = \Box(u \cdot u) - 2u \cdot \Box u$. Cette reformulation isole l'opérateur d'ondes agissant sur le produit, permettant au schéma de mieux approximer la structure nulle.
• Espaces de Bourgain discrets : L'analyse de convergence est menée dans des espaces de Bourgain discrets ($X^{s,b}_\tau$), qui sont les analogues discrets des espaces continus utilisés dans la théorie du bien-posé. Ces espaces utilisent la transformée de Fourier spatio-temporelle et des poids spécifiques impliquant le symbole de dérivée temporelle discret $d_\tau(\sigma) = (e^{i\sigma\tau}-1)/\tau$.
• Filtrage fréquentiel : Le schéma emploie un filtre fréquentiel $\Pi$ (un lissage à haute fréquence $\sim \tau^{-1/2}$) pour traiter les dérivées temporelles discrètes dans la non-linéarité. Bien que cette projection vers des fréquences plus basses puisse sembler restrictive, elle est nécessaire pour approximer efficacement les dérivées discrètes.
• Fractionnement de Lie : Le schéma sépare l'évolution linéaire des ondes ($L_\tau$) et le flot non linéaire. L'étape non linéaire est approchée en utilisant une étape d'Euler explicite avec la non-linéarité filtrée $N_\tau$.
• Estimations multilinéaires : Une réalisation technique centrale est le développement d'estimations bilinéaires et trilinéaires dans ces espaces discrets. Les auteurs analysent soigneusement la géométrie des interactions d'ondes dans l'espace de Fourier discret (en considérant la périodicité et la modulation) pour prouver que la structure nulle discrète produit les annulations nécessaires. Ils établissent des estimations de Strichartz pour l'opérateur d'ondes discret (s'appuyant sur les travaux de Ruff et Schnaubelt) et les utilisent pour contrôler les termes non linéaires.

Contributions clés

1. Préservation de la structure nulle discrète : L'article démontre qu'en reformulant la non-linéarité et en utilisant un fractionnement de Lie filtré spécifique, les annulations cruciales de la structure nulle peuvent être préservées au niveau discret, même en présence de dérivées temporelles.
2. Cadre de Bourgain discret : Les auteurs définissent rigoureusement les espaces de Bourgain discrets pour l'équation d'ondes et établissent les estimations linéaires et multilinéaires nécessaires (incluant les estimations de Strichartz et bilinéaires) requises pour l'analyse de faible régularité. Cela étend le cadre précédemment développé pour les équations de Schrödinger au contexte de l'équation d'ondes.
3. Preuve de convergence : Le résultat principal (Théorème 1.2) prouve la convergence du schéma pour des données initiales dans $H^s \times H^{s-1}$ pour tout $s > 3/2$ (la plage sous-critique complète en 3D). L'estimation d'erreur est donnée par :
   $$\|u(t_n) - u^n\|_{H^{s_1} \times H^{s_1-1}} \leq C \tau^{(s-s_1)/2}$$
   pour tout $3/2 < s_1 < s < 2$.
4. Traitement des dérivées temporelles : L'analyse aborde explicitement l'erreur introduite par les dérivées temporelles dans la non-linéarité, montrant que bien que cela soit la source d'erreur dominante, cela peut être contrôlé sous une condition CFL forte et la structure spécifique du schéma.

Résultats

• Théoriques : Le schéma converge vers la solution forte (Duhamel) pour toutes les données initiales sous-critiques. Le taux de convergence dépend de l'écart de régularité entre l'espace de la solution ($s$) et l'espace de mesure de l'erreur ($s_1$).
• Expériences numériques : Les auteurs présentent des simulations numériques sur un tore 3D.
  • Pour des données lisses ($s=3$), le schéma présente le taux de convergence prédit de $O(\tau)$ dans la norme $L^2 \times H^{-1}$.
  • Pour des données irrégulières ($H^{1.7} \times H^{0.7}$), le schéma converge dans $H^{1.6} \times H^{0.6}$ avec un taux de $\tau^{0.05}$, ce qui est plus lent que le cas lisse mais cohérent avec la prédiction théorique pour une faible régularité. La décroissance observée est notée comme étant plus forte que la borne inférieure théorique mais significativement plus lente que le cas lisse.

Signification et affirmations
L'article prétend fournir la première analyse de convergence d'un schéma de fractionnement pour l'équation des cartes d'ondes au seuil de régularité minimal ($s > d/2$). Les auteurs soulignent que la difficulté principale réside dans la préservation de la structure nulle et la gestion des dérivées temporelles dans la discrétisation.

La signification de ce travail est double :

1. Il comble le fossé entre la théorie continue du bien-posé à faible régularité (qui repose sur les structures nulles et les espaces de Bourgain) et l'analyse numérique, démontrant que ces outils analytiques sophistiqués peuvent être adaptés aux contextes discrets.
2. Il offre une méthode numérique viable pour les cartes d'ondes qui ne nécessite pas les hypothèses de grande régularité des méthodes basées sur la projection, permettant ainsi la simulation de données initiales plus irrégulières.

Les auteurs notent que bien que leur schéma ne préserve pas strictement la contrainte sphérique $u(t,x) \in S^2$ à chaque pas de temps, la convergence dans $H^s$ ($s > 3/2$) implique que les solutions s'approchent de la sphère uniformément lorsque le pas de temps $\tau \to 0$. Ils suggèrent également que leurs arguments peuvent probablement être généralisés aux dimensions $d \geq 3$ et potentiellement à $d \geq 2$ avec des modifications appropriées.