Stochastic numerical approximation for nonlinear Fokker-Planck equations with singular kernels

Résumé technique : Approximation numérique stochastique pour les équations de Fokker-Planck non linéaires à noyaux singuliers

Énoncé du problème
Cet article traite de l'approximation numérique des solutions d'équations aux dérivées partielles (EDP) de Fokker-Planck non linéaires de la forme :
$$\partial_t u = \Delta u - \nabla \cdot (u K * u)$$
où le noyau d'interaction $K$ est singulier. De telles équations apparaissent dans divers contextes physiques et biologiques, notamment le modèle de chimiotactisme de Keller-Segel, les potentiels de Coulomb en électrostatique/gravité, et la loi de Biot-Savart pour les équations de Navier-Stokes 2D.

Le défi principal réside dans le terme de dérive $K * u$, qui implique un noyau singulier. Les approches numériques standards, telles que l'application directe du schéma d'Euler-Maruyama (EM) à l'équation différentielle stochastique (EDS) de McKean-Vlasov associée, font face à deux difficultés :

1. Calcul de la densité : Le schéma EM nécessite d'évaluer la densité de la solution à chaque pas de temps, ce qui n'est pas directement disponible.
2. Singularité et bien-fondé : À mesure que le pas de temps $h \to 0$, les singularités du noyau d'interaction peuvent compromettre le bien-fondé du système de particules résultant.

Pour surmonter ces obstacles, l'auteur considère un système de $N$ particules en interaction approximant la solution via une mesure empirique. Cependant, en raison de la dérive singulière, le système est régularisé à l'aide d'une séquence de mollificateurs $(V^N)$ et d'une fonction de coupure lisse $F_A$ pour borner la dérive, garantissant ainsi la stabilité numérique et la bornitude théorique.

Méthodologie
L'article analyse la convergence du système de particules régularisé et discrétisé (Équation 1.4/2.2) vers la solution de l'équation de Fokker-Planck non linéaire. La méthodologie comprend :

1. Régularisation et discrétisation :

   • Le noyau singulier $K$ est convolué avec un mollificateur $V^N$ (paramètre d'échelle $\alpha$).
   • Une fonction de coupure $F_A$ est appliquée à la dérive pour assurer une bornitude uniforme, cohérente avec la borne $L^\infty$ de la dérive réelle $K * u_t$.
   • Le système est discrétisé en utilisant le schéma d'Euler-Maruyama avec un pas de temps $h$.

2. Stratégie d'analyse de convergence :
   L'analyse décompose l'erreur entre la mesure empirique mollifiée $\tilde{\mu}^{N,h}_t$ (et la densité d'une particule unique) et la solution réelle $u_t$ en plusieurs composantes :

   • Erreur initiale : La différence entre la mesure empirique mollifiée initiale et $u_0$.
   • Erreur de propagation du chaos : La différence entre le système de particules régularisé et l'EDS de McKean-Vlasov (non discrétisée).
   • Erreur de discrétisation : L'erreur introduite par le schéma EM (discrétisation temporelle).
   • Convolution stochastique : Un terme provenant du bruit dans le système de particules.

3. Outils techniques :

   • Approche par semi-groupes : L'auteur utilise la formulation faible de l'équation de Fokker-Planck impliquant le noyau de la chaleur $g_2(t, \cdot)$.
   • Formule d'Itô : Appliquée pour développer les densités du système de particules et de la solution de McKean-Vlasov autour du noyau de la chaleur.
   • Bornes de moments : L'auteur établit des bornes uniformes en $N$ pour la mesure empirique mollifiée dans les espaces $L^p$.
   • Lemme de Gronwall : Utilisé pour gérer les inégalités intégrales issues de la décomposition de l'erreur.
   • Hypothèses de régularité : Le noyau $K$ est supposé satisfaire des conditions d'intégrabilité spécifiques (Hypothèse AK) et une condition de régularité de type Hölder sur l'opérateur de convolution $K * f$ en termes de la norme $L^1 \cap L^r$.

Contributions clés et résultats
L'article dérive des taux de convergence explicites pour deux objets dans la limite de grand nombre de particules ($N \to \infty$) et la limite de petit pas de temps ($h \to 0$) :

1. Convergence de la mesure empirique mollifiée (Théorème 2.5) :
   L'article établit que la mesure empirique mollifiée $\tilde{\mu}^{N,h}_t$ converge vers la solution $u_t$ de l'équation de Fokker-Planck dans la norme $L^1 \cap L^r$. Le taux de convergence est donné par :
   $$\mathbb{E}\left[ \sup_{s \in [0,t]} \| \tilde{\mu}^{N,h}_s - u_s \|^m_{L^1 \cap L^r} \right]^{1/m} \leq C \left( \mathbb{E}[\| \tilde{\mu}^{N,h}_0 - u_0 \|^m_{L^1 \cap L^r}]^{1/m} + N^{-\rho + \epsilon} + N^{v_2} h^{v_3} \right)$$
   où les exposants dépendent des propriétés du noyau et du paramètre de mollification $\alpha$.

2. Convergence de la densité de particule (Théorème 2.6) :
   L'article fournit un taux de convergence ponctuelle de la densité $u^{i,N,h}_t(x)$ de la $i$-ème particule vers la densité de McKean-Vlasov $u_t(x)$. La borne inclut un terme de poids gaussien :
   $$|u^{i,N,h}_t(x) - u_t(x)| \leq C (\dots) \left( \int_{\mathbb{R}^d} g_{c/p}(t, x-z) u_0(dz) \right)^{1/p}$$
   Ce résultat quantifie l'erreur faible pour le schéma EM, même pour des fonctions de test pouvant être des masses de Dirac.

3. Taux explicites pour des noyaux spécifiques (Section 4) :
   L'auteur applique ses résultats généraux à des cas spécifiques :

   • Noyaux Lipschitz bornés : Récupérant des taux d'ordre $O(N^{-\frac{1}{d+2}} + h^{1/2})$.
   • Potentiels de Riesz et noyaux de Coulomb : Pour $d \geq 2$, le taux est $O(N^{-\frac{1}{2(d+1)}} + N^{\frac{d}{2(d+1)}} h^{1/2})$.
   • Noyaux de type Krylov-Röckner : Des taux similaires sont dérivés pour les noyaux satisfaisant $K \in L^p$ avec $p > d$.

Signification et affirmations
L'auteur affirme que, à sa connaissance, ce sont les premiers résultats fournissant des taux de convergence quantitatifs pour le schéma d'Euler-Maruyama appliqué à des systèmes de particules en interaction avec des noyaux d'interaction singuliers.

• Gestion des singularités : Contrairement aux travaux précédents qui nécessitent souvent la continuité Lipschitz de la dérive dans la métrique de Wasserstein (ce qui échoue pour les noyaux singuliers), ce travail utilise une structure de régularité spécifique de l'opérateur de convolution $K * f$ dans les espaces $L^1 \cap L^r$.
• Analyse de la discrétisation : Bien que la convergence quantitative pour le système de particules non discrétisé (temps continu) ait été récemment établie par Olivera et al. [35], cet article étend l'analyse au schéma EM discrétisé, quantifiant explicitement le compromis entre le nombre de particules $N$ et le pas de temps $h$.
• Robustesse : La méthode se montre robuste pour les noyaux qui ne satisfont pas la condition standard de Krylov-Röckner (par exemple, les noyaux de Keller-Segel hautement singuliers où $\alpha \geq 2$), à condition qu'ils satisfont les hypothèses structurelles spécifiques (AK) introduites dans l'article.

L'auteur souligne que les taux dérivés dépendent polynomialement de $N$ et $h$, et que l'analyse suit attentivement la dépendance vis-à-vis de la dimension $d$ et de la force de la singularité du noyau. Ce travail comble le fossé entre les résultats théoriques de propagation du chaos et la mise en œuvre numérique pratique pour les modèles d'interaction singuliers.