Hybrydowe metody rozwinięć dla ułamkowych nieliniowych układów matematycznych z operatorami pochodnych Erdélyi–Kober w teorii modelowania fal tsunami

Fale z długą pamięcią

Ogromne fale oceaniczne, takie jak tsunami, nie reagują tylko na to, co dzieje się tu i teraz; niosą także pewien rodzaj pamięci tego, co wydarzyło się wcześniej na dnie morskim i wzdłuż ich trasy. Artykuł bada, jak uchwycić tę pamięć w modelach matematycznych fal w płytkiej wodzie. Celem autorów jest stworzenie narzędzi, które potrafią wierniej opisywać ewolucję fal przypominających tsunami niż tradycyjne równania, bez twierdzenia, że bezpośrednio przewidują rzeczywiste katastrofy.

Dlaczego zwykłe równania falowe zawodzą

Standardowe równania falowe traktują wodę jak system reagujący natychmiast i lokalnie: to, co dzieje się w danym miejscu, zależy głównie od warunków właśnie tam i w tej chwili. Rzeczywiste oceany są jednak bardziej subtelne. Osady, złożone linie brzegowe i długotrwałe zaburzenia sprawiają, że bieżący ruch wody jest wpływany przez jej przeszłe zachowanie oraz przez zdarzenia odległe przestrzennie. Matematycy określają te efekty jako interakcje niezasobowe (nielokalne) i pamięć. Aby temu sprostać, badacze korzystają z rachunku „ułamkowego”, który pozwala na pochodne o niecałkowitym rzędzie i naturalnie uwzględnia historię oraz wpływy o długim zasięgu.

Nowy rodzaj rachunku dla morza

Autorzy koncentrują się na szczególnej rodzinie równań związanych z tsunami, znanej jako układ Whitham–Broer–Kaup (WBK), powszechnie używanej jako upraszczający model fal w płytkiej wodzie. Zastępują zwykłą pochodną czasową w tych równaniach pochodną ułamkową Erdélyi–Kober, wyspecjalizowanym operatorem, który potrafi zakodować efekt skalowania i pamięć. Mówiąc prościej, ta modyfikacja pozwala równaniom pamiętać, jak fala ewoluowała w czasie, zamiast reagować wyłącznie na najnowsze „pchnięcie”. W efekcie powstaje „ułamkowy” model WBK, w którym kluczowy parametr kontroluje, jak silnie przeszłe zdarzenia wpływają na kształt fali w chwili obecnej.

Hybrydowe metody ujarzmiania trudnych równań

Te bardziej realistyczne równania są jednocześnie trudniejsze do rozwiązania. Zamiast polegać wyłącznie na brutalnej mocy obliczeniowej, autorzy konstruują dwie metody półanalityczne, które łączą zalety rozwinięć szeregowych i iteracyjnych poprawek. Pierwsza, nazwana metodą nowego rozwinięcia iteracyjnego (ENIM), buduje rozwiązanie krok po kroku, wielokrotnie poprawiając początkowe przybliżenie z wykorzystaniem struktury równań i operatorów ułamkowych. Druga, metoda rozwinięcia homotopijno‑perturbacyjnego (EHPM), stopniowo przekształca prosty problem w pełny układ ułamkowy i śledzi, jak rozwiązanie się deformuje po drodze. W obu przypadkach profil fali i pozioma prędkość wody są reprezentowane jako szczególne szeregi potęgowe o ułamkowych wykładnikach, których współczynniki oblicza się, wykorzystując własności operatorów Erdélyi–Kober.

Co obliczenia ujawniają o falach przypominających tsunami

Aby przetestować swoje metody, autorzy stosują ENIM i EHPM do wzorcowej wersji ułamkowego układu WBK opisującego dynamikę fal tsunami w płytkiej wodzie. Gdy rząd ułamkowy ustawiono tak, by odpowiadał przypadkowi klasycznemu (efektywnie wyłączając pamięć), obie metody odtwarzają znane rozwiązanie dokładne z wysoką precyzją; EHPM konsekwentnie daje mniejsze błędy niż ENIM, często około połowy ich wielkości. Jeśli rząd ułamkowy zmniejszy się poniżej wartości klasycznej, fale modelowe zmieniają charakter: stają się gładsze i bardziej rozproszone, o niższych szczytach i szerszych profilach, co odzwierciedla silniejszą pamięć i efekty rozpraszania. W miarę zwiększania rzędu ułamkowego w kierunku wartości klasycznej, fale stają się ostrzejsze i zachowują się bardziej jak znane, wyraźnie zdefiniowane pulsacje pędzące w stronę brzegu.

Dlaczego to ma znaczenie dla przyszłego modelowania fal

Badanie konkluduje, że te hybrydowe metody rozwinięć są stabilnymi, wydajnymi i dokładnymi narzędziami do obsługi złożonych równań falowych uwzględniających pamięć. Chociaż wyniki nie są przeznaczone do bezpośredniego prognozowania rzeczywistych tsunami, pokazują, jak modele ułamkowe mogą płynnie łączyć zachowania dyfuzyjne, silnie zależne od historii, z ostrzejszymi falami opisywanymi przez równania tradycyjne. To czyni je obiecującymi elementami konstrukcyjnymi dla przyszłych badań matematycznych nad zjawiskami przypominającymi tsunami oraz innymi systemami, w których przeszłość wciąż kształtuje teraźniejszość.

Cytowanie: Damag, F.H., Saif, A., Alshammari, M. et al. Hybrid expansion methods for fractional non-linear mathematical systems with Erdelyi-Kober derivative operators in theory of tsunami wave modeling. Sci Rep 16, 10551 (2026). https://doi.org/10.1038/s41598-026-46268-5

Słowa kluczowe: modelowanie tsunami, rachunek ułamkowy, fale w płytkich wodach, metody półanalityczne, propagacja fal