A NURBS-based high-order panel method for three-dimensional radiation and diffraction problems with forward speed

A high-order Rankine panel method based on Non-Uniform Rational B-Spline (NURBS) is developed for solving the three-dimensional radiation and diffraction problems with forward speed. A NURBS surface is used to precisely represent the body geometry. Velocity potential on the body surface is described by B-spline after the source density distribution on the boundary surface is determined. A collocation approach is applied to numerical computation and the integral equations are evaluated by applying Gauss–Legendre quadrature. The m_j-terms are evaluated by a desingularized method which utilizes NURBS technique. In order to verify the method proposed, it is firstly applied to the unbounded flow problem of a sphere and spheroids. The numerical results are found to be in good agreement with analytical solutions. Then the method is used to solve the radiation and diffraction problems of a sphere and the diffraction problem of a spheroid moving with a forward speed beneath the free surface in frequency domain. The numerical results are satisfactory in comparison with the published analytical results and experimental results.     

Performance analysis of 3D hydrofoil under free surface

The purpose of the present paper is to develop a potential-based panel method for determining the steady potential flow about three-dimensional hydrofoil under free surface. The method uses constant-strength doublets and source density distribution over the foil body surface and thereby Dirichlet-type boundary condition is used instead of Neumann-type condition. On the undisturbed free surface source density is used to meet the free surface condition that is linearised in terms of double-body model approach and is discretised by a one-side, upstream, four-point finite difference operator. After solving the doublets on the foil and sources on the free surface, the numerical results of pressure, lift and resistance coefficients and also wave profiles can then be calculated for different Froude number and depth of submergence to demonstrate the influence of free surface and aspect ratio effects on performance of the hydrofoil.     

Time-domain simulation of wave–structure interaction based on multi-transmitting formula coupled with damping zone method for radiation boundary condition

Based on the Rankine source, this paper proposed a time-domain method for analyzing the three-dimensional wave–structure interaction problem in irregular wave. A stable integral form of the free-surface boundary condition (IFBC) is employed to update the velocity potential on the free surface. A multi-transmitting formula, with an artificial wave speed, is used to eliminate the wave reflection for radiation condition on the artificial boundary. An effective multi-transmitting formula, coupled with damping zone method, is further used to analyze the irregular wave diffraction at the artificial boundary. We investigate hydrodynamic forces on floating structure and compare our solution to the frequency-domain solution. It is shown that long time simulation can be done with high stability and the numerical results agree well with the solution obtained under the frequency domain. The efficiency of the proposed multi-transmitting formula and the coupled methods for radiation boundary make them promising candidates in studying the irregular water wave problem in time domain.     

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

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

Numerical and experimental study on seakeeping performance of ship in finite water depth

Vessels operating in shallow waters require careful observation of the finite-depth effect. In present study, a Rankine source method that includes the shallow water effect and double body steady flow effect is developed in frequency domain. In order to verify present numerical methods, two experiments were carried out respectively to measure the wave loads and free motions for ship advancing with forward speed in head regular waves. Numerical results are systematically compared with experiments and other solutions using the double body basis flow approach, the Neumann-Kelvin approach with simplified m-terms, and linearized free surface boundary conditions with double-body m-terms. Furthermore, the influence of water depths on added mass and damping coefficients, wave excitation forces, motions and unsteady wave patterns are deeply investigated. It is found that finite-depth effect is important and unsteady wave pattern in shallow water is dependent on both of the Brard number τ and depth Froude number F_h.     

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

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

On pressure forcing techniques for numerical modelling of unidirectional nonlinear water waves

A method of incorporating pressure forcing into a nonlinear potential flow wave model is presented. A semi-analytical pseudo-spectral method is used to calculate dynamic response of a water body exposed to evolving local pressure distribution. Surface slope coherent and slope proportional pressure functions are directly applied through a pressure term appearing in the dynamic free-surface boundary condition of a formulated initial boundary-value problem. First, a monochromatic pressure distribution is used to generate steady regular waves of permanent form. The pressure-induced wave motion exhibits stable harmonic structure for deepwater, transitional water and shallow water waves. In the next step, a more complex pressure system is used to initiate multi-component wave propagation. It is demonstrated that the proposed method provides well-posed initial conditions for studying various water wave scenarios within a framework of nonlinear potential flow solutions.     

High-speed submerged and surface piercing cavitating hydrofoils, including tandem case

The iterative method that is originally developed before both for two- and three-dimensional single cavitating hydrofoils moving with a constant speed under free surface is applied to the case of high-speed (Froude number up to 6.5) and some figures are given. The method is also extended to include the surface piercing hydrofoils (vertical struts) and the case of tandem hydrofoils into the calculations. The iterative nonlinear method based on the Green's theorem allows separating the cavitating hydrofoil problem(s) and the free surface problem. These two (or three in the case of tandem hydrofoil) problems are solved separately, with the effects of one on the other being accounted for in an iterative manner. The cavitating hydrofoil surface(s) and the free surface are modeled with constant strength dipole and constant strength source panels. The source strengths on the free surface are expressed in terms of perturbation potential by applying the linearized free surface conditions. No radiation condition is enforced for downstream and transverse boundaries. The cavitation number is expressed in terms of Froude number and the submergence depth of the hydrofoil from the free surface. An algebraic grid on the free surface has been described to get a smooth transition between the panels along the direction of uniform inflow and to have a long distance in the downstream direction depending on the wave-length (or Froude number) while keeping the number of panels fixed. First, the method is validated in the case of surface piercing hydrofoil. Then, the effects of high Froude number and the submergence depth of the hydrofoil from free surface on the results are discussed and some figures are given for interested engineers and designers. The method is later applied to the case of tandem hydrofoils and the effects of one hydrofoil on the other are discussed.     

Wave radiation and diffraction by a two-dimensional floating body with an opening near a side wall

The radiation and diffraction problem of a two-dimensional rectangular body with an opening floating on a semi-infinite fluid domain of finite water depth is analysed based on the linearized velocity potential theory through an analytical solution procedure. The expressions for potentials are obtained by the method of variation separation, in which the unknown coefficients are determined by the boundary condition and matching requirement on the interface. The effects of the position of the hole and the gap between the body and side wall on hydrodynamic characteristics are investigated. Some resonance is observed like piston motion in a moon pool and sloshing in a closed tank because of the existence of restricted fluid domains.     

A potential based panel method for 2-D hydrofoils

A potential based panel method for the hydrodynamic analysis of 2-D hydrofoils moving beneath the free surface with constant speed without considering cavitation is described. By applying Green's theorem and the Green function method, an integral equation for the perturbation velocity potential is obtained under the potential flow theory. Dirichlet type boundary condition is used instead of Neumann type boundary condition. The 2-D hydrofoil is approximated by line panels which have constant source strength and constant doublet strength distributions. The free surface condition is linearized and the method of images is used for satisfying this free surface condition. All the terms in fundamental solution (Green function) of perturbation potential are integrated over a line panel. Pressure distribution, lift, residual drag and free surface deformations are calculated for NACA4412, symmetric Joukowski and van de Vooren profile types of hydrofoil. The results of this method show good agreement with both experimental and numerical methods in the literature for the NACA4412 and symmetric Joukowski profile types. The lift and residual drag values of the van de Vooren profile are also presented. The effect of free surface is examined by a parametric variation of Froude number and depth of submergence.     

A quasi-cloaking phenomenon to reduce the wave drift force on an array of adjacent floating bodies

The mooring of offshore floating structures, such as offshore platforms, in large waves against drift forces and rotational moments is a challenging problem in offshore engineering. To accurately investigate such problems, called positioning problems, the time-averaged steady forces of the second order known as the wave drift forces must be taken into account. Fortunately, a cloaking phenomenon occurs under certain conditions and dramatically reduces the wave drift force acting on such a floating body, as previously reported by several researchers. In the diffraction problem of water waves, cloaking refers to the condition where there is no scattering in the form of radial outgoing waves. The reduction of wave drift force on a truncated cylinder with the occurrence of cloaking phenomenon has been numerically and experimentally confirmed. In this paper, the arrangement of several small circular cylinders at regular intervals in a circle concentric with a fixed floating body is considered as an effective means of reducing the wave drift force. Using a combination of a higher-order boundary element method (HOBEM) and wave interaction theory, the influences of the geometric parameters of the outer surrounding cylinders on the wave drift force and the total scattered-wave energy are systematically investigated and discussed. A quasi-cloaking phenomenon is first found and reported in the present study, which is beneficial and flexible for application in practical engineering. More than one quasi-cloaking trigger (where a trigger is an occurrence condition) can be found simply by varying the distance between the inner and outer floating bodies.     

Springing analysis of a seagoing vessel using fully coupled BEM-FEM in the time domain

The quasi-steady resonant vibration of a flexible seagoing vessel under resonant wave excitation force, called springing, is studied in this paper. A higher-order B-spline Rankine panel method is used to represent the effects of the fluid motion surrounding this flexible seagoing vessel, and a finite element formulation based on Vlasov beam is employed for structural response. The boundary integral equation and finite element equation, both for fluid and structural domains, are fully coupled with each other using an iterative implicit method in the time domain. Coupling between the two field equations is achieved by relying on fixed-point iteration with relaxation aided by Aitken's δ² process to maximize convergence speed. The steady-unsteady coupling term or m-term in the linearized body boundary condition derived by Timman and Newman is taken into account for accurate prediction of flexible body motion when forward speed is present. The 2nd derivative of basis potential in the m-term is obtained by modifying Nakos approach, which was originally developed using the Stokes theorem for rigid body ship motion problem. For the solution of the FE equation, instead of conventionally used modal superposition method, a direct integration scheme based on Newmark method is employed. It is believed that this technique is more attractive in the sense that it allows us free from the selection of optimum number of mode-shapes in the computation.     

Stability analysis of high-order finite-difference discretizations of the linearized forward-speed seakeeping problem

A high-order finite-difference method solution of the linearized, potential flow, seakeeping problem for a ship at steady forward speed was recently presented by Amini-Afshar et al. [1,2]. In this paper, we provide a detailed matrix-based eigenvalue stability analysis of this model, highlighting the sources of instability and the effects of possible remedies. In particular, we illustrate how both boundary treatment and grid stretching are important factors which are not typically captured by a von Neumann-type analysis. The new analysis shows that when grid stretching is used together with centered finite difference schemes, the method is generally unstable. The source of the instability can in some cases be traced to an effective downwinding of the convective terms. Stable solutions can be obtained either by introducing upwind-biased schemes for computing the convective derivatives on the free-surface, or by application of a mild filter at each time-step. A second source of instability is associated with the treatment of the convective derivatives of the free-surface elevation at points close to the domain boundaries. Here it is necessary to consider whether the surrounding fluid points lie in an upwind or a downwind direction. For upwinded points, ordinary one-sided differencing can be used, but for downwinded points we instead impose a Neumann-type boundary condition derived from the body and free-surface boundary conditions. As an example application to complement those already given in [1], [2], the method is applied to solve the steady wave resistance problem and comparison is made to reference solutions for a two-dimensional floating cylinder and a submerged sphere. Estimates of the wave resistance of the Wigley hull are also compared with experimental measurements.     

高阶边界元方法求解规则波在三维局部渗透海床上的传播

基于线性势流理论,利用高阶边界元法研究了规则波在三维局部渗透海床上的传播。根据Darcy渗透定律推导出渗透海床的控制方程,利用渗透海床顶部和海底处法向速度和压强连续条件得到渗透海床顶部满足的边界条件。根据绕射理论,利用满足自由水面条件的格林函数建立了求解渗透海床绕射势的边界积分方程,采用高阶边界元方法求解边界积分方程进而得到自由水面的绕射势和波浪在局部渗透海床上传播过程中幅值的变化情况。通过与已发表的波浪对圆柱形暗礁的时域全绕射结果对比,证明了本文建立的频域方法计算波幅的正确性和有效性。利用这一模型研究了三维矩形渗透海床区域上波浪的传播特性,并分析了入射波波长、海床渗透特性系数等参数对波浪传播的影响。     

A potential based panel method for 2-D hydrofoils

A potential based panel method for the hydrodynamic analysis of 2-D hydrofoils moving beneath the free surface with constant speed without considering cavitation is described. By applying Green's theorem and the Green function method, an integral equation for the perturbation velocity potential is obtained under the potential flow theory. Dirichlet type boundary condition is used instead of Neumann type boundary condition. The 2-D hydrofoil is approximated by line panels which have constant source strength and constant doublet strength distributions. The free surface condition is linearized and the method of images is used for satisfying this free surface condition. All the terms in fundamental solution (Green function) of perturbation potential are integrated over a line panel. Pressure distribution, lift, residual drag and free surface deformations are calculated for NACA4412, symmetric Joukowski and van de Vooren profile types of hydrofoil. The results of this method show good agreement with both experimental and numerical methods in the literature for the NACA4412 and symmetric Joukowski profile types. The lift and residual drag values of the van de Vooren profile are also presented. The effect of free surface is examined by a parametric variation of Froude number and depth of submergence.     

A panel-free method for time-domain analysis of the radiation problem

A panel-free method (PFM), based on the desingularized Green’s formulae proposed by Landweber and Macagno, has been developed to solve the radiation problem of a floating body in the time domain. The velocity potential due to a non-impulsive velocity is obtained by solving the boundary integral equation in terms of source strength distribution. The singularity in the Rankine source term of the time-dependent Green function is removed. The geometry of a body surface is mathematically represented by NURBS surfaces. The integral equation can be globally discretized over the body surface by Gaussian quadratures. No assumption is needed for certain degree of approximation of distributed source strength on the body surface. The accuracy of PFM was demonstrated by its application to a classical problem of uniform flow past a sphere. The response function of a hemisphere at zero speed was then computed by PFM. The computed response function, added-mass and damping coefficients are compared with other published results.     

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.     

Fully Nonlinear Time Domain Analysis for Hydrodynamic Performance of An Oscillating Wave Surge Converter

The hydrodynamic behaviour of an oscillating wave surge converter (OWSC) in large motion excited by nonlinear waves is investigated. The mechanism through which the wave energy is absorbed in the nonlinear system is analysed. The mathematical model used is based on the velocity potential theory together with the fully nonlinear boundary conditions on the moving body surface and deforming free surface. The problem is solved by the boundary element method. Numerical results are obtained to show how to adjust the mechanical properties of the OWSC to achieve the best efficiency in a given wave, together with the nonlinear effect of the wave height. Numerical results are also provided to show the behaviour of a given OWSC in waves of different frequencies and different heights.     

海洋中波浪下迎浪驳船运动对两矩形驳船间窄缝内水体共振的影响

A three-dimensional time-domain potential flow model is developed and applied to simulate the wave resonance in a gap between two side-by-side rectangular barges. A fourth-order predict-correct method is implemented to update free surface boundary conditions. The response of an up-wave barge is predicted by solving the motion equation with the Newmark-β method. Following the validation of the developed numerical model for wave radiation and diffraction around two side-by-side barges, the influence of up-wave barge motion on the gap surfaceresonance is investigated in two different locations of the up-wave barge relative to the back-wave barge at various frequencies. The results reveal that the freely floating up-wave barge significantly influences the resonance frequency and the resonance wave amplitude. Simultaneously, the up-wave barge located in the middle of the back-wave barge leads to a reduction in the resonance wave amplitude and motion response when compared with other configurations.     

有航速浮体时域分析的混合边界元法

给出了一种联合瞬态格林函数和Rankine源进行有航速浮体时域水动力分析的混合——边界元方法。在三维混合边界元方法中,通过一个匹配面将流体域划分为内域和外域,在内域中使用Rankine源以模拟直壁或非直壁船体及线性或非线性自由面条件,在外域中使用瞬态格林函数以满足自由面条件和远方辐射条件。使用该方法计算了一个有航速潜没圆球的波浪力,和解析结果的比较证明了该方法的正确性。进一步给出了一个有航速Wigley船的水动力结果,计算结果稳定,没有外传波向内反射的现象发生。