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 of 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.
more »
« less
Transition to turbulence in wind-drift layers
A light breeze rising over calm water initiates an intricate chain of events that culminates in a centimetres-deep turbulent shear layer capped by gravity–capillary ripples. At first, viscous stress accelerates a laminar wind-drift layer until small surface ripples appear. The surface ripples then catalyse the growth of a second instability in the wind-drift layer, which eventually sharpens into along-wind jets and downwelling plumes, before devolving into three-dimensional turbulence. In this paper, we compare laboratory experiments with simplified, wave-averaged numerical simulations of wind-drift layer evolution beneath monochromatic, constant-amplitude surface ripples seeded with random initial perturbations. Despite their simplicity, our simulations reproduce many aspects of the laboratory-based observations – including the growth, nonlinear development and turbulent breakdown the wave-catalysed instability – generally validating our wave-averaged model. But we also find that the simulated development of the wind-drift layer is disturbingly sensitive to the amplitude of the prescribed surface wave field, such that agreement is achieved through suspiciously careful tuning of the ripple amplitude. As a result of this sensitivity, we conclude that wave-averaged models should really describe the coupled evolution of the surface waves together with the flow beneath to be regarded as truly ‘predictive’.
more »
« less
- PAR ID:
- 10548472
- Publisher / Repository:
- Journal of Fluid Mechanics
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 976
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
null (Ed.)Abstract This study investigates the dynamics of velocity shear and Reynolds stress in the ocean surface boundary layer for idealized misaligned wind and wave fields using a large-eddy simulation (LES) model based on the Craik–Leibovich equations, which captures Langmuir turbulence (LT). To focus on the role of LT, the LES experiments omit the Coriolis force, which obscures a stress–current-relation analysis. Furthermore, a vertically uniform body force is imposed so that the volume-averaged Eulerian flow does not accelerate but is steady. All simulations are first spun-up without wind-wave misalignment to reach a fully developed stationary turbulent state. Then, a crosswind Stokes drift profile is abruptly imposed, which drives crosswind stresses and associated crosswind currents without generating volume-averaged crosswind currents. The flow evolves to a new stationary state, in which the crosswind Reynolds stress vanishes while the crosswind Eulerian shear and Stokes drift shear are still present, yielding a misalignment between Reynolds stress and Lagrangian shear (sum of Eulerian current and Stokes drift). A Reynolds stress budgets analysis reveals a balance between stress production and velocity–pressure gradient terms (VPG) that encloses crosswind Eulerian shear, demonstrating a complex relation between shear and stress. In addition, the misalignment between Reynolds stress and Eulerian shear generates a horizontal turbulent momentum flux (due to correlations of along-wind and crosswind turbulent velocities) that can be important in producing Reynolds stress (due to correlations of horizontal and vertical turbulent velocities). Thus, details of the Reynolds stress production by Eulerian and Stokes drift shear may be critical for driving upper-ocean currents and for accurate turbulence parameterizations in misaligned wind-wave conditions.more » « less
-
There has been an increase in recognition of the important role that the boundary layer turbulent flow structure has on wake recovery and concomitant wind farm efficiency. Most research thus far has focused on onshore wind farms, in which the ground surface is static. With the expected growth of offshore wind farms, there is increased interest in turbulent flow structures above wavy, moving surfaces and their effects on offshore wind farms. In this study, experiments are performed to analyze the turbulent structure above the waves in the wake of a fixed-bottom model wind farm, with special emphasis on the conditional averaged Reynolds stresses, using a quadrant analysis. Phase-averaged profiles show a correlation between the Reynolds shear stresses and the curvature of the waves. Using a quadrant analysis, Reynolds stress dependence on the wave phase is observed in the phase-dependent vertical position of the turbulence events. This trend is primarily seen in quadrants 1 and 3 (correlated outward and inward interactions). Quantification of the correlation between the Reynolds shear stress events and the surface waves provides insight into the turbulent flow mechanisms that influence wake recovery throughout the wake region and should be taken into consideration in wind turbine operation and placement.more » « less
-
Abstract Observations in the Pacific Equatorial Undercurrents (EUC) show that the nighttime deep-cycle turbulence (DCT) in the marginal-instability (MI) layer of the EUC exhibits seasonal variability that can modulate heat transport and sea surface temperature. Large-eddy simulations (LES), spanning a wide range of control parameters, are performed to identify the key processes that influence the turbulent heat flux at multiple time scales ranging from turbulent (minutes to hours) to daily to seasonal. The control parameters include wind stress, convective surface heat flux, shear magnitude, and thickness of the MI layer. In the LES, DCT occurs in discrete bursts during the night, exhibits high temporal variability within a burst, and modulates the mixed layer depth. At the daily time scale, turbulent heat flux generally increases with increasing wind stress, MI-layer shear, or nighttime convection. Convection is found to be important to mixing under weak wind, weak shear conditions. A parameterization for the daily averaged turbulent heat flux is developed from the LES suite to infer the variability of heat flux at the seasonal time scale. The LES-based parameterized heat flux, which takes into account the effects of all control parameters, exhibits a seasonal variability that is similar to the observed heat flux from theχ-pods.more » « less