On the convergence of explicit formulas for $L^2$ solutions to the Benjamin-Ono and continuum Calogero-Moser equations

1. Problématique et Contexte

L'article s'intéresse à l'approximation numérique de deux équations aux dérivées partielles (EDP) intégrables non locales définies sur le tore $\mathbb{T}$ :

1. L'équation de Benjamin-Ono (BO) : $\partial_t u = \partial_x (|\partial_x|u - u^2)$, modélisant les ondes internes dans un fluide stratifié.
2. L'équation de Calogero-Moser continue (CCM) : $i\partial_t u + \partial_x^2 u \pm \frac{2}{i} u \partial_x \Pi(|u|^2) = 0$, liée au système de particules de Calogero-Moser et à la modulation d'ondes (également connue sous le nom d'équation DNLS de Calogero-Sutherland).

Le défi principal :
La plupart des schémas numériques classiques pour ces équations nécessitent une régularité initiale élevée (souvent $H^s$ avec $s > 1$) pour garantir la convergence, ou ne convergent que dans des normes faibles (ex: $L^2([0,T], L^2_{loc})$) pour des données initiales peu régulières ($L^2$). Or, la régularité critique pour le bien-posé global de ces équations est précisément $L^2(\mathbb{T})$ (pour BO) et $L^2_+(\mathbb{T})$ (pour CCM, espace de Hardy). De plus, les schémas classiques souffrent souvent d'un coût computationnel qui dépend du temps final $T$ ou d'une dérive de l'erreur sur de longs intervalles de temps.

L'objectif est de concevoir et d'analyser des schémas numériques fully-discrete (discrétisés en temps et en espace) qui :

• Convergent pour des données initiales grossières ($u_0 \in L^2$).
• Préservent les structures géométriques (conservation de la masse, de la norme $L^2$).
• Soient exacts en temps (coût indépendant de $T$).

2. Méthodologie

Au lieu de discrétiser directement les équations différentielles, les auteurs exploitent la structure intégrable de ces systèmes via des formules explicites récentes découvertes par Gérard (pour BO) et Badreddine (pour CCM).

A. Utilisation des Formules Explicites

Les solutions $u(t)$ sont caractérisées par leurs coefficients de Fourier via des opérateurs de Lax :

• Pour BO : $\hat{u}(t, k) = \langle (e^{it(I+2L_{u_0})}S^*)^k \Pi u_0, 1 \rangle$
• Pour CCM : $\hat{u}(t, k) = \langle (e^{-it(I+2L_{u_0})}S^*)^k u_0, 1 \rangle$

Où :

• $L_{u_0}$ est l'opérateur de Lax (semi-borné, auto-adjoint).
• $S^*$ est l'opérateur de décalage (shift) dans l'espace de Hardy.
• $\Pi$ est le projecteur de Riesz-Szegő.

B. Conception du Schéma Numérique

Les auteurs généralisent le schéma précédent de [3] en introduisant une séquence de paramètres de troncature $n_K(k)$ dépendant de l'itération $k$.
Le schéma $u_K(t)$ est défini par :
$$\hat{u}_K(t, k) = \langle u_k, 1 \rangle, \quad \text{où } u_k = (e^{\pm it(I+2L_{n(k)})}S^*) u_{k-1}$$
avec $L_{n(k)}$ étant l'opérateur de Lax tronqué à la fréquence $n(k)$.

Innovation clé : La liberté dans le choix de la séquence $n_K(k)$ permet de concevoir des schémas qui préservent exactement la norme $L^2$ au niveau discret (en choisissant $n(k) \le K-k$), contrairement aux troncatures simples.

C. Analyse de Convergence

La preuve de convergence repose sur l'analyse spectrale des opérateurs de Lax tronqués $L_n$ :

1. Convergence de la résolvante : Démonstration que $L_n$ converge vers $L$ au sens de la résolvante en norme ($\|R_n(\kappa) - R(\kappa)\|_{op} \to 0$) pour $u_0 \in L^2$. Cela utilise des estimations de perturbations relatives de l'opérateur libre $-i\partial_x$.
2. Convergence des groupes unitaires : Utilisation du théorème de Stone-Weierstrass et de la densité pour montrer que les groupes unitaires $e^{itL_n}$ convergent uniformément vers $e^{itL}$ sur des intervalles de temps bornés et des compacts de $L^2$.
3. Passage de la convergence faible à forte :
   • On établit d'abord la convergence faible des coefficients de Fourier.
   • On utilise ensuite la conservation de la norme $L^2$ (ou sa borne uniforme) du schéma numérique (Proposition 2.4) pour élever la convergence faible en convergence forte dans $L^2$.

3. Résultats Principaux

Théorème de Convergence (Théorème 1.3)

Pour toute donnée initiale $u_0 \in L^2(\mathbb{T})$ (pour BO) ou $u_0 \in L^2_+(\mathbb{T})$ (pour CCM, avec $\|u_0\|_{L^2} < 1$ dans le cas focalisant), le schéma numérique $u_K(t)$ converge vers la solution exacte $u(t)$ dans la norme $L^2(\mathbb{T})$ :
$$\sup_{t \in [-T, T]} \|u_K(t) - u(t)\|_{L^2} \xrightarrow[K \to \infty]{} 0$$
Cette convergence est uniforme sur tout intervalle de temps borné $[-T, T]$.

Points forts de ce résultat :

• C'est le premier résultat d'erreur pour un schéma numérique convergeant dans $C([0, T], L^2)$ pour des données initiales purement $L^2$ (sans régularité Sobolev positive).
• Le coût computationnel est indépendant du temps final $T$ (les schémas sont exacts en temps, seule la discrétisation spatiale introduit l'erreur).
• La méthode ne nécessite pas d'estimations de régularité supplémentaires (comme les espaces de Bourgain) qui échouent souvent pour $s=0$.

Propriétés de Préservation de Structure

• Conservation de la masse (moyenne) : Le schéma préserve exactement $\int u_K dx = \int u_0 dx$.
• Préservation de la norme $L^2$ : En choisissant judicieusement la séquence de troncature $n(k)$ (ex: $n(k) \le K-k$), le schéma préserve exactement la norme $L^2$ de la donnée initiale tronquée, offrant une stabilité numérique supérieure.

Application : L'Effet Talbot

Les auteurs appliquent leur schéma à l'équation BO avec une condition initiale en forme d'onde carrée ($u_0(x) = \text{sgn}(x)$), qui appartient à $H^{1/2-}$ mais pas à $L^2_+$.

• Ils observent numériquement l'effet Talbot (reconstruction périodique de l'onde à des temps rationnels et comportement fractal/continu à des temps irrationnels) pour l'équation non linéaire.
• Cela fournit la première preuve numérique rigoureuse de l'effet Talbot pour BO, étayant l'hypothèse que la dynamique non linéaire préserve les singularités de la dynamique linéaire, bien que la preuve analytique rigoureuse de cet effet pour BO reste un problème ouvert.

4. Contributions Clés

1. Extension aux données $L^2$ : Généralisation des résultats de convergence aux régularités critiques, là où les méthodes classiques échouent.
2. Nouvelle classe de schémas : Introduction d'une famille de schémas paramétrés par $n_K(k)$ permettant de contrôler les propriétés de conservation (notamment la norme $L^2$) au niveau discret.
3. Méthodologie d'analyse : Développement d'une analyse de convergence basée sur la théorie des opérateurs (convergence de la résolvante, groupes unitaires) plutôt que sur des estimations d'énergie classiques ou des espaces de Bourgain discrets.
4. Efficacité computationnelle : Démonstration que l'exploitation des formules explicites permet des simulations à long terme avec une erreur contrôlée et un coût indépendant de $T$.

5. Signification et Impact

Ce travail représente une avancée majeure dans le domaine de l'analyse numérique des EDP intégrables non linéaires. Il démontre que l'exploitation directe de la structure intégrable (via les formules explicites et les opérateurs de Lax) offre une voie supérieure aux méthodes de discrétisation standard pour les problèmes de faible régularité.

• Théorique : Il comble un vide important en établissant la convergence forte en $L^2$ pour des données critiques, validant l'utilisation de ces schémas pour des solutions qui ne sont pas lisses.
• Pratique : Il permet l'étude numérique fiable de phénomènes complexes (comme l'effet Talbot ou les solutions avec singularités) sur des temps très longs, ce qui était auparavant prohibitif ou peu fiable avec les schémas classiques.
• Généralité : La méthode ouvre la porte à l'application de techniques similaires à d'autres équations intégrables possédant des formules explicites, bien que l'analyse spécifique dépende fortement de la structure des opérateurs de Lax (ce qui exclut, par exemple, l'équation de Szegő cubique dans ce cadre précis).