Fast computation of the transient motions of moving vessels in irregular ocean waves

A fast time-domain method is developed in this paper for the real-time prediction of the six degree of freedom motions of a vessel traveling in an irregular seaway in infinitely deep water. The fully coupled unsteady ship motion problem is solved by time-stepping the linearized boundary conditions on both the free surface and body surface. A velocity-based boundary integral method is then used to solve the Laplace equation at every time step for the fluid kinematics, while a scalar integral equation is solved for the total fluid pressure. The boundary integral equations are applied to both the physical fluid domain outside the body and a fictitious fluid region inside the body, enabling use of the fast Fourier transform method to evaluate the free surface integrals. The computational efficiency of the scheme is further improved through use of the method of images to eliminate source singularities on the free surface while retaining vortex/dipole singularities that decay more rapidly in space. The resulting numerical algorithm runs 2–3 times faster than real time on a standard desktop computer. Numerical predictions are compared to prior published results for the transient motions of a hemisphere and laboratory measurements of the motions of a free running vessel in oblique waves with good agreement.     

Nonlinear response of membranes to ocean waves using boundary and finite elements

The behavior of a highly deformable membrane to ocean waves was studied by coupling a nonlinear boundary element model of the fluid domain to a nonlinear finite element model of the membrane. The hydrodynamic loadings induced by water waves are computed assuming large body hydrodynamics and ideal fluid flow and then solving the transient diffraction/radiation problem. Either linear waves or finite amplitude waves can be assumed in the model and thus the nonlinear kinematic and dynamic free surface boundary conditions are solved iteratively. The nonlinear nature of the boundary condition requires a time domain solution. To implicitly include time in the governing field equation, Volterra's method was used. The approach is the same as the typical boundary element method for a fluid domain where the governing field equation is the starting point. The difference is that in Volterra's method the time derivative of the governing field equation becomes the starting point.The boundary element model was then coupled through an iterative process to a finite element model of membrane structures. The coupled model predicts the nonlinear interaction of nonlinear water waves with highly deformable bodies. To verify the coupled model a large scale test was conducted in the OH Hinsdale wave Research Laboratory at Oregon State University on a 3-ft-diameter fabric cylinder submerged in the wave tank. The model data verified the numerical prediction of the structure displacements and of the changes in the wave field.The boundary element model is an ideal modeling technique for modeling the fluid domain when the governing field equations is the Laplace equation. In this case the nonlinear boundary element model was coupled with a finite element model of membrane structures, but the model could have been coupled with other finite element models of more rigid structures, such as a pontoon floating breakwater.     

非线性波浪波面追踪的一种新模式

基于Laplace方程的Green积分表达式和波面BemouUi方程所建立的非线性波动数学模型,是一个时域上具有初始值的边值问题,而精确地追踪自由表面的波动位置,给出波面运动瞬时的波面高度和波面势函数,是建立时域内非线性波浪数值模式的基础。本文采用0-1混合型边界元剖分计算域边界并离散Laplace方程的Green积分表达式,采用有限元剖分自由水面并推导满足自由表面非线性边界条件的波面有限元方程,联立计算域内以节点波势函数和波面位置高度的时间增量为未知量的线性方程组,通过时步内的循环迭代,给出每个时步上的波面位置和波面势函数,从而建立了一种新的非线性波浪波面追踪模式。数值造波水槽内的波浪试验表明,其数值模拟结果具有良好的计算精度。     

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.     

Irregular points in 2-D free surface flows with surface tension for the wavemaker boundary value problem

Oscillatory, two dimensional free-surface, irrotational flows of an incompressible fluid are treated numerically. Rather than basing the treatment on perturbation expansions, the unknown fluid domain bounded by a free surface is mapped onto a rectangular domain. This has the effect of exchanging a boundary value problem for a linear, constant coefficient Laplace equation having nonlinear boundary conditions into a nonlinear elliptic equation specified on the known rectangular domain. In the known transformed domain, compatibility conditions at the corners are derived, the boundary value problem is solved and the results are transformed back to the original domain.     

Sloshing waves in electrified fluid layers in an excited rectangular container

This work aims to demonstrate the dynamical properties of the free interface separating two bounded layers of electrified Newtonian fluids inside a parametrically excited boxed basin. The mathematical model of the investigated problem is contacted with the linearized Navier–Stokes equation of viscous fluids and Maxwell equations together with the boundary conditions that are solved utilizing Laplace transform. Solutions of the governing equations in time domain are then numerically calculated by using Durbin's numerical inverse Laplace transform scheme. Graphical illustrations of the accomplished numerical results are provided to examine the influence of some selected parameters on the stability picture of the interfacial waves as well as the electric surface charge distribution. The numerical results proposed that the electric Euler number Eu, promotes the instability of the wave motion whatever the values of the frequency of oscillation and the destabilizing effect of Eu becomes weaker at larger values of the frequency; while a dual effect (in) stabilizing of the frequency of oscillations is reported whether the electric effects exist or not. This irregular influence of the frequency oscillation has also been monitored regarding to the forces exerted on the side wall of the container.     

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

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

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.     

Study on the Interaction of Water Waves with Semi-Circular Breakwater

—The present study investigates the interaction of steep waves with semi-circular breakwaterwith the complex plane＇s Cauchy boundary integral theorem.The boundary integral method is used totransform the calculation in fluid domain into its boundary alone.In the calculation the computation do-main is moved with the propagation of waves.A numerical solution is obtained for incident Stokes wavespassing the submerged obstacles.This method has been extended to the calculation of wave run-up on aslope for estimating wave overtopping.     

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

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

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

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

Time domain three-dimensional fully nonlinear computations of steady body–wave interaction problem

Three-dimensional fully nonlinear waves generated by moving disturbances with steady forward speed without motions are solved using a mixed Eulerian–Lagrangian method in terms of an indirect boundary integral method and a Runge–Kutta time marching approach which integrates the fully nonlinear free surface boundary conditions with respect to time.A moving computational window is used in the computations by truncating the fluid domain (the free surface) into a computational domain. The computational window maintains the computational domain and tracks the free surface profile by a node-shifting scheme applied within it. An implicit implement of far field condition is enforced automatically at the truncation boundary of the computational window.Numerical computations are applied to free surface waves generated by Wigley and Series 60 hulls for the steady problem. The present numerical results are presented and compared with existing linear theory, experimental measurements, and other numerical nonlinear computations. The comparisons show satisfactory agreements for these hydrodynamic problems.     

Time domain simulation of side-by-side floating bodies using a 3D numerical wave tank approach

The coupled system of two side-by-side fixed and/or floating bodies interacting with a large amplitude nonlinear wave is studied using a direct time domain solution method. The numerical model is based on a three-dimensional mixed Eulerian–Lagrangian (MEL) method under certain simplifying approximations permitting Rankine panel scheme to be implemented over a time-invariant boundary surface to solve the boundary value problem for the unknown velocity potentials. A 4th order Adams–Bashforth–Moulton scheme is used for time marching of rigid-body motion histories of the individual bodies and evolution of the free-surface including the gap region in which large resonant fluid motions occur. A systematic study has been carried out to evaluate the performance of the developed time domain method in simulating the forces and motions as well as the fluid motion in the gap region for the two body system under various arrangements and in different wave-headings. At first, the computed numerical results have been validated and verified with computational and experimental results available in literature for standard geometries such as vertical truncated cylinders and rectangular boxes. Secondly, effectiveness of the damping lid model which is introduced to suppress wave resonance in the gap region is investigated including its influence on maximum sway forces on fixed and floating rectangular barges in side-by-side configurations. Thirdly, comparative studies on absolute and relative motion response for two cases (two rectangular barges, and a FLNG-FPSO + shuttle tanker) in side-by-side arrangement are detailed to bring out the importance of nonlinearities arising due to steep nonlinear incident waves. Finally, coupled motions of the two-body system of an FPSO and a shuttle tanker floating in side-by-side configuration in a steep nonlinear wave field are studied in which the two bodies are connected through hawsers, and also the FPSO is moored to the ground. Additionally there is a fender between the two bodies.     

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.     

Numerical simulation of unsteady free-surface motions using a boundary integral formulation

A numerical time simulation method is described to solve fluid flow problems including unsteady free surface motion. The method is based on potential flow theory. At every time step, the problem is solved using a boundary integral formulation of the fluid domain. The linearized free surface conditions are integrated in time and the solution is marched forward. Computational results simulating the free surface motion for the cases of a linear progressive wave, wave propagating into a region of calm water and the wave maker motion are presented. Comparison with theoretical results demonstrate the feasibility of the proposed simulation scheme.     

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.     

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

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

New Derivation of Ordinary Differential Equations for Transient Free-Surface Green Functions

Based on the Laplace transform, a direct derivation of the ordinary differential equations for the three-dimensional transient free-surface Green function in marine hydrodynamics is presented. The results for the 3D Green function and all its spatial derivatives are a set of fourth-order ordinary differential equations, which are identical with that of Clement (1998). All of these results may be used to accelerate numerical computation for the time-domain boundary element method in marine hydrodynamics.