A Lagrangian model for boundary layer growth and bed shear stress in the swash zone

A new model for the boundary layer development and associated skin friction coefficients and shear stress within the swash zone is presented. The model is developed within a Lagrangian reference frame, following fluid trajectories, and can be applied to both laminar flow and smooth turbulent flow. The model is based on the momentum integral approach for steady, flat-plate boundary layers, with appropriate modifications to account for the unsteady flow regime and flow history. The model results are consistent with previous measurements of bed shear stress and skin friction coefficients within the swash zone. These indicate strong temporal and spatial variation throughout the swash cycle, and a clear distinction between the uprush and backwash phase. This variation has been previously attributed the unsteady flow regime and flow history effects, both of which are accounted for in the new model. Fluid particle trajectories and velocity are computed using the non-linear shallow water wave equations and the boundary layer growth across the entire swash zone is estimated. Predictions of the bed shear stress and skin friction coefficients agree reasonably well with direct bed shear stress measurements reported by Barnes et al. (Barnes, M.P., O’Donaghue, T., Alsina, J.M., Baldock, T.E., 2009. Direct bed shear stress measurements in bore-driven swash. Coastal Engineering 56 (8), 853–867) and, for a given flow velocity, give stresses which are consistent with the bias toward uprush sediment transport which has consistently been observed in measurements. The data and modelling suggest that the backwash boundary layer is initially laminar, which results in the late development of significant bed shear during the backwash, with a transition to a turbulent boundary layer later in the backwash. A new conceptual model for the boundary layer structure at the leading edge of the swash is proposed, which accounts for both the no-slip condition at the bed and the moving wet–dry interface. However, further development of the Lagrangian Boundary Layer Model is required in order to include bore-generated turbulence and to account for variable roughness and mobile beds.     

Flow under standing waves: Part 1. Shear stress distribution,energy flux and steady streaming

The conditions for energy flux, momentum flux and the resulting streaming velocity are analysed for standing waves formed in front of a fully reflecting wall. The exchange of energy between the outer wave motion and the near bed oscillatory boundary layer is considered, determining the horizontal energy flux inside and outside the boundary layer. The momentum balance, the mean shear stress and the resulting time averaged streaming velocities are determined. For a laminar bed boundary layer the analysis of the wave drift gives results similar to the original work of Longuet–Higgins from 1953. The work is extended to turbulent bed boundary layers by application of a numerical model. The similarities and differences between laminar and turbulent flow conditions are discussed, and quantitative results for the magnitude of the mean shear stress and drift velocity are presented. Full two-dimensional simulations of standing waves have also been made by application of a general purpose Navier–Stokes solver. The results agree well with those obtained by the boundary layer analysis. Wave reflection from a plane sloping wall is also investigated by using the same numerical model and by physical laboratory experiments. The phase shift of the reflected wave train is compared with theoretical and empirical models.     

A diapycnal diffusion algorithm for isopycnal ocean circulation models with special application to mixed layers

An algorithm is presented for solving the one-dimensional diffusion equation for density, written in terms of density (or a like surrogate) as the independent variable. The algorithm maintains nonnegative layer thicknesses, the premise of the transformation to density as the independent coordinate, under certain restrictions. Near-zero thickness layers can be maintained at the boundaries to accommodate future inflation in response to heating from the boundary. Layers can shrink to near-zero thickness in response to cooling from the boundary. A slight modification of the algorithm permits layers to have diffusion coefficients which differ by orders of magnitude. This provides a natural framework for a surface mixed layer in an isopycnal model, in which the mixed layer is distinguished as a zone of very high turbulent diffusivity overlying an ocean interior of much smaller turbulent diffusivity. The “mixed layer” may be an aggregation of several isopycnal layers rather than just one. A substantial jump in density at the mixed layer base can be represented by several near-zero thickness isopycnal layers. The specification of the thickness of the mixing zone, i.e., the mixed layer depth, is external to the algorithm. An illustration is given using a Kraus–Turner-type specification.     

Inleakage of a vortex-ring bunch onto a plane screen

A method is suggested for simulating axisymmetric laminar or turbulent flows formed during the motion of a vortex-ring bunch of given geometry and circulation toward a plane screen. Earlier, similar problems were simulated with the numerical solution of the Navier-Stokes equations for laminar flows. Turbulent flows have remained unconsidered until now. When a vortex ring approaches the screen, the secondary nonstationary flow is induced near the screen’s surface and this secondary flow causes the formation of the radial boundary layer (provided that air viscosity is taken into account). First, the medium spreads out from the critical point at the screen’s center with the negative pressure gradient along the radial coordinate and then detaches in the region of the positive pressure gradient. This radial wall flow and the corresponding boundary layer are considered in the quasi-stationary approximation. When the boundary layer detaches at successive instances, the flow is replenished with the radially moving secondary vortex rings whose circulations have the sign opposite to that of the circulation of the primary vortex ring. It is the interaction of the primary and secondary vortices that governs process dynamics, which differs substantially from that in the case when the formation of secondary vortices is disregarded. The suggested method is based on the method of discrete vortices (a perfect liquid) and the boundary-layer (laminar or turbulent) theory. During the development of the flow under investigation, the nonstationary ascending flow in the direction perpendicular to the screen’s plane is formed and then this flow decays and dissipates. Simulations for large Reynolds numbers corresponding to the formation of the turbulent boundary layer show that the velocity of ascending vortices in the plane of the initial vortex bunch is less than one-tenth of the initial velocity of the descending vortex ring. The boundary layer is introduced into calculations with the sole goal of determining the parameters of the secondary vortex rings formed during boundary-layer detachments. The interaction of the primary and secondary vortices is then considered within the framework of a perfect medium. Simulations for large Reynolds numbers corresponding to the formation of the turbulent boundary layer on the screen were correlated with the available data obtained in laboratory experiments for small Reynolds numbers. Qualitative agreement between the simulations and experiments is fairly satisfactory. The simulation for one combination of the circulation and vortex-ring geometry takes at most 10–15 min with the use of an average PC.     

Modification of the damping function in the k–ε model to analyse oscillatory boundary layers

A simple relationship has been developed between the wall coordinate y⁺ and Kolmogorov's length scale using direct numerical simulation (DNS) data for a steady boundary layer. This relationship is then utilized to modify two popular versions of low Reynolds number k–ε model. The modified models are used to analyse a transitional oscillatory boundary layer. A detailed comparison has been made by virtue of velocity profile, turbulent kinetic energy, Reynolds stress and wall shear stress with the available DNS data. It is observed that the low Reynolds number models used in the present study can predict the boundary layer properties in an excellent manner.     

Measurements of the wall pressure spectra on a full-scale experimental towed array

Turbulent wall pressure data acquired during tests of a full-scale experimental towed array over a range of tow speeds in straight tows and turns is presented. The experimental towed array contained a linear array of sensors mounted at the fluid–solid interface to measure the spectra of the wall pressure fluctuations due to the cylindrical turbulent boundary layer. The physics are dominated by the growth of a thick, high Reynolds number turbulent boundary layer at arc length Reynolds numbers as high as 9×10⁸. The measured wavenumber-frequency spectra, autospectra, cross-spectral decay and convection velocities are presented. A well-defined convective ridge exists in the wavenumber-frequency spectra obtained during straight tows and turns. Turns give rise to a complicated fluid–structure interaction problem, but do not lead to the separation of the turbulent boundary layer. As the array moves through a turn, flow-induced vibrations of the array are shown to dominate the spectra at low frequencies, with more rapid decay in the measured coherence occurring at higher frequencies. The use of tow speed as a velocity scale is shown to collapse autospectra and convection velocities.     

A new approach to oscillatory rough turbulent boundary layers

Velocity measurements have been performed in an oscillatory turbulent boundary layer over a rough wall, using a large oscillating water tunnel. These together with measurements by Kalkanis (1964) over an oscillating wall indicate the existence of universal wall and defect laws for velocity. A logarithmic overlap layer is predicted and observed as in a steady turbulent boundary layer, and this results in a new relationship between friction factor and relative boundary layer thickness. The phase lead of the defect velocity relative to the wall ditto seems to follow a universal law over the whole defect layer. A method is suggested for the calculation of the phase lead of wall shear stress over velocity in the free stream for large amplitude to roughness ratios. Apart from the inner layer, it is in principle possible to construct the velocity profiles in a turbulent oscillatory boundary layer at a rough wall, using the findings of this report. A review of experimental and theoretical investigations of the stability of the oscillatory boundary layer is also given.     

Reynolds number variation in oscillatory boundary layers: Part I. Purely oscillatory motion

A numerical model based upon a low Reynolds number turbulence closure is proposed to study Reynolds number variation in reciprocating oscillatory boundary layers. The model is used to compute the boundary layer for flow regimes ranging from smooth laminar to rough turbulent. Criteria for fully developed turbulence are derived for walls of the smooth and rough types. In particular, a new criterion to identify the rough turbulent regime is determined based on the time-averaged turbulence intensity. The reliability of the present model is assessed through comparisons with detailed experimental data collected by other investigators. The model globally improves upon standard high Reynolds number closures. Variation through the wave cycle of the main flow variables (ensemble-averaged velocity, shear stress, turbulent kinetic energy) is remarkably well-predicted for smooth walls. Predictions are satisfactory for rough walls as well. Yet, the turbulence level in the rough turbulent regime is overpredicted in the vicinity of the bed.     

Boundary layers and turbulence

We analyze the results of investigation of turbulent boundary layers typical of geophysical objects. It is shown that boundary layers of various nature are self-regulating sustems characterized by relatively slow evolutionary processes of formation accompanied by the growth of instabilities of different types and then replaced by the rapid development of instabilities and destruction of the boundary layers. This cycle is repeated with a certain quasiregular frequency. The destruction of the boundary layer is accompanied by the ejection of turbulent structures whose parameters are characterized by stable experimentally reproducible mean values. This mechanism is responsible for the process of exchange in the boundary layers. Translated by Peter V. Malyshev and Dmitry V. Malyshev     

Features of turbulent momentum and heat transfer in a stably stratified boundary layer over a rough surface

The features of the structure of a stable boundary layer over an urbanized surface are studied using a nonlocal model for the turbulent momentum and heat fluxes, which physically adequately takes into account the effect of buoyancy on turbulent transfer. The transformation of the structure of the boundary layer during transition from a state of convective mixing to a stable state is described by unified expressions for the turbulent momentum and heat fluxes. In some known schemes, different models are used for unstable and stable states. The model reproduces a stable dependence of the Prandtl number on the Richardson number and countergradient heat transfer in a strongly stable boundary layer. The results of numerical simulation are compared to the data of a laboratory experiment and the data obtained using the large-eddy simulation (LES) method.     

Turbulence-induced steady streaming in an oscillating boundary layer: On the reliability of turbulence closure models

The accuracy of several closure models of the Reynolds-Averaged Navier–Stokes Equations in predicting the characteristics of an oscillating turbulent wall boundary layer is analyzed. The analysis involves four low Reynolds number k − ε models and a k − ω model and it is carried out by comparing the model results both with experimental data and with data obtained by a Direct Numerical Simulation (DNS) of the Navier–Stokes equations. The boundary layer is generated by a spatially constant time-oscillating pressure gradient given by the sum of two harmonic components characterized by angular frequencies Ω^∗ and 2Ω^∗ respectively, which generates a steady streaming because of the asymmetry of turbulence intensity during the cycle. Thus the results are relevant to the boundary layer at the bottom of nonlinear sea waves. The attention is therefore focused on the accuracy of the models in reproducing the period averaged profiles of the hydrodynamic characteristics of the steady streaming. The instantaneous quantities, such as time development of the wall shear stress, profiles of the streamwise velocity, Reynolds stresses and turbulent kinetic energy are also considered and analyzed. The results shows that a model can be judged better or worse than other models depending on the specific flow characteristic under investigation. However, an approach has been adopted which allowed to rank the models according to their accuracy in predicting the values of the hydrodynamic quantities involved in the present study.     

The “slip law” of the free surface

A “slip law” connects the excess velocity or “slip” of a wind-blown water surface, relative to the motion in the middle of the mixed layer, to the wind stress, the wind-wave field, and buoyancy flux. An inner layer-outer layer model of the turbulent shear flow in the mixed layer is appropriate, as for a turbulent boundary layer or Ekman layer over a solid surface, allowing, however, for turbulent kinetic energy transfer from the air-side via breaking waves, and for Stokes drift. Asymptotic matching of the velocity distributions in inner and outer portions of the mixed layer yields a slip law of logarithmic form, akin to the drag law of a turbulent boundary layer. The dominant independent variable is the ratio of water-side roughness length to mixed layer depth or turbulent Ekman depth. Convection due to surface cooling is also an important influence, reducing surface slip. Water-side roughness length is a wind-wave property, varying with wind speed similarly to air-side roughness. Slip velocity is typically 20 times water-side friction velocity or 3% of wind speed, varying within a range of about 2 to 4.5%. A linearized model of turbulent kinetic energy distribution shows much higher values near the surface than in a wall layer. Nondimensional dissipation peaks at a value of about eight, a short distance below the surface.     

Mixing processes relevant to phytoplankton dynamics in lakes

Active turbulence in lakes is confined to the surface mixed layer, to boundary layers on the lake sides and bottom, and to turbulent patches in the interior. The density stratification present in most lakes fundamentally alters the pathways connecting external mechanical energy inputs, for example by wind, with its ultimate fate as dissipation to heat; the density stratification supports internal waves and intrusions that distribute the input energy throughout the lake. Intrusions may be viewed as internal waves with zero horizontal wavenumber and are formed each time localised mixing occurs in a stratified fluid. Intrusions are also formed in the epilimnion by differential heating or cooling and by differential deepening. The fraction of lake volume below the diurnal mixed layer that is subject to active turbulence is very small, probably of the order of 1% or less in small to medium‐sized lakes. By contrast, in the surface mixed layer, turbulence is less intermittent and maintains phytoplankton in suspension and controls their exposure to the underwater solar flux. Nutrient transport to individual cells depends not only on the cell Reynolds number but also on the Peclet number, which, if large, implies enhanced mass transfer above purely diffusive values.     

A note on a theoretical model of oscillatory smooth turbulent boundary layers

An analytical theory which describes the motion in an oscillatory smooth turbulent boundary layer using a two-layer time invariant eddy viscosity model is presented. The eddy viscosity in the inner layer increases quadratically with the height above the wall. In the outer layer the eddy viscosity is taken as a constant.     

Effect of surface wave breaking on the surface boundary layer of temperature in the Yellow Sea in summer

Near-surface enhancement of turbulent mixing and vertical mixing coefficient for temperature owing to the effect of surface wave breaking is investigated using a two-dimensional (2-D) ocean circulation model with a tidal boundary condition in an idealized shelf sea. On the basis of the 2-D simulation, the effect of surface wave breaking on surface boundary layer deepening in the Yellow Sea in summer is studied utilizing a 3-D ocean circulation model. A well-mixed temperature surface layer in the Yellow Sea can be successfully reconstructed when the effect of surface wave breaking is considered. The diagnostic analysis of the turbulent kinetic energy equation shows that turbulent mixing is enhanced greatly in the Yellow Sea in summer by surface wave breaking. In addition, the diagnostic analysis of momentum budget and temperature budget also show that surface wave breaking has an evident contribution to the turbulent mixing in the surface boundary layer. We therefore conclude that surface wave breaking is an important factor in determining the depth of the surface boundary layer of temperature in the Yellow Sea in summer.     

Axisymmetric viscoelastic response of flexible pipes in time domain

The challenges for determining the mechanical behavior of flexible pipes mainly arise from highly non-linear geometrical and material properties and complex contact interaction conditions between and within layers components. This paper develops an innovative model to investigate the linear viscoelastic behavior of flexible pipes under axisymmetric loads in time domain. The model is derived from an equivalent linear elastic axisymmetric model by invoking the elastic-viscoelastic correspondence principle. Analytical formulations that describe the behavior of the metallic helical layers based on a combination of differential geometry concepts and Clebsch–Kirchhoff equilibrium equations for initially curved slender elastic rods are presented. The elastic response of the homogenous polymeric cylindrical layers is also presented. The assemblage of both types of governing algebraic equations that approximate analytical solutions for force and moment distributions, deformations in each layer, as well as contact pressure between near layers, taking time-dependent characteristics of polymeric layers into account are provided and it is clear that the relationship between axial force and elongation is non-linear and encompasses a hysteretic response. Besides, the creep behavior in axial direction can also be found. Some insights into the differences in the behavior for several loading conditions are discussed by considering variable frequencies.     

Study of the Neutral Ekman Flow Using an Algebraic Reynolds Stress Model

A recently developed fully explicit algebraic model of Reynolds stress and turbulent heat flux in a thermally stratified planetary atmospheric boundary layer without stratification has been used for a numerical study of the Ekman turbulent boundary layer over a homogeneous rough surface for different dimensionless surface Rossby numbers. A comparative analysis has been conducted for a closure model of the transport term in the prognostic equation of turbulent kinetic energy dissipation including third-order moments. Dependences of the total wind rotation angle on the Rossby number have been obtained. The calculated vertical profiles of mean velocity, turbulent stress, turbulent kinetic energy, surface-friction velocity, and boundary-layer height agree satisfactorily with observational and earlier obtained LES data.     

Modeling wave�Ccurrent bottom boundary layers beneath shoaling and breaking waves

The boundary layer characteristics beneath waves transforming on a natural beach are affected by both waves and wave-induced currents, and their predictability is more difficult and challenging than for those observed over a seabed of uniform depth. In this research, a first-order boundary layer model is developed to investigate the characteristics of bottom boundary layers in a wave–current coexisting environment beneath shoaling and breaking waves. The main difference between the present modeling approach and previous methods is in the mathematical formulation for the mean horizontal pressure gradient term in the governing equations for the cross-shore wave-induced currents. This term is obtained from the wave-averaged momentum equation, and its magnitude depends on the balance between the wave excess momentum flux gradient and the hydrostatic pressure gradient due to spatial variations in the wave field of propagating waves and mean water level fluctuations. A turbulence closure scheme is used with a modified low Reynolds number k-ε model. The model was validated with two published experimental datasets for normally incident shoaling and breaking waves over a sloping seabed. For shoaling waves, model results agree well with data for the instantaneous velocity profiles, oscillatory wave amplitudes, and mean velocity profiles. For breaking waves, a good agreement is obtained between model and data for the vertical distribution of mean shear stress. In particular, the model reproduced the local onshore mean flow near the bottom beneath shoaling waves, and the vertically decreasing pattern of mean shear stress beneath breaking waves. These successful demonstrations for wave–current bottom boundary layers are attributed to a novel formulation of the mean pressure gradient incorporated in the present model. The proposed new formulation plays an important role in modeling the boundary layer characteristics beneath shoaling and breaking waves, and ensuring that the present model is applicable to nearshore sediment transport and morphology evolution.     

Boundary layer approach in the modeling of breaking solitary wave runup

The boundary layer is very important in the relation between wave motion and bed stress, such as sediment transport. It is a known fact that bed stress behavior is highly influenced by the boundary layer beneath the waves. Specifically, the boundary layer underneath wave runup is difficult to assess and thus, it has not yet been widely discussed, although its importance is significant. In this study, the shallow water equation (SWE) prediction of wave motion is improved by being coupled with the k–ω model, as opposed to the conventional empirical method, to approximate bed stress. Subsequently, the First Order Center Scheme and Monotonic Upstream Scheme of Conservation Laws (FORCE MUSCL), which is a finite volume shock-capturing scheme, is applied to extend the SWE range for breaking wave simulation. The proposed simultaneous coupling method (SCM) assumes the depth-averaged velocity from the SWE is equivalent to free stream velocity. In turn, free stream velocity is used to calculate a pressure gradient, which is then used by the k–ω model to approximate bed stress. Finally, this approximation is applied to the momentum equation in the SWE. Two experimental cases will be used to verify the SCM by comparing runup height, surface fluctuation, bed stress, and turbulent intensity values. The SCM shows good comparison to experimental data for all before-mentioned parameters. Further analysis shows that the wave Reynolds number increases as the wave propagates and that the turbulence behavior in the boundary layer gradually changes, such as the increase of turbulent intensity.