结合椭圆型缓坡方程模拟近岸波流场

波浪向近岸传播的过程中,由波浪破碎效应所产生的近岸波流场是近岸海域关键的水动力学因素之一.结合近岸波浪场的椭圆型缓坡方程和近岸波流场数学模型对近岸波浪场及由斜向入射波浪破碎后所形成的近岸波流场进行了数值模拟.计算中考虑到波浪向近岸传播中由于波浪的折射、绕射、反射等效应使局部复杂区域波向不易确定,采用结合椭圆型缓坡方程所给出的波浪辐射应力公式来计算波浪产生的辐射应力,在此基础上耦合椭圆型缓坡方程和近岸波流场数学模型对近岸波流场进行数值模拟,从而使模型综合考虑了波浪的折射、绕射、反射等效应且避免了对波向角的直接求解,可以应用于相对较复杂区域的近岸波流场模拟.     

非结构化网格下近岸波生流数值模拟

波浪破碎产生的近岸流是近岸海域关键的水动力因素之一。基于近岸波浪的椭圆型缓坡方程和二维近岸波生流方程,建立了非结构化网格下近岸波浪破碎形成的近岸流数值模型。数值模型中,在空间上采用有限体积法进行数值离散,在时间上采用欧拉向前格式数值离散。数值计算结果表明,该数值模型可以有效地模拟近岸波浪破碎产生的近岸流。     

斜坡上波浪破碎与越浪非静压数值模拟

为合理确定防浪建筑物的越浪量,基于含非静水压力梯度项的非线性浅水方程建立了近岸波浪越浪数值模型。通过采用域内造波、消波并结合波前静压假定的破碎模型,模拟了规则波和不规则波在斜坡上的波浪传播变形,并在此基础上进行了越浪量数值计算。数值计算结果与物理模型实验结果表明,非静压模型可合理地描述波浪破碎点位置、破碎后的波高、增减水以及斜坡上的堤后越浪量。数值模型具有较高的计算精度和计算效率,可为实际工程防浪建筑物越浪以及堤顶高程的设计提供一种新的数值研究手段。     

非结构化网格下大范围波生流模拟和应用

波浪破碎引起的沿岸流是近岸海域的关键水动力因素。利用基于缓坡方程得到的光程函数方程和波作用守恒方程建立了考虑绕射效应的大范围波浪传播模型,模型可以考虑流场的影响;将波浪模型计算得到的辐射应力、波浪紊动系数等参数添加到三维水动力模型中,得到大范围近岸波生流的计算模型。模型中流场和波浪可以共用计算网格,且可同步耦合计算;模型基于非结构化网格,可以拟合复杂岸线的变化。模型对波生沿岸流、环流和逆流进行了验证,同时对实际海域的波生流进行了计算,结果表明:该模型对近岸波浪破碎引起的波生流具有很好的精度和适用性,可用于实际工程的计算。     

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

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

基于抛物型缓坡方程模拟近岸植被区波浪传播

植被对波浪传播运动有重要影响。考虑近岸波浪在植被区传播中的折射、绕射、破碎及植被引起的波能耗损效应,基于抛物型缓坡方程建立了模拟近岸植被区波浪传播的数学模型,对模型进行了数值模拟验证,采用数值模拟试验分析了植被对波浪传播的影响。数值模拟结果表明,波浪在近岸植被区传播时,随着植被密度和植被高度的增加,波浪传播中的波高衰减增大,波能耗损增加;不同周期波浪在植被区传播中的波高衰减过程也明显不同。     

近岸波浪破碎区不规则波浪的数值模拟

基于近岸不规则波浪传播的抛物型缓坡方程和两类波浪破碎能量损耗因子,对近岸波浪破碎区不规则波浪的波高分布进行了数值模拟,并结合实验结果对数值模拟结果进行了验证分析,结果表明采用两类波浪破碎能量损耗因子所模拟的破碎区波高与实测值均吻合良好,波浪破碎能量损耗因子及波浪破碎指标对破碎区波浪波高分布影响较明显。     

潜堤上规则波辐射应力的数值研究

在近岸波浪相关研究中,辐射应力是波动在水体中引起的剩余动量流,是波浪运动的重要物理量.在波浪从深水逐渐传向浅水的过程中,波浪的非线性逐渐增强,甚至会发生破碎等剧烈变形,引起辐射应力的强烈变化,对次重力波生成等有重要贡献.应用OpenFOAM精细模拟波浪在潜堤上的传播,得出波浪运动的详细流场信息,计算了有波浪破碎情况下潜...     

波浪漫滩边界对波生流数值计算的影响

波浪漫滩是近岸波浪的小尺度运动,在实际海域的波生流数值计算中通常被忽略。本文基于Boussinesq方程的FUNWAVE模式,分别采用波浪漫滩边界、固壁边界、海绵边界进行Haller波浪港池物理模型实验的数值模拟,比较三种边界计算结果与实验观测数据的误差,检验波浪漫滩边界对波生流数值计算的影响;然后设计了多种周期、波高的波生流数值模拟试验,分析多种波浪入射条件下波浪漫滩边界对近岸波生流数值计算的影响。结果表明,波浪漫滩对邻近区域波生流有明显影响,漫滩边界下的波生流计算结果更接近实验观测值,在近岸波生流数值模型中引入波浪漫滩边界可以提高波生流计算精度。     

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

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

Numerical Simulation of Wave Height and Wave Set-Up in Nearshore Regions

Based on the time dependent mild slope equation including the effect of wave energy dissipation, an expression for the energy dissipation factor is derived in conjunction with the wave energy balance equation, and then a practical method for the simulation of wave height and wave set-up in nearshore regions is presented. The variation of the complex wave amplitude is numerically simulated by use of the parabolic mild slope equation including the effect of wave energy dissipation due to wave breaking. The components of wave radiation stress are calculated subsequently by new expressions for them according to the obtained complex wave amplitude, and then the depth-averaged equation is applied to the calculation of wave set-up due to wave breaking. Numerical results are in good agreement with experimental data, showing that the expression for the energy dissipation factor is reasonable and that the new method is effective for the simulation of wave set-up due to wave breaking in nearshore regions.     

Two-dimensional wave-current interaction with a locally enriched quadtree grid system

Wave-induced currents may drive nearshore mixing and transport processes, including coastal pollutant dispersion, littoral drift, and long-term morphological changes through beach erosion and accretion. In this study, a numerical model is newly developed to simulate wave climate and localized currents in complicated coastal environments. The model developed is based on a quadtree grid system. The two-dimensional hydrodynamic governing equations are solved by using an explicit Adams-Bashforth finite difference scheme. Effects of wave breaking, shoaling, refraction, diffraction, wave-current interaction, set-up and set-down, turbulent mixing, bed friction, and shoreline movement are incorporated in the model. Results are presented for set-up at a beach in a flume due to normally incident waves, and longshore currents generated by oblique waves on a plane beach.     

Numerical Simulation of Breaking Wave Based on Higher-Order Mild Slope Equation

The "surface roller" to simulate wave energy dissipation of wave breaking is introduced into the random wave model based on approximate parabolic mild slope equation in this paper to simulate the random wave transportation including diffraction, refraction and breaking in nearshore areas. The roller breaking random wave higher-order approximate parabolic equation model has been verified by the existing experimental data for a plane slope beach and a circular shoal, and the numerical results of random wave breaking model agree with the experimental data very well. This model can be applied to calculate random wave propagation from deep to shallow water in large areas near the shore over natu ral topography.     

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.     

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

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

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

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

Numerical solution of a mathematical model for water wavesin large coastal areas

on the evolution equation for water waves,a mathematical model for wave propagation in large mild-slope areas is derived.The model is solved by the finite difference method with the staggered grid system.The computational results are in good agreement with experimental data and show that the model can obtain better results with relatively coarser grids.The model can be used to simulate water wave propagation in large coastal areas and can be efficiently solved without much programming effort.     

近岸大区域水波数学模型及其数值求解

从水波发展方程出发,导出了大区域缓坡水波数学模型,并采用交错网格系统下的有限差分法对该数学模型进行了数值求解,计算结果表明该数学模型在粗网格下也能得到与实验结果比较一致的结果,从而表明该模型可用于较大区域的水波问题.该模型具有编程简单、求解比较快速、经济的优点.     

近岸沿岸流及污染物运动的数值模拟

基于双曲型缓坡方程和近岸浅水方程对近岸波浪斜向入射破碎所生成的沿岸流及污染物在沿岸波流作用下的运动进行了数值模拟,并对数值模拟结果进行了验证分析。数值模拟结果表明,在相近工况参数下,随着入射波高的增大,沿岸流流速和平均水面升高值均明显增大;随着岸坡坡度的增加,沿岸流流速和平均水面升高值明显增大;随着入射波浪周期的增大,平均水面升高值明显增大。在沿岸缓坡区域,由斜向入射波浪破碎所产生的沿岸流对污染物的运动起着不可忽略的影响。     

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.