用时域法求解二维二阶非线性水波

本文用二阶理论在时域范围内计算二维二阶非线性水波,一、二阶问题分别满足各自的自由表面条件和物面条件,流场内的速度势通过求解有限元方程得到,计算采用八结点四边形等参单元,采用人工阻尼来吸收反射波,对楔形体在水面上的振荡问题进行了计算,计算结果与有关文献相比符合较好。     

散射波的二阶幅射边界条件及其应用

本文借助于一个线性微分算子序列,对线性波浪绕射问题,推导出了散射波的一阶和二阶幅射边界条件,作为数值算例,在内域采用有限元,在有限距离幅射边界上分别采用一阶和二阶幅射条件及Sommerfeld条件,计算了作用在大尺度圆柱上的波浪荷载,结果表明,当采用二阶幅射条件时,精度得到很大的提高,因此,可以将内域划得很小,从而避免了以往采用Sommerfeld条件时,为获得足够的精度必须增加单元数目而导致计算量很大和占用很多计算机内存的弊病.     

浮式塔的动力计算

采用频域二阶摄动法研究了浮式塔在给定环境条件下的动力响应.计算了浮式塔在不规则波上的一阶运动响应函数和运动谱以及张力谱,同时进行了相应的模型试验.由结果的比较可以看出,试验结果与理论方法计算得到的结果吻合较好.计算结果表明,波浪的二阶作用对浮式塔的低频纵荡有影响.计算结果与试验结果的比较说明频域二阶摄动法可以用于浮式塔的动力分析.     

基于二阶斯托克斯波理论的辐射应力垂向分布

基于二阶斯托克斯波理论推导了辐射应力的垂向分布表达式,通过算例讨论了辐射应力在深水和有限水深条件下的垂向分布规律,并与基于微幅波理论的辐射应力进行了比较.结果表明,在波浪非线性不强时,基于二阶斯托克斯波理论的辐射应力与基于微幅波理论的辐射应力表达式计算结果接近;而当水深较浅波浪非线性较强时,基于二阶斯托克斯波理论的辐射应力在近表面处明显大于基于微幅波理论的辐射应力.采用二阶斯托克斯波理论推导的波浪辐射应力更为合理地反映了波浪非线性效应.     

三维二阶水波绕射问题的有限元时域计算

本文用有限元法配合时步处理来求解三维非线性水波的绕射问题,自由表面条件和物面条件都满足到二阶,采用人工阻尼区来吸收反射波,流场内的速度势通过求解有限元方程得到。对垂直圆柱体的绕射问题进行了计算,得到了自由表面波高时间历程和圆柱所受到的波浪力,计算结果和有关文献的理论计算结果进行了比较。     

有限水深二维绕射问题的二阶速度势和水动力

本文给出有限水深二维物体二阶绕射势在外域中的解析表达式，从而准确满足二阶绕射势的辐射条件。二阶绕射势在内域自由表面上的边界条件则由一阶势的数值微分求得。然后对内域用简单Ｇｒｅｅｎ函数法求得二阶绕射势。本文对二维浮体和潜体在不同水深和潜深情况下的绕射问题进行了计算，求得了二阶绕射势和二阶定常及倍频波浪力。讨论了水深和潜深对波浪力的影响以及二阶绕射势对非线性波浪力的贡献     

波浪产生和运动的二阶计算模型

在快速模拟波浪运动的谱方法基础上,引入造波边界,建立了模拟波浪产生和运动的二阶计算模型。采用摄动展开方法简化了带有造波边界的水波运动问题,将速度势分解,得到了满足造波边界和自由面边界的速度势的一般解,运用快速Fourier变换和时间积分,建立了模拟波浪产生和运动的数学模型。基于该模型,采用不同的数值造波条件,模拟了波浪的产生问题;考虑了波浪的初始运动问题;通过把数值结果与物理实验的比较,验证了波浪计算模型的有效性。     

水深对软刚臂单点系泊FPSO动力响应的影响

系泊系统的定位能力是浅水油田作业的软刚臂式单点系泊FPSO安全作业的重要保障,为研究不同水深/吃水比下单点系泊系统的受力性能,针对一艘16万吨级软刚臂单点系泊FPSO,在线性三维势流理论的基础上,基于多体动力学方法,建立FPSO-系泊腿-软刚臂的耦合模型,采用Newman近似法和Pinkster近似法分析了FPSO所受二阶波浪力,在时域内计算了不同水深/吃水比对系泊系统动力响应性能的影响。结果表明,随着水深/吃水比的增加,Newman近似法计算得到二阶波浪力先增大后减小,引起单点系泊系统载荷先增大后减小;而Pinkster近似法计算得到的二阶波浪力逐渐减下,引起单点载荷逐渐减下。在浅水条件下,Pinkster近似法具有较好的适用性,Newman近似法严重低估了FPSO所受的二阶波浪力;在深水条件下,Newman近似法能满足工程计算的要求;适用两种方法的临界水深/吃水比为1.64。     

三维物体上二阶波浪绕射的实时模拟

采用二阶时域理论对非线性波浪在任意三维物体周围的绕射问题进行了研究，对自由表面边界条件进行Taylor级数展开，应用摄动展开可以建立相应的边值问题，而且此边值问题的计算域不随时间变化，运用基于B－样条的边界元方法求解每一时刻的波浪场，二阶自由表面边界条件在时间上进行数值积分，在自由表面加了一个人工阻尼层以避免波浪的反射，速度势分解为已知的入射势和未知的散射势，初始条件采用二阶Stokes波浪场，通过加入物体表面边界条件，得到散射势在时间和空间上的发展，本文对圆柱所受规则波的二阶波浪力和波浪爬高进行了计算，数值结果表明此理论计算准确，效率高，数值稳定。     

双色入射波下二阶波浪力响应函数

应用边界元方法,对双频入射波在和频及差频下的二阶速度势做了完整的求解,通过物面积分计算了任意三维结构上的二阶波浪力的传递函数.对简单几何物体,与发表的结果做了对比,两者吻合良好,验证了本方法的正确性.应用这一方法还对复杂的张力腿平台模型做了实际计算,发现在低频和高频区二阶波浪力平方响应函数有着显著的能量分布.     

A composite numerical model for wave diffraction in a harbor with varying water depth

A composite numerical model is presented for computing the wave field in a harbor. The mild slope equation is discretized by a finite element method in the domain concerned. Out of the computational domain, the water depth is assumed to be constant. The boundary element method is applied to the outer boundary for dealing with the infinite boundary condition. Because the model satisfies strictly the infinite boundary condition, more accurate results can be obtained. The model is firstly applied to compute the wave diffraction in a narrow rectangular bay and the wave diffraction from a porous cylinder. The numerical results are compared with the analytical solution, experimental data and other numerical results. Good agreements are obtained. Then the model is applied to computing the wave diffraction in a square harbor with varying water depth. The effects of the water depth in the harbor and the incoming wave direction on the wave height distribution are discussed.     

A composite numerical model for wave diffraction in a harbor with varying water depth

A composite numerical model is presented for computing the wave field in a harbor. The mild slope equation is discretized by a finite element method in the domain concerned. Out of the computational domain, the water depth is assumed to be constant. The boundary element method is applied to the outer boundary for dealing with the infinite boundary condition. Because the model satisfies strictly the infinite boundary condition, more accurate results can be obtained. The model is firstly applied to compute the wave diffraction in a narrow rectangular bay and the wave diffraction from a porous cylinder. The numerical results are compared with the analytical solution, experimental data and other numerical results. Good agreements are obtained. Then the model is applied to computing the wave diffraction in a square harbor with varying water depth. The effects of the water depth in the harbor and the incoming wave direction on the wave height distribution are discussed.     

Second-Order Wave Diffraction Around 3-D Bodies by A Time-Domain Method

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.     

Time stepping solutions of the two-dimensional nonlinear wave radiation problem

The two-dimensional nonlinear time domain free surface flow problem is analysed using potential flow theory. The problem is solved by a time marching method. At each time step two numerical approaches are used. One is based on the boundary element method in the complex plane. The complex potential is assumed to vary linearly within each element and the solution is obtained by imposing the boundary conditions at the nodes of the elements. The other approach is based on the finite element formulation. Triangular elements and linear shape functions are used. The solution is obtained by the Galerkin method. Numerical results are obtained for the wave elevation generated by a vertical wave maker. Results are also provided for a circular cylinder oscillating below the free surface. For these cases the finite element method is found to provide substantially more efficient computations than the boundary element method using equivalent discretizations.     

Fully nonlinear simulations of wave resonance by an array of cylinders in vertical motions

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.     

A Finite Element Solution of Wave Forces on Submerged Horizontal Circular Cylinders

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.     

Fully Nonlinear Simulations of Wave Resonance by An Array of Cylinders in Vertical Motions

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.