Numerical study of the Serre-Green-Naghdi equations in 2D

S. Gavrilyuk,C. Klein
2024-03-01
Abstract:A detailed numerical study of solutions to the Serre-Green-Naghdi (SGN) equations in 2D with vanishing curl of the velocity field is presented. The transverse stability of line solitary waves, 1D solitary waves being exact solutions of the 2D equations independent of the second variable, is established numerically. The study of localized initial data as well as crossing 1D solitary waves does not give an indication of existence of stable structures in SGN solutions localized in two spatial dimensions. For the numerical experiments, an approach based on a Fourier spectral method with a Krylov subspace technique is applied.
Analysis of PDEs,Numerical Analysis
What problem does this paper attempt to address?
The problems that this paper attempts to solve mainly focus on two aspects: 1. **Research on the Transverse Stability of Line Solitary Waves**: - The paper studies the transverse stability of line solitary waves in the two - dimensional Serre - Green - Naghdi (SGN) equation by numerical methods. Line solitary waves are solutions of one - dimensional solitary waves in two - dimensional equations and do not depend on the second spatial variable. The author verifies the stability of these line solitary waves under transverse perturbations through numerical experiments. - Specifically, the author considers various perturbations of line solitary waves, including deformed line solitary waves and line solitary waves with Gaussian perturbations. Numerical results show that these perturbations will eventually evolve into unperturbed line solitary waves plus radiation, which provides strong evidence for the transverse stability of line solitary waves. 2. **Research on the Stability of Two - Dimensional Local Structures**: - The paper also explores the evolution of two - dimensional local initial data, especially intersecting line solitary waves and Gaussian - type initial data. The author finds that these initial data do not form stable two - dimensional local structures during the evolution process. - For intersecting line solitary waves, the author finds that these waves will eventually separate and not form a stable two - dimensional structure. For Gaussian - type initial data, numerical results show that these data will evolve into ring - like structures, and a depression will appear in the central area, further supporting the conclusion that the SGN equation has a defocusing effect in the two - dimensional case. ### Main Conjectures Based on the above research, the author proposes two main conjectures: 1. **Main Conjecture I**: Line solitary waves, as solutions of the two - dimensional SGN equation, are stable in the transverse direction. 2. **Main Conjecture II**: The SGN equation has no stable localized solitary waves in the two - dimensional case, and the equation has a defocusing effect. ### Research Methods To conduct these studies, the author uses a numerical method based on the Fourier spectral method and Krylov subspace techniques (GMRES) to solve elliptic equations. This method has high precision and high efficiency when dealing with rapidly decaying or periodic functions. The standard fourth - order explicit Runge - Kutta method is used for time integration. ### Conclusions Through detailed numerical experiments, the author provides strong numerical evidence for the transverse stability of line solitary waves and the defocusing effect of the SGN equation in the two - dimensional case. These results not only deepen the understanding of the SGN equation but also provide an important reference for further research in related fields.