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: 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 of the prediction, and the amplitude distribution shows remarkable agreement. Universal properties of solitary wave fission in other fluid dynamics problems are identified.  more » « less
Award ID(s):
1816934
PAR ID:
10177608
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
883
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Evolution of a solitary wave travelling along a submerged sill is studied. The disturbance from the sill creates a phase lag along the wave crest between the ambient water depth and the shallower depth over the sill. This phase lag causes wave diffraction between the different parts of the wave, which induces radiating waves off the edge of the sill. The radiating waves act as an outlet for wave energy, resulting in significant and continual amplitude reduction of the solitary wave. Findings from laboratory experiments are confirmed numerically by simulating a much longer propagation distance with different sill breadths. When the sill breadth is narrow, the solitary wave slowly attenuates by wave radiation, maintaining a quasi-steady wave pattern. This is not the case for a broader sill. The resulting phase lag on the sill continually changes the wave pattern and the attenuation rate is substantially greater than the rate for the case of the narrow sill. The significant energy radiation together with the continual change in the wave formation eventually leads to the complete annihilation of the solitary wave in a wave tank. We also report a wave-breaking process along the sill observed in laboratory experiments. This breaking is induced when the wave amplitude on the sill is smaller than the maximum amplitude of a solitary wave in a uniform depth. Also found is the wake-like formation of gravity–capillary waves behind the breaking crest forming on the sill. Other features associated with the breaking are presented. 
    more » « less
  2. 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
  3. Conduits generated by the buoyant dynamics between two miscible Stokes fluids with high viscosity contrast, a type of core–annular flow, exhibit a rich nonlinear wave dynamics. However, little is known about the fundamental wave dispersion properties of the medium. In the present work, a pump is used to inject a time-periodic flow that results in the excitation of propagating small- and large-amplitude periodic travelling waves along the conduit interface. This wavemaker problem is used as a means to measure the linear and nonlinear dispersion relations and corresponding periodic travelling wave profiles. Measurements are favourably compared with predictions from a fully nonlinear, long-wave model (the conduit equation) and the analytically computed linear dispersion relation for two-Stokes flow. A critical frequency is observed, marking the threshold between propagating and non-propagating (spatially decaying) waves. Measurements of wave profiles and the wavenumber–frequency dispersion relation quantitatively agree with wave solutions of the conduit equation. An upshift from the conduit equation's predicted critical frequency is observed and is explained by incorporating a weak recirculating flow into the full two-Stokes flow model. When the boundary condition corresponds to the temporal profile of a nonlinear periodic travelling wave solution of the conduit equation, weakly nonlinear and strongly nonlinear, cnoidal-type waves are observed that quantitatively agree with the conduit nonlinear dispersion relation and wave profiles. This wavemaker problem is an important precursor to the experimental investigation of more general boundary value problems in viscous fluid conduit nonlinear wave dynamics. 
    more » « less
  4. 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
  5. In this paper, nonlinearity associated with intense ultrasound is studied by using the one-dimensional motion of nonlinear shock wave in an ideal fluid. In nonlinear acoustics, the wave speed of different segments of a waveform is different, which causes distortion in the waveform and can result in the formation of a shock (discontinuity). Acoustic pressure of high-intensity waves causes particles in the ideal fluid to vibrate forward and backward, and this disturbance is of relatively large magnitude due to high-intensities, which leads to nonlinearity in the waveform. In this research, this vibration of fluid due to the intense ultrasonic wave is modeled as a fluid pushed by one complete cycle of piston. In a piston cycle, as it moves forward, it causes fluid particles to compress, which may lead to the formation of a shock (discontinuity). Then as the piston retracts, a forward-moving rarefaction, a smooth fan zone of continuously changing pressure, density, and velocity is generated. When the piston stops at the end of the cycle, another shock is sent forward into the medium. The variation in wave speed over the entire waveform is calculated by solving a Riemann problem. This study examined the interaction of shocks with a rarefaction. The flow field resulting from these interactions shows that the shock waves are attenuated to a Mach wave, and the pressure distribution within the flow field shows the initial wave is dissipated. The developed theory is applied to waves generated by 20 KHz, 500 KHz, and 2 MHz transducers with 50, 150, 500, and 1500 W power levels to explore the effect of frequency and power on the generation and decay of shock waves. This work enhances the understanding of the interactions of high-intensity ultrasonic waves with fluids. 
    more » « less