共查询到18条相似文献，搜索用时 109 毫秒
1.
建立基于四阶完全非线性Boussinesq水波方程的二维波浪传播数值模型。采用Kennedy等提出的涡粘方法模拟波浪破碎。在矩形网格上对控制方程进行离散,采用高精度的数值格式对离散方程进行数值求解。对规则波在具有三维特征地形上的传播过程进行了数值模拟,通过数值模拟结果与实验结果的对比,对所建立的波浪传播模型进行了验证。同时,为了考察非线性对波浪传播的影响,给出和上述模型具有同阶色散性、变浅作用性能但仅具有二阶完全非线性特征的波浪模型的数值结果。通过对比两个模型的数值结果以及实验数据,讨论非线性在波浪传播过程中的作用。研究结果表明,所建立的Boussinesq水波方程在深水范围内不但具有较精确的色散性和变浅作用性能,而且具有四阶完全非线性特征,适合模拟波浪在近岸水域的非线性运动。 相似文献
2.
为高效准确地对完全非线性波浪与二维固定结构物的相互作用进行模拟分析，建立了二维完全非线性时域耦合模型。耦合模型将计算域划分为靠近结构物的内域和远离结构物的外域，每个区域均采用满足完全非线性自由水面边界条件的波浪模型进行求解。在内域使用Laplace方程描述流体运动并采用高阶边界元法（BEM）对其进行求解；而在没有结构物的外域，波浪运动的控制方程为Irrotational Green-Naghdi（IGN）方程并采用有限元法（FEM）对其进行求解。内域和外域通过一段重叠区域进行耦合，从而实现模型间变量的传递。首先利用耦合模型分别对规则波的传播、直墙前立波的生成以及相关物理模型试验进行模拟，数值结果与精确解和试验结果的良好吻合验证了耦合模型耦合方式的合理性以及处理非线性问题的准确性；然后使用耦合模型模拟分析了波浪与固定结构物间的相互作用，并将结果与线性解析解以及完全非线性BEM模型的结果进行了对比分析，进一步证明了耦合模型的正确性与高效性。 相似文献
3.
4.
无限水深聚焦波完全非线性数值模拟 总被引：1，自引：1，他引：0
基于势流理论提出一种新的高阶边界元方法对无限水深的聚焦波浪进行完全非线性数值模拟.自由水面满足完全非线性边界条件,模拟波浪的非线性效果可以达到更高阶.利用镜像原理,建立一种全新的格林函数应用到无限水深的数值波浪水槽中,以致于两无限深水槽侧壁的积分可以被排除.为了产生相应的入射波和吸收出流波浪,一个由点源组成的造波装置被布置于计算域内,同时人工阻尼层被用来吸引出流波浪,由波浪聚焦的方法得到极限波浪.通过开展线性和完全非线性聚焦波浪的数值实验及与理论解对比,验证本数值模型可以用来模拟无限深水域的极限波浪,且在出流边界没有反射. 相似文献
5.
6.
基于高阶边界元的三维数值波浪港池 总被引：8，自引：1，他引：8
初步建立了一个基于高阶边界元的三维数值波浪港池，港池具有造波和消波功能。采用高阶边界元16节点四边形单元和基于二阶显式泰勒展开的混合欧拉-拉格朗日时间步进求解带自由表面的完全非线性势流方程。模型中对于影响数值精度的问题作了细致的处理。数值计算结果表明本港池可以用来模拟非线性波浪的传播，具有很高的数值精度和稳定性。 相似文献
7.
基于高阶边界元的三维数值波浪港池--波浪破碎的模拟 总被引：5，自引：1，他引：4
在势流理论的框架内,采用高阶边界元方法和混合欧拉-拉格朗日法,实现了对三维波浪破碎过程的数值模拟.数值模型使用可调节时间步长的基于二阶显式泰勒展开的混合欧拉-拉格郎日时间步进来求解自由表面的演化过程.在所使用的边界元方法中,采用16节点三次滑移四边形单元来表示,这种单元在单元内具有高阶的精度同时在单元之间具有良好的连续性.给出了孤立波的传播和周期性非线性波浪沿缓坡传播的计算结果,表明数值模型具有良好的稳定性. 相似文献
8.
9.
10.
本文基于具备间断捕捉能力的二阶全非线性Boussinesq数值模型,对规则波和随机波在礁坪地形上的传播变形进行了数值模拟。该模型采用高阶有限体积法和有限差分方法求解守恒格式的控制方程,将波浪破碎视为间断,同时采用静态重构技术处理了海岸动边界问题。重点针对礁坪上波浪传播过程中的波高空间分布和沿程衰减,礁坪上的平均水位变化,以及波浪能量频谱的移动和空间差异等典型水动力现象开展数值计算。将数值结果与实验结果对比,两者吻合情况良好,验证了模型具有良好的稳定性,具备模拟破碎波浪和海-岸动边界的能力,能较为准确地模拟波浪在礁坪地形上的传播过程中发生的各种水动力现象。 相似文献
11.
Simulation of Fully Nonlinear 3-D Numerical Wave Tank 总被引：6，自引：0，他引：6
A fully nonlinear numerical wave tank (NWT) has been simulated by use of a three-dimensional higher order bouodary element method (HOBEM) in the time domain. Within the frame of potential flow and the adoption of simply Rankine source, the resulting boundary integral equation is repeatedly solved at each time step and the fully nonlinear free surface boundary conditions are integrated with time to update its position and boundary values. A smooth technique is also adopted in order to eliminate the possible saw-tooth numerical instabilities. The incident wave at the uptank is given as theoretical wave in this paper. The outgoing waves are absorbed inside a damping zone by spatially varying artificial damping on the free surface at the wave tank end. The numerical results show that the NWT developed by these approaches has a high accuracy and good numerical stability. 相似文献
12.
Freak waves are generated based on the mechanism of wave focusing in a 2D numerical wave tank. To set up the nonlinear numerical wave tank, the Boundary Element Method is used to solve potential flow equations incorporated with fully nonlinear free surface boundary conditions. The nonlinear properties of freak waves, such as high frequency components and wave profile asymmetry, are discussed. The kinematic data, which can be useful for the evaluation of the wave forces exerted on structures to avoid underestimation of linear predictions, are obtained, and discussed, from the simulated results of freak waves. 相似文献
13.
Numerical Wave Channel with Absorbing Wave-Maker 总被引：2，自引：0，他引：2
Wang Yongxue
Professor State Key Laboratory of Coastal Offshore Engineering Dalian University of Technology Dalian 《中国海洋工程》1995,(2)
The numerical wave channel has been developed based on the volume of fluid method (VOF) in conjunction with the Navier-Stokes equations. The absorbing wave-maker boundary on the left side of the channel is presented by prescribing velocity reference to linear wave-maker theory. The principle of which is that the numerical wave-maker is designed to move in a way that generates the required incident wave and cancels out any reflected wave that reach it at the same time. On the right side of the channel, the open boundary is set to permit incident waves to be transmitted freely. The parametric studies have been carried out at a range of ratios of water depth to wave length d/ L from 0.124 to 0.219, with wave height in the front of paddle/water depth ratio (H0 / d) from 0.1 to 0.3. Wave height, wave pressure distribution along the channel and velocity field are obtained for both open boundary condition and reflective boundary condition at the other end of the channel. For a reflective case, it is shown that 相似文献
14.
A finite-difference scheme and a modified marker-and-cell (MAC) algorithm have been developed to investigate the interactions of fully nonlinear waves with two- or three-dimensional structures of arbitrary shape. The Navier–Stokes (NS) and continuity equations are solved in the computational domain and the boundary values are updated at each time step by the finite-difference time-marching scheme in the framework of a rectangular coordinate system. The fully nonlinear kinematic free-surface condition is implemented by the marker-density function (MDF) technique developed for two fluid layers.To demonstrate the capability and accuracy of the present method, the numerical simulation of backstep flows with free-surface, and the numerical tests of the MDF technique with limit functions are conducted. The 3D program was then applied to nonlinear wave interactions with conical gravity platforms of circular and octagonal cross-sections. The numerical prediction of maximum wave run-up on arctic structures is compared with the prediction of the Shore Protection Manual (SPM) method and those of linear and second-order diffraction analyses based on potential theory and boundary element method (BEM). Through this comparison, the effects of non-linearity and viscosity on wave loading and run-up are discussed. 相似文献
15.
非线性波浪波面追踪的一种新模式 总被引：1，自引：0，他引：1
基于Laplace方程的Green积分表达式和波面BemouUi方程所建立的非线性波动数学模型,是一个时域上具有初始值的边值问题,而精确地追踪自由表面的波动位置,给出波面运动瞬时的波面高度和波面势函数,是建立时域内非线性波浪数值模式的基础。本文采用0-1混合型边界元剖分计算域边界并离散Laplace方程的Green积分表达式,采用有限元剖分自由水面并推导满足自由表面非线性边界条件的波面有限元方程,联立计算域内以节点波势函数和波面位置高度的时间增量为未知量的线性方程组,通过时步内的循环迭代,给出每个时步上的波面位置和波面势函数,从而建立了一种新的非线性波浪波面追踪模式。数值造波水槽内的波浪试验表明,其数值模拟结果具有良好的计算精度。 相似文献
16.
In this paper, a numerical model is established for estimating the wave forces on a submerged horizontal circular cylinder. For predicting the wave motion, a set of two-dimensional Navier-Stokes equations is solved numerically with a finite element method. In order to track the moving non-linear wave surface boundary, the Navier-Stokes equations are discretized in a moving mesh system. After each computational time step, the mesh is modified according to the changed wave surface boundary. In order to stabilize the numerical procedure, a three-step finite element method is applied in the time integration. The water sloshing in a tank and wave propagation over a submerged bar are simulated for the first time to validate the present model. The computational results agree well with the analytical solution and the experimental data.Finally, the model is applied to the simulation of interaction between waves and a submerged horizontal circular cylinder.The effects of the KC number and the cylinder depth on the wave forces are studied. 相似文献
17.
基于线性势流理论,利用高阶边界元法研究了规则波在三维局部渗透海床上的传播。根据Darcy渗透定律推导出渗透海床的控制方程,利用渗透海床顶部和海底处法向速度和压强连续条件得到渗透海床顶部满足的边界条件。根据绕射理论,利用满足自由水面条件的格林函数建立了求解渗透海床绕射势的边界积分方程,采用高阶边界元方法求解边界积分方程进而得到自由水面的绕射势和波浪在局部渗透海床上传播过程中幅值的变化情况。通过与已发表的波浪对圆柱形暗礁的时域全绕射结果对比,证明了本文建立的频域方法计算波幅的正确性和有效性。利用这一模型研究了三维矩形渗透海床区域上波浪的传播特性,并分析了入射波波长、海床渗透特性系数等参数对波浪传播的影响。 相似文献
18.
A fully nonlinear numerical model based on a time-domain higher-order boundary element method (HOBEM) is founded to simulate the kinematics of extreme waves. In the model, the fully nonlinear free surface boundary conditions are satisfied and a semi-mixed Euler-Lagrange method is used to track free surface; a fourth-order Runga-Kutta technique is adopted to refresh the wave elevation and velocity potential on the free surface at each time step; an image Green function is used in the numerical wave tank so that the integrations on the lateral surfaces and bottom are excluded. The extreme waves are generated by the method of wave focusing. The physical experiments are carried out in a wave flume. On the horizontal velocity of the measured point, numerical solutions agree well with experimental results. The characteristics of the nonlinear extreme-wave kinematics and the velocity distribution are studied here. 相似文献