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: Evolution of interface singularities in shallow water equations with variable bottom topography
Abstract Wave front propagation with nontrivial bottom topography is studied within the formalism of hyperbolic long wave models. Evolution of nonsmooth initial data is examined, and, in particular, the splitting of singular points and their short time behavior is described. In the opposite limit of longer times, the local analysis of wave fronts is used to estimate the gradient catastrophe formation and how this is influenced by the topography. The limiting cases when the free surface intersects the bottom boundary, belonging to the so‐called “physical” and “nonphysical” vacuum classes, are examined. Solutions expressed by power series in the spatial variable lead to a hierarchy of ordinary differential equations for the time‐dependent series coefficients, which are shown to reveal basic differences between the two vacuum cases: for nonphysical vacuums, the equations of the hierarchy are recursive and linear past the first two pairs, whereas for physical vacuums, the hierarchy is nonrecursive, fully coupled, and nonlinear. The former case may admit solutions that are free of singularities for nonzero time intervals, whereas the latter is shown to develop nonstandard velocity shocks instantaneously. Polynomial bottom topographies simplify the hierarchy, as they contribute only a finite number of inhomogeneous forcing terms to the equations in the recursion relations. However, we show that truncation to finite‐dimensional systems and polynomial solutions is in general only possible for the case of a quadratic bottom profile. In this case, the system's evolution can reduce to, and is completely described by, a low‐dimensional dynamical system for the time‐dependent coefficients. This system encapsulates all the nonlinear properties of the solution for general power series initial data, and, in particular, governs the loss of regularity in finite times at the dry point. For the special case of parabolic bottom topographies, an exact, self‐similar solution class is introduced and studied to illustrate via closed‐form expressions the general results.  more » « less
Award ID(s):
1910824
PAR ID:
10446111
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
Wiley-Blackwell
Date Published:
Journal Name:
Studies in Applied Mathematics
Volume:
148
Issue:
4
ISSN:
0022-2526
Page Range / eLocation ID:
p. 1439-1476
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    This paper studies the local unstable manifold attached to an equilibrium solution of a system of delay differential equations (DDEs). Two main results are developed. The first is a general method for computing the formal Taylor series coefficients of a function parameterizing the unstable manifold. We derive linear systems of equations whose solutions are the Taylor coefficients, describe explicit formulas for assembling the linear equations for DDEs with polynomial nonlinearities. We also discuss a scheme for transforming non-polynomial DDEs into polynomial ones by appending auxiliary equations. The second main result is an a-posteriori theorem which—when combined with deliberate control of rounding errors— leads to mathematically rigorous computer assisted convergence results and error bounds for the truncated series. Our approach is based on the parameterization method for invariant manifolds and requires some mild non-resonance conditions between the unstable eigenvalues 
    more » « less
  2. null (Ed.)
    Abstract General rogue waves in (1+1)-dimensional three-wave resonant interaction systems are derived by the bilinear method. These solutions are divided into three families, which correspond to a simple root, two simple roots and a double root of a certain quartic equation arising from the dimension reduction, respectively. It is shown that while the first family of solutions associated with a simple root exists for all signs of the nonlinear coefficients in the three-wave interaction equations, the other two families of solutions associated with two simple roots and a double root can only exist in the so-called soliton-exchange case, where the nonlinear coefficients have certain signs. Many of these rogue wave solutions, such as those associated with two simple roots, the ones generated by a $$2\times 2$$ block determinant in the double-root case, and higher-order solutions associated with a simple root, are new solutions which have not been reported before. Technically, our bilinear derivation of rogue waves for the double-root case is achieved by a generalization to the previous dimension reduction procedure in the bilinear method, and this generalized procedure allows us to treat roots of arbitrary multiplicities. Dynamics of the derived rogue waves is also examined, and new rogue wave patterns are presented. Connection between these bilinear rogue waves and those derived earlier by Darboux transformation is also explained. 
    more » « less
  3. We report new rogue wave patterns in the nonlinear Schrödinger equation. These patterns include heart-shaped structures, fan-shaped sectors, and many others, that are formed by individual Peregrine waves. They appear when multiple internal parameters in the rogue wave solutions get large. Analytically, we show that these new patterns are described asymptotically by root structures of Adler–Moser polynomials through a dilation. Since Adler–Moser polynomials are generalizations of the Yablonskii–Vorob’ev polynomial hierarchy and contain free complex parameters, these new rogue patterns associated with Adler–Moser polynomials are much more diverse than previous rogue patterns associated with the Yablonskii–Vorob’ev polynomial hierarchy. We also compare analytical predictions of these patterns to true solutions and demonstrate good agreement between them. 
    more » « less
  4. Abstract We show that new types of rogue wave patterns exist in integrable systems, and these rogue patterns are described by root structures of Okamoto polynomial hierarchies. These rogue patterns arise when the τ functions of rogue wave solutions are determinants of Schur polynomials with index jumps of three, and an internal free parameter in these rogue waves gets large. We demonstrate these new rogue patterns in the Manakov system and the three‐wave resonant interaction system. For each system, we derive asymptotic predictions of its rogue patterns under a large internal parameter through Okamoto polynomial hierarchies. Unlike the previously reported rogue patterns associated with the Yablonskii–Vorob'ev hierarchy, a new feature in the present rogue patterns is that the mapping from the root structure of Okamoto‐hierarchy polynomials to the shape of the rogue pattern is linear only to the leading order, but becomes nonlinear to the next order. As a consequence, the current rogue patterns are often deformed, sometimes strongly deformed, from Okamoto‐hierarchy root structures, unless the underlying internal parameter is very large. Our analytical predictions of rogue patterns are compared to true solutions, and excellent agreement is observed, even when rogue patterns are strongly deformed from Okamoto‐hierarchy root structures. 
    more » « less
  5. We consider high-order strongly nonlinear long wave models expanded in a single small parameter measuring the ratio of the water depth to the characteristic wavelength. By examining its dispersion relation, the high-order system for the bottom velocity is found stable to all disturbances at any order of approximation. On the other hand, systems for other velocities can be unstable and even ill-posed, as signified by the unbounded maximum growth. Under the steady assumption, a new third-order solitary wave solution of the Euler equations is obtained using the high-order strongly nonlinear system and is expanded in an amplitude parameter, which is different from that used in weakly nonlinear theory. The third-order solution is shown to well describe various physical quantities induced by a finite-amplitude solitary wave, including the wave profile, horizontal velocity profile, particle velocity at the crest and bottom pressure. For numerical computations, the first- and second-order strongly nonlinear systems for the bottom velocity are considered. It is shown that finite difference schemes are unstable due to truncation errors introduced in approximating high-order spatial derivatives and, therefore, a more accurate spatial discretization scheme is necessary. Using a pseudo-spectral method based on finite Fourier series combined with an iterative scheme for the inversion of a non-local operator, the strongly nonlinear systems are solved numerically for the propagation of a single solitary wave and the head-on collision of two counter-propagating solitary waves of finite amplitudes, and the results are compared with previous laboratory measurements. 
    more » « less