共查询到18条相似文献，搜索用时 187 毫秒
1.
基于线性势流理论,利用高阶边界元法研究了规则波在三维局部渗透海床上的传播。根据Darcy渗透定律推导出渗透海床的控制方程,利用渗透海床顶部和海底处法向速度和压强连续条件得到渗透海床顶部满足的边界条件。根据绕射理论,利用满足自由水面条件的格林函数建立了求解渗透海床绕射势的边界积分方程,采用高阶边界元方法求解边界积分方程进而得到自由水面的绕射势和波浪在局部渗透海床上传播过程中幅值的变化情况。通过与已发表的波浪对圆柱形暗礁的时域全绕射结果对比,证明了本文建立的频域方法计算波幅的正确性和有效性。利用这一模型研究了三维矩形渗透海床区域上波浪的传播特性,并分析了入射波波长、海床渗透特性系数等参数对波浪传播的影响。 相似文献
2.
本文基于作者在文献中给出的二阶绕射问题的无穷远场条件及无穷远场解,构造了一个适用于二阶绕射问题的无限元。无限元中沿径向(r方向)的插值与二阶绕射问题的远场高阶渐近解的r方向变化规律相同。本文利用内域有限元与外域无限元协调匹配的方法,在二阶逼近的意义下求解了非线性圆柱体绕射问题,给出了波力与波浪爬高的解。本文求出的由二阶势在物面上积分所表达的二阶力部分与Taylor与Hung用一阶势及一个假想的辐射势在自由表面上积分所表达的对应的二阶力部分完全吻合。本文给出的波力与波浪爬高计算结果分别与Chakrabati的波力实验及Isaacson的波浪爬高实验进行了比较,结果吻合得相当好。 相似文献
3.
本文给出有限水深二维物体二阶绕射势在外域中的解析表达式，从而准确满足二阶绕射势的辐射条件。二阶绕射势在内域自由表面上的边界条件则由一阶势的数值微分求得。然后对内域用简单Ｇｒｅｅｎ函数法求得二阶绕射势。本文对二维浮体和潜体在不同水深和潜深情况下的绕射问题进行了计算，求得了二阶绕射势和二阶定常及倍频波浪力。讨论了水深和潜深对波浪力的影响以及二阶绕射势对非线性波浪力的贡献 相似文献
4.
采用二阶时域理论对非线性波浪在任意三维物体周围的绕射问题进行了研究，对自由表面边界条件进行Taylor级数展开，应用摄动展开可以建立相应的边值问题，而且此边值问题的计算域不随时间变化，运用基于B－样条的边界元方法求解每一时刻的波浪场，二阶自由表面边界条件在时间上进行数值积分，在自由表面加了一个人工阻尼层以避免波浪的反射，速度势分解为已知的入射势和未知的散射势，初始条件采用二阶Stokes波浪场，通过加入物体表面边界条件，得到散射势在时间和空间上的发展，本文对圆柱所受规则波的二阶波浪力和波浪爬高进行了计算，数值结果表明此理论计算准确，效率高，数值稳定。 相似文献
5.
本文用二阶理论在时域范围内计算二维二阶非线性水波,一、二阶问题分别满足各自的自由表面条件和物面条件,流场内的速度势通过求解有限元方程得到,计算采用八结点四边形等参单元,采用人工阻尼来吸收反射波,对楔形体在水面上的振荡问题进行了计算,计算结果与有关文献相比符合较好。 相似文献
6.
抛物型缓坡方程的变分及数值模拟 总被引：1，自引：0，他引：1
对线性水波的折射一绕射问题应用变分原理,对非等深、具有缓坡和不连续的底被导出了一种修改的抛物型缓坡方程近似模型,可预测三维地形上波浪的折射一绕射。同抛物型缓坡方程的线性方程进行了对比。通过数值模拟方法进行数值求解,表明本方法可用于地形条件下的波浪折射一绕射问题。 相似文献
7.
8.
水下圆形浅滩附近波浪绕射的计算 总被引：2，自引：0，他引：2
采用波数矢量无旋和波能守恒方程对圆形浅滩附近水域波浪绕射进行了数值计算,计算模型中采用Ｂａｔｔｊｅｓ关系与波数矢量无旋,波能量守恒方程一起联合求解圆形浅滩附近水域波浪折射影响下的波浪要素。本文的数值计算模型对圆形浅滩水域波浪折射绕射现象的验证结果表明,计算所得结果与试验结果是吻合的,数学模型是可靠和合理的,具有实用价值。 相似文献
9.
坐底式平台的沉垫,由于其尺度较大,一般可用绕射问题求解方法中的半解析法和数值计算法,比较准确地获得作用其上的波浪力值。然而在工程实际中,尤其是在平台方案论证和初步设计时,常需要能较快地估算出平台的波浪载荷。因此,本文应用F—K法进行求解,推导出实用的沉垫水平与垂直波浪力的工程计算公式。并且根据一系列模型试验来确定附加的绕射效应和质量效应,据此编制程序的计算结果与国外发表的有关资料相当接近。文中给出了一个平台的计算实例。此外,本文还对沉垫垂直波浪力对平台坐底稳性的影响进行了分析。 相似文献
10.
11.
12.
The finite element method (FEM) is employed to analyze the resonant oscillations of the liquid confined within multiple or an array of floating bodies with fully nonlinear boundary conditions on the free surface and the body surface in two dimensions. The velocity potentials at each time step are obtained through the FEM with 8-node quadratic shape functions. The finite element linear system is solved by the conjugate gradient (CG) method with a symmetric successive overelaxlation (SSOR) preconditioner. The waves at the open boundary are absorbed by the combination of the damping zone method and the Sommerfeld-Orlanski equation. Numerical examples are given by an array of floating wedge- shaped cylinders and rectangular cylinders. Results are provided for heave motions including wave elevations, profiles and hydrodynamic forces. Comparisons are made in several cases with the results obtained from the second order solution in the time domain. It is found that the wave amplitude in the middle region of the array is larger than those in other places, and the hydrodynamic force on a cylinder increases with the cylinder closing to the middle of the array. 相似文献
13.
A time-domain method is applied to simulate nonlinear wave diffraction around a surface piercing 3-D arbitrary body. The method involves the application of Taylor series expansions and the use of perturbation procedure to establish the corresponding boundary value problems with respect to a time-independent fluid domain. A boundary element method based on B-spline expansion is used to calculate the wave field at each time step, and the free surface boundary condition is satisfied to the second order of wave steepness by a numerical integration in time. An artificial damping layer is adopted on the free surface for the removal of wave reflection from the outer boundary. As an illustration, the method is used to compute the second-order wave forces and run-up on a surface-piercing circular cylinder. The present method is found to be accurate, computationally efficient, and numerically stable. 相似文献
14.
The finite element method(FEM) is employed to analyze the resonant oscillations of the liquid confined within multiple or an array of floating bodies with fully nonlinear boundary conditions on the free surface and the body surface in two dimensions.The velocity potentials at each time step are obtained through the FEM with 8-node quadratic shape functions.The finite element linear system is solved by the conjugate gradient(CG) method with a symmetric successive overelaxlation(SSOR) preconditioner.The waves at the open boundary are absorbed by the combination of the damping zone method and the Sommerfeld-Orlanski equation.Numerical examples are given by an array of floating wedgeshaped cylinders and rectangular cylinders.Results are provided for heave motions including wave elevations,profiles and hydrodynamic forces.Comparisons are made in several cases with the results obtained from the second order solution in the time domain.It is found that the wave amplitude in the middle region of the array is larger than those in other places,and the hydrodynamic force on a cylinder increases with the cylinder closing to the middle of the array. 相似文献
15.
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. 相似文献
16.
A full time-domain analysis program is developed for the coupled dynamic analysis of offshore structures. For the hydrodynamic loads, a time domain second order method is developed. In this approach, Taylor series expansions are applied to the body surface and free-surface boundary conditions, and the Stokes perturbation procedure is then used to establish the corresponding boundary value problems with time-independent boundaries. A higher-order boundary element method (HOBEM) is developed to calculate the velocity potential of the resulting flow field at each time step. The free-surface boundary condition is satisfied to the second order by fourth order Adams–Bashforth–Moultn method. An artificial damping layer is adopted on the free surface to avoid the wave reflection. The mooring-line/tendon/riser dynamics are based on the rod theory and the finite element method (FEM), with the governing equations described in a global coordinate system. In the coupled dynamic analysis, the motion equation for the hull and dynamic equations for mooring-lines/tendons/risers are solved simultaneously using the Newmark method. The coupled analysis program is applied for a truss Spar motion response simulation. Numerical results including motions and tensions at the top of mooring-lines/risers are presented, and some significant conclusions are derived. 相似文献
17.
利用完全非线性数值波浪水槽技术研究水下平板与波浪的相互作用。假定水下平板厚度极薄、刚性,位于有限水深并且非常接近自由水面。应用四阶龙格库塔方法追踪每一时刻的波面形状,采用阻尼层来吸收反射波以保证算法的稳定性,同时引入平滑和重组的方法抑制自由表面控制点的较高梯度。通过对波浪与浮动圆柱相互作用的数值模拟证实了数值波浪水槽方法的有效性,计算结果与线性理论吻合良好。在波浪数值水槽方法中引入造波板模拟波浪产生并与水下平板发生相互作用,应用傅立叶解析方法对波面变形、波浪力作了分析。结果表明在板非常接近自由水面的情况下会表现出现很强的非线性,揭示了线性理论的局限性。 相似文献
18.
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. 相似文献