完全非线性波浪与二维固定结构物作用的IGN-BEM耦合分析模型

为高效准确地对完全非线性波浪与二维固定结构物的相互作用进行模拟分析,建立了二维完全非线性时域耦合模型。耦合模型将计算域划分为靠近结构物的内域和远离结构物的外域,每个区域均采用满足完全非线性自由水面边界条件的波浪模型进行求解。在内域使用Laplace方程描述流体运动并采用高阶边界元法（BEM）对其进行求解;而在没有结构物的外域,波浪运动的控制方程为Irrotational Green-Naghdi（IGN）方程并采用有限元法（FEM）对其进行求解。内域和外域通过一段重叠区域进行耦合,从而实现模型间变量的传递。首先利用耦合模型分别对规则波的传播、直墙前立波的生成以及相关物理模型试验进行模拟,数值结果与精确解和试验结果的良好吻合验证了耦合模型耦合方式的合理性以及处理非线性问题的准确性;然后使用耦合模型模拟分析了波浪与固定结构物间的相互作用,并将结果与线性解析解以及完全非线性BEM模型的结果进行了对比分析,进一步证明了耦合模型的正确性与高效性。     

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.     

三维完全非线性波浪水槽的数值模拟

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

An Efficient Model for Transient Surface Waves in Both Finite and Infinite Water Depth

A numerical model is developed to simulate fully nonlinear extreme waves in finite and infinite water-depth wave tanks. A semi-mixed Eulerian-Lagrangian formulation is adopted and a higher-order boundary element method in conjunction with an image Green function is used for the fluid domain. The boundary values on the free surface are updated at each time step by a fourth-order Runga-Kutta time-marching scheme at each time step. Input wave characteristics are specified at the upstream boundary by an appropriate wave theory. At the downstream boundary, an artificial damping zone is used to prevent wave reflection back into the computational domain. Using the image Green function in the whole fluid domain, the integrations on the two lateral walls and bottom are excluded. The simulation results on extreme wave elevations in finite and infinite water-depths are compared with experimental results and second-order analytical solutions respectively. The wave kinematics is also discussed in the present study.     

An Efficient Model for Transient Surface Waves in Both Finite and Infinite Water Depths

A numerical model is developed to simulate fully nonlinear extreme waves in finite and infinite water-depth wave tanks. A semi-mixed Eulerian-Lagrangian formulation is adopted and a higher-order boundary element method in conjunction with an image Green function is used for the fluid domain. The boundary values on the free surface are updated at each time step by a fourth-order Runga-Kutta time-marching scheme at each time step. Input wave characteristics are specified at the upstream boundary by an appropr...     

Wave Refraction-Diffraction in the Presence of a Current

-Wave refraction-diffraction due to a large ocean structure and topography in the presence of a 'current are studied numerically. The mathematical model is the mild-slope equation developed by Kirby (1984). This equation is solved using a finite and boundary element method. The physical domain is devid-ed into two regions: a slowly varying topography region and a constant water depth region. For waves propagating in the constant water depth region, without current interfering, the mild- slope equation is then reduced to the Helmholtz equation which is solved by boundary element method. In varying topography region, this equation will be solved by finite element method. Conservation of mass and energy flux of the fluid between these two regions is required for composition of these two numerical methods. The numerical scheme proposed here is capable of dealing with water wave problems of different water depths with the main characters of these two methods.     

Simulation of nearshore wave processes by a depth-integrated non-hydrostatic finite element model

This paper presents CCHE2D-NHWAVE, a depth-integrated non-hydrostatic finite element model for simulating nearshore wave processes. The governing equations are a depth-integrated vertical momentum equation and the shallow water equations including extra non-hydrostatic pressure terms, which enable the model to simulate relatively short wave motions, where both frequency dispersion and nonlinear effects play important roles. A special type of finite element method, which was previously developed for a well-validated depth-integrated free surface flow model CCHE2D, is used to solve the governing equations on a partially staggered grid using a pressure projection method. To resolve discontinuous flows, involving breaking waves and hydraulic jumps, a momentum conservation advection scheme is developed based on the partially staggered grid. In addition, a simple and efficient wetting and drying algorithm is implemented to deal with the moving shoreline. The model is first verified by analytical solutions, and then validated by a series of laboratory experiments. The comparison shows that the developed wave model without the use of any empirical parameters is capable of accurately simulating a wide range of nearshore wave processes, including propagation, breaking, and run-up of nonlinear dispersive waves and transformation and inundation of tsunami waves.     

浅水中浮体在波浪作用下运动的三维时域计算模型

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

Coupled dynamic analysis for wave interaction with a truss spar and its mooring line/riser system in time domain

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.     

Internal wave generation in an improved two-dimensional Boussinesq model

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.     

Internal wave generation in an improved two-dimensional Boussinesq model

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.     

系泊船非线性波浪力时域计算:二维模型

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

Transient Hydroelastic Response of VLFS by FEM with Impedance Boundary Conditions in Time Domain

A time-dependent finite element method (FEM) is developed to analyze the transient hydroelastie responses of very large floating structures (VLFS) subjected to dynamic loads. The hydrodynamic problem is formulated based on the linear theory of fluid and the structural response is analyzed based on the thin plate theory. The FEM truncates the unbounded fluid domain by introducing an artificial boundary surface, thus defining a finite computational domain. At this boundary surface an impedance boundary conditions are applied so that no wave reflections occur. In the proposed scheme, all of the procedures are processed directly in time domain, which is efficient for nonlinear analyses of structure floating on unbounded fluid. Numerical results indicate acceptable accuracy of the proposed method.     

波浪沿斜面传播过程的数值模拟

精确模拟非线性波沿斜面传播过程非常困难,为此论文从势函数的边界积分方程出发,建立了一种时域内二维波浪模拟的数值模型,主要用来模拟完全非线性波浪的传播变形过程。论文的数值模型使用高阶二维边界元方法,采用可调节时间步长的基于二阶显式泰勒展开的混合欧拉-拉格郎日时间步进来求解带自由表面的线性或完全非线性波浪传播问题。在计算区域一端造出线性或非线性的周期性波浪,另一端采用消除反射波的人工粘性吸收边界。通过与现有理论比较证明了论文数值方法所得结果是准确可靠的。     

Short-crested waves interaction with a concentric porous cylinder system with partially porous outer cylinder

In this paper, based on the linear wave theory, the interaction of short-crested waves with a concentric dual cylindrical system with a partially porous outer cylinder is studied by using the scaled boundary finite element method (SBFEM), which is a novel semi-analytical method with the advantages of combining the finite element method (FEM) with the boundary element method (BEM). The whole solution domain is divided into one unbounded sub-domain and one bounded sub-domain by the exterior cylinder. By weakening the governing differential equation in the circumferential direction, the SBFEM equations for both domains can be solved analytically in the radial direction. Only the boundary on the circumference of the exterior porous cylinder is discretized with curved surface finite elements. Meanwhile, by introducing a variable porous-effect parameter G, non-homogeneous materials caused by the complex configuration of the exterior cylinder are modeled without additional efforts. Comparisons clearly demonstrate the excellent accuracy and computational efficiency associated with the present SBFEM. The effects of the wide range wave parameters and the structure configuration are examined. This parametric study will help determine the various hydrodynamic effects of the concentric porous cylindrical structure.     

波浪作用下海床动力反应的数值分析

近海和离岸建筑物的基通常处于连续不断的小风浪作用之下，可将其变形视为弹性。基于二维广义Biot理论，提出了线性或非线性波浪作用下饱和弹性海床动力应应的时域有限元数值解法。静力平衡条件和Biot方程组成的边值方程可视为其特例，在比较算例中，数值计算得到的孔压和有效应力幅值沿海床深度的分布与解析解十分吻合。土骨架和孔隙流体的加速度对海床动力反应的影响很小。具体算例表明，线性波沿缓坡海床传播时，土层中超静孔压和有效应力幅值随之增大，有可能发生滑动坡坏。所提出的数值解法能够灵活地处理非线性波浪荷载，海底复杂地质条件和波浪沿缓坡传播等复杂情况。     

Three-dimensional numerical simulation of oblique wave action on vertical walls

1 Introduction A vertical wall is one of the typical maritime sheltering structures in the coastal region, which is mainly subjected to wave forces. It has been a com- mon engineering assumption adopted for design pur- poses that normally incident wave fo…     

Time-domain simulation of berthing problem between FPSO and shuttle tanker in waves

A time-domain simulation method based on potential flow model has been developed to investigate the berthing problem between two floating bodies in wave. The boundary value problem is formulated with respect to an earth-fixed coordinate system because the relative positions of the two vessels continuously change during the berthing operation. The classical finite element method is used to solve the Laplace equation in the fluid domain with moving boundary. The linearized free-surface boundary conditions are integrated in time by applying 4th-order Adams–Bashforth–Moulton method. A simple re-mesh algorithm with local and global mesh systems is introduced to update mesh by considering large horizontal movement of the berthing vessel. The developed numerical method is used to investigate the berthing problem between a FPSO and shuttle tanker in waves. The focus is on the wave-induced motion response during the berthing process. The characteristics of the motion responses in berthing operation are examined with various wave frequencies, berthing speeds and wave headings.     

近岸波浪折射-绕射-破波耗散联合模式的有限元数值研究

建立了近岸波浪折射-绕射-破波耗散的有限元数值模式。采用的有限元方法为改进的混合元法,其中外域开边界条件得到改进,内域有限元为伽廖金有限元。用理论解检验了所建立的数值模式,并将该模式应用到一个模型港湾。