复杂地形下海浪数值模式的特征线计算格式

LAGFD-WAM海浪数值模式是一种第三代海浪数值模式，通过求解波数谱平衡方程，并考虑风输入、波浪破碎耗散、底摩擦耗散、波波非线性相互作用和波流相互作用等源函数，模拟波数空间下的海浪方向谱，并依此获得海浪的波高、周期和平均波向。该模式的一个显著特点是采用特征线嵌入格式求解海浪的传播。在进行浅水区域的海浪模拟时，特征线嵌入格式的数值计算方案是否合理对海浪数值模拟结果产生直接的影响。为此LAGFD-WAM海浪数值模式提出了一种新的特征线混合数值计算格式，并应用于浅水海浪数值模拟。结果表明，采用该计算方法，能够使数值模拟结果与实测结果很好符合。     

LAGFD-WAM numerical wave model——Ⅰ. Basic physical model

The LAGFD-WAM wave model is a third generation wave model. In the present paper the physical aspect of the model was shown in great detail including energy spectrum balance equation, complicated characteristics equations and source functions.     

球坐标系下MASNUM海浪数值模式的建立及其应用

为开展海浪对海洋上混合层的搅拌混合作用及其对海气界面通量的影响等研究,在LAGFD WAM区域海浪数值模式基础上建立了球坐标系下的全球海浪数值模式.重点导出了球坐标系下的海浪能量谱平衡方程及其复杂特征线方程,该组方程包含了背景流场对波动传播的调整、波动沿大圆传播的折射等.数值积分则采用复杂特征线嵌入计算格式.初步数值模拟结果表明,该海浪全球数值模式能够较为精确地刻画海浪的动力过程.     

利用海浪谱共轭方程进行东海敏感性区域分布研究

近岸局部范围内的波浪是由远处洋面上生成并传播过来的，它的生成发展源地(敏感性海域)对于所关注范围的波浪起决定性影响。共轭方程中的模式变量反映了物理参量的梯度变化，利用它可以达到精确确定敏感海域具体位置的目的。本文在LAGFD-WAM海浪模式基础上建立了海浪谱共轭方程，提出了通过在关注点邻域加一小扰动作为逆向积分海浪谱共轭方程的驱动项，来计算出梯度极值分布的思想。将这一思想应用于东中国海区域，通过一系列试验来确定舟山及长江口附近海域的敏感性区域位置。试验结果符合物理意义，梯度极值分布结果为下一步的多源卫星观测资料复合分析及最优实测方案设计奠定了基础。     

基于MPI的LAGFD-WAM海浪数值模式并行算法研究

在LAGFD-WAM海浪串行数值模式基础上,利用MPI信息传递机制实现其并行化。通过对模拟区域合理划分,对数据采取分块加载,实现了各个节点的负载平衡;通过对算法的改进实现了粗粒度计算,大幅度减少了通信量,从而提高了程序的执行效率。对串行计算和并行计算的效率比较表明,本文建立的方法能够得到较高的加速比。对全球海浪模式,加速比和CPU数目大致呈线性关系。对高分辨率的区域海浪模式,在128 CPU条件下加速比可以达到91.9。     

卫星高度计波高数据同化对西北太平洋海浪数值预报的影响评估

通过数值试验验证卫星高度计波高数据同化对西北太平洋3 d海浪预报的改进效果。驱动海浪模式的强迫场采用国家海洋环境预报中心基于MM5模式预报的风场,波高数据同化使用的观测数据是Jason-1卫星高度计有效波高。用最优插值数据同化方法获得海浪有效波高的最优估计并重构相应的海浪方向谱,以此为初始场进行为期3 d的数值预报试验。与没有同化的预报进行了比较和分析,结果表明卫星高度计海浪数据同化对0~72 h预报有不同程度的明显改善,改进程度随预报时效的增加而减少。     

Experimental Research on Wave Transmission over Submerged Rubble Mound Breakwaters

This paper discusses some previous, and presents some new experimental results on wave transmission over submerged breakwaters. The objective of this study is to evaluate wave transmission coefficient and develop a two-dimensional （2D） model as an improvement to the existing wave transmission coefficient models. Factors which affect wave transmission over stbmerged breakwaters are discussed through a series of laboratory experiments. Basic recommendations for evaluation and design of submerged rubble-monud breakwaters are presented. From the test results, a calculation formula of wave transmission coefficient is proposed.     

Assimilation of altimeter wave height data into a spectral wave model using statistical interpolation

A scheme for the assimilation of altimeter wave height data into a second generation wave prediction model is developed. The scheme modifies the model wave spectrum so as to be consistent with the satellite observed values of significant wave height. This modification is achieved, however, so that the adjusted wave field is consistent with the model physics and the forcing wind field. In this manner the modifications to the wave field persist in the model, thus yielding long term improvements in model performance. In addition, a statistical interpolation scheme is used to ensure that maximum use is made of the point observations made by the satellite. In this manner, not only points directly beneath the satellite track are updated. Points adjacent to the track are also modified, the extent of this modification depending on the spatial correlation of the wave field. The scheme is applied to a computationally efficient second generation wave model. The improvement in performance is significant, demonstrating that assimilation may be a very efficient alternative to proceeding to more sophisticated and expensive third generation models. This is particularly true where the forcing wind field may be of poor quality. The results also demonstrate that with the addition of assimilation, relatively small computational grids can be utilized. Swell generated external to the grid will be included through the assimilation cycle.     

Iceberg drift affected by wave action

A mathematical model to predict iceberg drift pattern has been developed, which includes the wave drift force, in addition to the other conventional force components such as forces due to wind, current, Coriolis effect, and geostrophic effect. Trajectories of two icebergs were computed first with the wave effect then without the wave effect. All were compared with the observed results from the field. The model with the wave effect shows a significant improvement in the correlation.     

Random breaking waves: Field verification of a wave-by-wave algorithm for engineering application

A simple numerical algorithm is utilized to test the wave-by-wave approach to modeling random wave transformation caused by shoaling, breaking and reforming across nearshore profiles of arbitrary shape. Histograms of wave height and associated statistics predicted by the model show favorable comparison to measurements from the nine tests of the DUCK'85 field experiment, with a root mean square error less than 10% for most wave statistics. Although suitable for use in many engineering applications, the model does overpredict the number of waves that reform as they pass over a trough formation, which is one facet that warrants improvement. Using the model predictions as a point of reference, the qualitative effects of (1) an offshore wind, (2) wave reflection and (3) surf beat are identified and discussed.     

The synchronous coupling of a third-generation wave model and a two-dimensional storm surge model

The interaction of waves and currents is studied by the dynamical coupling of a third generation wave model and a two-dimensional storm surge model. The coupling process of the two models is implemented synchronously. To estimate the effects of waves on the generation of storm surges, the theory of Janssen is used. The effects of the wave radiation stress on surge levels and the effects of storm-induced currents on waves are also investigated.The coupled wave and storm surge models have been tested by hindcasting two storm events in the northern South China Sea. The use of the Simth and Banke stress relation underestimates the surges by 10%. The inclusion of the radiation stress improves the accuracy of the computed results slightly by 2%. The introduction of a wave-dependent surface drag gives a significant improvement. The storm-induced currents clearly affect the wave characteristics at the peak stage. However, as far as the prediction of wave height is concerned, it is better not to consider the wave radiation stress in the storm surge model unless this is accompanied by a wave-dependent surface drag.     

On Radiation Boundary Conditions and Wave Transformation Across the Surf Zone

The purpose of this paper is to extend the validity of Li's parabolic model (1994) by incorporating a combined energy factor in the mild-slope equation and by improving the traditional radiation boundary conditions. With wave breaking and energy dissipation expressed in a direct form in the equation, the proposed model could provide an efficient numerical scheme and accurate predictions of wave transformation across the surf zone. The radiation boundary conditions are iterated in the model without use of approximations. The numerical predictions for wave height distributions across the surf zone are compared with experimental data over typical beach profiles. In addition, tests of waves scattering around a circular pile show that the proposed model could also provide reasonable improvement on the radiation boundary conditions for large incident angles of waves.     

Numerical Study of Wave Diffraction Effect Introduced in the SWAN Model

New version of SWAN model includes the wave diffraction effect which is the main improvement compared with the previous versions.Experimental data collected in the wave basin of the University of Delaware were used to test its performance.Wave heights were compared in the four cases(with different wave energies and directional spreading spectra).The results agreed well with the measurements,especially for the broad directional spectra cases.The effect of wave diffraction was analyzed by switching on/off the corresponding term.By introducing the diffraction term,the distributions of wave height and wave direction were smoothed,especially obvious for the narrow spectrum cases.Compared with the calculations without diffraction,the model with diffraction effect gave better results.     

Experimental Research on Wave Transmission over Submerged Rubble-Mound Breakwaters

This paper discusses some previous, and presents some new experimental results on wave transmission oversubmerged breakwaters. The objective of this study is to evaluate wave transmission coefficient and develop two-dimensional (2D) model as an improvement to the existing wave transmission coefficient models. Factors which affect wave transmission over submerged breakwaters are discussed through series of laboratory experiments. Basic recommendations for evaluation and design of submerged rubble-monud breakwaters are presented.From the test results, calculation formula of wave transmission coefficient is proposed.     

CFD modelling of a small–scale fixed multi–chamber OWC device

Wave Energy Converters (WECs) have excellent potential as a source of renewable energy that is yet to be commercially realised. Recent attention has focused on the installation of Oscillating Water Column (OWC) devices as a part of harbor walls to provide advantages of cost–sharing structures and proximity of power generation facilities to existing infrastructure. In this paper, an incompressible three–dimensional CFD model is constructed to simulate a fixed Multi–Chamber OWC (MC–OWC) device. The CFD model is validated; the simulation results are found to be in good agreement with experimental results obtained from a scale physical model tested in a wave tank. The validated CFD model is then used for a benchmark study of 96 numerical tests. These investigate the effects of the PTO damping caused by the power take–off (PTO) system on device performance. The performance is assessed for a range of regular wave heights and periods. The results demonstrate that a PTO system with an intermediate damping can be used for all chambers in the MC–OWC device for most wave period ranges, except for the long wave periods. These require a higher PTO damping. An increased incident wave height reduces the device capture width ratio, but there is a noticeable improvement for long wave periods.     

SWAN模型对实验室小尺度浅水波浪模拟研究

基于浅水斜坡地形的物理模型试验数据,考察SWAN模型对实验室小尺度浅水波浪的模拟效果,进而检验其浅水项的模拟精度。模拟中采用直接输入初始测点的实测海浪谱进行造波,重点考察浅水中三波相互作用和变浅破碎两个源项,对不同工况下,SWAN模式在水深条件变化下的有效波高、谱平均周期、海浪谱演化的模拟能力进行研究。研究表明:模拟的有效波高较符合实测波浪的增长和衰减,但谱平均周期计算值明显偏小;海浪谱的能量转移机制同实测有较大区别,频谱模拟结果出现高频高估、低频低估现象。对两个源项进行对比分析得出三波相互作用对海浪谱的能量转换影响远大于变浅破碎耗散。想要提高近岸区谱平均周期和海浪谱的模拟精度则SWAN模型中三波非线性项的计算精确度仍需更多研究和改进。     

Wave characteristics and extreme parameters in the Bohai Sea

This paper is aimed at the whole Bohai Sea,as the complement and improvement of wave characteristics and extreme parameters.Wave fields were simulated in the Bohai Sea by using wave model SWAN from 1985 to 2004.The input data based on the hindcast of high-resolution wind fields from RAMS and water level fields from POM,which have been tested and verified well.Comparisons of significant wave heights between simulation and station observations show a good agreement in general.By statistical analysis,the wave characteristics such as significant wave heights, dominant wave directions and their seasonal variations are discussed.In addition,main wave extreme parameters and directional extreme values particularly for 100-year return period are investigated.     

Prediction of swash motion and run-up including the effects of swash interaction

Modifications to a model describing swash motion based on solutions to the non-linear shallow water equations were made to account for interaction between up-rush and back-wash at the still water shoreline and within the swash zone. Inputs to the model are wave heights and arrival times at the still water shoreline. The model was tested against wave groups representing idealized vessel-generated wave trains run in a small wave tank experiment. Accounting for swash interaction markedly improved results with respect to the maximum run-up length for cases with rather gentle foreshore slopes (tanβ=0.07). For the case with a steep foreshore slope (tanβ=0.20) there was very little improvement compared to model results if swash interaction was not accounted for. In addition, an equation was developed to predict the onset and degree of swash interaction including the effects of bed friction.     

Simulation of wave overtopping by an incompressible SPH model

The paper presents an incompressible Smoothed Particle Hydrodynamics (SPH) model to investigate the wave overtopping of coastal structures. The SPH method is a grid-less Lagrangian approach which is capable of tracking the large deformations of the free surface with good accuracy. The incompressible algorithm of the model is implemented by enforcing the constant particle density in the pressure projection. The SPH model is employed to reproduce a transient wave overtopping over a fixed horizontal deck and the regular/irregular waves overtopping of a sloping seawall. The computations are validated against the experimental and numerical data and a good agreement is observed. The SPH modelling is shown to provide a promising tool to predict the overtopping characteristics of different waves. The present model is expected to be of practical purpose if further improvement in the spatial resolution and CPU time can be adequately made.     

Uncertainty Analysis of Wind-Wave Predictions in Lake Michigan

With all the improvement in wave and hydrodynamics numerical models, the question rises in our mind that how the accuracy of the forcing functions and their input can affect the results. In this paper, a commonly used numerical third-generation wave model, SWAN is applied to predict waves in Lake Michigan. Wind data are analyzed to determine wind variation frequency over Lake Michigan. Wave predictions uncertainty due to wind local effects are compared during a period where wind has a fairly constant speed and direction over the northern and southern basins. The study shows that despite model calibration in Lake Michigan area, the model deficiency arises from ignoring wind effects in small scales. Wave prediction also emphasizes that small scale turbulence in meteorological forces can increase prediction errors by 38%. Wave frequency and coherence analysis show that both models can predict the wave variation time scale with the same accuracy. Insufficient number of meteorological stations can result in neglecting local wind effects and discrepancies in current predictions. The uncertainty of wave numerical models due to input uncertainties and model principals should be taken into account for design risk factors.