有限水深三维脉动源格林函数数值算法研究

由于柯西主值积分的奇异性和贝塞尔函数的振荡性,有限水深情况下复杂格林函数及其导数的精确数值求解一直是浅水中波浪水动力计算的难点,因此寻找格林函数的精确数值解显得非常重要.通过对格林函数奇异项进行变形推导,文中给出了一种去掉了奇点的新积分形式.另外通过改进前人推导的Gauss-Laguerre积分方法,给出了一种改进的新Gauss-Laguerre积分方法.格林函数及其导数的数值结果显示文中给出的两种新方法可以有效地计算复杂格林函数及其导数值.最后对这两种方法、级数解以及传统的Gauss-Laguerre积分方法的计算精度和效率进行了比较研究,结果显示文中给出的两种方法的计算精度高于传统的Gauss-Laguerre积分方法,但其计算效率低于级数解.但在接近于0的近场附近级数解失真,而文中提出的改进的新Gauss-Laguerre积分方法可以获得正确结果.同时当频率和水深均较大时级数解失真,而文中提出的方法也可以获得正确结果.最后针对这些方法的优缺点,该文提出了建议的策略用于计算有限水深格林函数.     

时域有限水深格林函数及其导数的数值计算

精确快速地求取有限水深格林函数及它的导数是应用时域法求解有限水深水动力学问题的关键.通过对有限水深格林函数在不同区域渐近特性的分析,给出了使用多维Chebyshev多项式和渐近展开式快速近似计算有限水深格林函数及它的导数的方法.近似方法与直接积分方法的比较表明近似方法的计算工作量大幅减少,可以给出足够精度的计算结果.     

有限水深复合格林函数的数值计算

有限水深格林函数值的正确求解对处于浅水的海洋结构物运动响应和载荷计算起着非常关键的作用.本文根据海洋结构物的物面对称性,在有限水深复合格林函数表达式的基础上推导出了格林函数偏导数的数值计算表达式,并对格林函数的奇异性进行了分析.采用抛物线积分法对其进行了积分,为有限水深海洋结构物的水弹性分析奠定了基础.计算结果表明,利用抛物线积分法所求得的运动响应与挪威DNV船级社SESAM软件的计算结果吻合得很好.     

Chebyshev多项式在频域无限水深格林函数的高效数值逼近中的应用（英文）

文章对频域无限水深无航速格林函数及其导数的高效计算进行了新的研究。该研究针对含有贝塞尔函数的半无穷限柯西主值积分的格林函数的经典表达形式,将整个计算内容分为计算具有简单解析表达形式和缓变剩余函数两部分。为了逼近剩余函数,首先在改进后的三个计算域内得到经济高效的Chebyshev多项式,然后利用Clenshaw法进行该逼近多项式的快速双重求和从而得到格林函数本身。为了达到6D计算精度,其导数的计算也利用了几乎相同的多项式,确切地说是仅需额外增加或减少几项。数值结果表明,文中结果与著名HydroStar的吻合程度很好,而且文中算法与后者的分析效率几乎相同。文中的算法可以提供的数值精度与分析效率完全能够胜任实际应用。     

有限水深频域格林函数的分区快速算法

研究波浪与海洋结构物相互作用问题的关键是如何准确、有效地计算格林函数。本文在Noblesse的无限水深频域格林函数的分区算法的基础上,将分区思想扩展应用到有限水深频域格林函数及其导数的求解中。在计算过程中,对格林函数的空间、时间坐标实行无量纲化,提出积分变量,并且成功推导出符合有限水深边界条件的实函数R_0(h,v)和R_1(h,v)表达式。根据近远场方法,将h-v坐标空间划分为五个区域,分别采取多项式展开、递增级数、泰勒级数与Haskind积分实现格林函数及其导数积分项的快速计算。通过将本文各分区计算结果与经典算法结果进行比较可得本文方法的计算精度较高,至少达到5D。最后,以在波浪中自由运动的浅水漂浮方箱作为算例,采用本文边界元算法分析方箱的动力响应。通过比较发现本文结果与经典算法结果吻合较好,并且与实际物理过程相符,说明本文算法在计算海洋结构物的浅水波浪运动问题时具有良好的工程精度和有效性。     

海洋结构物时II边界积分方程的改进方法

文章基于三维时域格林函数理论,提出了一种时域边界积分方程的改进计算方法,用于预报海洋结构物在波浪中运动受到的时域波浪力.该方法通过利用物体内部构造的满足拉普拉斯方程的内部解,在时域边界积分方程中只需要计及时域格林函数记忆项本身,从而有效地避免了求解具有高震荡积分特性的时域格林函数记忆项的法向导数.对潜体和浮体作了计算,结果表明,这种方法十分有效,可以有效节省内存使用量和计算机空间.通过与相关已发表文献的结果进行比较,计算结果令人满意.     

三维移动源格林函数的新表达式及其数值积分研究

文章以三维移动脉动源格林函数为基础,推导了频率为零有航速时该函数的数学表达式,即定常移动兴波源格林函数,该函数包含Rankine源及其关于静水面的镜像源项、近场扰动项和远场传播项。与一般的移动源格林函数相比,该表达式不存在无穷间断点带来的奇异性,且积分区间变为一般表达式积分区间的一半,在积分计算中具有较大的数值优势。依据近场扰动项和远场传播项的函数特性,提出采用变步长自适应Simpson法和数学变换方法以提高数值积分计算的稳定性,构建了一套完备的兼顾积分效率和精度的数值积分方法。分析新引入的复变函数特性发现,其伪奇异性与Kelvin源传播波的传播范围相对应,且伪奇异性点的个数最多为两个,积分计算过程中可通过寻根公式快速直接求解伪奇异点。数值计算结果表明,文中提出的方法和计算程序可靠,适用于不同航速和任意源点、场点位置条件下的移动源格林函数及其偏导数的数值计算。     

FPSO型采油平台附加质量与阻尼力的分析研究

文章应用边界积分方程的方法研究了FPSO型采油平台的附加质量和阻尼力问题.计算中应用了有限水深格林函数并采用高阶边界元对边界积分方程进行离散计算.研究揭示了FPSO在奇模态运动(纵荡-垂荡-纵摇)和偶模态运动(横荡-横摇-首摇)时附加质量和阻尼力的变化规律.附加质量和阻尼力的计算结果在高频情况下仍具有很好的收敛精度,可作为FPSO型采油平台在设计和使用过程中的荷载分析依据.     

浮体浅水波浪载荷数值计算方法研究

[目的]在浅水环境下,对波浪中浮体的响应进行求解时,其主要难点在于对有限水深格林函数及其偏导数的准确求解和快速计算。[方法]为此,利用改进的Gauss-Laguerre积分法,提出一种可精确计算有限水深格林函数及其偏导数的方法,结合循环矩阵原理,给出对称性的处理方法和简化的级数求解公式。并将所提方法与其他商业软件计算结果进行对比分析。[结果]结果表明,所提计算方法精确度较高。[结论]该方法可用于准确评估浅水中的浮体运动和波浪载荷。     

基于垂向积分形式时域格林函数的多域高阶面元法及其在船舶波阻增加计算中的应用

基于三维时域势流理论,提出基于垂向积分形式的时域格林函数的多域高阶面元法,计算船舶在波浪中的运动响应和阻力增加.在外域采用时域自由面格林函数法,将控制面时域格林函数面元积分分解为垂向积分和水平积分,推导了垂向积分形式核函数新的表达式,对水平积分采用数值离散方法,提高了计算效率和收敛性.在内域采用满足非线性自由面条件的R...     

经过海峡的双层流体由底部波动产生的斜向波散射问题研究(英文)

The problem of oblique wave(internal wave) propagation over a small deformation in a channel flow consisting of two layers was considered.The upper fluid was assumed to be bounded above by a rigid lid,which is an approximation for the free surface,and the lower one was bounded below by an impermeable bottom surface having a small deformation;the channel was unbounded in the horizontal directions.Assuming irrotational motion,the perturbation technique was employed to calculate the first-order corrections of the velocity potential in the two fluids by using Green’s integral theorem suitably with the introduction of appropriate Green’s functions.Those functions help in calculating the reflection and transmission coefficients in terms of integrals involving the shape function representing the bottom deformation.Three-dimensional linear water wave theory was utilized for formulating the relevant boundary value problem.Two special examples of bottom deformation were considered to validate the results.Consideration of a patch of sinusoidal ripples(having the same wave number) shows that the reflection coefficient is an oscillatory function of the ratio of twice the x-component of the wave number to the ripple wave number.When this ratio approaches one,the theory predicts a resonant interaction between the bed and the interface,and the reflection coefficient becomes a multiple of the number of ripples.High reflection of incident wave energy occurs if this number is large.Similar results were observed for a patch of sinusoidal ripples having different wave numbers.It was also observed that for small angles of incidence,the reflected energy is greater compared to other angles of incidence up to.These theoretical observations are supported by graphical results.     

半无限惯性表面产生的波浪散射问题研究的另一种方法(英文)

A new method to solve the boundary value problem arising in the study of scattering of two-dimensional surface water waves by a discontinuity in the surface boundary conditions is presented in this paper. The discontinuity arises due to the floating of two semi-infinite inertial surfaces of different surface densities. Applying Green’s second identity to the potential functions and appropriate Green’s functions, this problem is reduced to solving two coupled Fredholm integral equations with regular kernels. The solutions to these integral equations are used to determine the reflection and the transmission coefficients. The results for the reflection coefficient are presented graphically and are compared to those obtained earlier using other research methods. It is observed from the graphs that the results computed from the present analysis match exactly with the previous results.     

Choquet-like积分的一致自连续性和一致逆自连续性

任意给定非负单调集函数和非负可测函数,Choquet-like积分将确定一个新的集函数。在本文中,我们讨论了一类Choquet-like积分确定的集函数的性质。这类Choquet-like积分中伪乘的单位元不是相应伪加的幂等元。通过引入次可加性、超可加性、一致自连续和一致逆自连续概念的变体形式,我们证明了这类Choquet-like积分保持这些结构特征。这些结果将丰富Choquet-like积分的性质。     

一种三维时域格林函数计算方法

如何高效地计算出三维时域格林函数是线性时域方法能否有效地求解船舶水动力问题的关键.本文在前人工作的基础上,采用了一种半解析的数值求积方法,对三维时域格林函数及其导数作了计算,这种方法十分有效,计算结果令人满意.     

The interaction of oblique flexural gravity waves with a small bottom deformation on a porous ocean-bed: Green’s function approach

The interaction of oblique incident water waves with a small bottom deformation on a porous ocean-bed is examined analytically here within the framework of linear water wave theory. The upper surface of the ocean is assumed to be covered by an infinitely extended thin uniform elastic plate, while the lower surface is bounded by a porous bottom surface having a small deformation. By employing a simplified perturbation analysis, involving a small parameter δ(=1), which measures the smallness of the deformation, the governing Boundary Value Problem (BVP) is reduced to a simpler BVP for the first-order correction of the potential function. This BVP is solved using a method based on Green’s integral theorem with the introduction of suitable Green’s function to obtain the first-order potential, and this potential function is then utilized to calculate the first-order reflection and transmission coefficients in terms of integrals involving the shape function c(x) representing the bottom deformation. Consideration of a patch of sinusoidal ripples shows that when the quotient of twice the component of the incident field wave number propagating just below the elastic plate and the ripple wave number approaches one, the theory predicts a resonant interaction between the bed and the surface below the elastic plate. Again, for small angles of incidence, the reflected wave energy is more as compared to the other angles of incidence. It is also observed that the reflected wave energy is somewhat sensitive to the changes in the flexural rigidity of the elastic plate, the porosity of the bed and the ripple wave numbers. The main advantage of the present study is that the results for the values of reflection and transmission coefficients obtained are found to satisfy the energy-balance relation almost accurately.     

Scattering of oblique surface water waves by thin vertical barrier over undulating bed topography

The present study deals with the scattering of oblique surface water waves by small undulation on the bottom in the presence of a thin vertical barrier. Here, three different configurations of vertical barriers are investigated. Perturbation analysis is employed to determine the physical quantities, namely, the reflection and transmission coefficients. In this analysis, many different Boundary Value Problems (BVPs) are obtained out of which the first two bvps are considered. The zeroth order bvp is solved with the aid of eigenfunction expansion method. The first order reflection and transmission coefficients are derived in terms of the integrals by the method of the Green’s integral theorem. The variation of these coefficients is plotted and analyzed for different physical parameters. Furthermore, the energy balance relation, an important relation in the study of water wave scattering, is derived and checked for assuring the correctness of the numerical results for the present problem.     

三维频域格林函数的高效率计算方法研究

传统的频域格林函数计算方法主要采用级数或渐进展开式对其进行数值逼近,但这种方法需要非常细致的分区,计算过程复杂,循环次数也较多,极大地影响计算效率.将高斯积分引入到频域格林函数及其导数的数值计算,并将其计算结果与已有参考文献进行对比,证明该方法在满足足够精度的基础上,使计算过程大大简化,减少了分区,提高了计算效率.     

水线源强消弃法在船舶定常兴波计算中的应用

当采用等强度源汇分布法求解线性船舶水动力学问题时,船体的水线源强必须等于零,基于这一假定,杜双兴在基于边界元法的数值计算中,令水线附近面元上的源强为零,且水线积分项也自动消失,称这种方法为"水线源强消弃法".本文应用Kelyin移动兴波源格林函数,研究水线源强消弃法在定常航行船舶的兴波计算中的应用,首先以潜航椭球体为例,对Kelvin移动兴波源格林函数进行了计算和验证,然后分别采用常规面元法和水线源强消弃法,对Wigley船型的兴波阻力、波形及船体压力分布进行了计算及试验比较,另外为了考查水线源强消弃法的工程实用性,最后对一艘SWATH船型的兴波阻力进行了计算.结果表明:与考虑水线单元和水线积分的常规面元法相比,水线源强消弃法能够较好地消除水线积分时格林函数的奇异性所引起的误差,给出的预报结果与试验值吻合更好;另外水线源强消弃法的预报结果受网格划分形式、忽略水线单元尺寸大小的影响,文中的算例表明:忽略10%-15%吃水以上的水线单元能够给出较为合理的预报结果.     

有限水深中垂直下潜弹性薄板的水波散射(英文)

The problem of water wave scattering by a thin vertical elastic plate submerged in uniform finite depth water is investigated here.The boundary condition on the elastic plate is derived from the Bernoulli-Euler equation of motion satisfied by the plate.Using the Green’s function technique,from this boundary condition,the normal velocity of the plate is expressed in terms of the difference between the velocity potentials(unknown)across the plate.The two ends of the plate are either clamped or free.The reflection and transmission coefficients are obtained in terms of the integrals’involving combinations of the unknown velocity potential on the two sides of the plate,which satisfy three simultaneous integral equations and are solved numerically.These coefficients are computed numerically for various values of different parameters and depicted graphically against the wave number in a number of figures.     

A note on the solution of water wave scattering problem involving small deformation on a porous channel-bed

The solution of water wave scattering problem involving small deformation on a porous bed in a channel, where the upper surface is bounded above by an infinitely extent rigid horizontal surface, is studied here within the framework of linearized water wave theory. In such a situation, there exists only one mode of waves propagating on the porous surface. A simplified perturbation analysis, involving a small parameter ε ( ? 1), which measures the smallness of the deformation, is employed to reduce the governing Boundary Value Problem (BVP) to a simpler BVP for the first-order correction of the potential function. The first-order potential function and, hence, the first-order reflection and transmission coefficients are obtained by the method based on Fourier transform technique as well as Green’s integral theorem with the introduction of appropriate Green’s function. Two special examples of bottom deformation: the exponentially damped deformation and the sinusoidal ripple bed, are considered to validate the results. For the particular example of a patch of sinusoidal ripples, the resonant interaction between the bed and the upper surface of the fluid is attained in the neighborhood of a singularity, when the ripples wavenumbers of the bottom deformation become approximately twice the components of the incident field wavenumber along the positive x-direction. Also, the main advantage of the present study is that the results for the values of reflection and transmission coefficients are found to satisfy the energy-balance relation almost accurately.