Numerical Study of Nonlinear Dispersive Wave Models with SpecTraVVave

In nonlinear dispersive evolution equations, the competing effects of nonlinearity and dispersion make a number of interesting phenomena possible. In the current work, the focus is on the numerical approximation of traveling-wave solutions of such equations. We describe our efforts to write a dedicated Python code which is able to compute traveling-wave solutions of nonlinear dispersive equations of the general form \begin{equation*} u_t + [f(u)]_{x} + \mathcal{L} u_x = 0, \end{equation*} where $\mathcal{L}$ is a self-adjoint operator, and $f$ is a real-valued function with $f(0) = 0$. The SpectraVVave code uses a continuation method coupled with a spectral projection to compute approximations of steady symmetric solutions of this equation. The code is used in a number of situations to gain an understanding of traveling-wave solutions. The first case is the Whitham equation, where numerical evidence points to the conclusion that the main bifurcation branch features three distinct points of interest, namely a turning point, a point of stability inversion, and a terminal point which corresponds to a cusped wave. The second case is the so-called modified Benjamin-Ono equation where the interaction of two solitary waves is investigated. It is found that is possible for two solitary waves to interact in such a way that the smaller wave is annihilated. The third case concerns the Benjamin equation which features two competing dispersive operators. In this case, it is found that bifurcation curves of periodic traveling-wave solutions may cross and connect high up on the branch in the nonlinear regime

Conservation laws, exact travelling waves and modulation instability for an extended nonlinear Schr\"odinger equation

We study various properties of solutions of an extended nonlinear Schr\"{o}dinger (ENLS) equation, which arises in the context of geometric evolution problems -- including vortex filament dynamics -- and governs propagation of short pulses in optical fibers and nonlinear metamaterials. For the periodic initial-boundary value problem, we derive conservation laws satisfied by local in time, weak $H^2$ (distributional) solutions, and establish global existence of such weak solutions. The derivation is obtained by a regularization scheme under a balance condition on the coefficients of the linear and nonlinear terms -- namely, the Hirota limit of the considered ENLS model. Next, we investigate conditions for the existence of traveling wave solutions, focusing on the case of bright and dark solitons. The balance condition on the coefficients is found to be essential for the existence of exact analytical soliton solutions; furthermore, we obtain conditions which define parameter regimes for the existence of traveling solitons for various linear dispersion strengths. Finally, we study the modulational instability of plane waves of the ENLS equation, and identify important differences between the ENLS case and the corresponding NLS counterpart. The analytical results are corroborated by numerical simulations, which reveal notable differences between the bright and the dark soliton propagation dynamics, and are in excellent agreement with the analytical predictions of the modulation instability analysis.Comment: 27 pages, 5 figures. To be published in Journal of Physics A: Mathematical and Theoretica