We study numerically the small dispersion limit for the Korteweg-de Vries
(KdV) equation $u_t+6uu_x+\epsilon^{2}u_{xxx}=0$ for $\epsilon\ll1$ and give a
quantitative comparison of the numerical solution with various asymptotic
formulae for small $\epsilon$ in the whole $(x,t)$-plane. The matching of the
asymptotic solutions is studied numerically.
We numerically investigate transverse stability and instability of so-called
cnoidal waves, i.e., periodic traveling wave solutions of the Korteweg-de Vries
equation, under the time-evolution of the Kadomtsev-Petviashvili equation. In
particular, we find that in KP-I small amplitude cnoidal waves are stable (at
least for spatially localized perturbations) and only become unstable above a
certain threshold. In contrast to that, KP-II is found to be stable for all
amplitudes, or, equivalently, wave speeds.
Purely dispersive partial differential equations as the Korteweg-de Vries
equation, the nonlinear Schr\"odinger equation and higher dimensional
generalizations thereof can have solutions which develop a zone of rapid
modulated oscillations in the region where the corresponding dispersionless
equations have shocks or blow-up. To numerically study such phenomena, fourth
order time-stepping in combination with spectral methods is beneficial to
resolve the steep gradients in the oscillatory region.
An efficient algorithm for computing the branching structure of a compact
Riemann surface defined via an algebraic curve is presented. Generators of the
fundamental group of the base of the ramified covering punctured at the
discriminant points of the curve are constructed via a minimal spanning tree of
the discriminant points. This leads to paths of minimal length between the
points, which is important for a later stage where these paths are used as
integration contours to compute periods of the surface.
Physically meaningful periodic solutions to certain integrable partial
differential equations are given in terms of multi-dimensional theta functions
associated to real Riemann surfaces. Typical analytical problems in the
numerical evaluation of these solutions are studied. In the case of
hyperelliptic surfaces efficient algorithms exist even for almost degenerate
surfaces. This allows the numerical study of solitonic limits.
The small dispersion limit of solutions to the Camassa-Holm (CH) equation is
characterized by the appearance of a zone of rapid modulated oscillations. An
asymptotic description of these oscillations is given, for short times, by the
one-phase solution to the CH equation, where the branch points of the
corresponding elliptic curve depend on the physical coordinates via the Whitham
equations. We present a conjecture for the phase of the asymptotic solution. A
numerical study of this limit for smooth hump-like initial data provides strong
evidence for the validity of this conjecture.
The small dispersion limit of solutions to the Camassa-Holm (CH) equation is
characterized by the appearance of a zone of rapid modulated oscillations. An
asymptotic description of these oscillations is given, for short times, by the
one-phase solution to the CH equation, where the branch points of the
corresponding elliptic curve depend on the physical coordinates via the Whitham
equations. We present a conjecture for the phase of the asymptotic solution. A
numerical study of this limit for smooth hump-like initial data provides strong
evidence for the validity of this conjecture.