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: Wave-induced shallow-water monopolar vortex: large-scale experiments
Numerous field observations of tsunami-induced eddies in ports and harbours have been reported for recent tsunami events. We examine the evolution of a turbulent shallow-water monopolar vortex generated by a long wave through a series of large-scale experiments in a rectangular wave basin. A leading-elevation asymmetric wave is guided through a narrow channel to form a flow separation region on the lee side of a straight vertical breakwater, which coupled with the transient flow leads to the formation of a monopolar turbulent coherent structure (TCS). The vortex flow after detachment from the trailing jet is fully turbulent ( $$Re_h \sim O(10^{4}\text {--}10^{5}$$ )) for the remainder of the experimental duration. The free surface velocity field was extracted through particle tracking velocimetry over several experimental trials. The first-order model proposed by Seol & Jirka ( J. Fluid Mech. , vol. 665, 2010, pp. 274–299) to predict the decay and spatial growth of shallow-water vortices fits the experimental data well. Bottom friction is predicted to induce a $$t^{-1}$$ azimuthal velocity decay and turbulent viscous diffusion results in a $$\sqrt {t}$$ bulk vortex radial growth, where $$t$$ represents time. The azimuthal velocity, vorticity and free surface elevation profiles are well described through an idealised geophysical vortex. Kinematic free surface boundary conditions predict weak upwelling in the TCS-centre, followed by a zone of downwelling in a recirculation pattern along the water column. The vertical confinement of the flow is quantified through the ratio of kinetic energy contained in the secondary and primary surface velocity fields and a transition point towards a quasi-two-dimensional flow is identified.  more » « less
Award ID(s):
1661052
PAR ID:
10221501
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
910
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Generating mechanisms and parameterizations for enhanced turbulence in the wake of a seamount in the path of the Kuroshio are investigated. Full-depth profiles of finescale temperature, salinity, horizontal velocity, and microscale thermal-variance dissipation rate up- and downstream of the ∼10-km-wide seamount were measured with EM-APEX profiling floats and ADCP moorings. Energetic turbulent kinetic energy dissipation ratesand diapycnal diffusivitiesabove the seamount flanks extend at least 20 km downstream. This extended turbulent wake length is inconsistent with isotropic turbulence, which is expected to decay in less than 100 m based on turbulence decay time ofN−1∼ 100 s and the 0.5 m s−1Kuroshio flow speed. Thus, the turbulent wake must be maintained by continuous replenishment which might arise from (i) nonlinear instability of a marginally unstable vortex wake, (ii) anisotropic stratified turbulence with expected downstream decay scales of 10–100 km, and/or (iii) lee-wave critical-layer trapping at the base of the Kuroshio. Three turbulence parameterizations operating on different scales, (i) finescale, (ii) large-eddy, and (iii) reduced-shear, are tested. Averageεvertical profiles are well reproduced by all three parameterizations. Vertical wavenumber spectra for shear and strain are saturated over 10–100 m vertical wavelengths comparable to water depth with spectral levels independent ofεand spectral slopes of −1, indicating that the wake flows are strongly nonlinear. In contrast, vertical divergence spectral levels increase withε. 
    more » « less
  2. Accurate characterization of the response of coastal structures when subjected to tsunami-like waves is important for structural engineering assessment and design. The weakly-compressive Smoothed Particle Hydrodynamic (SPH) model can theoretically investigate such phenomena in both horizontal and vertical directions. Yet, the convergence of the solutions is sensitive to physical and numerical parameters used in the modeling. In this paper, multiple three- and two-dimensional SPH models are used to study the numerical convergence of free-surface elevation solutions for various initial inter-particle distances, domain locations along the flume and vicinity of the structure, and unbroken and broken wave flow conditions. The results are used to infer on the trade-offs between the accuracy of the SPH solutions and computational costs of the simulations, including computing time and data storage requirements. Two-dimensional models and an approximate ratio of ten particles per wave height can reasonably predict the nonturbulent unbroken wave case. The broken wave case requires three-dimensional models and four times the ratio of particles per wave height. A correlation between experimental and numerical results is then performed, showing adequacy of the free-surface elevation converged SPH models to capture global force responses. The distribution of horizontal and vertical pressures exerted on the elevated structure are characterized and compared with an analytical equation derived from the experimental dataset, highlighting the symbiotic relationship between experimental data, for calibration of the models, and numerical insights, for physical setup design. For example, additional instruments should be placed at strategic locations in future experimental programs to further validate numerical local responses, such as pressures near the edges and corners of structures. Such insights are important to support future work and development of updated US and European guidelines for the design of overland built infrastructures. 
    more » « less
  3. A Computational Fluid Dynamic (CFD) model study of wave and structure interactions on an elevated residential building under various air gap and surge/wave conditions was performed using the olaFlow, an open-source program using the OpenFOAM (Open-source Fields Operation And Manipulation) platform. The numerical model results, including free surface elevation, wave velocity, and vertical pressures on the underside of the elevated structure, showed a good agreement with the measured time-series data from the 1:6 scale hydraulic experiment (Duncan et al., 2021). The numerical simulations were used to extend the physical model tests by computing the vertical distribution of the pressure and resulting wave-induced horizontal forces/pressures, which were not measured in the physical model studies. The simulated results indicate that the pattern of pressure distributions at the frontal face of the elevated structure was controlled by water depth and wavebreaking types (nonbreaking, breaking, and broken waves). The wave induced-vertical force on the elevated structure strongly depends on wave height and the air gap, which is a net elevation from the still water level to the bottom of the structure, but the horizontal force shows complicated patterns due to the varied surge levels (flow depth), wave heights and air gaps. The new dimensionless parameter, α′/h, comprised of the air gap, incident wave height, and flow depth, is introduced and utilized to predict the horizontal forces on the elevated structure. 
    more » « less
  4. Measurements collected by a REMUS 600 AUV off the coast of southern California demonstrate large-scale coherent wave-driven vortices, consistent with Langmuir turbulence (LT), played a dominant role in structuring turbulent dissipation within the oceanic surface boundary layer. During a 10-hour period with sustained wind speeds of 10 m/s, Langmuir circulations were limited to the upper third of the surface mixed layer by persistent stratification within the water column. The ensemble-averaged circulation, calculated using conditional averaging of AD2CP velocity profiles using elevated backscattering intensity associated with subsurface bubble clouds, indicates that LT vortex pairs were characterized by an energetic downwelling zone flanked by broader, weaker upwelling regions with vertical velocity magnitudes similar to previous numerical studies of LT. Horizontally-distributed microstructure estimates of turbulent kinetic energy dissipation rates were lognormally-distributed near the surface in the wave mixing layer with the majority of values falling between wall layer scaling and wave transport layer scaling. Partitioning dissipation rates between downwelling centers and ambient conditions suggests that LT may play a dominant role in elevating dissipation rates in the OSBL by redistributing wave-breaking turbulence. 
    more » « less
  5. This study utilizes a large-eddy simulation (LES) approach to systematically assess the directional variability of wave-driven Langmuir turbulence (LT) in the ocean surface boundary layer (OSBL) under tropical cyclones (TCs). The Stokes drift vector, which drives LT through the Craik–Leibovich vortex force, is obtained through spectral wave simulations. LT’s direction is identified by horizontally elongated turbulent structures and objectively determined from horizontal autocorrelations of vertical velocities. In spite of a TC’s complex forcing with great wind and wave misalignments, this study finds that LT is approximately aligned with the wind. This is because the Reynolds stress and the depth-averaged Lagrangian shear (Eulerian plus Stokes drift shear) that are key in determining the LT intensity (determined by normalized depth-averaged vertical velocity variances) and direction are also approximately aligned with the wind relatively close to the surface. A scaling analysis of the momentum budget suggests that the Reynolds stress is approximately constant over a near-surface layer with predominant production of turbulent kinetic energy by Stokes drift shear, which is confirmed from the LES results. In this layer, Stokes drift shear, which dominates the Lagrangian shear, is aligned with the wind because of relatively short, wind-driven waves. On the contrary, Stokes drift exhibits considerable amount of misalignments with the wind. This wind–wave misalignment reduces LT intensity, consistent with a simple turbulent kinetic energy model. Our analysis shows that both the Reynolds stress and LT are aligned with the wind for different reasons: the former is dictated by the momentum budget, while the latter is controlled by wind-forced waves. 
    more » « less