Approximations of dispersive PDEs in the presence of low-regularity randomness

Technische Samenvatting: Approximaties van Dispersieve PDE's in de Aanwezigheid van Laag-Reguliere Willekeur

Probleemstelling
Dit artikel behandelt de numerieke approximatie van de verwachting van het kwadraat van de modulus van de $k$-de Fouriercoëfficiënt van de oplossing van nietlineaire dispersieve partiële differentiaalvergelijkingen (PDE's) met willekeurige initiële data. Specifiek beschouwen de auteurs vergelijkingen van de vorm:
$$i\partial_t u(t, x) + L(\nabla)u(t, x) = |\nabla|^\alpha p(u(t, x), \bar{u}(t, x))$$
met willekeurige initiële data $u(0, x) = v_\eta(x)$, waarbij de willekeur wordt gemodelleerd door een familie van standaard complexe Gaussische variabelen. De doelgrootheid is $E(|u_k(t, v_\eta)|^2)$, die een centrale rol speelt in de studie van golfturbulentie en de afleiding van de golfkinetische vergelijking (WKE).

Klassieke numerieke benaderingen voor dergelijke problemen discretiseren doorgaans de PDE zelf en berekenen vervolgens de verwachting, of discretiseren de oplossing $u$ en evalueren vervolgens $E(|u|^2)$. Deze methoden vereisen vaak hoge regulariteitsveronderstellingen op de initiële data om optimale convergentiesnelheden te bereiken, aangezien de lokale foutbounds afhangen van hogere-orde afgeleiden van de oplossing. De auteurs merken op dat bestaande low-regularity integratoren (bijv. resonantie-gebaseerde schema's) zijn ontwikkeld voor deterministische PDE's, maar nog niet volledig zijn aangepast aan de specifieke structuur van de verwachting van het tweede moment in de aanwezigheid van willekeur.

Methodologie
De auteurs stellen een nieuwe klasse numerieke schema's voor die de verwachting direct discretiseren, in plaats van de PDE of het oplossingspad. Deze aanpak maakt gebruik van de algebraïsche structuur van het probleem om resonantie-annuleringen uit te buiten en regulariteit te winnen op het discrete niveau.

1. Combinatorisch Kader (Paired Decorated Forests):
   De kerninnovatie is de introductie van een nieuwe combinatorische structuur genaamd paired decorated forests. Terwijl eerder werk over low-regularity integratoren (specifiek [13]) gedecoreerde bomen gebruikte om geïtereerde integralen van de oplossing te representeren, representeert dit werk het tweede moment $E(|u_k|^2)$ als een product van twee dergelijke bomen.

   • De bladeren van deze bomen komen overeen met de frequenties van de willekeurige initiële data.
   • De koppeling van bladeren tussen de twee bomen wordt bepaald door de Wick-stelling toegepast op de Gaussische willekeurige variabelen.
   • Deze structuur, aangeduid als $F = T_1 \cdot T_2$, codeert de specifieke koppelingen van frequenties die optreden bij het berekenen van de verwachting van het product van geïtereerde integralen.

2. Resonantie-gebaseerde Discretisatie:
   Het schema rust op een resonantie-gebaseerde discretisatie van de oscillerende geïtereerde integralen die in de boom-expansie voorkomen.

   • De auteurs definiëren een karakter $\Pi$ die gedecoreerde bossen naar oscillerende integralen mapt.
   • Ze introduceren een benaderingsoperator $\Pi_{n,r}$ (waarbij $n$ de veronderstelde regulariteit van de initiële data is en $r$ de orde van het schema) die een Taylor-expansie uitvoert op de "lagere orde" delen van de resonantie, terwijl de "dominante" delen exact worden geïntegreerd.
   • Cruciaal is dat omdat de berekening een product van twee Fouriercoëfficiënten betreft (bijv. $v_k \bar{v}_k$), de regulariteitsvereisten worden verdeeld. De auteurs argumenteren dat termen zoals $k^2 |v_k|^2$ geïnterpreteerd kunnen worden als convoluties ($\nabla v * \nabla \tilde{v}$), die minder afgeleiden vereisen dan de klassieke fouttermen die $k^2 v_k$ bevatten.

3. Recursieve Foutanalyse:
   Het artikel biedt een rigoureuze recursieve definitie voor de lokale foutterm $L^r_{low}(F, n)$. Deze definitie houdt rekening met de specifieke koppelingsstructuur van de bossen, die verschilt van de enkelvoudige boom-casus in [13] omdat de frequenties op de bladeren van $T_1$ en $T_2$ niet noodzakelijkerwijs disjoint zijn. De foutbound wordt getoond te zijn van de vorm $O(\tau^{r+2} L^r_{low})$, waarbij $L^r_{low}$ slechts lagere-orde afgeleiden bevat vergeleken met klassieke methoden.

Belangrijkste Bijdragen

• Nieuwe Combinatorische Structuur: De introductie van "paired decorated forests" om de Wick-contracties die nodig zijn voor het berekenen van het tweede moment van de oplossingen van dispersieve PDE's met willekeurige initiële data systematisch te behandelen.
• Directe Discretisatie van de Verwachting: Een framework dat de verwachting $E(|u_k|^2)$ direct discretiseert. Dit stelt het schema in staat om de verbeterde resonantiestructuur en de inherente winst in regulariteit in de tweede-momentberekening uit te buiten, die verloren gaan als men eerst de oplossing discretiseert.
• Low-Regularity Schema's: De afleiding van algemene low-regularity numerieke schema's voor het tweede moment die aanzienlijk minder regulariteit vereisen op de initiële data dan klassieke splitting-methoden of exponentiële integratoren.
• Stabilisatietechnieken: De voorstel van filterfuncties (bijv. sinc-gebaseerde filters) om de hogere-orde schema's, die ruimtelijke afgeleiden bevatten, te stabiliseren zonder restrictieve CFL-condities op te leggen.

Resultaten
De auteurs passen hun algemene framework toe op twee fundamentele voorbeelden:

1. Nietlineaire Schrödinger (NLS) Vergelijking:

   • Ze construeren eerste-orde en tweede-orde schema's voor het tweede moment.
   • Het eerste-orde schema vereist slechts één ruimtelijke afgeleide op de initiële data (vergeleken met twee voor klassieke methoden).
   • Het tweede-orde schema is expliciet afgeleid in zowel Fourier- als fysieke ruimte, waarbij termen als $\nabla u * \nabla \tilde{u}$ betrokken zijn.
   • Lokale foutanalyse bevestigt de theoretische winst in regulariteit, waarbij fouten van orde $O(\tau^2 |\nabla| v)$ en $O(\tau^3 |\nabla|^2 v)$ worden getoond respectievelijk.

2. Korteweg–de Vries (KdV) Vergelijking:

   • Vergelijkbare eerste- en tweede-orde schema's zijn afgeleid.
   • De analyse benadrukt hoe de specifieke resonantiestructuur van KdV leidt tot annuleringen die de discretisatie vereenvoudigen.
   • Stabilisatiefilters worden geïntroduceerd om de ruimtelijke afgeleiden in het tweede-orde schema te behandelen.

Betekenis en Claims
Het artikel claimt dat door de focus te verleggen van het discretiseren van de PDE naar het discretiseren van de verwachting, men "zwaar kan profiteren van de optimale resonantiestructuur en de onderliggende winst in regulariteit op het eindige dimensionale (discrete) niveau."

De auteurs benadrukken dat hun aanpak onderscheidt van klassieke methoden omdat deze de PDE zelf niet discretiseert. In plaats daarvan benadert het de verwachting met behulp van de specifieke combinatorische eigenschappen van Gaussische willekeurige variabelen (Wick-stelling). Dit maakt de constructie van schema's mogelijk die robuust zijn voor ruwe initiële data (lage regulariteit), een regime waarin klassieke methoden vaak falen of excessieve regulariteitsveronderstellingen vereisen.

Het werk wordt gepresenteerd als een generalisatie van de resonantie-gebaseerde discretisatie geïntroduceerd in [13] naar de stochastische setting van willekeurige initiële data. De auteurs merken op dat hoewel de afgeleide foutbounds formeel zijn, ze rigoureus gemaakt kunnen worden met behulp van klassieke Sobolev-ruimtes of discrete Strichartz-schattingen, waarbij zij verwijzen naar eerder werk ([55, 2, 39, 50]) voor dergelijke rigoureuze bewijzen. Het artikel claimt niet zelf de afleiding van de golfkinetische vergelijking op te lossen, maar biedt de noodzakelijke numerieke instrumenten om de grootheden ($E(|u_k|^2)$) te benaderen die centraal staan in dergelijke afleidingen.