skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Stability of smooth solitary waves under intensity-dependent dispersion
Abstract The nonlinear Schrödinger (NLS) equation in one dimension is considered in the presence of an intensity-dependent dispersion term. We study bright solitary waves with smooth profiles that extend from the limit where the dependence of the dispersion coefficient on the wave intensity is negligible to the limit where the solitary wave becomes singular due to vanishing dispersion coefficient. We analyse and numerically explore the stability for such smooth solitary waves, showing with the help of numerical approximations that the family of solitary waves becomes unstable in an intermediate region between the two limits, while being stable in both limits. This bistability, which has also been observed in other NLS equations with generalized nonlinearity, brings about interesting dynamical transitions from one stable branch to another stable branch, which are explored in direct numerical simulations of the NLS equation with the intensity-dependent dispersion term.  more » « less
Award ID(s):
2110030 2204702
PAR ID:
10634582
Author(s) / Creator(s):
; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
IMA Journal of Applied Mathematics
Volume:
89
Issue:
6
ISSN:
0272-4960
Page Range / eLocation ID:
989 to 1005
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We establish dispersive estimates and local decay estimates for the time evolution of non-self-adjoint matrix Schrödinger operators with threshold resonances in one space dimension. In particular, we show that the decay rates in the weighted setting are the same as in the regular case after subtracting a finite rank operator corresponding to the threshold resonances. Such matrix Schrödinger operators naturally arise from linearizing a focusing nonlinear Schrödinger equation around a solitary wave. It is known that the linearized operator for the 1D focusing cubic NLS equation exhibits a threshold resonance. We also include an observation of a favorable structure in the quadratic nonlinearity of the evolution equation for perturbations of solitary waves of the 1D focusing cubic NLS equation. 
    more » « less
  2. Abstract Whitham modulation equations are derived for the nonlinear Schrödinger equation in the plane ((2+1)‐dimensional nonlinear Schrödinger [2d NLS]) with small dispersion. The modulation equations are obtained in terms of both physical and Riemann‐type variables; the latter yields equations of hydrodynamic type. The complete 2d NLS Whitham system consists of six dynamical equations in evolutionary form and two constraints. As an application, we determine the linear stability of one‐dimensional traveling waves. In both the elliptic and hyperbolic cases, the traveling waves are found to be unstable. This result is consistent with previous investigations of stability by other methods and is supported by direct numerical calculations. 
    more » « less
  3. The generalized nonlinear Schr\"odinger equation with full dispersion (FDNLS) is considered in the semiclassical regime. The Whitham modulation equations are obtained for the FDNLS equation with general linear dispersion and a generalized, local nonlinearity. Assuming the existence of a four-parameter family of two-phase solutions, a multiple-scales approach yields a system of four independent, first-order, quasi-linear conservation laws of hydrodynamic type that correspond to the slow evolution of the two wavenumbers, mass, and momentum of modulated periodic traveling waves. The modulation equations are further analyzed in the dispersionless and weakly nonlinear regimes. The ill-posedness of the dispersionless equations corresponds to the classical criterion for modulational instability (MI). For modulations of linear waves, ill-posedness coincides with the generalized MI criterion, recently identified by Amiranashvili and Tobisch [New J. Phys., 21 (2019), 033029]. A new instability index is identified by the transition from real to complex characteristics for the weakly nonlinear modulation equations. This instability is associated with long wavelength modulations of nonlinear two-phase wavetrains and can exist even when the corresponding one-phase wavetrain is stable according to the generalized MI criterion. Another interpretation is that while infinitesimal perturbations of a periodic wave may not grow, small but finite amplitude perturbations may grow, hence this index identifies a nonlinear instability mechanism for one-phase waves. Classifications of instability indices for multiple FDNLS equations with higher-order dispersion, including applications to finite-depth water waves and the discrete NLS equation, are presented and compared with direct numerical simulations. 
    more » « less
  4. The YBJ equation (Young & Ben Jelloul, J. Marine Res. , vol. 55, 1997, pp. 735–766) provides a phase-averaged description of the propagation of near-inertial waves (NIWs) through a geostrophic flow. YBJ is obtained via an asymptotic expansion based on the limit $$\mathit{Bu}\rightarrow 0$$ , where $$\mathit{Bu}$$ is the Burger number of the NIWs. Here we develop an improved version, the YBJ + equation. In common with an earlier improvement proposed by Thomas, Smith & Bühler ( J. Fluid Mech. , vol. 817, 2017, pp. 406–438), YBJ + has a dispersion relation that is second-order accurate in $$\mathit{Bu}$$ . (YBJ is first-order accurate.) Thus both improvements have the same formal justification. But the dispersion relation of YBJ + is a Padé approximant to the exact dispersion relation and with $$\mathit{Bu}$$ of order unity this is significantly more accurate than the power-series approximation of Thomas et al. (2017). Moreover, in the limit of high horizontal wavenumber $$k\rightarrow \infty$$ , the wave frequency of YBJ + asymptotes to twice the inertial frequency $2f$ . This enables solution of YBJ + with explicit time-stepping schemes using a time step determined by stable integration of oscillations with frequency $2f$ . Other phase-averaged equations have dispersion relations with frequency increasing as $$k^{2}$$ (YBJ) or $$k^{4}$$ (Thomas et al. 2017): in these cases stable integration with an explicit scheme becomes impractical with increasing horizontal resolution. The YBJ + equation is tested by comparing its numerical solutions with those of the Boussinesq and YBJ equations. In virtually all cases, YBJ + is more accurate than YBJ. The error, however, does not go rapidly to zero as the Burger number characterizing the initial condition is reduced: advection and refraction by geostrophic eddies reduces in the initial length scale of NIWs so that $$\mathit{Bu}$$ increases with time. This increase, if unchecked, would destroy the approximation. We show, however, that dispersion limits the damage by confining most of the wave energy to low  $$\mathit{Bu}$$ . In other words, advection and refraction by geostrophic flows does not result in a strong transfer of initially near-inertial energy out of the near-inertial frequency band. 
    more » « less
  5. We consider a strongly nonlinear long wave model for large amplitude internal waves in a three-layer flow between two rigid boundaries. The model extends the two-layer Miyata–Choi–Camassa (MCC) model (Miyata, Proceedings of the IUTAM Symposium on Nonlinear Water Waves , eds. H. Horikawa & H. Maruo, 1988, pp. 399–406; Choi & Camassa, J. Fluid Mech. , vol. 396, 1999, pp. 1–36) and is able to describe the propagation of long internal waves of both the first and second baroclinic modes. Solitary-wave solutions of the model are shown to be governed by a Hamiltonian system with two degrees of freedom. Emphasis is given to the solitary waves of the second baroclinic mode (mode 2) and their strongly nonlinear characteristics that fail to be captured by weakly nonlinear models. In certain asymptotic limits relevant to oceanic applications and previous laboratory experiments, it is shown that large amplitude mode-2 waves with single-hump profiles can be described by the solitary-wave solutions of the MCC model, originally developed for mode-1 waves in a two-layer system. In other cases, however, e.g. when the density stratification is weak and the density transition layer is thin, the richness of the dynamical system with two degrees of freedom becomes apparent and new classes of mode-2 solitary-wave solutions of large amplitudes, characterized by multi-humped wave profiles, can be found. In contrast with the classical solitary-wave solutions described by the MCC equation, such multi-humped solutions cannot exist for a continuum set of wave speeds for a given layer configuration. Our analytical predictions based on asymptotic theory are then corroborated by a numerical study of the original Hamiltonian system. 
    more » « less