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

Technische Samenvatting: Stochastische Numerieke Approximatie voor Nietlineaire Fokker-Planck-vergelijkingen met Singuliere Kernels

Probleemstelling
Dit artikel behandelt de numerieke approximatie van oplossingen van nietlineaire Fokker-Planck partiële differentiaalvergelijkingen (PDE's) van de vorm:
$$\partial_t u = \Delta u - \nabla \cdot (u K * u)$$
waarbij de interactiekernel $K$ singulier is. Dergelijke vergelijkingen komen voor in diverse fysische en biologische contexten, waaronder het Keller-Segel chemotaxis-model, Coulomb-potentialen in elektrostatica/zwaartekracht, en de Biot-Savart-wet voor 2D Navier-Stokes-vergelijkingen.

De primaire uitdaging ligt in de driftterm $K * u$, die een singuliere kernel bevat. Standaard numerieke benaderingen, zoals het direct toepassen van het Euler-Maruyama (EM) schema op de bijbehorende McKean-Vlasov stochastische differentiële vergelijking (SDE), kampen met twee moeilijkheden:

1. Dichtheidsberekening: Het EM-schema vereist de evaluatie van de dichtheid van de oplossing in elk tijdstap, wat niet direct beschikbaar is.
2. Singulariteit en Welbedaardheid: Naarmate de tijdstap $h \to 0$ gaat, kunnen de singulariteiten van de interactiekernel de welbedaardheid van het resulterende deeltjessysteem in gevaar brengen.

Om deze problemen te overwinnen, overweegt de auteur een systeem van $N$ interagerende deeltjes die de oplossing benaderen via een empirische maat. Vanwege de singuliere drift wordt het systeem echter geregulariseerd met behulp van een mollifier-sequentie $(V^N)$ en een gladde afsnijfunctie $F_A$ om de drift te begrenzen, wat numerieke stabiliteit en theoretische begrensdheid garandeert.

Methodologie
Het artikel analyseert de convergentie van het geregulariseerde, gediscretiseerde deeltjessysteem (Vergelijking 1.4/2.2) naar de oplossing van de nietlineaire Fokker-Planck-vergelijking. De methodologie omvat:

1. Regulering en Discretisatie:

   • De singuliere kernel $K$ wordt geconvolveerd met een mollifier $V^N$ (schalingsparameter $\alpha$).
   • Een afsnijfunctie $F_A$ wordt toegepast op de drift om uniforme begrensdheid te waarborgen, consistent met de $L^\infty$-bound van de ware drift $K * u_t$.
   • Het systeem wordt gediscretiseerd met behulp van het Euler-Maruyama-schema met tijdstap $h$.

2. Convergentie-analyse Strategie:
   De analyse ontleedt de fout tussen de gemollificeerde empirische maat $\tilde{\mu}^{N,h}_t$ (en de dichtheid van een enkel deeltje) en de ware oplossing $u_t$ in verschillende componenten:

   • Initiële Fout: Het verschil tussen de initiële gemollificeerde empirische maat en $u_0$.
   • Propagatie van Chaos Fout: Het verschil tussen het geregulariseerde deeltjessysteem en de McKean-Vlasov SDE (niet-gediscretiseerd).
   • Discretisatiefout: De fout geïntroduceerd door het EM-schema (tijddiscretisatie).
   • Stochastische Convolutie: Een term die voortkomt uit de ruis in het deeltjessysteem.

3. Technische Instrumenten:

   • Semigroep-benadering: De auteur maakt gebruik van de milde formulering van de Fokker-Planck-vergelijking met de warmtekernel $g_2(t, \cdot)$.
   • Itô's Formule: Toegepast om de dichtheden van het deeltjessysteem en de McKean-Vlasov-oplossing rond de warmtekernel uit te breiden.
   • Momenten-bounds: De auteur stelt uniforme bounds vast in $N$ voor de gemollificeerde empirische maat in $L^p$-ruimten.
   • Gronwall's Lemma: Gebruikt om de integraalongelijkheid te behandelen die voortvloeit uit de foutontleding.
   • Regulariteitsveronderstellingen: De kernel $K$ wordt aangenomen te voldoen aan specifieke integrabiliteitsvoorwaarden (Aanname AK) en een Hölder-type regulariteitsvoorwaarde op de convolutieoperator $K * f$ in termen van de $L^1 \cap L^r$ norm.

Belangrijkste Bijdragen en Resultaten
Het artikel leidt expliciete convergentiesnelheden af voor twee objecten in de limiet van een groot aantal deeltjes ($N \to \infty$) en een kleine tijdstap ($h \to 0$):

1. Convergentie van de Gemollificeerde Empirische Maat (Theorema 2.5):
   Het artikel stelt vast dat de gemollificeerde empirische maat $\tilde{\mu}^{N,h}_t$ convergeert naar de oplossing $u_t$ van de Fokker-Planck-vergelijking in de $L^1 \cap L^r$ norm. De convergentiesnelheid wordt gegeven door:
   $$\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)$$
   waarbij de exponenten afhangen van de eigenschappen van de kernel en de mollificatieparameter $\alpha$.

2. Convergentie van de Deeltjesdichtheid (Theorema 2.6):
   Het artikel geeft een puntgewijze convergentiesnelheid voor de dichtheid $u^{i,N,h}_t(x)$ van het $i$-de deeltje naar de McKean-Vlasov-dichtheid $u_t(x)$. De bound bevat een Gaussische wegingsterm:
   $$|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}$$
   Dit resultaat kwantificeert de zwakke fout voor het EM-schema, zelfs voor testfuncties die Dirac-massa's kunnen zijn.

3. Expliciete Snelheden voor Specifieke Kernels (Sectie 4):
   De auteur past de algemene resultaten toe op specifieke gevallen:

   • Begrensde Lipschitz-kernels: Herstel van snelheden van de orde $O(N^{-\frac{1}{d+2}} + h^{1/2})$.
   • Riesz-potentialen en Coulomb-kernels: Voor $d \geq 2$, is de snelheid $O(N^{-\frac{1}{2(d+1)}} + N^{\frac{d}{2(d+1)}} h^{1/2})$.
   • Krylov-Röckner Type Kernels: Vergelijkbare snelheden worden afgeleid voor kernels die voldoen aan $K \in L^p$ met $p > d$.

Betekenis en Claims
Het artikel claimt dat, naar beste weten van de auteur, dit de eerste resultaten zijn die kwantitatieve convergentiesnelheden bieden voor het Euler-Maruyama-schema toegepast op interagerende deeltjelsystemen met singuliere interactiekernels.

• Omgaan met Singulariteiten: In tegenstelling tot eerdere werken die vaak Lipschitz-continuïteit van de drift in de Wasserstein-metriek vereisen (wat faalt voor singuliere kernels), maakt dit werk gebruik van een specifieke regulariteitsstructuur van de convolutieoperator $K * f$ in $L^1 \cap L^r$-ruimten.
• Discretisatie-analyse: Hoewel kwantitatieve convergentie voor het niet-gediscretiseerde (continue tijd) deeltjessysteem recentelijk is vastgesteld door Olivera et al. [35], breidt dit artikel de analyse uit naar het gediscretiseerde EM-schema, waarbij de afruil tussen het aantal deeltjes $N$ en de tijdstap $h$ expliciet wordt gekwantificeerd.
• Robuustheid: De methode wordt getoond robuust te zijn voor kernels die niet voldoen aan de standaard Krylov-Röckner-voorwaarde (bijv. zeer singuliere Keller-Segel kernels waar $\alpha \geq 2$), mits zij voldoen aan de specifieke structurele aannames (AK) die in het artikel worden geïntroduceerd.

De auteur benadrukt dat de afgeleide snelheden polynomiaal afhankelijk zijn van $N$ en $h$, en dat de analyse zorgvuldig de afhankelijkheid van de dimensie $d$ en de singulariteitssterkte van de kernel bijhoudt. Het werk overbrugt de kloof tussen theoretische resultaten over propagatie van chaos en praktische numerieke implementatie voor singuliere interactiemodellen.