Direct bed shear stress measurements in bore-driven swash

Direct measurements of bed shear in the swash zone are presented. The data were obtained using a shear plate in medium and large-scale laboratory bore-driven swash and cover a wide range of bed roughness. Data were obtained across the full width of the swash zone and are contrasted with data from the inner surf zone. Estimates of the flow velocities through the full swash cycle were obtained through numerical modelling and calibrated against measured velocity data. The measured stresses and calculated flow velocities were subsequently used to back-calculate instantaneous local skin friction coefficients using the quadratic drag law. The data show rapid temporal variation of the bed shear stress through the leading edge of the uprush, which is typically two–four times greater than the backwash shear stresses at corresponding flow velocity. The measurements indicate strong temporal variation in the skin friction coefficient, particularly in the backwash. The general behaviour of the skin friction coefficient with Reynolds number is consistent with classical theory for certain stages of the swash cycle. A spatial variation in skin friction coefficient is also identified, which is greatest across the surf-swash boundary and likely related to variations in local turbulent intensities. Skin friction coefficients during the uprush are approximately twice those in the backwash at corresponding Reynolds number and cross-shore location. It is suggested that this is a result of the no-slip condition at the tip leading to a continually developing leading edge and boundary layer, into which high velocity fluid and momentum are constantly injected from the flow behind and above the tip region. Finally, the measured stress data are used to determine the asymmetry and cross-shore variation in potential sediment transport predicted by three forms of sediment transport formulae.     

RANS v2–f simulation of a swash event: Detailed flow structure

The flow structure of a swash event over a uniform slope is studied using a RANS-VOF numerical model coupled with a v²–f turbulence closure. The model is compared with experimental data of recent laboratory experiments. The ability of the turbulence modelling for simulating swash flow and the evolution of the computed bed shear stress during run-up and run-down are investigated. The agreement between numerical results and measured data, such as water depth, depth-averaged velocity and bed shear stress is very good during run-up. Main discrepancies are found during run-down. The paper also examines the aeration of the water layer in the swash flow, taking advantage of the PLIC method for computing the air–water interfaces. Air is continuously entrapped in the swash front and released at its rear during run-up. A detailed analysis indicates that the flow reversal is initiated near the bottom at the outer boundary of the swash zone and progresses landward. The study highlights the asymmetry between run-up and run-down. During run-up, the swash front propagation determines the turbulence properties and the bed shear stress profile on the beach, whereas the flow properties are more homogeneously distributed in the swash area during run-down.     

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.     

Suspended sediment transport in the swash zone of a dissipative beach

Simultaneous high frequency field measurements of water depth, flow velocity and suspended sediment concentration were made at three fixed locations across the high tide swash and inner surf zones of a dissipative beach. The dominant period of the swash motion was 30–50 s and the results are representative of infragravity swash motion. Suspended sediment concentrations, loads and transport rates in the swash zone were almost one order of magnitude greater than in the inner surf zone. The vertical velocity gradient near the bed and the resulting bed shear stress at the start of the uprush was significantly larger than that at the end of the backwash, despite similar flow velocities. This suggests that the bed friction during the uprush was approximately twice that during the backwash.The suspended sediment profile in the swash zone can be described reasonably well by an exponential shape with a mixing length scale of 0.02–0.03 m. The suspended sediment transport flux measured in the swash zone was related to the bed shear stress through the Shields parameter. If the bed shear stress is derived from the vertical velocity gradient, the proportionality coefficient between shear stress and sediment transport rate is similar for the uprush and the backwash. If the bed shear stress is estimated using the free-stream flow velocity and a constant friction factor, the proportionality factor for the uprush is approximately twice that of the backwash. It is suggested that the uprush is a more efficient transporter of sediment than the backwash, because the larger friction factor during the uprush causes larger bed shear stresses for a given free-stream velocity. This increased transport competency of the uprush is necessary for maintaining the beach, otherwise the comparable strength and greater duration of the backwash would progressively remove sediment from the beach.     

Combined tidal and wind driven flows and bedload transport over a flat bottom

The combined tidal and wind driven flow and resulting sediment transport in the ocean over a flat bottom at intermediate water depth has been investigated, using a simple one dimensional two-equation turbulence closure model. This model has been verified against field measurements of a tidal flow in the Celtic Sea. The tidal velocity ellipses and the time series of the horizontal velocity components at given elevations above the bottom are well predicted through the water column although there are some deviations between the predicted and measured velocities near the bottom due to the uncertainty of the bottom roughness. For the combined tidal and wind driven flows the velocity profiles, turbulent kinetic energy profiles and surface particle trajectories are predicted for weak and strong winds. Furthermore, the bottom shear stress and the resulting bedload transport have been predicted; the parts of the particle trajectories in the close vicinity of the bottom where the bedload transport exists are displayed. Finally, the direction and magnitude of the surface drift, the depth-averaged mean velocity and the mean bedload transport are given, and the effect of the bottom roughness on the sea surface drift is investigated.     

Validation of volume continuity method for estimation of cross-shore swash flow velocity

The use of flow volume continuity to estimate depth-averaged cross-shore swash flow velocity is validated using appropriate field data. Swash surface and beach face elevation were measured by a network of ultrasonic altimeters mounted through the swash zone of a sandy beach. Estimates of cross-shore flow velocity derived from these data compare well with concurrent measurements obtained from a vertical array of four electromagnetic current meters (EMCMs). Significantly, the continuity technique allows swash velocity to be estimated throughout the full swash cycle. This is not possible using fixed current meters and the new measurements indicate that swash velocity is more negatively skewed than previously observed.     

一个湍能封闭模型对风暴潮海底底应力的应用研究

将一个三维湍能封闭模型应用于开阔海区的风暴潮，通过数值计算探讨了Ｔａｙｌｏｒ底摩擦二次率的拖曳系数随空间的分布及拖曳系数与水深、海底粗糙度、风向和风速等因素的关系。本文对底摩擦二次率的可靠性做了评价。     

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.     

剪切流中界面波的二阶Stokes 波解

考虑了均匀剪切流场中强非线性界面波,建立了基于任意水深处速度而不是通常所用的平均速度为速度变量的模型,分析了其色散关系,并求得了各层速度的二阶渐近解和界面内波波面位移的二阶Stokes解,揭示了波流之间的非线性相互作用和界面波解之间的非线性相互作用。     

基于FVCOM模式钱塘江河口涌潮三维数值模拟研究

应用FVCOM模型进行钱塘江涌潮的二维和三维数值模拟,并使用2007年10月大潮期间的现场观测结果进行模型验证。结果表明,平面二维模拟结果与实测值吻合较好,与基于Godunov算法求解钱塘江涌潮二维特征的结果基本一致。三维计算结果再现了流速在垂向上的差异,以及在落潮转向涨潮时刻,涨潮先从底部开始,垂向上从底层向表层逐步从落潮流转为涨潮流。敏感性分析显示模型网格大小和曼宁系数的选取对模拟结果有较大的影响。     

The influence of seaward boundary conditions on swash zone hydrodynamics

The influence of the seaward boundary condition on the internal swash hydrodynamics is investigated. New numerical solutions of the characteristics form of the nonlinear shallow-water equations are presented and applied to describe the swash hydrodynamics forced by breaking wave run-up on a plane beach. The solutions depend on the specification of characteristic variables on the seaward boundary of the swash zone, equivalent to prescribing the flow depth or the flow velocity. It is shown that the analytical solution of Shen and Meyer [Shen, M.C., Meyer, R.E., 1963. Climb of a bore on a beach. Part 3. Runup. J. Fluid Mech. 16, 113–125] is a special case of the many possible solutions that can describe the swash flow, but one that does not appear appropriate for practical application for real waves. The physical significance of the boundary conditions is shown by writing the volume and momentum fluxes in terms of the characteristic variables. Results are presented that illustrate the dependence of internal flow depth and velocity on the boundary condition. This implies that the internal swash hydrodynamics depend on the shape and wavelength of the incident bore, which differs from the hydrodynamic similarity inherent in the analytical solution. A solution appropriate for long bores is compared to laboratory data to illustrate the difference from the analytical solution. The results are important in terms of determining overwash flows, flow forces and sediment dynamics in the run-up zone.     

均匀剪切流场中的强非线性波及其色散关系

考虑了剪切流场中强非线性波,唯一的假设是水深与特征波长之比是小量,建立了基于任意水深处速度而不是通常所用的平均速度为速度变量的模型。不仅改进了色散关系,使模型的水深适用范围更大,而且由于整个推导过程对波的振幅没有做任何假设,因而所获得的模型可以用于任何振幅的波。     

The role of bore collapse and local shear stresses on the spatial distribution of sediment load in the uprush of an intermediate-state beach

Spatial variations in sediment load in the swash uprush and textural properties of sediment in transport were evaluated to investigate the mechanisms responsible for sediment transport during wave uprush. Four streamer traps were deployed at 2.0-m intervals across the swash zone of a sheltered, microtidal sandy beach at Port Beach, Western Australia, over a 4-day period. During these trapping experiments, offshore significant wave heights were 0.3–0.5 m and wave periods were about 10 s. The average width of the uprush zone was 6.9 m and the average uprush duration was 5.9 s. Cross-shore distributions of sediment load for 70 uprush events reveal a maximum in sediment load landward of the base of the swash (at about 20% of swash width) during single events and a maximum closer to mid-swash (at about 40% of swash width) during multiple events characterized by swash interactions. Settling velocity distributions of trap samples during individual uprush events are similar to distributions found on the beach surface, with the lowest settling velocities (finest sediments) near the base of the swash zone and maximum settling velocities (coarsest sediments) around the mid-swash position. It was found that sediment transport during wave uprush occurs through two distinct mechanisms: (1) sediment entrainment during bore collapse seaward of the base of the swash zone and subsequent advection of this bore-entrained sediment up the beach by wave uprush; and (2) in situ sediment entrainment and transport induced by local shear stresses during wave uprush. Both mechanisms are considered important, but the first mechanism is considered most significant during the early stages of wave uprush when sediment is transported mainly in suspension, while the second mechanism is likely to dominate the mid- to later stages of wave uprush when sediment is transported mainly by sheet flow. The relative importance of the two mechanisms will vary between different beaches with the morphodynamic state of the beach (reflective versus dissipative) expected to play a major role.     

Experimental study of turbulent oscillatory boundary layers in an oscillating water tunnel

A high-quality experimental study including a large number of tests which correspond to full-scale coastal boundary layer flows is conducted using an oscillating water tunnel for flow generations and a Particle Image Velocimetry system for velocity measurements. Tests are performed for sinusoidal, Stokes and forward-leaning waves over three fixed bottom roughness configurations, i.e. smooth, “sandpaper” and ceramic-marble bottoms. The experimental results suggest that the logarithmic profile can accurately represent the boundary layer flows in the very near-bottom region, so the log-profile fitting analysis can give highly accurate determinations of the theoretical bottom location and the bottom roughness. The first-harmonic velocities of both sinusoidal and nonlinear waves, as well as the second-harmonic velocities of nonlinear waves, exhibit similar patterns of vertical variation. Two dimensionless characteristic boundary layer thicknesses, the elevation of 1% velocity deficit and the elevation of maximum amplitude, are found to have power-law dependencies on the relative roughness for rough bottom tests. A weak boundary layer streaming embedded in nonlinear waves and a small but meaningful third-harmonic velocity embedded in sinusoidal waves are observed. They can be only explained by the effect of a time-varying turbulent eddy viscosity. The measured period-averaged vertical velocities suggest the presence of Prandtl's secondary flows of the second kind in the test channel. Among the three methods to infer bottom shear stress from velocity measurements, the Reynolds stress method underestimates shear stress due to missed turbulent eddies, and the momentum integral method also significantly underestimates bottom shear stress for rough bottom tests due to secondary flows, so only the log-profile fitting method is considered to yield the correct estimate. The obtained bottom shear stresses are analyzed to give the maximum and the first three harmonics, and the results are used to validate some existing theoretical models.     

Measurement and modeling of bed shear stress under solitary waves

Direct measurements of bed shear stresses (using a shear cell apparatus) generated by non-breaking solitary waves are presented. The measurements were carried out over a smooth bed in laminar and transitional flow regimes (~ 10⁴ < R_e < ~ 10⁵). Measurements were carried out where the wave height to water depth (h/d) ratio varied between 0.12 and 0.68; maximum near bed velocity varied between 0.16 m/s and 0.51 m/s and the maximum total shear stress (sum of skin shear stress and Froude–Krylov force) varied between 0.386 Pa and 2.06 Pa. The total stress is important in determining the stability of submarine sediment and in sheet flow regimes. Analytical modeling was carried out to predict total and skin shear stresses using convolution integration methods forced with the free stream velocity and incorporating a range of eddy viscosity models. Wave friction factors were estimated from skin shear stress at different instances over the wave (viz., time of maximum positive total shear stress, maximum skin shear stress and at the time of maximum velocity) using both the maximum velocity and the instantaneous velocity at that phase of the wave cycle. Similarly, force coefficients obtained from total stress were estimated at time of maximum positive and negative total stress and at maximum velocity. Maximum positive total shear stress was approximately 1.5 times larger than minimum negative total stress. Modeled and measured positive bed shear stresses are well correlated using the best convolution model, but the model underestimates the data by about 4%. Friction factors are dependent on the choice of normalizing using the maximum velocity, as is conventional, or the instantaneous velocity. These differ because the stress is not in phase with the velocity in general. Friction factors are consistent with previous data for monochromatic waves, and vary inversely with the square-root of the Reynolds number. The total shear stress leads the free stream fluid velocity by approximately 50°, whereas the skin friction shear stress leads by about 30°, which is similar to that reported by earlier researchers.     

Numerical Simulation of Wave Scour Around A Large-Scale Circular Cylinder

A numerical model is developed for estimation of local scour around a large circular cylinder under vvave action. The model includes wave diffraction around structures, bed shear stress calculation inside the vvave boundary layer and topo-graphical change model. The vvave model is based on the improved Boussinesq equations for varying depth. The vvave boundary layer is calculaled by solving the integrated momentum equation over the boundary layer. The bed shear stress due to streaming, an important factor affecting the sediment transport around a large-scale cylinder, is calculated. The Lagrangian drift velocity is included in calculation of the suspended sediment transport rates. The model is implemented by a finite element method and the results from the present model, which agree well with experimental data, are com-pared vvith those from other methods.     

On the transport of suspended sediment by a swash event on a plane beach

We develop solutions for the transport of suspended sediment by a single swash event following the collapse of a bore on a plane beach, and we investigate the morphodynamical role that such transport may play. Although the intrinsic asymmetry between uprush and backwash velocities tends to encourage the export of sediment, we find that swash events may be effective in distributing across the swash zone much or all of the sediment mobilised by bore collapse; additionally, settling lag effects may promote a weak onshore movement of sediment. We quantify both effects in terms of the properties of the sediment and of the swash event, and comment on the relationship between our findings and recent field studies of swash zone sediment transport.