We present a method of detecting bifurcations by locating zeros of a signed version of the smallest singular value of the Jacobian. This enables the use of quadratically convergent root-bracketing techniques or Chebyshev interpolation to locate bifurcation points. Only positive singular values have to be computed, though the method relies on the existence of an analytic or smooth singular value decomposition (SVD). The sign of the determinant of the Jacobian, computed as part of the bidiagonal reduction in the SVD algorithm, eliminates slope discontinuities at the zeros of the smallest singular value. We use the method to search for spatially quasi-periodic traveling water waves that bifurcate from large-amplitude periodic waves. The water wave equations are formulated in a conformal mapping framework to facilitate the computation of the quasi-periodic Dirichlet-Neumann operator. We find examples of pure gravity waves with zero surface tension and overhanging gravity-capillary waves. In both cases, the waves have two spatial quasi-periods whose ratio is irrational. We follow the secondary branches via numerical continuation beyond the realm of linearization about solutions on the primary branch to obtain traveling water waves that extend over the real line with no two crests or troughs of exactly the same shape. The pure gravity wave problem is of relevance to ocean waves, where capillary effects can be neglected. Such waves can only exist through secondary bifurcation as they do not persist to zero amplitude. The gravity-capillary wave problem demonstrates the effectiveness of using the signed smallest singular value as a test function for multi-parameter bifurcation problems. This test function becomes mesh independent once the mesh is fine enough.
more »
« less
Parasitic Capillary Waves on Small-Amplitude Gravity Waves with a Linear Shear Current
This paper describes a numerical investigation of ripples generated on the front face of deep-water gravity waves progressing on a vertically sheared current with the linearly changing horizontal velocity distribution, namely parasitic capillary waves with a linear shear current. A method of fully nonlinear computation using conformal mapping of the flow domain onto the lower half of a complex plane enables us to obtain highly accurate solutions for this phenomenon with the wide range of parameters. Numerical examples demonstrated that, in the presence of a linear shear current, the curvature of surface of underlying gravity waves depends on the shear strength, the wave energy can be transferred from gravity waves to capillary waves and parasitic capillary waves can be generated even if the wave amplitude is very small. In addition, it is shown that an approximate model valid for small-amplitude gravity waves in a linear shear current can reasonably well reproduce the generation of parasitic capillary waves.
more »
« less
- Award ID(s):
- 2108524
- PAR ID:
- 10345106
- Date Published:
- Journal Name:
- Journal of Marine Science and Engineering
- Volume:
- 9
- Issue:
- 11
- ISSN:
- 2077-1312
- Page Range / eLocation ID:
- 1217
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We present a numerical study of spatially quasi-periodic gravity-capillary waves of finite depth in both the initial value problem and travelling wave settings. We adopt a quasi-periodic conformal mapping formulation of the Euler equations, where one-dimensional quasi-periodic functions are represented by periodic functions on a higher-dimensional torus. We compute the time evolution of free surface waves in the presence of a background flow and a quasi-periodic bottom boundary and observe the formation of quasi-periodic patterns on the free surface. Two types of quasi-periodic travelling waves are computed: small-amplitude waves bifurcating from the zero-amplitude solution and larger-amplitude waves bifurcating from finite-amplitude periodic travelling waves. We derive weakly nonlinear approximations of the first type and investigate the associated small-divisor problem. We find that waves of the second type exhibit striking nonlinear behaviour, e.g. the peaks and troughs are shifted non-periodically from the corresponding periodic waves due to the activation of quasi-periodic modes.more » « less
-
We explore basic mechanisms for the instability of finite-amplitude interfacial gravity waves through a two-dimensional linear stability analysis of the periodic and irrotational plane motion of the interface between two unbounded homogeneous fluids of different density in the absence of background currents. The flow domains are conformally mapped into two half-planes, where the time-varying interface is always mapped onto the real axis. This unsteady conformal mapping technique with a suitable representation of the interface reduces the linear stability problem to a generalized eigenvalue problem, and allows us to accurately compute the growth rates of unstable disturbances superimposed on steady waves for a wide range of parameters. Numerical results show that the wave-induced Kelvin–Helmholtz (KH) instability due to the tangential velocity jump across the interface produced by the steady waves is the major instability mechanism. Any disturbances whose dominant wavenumbers are greater than a critical value grow exponentially. This critical wavenumber that depends on the steady wave steepness and the density ratio can be approximated by a local KH stability analysis, where the spatial variation of the wave-induced currents is neglected. It is shown, however, that the growth rates need to be found numerically with care and the successive collisions of eigenvalues result in the wave-induced KH instability. In addition, the present study extends the previous results for the small-wavenumber instability, such as modulational instability, of relatively small-amplitude steady waves to finite-amplitude ones.more » « less
-
Abstract. The design and implementation of boundary conditions for the robust generation and simulation of periodic finite-amplitude internal waves is examined in a quasi two-layer continuous stratification using a spectral-element-method-based incompressible flow solver. The commonly used Eulerian approach develops spurious, and potentially catastrophic small-scale numerical features near the wave-generating boundary in a non-linear stratification when the parameter A/(δc) is sufficiently larger than unity; A and δ are measures of the maximum wave-induced vertical velocity and pycnocline thickness, respectively, and c is the linear wave propagation speed. To this end, an Euler–Lagrange approach is developed and implemented to generate robust high-amplitude periodic deep-water internal waves. Central to this approach is to take into account the wave-induced (isopycnal) displacement of the pycnocline in both the vertical and (effectively) upstream directions. With amplitudes not restricted by the limits of linear theory, the Euler–Lagrange-generated waves maintain their structural integrity as they propagate away from the source. The advantages of the high-accuracy numerical method, whose minimal numerical dissipation cannot damp the above near-source spurious numerical features of the purely Eulerian case, can still be preserved and leveraged further along the wave propagation path through the robust reproduction of the non-linear adjustments of the waveform. The near- and far-source robustness of the optimized Euler–Lagrange approach is demonstrated for finite-amplitude waves in a sharp quasi two-layer continuous stratification representative of seasonally stratified lakes. The findings of this study provide an enabling framework for two-dimensional simulations of internal swash zones driven by well-developed non-linear internal waves and, ultimately, the accompanying turbulence-resolving three-dimensional simulations.more » « less
-
These files contain data supporting all results reported in Lloret et al. "A robust numerical method for the generation and propagation of periodic finite-amplitude internal waves in natural waters using high-accuracy simulations". In Lloret et al. we found: The design and implementation of boundary conditions for the robust generation and simulation of periodic finite-amplitude internal waves is examined in a quasi two-layer continuous stratification using a spectralelement-method-based incompressible flow solver. The commonly used Eulerian approach develops spurious, and potentially catastrophic small-scale numerical features near the wave-generating boundary in a non-linear stratification when the parameter A/(δc) is sufficiently larger than unity; A and δ are measures of the maximum wave-induced vertical velocity and pycnocline thickness, respectively, and c is the linear wave propagation speed. To this end, an Euler–Lagrange approach is developed and implemented to generate robust high-amplitude periodic deep-water internal waves. Central to this approach is to take into account the wave- induced (isopycnal) displacement of the pycnocline in both the vertical and (effectively) upstream directions. With amplitudes not restricted by the limits of linear theory, the Euler–Lagrange-generated waves maintain their structural integrity as they propagate away from the source. The advantages of the high-accuracy numerical method, whose minimal numerical dissipation cannot damp the above near-source spurious numerical features of the purely Eulerian case, can still be preserved and leveraged further along the wave propagation path through the robust reproduction of the non-linear adjustments of the waveform. The near- and far-source robustness of the optimized Euler–Lagrange approach is demonstrated for finite-amplitude waves in a sharp quasi two- layer continuous stratification representative of seasonally stratified lakes. The findings of this study provide an enabling framework for two-dimensional simulations of internal swash zones driven by well-developed non- linear internal waves and, ultimately, the accompanying turbulence-resolving three-dimensional simulations. Please cite as: Lloret, P., Diamessis, P., Stastna, M., & Thomsen, G. N. (2024). Data and scripts from: A robust numerical method for the generation and propagation of periodic finite-amplitude internal waves in natural waters using high-accuracy simulations [Data set]. Cornell University eCommons Repository. https://doi.org/10.7298/5VKW-0303more » « less
An official website of the United States government

