skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Friday, September 29 until 11:59 PM ET on Saturday, September 30 due to maintenance. We apologize for the inconvenience.

Title: Strongly nonlinear effects on internal solitary waves in three-layer flows
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
Award ID(s):
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The flow of tides over rough bathymetry in the deep ocean generates baroclinic motion including internal waves and bottom‐trapped tides. The stresses generated by this motion feedback on the amplitude and phase of the large‐scale tide. Quantifying the stresses associated with tidal flow over abyssal hills is especially important, as this scale of bathymetry is often unresolved in global baroclinic tide models, and the stresses must therefore be parameterized. Here, we extend the previous theoretical work of the authors to determine the amplitude, phasing, and vertical location of the stresses exerted on a flow driven by a time‐periodic body force when it encounters rough bathymetry. The theory compares favorably with a suite of fully nonlinear numerical simulations. It is shown that all topographic stresses are applied directly above the bathymetry, leading to a two‐layer baroclinic flow, with the near‐bottom spatial‐mean flow (the benthic tide) strongly modified by topographic stresses, and the flow at height unperturbed by the presence of topography. Our results provide a framework to improve baroclinic tide models by (i) providing a simple parameterization for the subinertial stress which is currently not included in any models, (ii) establishing that parameterized stresses should be applied in the diffusive boundary layer directly above the topography, independent of where internal tides may dissipate, and (iii) identifying a minimum resolution of ∼10 km for baroclinic tidal models to adequately capture wave resonance effects that can significantly impact the magnitude of the benthic tide.

    more » « less
  2. Abstract

    This paper considers two-dimensional steady solitary waves with constant vorticity propagating under the influence of gravity over an impermeable flat bed. Unlike in previous works on solitary waves, we allow for both internal stagnation points and overhanging wave profiles. Using analytic global bifurcation theory, we construct continuous curves of large-amplitude solutions. Along these curves, either the wave amplitude approaches the maximum possible value, the dimensionless wave speed becomes unbounded, or a singularity develops in a conformal map describing the fluid domain. This is stronger than what one would expect from a straightforward generalization of existing results for periodic waves. We also show that an arbitrary solitary wave of elevation with constant vorticity must be supercritical. The existence proof relies on a novel reformulation of the problem as an elliptic system for two scalar functions in a fixed domain, one describing the conformal map of the fluid region and the other the flow beneath the wave.

    more » « less
  3. The interaction of localised solitary waves with large-scale, time-varying dispersive mean flows subject to non-convex flux is studied in the framework of the modified Korteweg–de Vries (mKdV) equation, a canonical model for internal gravity wave propagation and potential vorticity fronts in stratified fluids. The effect of large amplitude, dynamically evolving mean flows on the propagation of localised waves – essentially ‘soliton steering’ by the mean flow – is considered. A recent theoretical and experimental study of this new type of dynamic soliton–mean flow interaction for convex flux has revealed two scenarios where the soliton either transmits through the varying mean flow or remains trapped inside it. In this paper, it is demonstrated that the presence of a non-convex cubic hydrodynamic flux introduces significant modifications to the scenarios for transmission and trapping. A reduced set of Whitham modulation equations is used to formulate a general mathematical framework for soliton–mean flow interaction with non-convex flux. Solitary wave trapping is stated in terms of crossing modulation characteristics. Non-convexity and positive dispersion – common for stratified fluids – imply the existence of localised, sharp transition fronts (kinks). Kinks play dual roles as a mean flow and a wave, imparting polarity reversal to solitons and dispersive mean flows, respectively. Numerical simulations of the mKdV equation agree with modulation theory predictions. The mathematical framework developed is general, not restricted to completely integrable equations like mKdV, enabling application beyond the mKdV setting to other fluid dynamic contexts subject to non-convex flux such as strongly nonlinear internal wave propagation that is prevalent in the ocean. 
    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. Abstract

    A Hamiltonian reduction approach is defined, studied, and finally used to derive asymptotic models of internal wave propagation in density stratified fluids in two-dimensional domains. Beginning with the general Hamiltonian formalism of Benjamin (1986J. Fluid Mech.165445–74) for an ideal, stably stratified Euler fluid, the corresponding structure is systematically reduced to the setup of two homogeneous fluids under gravity, separated by an interface and confined between two infinite horizontal plates. A long-wave, small-amplitude asymptotics is then used to obtain a simplified model that encapsulates most of the known properties of the dynamics of such systems, such as bidirectional wave propagation and maximal amplitude travelling waves in the form of fronts. Further reductions, and in particular devising an asymptotic extension of Dirac’s theory of Hamiltonian constraints, lead to the completely integrable evolution equations previously considered in the literature for limiting forms of the dynamics of stratified fluids. To assess the performance of the asymptotic models, special solutions are studied and compared with those of the parent equations

    more » « less