二维水波辐射问题的二阶时域解

用二阶理论在时域范围内研究二维二阶非线性水波的辐射问题,一阶、二阶问题分别满足各自的自由表面条件和物面条件,采用人工阻尼来吸收反射波.计算采用有限元法,对水平半圆柱体在水面的振荡问题进行了计算,计算了波高时间历程和水动力,计算结果与有关文献中的数据相比符合较好.     

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

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

浮式塔的动力计算

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

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

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

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

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

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

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

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

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

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

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

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

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

非线性波浪绕射问题的有限元-无限元解法

本文基于作者在文献中给出的二阶绕射问题的无穷远场条件及无穷远场解,构造了一个适用于二阶绕射问题的无限元。无限元中沿径向(r方向)的插值与二阶绕射问题的远场高阶渐近解的r方向变化规律相同。本文利用内域有限元与外域无限元协调匹配的方法,在二阶逼近的意义下求解了非线性圆柱体绕射问题,给出了波力与波浪爬高的解。本文求出的由二阶势在物面上积分所表达的二阶力部分与Taylor与Hung用一阶势及一个假想的辐射势在自由表面上积分所表达的对应的二阶力部分完全吻合。本文给出的波力与波浪爬高计算结果分别与Chakrabati的波力实验及Isaacson的波浪爬高实验进行了比较,结果吻合得相当好。     

楔形体在波浪中自由入水的数值模拟

物体入水时波浪的影响不可忽略,基于流体力学模型采用VOF法,并利用自定义函数,模拟了楔形体的自由入水过程;同时结合推波板原理及海绵层消波理论实现了数值水槽的造消波,完成了波浪中楔形体自由入水的模拟,计算了楔形体入水时所受的水作用力、自由液面变化及物面压强分布等,研究了不同波高、周期以及在波浪不同位置入水时对楔形体的影响。结果表明:本文建立的数值模型可很好地模拟楔形体入水造成的射流及空泡的形成发展过程,波浪对楔形体入水的影响主要由波浪内部流场变化及表面波形决定,在波浪不同位置处入水对楔形体受力及入水形态均有较大影响。     

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.     

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.     

Interactions between nonlinear water waves and non-wall-sided 3D structures

Interactions between water waves and non-wall-sided cylinders are analyzed based on velocity potential theory with fully nonlinear boundary conditions on the free surface and the body surface. The finite element method (FEM) is adopted together with a 3D mesh generated through an extension of a 2D Delaunay grid on a horizontal plane along the depth. The linear matrix equation for the velocity potential is constructed by imposing the governing equation and boundary conditions through the Galerkin method and is solved through an iterative method. By imposing the gradient of the potential equal to the velocity, the Galerkin method is used again to obtain the velocity field in the fluid domain. Simulations are made for bottom mounted and truncated cylinders with flare in a numerical tank. Periodic waves and wave groups are generated by a piston type wave maker mounted on one end of the tank. Results are obtained for forces, wave profiles and wave runups. Further simulations are made for a cylinder with flare subjected to forced motion in otherwise still open water. Results are provided for surge and heave motion in different amplitudes, and for a body moving in a circular path in the horizontal plane. Comparisons are made in several cases with the results obtained from the second order solution in the time domain.     

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.     

A Numerical Model for Edge Waves on A Compound Slope

An edge wave is a kind of surface gravity wave basically travelling along a shoaling beach. Based on the periodic assumption in the longshore direction, a second order ordinary differential equation is obtained for numerical simulation of the cross-shore surface elevation. Given parameters at the shoreline, a cross-shore elevation profile is obtained through integration with fourth-order Runge-Kutta technique. For a compound slope, a longshore wavenumber is obtained by following a geometrical approach and solving a transcendental equation with an asymptotic method. Numerical results on uniform and compound sloping beaches with different wave periods, slope angles, modes and turning point positions are presented. Some special scenarios, which cannot be predicted by analytical models are also discussed.     

Water entry of a finite width wedge near a floating body

The problem of a two-dimensional finite-width wedge entering water near a freely floating body is considered through the velocity potential theory for the incompressible liquid with the fully nonlinear boundary conditions on the free surface. The problem is solved by using the boundary element method in the time domain. The numerical process is divided into two phases based on whether the interaction between the wedge and floating body is significant. In the first phase, when the single wedge enters water at initial stage, only a small part near its tip is in the fluid, the problem is studied in a stretched coordinate system and the presence of the floating body has no major effect. In the second phase, the disturbance by water entry of the wedge has reached the floating body, and both are considered together in the physical system. The auxiliary function method is adopted to decouple the nonlinear mutual dependence between the motions of the wedge and floating body, both in three degrees of freedom, and the fluid flow, as well as the interaction effects between them. Case studies are undertaken for a wedge entering water in forced or free fall motion, vertically or obliquely. Results are provided for the accelerations, velocities, pressure distribution and free surface deformation, and the interaction effects are discussed.     

Similarity solution of oblique impact of wedge-shaped water column on wedged coastal structures

The initial stage of plunging wave impact obliquely on coastal structures is analysed. The problem is modelled through an oblique collision of an asymmetrical water wedge and an asymmetrical solid wedge. The gravity effect on the flow is ignored based on the assumption that the ratio of the incoming speed of the wave to the acceleration due to gravity is much larger than the time scale of interest. Similarity solution method based on the velocity potential theory is then used. The problem of this similarity flow is solved by a boundary element method through the Cauchy theorem in the complex plan. Results for the wave elevation and pressure distribution are provided, including the forces and moments, effects of different impact angles and the effects of oblique impact are investigated. In particular, negative pressure near the tip of the solid wedge is observed and its implications are discussed.