Characterization of singular flows of zeroth-order pseudo-differential operators via elliptic eigenfunctions: a numerical study
The propagation of internal gravity waves in stratified media, such as those found in ocean basins and lakes, leads to the development of geometrical patterns called "attractors". These structures accumulate much of the wave energy and make the fluid flow highly singular. In more analytical terms, the cause of this phenomenon has been attributed to the presence of a continuous spectrum in some nonlocal zeroth-order pseudo-differential operators. In this work, we analyze the generation of these attractors from a numerical analysis perspective. First, we propose a high-order pseudo-spectral method to solve the evolution problem (whose long-term behaviour is known to be not square-integrable). Then, we use similar tools to discretize the corresponding eigenvalue problem. Since the eigenvalues are embedded in a continuous spectrum, we compute them using viscous approximations. Finally, we explore the effect that the embedded eigenmodes have on the long-term evolution of the system.
- Clone the repository:
git clone [email protected]:javieralmonacid/zeroth-order-operators.git
- Open Matlab and add the
tools
folder to the path:addpath tools
- Run the scripts (most of the codes can be run independently from each other). Each script produces a .mat files to avoid doing the main computation every time the script is executed.
Just cite the main article:
J. A. Almonacid and N. Nigam. Characterization of singular flows of zeroth-order pseudo-differential operators via elliptic eigenfunctions: a numerical study. Journal of Computational and Applied Mathematics 438 (2024), 115510.
Available on the JCAM website and arXiv.
- ColorBrewer: Attractive and Distinctive Colormaps (GitHub, Mathworks File Exchange),
sptoeplitz
(Mathworks File Exchange),textprogressbar
(Mathworks File Exchange).