An official website of the United States government
A strongly nonlinear long-wave approximation is adopted to obtain a high-order model for large-amplitude long internal waves in a two-layer system by assuming the water depth is much smaller than the typical wavelength. When truncated at the first order, the model can be reduced to the regularized strongly nonlinear model of Choiet al.(J. Fluid Mech., vol. 629, 2009, pp. 73–85), which lessens the Kelvin–Helmholtz instability excited by the tangential velocity jump across the interface in the inviscid Miyata–Choi–Camassa (MCC) equations. Using the second-order model, the next-order correction to the internal solitary wave solution of the MCC equations is found and its validity is examined with the Euler solution in terms of the wave profile, the effective wavelength and the velocity profile. It is shown that the correction greatly improves the comparison with the Euler solution for the whole range of wave amplitudes and no further correction is necessary for practical applications. Based on a local stability analysis, the region of stability for the second-order long-wave model is identified in the physical parameter space so that the efficient numerical scheme developed for the first-order model can be used for the second-order model.
more »
« less
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):
- 1634939
- PAR ID:
- 10128452
- 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
-
-
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
-
Abstract This paper evaluates the experimental generation of internal solitary waves (ISWs) in a miscible two-layer system with a free surface using a jet-array wavemaker (JAW). Unlike traditional gate-release experiments, the JAW system generates ISWs by forcing a prescribed vertical distribution of mass flux. Experiments examine three different layer-depth ratios, with ISW amplitudes up to the maximum allowed by the extended Korteweg-de Vries (eKdV) solution. Phase speeds and wave profiles are captured via planar laser-induced fluorescence and the velocity field is measured synchronously using particle imaging velocimetry. Measured properties are directly compared with the eKdV predictions. As expected, small- and intermediate-amplitude waves match well with the corresponding eKdV solutions, with errors in amplitude and phase speed below 10%. For large waves with amplitudes approaching the maximum allowed by the eKdV solution, the phase speed and the velocity profiles resemble the eKdV solution while the wave profiles are distorted following the trough. This can potentially be attributed to Kelvin-Helmholtz instabilities forming at the pycnocline. Larger errors are generally observed when the local Richardson number at the JAW inlet exceeds the threshold for instability.more » « less
-
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
-
The Degasperis-Procesi (DP) equation is an integrable Camassa-Holm-type model which is an asymptotic approximation for the unidirectional propagation of shallow water waves. This work establishes the orbital stability of localized smooth solitary waves to the DP equation on the real line, extending our previous work on their spectral stability [J. Math. Pures Appl. (9) 142 (2020), pp. 298–314]. The main difficulty stems from the fact that the natural energy space is a subspace of L 3 L^3 , but the translation symmetry for the DP equation gives rise to a conserved quantity equivalent to the L 2 L^2 -norm, resulting in L 3 L^3 higher-order nonlinear terms in the augmented Hamiltonian. But the usual coercivity estimate is in terms of L 2 L^2 norm for DP equation, which cannot be used to control the L 3 L^3 higher order term directly. The remedy is to observe that, given a sufficiently smooth initial condition satisfying some mild constraint, the L ∞ L^\infty orbital norm of the perturbation is bounded above by a function of its L 2 L^2 orbital norm, yielding the higher order control and the orbital stability in the L 2 ∩ L ∞ L^2\cap L^\infty space.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2019.795
