Fourier pseudospectral method with Kepler mapping for travelling waves with discontinuous slope: application to corner waves of the Ostrovsky-Hunter equation and equatorial Kelvin waves in the four-mode approximation

Authors:
John P. Boyd
Affiliations:
Department of Atmospheric, Oceanic and Space Science and Laboratory for Scientific Computation, University of Michigan, Ann Arbor, MI
Venue:
Applied Mathematics and Computation
Year:
2006

Citing 3
Cited 1

Solitons from sine waves: analytical and numerical methods for non-integrable solitary and cnoidal waves

Physica D
Weakly nonlinear wavepackets in the Korteweg—de Vries equation: the KdV/NLS connection

Mathematics and Computers in Simulation - IMACS sponsored special issue on nonlinear waves: computation and theory
A Legendre-pseudospectral method for computing travelling waves with corners (slope discontinuities) in one space dimension with application to Whitham's equation family

Journal of Computational Physics

Why Newton's method is hard for travelling waves: Small denominators, KAM theory, Arnold's linear Fourier problem, non-uniqueness, constraints and erratic failure

Mathematics and Computers in Simulation

Quantified Score

Hi-index	0.48

Visualization

Abstract

Many species of travelling waves have a single branch of solutions which ends at finite amplitude with a singular wave whose slope is discontinuous, a so-called "corner wave". Fourier pseudospectral methods converge exponentially fast with N, the truncation of the series, for the smooth waves, but the error falls only as O(1/N) for a wave with a slope discontinuity. We show that the error rate can be accelerated to O(1/N3) by making the change-of-coordinate ("Kepler map") x = z - sin(z). Unfortunately, there is a subtlety: the end-in-a-singular-solution bifurcation is possible only for infinite-dimensional systems. We show that the truncated Fourier pseudospectral approximation has roots on both sides of the corner wave. The bifurcation can be detected by (i) a step-function-like jump in the residual of the differential equation at the corner wave and (ii) by observing the slope of u(x), which becomes discontinuous only at the corner wave itself. These concepts are illustrated for the equatorially trapped ocean Kelvin wave in the so-called "four-latitudinal-mode" approximation. However, corner waves arise in many species of waves and the concepts explained here are applicable to all.