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.     

Measurement of wave-by-wave bed-levels in the swash zone

A technique is described to observe and quantify wave-by-wave bed-level changes in the swash zone. The ultrasonic instrument system is non-contact with the beach face surface being measured and the sensors remain outside of the fluid flows causing sediment movement. Sensor resolution combined with the electronic noise inherent within a digital network data-logging system results in a (conservative) measurement accuracy of ± 1 mm, equating to a couple of sand grain diameters in height. Illustrative field results demonstrate the practical use of the instrumentation, and a simple data pre-processing method to separate swashes and intervening bed-level ‘events’ is discussed. These example data reveal rather complex fluctuations of the bed observed over time periods of minutes to hours. Rather strikingly, gross bed-level changes per wave are revealed to be up to many orders of magnitude larger than the observed net rate of beach face evolution. It is outlined how observations of successive bed-level changes at multiple locations within a dense grid, combined with a consideration of sediment continuity, will now enable the total net sediment transported per uprush–backwash to be quantified.     

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.     

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.     

A numerical study of swash flows generated by bores

The dynamic processes of bore propagation over a uniform slope are studied numerically using a 2-D Reynolds Averaged Navier–Stokes (RANS) solver, coupled to a non-linear k − ε turbulence closure and a volume of fluid (VOF) method. The dam-break mechanism is used to generate bores in a constant depth region. Present numerical results for the ensemble-averaged flow field are compared with existing experimental data as well as theoretical and numerical results based on non-linear shallow water (NSW) equations. Reasonable agreement between the present numerical solutions and experimental data is observed. Using the numerical results, small-scale bore behaviors and flow features, such as the bore collapse process near the still-water shoreline, the ‘mini-collapse’ during the runup phase and the ‘back-wash bore’ in the down-rush phase, are described. In the case of a strong bore, the evolution of the averaged turbulence kinetic energy (TKE) over the swash zone consists of two phases: in the region near the still-water shoreline, the production and the dissipation of TKE are roughly in balance; in the region farther landwards of the still-water shoreline, the TKE decay rate is very close to that of homogeneous grid turbulence. On the other hand, in the case of a weak bore, the bore collapse generated turbulence is confined near the bottom boundary layer and the TKE decays at a much slower rate.     

Numerical study on wave-induced filtration flow across the beach face and its effects on swash zone sediment transport

A numerical model, coupling an analysis of beach groundwater flow with an analysis of swash wave motion over a uniform slope, is presented. Model calculations are performed to investigate the variations of swash-induced filtration flows across the beach face for different input parameters. Swash zone sediment transport under the influence of such filtration flow across the beach face is investigated through modification of effective weight of sediment particle and modification of swash boundary layer thickness. These effects are quantified based on a bed load transport model with a modified Shields parameter.     

Laboratory study of combined wave overtopping and storm surge overflow of a levee

Combined wave overtopping and storm surge overflow of a levee with a trapezoidal cross section was studied in a two-dimensional laboratory wave/flow flume at a nominal prototype-to-model length scale of 25-to-1. The goal of this study was to develop design guidance in the aftermath of Hurricane Katrina. Time series of water depth at two locations on the levee crown and flow thickness at five locations on the landward-side slope were measured along with horizontal velocity near the landward edge of the crown. New equations are presented for average overtopping discharge, distribution of instantaneous discharge, and distribution of individual wave volumes. Equations are also given for mean flow thickness, RMS wave height, mean velocity, and velocity of the wave front down the landward-side slope.     

Lagrangian measurements and modelling of fluid advection in the inner surf and swash zones

New laboratory and field data are presented on fluid advection into the swash zone. The data illustrate the region of the inner surf zone from which sediment can be directly advected into the swash zone during a single uprush, which is termed the advection length. Experiments were conducted by particle tracking in a Lagrangian reference frame, and were performed for monochromatic breaking waves, solitary bores, non-breaking solitary waves and field conditions. The advection length is normalised by the run-up length to give an advection ratio, A, and different advection ratios are identified on the basis of the experimental data. The data show that fluid enters the swash zone from a region of the inner surf zone that can extend a distance seaward of the bore collapse location that is approximately equal to half of the run-up length. This region is about eight times wider than the region predicted by the classical swash solution of Shen and Meyer [Shen, M.C., Meyer, R.E., 1963. Climb of a bore on a beach. Part 3. Runup. Journal of Fluid Mechanics 16, 113–125], as illustrated by Pritchard and Hogg [Pritchard, D., Hogg, A.J., 2005. On the transport of suspended sediment by a swash event on a plane beach. Coastal Engineering 52, 1–23]. Measured advection ratios for periodic waves show no significant trend with Iribarren number, consistent with self-similarity in typical swash flows. The data are compared to recent characteristic solutions of the non-linear shallow water wave (NLSW) equations and both finite difference and finite volume solutions of the NLSW equations.     

Alongshore fluid motions in the swash zone of a sandy and gravel beach

Field experiments were conducted on a low-gradient, high-energy sandy beach (Truc Vert, France) and a steep, low-energy gravel beach (Slapton, UK) to examine alongshore-directed currents within the swash zone. At Truc Vert, data were collected over 33 tidal cycles with offshore significant wave heights of 1–4 m and periods of 5–12 s. At Slapton data were collected during 12 tides with wave heights of 0.3–1 m and periods of 4–9 s. The swash motion was predominantly at infragravity frequencies at Truc Vert and incident frequencies at Slapton.     

Numerical study of combined overflow and wave overtopping over a smooth impermeable seawall

A numerical wave flume is used to investigate the discharge characteristics of combined overflow and wave overtopping of impermeable seawalls. The numerical procedure computes solutions to the Reynolds-averaged Navier–Stokes equations and includes the generation of an irregular train of waves, the simulation of wave breaking and interaction with a sloping, impermeable wall. The numerical model is first tested against published experimental observations, approximate analytical solutions and empirical design formulae for the cases of pure overflow and pure overtopping. A sequence of numerical experiments simulating combined overflow and overtopping are described. The results are used to determine empirical discharge formulae of the form used in current practice.     

Large-scale experiments on beach profile evolution and surf and swash zone sediment transport induced by long waves,wave groups and random waves

New large-scale laboratory data are presented on the influence of long waves, bichromatic wave groups and random waves on sediment transport in the surf and swash zones. Physical model testing was performed in the large-scale CIEM wave flume at UPC, Barcelona, as part of the SUSCO (swash zone response under grouping storm conditions) experiment in the Hydralab III program (Vicinanza et al., 2010). Fourteen different wave conditions were used, encompassing monochromatic waves, bichromatic wave groups and random waves. The experiments were designed specifically to compare variations in beach profile evolution between monochromatic waves and unsteady waves with the same mean energy flux. Each test commenced with approximately the same initial profile. The monochromatic conditions were perturbed with free long waves, and then subsequently substituted with bichromatic wave groups with different bandwidth and with random waves with varying groupiness. Beach profile measurements were made at half-hourly and hourly intervals, from which net cross-shore transport rates were calculated for the different wave conditions. Pairs of experiments with slightly different bandwidth or wave grouping show very similar net cross-shore sediment transport patterns, giving high confidence to the data set. Consistent with recent small-scale experiments, the data clearly show that in comparison to monochromatic conditions the bichromatic wave groups reduce onshore transport during accretive conditions and increase offshore transport during erosive conditions. The random waves have a similar influence to the bichromatic wave groups, promoting offshore transport, in comparison to the monochromatic conditions. The data also indicate that the free long waves promote onshore transport, but the conclusions are more tentative as a result of a few errors in the test schedule and modifications to the setup which reduced testing time. The experiments suggest that the inclusion of long wave and wave group sediment transport is important for improved near-shore morphological modeling of cross-shore beach profile evolution, and they provide a very comprehensive and controlled series of tests for evaluating numerical models. It is suggested that the large change in the beach response between monochromatic conditions and wave group conditions is a result of the increased significant and maximum wave heights in the wave groups, as much as the presence of the forced and free long waves induced by the groupiness. The equilibrium state model concept can provide a heuristic explanation of the influence of the wave groups on the bulk beach profile response if their effective relative fall velocity is larger than that of monochromatic waves with the same incident energy flux.     

Laboratory experiment and numerical simulation on authigenic mineral formation induced by seabed methane seeps

In natural gas hydrate marine settings, cold seepage of methane fluid is a widely observed phenomenon, where authigenic minerals serve as an indication of potential gas hydrate-bearing reservoirs at depth. In this study, based on the data from the site HD196 near Dongsha Island, northern continental slope of South China Sea, laboratory experiments and numerical simulation studies were conducted to investigate the biogeochemical processes of authigenic mineral formation induced by methane seepage. The bioreactor experimental results show that in response to methane flux, pH increased to 8.5 after 20 days of reaction, and Eh declined rapidly first and then remained unchanged at about 100 mV. The decrease in SO₄²⁻, HS⁻ and HCO₃⁻ concentrations indicated the occurrence of anaerobic oxidation of methane coupled with sulfate reduction (AOM-SR). The depletion of Fe²⁺ implied the formation of iron-bearing minerals, e.g., iron sulfides. Carbonate minerals were also identified in the experimental phase of this study. Most iron sulfides existed as massive pieces, and in some cases as spherical or rod-shape pieces. The calcium carbonates were observed as blocky pieces. Numerical simulations were also performed to reproduce the biogeochemical reactions that occurred in the reactor experiments. Based on experimental data, kinetic parameters associated with the observed reactions were calibrated. The model simulated results are general consistent with those obtained by the experiments conducted in this study. The combination of simulation and experimental studies provided a powerful tool to investigate the biogeochemical processes in the methane leakage environment at different temporal and spatial scales. This study gave a new perspective to understand the formation of cold seep authigenic minerals in marine sediments, and was significant for future investigations on the effects of hydrate decomposition.     

Topographic steering of dense overflows: Laboratory experiments with V-shaped ridges and canyons

Topographic corrugations such as canyons and ridges cross-cutting the path of a dense plume may effectively steer all or part of the plume downslope. Here, topographically steered flows are investigated experimentally, as laminar, dense gravity currents are observed to impinge on and flow along sloping, V-shaped canyons and ridges. Ridges, as well as canyons, were observed to steer the dense water downslope. A dynamical regime, in which the along-slope transport is balanced by a return flow in the Ekman layer to maintain a geostrophically balanced downslope flow along the corrugation, has been proposed. Results from a previously published analytical model describing such flows are compared with the laboratory experiments. The response of the flow to variations in four governing parameters (slope, rotation, volume flux and reduced gravity) is generally described well by the model and results agree qualitatively, although theory slightly underestimates the dense layer thickness. Vertical velocity profiles resolving the Ekman spiral were obtained using a laser Doppler velocimeter and they showed the secondary, transverse circulation superimposed on the primary, downslope flow. A particle flowing down the canyon/along the ridge can be expected to follow a helix-like path, and dye released within the dense layer showed this. The experiments support the analytical model and the dynamical regime proposed for topographically steered flows. The gravity current split in two when the transport capacity of the corrugation was exceeded; one part continued along the slope and the other flowed downslope along the corrugation.     

Field investigation on the results of non-linear shallow water equations in the Swash Zone

For the study of the cross-shore wave-induced hydrodynamics in the swash zone, a numerical model is developed based on the one-dimensional non-linear shallow water (NSW) equations for prediction of hydrodynamic parameters in the swash zone. In order to evaluate the accuracy of the outputs of the numerical model, the model's predictions in terms of water surface elevations and cross-shore velocities, are compared to field data from full-scale experiments conducted on three sites with different beach slope; mild and steep, several bed particle sizes and under various incident wave conditions. The quantitative and qualitative comparison of the results of the numerical model and the full-scale data reveals that the model can generally predict many aspects of the flow in the surf and swash zone on both types of beach. The accuracy is adequate for application in a sediment transport study. Considering the time-history and probability distribution of water surface elevation, the model is generally more accurate on steep beaches than on the mild beach. The model can adequately simulate the dominant frequency across the beach and saturation of higher frequencies on both mild and steep beaches for various incident wave energy characteristics. With regard to the horizontal (cross-shore) velocity, the sawtooth shape of time-history and negative acceleration of water are well predicted by the model for both mild and steep beaches. Due to the uncertainties in maximum and minimum values of velocity data, clear judgement about the accuracy of the numerical model in this matter was not possible. However, the comparison of the minimum velocities (offshore direction) revealed that the application of friction factors below the range which is suggested by literature best match the data.     

Forces on a pitching plate: An experimental and numerical study

A flat plate in pitching motion is considered as a fundamental source of locomotion in the general context of marine propulsion. The experimental as well as numerical investigation is carried out at a relatively small Reynold number of 2000 based on the plate length c and the inflow velocity U_∞. The plate oscillates sinusoidally in pitch about its 1/3 − c axis and the peak to peak amplitude of motion is 20°. The reduced frequency of oscillation k = πfc/U_∞ is considered as a key parameter and it may vary between 1 and 5. The underlying fluid-structure problem is numerically solved using a compact finite-differences Navier–Stokes solution procedure and the numerical solution is compared with Particle Image Velocimetry (PIV) measurements of the flow field around the pitching foil experimental device mounted in a water-channel. A good agreement is found between the numerical and experimental results and the threshold oscillation frequency beyond which the wake exhibits a reverse von Kármán street pattern is determined. Above threshold, the mean velocity in the wake exhibits jet-like profiles with velocity excess, which is generally considered as the footprint of thrust production. The forces exerted on the plate are extracted from the numerical simulation results and it is shown, that reliable predictions for possible thrust production can be inferred from a conventional experimental control volume analysis, only when besides the wake's mean flow the contributions from the velocity fluctuation and the pressure term are taken into account.     

A numerical study on dynamic properties of the gravity cage in combined wave-current flow

Recent work in the area of open-ocean aquaculture-system dynamics has focused either on the response of fish cages in waves or the steady drag response from ocean currents, not on them combined. In reality, however, the forces bearing on these open-ocean structures are a nonlinear, multidirectional combination of both waves and current profiles. In this paper, a numerical model has been developed to simulate the dynamic response of the gravity cage to waves combined with currents. When current flows are combined with regular waves, gravity-cage motion response (including heave, surge, and pitch) and mooring-line forces have been calculated. To examine the validity of simulated results, a series of physical model tests have been carried out. The results of our numerical simulation are all in close agreement with the experimental data.     

Unified numerical study of shallow foundation on structured soft clay under unconsolidated and consolidated-undrained loadings

Abstract

Based on a new elasto-plastic constitutive model, this paper presents a soil–water coupled numerical prediction of the bearing capacity for shallow foundation constructed on Ballina soft clay for unconsolidated undrained (UU) and consolidated undrained (CU) conditions. This elasto-plastic constitutive Shanghai model has an advantage of describing the mechanical behaviour of over-consolidated and structured soil under different loading and drainage conditions, by using one set of material parameter. In this paper, the Shanghai model used for both UU and CU conditions has the same initial parameters obtained from laboratory test results. The loading conditions and consolidation stages vary based on construction details. The predicted bearing pressure-settlement responses for UU and CU, approves the field observation. The phenomenon of gaining the bearing capacity due to consolidation is captured and explained by the use of soil–water coupled numerical analysis with a new elasto-plastic model. The stress strain behaviour, stress paths and the decay of the structure of elements at different depths presented in this study, reveal the mechanism for the difference between UU and CU conditions to understand the foundation behaviour. Effect of the initial degree of soil structure on the bearing capacity is also addressed. Overall, this approach provides the integrated solution for the shallow foundation design problems under short and long-term loadings.     

Deterministic combination of numerical and physical coastal wave models

A deterministic combination of numerical and physical models for coastal waves is developed. In the combined model, a Boussinesq model MIKE 21 BW is applied for the numerical wave computations. A piston-type 2D or 3D wavemaker and the associated control system with active wave absorption provides the interface between the numerical and physical models. The link between numerical and physical models is given by an ad hoc unified wave generation theory which is devised in the study. This wave generation theory accounts for linear dispersion and shallow water non-linearity. Local wave phenomena (evanescent modes) near the wavemaker are taken into account. With this approach, the data transfer between the two models is thus on a deterministic level with detailed wave information transmitted along the wavemaker.