Effect of a submerged plate on the near-bed dynamics underincoming waves in deep water conditions

It is well known that wave induced bottom oscillations become more and more negligible when the water depth exceeds half the wavelength of the surface gravity wave. However, it was experimentally demonstrated for regular waves that the bottom pressure oscillations at both first and second wave harmonic frequencies could be significant even for incoming waves propagating in deep water condition in the presence of a submerged plate [16]. For a water depth h of about the wavelength of the wave, measurements under the plate (depth immersion of top of plate h/6, length h/2) have shown bottom pressure variations at the wave frequency, up to thirty times larger than the pressure expected in the absence of the plate. In this paper, not only regular but also irregular wave are studied together with wave following current conditions. This behavior is numerically verified by use of a classical linear theory of waves. The wave bottom effect is explained through the role of evanescent modes and horizontally oscillating water column under the plate which still exist whatever the water depth. Such a model, which allows the calculation of the velocity fields, has shown that not only the bottom pressure but also the near bed fluid velocity are enhanced. Two maxima are observed on both sides of the location of the plate, at a distance of the plate increasing with the water depth. The possible impact of such near bed dynamics is then discussed for field conditions thanks to a scaling based on a Froude similarity. It is demonstrated that these structures may have a significant impact at the sea bed even in very deep water conditions, possibly enhanced in the presence of current.     

Computer vision applied to wave flume measurements

A new measurement system for laboratory wave flumes is presented, based on the analysis of digital images by means of computer vision techniques. Unlike conventional wave gauges, the system detects the motions of the free surface along the flume section of interest as opposed to a point. A further advantage lies in that no sensors are necessary within the flume. The system is particularly useful in the case of wave flume tests of floating bodies, such as an oil boom—the motions of the model boom are measured without any sensors in contact with it, which otherwise might place constraints on its displacements. The new method is validated based on free surface measurements carried out with a conventional wave gauge, with excellent results.     

Investigations on 3D effects and correlation between wave height and lip submergence of an offshore stationary OWC wave energy converter

Understanding the hydrodynamic interactions between ocean waves and the oscillating water column (OWC) wave energy converter is crucial for improving the device performance. Most previous relevant studies have focused on testing onshore and offshore OWCs using 2D models and wave flumes. Conversely, this paper provides experimental results for a 3D offshore stationary OWC device subjected to regular waves of different heights and periods under a constant power take–off (PTO) damping simulated by an orifice plate of fixed diameter. In addition, a 3D computational fluid dynamics (CFD) model based on the RANS equations and volume of fluid (VOF) surface capturing scheme was developed and validated against the experimental data. Following the validation stage, an extensive campaign of computational tests was performed to (1) discover the impact of testing such an offshore OWC in a 2D domain or a wave flume on device efficiency and (2) investigate the correlation between the incoming wave height and the OWC front wall draught for a maximum efficiency via testing several front lip draughts for two different rear lip draughts under two wave heights and a constant PTO damping. It is found that the 2D and wave flume modelling of an offshore OWC significantly overestimate the overall power extraction efficiency, especially for wave frequencies higher than the chamber resonant frequency. Furthermore, a front lip submergence equal to the wave amplitude affords maximum efficiency whilst preventing air leakage, hence it is recommended that the front lip draught is minimized.     

Laboratory wave reflection analysis using co-located gages

Most coastal and ocean engineering laboratories employ techniques that use two or three spatially separated wave gages to estimate reflection of irregular waves in two-dimensional wave flumes. This paper presents a frequency domain method for separating incident and reflected wave spectra from co-located gages (gages located on the same vertical line). The technique is based on linear wave theory, and it can be applied to time series of sea surface elevation and horizontal water velocity collected in a vertical array, or it can be used with horizontal and vertical water velocity time series collected at the same point in the water column. Application of the method is limited to those frequencies showing good coherence between time series signals. Outside the range of good coherence, gross inaccuracies occur.The utility of the co-located gage method is illustrated using water velocity data collected in a wave flume with a laser Doppler velocimeter, and the method is validated for the case of complete reflection by a vertical wall. Side-by-side comparison to the spatially-separated wave gage method of Goda and Suzuki (1976) exhibited close agreement for a variety of irregular wave trains being weakly reflected by a mild sloping beach. The co-located gage method is useful in situations where there are spatial variations in the wavelength, such as on a mildly sloping bottom, or in the region close to highly reflective structures where errors arising from spatial variations in characteristic wave parameters would corrupt estimates made using the spatially-separated wave gage method.     

波浪能振荡浮子转换方式实验装置的研究与设计

针对阵列式振荡浮子波浪能转换设备,提出一种在实验室进行测试研究的实验装置。该实验装置由一组阵列柱形水槽组成,水槽液面在水动力系统的驱动下,按照计算机设定的函数做升降运动;对该装置模型——柱形阵列水槽的假设,依据波浪理论进行了论证;对装置的工作原理、结构设计、技术参数、控制电路和操作运行等做了较详细的说明;对实验装置实际运行和测试情况做了介绍。结果表明:该实验装置设计合理可行,模拟波浪高度,远高于现有的造波水池,构成了一种专用的模拟造波装置;该装置为不同形式振荡浮子式波浪能转换设备的研制提供了一种多功能、大振幅、高效灵活的实验平台。     

Contribution of hydrodynamic conditions during shallow water stages to the sediment balance on a tidal flat: Mont-Saint-Michel Bay, Normandy, France

Field measurements were conducted in Mont-Saint-Michel Bay, a megatidal embayment (spring tidal range of 15 m), in order to monitor, over the course of a tidal cycle, sediment transport variability due to waves and tides on the upper part of a tidal flat characterised by shallow water depths. Sensors used to measure currents, water depth and turbidity were installed just above the bed (0.04 m). Two experiments were conducted under contrasting hydrodynamic conditions. The results highlight wave activity over the tidal flat even though observed wind waves were largely dissipated due to the very shallow water depths. Very high suspended sediment concentrations (up to 6 kg/m³) were recorded in the presence of wave activity at the beginning of the local flood, when significant sediment transport occurred, up to 7 times as much as under conditions of no wave activity. This influence may be attributed to the direct action of waves on bed sediments, to wave-induced liquefaction, and to the erosive action of waves on tidal channel banks. The sediment composition, comprising a clay fraction of 2-5%, may also enhance sediment transport by reducing critical shear stress through the sand lubrication effect. The results also show that antecedent meteorological conditions play an important role in suspended sediment transport on the tidal flat. Total sediment flux directions show a net transport towards the inner part of the bay that contributes to deposition over the adjacent salt marshes, and this tendency also prevails during strong wave conditions. Such sediment transport is characterised by significant variability over the course of the tidal cycle. During fair and moderate weather conditions, 83% and 71% of the total flux was observed, respectively, over only 11% and 28% of the duration of the local tidal cycle and with water depths between 0.04 and 0.3 m. These results suggest that in order to improve our understanding of sediment budgets in this type of coastal environment, it is essential to record data just at the beginning and at the end of tidal submergence close to the bed.     

Magnitude and cross-shore distribution of bed return flow measured on natural beaches

Field measurements of cross-shore currents 0.25 m from the bed were made on two natural beaches under a range of incident wave conditions. The results indicated the presence of a relatively strong, offshore-directed mean current, both within and seaward of the surf zone. Typical velocities within the surf zone were of the order of 0.2–0.3 m/s. This bed return flow, or “undertow”, represents a mass conservation response, returning water seaward that was initially transported onshore in the upper water column, primarily above the trough of the incident waves. The measurements demonstrated that the bed return flow velocity increases with the incident wave height. In addition, the crossshore distribution of the bed return flow is characterised by a mid-surf zone maximum, which exhibits a strong decrease in velocity towards the shoreline and a more gradual decay in the offshore direction. Several bed return flow models based on mass continuity were formulated to predict the cross-shore distribution of the bed return flow under an irregular wave field and were compared with the field data. Best agreement was obtained using shallow water linear wave theory, after including the mass transport associated with unbroken waves. The contribution of the unbroken waves enables net offshore-directed bottom currents to persist outside the region of breaking waves, providing a mechanism, other than rip currents, to transport sediment offshore beyond the surf zone.     

Detection of free gas and gas hydrate based on 3D seismic data and cone penetration testing: An example from the Nigerian Continental Slope

We present a new method to characterize free gas, gas hydrates and carbonate concretions occurrence which are considered as high-risk factors for sub-sea developments in the Niger delta. This method is based on the combination of 3D seismic data to the geotechnical site characterizations using piezocone CPTU tests (Cone Penetration Test with additional measurement of the pore water pressure). A special processing of the 3D seismic data has enabled the determination of the interval compressional velocity. Using the effective-medium theory, velocity anomalies (negative and positive) within the first 15 m were translated in gas hydrate and free gas distribution. The calibration of the P wave velocity anomalies was done thanks to in-situ geotechnical testing carried out during two oceanographic surveys (2003 and 2004). Comparison between in-situ testing, recovered cores and the prediction of the gas and the gas hydrate distribution based on the compressional wave velocity have shown that 3D seismic data is a valuable tool to identify heterogeneous areas but the use of the piezocone was essential to discriminate between gas hydrate occurrences and carbonate concretions' presence. Furthermore, in-situ compressional wave velocity (V_p) measurements have clearly demonstrated what it was suspected from the 3D seismic data, the co-existence in the study area between gas hydrate and free gas.     

Wave-formed structures and paleoenvironmental reconstruction

Wave-formed sedimentary structures can be powerful interpretive tools because they reflect not only the velocity and direction of the oscillatory currents, but also the length of the horizontal component of orbital motion and the presence of velocity asymmetry within the flow. Several of these aspects can be related through standard wave theories to combinations of wave dimensions and water depth that have definable natural limits. For a particular grain size, threshold of particle movement and that of conversion from a rippled to flat bed indicate flow-velocity limits. The ratio of ripple spacing to grain size provides an estimate of the length of the near-bottom orbital motion. The degree of velocity asymmetry is related to the asymmetry of the bedforms, though it presently cannot be estimated with confidence. A plot of water depth versus wave height (h—H diagram) provides a convenient approach for showing the combination of wave parameters and water depths capable of generating any particular structure in sand of a given grain size. Natural limits on wave height and inferences or assumptions regarding either water depth or wave period based on geologic evidence allow refinement of the paleoenvironmental reconstruction. The assumptions and the degree of approximation involved in the different techniques impose significant constraints. Inferences based on wave-formed structures are most reliable when they are drawn in the context of other evidence such as the association of sedimentary features or progradational sequences.     

Surface gravity wave interaction with elastic bottom

In the present paper, a hydroelastic model is developed to deal with surface gravity wave interaction with an elastic bed based on the small amplitude water wave theory and plate deflection in finite water depth. The elastic bottom bed is modelled as a thin elastic plate and is based on the Euler-Bernoulli beam equation. The wave characteristics in the presence of the elastic bed is analyzed in both the cases of deep and shallow water waves. Further, the linearized long wave equation is generalized to include bottom flexibility. A generalized expansion formula for the velocity potential is derived to deal with the boundary value problems associated with surface gravity waves having an elastic bed. The utility of the expansion formula is illustrated by demonstrating specific physical problems which will play significant role in the analysis of wave structure interaction problems. Behavior of the wave spectra are discussed in the case of closed basin having a free surface and an elastic bottom topography.     

Predictability of wave-induced net sediment transport using the conventional 1DV RANS diffusion model

Based on a large database of laboratory experiments, the predictability of the conventional one-dimensional vertical Reynolds-averaged Navier–Stokes (RANS) diffusion model is systematically investigated with respect to wave-induced net sediment transport. The predicted net sediment transport rates are compared with the measured data of 176 physical experiments in wave flumes and oscillating water tunnels, covering a wide range of wave conditions (surface, skewed, and asymmetric waves with and without currents), sediment conditions (fine, medium, and coarse sands with median grain diameters ranging from 0.13 to 0.97 mm) and bed forms (flat beds and rippled beds), corresponding to various sediment dynamic regions in the near-shore area. Comparisons show that the majority (73 %) of predictions on a flat bed are within a factor 2 of the measurements. The model behaves much better for medium/coarse sand than for fine sand. The model generally underpredicts the transport rates beneath asymmetric waves and overpredicts the fine sand transport beneath skewed waves. Nevertheless, the model behaves well in reproducing the transport rates under surface waves. A detailed discussion and a quantitative measure of the overall model performance are made. The poor model predictability for fine sand cases is mainly due to the underestimation of unsteady phase-lag effect. It is revealed that the model predictability can be significantly improved by implementing alternative bedload formulas and incorporating more physical processes (mobile-bed roughness, hindered settling, and turbulence damping).     

An experimental study of rip channel flow

A laboratory study of the flow over a bar with a single rip channel has been performed. First, the well-known pattern of a bar circulation cell with a strong offshore-directed current out through the rip channel and a weaker onshore-directed return flow over the bar is documented. Then measurements of the three-dimensional structure of the flow in the area where the rip channel, the bar and the trough meet and well inside the rip channel are presented. These measurements reveal that 3D effects play an important role, and that a depth-integrated viewpoint may not always be sufficient for predicting the flow in the near bed region. Particle-tracking experiments illustrate the near bed flow pattern over the entire area. These demonstrate how the overall trajectory pattern changes as a function of the distance of wave breaking from the bar crest: For some conditions, the rip current is fed from the trough and for other conditions it is fed directly from the bar. Both the 3D measurements and the trajectory tests show the existence of a weaker onshore-directed near-bed drift in the area where the rip current ceases. Finally, in a series of sensitivity tests, measurements of the rip current intensity for different wave climate and water level conditions reveal a strong correlation between the rip current intensity and the wave height (both normalized).     

对应分析方法在粤西水东湾动力沉积环境分析的应用

本文以粤西水东湾表层沉积物样品的粒度分布为基本数据,使用对应分析方法并结合地形与波流特征,分析了海湾动力沉积环境。分析结果表明,海湾现代沉积环境可划分为4个动力沉积区。     

风浪和涌浪分离方法的比较

海浪通常以风浪和涌浪混合的形式存在。文中利用模型试验和实测资料,对目前提出的一种二维谱风涌浪分离方法(2D法)和3种一维谱风涌浪分离方法(PM法、WH法、JP法)进行了检验,分析发现:2D法给出的结果整体而言最为可靠,与2D法相比,PM法明显高估了风浪成分,WH法低风速时高估了风浪,高风速时跟2D法比较接近,而JP法在整体上高估了风浪成分。通过调整分割频率的比例系数,改进了PM法,改进后的PM法给出的分离结果与2D法最为一致。     

Implementation and modification of a three-dimensional radiation stress formulation for surf zone and rip-current applications

Regional Ocean Modeling System (ROMS v 3.0), a three-dimensional numerical ocean model, was previously enhanced for shallow water applications by including wave-induced radiation stress forcing provided through coupling to wave propagation models (SWAN, REF/DIF). This enhancement made it suitable for surf zone applications as demonstrated using examples of obliquely incident waves on a planar beach and rip current formation in longshore bar trough morphology (Haas and Warner, 2009). In this contribution, we present an update to the coupled model which implements a wave roller model and also a modified method of the radiation stress term based on Mellor (2008, 2011a,b,in press) that includes a vertical distribution which better simulates non-conservative (i.e., wave breaking) processes and appears to be more appropriate for sigma coordinates in very shallow waters where wave breaking conditions dominate. The improvements of the modified model are shown through simulations of several cases that include: (a) obliquely incident spectral waves on a planar beach; (b) obliquely incident spectral waves on a natural barred beach (DUCK'94 experiment); (c) alongshore variable offshore wave forcing on a planar beach; (d) alongshore varying bathymetry with constant offshore wave forcing; and (e) nearshore barred morphology with rip-channels. Quantitative and qualitative comparisons to previous analytical, numerical, laboratory studies and field measurements show that the modified model replicates surf zone recirculation patterns (onshore drift at the surface and undertow at the bottom) more accurately than previous formulations based on radiation stress (Haas and Warner, 2009). The results of the model and test cases are further explored for identifying the forces operating in rip current development and the potential implication for sediment transport and rip channel development. Also, model analysis showed that rip current strength is higher when waves approach at angles of 5° to 10° in comparison to normally incident waves.     

Experimental determination of wavelengths in the presence of a mean current

This paper describes a simple method for determining the wavelength of small amplitude waves under laboratory conditions where reflected wave components are present both with and without a mean current flow superimposed. It assumes a locally horizontal bed but requires no a priori assumption concerning the form of the dispersion relation with a coexisting current. Synchronous measurements of the water surface recorded along any straight line are analysed to yield Fourier coefficients at each location. It is then shown that for all practical conditions excluding a perfect standing wave, the average rate of change of wave phase in the chosen direction can be related directly to the component of incident wave number in that direction, irrespective of reflection coefficient or relative current strength. The technique has been applied to regular and bichromatic waves in a flume with an absorbing wave generator, and can also be applied in 3-D wave basins where waves and currents intersect at arbitrary angles. In combined wave–current experiments, by assuming the linear dispersion relation, it is also possible to estimate the effective current velocity.     

A numerical scheme for morphological bed level calculations

The typical equation for bed level change in sediment transport in river, estuary and near shore systems is based on conservation of sediment mass. It is generally a nonlinear conservation equation for bed level. The physics here are similar to shallow water wave equations and gas dynamics equation which will develop shock waves in many circumstances. Many state-of-art morphological models use classical lower order Lax–Wendroff or modified Lax–Wendroff schemes for morphology which are not very stable for long time sediment transport processes simulation. Filtering or artificial diffusion are often added to achieve stability. In this paper, several shock capturing schemes are discussed for simulating bed level change with different accuracy and stability behaviors. The conclusion is in favor of a fifth order Euler-WENO scheme which is introduced to sediment transport simulations here over other schemes. The Euler-WENO scheme is shown to have significant advantages over schemes with artificial viscosity and filtering processes, hence is highly recommended especially for phase-resolving sediment transport models.     

Approximate Stream Function wavemaker theory for highly non-linear waves in wave flumes

An approximate Stream Function wavemaker theory for highly non-linear regular waves in flumes is presented. This theory is based on an ad hoc unified wave-generation method that combines linear fully dispersive wavemaker theory and wave generation for non-linear shallow water waves. This is done by applying a dispersion correction to the paddle position obtained for non-linear long waves. The method is validated by a number of wave flume experiments while comparing with results of linear wavemaker theory, second-order wavemaker theory and Cnoidal wavemaker theory within its range of application.     

Breaking solitary wave evolution over a porous underwater step

Solitary wave evolution over a shelf including porous damping is investigated using Volume-Averaged Reynolds Averaged Navier–Stokes equations. Porous media induced damping is determined based on empirical formulations for relevant parameters, and numerical results are compared with experimental information available in the literature. The aim of this work is to investigate the effect of wave damping on soliton disintegration and evolution along the step for both breaking and non-breaking solitary waves. The influence of several parameters such as geometrical configuration (step height and still water level), porous media properties (porosity and nominal diameter) or solitary wave characteristics (wave height) is analyzed. Numerical simulations show the porous bed induced wave damping is able to modify wave evolution along the step. Step height is observed as a relevant parameter to influence wave evolution. Depth ratio upstream and downstream of the edge appears to be the more relevant parameter in the transmission and reflection coefficients than porosity or the ratio of wave height–water depth. Porous step also modifies the fission and the solitary wave disintegration process although the number of solitons is observed to be the same in both porous and impermeable steps. In the absence of breaking, porous bed triggers a faster fission of the incident wave into a second and a third soliton, and the leading and the second soliton reduces their amplitude while propagating. This decrement is observed to increase with porosity. Moreover, the second soliton is released before on an impermeable step. Breaking process is observed to dominate over the wave dissipation at the porous bottom. Fission is first produced on a porous bed revealing a clear influence of the bottom characteristics on the soliton generation. The amplitude of the second and third solitons is very similar in both impermeable and porous steps but they evolved differently due to the effect of bed damping.     

On the accurate and efficient calibration of a 3D wave basin

This paper describes the calibration procedure adopted for the new 3D wave basin located in the Hydrodynamics Laboratory at Imperial College London. Unlike traditional calibrations, based on observations of regular wave trains, the method described herein uses a focused wave approach. Such waves, produced by the constructive interference of freely propagating wave components, have led to a number of recent advances in theoretical wave modelling in which it was essential to know the underlying linear components. In the context of a laboratory study, similar advantages can be realised provided the linear wave components generated by the wave paddles are well defined. This, in turn, can only be achieved if the calibration is sufficiently accurate. The current study provides a calibration based upon a realistic JONSWAP spectrum, describes the details of the methodology employed, and highlights how the application of focused wave techniques eliminates spurious calibration effects due to unwanted reflections from the boundaries of the basin. The final calibration is verified through the generation of test cases, involving linear and nonlinear, unidirectional and directionally spread waves. These confirm both the accuracy of the calibration and the suitability of the methods employed.