Turbulent processes in the ocean surface boundary layer (OSBL) play a key role in weather and climate systems. This study explores a Lagrangian analysis of wave-driven OSBL turbulence, based on a large-eddy simulation (LES) model coupled to a Lagrangian stochastic model (LSM). Langmuir turbulence (LT) is captured by Craik–Leibovich wave forcing that generates LT through the Craik–Leibovich type 2 (CL2) mechanism. Breaking wave (BW) effects are modeled by a surface turbulent kinetic energy flux that is constrained by wind energy input to surface waves. Unresolved LES subgrid-scale (SGS) motions are simulated with the LSM to be energetically consistent with the SGS model of the LES. With LT, Lagrangian autocorrelations of velocities reveal three distinct turbulent time scales: an integral, a dispersive mixing, and a coherent structure time. Coherent structures due to LT result in relatively narrow peaks of Lagrangian frequency velocity spectra. With and without waves, the high-frequency spectral tail is consistent with expectations for the inertial subrange, but BWs substantially increase spectral levels at high frequencies. Consistently, over short times, particle-pair dispersion results agree with the Richardson–Obukhov law, and near-surface dispersion is significantly enhanced because of BWs. Over longer times, our dispersion results are consistent with Taylor dispersion. In this case, turbulent diffusivities are substantially larger with LT in the crosswind direction, but reduced in the along-wind direction because of enhanced turbulent transport by LT that reduces mean Eulerian shear. Our results indicate that the Lagrangian analysis framework is effective and physically intuitive to characterize OSBL turbulence.
We examine how Eulerian statistics of wave breaking and associated turbulence dissipation rates in a field of intermittent events compare with those obtained from sparse Lagrangian sampling by surface following drifters. We use a polydisperse two-fluid model with large-eddy simulation (LES) resolution and volume-of-fluid surface reconstruction (VOF) to simulate the generation and evolution of turbulence and bubbles beneath short-crested wave breaking events in deep water. Bubble contributions to dissipation and momentum transfer between the water and air phases are considered. Eulerian statistics are obtained from the numerical results, which are available on a fixed grid. Next, we sample the LES/VOF model results with a large number of virtual surface-following drifters that are initially distributed in the numerical domain, regularly or irregularly, before each breaking event. Time-averaged Lagrangian statistics are obtained using the time series sampled by the virtual drifters. We show that convergence of statistics occurs for signals that have minimum length of approximately 1000–3000 wave periods with randomly spaced observations in time and space relative to three-dimensional breaking events. We further show important effects of (i) extent of measurements over depth and (ii) obscuration of velocity measurements due to entrained bubbles, which are the two typical challenges in most of the available in situ observations of upper ocean wave breaking turbulence. An empirical correction factor is developed and applied to the previous observations of Thomson et al. Applying the new correction factor to the observations noticeably improves the inferred energy balance of wind input rates and turbulence dissipation rates. Finally, both our simulation results and the corrected observations suggested that the total wave breaking dissipation rates have a nearly linear relation with active whitecap coverage.
more » « less- NSF-PAR ID:
- 10140230
- Publisher / Repository:
- American Meteorological Society
- Date Published:
- Journal Name:
- Journal of Physical Oceanography
- Volume:
- 50
- Issue:
- 4
- ISSN:
- 0022-3670
- Page Range / eLocation ID:
- p. 867-885
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
more » « less
Abstract In shallow coastal oceans, turbulent flows driven by surface winds and waves and constrained by a solid bottom disperse particles. This work examines the mechanisms driving horizontal and vertical dispersion of buoyant and sinking particles for times much greater than turbulent integral time scales. Turbulent fields are modeled using a wind‐stress driven large eddy simulation (LES), incorporating wave‐driven Langmuir turbulence, surface breaking wave turbulent kinetic energy inputs, and a solid bottom boundary. A Lagrangian stochastic model is paired to the LES to incorporate Lagrangian particle tracking. Within a subset of intermediate buoyant rise velocities, particles experience synergistic vertical mixing in which breaking waves (BW) inject particles into Langmuir downwelling velocities sufficient to drive deep mixing. Along‐wind dispersion is controlled by vertical shear in mean along‐wind velocities. Wind and bottom friction‐driven vertical shear enhances dispersion of buoyant and sinking particles, while energetic turbulent mixing, such as from BW, dampens shear dispersion. Strongly rising and sinking particles trapped at the ocean surface and bottom, respectively, experience no vertical shear, resulting in low rates of along‐wind dispersion. Crosswind dispersion is shaped by particle advection in wind‐aligned fields of counter‐rotating Langmuir and Couette roll cells. Langmuir cells enhance crosswind dispersion in neutrally to intermediately buoyant particles through enhanced cell hopping. Surface trapping restricts particles to Langmuir convergence regions, strongly inhibiting crosswind dispersion. In shallow coastal systems, particle dispersion depends heavily on particle buoyancy and wave‐dependent turbulent effects.
-
null (Ed.)The focus of this paper is a numerical simulation study of the flow dynamics in a periodic porous medium to analyse the physics of a symmetry-breaking phenomenon, which causes a deviation in the direction of the macroscale flow from that of the applied pressure gradient. The phenomenon is prominent in the range of porosity from 0.43 to 0.72 for circular solid obstacles. It is the result of the flow instabilities formed when the surface forces on the solid obstacles compete with the inertial force of the fluid flow in the turbulent regime. We report the origin and mechanism of the symmetry-breaking phenomenon in periodic porous media. Large-eddy simulation (LES) is used to simulate turbulent flow in a homogeneous porous medium consisting of a periodic, square lattice arrangement of cylindrical solid obstacles. Direct numerical simulation is used to simulate the transient stages during symmetry breakdown and also to validate the LES method. Quantitative and qualitative observations are made from the following approaches: (1) macroscale momentum budget and (2) two- and three-dimensional flow visualization. The phenomenon draws its roots from the amplification of a flow instability that emerges from the vortex shedding process. The symmetry-breaking phenomenon is a pitchfork bifurcation that can exhibit multiple modes depending on the local vortex shedding process. The phenomenon is observed to be sensitive to the porosity, solid obstacle shape and Reynolds number. It is a source of macroscale turbulence anisotropy in porous media for symmetric solid-obstacle geometries. In the macroscale, the principal axis of the Reynolds stress tensor is not aligned with any of the geometric axes of symmetry, nor with the direction of flow. Thus, symmetry breaking in porous media involves unresolved flow physics that should be taken into consideration while modelling flow inhomogeneity in the macroscale.more » « less
-
null (Ed.)We use an online database of a turbulent channel-flow simulation at $Re_\tau =1000$ (Graham et al. J. Turbul. , vol. 17, issue 2, 2016, pp. 181–215) to determine the origin of vorticity in the near-wall buffer layer. Following an experimental study of Sheng et al. ( J. Fluid Mech. , vol. 633, 2009, pp.17–60), we identify typical ‘ejection’ and ‘sweep’ events in the buffer layer by local minima/maxima of the wall stress. In contrast to their conjecture, however, we find that vortex lifting from the wall is not a discrete event requiring $\sim$ 1 viscous time and $\sim$ 10 wall units, but is instead a distributed process over a space–time region at least $1\sim 2$ orders of magnitude larger in extent. To reach this conclusion, we exploit a rigorous mathematical theory of vorticity dynamics for Navier–Stokes solutions, in terms of stochastic Lagrangian flows and stochastic Cauchy invariants, conserved on average backward in time. This theory yields exact expressions for vorticity inside the flow domain in terms of vorticity at the wall, as transported by viscous diffusion and by nonlinear advection, stretching and rotation. We show that Lagrangian chaos observed in the buffer layer can be reconciled with saturated vorticity magnitude by ‘virtual reconnection’: although the Eulerian vorticity field in the viscous sublayer has a single sign of spanwise component, opposite signs of Lagrangian vorticity evolve by rotation and cancel by viscous destruction. Our analysis reveals many unifying features of classical fluids and quantum superfluids. We argue that ‘bundles’ of quantized vortices in superfluid turbulence will also exhibit stochastic Lagrangian dynamics and satisfy stochastic conservation laws resulting from particle relabelling symmetry.more » « less
-
more » « less
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.
