Solitary waves of moderate amplitude in the SSGGN equations: the extended KdV-Whitham approximation

Dit artikel toont aan dat de extended KdV-Whitham-benadering en de ruimtelijke formulering van de extended KdV-vergelijking geschikte regularisaties zijn voor het modelleren van solitaire golven van gematigde amplitude, waarbij numerieke vergelijkingen met het SSGGN-systeem aantonen dat deze modellen de niet-convexe dispersieproblemen van de standaard extended KdV-vergelijking oplossen.

De kern

Stel je voor dat je naar een meer kijkt en een grote, perfecte golf ziet die over het water glijdt. Deze golf is een "soliton": een eenzame golfdop die zijn vorm behoudt terwijl hij reist, zonder uit elkaar te vallen. Wiskundigen hebben al eeuwenlang geprobeerd om deze golven met formules te beschrijven.

Deze paper is als het ware een reparatiehandleiding voor de wiskundige modellen die we gebruiken om deze golven te voorspellen, vooral wanneer ze niet heel klein en zwak zijn, maar "matig groot" en krachtig.

Hier is de uitleg in simpele taal, met een paar creatieve vergelijkingen:

1. Het oude gereedschap (De KdV-vergelijking)

Vroeger gebruikten wetenschappers een simpele formule, de KdV-vergelijking. Dit is als een fiets: perfect voor een ritje in het park op een zonnige dag (kleine, zwakke golven). Maar als je de fiets wilt gebruiken om een zware berg af te racen (grote, krachtige golven), begint hij te haperen. Hij kan de complexiteit van de echte natuur niet meer volgen.

2. De nieuwe, zwaardere fiets (De eKdV-vergelijking)

Om grotere golven te kunnen beschrijven, hebben de auteurs een geavanceerdere versie bedacht: de eKdV-vergelijking. Dit is als een sportfiets met een turbo. Hij is veel krachtiger en kan steilere hellingen aan.

Maar er zit een addertje onder het gras:
Wanneer je deze "turbo-fiets" gebruikt in een bepaalde richting (de "slow time" richting), begint hij vreemde trillingen te maken die niet in de echte natuur voorkomen.

• De analogie: Stel je voor dat je een auto bestuurt en plotseling een spooktrein voor je verschijnt die niet echt bestaat. In de wiskunde noemen we dit resonante straling. De formule "hallucineert" een rij kleine golfjes voor de grote soliton, terwijl de echte watergolf dat niet doet. De formule is dus te snel en te gevoelig voor de verkeerde frequenties.

3. De oplossing: Twee manieren om het te fixen

De auteurs tonen aan dat je deze "hallucinaties" kunt stoppen op twee manieren:

Optie A: De route veranderen (Slow Space)
In plaats van de tijd als hoofdpersoon te nemen, kijken we naar de ruimte.

• De analogie: Het is alsof je in plaats van te rennen met de auto, de auto op een trein zet. De trein volgt een spoor dat perfect is afgestemd op de echte natuur. Door de berekening op deze manier te doen, verdwijnen de "spooktreinen" (de resonante straling) volledig. De formule werkt nu weer perfect.

Optie B: De motor aanpassen (eKdV-Whitham)
Als je toch de "turbo-fiets" (de tijd-variant) wilt gebruiken, moet je de motor (de dispersie) vervangen door die van de echte natuur.

• De analogie: Stel je voor dat je een robot bouwt die watergolven nabootst. De robot heeft een eigen brein (de eKdV-formule) dat soms fouten maakt. De auteurs zeggen: "Laat de robot zijn eigen brein houden voor het denken (de niet-lineaire delen), maar vervang zijn oren door de oren van een echte vis (de dispersie van het echte water)."
• Ze noemen dit de eKdV-Whitham benadering. De robot luistert nu naar de echte natuur voor hoe golven zich verspreiden, maar denkt nog steeds als een slimme wiskundige. Hierdoor verdwijnen de hallucinaties en blijft de robot trouw aan de realiteit, zelfs bij grote golven.

4. Het testrijden (De resultaten)

De auteurs hebben hun nieuwe modellen getest in een computer-simulatie:

• Bij positieve golven (golven die omhoog komen): De nieuwe "Whitham-motor" werkt fantastisch. Geen spooktreinen, alleen een perfecte soliton.
• Bij negatieve golven (golven die een put vormen, wat bij water soms gebeurt): De oude turbo-fiets (eKdV) faalt volledig en maakt een puinhoop van de golf. De nieuwe "Whitham-motor" daarentegen is zo goed dat hij ononderscheidbaar is van de echte natuur.

5. Hoe weet je welke formule je moet kiezen?

De paper geeft ook een slimme tip voor de toekomst. Voordat je een dure simulatie draait, kun je kijken naar de energie van de golf.

• De analogie: Het is alsof je een weegschaal gebruikt. Als je ziet dat de energie van de golf voornamelijk in de grote golf zit, kun je de simpele "fiets" (eKdV) gebruiken. Maar als je ziet dat er veel energie in de "schuimkopjes" (straling) zit, moet je direct overschakelen op de "trein" of de "Whitham-motor". Je kunt dit voorspellen zonder de hele simulatie te draaien.

Conclusie

Kortom: De auteurs hebben een manier gevonden om de wiskundige modellen voor watergolven te "repareren". Ze hebben laten zien dat je ofwel de berekeningsmethode moet veranderen (de route), ofwel de "oren" van je model moet vervangen door die van de echte natuur (Whitham). Hierdoor kunnen we nu veel nauwkeuriger voorspellen hoe grote, krachtige golven zich gedragen, zonder dat de computer "hallucineert" met valse golven.

Het is een beetje alsof ze de GPS van een schip hebben geüpdatet zodat hij niet meer vastloopt in de mist, maar altijd de echte kustlijn volgt, zelfs in stormachtig weer.

Technische samenvatting

Titel: Solitaire golven van gemiddelde amplitude in de SSGGN-vergelijkingen: de uitgebreide KdV-Whitham benadering

Auteurs: Benjamin Martin, Dmitri Tseluiko, en Karima Khusnutdinova (Loughborough University, VK)

---

1. Probleemstelling

Het artikel onderzoekt de evolutie van langdurige, matig niet-lineaire oppervlaktewatergolven. Hoewel de klassieke Korteweg-de Vries (KdV)-vergelijking goed werkt voor zwak niet-lineaire golven, verliest deze nauwkeurigheid naarmate de amplitude toeneemt (gemiddelde niet-lineariteit). De uitgebreide KdV-vergelijking (eKdV) is ontwikkeld om deze regio te dekken door hogere-orde termen in de asymptotische expansie mee te nemen.

Een kritiek punt in de bestaande literatuur is dat de eKdV-vergelijking, wanneer deze wordt geformuleerd in trage tijd ($T = \epsilon t$), een artefact vertoont dat niet voorkomt in het onderliggende fysische systeem: resonante straling. Dit fenomeen ontstaat doordat de lineaire dispersierelatie van de eKdV-vergelijking in trage tijd niet-convex is, wat resonantie toelaat tussen golven met grote en kleine golfgetallen. Dit leidt tot een onfysische modulaire golftraine voor de soliton, wat de nauwkeurigheid van het model beperkt voor golven met gemiddelde amplitude.

Het onderliggende "parent"-systeem, de 1D Serre–Su–Gardner–Green–Naghdi (SSGGN) vergelijkingen, is niet-integrabel maar biedt een nauwkeurige asymptotische benadering van de Euler-vergelijkingen zonder de aanname van kleine amplitude, en vertoont geen dergelijke resonante straling.

2. Methodologie

De auteurs hanteren een combinatie van analytische afleidingen, asymptotische analyse en numerieke simulaties:

• Afleiding van modellen: De eKdV-vergelijking wordt afgeleid uit de 1D SSGGN-vergelijkingen (en de Boussinesq-Peregrine vergelijkingen) via asymptotische meervoudige schaal-expansies. Dit gebeurt in twee variabelenstelsels:
  1. Trage tijd ($T, \xi$): Waar $\xi = x - c_0 t$.
  2. Trage ruimte ($X, \xi$): Waar $X = \epsilon x$.
• Analyse van dispersierelaties: De lineaire dispersierelaties van de KdV, eKdV (in beide formuleringen) en de SSGGN worden vergeleken om de convexiteit en het optreden van resonantie te analyseren.
• Regularisatie (eKdVW): Om het probleem van resonante straling in de trage tijd-formulering op te lossen, introduceren de auteurs de uitgebreide KdV-Whitham (eKdVW) benadering. Hierbij worden de niet-lineaire termen van de eKdV-vergelijking behouden, maar wordt de lineaire dispersierelatie vervangen door die van het onderliggende SSGGN-systeem. Dit resulteert in een integro-differentiaalvergelijking.
• Asymptotische oplossingen: Er worden soliton-oplossingen voor de eKdV-vergelijking afgeleid met behulp van de Kodama-Fokas-Liu near-identity transformatie (NIT) die de eKdV-vergelijking mapt naar de KdV-vergelijking. Dit levert een oplossing op die nauwkeurig is tot $O(\epsilon^2)$. Dit wordt vergeleken met oplossingen afgeleid van de verbeterde Gardner-vergelijking.
• Numerieke simulaties: Directe numerieke oplossingen van de SSGGN-vergelijkingen worden vergeleken met de gereduceerde modellen (KdV, eKdV, Gardner, eKdVW). De simulaties gebruiken een pseudospectrale methode voor ruimtelijke afgeleiden en een 4e-orde Runge-Kutta methode voor tijdsintegratie.
• Voorspelling van modelkeuze: Er wordt een methode ontwikkeld om a priori te bepalen welk model het beste werkt, gebaseerd op de Inverse Scattering Transform (IST) van de KdV-vergelijking en de behoudswetten (massa, impuls, energie).

3. Belangrijkste Bijdragen en Resultaten

A. Het probleem met resonante straling

De studie bevestigt dat de eKdV-vergelijking in trage tijd ($T$) een niet-convexe dispersierelatie heeft, wat leidt tot resonante straling voor de soliton. Dit fenomeen is afwezig in de SSGGN-vergelijkingen en in de eKdV-vergelijking in trage ruimte ($X$), waar de dispersierelatie convex blijft.

B. De eKdVW-benadering als oplossing

De introductie van de eKdVW-vergelijking (waarbij de dispersierelatie van SSGGN wordt gebruikt) regulariseert het eKdV-model effectief:

• Voor positieve amplitude solitons elimineert de eKdVW de resonante straling volledig en behoudt hij de nauwkeurigheid van de eKdV voor de soliton zelf, met slechts een minimale fout in fase en amplitude.
• Voor negatieve amplitude golven (waar geen solitons bestaan en de golf uit elkaar valt in een dispersieve trein) presteert de eKdVW aanzienlijk beter dan de standaard eKdV. De standaard eKdV faalt hier volledig door sterke resonantie, terwijl de eKdVW bijna identiek is aan de parent SSGGN-simulaties.

C. Vergelijking van soliton-oplossingen

De auteurs vergelijken twee asymptotische benaderingen voor de eKdV-soliton:

1. Een oplossing afgeleid via NIT naar de KdV-vergelijking.
2. Een oplossing afgeleid via NIT naar de verbeterde Gardner-vergelijking.
   Resultaat: De oplossing gebaseerd op de KdV-NIT levert de beste nauwkeurigheid voor oppervlaktewatergolven, zelfs bij matige amplitudes. De Gardner-benadering wordt beperkt door complexwordende amplitudes bij hogere waarden van $\epsilon$.

D. Voorspelling van het juiste model

Een belangrijke bijdrage is de ontwikkeling van een criterium om te voorspellen welk model (eKdV of eKdVW) het beste is voor een specifieke initiële voorwaarde:

• Door de IST toe te passen op de initiële KdV-data, kan men het aantal solitons en de hoeveelheid stralingsenergie (massa, impuls, energie) bepalen.
• Als de initiële voorwaarde evolueert naar puur solitons (weinig straling), presteert de eKdV het beste.
• Als er een significant deel van de energie in de straling zit (dispersieve golven), presteert de eKdVW (of KdVW) het beste omdat deze de dispersie correcter modelleert.

E. Numerieke Validatie

Numerieke tests tonen aan dat:

• De eKdVW-vergelijking computatieel veel goedkoper is dan de volledige SSGGN-vergelijkingen, maar een vergelijkbare nauwkeurigheid biedt voor gematigde amplitudes.
• De trage ruimte-formulering van eKdV ook een geldige regularisatie is, maar de eKdVW in trage tijd biedt meer flexibiliteit voor het gebruik van dezelfde initiële voorwaarden als het parent-model.

4. Significance en Conclusie

Dit onderzoek biedt een robuuste methode om de beperkingen van de eKdV-vergelijking voor oppervlaktewatergolven van gemiddelde amplitude te overwinnen. De belangrijkste inzichten zijn:

1. Regularisatie is noodzakelijk: De standaard eKdV in trage tijd is ongeschikt voor nauwkeurige modellering van dispersieve golven door kunstmatige resonantie.
2. eKdVW als superieur alternatief: De uitgebreide KdV-Whitham benadering combineert de voordelen van de eKdV (niet-lineariteit) met de nauwkeurige dispersie van het parent-systeem, waardoor het een universeel en nauwkeurig model wordt voor zowel solitons als dispersieve golven.
3. A priori selectie: De auteurs bieden een praktische tool (gebaseerd op IST en behoudswetten) om onderzoekers te helpen het meest geschikte gereduceerde model te kiezen voordat ze dure simulaties uitvoeren.

De studie bevestigt dat het vervangen van de lineaire dispersierelatie door die van het parent-systeem een krachtige techniek is om de geldigheidsgebieden van asymptotische golvenmodellen te verruimen zonder de complexiteit van de volledige niet-lineaire vergelijkingen te hoeven oplossen.