A new ordinary differential equation for the evaluation of the frequency-domain Green function

Clément (2013) derived a second order ordinary differential equation (ODE) satisfied by the free-surface Green function in the frequency domain. Since then, similar ODEs for the gradient of the Green function have been developed. Unfortunately, all these ODEs degenerate at zero frequency. Therefore, it is not possible to initialize the numerical solution of these ODEs from this zero frequency. Alternative methods based on the shifting of the initial condition to frequencies strictly greater than zero have then been developed.The present paper describes an alternative approach to address this issue. It involves a new function which is the solution of a modified ODE which can be solved from the zero frequency.Finally, comparisons with evaluations of the Green function using the classical direct integration method are provided. They show that the new ODE can provide accurate estimates of the Green function.     

On the evaluation of time-domain Green function

An analytical method has been developed to evaluate the wave part of the time-domain Green function and its derivatives. Based on Taylor series expansion, the Green function is obtained by solving a fourth-order ordinary differential equation. The method accelerates the convergence of the summation of an infinite series in the numerical computation. The accuracy of this method was demonstrated by its comparison with other method and its application to solve the radiation problem of a floating hemisphere using a panel-free method. The computed hydrodynamic coefficients agree well with the analytical solutions.     

Extension of the Frequency-Domain pFFT Method for Wave Structure Interaction in Finite Depth

To analyze wave interaction with a large scale body in the frequency domain, a precorrected Fast Fourier Transform (pFFT) method has been proposed for infinite depth problems with the deep water Green function, as it can form a matrix with Toeplitz and Hankel properties. In this paper, a method is proposed to decompose the finite depth Green function into two terms, which can form matrices with the Toeplitz and a Hankel properties respectively. Then, a pFFT method for finite depth problems is developed. Based on the pFFT method, a numerical code pFFT-HOBEM is developed with the discretization of high order elements. The model is validated, and examinations on the computing efficiency and memory requirement of the new method have also been carried out. It shows that the new method has the same advantages as that for infinite depth.     

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.     

Numerical simulation of non-linear regular and focused waves in an infinite water-depth

Inviscid three-dimensional free surface wave motions are simulated using a novel quadratic higher order boundary element model (HOBEM) based on potential theory for irrotational, incompressible fluid flow in an infinite water-depth. The free surface boundary conditions are fully non-linear. Based on the use of images, a channel Green function is developed and applied to the present model so that two lateral surfaces of an infinite-depth wave tank can be excluded from the calculation domain. In order to generate incident waves and dissipate outgoing waves, a non-reflective wave generator, composed of a series of vertically aligned point sources in the computational domain, is used in conjunction with upstream and downstream damping layers. Numerical experiments are carried out, with linear and fully non-linear, regular and focused waves. It can be seen from the results that the present approach is effective in generating a specified wave profile in an infinite water-depth without reflection at the open boundaries, and fully non-linear numerical simulations compare well with theoretical solutions. The present numerical technique is aimed at efficient modelling of the non-linear wave interactions with ocean structures in deep water.     

Interaction of uniform current with a circular cylinder submerged below an ice sheet

The problem of a uniform current passing through a circular cylinder submerged below an ice sheet is considered. The fluid flow is described by the linearized velocity potential theory, while the ice sheet is modelled through a thin elastic plate floating on the water surface. The Green function due to a source is first derived, which satisfies all the boundary conditions apart from that on the body surface. Through differentiating the Green function with respect to the source position, the multipoles are obtained. This allows the disturbed velocity potential to be constructed in the form of an infinite series with unknown coefficients which are obtained from the boundary condition. The result shows that there is a critical Froude number which depends on the physical properties of the ice sheet. Below this number there will be no flexural waves propagating to infinity and above this number there will be two waves, one on each side of the body. When the depth based Froude number is larger than 1, there will always be a wave at far upstream of the body. This is similar to those noticed in the related problem and is different from that in the free surface problem without ice sheet. Various results are provided, including the properties of the dispersion equation, resistance and lift, ice sheet deflection, and their physical features are discussed.     

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...     

On the radiation and diffraction of linear water waves by an infinitely long rectangular structure submerged in oblique seas

The radiation and diffraction of linear water waves by an infinitely long rectangular structure submerged in oblique seas of finite depth is investigated. The analytical expressions for the radiated and diffracted potentials are derived as infinite series by use of the method of separation of variables. The unknown coefficients in the series are determined by the eigenfunction expansion matching method. The expressions for wave forces, hydrodynamic coefficients and reflection and transmission coefficients are given and verified by the boundary element method. Using the present analytical solution, the hydrodynamic influences of the angle of incidence, the submergence, the width and the thickness of the structure on the wave forces, hydrodynamic coefficients, and reflection and transmission coefficients are discussed in detail.     

有限水深中零航速浮体波浪力计算奇点积分的数值处理方法

使用三维源汇分布法Ⅲ计算有限水深中零航速浮体所受到的波浪力,对两种不同形式的格林函数中所共同存在的奇点问题分别进行了处理。公式推导表明,使用级数形式的格林函数可以使计算更加快捷。最后,对不同尺度的圆柱体进行了验算,在对计算结果与解析解进行了比较之后,工程计算也证明选择格林函数级数计算公式是更令人满意的方案。     

Study on characteristics of 3-D translating-pulsating source green function of deep-water havelock form and its fast integration method

The singularities, oscillatory performances and the contributing factors to the 3-D translating-pulsating source Green function of deep-water Havelock form which consists of a local disturbance part and a far-field wave-like part, are analyzed systematically. Relative numerical integral methods about the two parts are presented in this paper. An improved method based on LOBATTO rule is used to eliminate singularities caused respectively by infinite discontinuity and jump discontinuous node from the local disturbance part function, which makes the improvement of calculation efficiency and accuracy possible. And variable substitution is applied to remove the singularity existing at the end of the integral interval of the far-field wave-like part function. Two auxiliary techniques such as valid interval calculation and local refinement of integral steps technique in narrow zones near false singularities are applied so as to avoid unnecessary integration of invalid interval and improve integral accordance. Numerical test results have proved the efficiency and accuracy in these integral methods that thus can be applied to calculate hydrodynamic performance of floating structures moving in waves.     

Radiation and diffraction of linear water waves by an infinitely long submerged rectangular structure parallel to a vertical wall

The radiation and the diffraction of linear water waves by an infinitely long floating rectangular structure submerged in water of finite depth with leeward boundary being a vertical wall are analyzed in this paper by using the method of separation of variables. Analytical expressions for the radiated and diffracted potentials are derived as infinite series with unknown coefficients determined by the eigenfunction expansion matching method. The expressions for wave forces and hydrodynamic coefficients are given. A comparison is made between the results obtained by the present analytical solution and those obtained by the boundary element method. By using the present analytical solution, the hydrodynamic influences of the submergence, the width, the thickness of the structure, and the distance between the structure and the wall on the wave forces and hydrodynamic coefficients are discussed in detail.     

无限水深聚焦波完全非线性数值模拟

基于势流理论提出一种新的高阶边界元方法对无限水深的聚焦波浪进行完全非线性数值模拟.自由水面满足完全非线性边界条件,模拟波浪的非线性效果可以达到更高阶.利用镜像原理,建立一种全新的格林函数应用到无限水深的数值波浪水槽中,以致于两无限深水槽侧壁的积分可以被排除.为了产生相应的入射波和吸收出流波浪,一个由点源组成的造波装置被布置于计算域内,同时人工阻尼层被用来吸引出流波浪,由波浪聚焦的方法得到极限波浪.通过开展线性和完全非线性聚焦波浪的数值实验及与理论解对比,验证本数值模型可以用来模拟无限深水域的极限波浪,且在出流边界没有反射.     

On the side wall effects upon bodies of arbitrary geometry in wave tanks

A numerical model based on using a tank Green function, has been developed to compute the side wall effects on first- and second-order loadings upon bodies of arbitrary geometry in wave tanks. This tank Green function (TGF) is composed of a finite series of open-sea Green functions and an asymptotic part represented by two single integrals whose kernels decrease exponentially with the integral variable. This consistent expression of the TGF permits one to highlight the side wall effects and to give some criteria for the choice of tank width and the measurement duration to limit the reflection of diffraction and radiation waves.

The efficiency of the developed model is shown in the application to hemispheres and a box-shaped barge placed in the center of the wave tanks. The numerical results explain well the irregularities in the experimental measurements and show that the side walls have important effects on the first-order quantities. These effects are much more pronounced on the second-order drift loads.     

Analysis of Wave Loads on A Semi-Submersible Platform

For the global and structural fatigue strength analysis of a semi-submersible platform, wave loads under design conditions are calculated by use of the three-dimensional boundary element method. Methods for calculating the forward-speed free-surface Green function are discussed and a computer program with this Green function is developed. According to the special rules, the wave loads under several typical design conditions of the platform are calculated. The maximum vertical bending moment, torsion moment and horizontal split force are determined from a series of contour maps of wave loads for the wave period of 5 to 18 seconds at a certain interval and the wave phase of 0°to 360°at a certain interval. The wave height is determined by the function of wave period with a given exceedance probability. The maximum wave loads under the combination of wave parameters are used as the input of hydrodynamic pressure in the three-dimensional finite element analysis process. The transfer functions of wave loads     

An improved model for the dielectric constant of sea water at microwave frequencies

The advent of precision microwave radiometry has placed a stringent requirement on the accuracy with which the dielectric constant of sea water must be known. To this end, measurements of the dielectric constant have been conducted atS-band andL-band with a quoted uncertainty of tenths of a percent. These and earlier results are critically examined, and expressions are developed which will yield computations of brightness temperature having an error of no more than 0.3 K for an undisturbed sea at frequencies lower thanX-band. At the higher microwave and millimeter wave frequencies, the accuracy is in question because of uncertainties in the relaxation time and the dielectric constant at infinite frequency.     

预修正快速傅里叶变换高阶边界元方法在波浪对物体作用问题中的应用

边界元方法被广泛应用于波浪对海上婕筑物作用领域,但由于传统边界元方法的存储量和计算量均为未知量的平方量级,很难满足大范围多未知量问题的计算需要.本文基于高阶边界元方法,应用预修正快速傅里叶变换方法,使计算量与存储量分别降低至O(NlogN)和O(N)量级,并得到一个连续的压强分布以适应结构设计的要求,同时可以通过使用满...     

Hydrodynamic coefficients of a two-dimensional,truncated rectangular floating structure in shallow water

The paper seeks to examine hydrodynamic coefficients of a rectangular structure in shallow water and to establish analytical formulae for fast computations. A two-dimensional rectangular profile is considered with the under-bottom clearance assumed to be small compared with structure dimensions and the water depth. Following the method of matched asymptotic expansions, the radiation problem is solved under assumptions of the linear wave theory, by matching two ‘outer’ flows with the ‘inner’ flow near the structure edge. Closed asymptotic formulae are obtained for all hydrodynamic coefficients for heave, sway and roll motions. The zero and infinite frequency values of the added mass are examined and formulae are derived intended for quick engineering estimations. Numerical results compare well with those published in literature, and the approach is shown to be consistent with known fundamental relations in the body–wave interaction theory.     

Fast Evaluation of Time-Domain Green Function for Finite Water Depth

1 .IntroductionShipsrestrainedbycablesandfendersinfrontofdocksundergolargeamplitudenon harmonicmotionsinwaves.Forthiskindofnon harmonicproblem ,atime domainmethodmustbeapplied .LinandYue ( 1 990 )usedanintegralequationwiththetime domainGreenfunctionforinfinitewaterdepthtocomputetheshipmotionindeepwater.Butforthepresentproblem ,theintegralequationwiththetime domainGreenfunctionforfinitewaterdepthmustbeapplied .TheGreenfunctionisafieldwithacertainboundaryandinitialconditionsproducedbyasourceat…     

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.