On the modeling of wave propagation on non-uniform currents and depth

By transforming two different time-dependent hyperbolic mild slope equations with dissipation term for wave propagation on non-uniform currents into wave-action conservation equation and eikonal equation, respectively, shown are the different effects of dissipation term on the eikonal equation in the two different mild slope equations. The performances of intrinsic frequency and wave number are also discussed. Thus the suitable mathematical model is chosen in which the wave number vector and intrinsic frequency are expressed both more rigorously and completely. By using the perturbation method, an extended evolution equation, which is of time-dependent parabolic type, is developed from the time-dependent hyperbolic mild slope equation which exists in the suitable mathematical model, and solved by using the alternating direction implicit (ADI) method. Presented is the numerical model for wave propagation and transformation on non-uniform currents in water of slowly varying topography. From the comparisons of the numerical solutions with the theoretical solutions of two examples of wave propagation, respectively, the results show that the numerical solutions are in good agreement with the exact ones. Calculating the interactions between incident wave and current on a sloping beach [Arthur, R.S., 1950. Refraction of shallow water waves. The combined effects of currents and underwater topography. EOS Transactions, August 31, 549–552], the differences of wave number vector between refraction and combined refraction–diffraction of waves are discussed quantitatively, while the effects of different methods of calculating wave number vector on numerical results are shown.     

A Eulerian–Lagrangian method for the simulation of wave propagation

A Eulerian–Lagrangian method (ELM) is employed for the simulation of wave propagation in the present research. The wave action conservation equation, instead of the wave energy balance equation, is used. The wave action is conservative and the action flux remains constant along the wave rays. The ELM correctly accounts for this physical characteristic of wave propagation and integrates the wave action spectrum along the wave rays. Thus, the total derivative for wave action spectrum may be introduced into the numerical scheme and the complicated partial differential wave action balance equation is simplified into an ordinary differential equation. A number of test cases on wave propagation are carried out and show that the present method is stable, accurate and efficient. The results are compared with analytical solutions and/or other computed results. It is shown that the ELM is superior to the first-order upwind method in accuracy, stability and efficiency and may better reflect the complicated dynamics due to the complicated bathymetry features in shallow water areas.     

单点高频地波雷达资料估算潮能耗散的方法（英文）

使用嵊泗站所布设的地波雷达观测获取的径向流数据,以及嵊泗、芦潮港、岱山3个潮汐观测站水位资料,采用两点近似投影方法反演流场全矢量流速,并用T_tide程序计算调和常数,分别计算O1、K1、M2、S2各分潮流速场及迟角场,并计算各点上的潮能通量及潮能耗散,得到嵊泗岛以西杭州湾口区域潮能耗散同地形存在良好对应关系,充分证明了采用地波雷达观测数据进行潮能耗散计算这一方法的可行性,供相关工作者作进一步研究和讨论。     

单点高频地波雷达资料估算潮能耗散的方法

使用嵊泗站所布设的地波雷达观测获取的径向流数据,以及嵊泗、芦潮港、岱山3个潮汐观测站水位资料,采用两点近似投影方法反演流场全矢量流速,并用T-Tide程序计算调和常数,分别计算O1、K1、M2、S2各分潮流速场及迟角场,并计算各点上的潮能通量及潮能耗散,得到嵊泗岛以西杭州湾口区域潮能耗散同地形存在良好对应关系,充分证明了采用地波雷达观测数据进行潮能耗散计算这一方法的可行性,供相关工作者作进一步研究和讨论.     

Near-optimal control of a wave energy device in irregular waves with deterministic-model driven incident wave prediction

This paper investigates wave-by-wave control of a wave energy converter using incident wave prediction based on up-wave surface elevation measurement. The goal of control is to approach the hydrodynamically optimum velocity leading to optimum power absorption. This work aims to study the gains in energy conversion from a deterministic wave propagation model that accounts for a range of group velocities in deriving the prediction. The up-wave measurement distance is assumed to be small enough to allow a deterministic propagation model, and further, both wave propagation and device response are assumed to be linear. For deep water conditions and long-crested waves, the propagation process is also described using an impulse response function (e.g. [1]). Approximate low and high frequency limits for realistic band-limited spectra are used to compute the corresponding group velocity limits. The prediction time into the future is based on the device impulse response function needed for the evaluation of the control force. The up-wave distance and the duration of measurement are then determined using the group velocity limits above.A 2-body axisymmetric heaving device is considered, for which power capture is through the relative heave oscillation between the two co-axial bodies. The power take-off is assumed to be linear and ideal as well as capable of applying the necessary resistive and reactive load components on the relative heave oscillation. The predicted wave profile is used along with device impulse response functions to compute the actuator force components at each instant. Calculations are carried out in irregular waves generated using a number of uni-modal wave spectra over a range of energy periods and significant wave heights. Results are compared with previous studies based on the use of instantaneous up-wave wave-profile measurements, both without and with oscillation constraints imposed. Considerable improvements in power capture are observed with the present approach over the range of wave conditions studied.     

An Improved Coupling of Numerical and Physical Models for Simulating Wave Propagation简

An improved coupling of numerical and physical models for simulating 2D wave propagation is developed in this paper. In the proposed model, an unstructured finite element model （FEM） based Boussinesq equations is applied for the numerical wave simulation, and a 2D piston-type wavemaker is used for the physical wave generation. An innovative scheme combining fourth-order Lagrange interpolation and Runge-Kutta scheme is described for solving the coupling equation. A Transfer function modulation method is presented to minimize the errors induced from the hydrodynamic invalidity of the coupling model and/or the mechanical capability of the wavemaker in area where nonlinearities or dispersion predominate. The overall performance and applicability of the coupling model has been experimentally validated by accounting for both regular and irregular waves and varying bathymetry. Experimental results show that the proposed numerical scheme and transfer function modulation method are efficient for the data transfer from the numerical model to the physical model up to a deterministic level.     

Generation and propagation of transient nonlinear waves in a wave flume

A semi-analytical nonlinear wavemaker model is derived to predict the generation and propagation of transient nonlinear waves in a wave flume. The solution is very efficient and is achieved by applying eigenfunction expansions and FFT. The model is applied to study the effect of the wavemaker and its motion on the generation and propagation of nonlinear waves. The results indicate that the linear wavemaker theory may be applied to predict only the generation of waves of low steepness for which the nonlinear terms in the kinematic wavemaker boundary condition and free-surface boundary conditions are of secondary importance. For waves of moderate steepness and steep waves these nonlinear terms have substantial effects on wave profile and wave spectrum just after the wavemaker. A wave spectrum corresponding to a sinusoidally moving wavemaker possesses a multi-peak form with substantial nonlinear components, which disturbs or may even exclude physical modeling in wave flumes. The analysis shows that the widely recognized weakly nonlinear wavemaker theory may only be applied to describe the generation and propagation of waves of low steepness. This is subject to further restrictions in shallow and deep waters because the kinematic wavemaker boundary condition as well as the nonlinear interaction of wave components and the evolution of wave energy spectrum is not properly described by weakly nonlinear wavemaker theory. Laboratory experiments were conducted in a wave flume to verify the nonlinear wavemaker model. The comparisons show a reasonable agreement between predicted and measured free-surface elevation and the corresponding amplitudes of Fourier series. A reasonable agreement between theoretical results and experimental data is observed even for fairly steep waves.     

Effect of sampling rate on wave height statistics

This paper examines the effects of digital processing, and therefore discretisation or sampling, of sea surface elevations that are, in reality, continuous. Using random linear wave theory, probability distributions for the measured (as opposed to actual) wave amplitudes and heights have been obtained which are conditional on the sampling rate. It is shown that with low sampling rates there are significant departures from the usual Rayleigh distribution. Rates of 1 Hz or less may lead to significant underestimation of the probability of very large waves. An analysis of full-scale measurements obtained from a platform in the North Sea supports these results.     

Simulation of wave propagation over a submerged bar using the VOF method with a two-equation k– turbulence modeling

A two-equation k– turbulence model is used in this paper to simulate the propagation of cnoidal waves over a submerged bar, where the free surface is handled by the volume-of-fluid (VOF) method. Using a VOF partial-cell variable and a donor–acceptor method, the model is capable of treating irregular boundaries, including arbitrary bottom topography and internal obstacles, where the no-slip condition is satisfied. The model also allows the viscous sublayer to be modeled by a wall function approximation implemented in the grid nodes that are immediately adjacent to a wall boundary. The numerical model applied to the propagation of cnoidal waves over a submerged bar can produce results that are in general agreement with some laboratory measurements. Some remarks arising from the comparison between the computational and experimental results are presented.     

近岸流对波浪传播影响的数值分析

使用近岸波浪模型SWAN计算存在沿岸流和离岸流时的近岸波浪传播。先设离岸流u=0m/s,模拟均匀、非均匀沿岸流的流速和梯度对波高传播的影响;再设沿岸流v=0.5m/s,模拟均匀、非均匀离岸流的流速和梯度对波能高传播的影响。从模拟中得到,近岸波浪传播受沿岸流、离岸流的流速和梯度影响时,波高的变化规律。     

Measured and modeled wave overtopping on a natural beach

The rate of wave overtopping of a barrier beach is measured and modeled. Unique rate of wave overtopping field data are obtained from the measure of the Carmel River, California, lagoon filling during a time when the lagoon is closed-off with no river inflow. Volume changes are based on measured lagoon height changes applied to a measured hypsometric curve. Wave heights and periods are obtained from directional wave spectra data in 15 m fronting the beach. Beach morphology was measured by GPS walking surveys. Three empirical overtopping models by Van der Meer and Janssen (1995), Hedges and Reis (1998) and Pullen et al. (2007) with differing parameterizations on wave height, period and beach slope and calibrated using extensive laboratory data obtained over plane, impermeable beaches are applied in a quasi-2D manner and compared with the field observations. Three overtopping events are considered when morphology data were available less than 2 weeks prior to the event. The models are tuned to fit the data using a reduction factor to account for beach permeability, berm characteristics, non-normal wave incidence and surface roughness influence. In addition, the run-up model by Stockdon et al. (2006) based on field data is examined and found to underestimate run-up as the calculated values were too small to predict any of the observed overtopping. The three overtopping models performed similarly well with values of 0.72–0.87 for the two narrow-banded wave cases, with an average reduction factor of 0.78. The European model (Pullen et. al., 2007) performed best overall and in particular for the case of the broad-banded, double peaked wave spectrum.     

Numerical study for vegetation effects on coastal wave propagation by using nonlinear Boussinesq model

The vegetation has important impacts on coastal wave propagation. In the paper, the sensitivities of coastal wave attenuation due to vegetation to incident wave height, wave period and water depth, as well as vegetation configurations are numerically studied by using the fully nonlinear Boussinesq model. The model is based on the implementation of drag resistances due to vegetation in the fully nonlinear Boussinesq equation where the drag resistance is provided by the Morison’s formulation for rigid structure induced drag stresses. The model is firstly validated by comparing with the experimental results for wave propagation in vegetation zones. Subsequently, the model is used to simulate waves with different height, period propagating on vegetation zones with different water depth and vegetation configurations. The sensitivities of wave attenuation to incident wave height, wave period, water depth, as well as vegetation configurations are investigated based on the numerical results. The numerical results indicate that wave height attenuation due to vegetation is sensitive to incident wave height, wave period, water depth, as well as vegetation configurations, and attenuation ratio of wave height is increased monotonically with increases of incident wave height and decreases of water depth, while it is complex for wave period. Moreover, more vegetation segments can strengthen the interaction of vegetation and wave in a certain range.     

Effects of wave breaking on wave statistics for deep-water random wave train

The experimental studies of the breaking effects on wave statistics for deep-water random waves are presented. It is especially focused on the behavior of kurtosis of surface elevations due to wave breaking. Wave breaking suppresses the maximum limit of kurtosis of the surface elevation, although skewness depends on characteristic wave steepness. The mean instantaneous wave steepness of breaking waves defined using the zero-down-crossing method was much lower than expected from the Stokes waves.     

海浪特性及其传播变形研究新进展

综述第27届国家海岸工程会议（ICCE2000）上有关海浪特性及其传播变形研究的新进展，内容包括海浪的时域特性，海浪的特征，海浪的传播变形和波浪的破碎，反射等四部分。     

有限长正弦地形上的波浪传播问题研究

应用基于势流理论的频域边界元方法,研究了波浪在有限长正弦周期地形上的传播问题.开展了网格收敛性分析,并与物理模型试验结果做了对比,验证了数值方法的准确性.研究了正弦地形长度和幅度对水波反射率的影响,以及波浪在周期性地形上传播时的沿程波幅变化.研究发现,波浪反射系数随波数的变化及波浪幅值沿程的变化均呈包络状;当频率处于布拉格共振区域边缘且反射率为零时,周期地形上会出现极大的波幅;随海床波幅增大,第一布拉格共振带变宽且向低频移动.     

Wave energy assessment based on trivariate distribution of significant wave height,mean period and direction

Based on the Vine copula theory, a trivariate statistical model of significant wave height, characterized wave period and mean wave direction was constructed. To maintain the properties of the different types of variables, a special copula function was derived from the model developed by Johnson and Wehrly based on the maximum entropy principle. It was then combined with the Archimedean copulas to construct the proposed model. An effective algorithm for generating corresponding joint pseudo-random numbers was also developed. Statistical analysis of hindcast data for the significant wave height, mean wave period, and direction, which were collected from an observation point in the North Atlantic every three hours from 1997 to 2001, was performed. The marginal distributions of the significant wave height and mean wave period were fitted by a modified maximum entropy distribution, and the mean wave direction was fitted by a mixture of von Mises distributions. It was shown that the proposed model is a good fit for the data. The seasonal wave energy resources in the target area were assessed using the model estimates. Histograms of the directional wave energy, wave energy roses, and scatter and energy diagrams were presented.     

Improvement on open boundaries on a time dependent numerical model of wave propagation in the nearshore region

An improvement on the simulation of outgoing waves on a time dependent numerical model for water wave propagation in the nearshore region is presented. The governing equations consist of a system of first order partial differential equations (PDEs), the equation of continuity and the equation of motion. A comparative study of first order radiation boundary conditions (BCs) and first order radiation BCs combined with sponge layers is presented for cases where outgoing waves leave the numerical domain of calculation through the open boundary. A reduction of spurious reflections from the numerical open boundaries can be obtained with an irrelevant increase in terms of computational cost.     

网衣对波浪传播影响的三维数值模拟

为研究网衣对波浪传播的影响,采用多孔介质模型模拟网衣,建立了模拟网衣与波浪相互作用的三维数值波浪水槽模型。基于该数值模型研究了波浪经过网衣作用后波高衰减的变化规律,并与物理模型试验结果进行比较。通过对比数值模拟和物理模型试验的结果,证明了该三维数值模型模拟网衣对波浪传播影响的可行性。     

基于改进缓坡方程的波浪传播数值模拟

用变分原理导出考虑底坡一阶导数平方项和二阶曲率项影响的缓坡方程,对传统缓坡方程作了改进,提高波浪在海底地形变化剧烈、水深较浅时数值模拟精度。数值计算与已有实验室试验资料比较表明,该模型可以较好地模拟有剧烈变化的海底地形的波浪传播,比传统缓坡方程模型计算结果在精度上有明显提高。