skip to main content

Title: Dynamic soliton–mean flow interaction with non-convex flux
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 more » 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. « less
; ;
Award ID(s):
Publication Date:
Journal Name:
Journal of Fluid Mechanics
Sponsoring Org:
National Science Foundation
More Like this
  1. A new type of wave–mean flow interaction is identified and studied in which a small-amplitude, linear, dispersive modulated wave propagates through an evolving, nonlinear, large-scale fluid state such as an expansion (rarefaction) wave or a dispersive shock wave (undular bore). The Korteweg–de Vries (KdV) equation is considered as a prototypical example of dynamic wavepacket–mean flow interaction. Modulation equations are derived for the coupling between linear wave modulations and a nonlinear mean flow. These equations admit a particular class of solutions that describe the transmission or trapping of a linear wavepacket by an unsteady hydrodynamic state. Two adiabatic invariants of motion are identified that determine the transmission, trapping conditions and show that wavepackets incident upon smooth expansion waves or compressive, rapidly oscillating dispersive shock waves exhibit so-called hydrodynamic reciprocity recently described in Maiden et al.  ( Phys. Rev. Lett. , vol. 120, 2018, 144101) in the context of hydrodynamic soliton tunnelling. The modulation theory results are in excellent agreement with direct numerical simulations of full KdV dynamics. The integrability of the KdV equation is not invoked so these results can be extended to other nonlinear dispersive fluid mechanic models.
  2. Abstract The interaction of an oblique line soliton with a one-dimensional dynamic mean flow is analyzed using the Kadomtsev–Petviashvili II (KPII) equation. Building upon previous studies that examined the transmission or trapping of a soliton by a slowly varying rarefaction or oscillatory dispersive shock wave (DSW) in one space and one time dimension, this paper allows for the incident soliton to approach the changing mean flow at a nonzero oblique angle. By deriving invariant quantities of the soliton–mean flow modulation equations—a system of three (1 + 1)-dimensional quasilinear, hyperbolic equations for the soliton and mean flow parameters—and positing the initial configuration as a Riemann problem in the modulation variables, it is possible to derive quantitative predictions regarding the evolution of the line soliton within the mean flow. It is found that the interaction between an oblique soliton and a changing mean flow leads to several novel features not observed in the (1 + 1)-dimensional reduced problem. Many of these interesting dynamics arise from the unique structure of the modulation equations that are nonstrictly hyperbolic, including a well-defined multivalued solution interpreted as a solution of the (2 + 1)-dimensional soliton–mean modulation equations, in which the soliton interacts with the mean flowmore »and then wraps around to interact with it again. Finally, it is shown that the oblique interactions between solitons and DSW solutions for the mean flow give rise to all three possible types of two-soliton solutions of the KPII equation. The analytical findings are quantitatively supported by direct numerical simulations.« less
  3. This paper presents a theoretical and experimental study of the long-standing fluid mechanics problem involving the temporal resolution of a large localised initial disturbance into a sequence of solitary waves. This problem is of fundamental importance in a range of applications, including tsunami and internal ocean wave modelling. This study is performed in the context of the viscous fluid conduit system – the driven, cylindrical, free interface between two miscible Stokes fluids with high viscosity contrast. Owing to buoyancy-induced nonlinear self-steepening balanced by stress-induced interfacial dispersion, the disturbance evolves into a slowly modulated wavetrain and further into a sequence of solitary waves. An extension of Whitham modulation theory, termed the solitary wave resolution method, is used to resolve the fission of an initial disturbance into solitary waves. The developed theory predicts the relationship between the initial disturbance’s profile, the number of emergent solitary waves and their amplitude distribution, quantifying an extension of the well-known soliton resolution conjecture from integrable systems to non-integrable systems that often provide a more accurate modelling of physical systems. The theoretical predictions for the fluid conduit system are confirmed both numerically and experimentally. The number of observed solitary waves is consistently within one to two wavesmore »of the prediction, and the amplitude distribution shows remarkable agreement. Universal properties of solitary wave fission in other fluid dynamics problems are identified.« less
  4. Abstract

    Large-amplitude internal solitary wave (ISW) shoaling, breaking, and run-up was tracked continuously by a dense and rapidly sampling array spanning depths from 500 m to shore near Dongsha Atoll in the South China Sea. Incident ISW amplitudes ranged between 78 and 146 m with propagation speeds between 1.40 and 2.38 m s−1. The ratio between wave amplitude and a critical amplitudeA0controlled breaking type and was related to wave speedcpand depth. Fissioning ISWs generated larger trailing elevation waves when the thermocline was deep and evolved into onshore propagating bores in depths near 100 m. Collapsing ISWs contained significant mixing and little upslope bore propagation. Bores contained significant onshore near-bottom kinetic and potential energy flux and significant offshore rundown and relaxation phases before and after the bore front passage, respectively. Bores on the shallow forereef drove bottom temperature variation in excess of 10°C and near-bottom cross-shore currents in excess of 0.4 m s−1. Bores decelerated upslope, consistent with upslope two-layer gravity current theory, though run-up extentXrwas offshore of the predicted gravity current location. Background stratification affected the bore run-up, withXrfarther offshore when the Korteweg–de Vries nonlinearity coefficientαwas negative. Fronts associated with the shoaling local internal tide, but equal in magnitude tomore »the soliton-generated bores, were observed onshore of 20-m depth.

    « less
  5. Simulations of strongly stratified turbulence often exhibit coherent large-scale structures called vertically sheared horizontal flows (VSHFs). VSHFs emerge in both two-dimensional (2D) and three-dimensional (3D) stratified turbulence with similar vertical structure. The mechanism responsible for VSHF formation is not fully understood. In this work, the formation and equilibration of VSHFs in a 2D Boussinesq model of stratified turbulence is studied using statistical state dynamics (SSD). In SSD, equations of motion are expressed directly in the statistical variables of the turbulent state. Restriction to 2D turbulence facilitates application of an analytically and computationally attractive implementation of SSD referred to as S3T, in which the SSD is expressed by coupling the equation for the horizontal mean structure with the equation for the ensemble mean perturbation covariance. This second-order SSD produces accurate statistics, through second order, when compared with fully nonlinear simulations. In particular, S3T captures the spontaneous emergence of the VSHF and associated density layers seen in simulations of turbulence maintained by homogeneous large-scale stochastic excitation. An advantage of the S3T system is that the VSHF formation mechanism, which is wave–mean flow interaction between the emergent VSHF and the stochastically excited large-scale gravity waves, is analytically understood in the S3T system. Comparisonmore »with fully nonlinear simulations verifies that S3T solutions accurately predict the scale selection, dependence on stochastic excitation strength, and nonlinear equilibrium structure of the VSHF. These results constitute a theory for VSHF formation applicable to interpreting simulations and observations of geophysical examples of turbulent jets such as the ocean’s equatorial deep jets.« less