首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 718 毫秒
1.
苏高飞  勾莹  滕斌 《海洋工程》2023,41(3):1-13
为高效准确地对完全非线性波浪与二维固定结构物的相互作用进行模拟分析,建立了二维完全非线性时域耦合模型。耦合模型将计算域划分为靠近结构物的内域和远离结构物的外域,每个区域均采用满足完全非线性自由水面边界条件的波浪模型进行求解。在内域使用Laplace方程描述流体运动并采用高阶边界元法(BEM)对其进行求解;而在没有结构物的外域,波浪运动的控制方程为Irrotational Green-Naghdi(IGN)方程并采用有限元法(FEM)对其进行求解。内域和外域通过一段重叠区域进行耦合,从而实现模型间变量的传递。首先利用耦合模型分别对规则波的传播、直墙前立波的生成以及相关物理模型试验进行模拟,数值结果与精确解和试验结果的良好吻合验证了耦合模型耦合方式的合理性以及处理非线性问题的准确性;然后使用耦合模型模拟分析了波浪与固定结构物间的相互作用,并将结果与线性解析解以及完全非线性BEM模型的结果进行了对比分析,进一步证明了耦合模型的正确性与高效性。  相似文献   

2.
系泊船非线性波浪力时域计算: 二维模型   总被引:6,自引:1,他引:6  
王大国  邹志利 《海洋学报》2004,26(2):104-117
为找到具有工程实用价值的港口系泊船波浪力的时域计算方法,建立了在港口中存在系泊船时非线性波浪力时域计算的垂直二维耦合模型:用Boussinesq方程计算船的两侧的外域,用欧拉方程计算船底面下的内域,两域在交界面处的连接条件是流量连续和压力相等.将复平面内的边界元方法应用于所研究问题,对耦合模型进行了验证.进行了相关模型实验,实验结果与数值计算结果比较表明这两种数值计算模型都具有满意的精度,但耦合模型的计算效率要远远高于边界元方法的计算效率.本耦合模型的数学处理简单,可适用于工程计算.  相似文献   

3.
修正型缓坡方程的有限元模型   总被引:1,自引:1,他引:0  
倪云林  滕斌  丛龙飞 《海洋学报》2017,39(1):104-110
与缓坡方程相比,修正型缓坡方程增加了地形曲率项和坡度平方项,从而提高了数值求解的复杂性。本文将计算域划分为内域和外域,内域为水深变化区域,使用修正型缓坡方程,其中的地形曲率项和坡度平方项可用有限单元各节点的水深信息和单元插值函数表示,外域为水深恒定区,速度势满足Helmholtz方程,通过内外域的边界匹配建立有限元方程,并用高斯消去法求解。进而分别模拟了波浪传过Homma岛和圆形浅滩的变形,其结果与相关的解析解和实验数据吻合良好,证明了本文有限元模型的正确性。同时,通过与实验数据的对比也明显看出,在地形坡度较陡的情况下,修正型缓坡方程较缓坡方程具有更高的计算精度。  相似文献   

4.
波浪与大孔隙多孔介质相互作用的耦合数学模型   总被引:1,自引:0,他引:1       下载免费PDF全文
建立了波浪与大孔隙多孔介质相互作用的耦合数学模型,波浪域的控制方程为雷诺时均方程和k-ε紊流模型。对于计算域的入射波采用推板式造波,它可以是线性波、椭圆余弦波和孤立波。采用PLIC-VOF法追踪波浪自由表面。对于多孔介质内的孔隙流场采用非线性Forchheimer方程,两区域共享连续方程,最后导出的波浪域与孔隙流域的压力修正方程具有完全相同的形式,利用这个方程能够同时而不是分别求解波浪场和孔隙流场,避免了在内部边界上给定匹配条件,实现了波浪场与孔隙流场的同步耦合。波浪与粗颗粒海床、平底床面上抛石潜堤及斜坡上抛石潜堤相互作用的验证计算结果表明该模型可用于研究波浪与大孔隙多孔介质相互作用的问题。  相似文献   

5.
滕斌  于梅 《海洋与湖沼》2022,53(4):822-829
针对无限水深下波浪与二维水面物体相互作用问题,传统的波浪格林函数形式复杂、计算缓慢,为了提高计算效率和计算精度,将流域分为物体周围的内域及远离物体的外域,内域采用简单格林函数法,外域采用多极子展开方法,通过内外域边界匹配,耦合求解得到流域中任意一点的速度势,并可计算物体在波浪作用下的波浪激振力、附加质量、辐射阻尼及透射和反射系数。应用该方法计算了二维水面漂浮半圆和水面漂浮方箱的算例,数值计算结果表明,该方法可以方便、准确、快速地计算无限水深下波浪与任意漂浮物体的作用问题。  相似文献   

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

7.
三维完全非线性波浪水槽的数值模拟   总被引:7,自引:0,他引:7  
用有限元求解拉普拉斯方程,建立了三维完全非线性数值波浪水槽.跟踪流体自由表面的方法为满足完全非线性自由表面条件的半拉格朗日法,对离散单元采用20节点的六面体二次等参数单元.并把数值计算结果与水面初始升高产生箱体内流体运动解析解和二阶斯托克斯波理论解进行了对比,结果表明该模型是稳定的、守恒的,能精确模拟非线性波浪的产生和传播.  相似文献   

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

9.
北部湾三维潮流数值模拟   总被引:11,自引:7,他引:11  
模型以经Sigma坐标变换后具有自由表面的三维非线性Navier-Stokes方程为基本方程,以分裂算子法剖分动量方程,用全隐差分格式求解连续方程得到自由表面,最后计算完整的速度场.此外,采用紊流封闭模型求解垂直方向紊动粘滞系数,准确获得了摩擦影响层的潮流结构.应用本模型计算了北部湾的潮波运动,重现了K1、M2分潮潮波系统并揭示了潮流的空间结构特征.  相似文献   

10.
港口中系泊船在波浪作用下运动问题的本质是浅水波浪与浮体的相互作用。与深水情况不同,浅水问题应当考虑水底、水域边界的影响及浅水波浪自身的特性,单一模型很难实现该模拟过程。为此,建立了Boussinesq方程计算入射波和Laplace方程计算散射波的全时域组合计算模型。有限元法求解的Boussinesq方程能使入射波充分考虑到水底、水域边界的影响和浅水波浪的特性;散射波被线性化,采用边界元法求解,并以浮体运动时的物面条件为入射波和散射波求解的匹配条件。该方法为完全的时域方法,计算网格不随时间变动,计算过程较为方便。通过与实验及其他数值方法的结果进行比较,验证了本模型对非线性波面、浮体的运动都有比较理想的计算结果,显示了本模型对非线性问题具有较好的计算能力。  相似文献   

11.
A 2-D time-domain numerical coupled model is developed to obtain an efficient method for nonlinear wave forces on a fixed box-shaped ship in a harbor.The domain is divided into an inner domain and an outer domain.The inner domain is the area beneath the ship and the flow is described by the simplified Euler equations.The other area is the outer domain and the flow is defined by the higher-order Boussinesq equations in order to consider the nonlinearity of the wave motions.Along the interface boundaries between the inner domain and the outer domain,the volume flux is assumed to be continuous and the wave pressures are equal.Relevant physical experiment is conducted to validate the present model.It is shown that the numerical results agree with the experimental data.Compared with the coupled model with the flow in the inner domain governed by the Laplace equation,the present coupled model is more efficient and its solution procedure is more simple,which is particularly useful for the study on the effect of the nonlinear wave forces on a fixed box-shaped ship in a large harbor.  相似文献   

12.
A 3-D time-domain numerical coupled model is developed to obtain an efficient method for nonlinear waves acting on a box-shaped ship fixed in a harbor. The domain is divided into the inner domain and the outer domain. The inner domain is the area beneath the ship and the flow is described by the simplified Euler equations. The remaining area is the outer domain and the flow is defined by the higher-order Boussinesq equations in order to consider the nonlinearity of the wave motions. Along the interface boundaries between the inner domain and the outer domain, the volume flux is assumed to be continuous and the wave pressures are equal. Relevant physical experiment is conducted to validate the present model and it is shown that the numerical results agree with the experimental data. Compared the coupled model with the flow in the inner domain governed by the Laplace equation, the present coupled model is more efficient and its solution procedure is simpler, which is particularly useful for the study on the effect of the nonlinear waves acting on a fixed box-shaped ship in a large harbor.  相似文献   

13.
非线性弱色散波内部流场的重构   总被引:1,自引:0,他引:1  
基于势流理论和级数直接求逆方法,本文建立了基于Bousinesq方程或Green-Naghdi方程给出的水深平均流速或某特征流速及波面信息重构非线性弱色散波内部流场的算法。以Bousinesq方程的孤立波解为例,用本反演方法计算了孤立波的表面水平流速及底部水平流速。结果表明本算法是有效的。本反演算法可用于获取非线性弱色散波的内部流场的详细信息。  相似文献   

14.
基于同量阶迭代法,在保留同阶面的前提下,对林建国等(1998a)得到的二阶Boussinesq类方程进行了求解,得到了与其量阶相对应的取立波解,并春与Euler方程的二阶孤立波解进行了比较,结果显示,本文解比传统Boussinesq方程的孤立波解有明显的改善,扩大了孤立的适用范围。  相似文献   

15.
A fully nonlinear domain decomposed solver is proposed for efficient computations of wave loads on surface piercing structures in the time domain. A fully nonlinear potential flow solver was combined with a fully nonlinear Navier–Stokes/VOF solver via generalized coupling zones of arbitrary shape. Sensitivity tests of the extent of the inner Navier–Stokes/VOF domain were carried out. Numerical computations of wave loads on surface piercing circular cylinders at intermediate water depths are presented. Four different test cases of increasing complexity were considered; 1) weakly nonlinear regular waves on a sloping bed, 2) phase-focused irregular waves on a flat bed, 3) irregular waves on a sloping bed and 4) multidirectional irregular waves on a sloping bed. For all cases, the free surface elevation and the inline force were successfully compared against experimental measurements.  相似文献   

16.
A new form of generalized Boussinesq equations for varying water depth   总被引:1,自引:0,他引:1  
M. Zhao  B. Teng  L. Cheng 《Ocean Engineering》2004,31(16):597-2072
A new set of equations of motion for wave propagation in water with varying depth is derived in this study. The equations expressed by the velocity potentials and the wave surface elevations include first-order non-linearity of waves and have the same dispersion characteristic to the extended Boussinesq equations. Compared to the extended Boussinesq equations, the equations have only two unknown scalars and do not contain spatial derivatives with an order higher than 2. The wave equations are solved by a finite element method. Fourth-order predictor–corrector method is applied in the time integration and a damping layer is applied at the open boundary for absorbing the outgoing waves. The model is applied to several examples of wave propagation in variable water depth. The computational results are compared with experimental data and other numerical results available in literature. The comparison demonstrates that the new form of the equations is capable of calculating wave transformation from relative deep water to shallow water.  相似文献   

17.
The three-dimensional scattering of cnoidal waves by cylinder arrays are studied numerically by using the generalized Boussinesq equations. The boundary-fitted coordinate transformation and a dual-grid technique are used to simplify the finite-difference computation. Also, a set of open boundary conditions and an incident cnoidal wave are incorporated for time-domain simulation. The free-surface elevation and hydrodynamic forces on each cylinder are calculated to illustrate the evolution of nonlinear waves and their interactions with large cylinder arrays. Comparisons are made between the present nonlinear wave loads and those obtained from linear diffraction theory. The sheltering role played by the neighboring cylinders and the feature of wave interference are discussed.  相似文献   

18.
Internal wave generation in an improved two-dimensional Boussinesq model   总被引:1,自引:0,他引:1  
A set of Boussinesq-type equations with improved linear frequency dispersion in deeper water is solved numerically using a fourth order accurate predictor-corrector method. The model can be used to simulate the evolution of relatively long, weakly nonlinear waves in water of constant or variable depth provided the bed slope is of the same order of magnitude as the frequency dispersion parameter. By performing a linearized stability analysis, the phase and amplitude portraits of the numerical schemes are quantified, providing important information on practical grid resolutions in time and space. In contrast to previous models of the same kind, the incident wave field is generated inside the fluid domain by considering the scattered wave field in one part of the fluid domain and the total wave field in the other. Consequently, waves leaving the fluid domain are absorbed almost perfectly in the boundary regions by employment of damping terms in the mass and momentum equations. Additionally, the form of the incident regular wave field is computed by a Fourier approximation method which satisfies the governing equations accurately in water of constant depth. Since the Fourier approximation method requires an Eulerian mean current below wave trough level or a net mass transport velocity to be specified, the method can be used to study the interaction of waves and currents in closed as well as open basins. Several computational examples are given. These illustrate the potential of the wave generation method and the capability of the developed model.  相似文献   

19.
基于物理模型试验,探究畸形波和不规则波作用下浮体系泊张力差异问题。讨论相对波高、相对周期和畸形波参数α1对系泊张力的影响。结果显示:畸形波参数α1和浮体系泊张力显著相关。在α1=2.0~2.83范围内,畸形波作用下迎浪侧系泊张力最大值可达不规则波作用的1.9倍。在相对波高Hs/d=0.032~0.097范围内,畸形波作用下迎浪侧系泊张力最大值显著大于不规则波的作用结果,但畸形波和不规则波对应的1/3值及平均值几乎一致。就相对周期影响而言,迎浪侧系泊张力最大差别出现在谱峰周期Tp0p范围内。频域方面采用小波分析方法讨论畸形波和不规则波作用下浮体系泊张力时频谱特征,两种波浪作用下系泊张力时频特征有显著差别。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号