Craik Leibovich Research Articles

Scaling and Flow Structure of Langmuir Turbulence in Inertial Frames

Abstract This paper provides a framework that unifies the characteristics of Langmuir turbulence, including the vortex force effect, velocity scalings, vertical flow structure, and crosswind spacing between surface streaks. The widely accepted CL2 mechanism is extended to explain the observed maximum alongwind velocity and downwelling velocity below the surface. Balancing the extended mechanism in the Craik–Leibovich equations, the scalings for the alongwind velocity u, crosswind velocity υ, and vertical velocity w are formulated as Here, Uf is the friction velocity, Us is the Stokes drift on the surface, and La = (Uf/Us)1/2 is the Langmuir number. Simulations using the Stratified Ocean Model with Adaptive Refinement in large-eddy simulation (LES-SOMAR) mode validate the scalings and reveal physical similarity for velocity and crosswind spacing. The horizontally averaged velocity along the wind on the surface grows with time, whereas υ/V and w/W are confined. The root-mean-square (rms) of w peaks at wrms/W ≈ 0.85 at a depth of 1.3Zs, where Zs is the e-folding scale of the Stokes drift. The crosswind spacing Lp grows linearly with time but is finally limited by the depth of the water H, with maximum Lp/H = 3.3. This framework agrees with measurement collected in six different field campaigns.

Langmuir turbulence in suspended kelp farms

This study investigates the influence of suspended kelp farms on ocean mixed layer hydrodynamics in the presence of currents and waves. We use the large eddy simulation method, where the wave effect is incorporated by solving the wave-averaged equations. Distinct Langmuir circulation patterns are generated within various suspended farm configurations, including horizontally uniform kelp blocks and spaced kelp rows. Intensified turbulence arises from the farm-generated Langmuir circulation, as opposed to the standard Langmuir turbulence observed without a farm. The creation of Langmuir circulation within the farm is attributed to two primary factors depending on farm configuration: (i) enhanced vertical shear due to kelp frond area density variability, and (ii) enhanced lateral shear due to canopy discontinuity at lateral edges of spaced rows. Both enhanced vertical and lateral shear of streamwise velocity, representing the lateral and vertical vorticity components, respectively, can be tilted into downstream vorticity to create Langmuir circulation. This vorticity tilting is driven by the Craik–Leibovich vortex force associated with the Stokes drift of surface gravity waves. In addition to the farm-generated Langmuir turbulence, canopy shear layer turbulence is created at the farm bottom edge due to drag discontinuity. The intensity of different types of turbulence depends on both kelp frond area density and the geometric configuration of the farm. The farm-generated turbulence has substantial consequences for nutrient supply and kelp growth. These findings also underscore the significance of the presence of obstacle structures in modifying ocean mixed layer characteristics.

A numerical study on Langmuir circulations and coherent vortical structures beneath surface waves

Langmuir circulations (LCs) arise through the interaction between the Lagrangian drift of the surface waves and the wind-driven shear layer. Quasi-streamwise vortices (QSVs) also form in the turbulent shear layer next to a flat surface. Both vortical structures manifest themselves by inducing wind-aligned streaks on the surface. In this study, numerical simulations of a stress-driven turbulent shear layer bounded by monochromatic surface waves are conducted to reveal the vortical structures of LCs and QSVs, and their interactions. The LC structure is educed from conditional averaging guided by the signatures of predominant streaks obtained from empirical mode decomposition; the width of the averaged LC pair is found to be comparable to the most unstable wavelength of the Craik–Leibovich equation. Coherent vortical structures (CVSs) are identified using a detection criterion based on local analysis of the velocity-gradient tensor and their topological geometry; QSVs accumulated beneath the windward surface are found to dominate the distribution. Employing the variable-interval spatial average to the identified QSVs further reveals that QSVs tend to form in the edge vicinity of the surface streaks induced by the LCs. The transport budgets of streamwise enstrophy are examined to reveal the interaction. It is found that QSVs perturb the streaks resulting in a localized streamwise gradient of the spanwise velocity, that is, vertical vorticity. The vertical shear tilts the vertical vorticity, therefore enhancing streamwise enstrophy production and the formation of QSVs. The results highlight the differences in the CVSs between the Langmuir turbulence and the wall turbulence.

Stokes drift and its discontents.

The Stokes velocity , defined approximately by Stokes (1847, Trans. Camb. Philos. Soc., 8, 441–455.), and exactly via the Generalized Lagrangian Mean, is divergent even in an incompressible fluid. We show that the Stokes velocity can be naturally decomposed into a solenoidal component, , and a remainder that is small for waves with slowly varying amplitudes. We further show that arises as the sole Stokes velocity when the Lagrangian mean flow is suitably redefined to ensure its exact incompressibility. The construction is an application of Soward & Roberts’s glm theory (2010, J. Fluid Mech., 661, 45–72. (doi:10.1017/S0022112010002867)) which we specialize to surface gravity waves and implement effectively using a Lie series expansion. We further show that the corresponding Lagrangian-mean momentum equation is formally identical to the Craik–Leibovich (CL) equation with replacing , and we discuss the form of the Stokes pumping associated with both and .This article is part of the theme issue ‘Mathematical problems in physical fluid dynamics (part 1)’.

Generation of attached Langmuir circulations by a suspended macroalgal farm

Abstract

Wind- and Wave-Driven Reynolds Stress and Velocity Shear in the Upper Ocean for Idealized Misaligned Wind-Wave Conditions

AbstractThis 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.

Numerical study of effect of wave phase on Reynolds stresses and turbulent kinetic energy in Langmuir turbulence

Abstract

Wave-resolving simulations of viscous wave attenuation effects on Langmuir circulation

Interaction processes of attenuating surface waves and underlying currents are explicitly simulated with a nonhydrostatic free surface numerical model. The model, designed with special attention to the numerical conservation properties of energy and momentum, is capable of accurately capturing the effect of the virtual wave stress, a viscous momentum transfer from waves to current. Simple experiments with the present model show that the Langmuir circulations in the wave-resolving simulations tend to be stronger than the ones in the Craik–Leibovich equation with corresponding Stokes drift profiles. The underestimation in the Craik–Leibovich equation can be corrected by addressing momentum transfer from attenuating waves to currents as the virtual wave stress at the surface.

Mutual Interaction between Surface Waves and Langmuir Circulations Observed in Wave-Resolving Numerical Simulations

AbstractWave-resolving simulations of monochromatic surface waves and Langmuir circulations (LCs) under an idealized condition are performed to investigate the dynamics of wave–current mutual interaction. When the Froude number (the ratio of the friction velocity of wind stress imposed at the surface and wave phase speed) is large, waves become refracted by the downwind jet associated with LCs and become amplitude modulated in the crosswind direction. In such cases, the simulations using the Craik–Leibovich (CL) equation with a prescribed horizontally uniform Stokes drift profile are found to underestimate the intensity of LCs. Vorticity budget analysis reveals that horizontal shear of Stokes drift induced by the wave modulation tilts the wind-driven vorticity to the downwind direction, intensifying the LCs that caused the waves to be modulated. Such an effect is not reproduced in the CL equation unless the Stokes drift of the waves modulated by LCs is prescribed. This intensification mechanism is similar to the CL1 mechanism in that the horizontal shear of the Stokes drift plays a key role, but it is more likely to occur because the shear in this interaction is automatically generated by the LCs whereas the shear in the CL1 mechanism is retained only when a particular phase relation between two crossing waves is kept locked for many periods.

Wave–vortex interactions and effective mean forces: three basic problems

Three examples of wave–vortex interaction are studied, in analytically tractable weak refraction regimes with attention to the mean recoil forces, local and remote, that are associated with refractive changes in wave pseudomomentum fluxes. Wave-induced mean forces of this kind can be persistent, with cumulative effects, even in the absence of wave dissipation. In each example, a single wavetrain propagates past a single vortex. In the first two examples, in a two-dimensional, non-rotating acoustic or shallow-water setting, the focus is on whether or not the wavetrain overlaps the vortex core. In the overlapping case, the recoil has a local contribution given by the Craik–Leibovich force on the vortex core, the vector product of Stokes drift and mean vorticity. (For a quantum vortex this contribution is called the Iordanskii force arising from the Aharonov–Bohm effect on a phonon current.) However, in all except one special limiting case there are additional “remote” contributions, mediated by Stokes-drift-induced return flows that can intersect the vortex core well away from locations where the waves are refracted. The third example is a non-overlapping, remote-recoil-only example in a rapidly rotating frame, in which the waves are deep-water gravity waves and the mean flow obeys shallow-water quasigeostrophic dynamics. Contrary to what might at first be thought, the Ursell “anti-Stokes flow” induced by the rotation – an Eulerian-mean flow tending to cancel the Stokes drift – fails to suppress remote recoil. There are nontrivial open questions about extending these results to regimes of stronger refraction, especially regarding the scope of the “pseudomomentum rule” for the wave-induced recoil forces.

Large-eddy simulation of small-scale Langmuir circulation and scalar transport

Large-eddy simulation (LES) of a wind- and wave-forced water column based on the Craik–Leibovich (C–L) vortex force is used to understand the structure of small-scale Langmuir circulation (LC) and associated Langmuir turbulence. The LES also serves to understand the role of the turbulence in determining molecular diffusive scalar flux from a scalar-saturated air side to the water side and the turbulent vertical scalar flux in the water side. Previous laboratory experiments have revealed that small-scale LC beneath an initially quiescent air–water interface appears shortly after the initiation of wind-driven gravity–capillary waves and provides the laminar–turbulent transition in wind speeds between 3 and . The LES reveals Langmuir turbulence characterized by multiple scales ranging from small bursting eddies at the surface that coalesce to give rise to larger (centimetre-scale) LC over time. It is observed that the smaller scales account for the bulk of the near-surface turbulent vertical scalar flux. Although the contribution of the larger (centimetre-scale) LC to the near-surface turbulent flux increases over time as these scales emerge and become more coherent, the contribution of the smaller scales remains dominant. The growing LC scales lead to increased vertical scalar transport at depths below the interface and thus greater scalar transfer efficiency. Simulations were performed with a fixed wind stress corresponding to a wind speed but with different wave parameters (wavelength and amplitude) in the C–L vortex force. It is observed that longer wavelengths lead to more coherent, larger centimetre-scale LC providing greater contribution to the turbulent vertical scalar flux away from the surface. In all cases, the molecular diffusive scalar flux at the water surface relaxes to the same statistically steady value after transition to Langmuir turbulence occurs, despite the different wave parameters in the C–L vortex force across the simulations. This implies that the small-scale turbulence intensity and the molecular diffusive scalar flux at the surface scale with the wind shear and not with the wave parameters. Furthermore, it is seen that the Langmuir (wave) forcing (provided by the C–L vortex force) is necessary to trigger the turbulence that induces elevated molecular diffusive scalar flux at the water surface relative to wind-driven flow without wave forcing.

Deciphering the Role of Small-Scale Inhomogeneity on Geophysical Flow Structuration: A Stochastic Approach

AbstractAn important open question in fluid dynamics concerns the effect of small scales in structuring a fluid flow. In oceanic or atmospheric flows, this is aptly captured in wave–current interactions through the study of the well-known Langmuir secondary circulation. Such wave–current interactions are described by the Craik–Leibovich system, in which the action of a wave-induced velocity, the Stokes drift, produces a so-called “vortex force” that causes streaking in the flow. In this work, we show that these results can be generalized as a generic effect of the spatial inhomogeneity of the statistical properties of the small-scale flow components. As demonstrated, this is well captured through a stochastic representation of the flow.

Wind–Wave Misalignment Effects on Langmuir Turbulence in Tropical Cyclone Conditions

AbstractThis 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.

On the physical mechanisms of the two-way coupling between a surface wave field and a circulation consisting of a roll and streak

The governing equations of a surface wave field and a coexisting roll–streak circulation typical of Langmuir circulations or submesoscale frontal circulations are derived to better describe their two-way interactions. The gradients and vertical velocities of the roll–streak circulation induce wave refraction, amplitude modulation and higher-order waves. These changes then produce wave–wave nonlinear forces and divergence of the wave-induced mass transport, both of which in turn affect the circulation. To accurately represent these processes, both a wave theory and a wave-averaged theory are developed without relying on any extrapolation, any spatiotemporal mapping or an approximation that treats the wave-induced mass divergence as being concentrated at the surface. This wave theory finds seven types of current-induced higher-order wave motions. It also determines the wave dynamics such as the governing equation of the wave action density valid in the presence of the complex circulation. The evolution of the wave action density is clearly affected by the upwelling or downwelling. The new wave-averaged theory presents the governing equations of the wave-averaged circulation which satisfies the wave-averaged mass conservation. This circulation is different from the circulation considered to satisfy the mass conservation in the Craik–Leibovich theory, and the difference becomes critical when the wave field evolves due to refraction. In this case, compared to the Craik–Leibovich theory, long waves are more important and also the rolls are more weakly forced.

Wave–vortex interactions, remote recoil, the Aharonov–Bohm effect and the Craik–Leibovich equation

Three examples of non-dissipative yet cumulative interaction between a single wavetrain and a single vortex are analysed, with a focus on effective recoil forces, local and remote. Local recoil occurs when the wavetrain overlaps the vortex core. All three examples comply with the pseudomomentum rule. The first two examples are two-dimensional and non-rotating (shallow water or gas dynamical). The third is rotating, with deep-water gravity waves inducing an Ursell ‘anti-Stokes flow’. The Froude or Mach number, and the Rossby number in the third example, are assumed small. Remote recoil is all or part of the interaction in all three examples, except in one special limiting case. That case is found only within a severely restricted parameter regime and is the only case in which, exceptionally, the effective recoil force can be regarded as purely local and identifiable with the celebrated Craik–Leibovich vortex force – which corresponds, in the quantum fluids literature, to the Iordanskii force due to a phonon current incident on a vortex. Another peculiarity of that exceptional case is that the only significant wave refraction effect is the Aharonov–Bohm topological phase jump.

On Langmuir circulation in 1 : 2 and 1 : 3 resonance

This paper is concerned with the nonlinear dynamics of spanwise periodic longitudinal vortex modes (Langmuir circulation (LC)) that arise through the instability of two-dimensional periodic flows (waves) in a non-stratified uniformly sheared layer of finite depth. Of particular interest is the excitation of the vortex modes either in the absence of interaction or in resonance, as described by nonlinear amplitude equations built upon the mean field Craik–Leibovich (CL) equations. Since Y-junctions in the surface footprints of Langmuir circulation indicate sporadic increases (doubling) in spacing as they evolve to the scale of sports stadiums, interest is focused on bifurcations that instigate such changes. To that end, surface patterns arising from the linear and nonlinear excitation of the vortex modes are explored, subject to two parameters: a Rayleigh number ${\mathcal{R}}$ present in the CL equations and a symmetry breaking parameter $\unicode[STIX]{x1D6FE}$ in the mixed free surface boundary conditions that relax to those at the layer bottom where $\unicode[STIX]{x1D6FE}=0$. Looking first to linear instability, it is found as $\unicode[STIX]{x1D6FE}$ increases from zero to unity, that the neutral curves evolve from asymmetric near onset to almost symmetric. The nonlinear dynamics of single modes is then studied via an amplitude equation of Ginzburg–Landau type. While typically of cubic order when the bifurcation is supercritical (as it is here) and the neutral curves are parabolic, the Ginzburg–Landau equation must instead here be of quartic order to recover the asymmetry in the neutral curves. This equation is then subjected to an Eckhaus instability analysis, which indicates that linearly unstable subharmonics mostly reside outside the Eckhaus boundary, thereby excluding them as candidates for excitation. The surface pattern is then largely unchanged from its linear counterpart, although the character of the pattern does change when $\unicode[STIX]{x1D6FE}\ll 1$ as a result of symmetry breaking. Attention is then turned to strong resonance between the least stable linear mode and a sub-harmonic of it, as described by coupled nonlinear amplitude equations of Stuart-Landau type. Both 1 : 2 and 1 : 3 resonant interactions are considered. Phase plots and bifurcation diagrams are employed to reveal classes of solution that can occur. Dominant over much of the ${\mathcal{R}}$-$\unicode[STIX]{x1D6FE}$ range considered are non-travelling pure- and mixed-mode equilibrium solutions that act singly or together. To wit, pure modes solutions alone act to realise windrows with spacings in accord with linear theory, while bistability can realise Y-junctions and, depending upon initial conditions, double or even triple the dominant spacing of LC.

Reply to “Comments on ‘A Wave-Resolving Simulation of Langmuir Circulations with a Nonhydrostatic Free-Surface Model: Comparison with Craik–Leibovich Theory and an Alternative Eulerian View of the Driving Mechanism’”

AbstractFujiwara et al. explicitly simulated Langmuir circulations using a wave-resolving simulation (WRS) technique and found that the residual wave effect on vorticity was well represented by the vortex force of the Craik–Leibovich (CL) equation, at least in the simulated situation. In response to the simulation results, Mellor has proposed a view that ubiquitous applicability of the CL formulation is still questionable and that the three-dimensional radiation stress (3DRS) formulation that he has derived encompasses both of the vortex force effect and an effect that is lower order in terms of wave steepness. Here, these opinions are discussed in terms of the approximations used in the wave-averaged formulations. The asymptotic expansion of the Eulerian-averaged momentum equation allows the separate discussion of two different wave effects: pressure correction and torque. It is argued that the approximation adopted in Mellor’s 3DRS formulation is presumably not accurate enough to properly parameterize the wave torque effect, and possible approaches to examine its performance are proposed. We agree with the view that the applicability of the CL formulation needs further investigation. WRS will be a helpful tool for this purpose.

Comments on “A Wave-Resolving Simulation of Langmuir Circulations with a Nonhydrostatic Free-Surface Model: Comparison with Craik–Leibovich Theory and an Alternative Eulerian View of the Driving Mechanism”

AbstractThe results of the subject paper are reviewed wherein credible Langmuir cells are produced by a numerical solution of the primitive fluid dynamic equations with a free surface. Whereas it is a major achievement, the claim that the same results support the general application of the so-called vortex force equations is challenged.

An inviscid analysis of the Prandtl azimuthal mass transport during swirl-type sloshing

An inviscid analytical theory of a slow steady liquid mass rotation during the swirl-type sloshing in a vertical circular cylindrical tank with a fairly deep depth is proposed by utilising the asymptotic steady-state wave solution by Faltinsen et al. (J. Fluid Mech., vol. 804, 2016, pp. 608–645). The tank performs a periodic horizontal motion with the forcing frequency close to the lowest natural sloshing frequency. The azimuthal mass transport (first observed in experiments by Prandtl (Z. Angew. Math. Mech., vol. 29(1/2), 1949, pp. 8–9)) is associated with the summarised effect of a vortical Eulerian-mean flow, which, as we show, is governed by the inviscid Craik–Leibovich equation, and an azimuthal non-Eulerian mean. Suggesting the mass-transport velocity tends to zero when approaching the vertical wall (supported by existing experiments) leads to a unique non-trivial solution of the Craik–Leibovich boundary problem and, thereby, gives an analytical expression for the summarised mass-transport velocity within the framework of the inviscid hydrodynamic model. The analytical solution is validated by comparing it with suitable experimental data.

Lagrangian Investigation of Wave-Driven Turbulence in the Ocean Surface Boundary Layer

AbstractTurbulent 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.