基于四阶完全非线性Boussinesq水波方程二维波浪传播数值模型

建立基于四阶完全非线性Boussinesq水波方程的二维波浪传播数值模型。采用Kennedy等提出的涡粘方法模拟波浪破碎。在矩形网格上对控制方程进行离散,采用高精度的数值格式对离散方程进行数值求解。对规则波在具有三维特征地形上的传播过程进行了数值模拟,通过数值模拟结果与实验结果的对比,对所建立的波浪传播模型进行了验证。同时,为了考察非线性对波浪传播的影响,给出和上述模型具有同阶色散性、变浅作用性能但仅具有二阶完全非线性特征的波浪模型的数值结果。通过对比两个模型的数值结果以及实验数据,讨论非线性在波浪传播过程中的作用。研究结果表明,所建立的Boussinesq水波方程在深水范围内不但具有较精确的色散性和变浅作用性能,而且具有四阶完全非线性特征,适合模拟波浪在近岸水域的非线性运动。     

抛物型缓坡方程的变分及数值模拟

对线性水波的折射一绕射问题应用变分原理,对非等深、具有缓坡和不连续的底被导出了一种修改的抛物型缓坡方程近似模型,可预测三维地形上波浪的折射一绕射。同抛物型缓坡方程的线性方程进行了对比。通过数值模拟方法进行数值求解,表明本方法可用于地形条件下的波浪折射一绕射问题。     

数值波浪水槽消波器的改进与检验

数值波浪水槽技术是一种新兴的水波动力学数值模拟技术, 它能够实现水波现象的精确模拟。作为数值波浪水槽的一个重要功能, 消波技术被用于消除传入波浪在水槽末端的反射效应, 以防止反射波对有效实验区域的污染, 从而保证特定要求的水波实验的精确性。目前被广泛采用的消波技术可分为阻尼消波区和主动消波器两类。     

潜堤上破碎波浪传播变形的数值模型及其验证

为模拟潜堤上破碎波浪传播时产生能量的耗散这一特性,在改进的具有四阶色散的Boussinesq水波方程中中入二阶紊动粘性项,建立了考虑波浪破碎的水波数学模型.在非交错网格下建立了有限差分数值模型,并利用三阶Adams-Bash forth格式预报、四阶Adams-Mouton格式校正对数值模型进行求解.通过数值试验,模拟...     

C网格嵌套技术及其在海洋波动传播模拟中的应用

采用粗细嵌套的ArakawaC网格模拟Klevin波和Rossby波沿赤道传播的过程，研究在粗细网格嵌套边界产生数值振荡和反射的原因及其消减方法，得到合理的并能用于复杂海洋模式的ArakawaC网格下的嵌套方案。数值实验结果表明：波形在粗细网格边界产生数值振荡的原因是波形在不同分辨率下的形态有差异，在粗网格下波形趋向于平坦化，而细网格下趋向于锐化。采用双向嵌套和粗细网格交界处加松弛的方法可以有效地消减数值振荡和反射。     

波浪沿斜面传播过程的数值模拟

精确模拟非线性波沿斜面传播过程非常困难,为此论文从势函数的边界积分方程出发,建立了一种时域内二维波浪模拟的数值模型,主要用来模拟完全非线性波浪的传播变形过程。论文的数值模型使用高阶二维边界元方法,采用可调节时间步长的基于二阶显式泰勒展开的混合欧拉-拉格郎日时间步进来求解带自由表面的线性或完全非线性波浪传播问题。在计算区域一端造出线性或非线性的周期性波浪,另一端采用消除反射波的人工粘性吸收边界。通过与现有理论比较证明了论文数值方法所得结果是准确可靠的。     

数值波浪水槽中完全非线性浅水波仿真特性研究

应用基于势流理论的时域高阶边界元方法,建立一个完全非线性的三维数值波浪水槽,通过实时模拟推板造波运动的方式产生波浪。通过混合欧拉-拉格朗日方法和四阶Runge-Kutta方法更新自由水面和造波板的瞬时位置。利用所建模型分别模拟了有限水深波和浅水波,与试验结果、相关文献结果和浅水理论结果吻合较好,且波浪能够稳定传播。系统地讨论造波板的运动圆频率、振幅和水深等对波浪传播和波浪特性的影响,并对波浪的非线性特性进行分析,研究发现造波板运动频率、运动振幅以及水深均将对波浪形态和波浪非线性产生显著影响。结果为真实水槽造波机的运动控制以及波浪生成试验提供了依据,便于实验室设置更合理的参数来准确模拟不同条件下的波浪。     

一种新型三维水流数值模型

以不可压缩流体的N-S方程为基本控制方程,用快速粒子level set方法(FPLS)追踪自由表面,提出了一种新的三维水流数值模型。在自由表面处应用虚拟压力法来封闭压力泊松方程,同时用速度等值外插的方法构造自由表面外侧的虚拟速度分布。通过模拟水波振荡、水柱崩塌、水滴滴落和空箱注水过程证明了模型的有效性。     

欧拉模型模拟波浪的抑制数值耗散效应的一种简单方法

数值波浪水槽是研究波浪及波浪与结构相互作用的常用工具,可在真实尺度下产生波浪,并提供流场的详细数据。然而,大部分数值波浪水槽都存在数值耗散和数值色散问题,数值耗散使波能缓慢消散,数值弥散在波传播过程中使波频移。本文在有限差分法（FDM）求解欧拉方程的基础上,提出了一种抑制数值耗散效应的简单方法,考虑阻尼项的影响,对波的传播解进行了解析求解。该方法的主要思想是在动量方程中附加一个源项,其强度由数值阻尼效应的强度决定。本文通过对规则线性波、Stokes波和不规则波的数值模拟,验证了该方法的有效性。结果表明,本文方法可有效减小数值波浪水槽中存在的数值耗散现象。     

近岸波浪在刚性植被区域传播的数值模型

基于扩展型Boussinesq水波方程,建立了波浪在刚性植被覆盖的近岸海域传播的数值模型。通过在动量方程源项中引入拖曳阻力项考虑植被对波浪的衰减作用。控制方程采用有限差分和有限体积混合格式求解,模型稳定性强,具备间断捕捉能力,能有效模拟近岸区域波浪的传播变形、破碎和处理海岸动边界问题。利用所建立模型对典型物理模型实验进行模拟,计算结果与实验结果吻合良好,表明模型可用于波浪在刚性植被覆盖海域的数值计算。     

A numerical model for predicting the jump of lift force on an oscillating cylinder

The fluid force coefficients on a transversely oscillating cylinder are calculated by applying two-dimensional large eddy simulation method. Considering the “jump” phenomenon of the amplitude of lift coefficient is harmful to the security of the submarine slender structures, the characteristics of this “jump” are dissertated concretely. By comparing with experiment results, we establish a numerical model for predicting the jump of lift force on an oscillating cylinder, providing consultation for revising the hydrodynamic parameters and checking the fatigue life scale design of submarine slender cylindrical structures.     

Solitary wave propagation influenced by submerged breakwater

The form of Boussinesq equation derived by Nwogu （1993） using velocity at an arbitrary distance and surface elevation as variables is used to simulate wave surface elevation changes. In the numerical experiment, water depth was divided into five layers with six layer interfaces to simulate velocity at each layer interface. Besides, a physical experiment was carried out to validate numerical model and study solitary wave propagation.“Water column collapsing”method （WCCM） was used to generate solitary wave. A series of wave gauges around an impervious breakwater were set-up in the flume to measure the solitary wave shoaling, run-up, and breaking processes. The results show that the measured data and simulated data are in good agreement. Moreover, simulated and measured surface elevations were analyzed by the wavelet transform method. It shows that different wave frequencies stratified in the wavelet amplitude spectrum. Finally, horizontal and vertical velocities of each layer interface were analyzed in the process of solitary wave propagation through submerged breakwater.     

Comparative study on the accuracy of solitary wave generations in an ISPH-based numerical wave flume

In the present study, six solitary wave generations by different mathematical approximations are investigated using a piston type wave maker at dimensionless amplitudes ranging from 0.1 to 0.6 and two water depths. Incompressible smoothed particle hydrodynamics is used to simulate solitary wave propagation along the fixed depth channel. The present numerical results are compared with analytical results and experimental data in terms of free surface displacements, fluid particle velocity, phase speed, paddle motion, etc. The present mesh-free numerical results of wave profile variations over time proved that “Rayleigh” has the lowest relative wave height variation. However, its solitary wave has notable phase lead, while “Third order” and “Ninth order” have the least wave lags. Furthermore, the record of present numerical free surface elevation at different distances and the loss of amplitude of the main pulse showed that regarding both of them, “Ninth order” has supremacy over five others. Considering the numerical velocity components of generated solitary wave, “Third order” and “Ninth order” trace analytical results more accurately than other four ones, whereas “Rayleigh” is the most accurate one in predicting the maximum runup. Finally, the paddle motion, its velocity, and displacement, as well as phase speed and outskirts decay coefficient are also compared and discussed intensely.     

三维斜压台风模式 Ⅱ.预报试验

一种斜压多重移动套网格台网模式在国家海洋环境预报中心已被应用于海洋环境预报。本文第一部分已描述了模式方程组和数值方法。本文继续概述模式网格、变分调整初始化和预报试验结果。最外粗网格域固定,内部细网格域随台风中心轨迹移动。模式中,粗细网格变量采用双向耦合。平衡方程和方程,理想台风场和变分调整方案被用于台风模式初始化。一种简单而有效的资料同化方法,即用第6h台风报和弱约束变分原理调整初始场,被应用于提高预报结果。最后本文给出预报试验结果。预报误差统计显示本模式对台风路径预报具有相当能力,而且可以提供海面风和气压场较好的预报。本模式已经与海浪模式联结,得到满意的波高预报结果。     

振荡水柱波能发电装置气室的三维数值模拟研究

振荡水柱装置是目前世界上应用最为广泛的岸式波能发电装置.气室作为该装置的主要结构可将入射波浪的能量转换为往复振荡的气流动能,是完成能量一次转换的关键结构.为了建立用于考察入射波浪、气室内的波面振荡变化,准确预测气室工作性能的三维数值模拟模型,构建了基于VOF模型的三维数值波浪水槽.通过与物理模型试验的结果对比发现,该模...     

The morphodynamic modelling of tidal sand waves on the shoreface

An artificial sand wave on the Dutch shoreface of the North Sea has been studied in conditions with relatively strong tidal currents in the range of 0.5 to 1 m/s and sediments in the medium sand size range of 0.2 to 0.5 mm. The sand wave is perpendicular to the tidal current and has a maximum height and length of the order of 5 m and 1 km, respectively. The sand wave is dynamically active and shows migration rates of the order of a few metres per year. A numerical morphodynamic model (DELFT3D model) has been used to simulate the morphological behaviour of the sand wave in the North Sea. This model approach is based on the numerical solution of the three-dimensional shallow water equations in combination with a surface wave propagation model (wind waves) and the advection–diffusion equation for the sediment particles with online bed updating after each time step. The model results show that the sand wave grows in the case of dominant bed-load transport (weak tidal currents; relatively coarse sediment; small roughness height; low waves) and that the sand wave decays in the case of dominant suspended transport (strong currents, relatively fine sediment, large roughness height; storm waves).     

基于激波捕捉类Boussinesq模型对岛礁地形上波浪传播的数值模拟

为了探究激波捕捉类Boussinesq模型在模拟岛礁地形上规则波和不规则波传播的可行性,采用基于完全非线性Boussinesq方程并具有激波捕捉能力的数值模型Funwave-TVD对规则波和不规则波在岛礁地形上的传播进行了数值模拟,通过与试验数据对比,分析模型中空间步长的影响,验证模型在模拟波高、平均水位分布以及波谱空间演变的能力,结果表明:采用合适的空间步长,模型能较好地模拟规则波和不规则波在岛礁地形上的传播和演化过程。对于规则波,较小的空间步长可改善破碎点处波高峰值的预测,并能更好地预测波浪破碎后波高的空间分布,相比结合经验破碎的Boussinesq模型,Funwave-TVD能更好地模拟规则波在岛礁地形上的破碎,以及破碎以后行进涌波的再生成过程;对于不规则波,Funwave-TVD总体而言能较好地模拟涌浪有效波高、次重力波的生成及空间演化和平均水位,但会低估礁坪上次重力波波高,较粗的空间步长也会低估礁坪上涌浪有效波高。     

On the tuning of a wave-energy driven oscillating-water-column seawater pump to polychromatic waves

Performance of wave-energy devices of the oscillating water column (OWC) type is greatly enhanced when a resonant condition with the forcing waves is maintained. The natural frequency of such systems can in general be tuned to resonate with a given wave forcing frequency. In this paper we address the tuning of an OWC sea-water pump to polychromatic waves. We report results of wave tank experiments, which were conducted with a scale model of the pump. Also, a numerical solution for the pump equations, which were proven in previous work to successfully describe its behavior when driven by monochromatic waves, is tested with various polychromatic wave spectra. Results of the numerical model forced by the wave trains measured in the wave tank experiments are used to develop a tuning criterion for the sea-water pump.     

Numerical study of vegetation damping effects on solitary wave run-up using the nonlinear shallow water equations

Vegetation damping effects on propagating water waves have been investigated by many researchers. This paper investigates the effects of damping due to vegetation on solitary water wave run-up via numerical simulation. The numerical model is based on an implementation of Morison's formulation for vegetation induced inertia and drag stresses in the nonlinear shallow water equations. The numerical model is solved via a finite volume method on a Cartesian cut cell mesh. The accuracy of the numerical scheme and the effects of the vegetation terms in the present model are validated by comparison with experiment results. The model is then applied to simulate a solitary wave propagating on a plane slope with vegetation. The sensitivity of solitary wave run-up to plant height, diameter and stem density is investigated by comparison of the numerical results for different patterns of vegetation. The numerical results show that vegetation can effectively reduce solitary wave propagation velocity and that solitary wave run-up is decreased with increase of plant height in water and also diameter and stem density.     

四叉树网格下的椭圆型缓坡方程数值模型研究

波浪是近岸海域关键的水动力因素之一。考虑到近岸地形复杂、波浪演化显著的特点,建立了四叉树网格体系下的椭圆型缓坡方程数值模型,采用有限体积法对模型进行数值离散,应用GPBiCG(m, n)算法求解离散后的控制方程。模型中根据波浪波长布局计算网格,生成多层次四叉树网格,对复杂计算域有较好的适应性,并且在离散和方程求解中无需引入形函数、不产生复杂的交叉项,节约了存储空间和计算时间。将模型成功应用于物理模型实验及Acapulco海湾的波浪场数值模拟,结果表明该模型能够准确、高效地模拟近岸波浪场,可为近岸波浪场的模拟提供一定的理论和技术支持。