skip to main content

Title: Solitary wave fission of a large disturbance in a viscous fluid conduit
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 waves more » of the prediction, and the amplitude distribution shows remarkable agreement. Universal properties of solitary wave fission in other fluid dynamics problems are identified. « less
; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Journal of Fluid Mechanics
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 andmore »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
  2. 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 classesmore »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.« less
  3. 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.

  4. 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 formore »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.« less
  5. We investigate wind wave growth by direct numerical simulations solving for the two-phase Navier–Stokes equations. We consider the ratio of the wave speed $c$ to the wind friction velocity $u_*$ from $c/u_*= 2$ to 8, i.e. in the slow to intermediate wave regime; and initial wave steepness $ak$ from 0.1 to 0.3; the two being varied independently. The turbulent wind and the travelling, nearly monochromatic waves are fully coupled without any subgrid-scale models. The wall friction Reynolds number is 720. The novel fully coupled approach captures the simultaneous evolution of the wave amplitude and shape, together with the underwater boundary layer (drift current), up to wave breaking. The wave energy growth computed from the time-dependent surface elevation is in quantitative agreement with that computed from the surface pressure distribution, which confirms the leading role of the pressure forcing for finite amplitude gravity waves. The phase shift and the amplitude of the principal mode of surface pressure distribution are systematically reported, to provide direct evidence for possible wind wave growth theories. Intermittent and localised airflow separation is observed for steep waves with small wave age, but its effect on setting the phase-averaged pressure distribution is not drastically different from that ofmore »non-separated sheltering. We find that the wave form drag force is not a strong function of wave age but closely related to wave steepness. In addition, the history of wind wave coupling can affect the wave form drag, due to the wave crest shape and other complex coupling effects. The normalised wave growth rate we obtain agrees with previous studies. We make an effort to clarify various commonly adopted underlying assumptions, and to reconcile the scattering of the data between different previous theoretical, numerical and experimental results, as we revisit this longstanding problem with new numerical evidence.« less