礁坪上波浪传播的数值模拟

本文基于具备间断捕捉能力的二阶全非线性Boussinesq数值模型,对规则波和随机波在礁坪地形上的传播变形进行了数值模拟。该模型采用高阶有限体积法和有限差分方法求解守恒格式的控制方程,将波浪破碎视为间断,同时采用静态重构技术处理了海岸动边界问题。重点针对礁坪上波浪传播过程中的波高空间分布和沿程衰减,礁坪上的平均水位变化,以及波浪能量频谱的移动和空间差异等典型水动力现象开展数值计算。将数值结果与实验结果对比,两者吻合情况良好,验证了模型具有良好的稳定性,具备模拟破碎波浪和海-岸动边界的能力,能较为准确地模拟波浪在礁坪地形上的传播过程中发生的各种水动力现象。     

热带海气相互作用对大气BSISO年际振荡的影响

利用卫星观测OLR资料以及海气耦合数值模拟试验结果,从每年波-频分析结果提取了各种传播模态的强度指数序列,分析了热带北半球夏季季节内振荡(BSISO)各种传播模态的年际变化谱特征,探讨了热带各海区海气相互作用对其影响。主要结果如下:赤道外西传波和印度洋北传波以准2 a为显著振荡周期,赤道东传波、南海北传波和西太平洋北传波则都包含准2 a和准5 a两种周期,南海北传波是5种指数中惟一以准5 a为最主要周期振荡的模态。热带印度洋、西太平洋、东太平洋各海区海气相互作用对各指数准2 a振荡、准5 a振荡既有加强作用,也有削弱作用。各海区比较而言,对赤道东传波准2 a和准5 a振荡、南海北传波准2 a和准5 a振荡起最大加强作用的是西太平洋海区海气相互作用;对赤道外西传波准2 a振荡、西太平洋北传波准2 a和准5 a振荡起最大加强作用的是印度洋海区海气相互作用。     

热带海气相互作用对大气BSISO年际振荡的影响

利用卫星观测OLR资料以及海气耦合数值模拟试验结果,从每年波-频分析结果提取了各种传播模态的强度指数序列,分析了热带北半球夏季季节内振荡(BSISO)各种传播模态的年际变化谱特征,探讨了热带各海区海气相互作用对其影响。主要结果如下:赤道外西传波和印度洋北传波以准2 a为显著振荡周期,赤道东传波、南海北传波和西太平洋北传波则都包含准2 a和准5 a两种周期,南海北传波是5种指数中惟一以准5 a为最主要周期振荡的模态。热带印度洋、西太平洋、东太平洋各海区海气相互作用对各指数准2 a振荡、准5 a振荡既有加强作用,也有削弱作用。各海区比较而言,对赤道东传波准2 a和准5 a振荡、南海北传波准2 a和准5 a振荡起最大加强作用的是西太平洋海区海气相互作用;对赤道外西传波准2 a振荡、西太平洋北传波准2 a和准5 a振荡起最大加强作用的是印度洋海区海气相互作用。     

Application of desingularized approach to water wave propagation over three-dimensional topography

A numerical approach based on desingularized boundary element method and mixed Eulerian–Lagrangian formulation [Zhang et al., 2006. Wave propagation in a fully nonlinear numerical wave tank: a desingularized method. Ocean Engineering 33, 2310–2331] is extended to solve the water wave propagation over arbitrary topography in a three-dimensional wave tank. A robust damping layer applicable for regular and irregular incident waves is employed to minimize the outgoing wave reflection back into the wave tank. Numerical results on the propagation of regular and irregular incident waves over the flat bottom and linear incident waves over an elliptical shoal show good concurrence with the corresponding analytical solutions and experimental data.     

用含高色散性Boussinesq方程模拟潜堤上的波浪变形

基于一种高阶Boussiensq方程(刘忠波等,2004),采用预报-校正格式的有限差分法对该方程进行了数值离散,建立了数值模型。针对动量方程中三阶项的差分形式,采用了迎风格式和五点格式。通过数值模拟常水深下不同周期波浪传播变形,指出迎风格式在计算小周期波浪时存在的问题。为进一步验证数值模型的适用性,模拟了淹没潜堤上的传播变形。从数值结果与实验值的对比结果上看,该数值模型能较好地模拟波浪变形,可用于模拟实际中的波浪场问题。     

斜坡平台波浪破碎数值模拟

波浪在斜坡地形上破碎,破波后稳定波高多采用物理模型试验方法进行研究,利用近岸波浪传播变形的抛物型缓坡方程和波能流平衡方程,导出了适用于斜坡上波浪破碎的数值模拟方法。首先根据波能流平衡方程和缓坡方程基本型式分析波浪在破波带内的波能变化和衰减率,推导了波浪传播模型中波能衰减因子和破波能量流衰减因子之间的关系;其次,利用陡坡地形上的高阶抛物型缓坡方程建立了波浪传播和波浪破碎数学模型;最后,根据物理模型试验实测数据对数值模拟的效果进行验证。数值计算与试验资料比较表明,该模型可以较好地模拟斜坡地形的波浪传播波高变化。     

Beam propagation of tidal internal waves over a submarine slope of the Mascarene Ridge

The generation of internal tides over the Mascarene Ridge is studied on the basis of moored measurements and numerical modeling. The beam structure of the internal wave propagation over a submarine ridge is analyzed. The dependence of the beam propagation of the perturbations on the steepness of the slope, the depth of the ridge crest, and the stratification is studied.     

非线性波传播的新型数值模拟模型及其实验验证

以一种新型的Boussinesq型方程为控制方程组,采用五阶Runge-Kutta-England格式离散时间积分,采用七点差分格式离散空间导数,并通过采用恰当的出流边界条件,从而建立了非线性波传播的新型数值模拟模型.通过对均匀水深水域内波浪传播的数值模拟说明,模型能较好地模拟大水深水域和强非线性波的传播.通过设置不同的入射波参数来进行潜堤地形上波浪传播的物理模型实验,并将数值解与物理模型实验结果进行了比较.     

Practical modified scheme of linear shallow-water equations for distant propagation of tsunamis

A simple but practical numerical model describing a distant propagation of tsunamis is newly proposed by introducing an additional term to the existing modified scheme. The numerical dispersion of the proposed model is manipulated to replace the physical dispersion of the linear Boussinesq equations without any limitation. The new model developed in this study is applied to propagation of a Gaussian hump over a constant water depth and the predicted free surface displacements are compared with available analytical solutions. A very reasonable agreement is observed.     

适于模拟不规则水域波浪的缓坡方程两种数值模型比较

本文分析比较了适于不规则水域波浪模拟的椭圆型缓坡方程两种数值模型。两种数值模型均采用有限体积法离散,分别基于四叉树网格和非结构化三角形网格建立。首先结合近岸缓坡地形上波浪传播的经典物理模型实验对两种数值模型分别进行了验证,并结合计算结果对比分析了两种模型的计算精度和效率。计算结果表明,两种数值模型均可有效地模拟近岸波浪的传播变形;相对非结构化三角形网格下的模型,基于四叉树网格建立的数值模型在数值离散和求解过程中无需引入形函数、不产生复杂的交叉项,离散简单,易于程序实现,且节约计算存储空间,计算效率高。     

Numerical simulation of wave transformation and runup incorporating porous media wave absorber and turbulence models

A numerical wave tank is first established using the Navier–Stokes equations and the VOF method assuming laminar flow. The standard k–ε, realizable k–ε and RNG k–ε turbulent models are then incorporated to the numerical tank. An effective numerical method for wave absorption utilizing the energy-dissipating property of porous media is also included. To validate the accuracy of the proposed models, the propagation of a solitary wave, where analytical solution is available for comparison, is first simulated. This is followed by the simulation of irregular wave runup on a composite seawall, wave propagation over submerged bars and wave refraction and diffraction over an elliptic shoal, where experimental data are available for comparison. All computed results agree well with either the analytical solution or the experimental data.     

A Finite-Difference Approach to the Time-Dependent Mild-Slope Equation

A finite-difference approach is used to develop a time-dependent mild-slope equation incorporating the effects of bottom dissipation and nonlinearity.The Euler predictor-corrector method and the three-point finite-difference method with varying spatial steps are adopted to discretize the time derivatives and the two-dimensional horizontal ones,respectively,thus leading both the time and spatial derivatives to the second-order accuracy.The boundary conditions for the present model are treated on the basis of the general conditions for open and fixed boundaries with an arbitrary reflection coefficient and phase shift.Both the linear and nonlinear versions of the numerical model are applied to the wave propagation and transformation over an elliptic shoal on a sloping beach,respectively,and the linear version is applied to the simulation of wave propagation in a fully open rectangular harbor.From comparison of numerical results with theoretical or experimental ones,it is found that they are in reasonable agreement.     

Numerical modelling of tsunami propagation in the open Black Sea

This paper focuses on the numerical modelling of tsunami propagation in the open Black Sea. Two types of numerical models are discussed: a model for the radial propagation of long waves and an evolutionary finite-difference prognostic model. Experimentally derived numerical data on the model source of tsunamis are reported. Translated by Vladimir A. Puchkin.     

An explicit finite difference model for simulating weakly nonlinear and weakly dispersive waves over slowly varying water depth

In this paper, a modified leap-frog finite difference (FD) scheme is developed to solve Non linear Shallow Water Equations (NSWE). By adjusting the FD mesh system and modifying the leap-frog algorithm, numerical dispersion is manipulated to mimic physical frequency dispersion for water wave propagation. The resulting numerical scheme is suitable for weakly nonlinear and weakly dispersive waves propagating over a slowly varying water depth. Numerical studies demonstrate that the results of the new numerical scheme agree well with those obtained by directly solving Boussinesq-type models for both long distance propagation, shoaling and re-fraction over a slowly varying bathymetry. Most importantly, the new algorithm is much more computationally efficient than existing Boussinesq-type models, making it an excellent alternative tool for simulating tsunami waves when the frequency dispersion needs to be considered.     

Experimental verification of horizontal two-dimensional modified mild-slope equation model

In order to verify modified mild-slope equation models in a horizontal two-dimensional space, a hydraulic experiment is made for surface wave propagation over a circular shoal on which water depth varies substantially. A horizontal two-dimensional numerical model is also constructed based on the hyperbolic equations that have been developed from the modified mild-slope equation to account for the substantial depth variation. Comparison between experimental measurements and numerical results shows that the modified mild-slope equation model is capable of producing accurate results for wave propagation in a region where water depth varies substantially, while the conventional mild-slope equation model gives large errors as the mild-slope assumption is violated.     

Numerical Model of Internal Solitary Wave Evolution on Impermeable Variable Seabed in A Stratified Two-Layer Fluid System

1.IntroductionIn mathematics and physics,a soliton is a self-reinforcing solitary wave caused by nonlinear ef-fectsinthe medium.Solitons are found in many physical phenomena(Chou and Shih,1996;Chouand Quyang,1999;Chouet al.,2003;Chenet al.,2004;Wang,2004;…     

Combined refraction and diffraction of short waves using the finite element method

A two-dimensional finite element numerical model is presented that calculates combined refraction and diffraction of short waves. The wave equation solved governs the propagation of periodic, small amplitude surface gravity waves over a variable depth seabed of mild slope. An efficient computational scheme is employed that allows the solution of practical problems that typically require large computational grids. Comparisons are presented between the finite element model calculations and an analytical solution, a two-dimensional numerical solution, a three-dimensional numerical solution, and measurements from a hydraulic model.     

Numerical simulation of wave overtopping at coastal dikes and low-crested structures by means of a shock-capturing Boussinesq model

In this paper, a shock-capturing numerical model, based on the combined solution of Boussinesq and nonlinear shallow water equations is applied to the simulation of wave runup, overtopping and wave train propagation over impermeable, emerged and low-crested, structures. In order to improve the performances of the scheme at wet–dry interfaces, a numerical treatment is introduced to provide well-balancing of advective fluxes and source terms. One- and two-dimensional test cases are presented to validate the performances of the model under both regular and irregular wave conditions.     

Numerical study on the characteristics of flow field and wave propagation near submerged breakwater on slope

In this study,characteristics of flow field and wave propagation near submerged breakwater on a sloping bed are investigated with numerical model. The governing equations of the vertical twodimensional model are Reynolds Averaged Navier Stokes equations. The Reynolds stress terms are closed by a nonlinear k ε turbulence transportation model. The free surface is traced through the PILC-VOF method. The proposed numerical model is verified with experimental results. The numerical result shows that the wave profile may become more asymmetrical when wave propa-gates over breakwater. When wave crest propagates over breakwater,the anticlockwise vortex may generate. On the contrary,when wave hollow propagates over breakwater,the clockwise vortex may generate. Meanwhile,the influenced zone of vortex created by wave crest is larger than that created by wave hollow. All the maximum values of the turbulent kinetic energy,turbulent dissi-pation and eddy viscosity occur on the top of breakwater. Both the turbulent dissipation and eddy viscosity increase as the turbulent kinetic energy increases. Wave energy may rapidly decrease near the breakwater because turbulent dissipation increases and energy in lower harmonics is transferred into higher harmonics.