满足自由水面的水槽格林函数的研究

研究满足自由水面的水槽格林函数的两种数值计算方法,一是Linton提出的快速收敛数值计算算法,二是对John提出的开阔海域中格林函数关于水槽两侧壁作镜像的数值计算方法。研究表明,Linton提出的快速收敛数值计算算法在特定情况下不适用,因此建议采用镜像格林函数来计算常水深和常宽度水槽中波浪对结构的作用问题。     

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

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

船舶有航速辐射问题的新实用数值方法

为解决船舶在波浪中以恒定速度运动的辐射问题,基于势流理论,采用三维瞬时自由面格林函数方法对船舶水动力问题进行求解。提出三维瞬时自由面格林函数波动项的新数值计算方法,并在大时间区域内基于精细积分法对三维瞬时自由面格林函数的波动项进行插值调用,最终提高波动时域格林函数的数值计算精度。对Wigley船型的水动力系数进行数值计算,并将附加质量和阻尼系数的数值计算结果与试验结果和参考文献进行对比,验证新数值计算方法的有效性。     

基于Rankine源和Kelvin源格林函数求解兴波阻力的复合算法

[目的]运用边界元法计算船舶兴波阻力基本上是先求解船体附近的速度分布,然后采用伯努利方程进行压力积分,其计算过程复杂,且误差非常大。[方法]提出一种可快速计算船舶兴波阻力的复合算法,利用Rankine源格林函数求解船体表面源强,结合Lagally定理进行受力计算,并基于Kelvin源格林函数求解船舶兴波阻力。运用该算法对Wigley船的兴波阻力进行计算。[结果]计算结果表明,所用算法相较于运用线性兴波阻力中的薄船理论得到的结果精度更高,而且与完全使用Kelvin源格林函数的算法相比效率也更高。[结论]所用算法可在计算兴波阻力时作为精度与效率之间的一种折中方法。     

深水波频域格林函数及其高阶导数算法研究

文章研究了无限水深自由面格林函数一种新的算法,该算法将计算域重新划分为五个子区域,在每一个子区域内只涉及级数展开,没有使用插值表及数值积分,与现有格林函数算法相比,不仅给出一阶导数值,而且能高效地计算二阶导数,这在格林函数微分方程求解及高阶面元法中非常有用.数值结果表明,该方法计算速度快且一阶导数和二阶导数的计算精度达到10-9.     

基于机器学习的三维频域格林函数递推预报与应用

频域自由面格林函数法在船体水动力的求解中具有重要应用价值,求解的关键在于自由面格林函数本身及其偏导数的计算。论文分析了自由面格林函数与其偏导数之间的关系,提出了基于机器学习的格林函数及其偏导数的递推预报方法,称为脉动源递推预报方法。研究了在不同参数设置下用该方法预报格林函数的精度和效率,并将其应用于S175船的水动力计算。计算得到的自由面格林函数的精度在10^-3～10^-6;计算效率优于数值方法,与解析方法相当;水动力计算结果准确,计算方法具有实用价值。     

基于CFD实尺2339标准箱集装箱船的阻力分析和实船验证

基于计算流体力学(CFD)方法对2339标准箱支线型集装箱船进行实尺数值模拟阻力预报。首先建立NUMECA/HEXPRESS全六面体非结构网格,采用湍流K-Omega(SST)-Menter模型进行数值模拟,计算其在压载工况下不同傅汝德数时的阻力;然后将CFD数值计算阻力结果同实船试航阻力进行比较与分析。结果表明:CFD数值模拟同实船以及实验结果趋势一致,误差较小。论文所述CFD分析方法对船舶的线型设计、优化和航速预报,具有一定的指导与参考意义。     

基于SHIPFLOW的KCS船型兴波阻力数值计算

本文以韩国船舶与海洋工程研究所(KRISO)的集装箱船标模KCS为研究对象,利用CFD软件SHIPFLOW对其进行兴波阻力性能数值模拟,比较网格质量、傅汝德数对兴波阻力性能计算结果的影响。结果表明:傅汝德数相同,网格质量不同时,计算结果不同;网格质量相同,傅汝德数不同,计算结果与实验值有所差异,但是总体趋势与实验值较为一致。本文对不同网格质量、不同傅汝德数下的KCS船型作了兴波阻力的数值模拟,模拟结果具有一定可靠性,为基于兴波阻力的船型优化提供兴波阻力性能数值计算支撑。     

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

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

近水面水翼影响的船舶兴波时域计算研究

近水面水翼的兴波由升力线的自由面兴波来表达,船体的兴波考虑升力线兴波的干扰影响,三维时域格林函数法用于船体的兴波分析计算.以Wigley船型为数值计算例,考虑近水面水翼的影响,分别对单体船和双体船进行了船舶兴波计算研究,提供了兴波波形和兴波阻力的计算结果.以系列60船型为数值计算例,把计算的结果与发表的有关结果进行对比,具有较好的吻合.本文还对船舶自由表面兴波波形进行了时域数值仿真.     

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

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

三维脉动源格林函数的一种简化算法

在三维脉动源格林函数的数值计算中,由于格林函数中通常含有形式复杂的函数,如Bessel函数,会导致计算比较困难。本文从三维脉动源格林函数中所包含的三种Bessel函数积分形式入手,采用了可变误差多面体算法来拟和形式复杂的函数,提出了三维脉动源格林函数的一种简化算法。采用这一简化算法大大减少了计算量,并且精度可以满足工程需要。     

基于Rankine源函数的船舶兴波阻力计算

本文基于格林函数的Rankine源函数对船舶运动中的兴波阻力展开计算,由此提升直接法的计算精度,保障船舶的航行安全。     

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

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

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

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

极地船连续破冰阻力数值模拟研究

连续破冰是自破冰型极地船舶中最常用的破冰模式,连续破冰阻力精确计算是极地船舶航行性能预报和船型优化设计的基础。论文提出一种基于预设冰网格技术的连续破冰阻力数值模拟方法,预先对船体水线附近的线型进行分层离散及对冰场进行方形网格离散,通过几何计算和受力分析数值模拟连续破冰进程及产生的破冰阻力分量。进一步考虑连续破冰阻力与运动耦合效应,在时域内建立船体6自由度运动微分方程,数值求解获得连续破冰阻力时间历程结果。利用该数值模拟方法计算获得了一艘PC3级极地油船连续破冰阻力时间历程结果,并将冰阻力有义值与对应的模型试验结果进行了对比。研究表明,数值模拟获得的连续破冰阻力时历结果能够反映连续破冰进程,数值模拟获得的冰阻力有义值与模型试验结果的误差在12%以内。论文的数值模拟方法能够高效、精确地模拟连续破冰进程及冰阻力,具有一定的工程应用价值。     

基于精细时程积分方法的时域格林函数计算

文章提出了一种精确计算三维时域格林函数记忆项本身及其导数的方法。基于精细时程方法求解时域格林函数及其导数所满足的定常非齐次常微分方程组,得到时域格林函数数值解。文中通过与解析方法和其他方法的比较,验证了该方法的精确性、稳定性和高效性。     

基于时域格林函数法的船舶有航速问题数值求解

为了解决船舶以恒定速度在波浪中运动的水动力问题,文章开发了基于瞬时格林函数方法的三维时域计算程序。从数值精度和计算效率的角度出发,采用了精细时程法和一些数学上的推导计算时域格林函数的波动部分及其导数。并对Wigley型船进行了数值计算,通过水动力系数与试验数据及文献结果的比较,验证了三维时域计算程序的可靠性。     

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

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

基于Kelvin源格林函数的船舶兴波阻力预报

实现了基于Kelvin源格林函数的实船模型兴波阻力预报,采用自适应数值积分算法处理被积函数的高频振荡。首先介绍了计算所用的数学模型及数值方法,而后将该方法应用于Wigley数学船型,初步验证了方法和程序的正确性.在此基础上计算系列60船型兴波阻力系数,对比试验结果分析了面元上Gauss节点数与数值积分算法选择对计算精度的影响,得出不同数值计算方案优劣的定性结论。针对方尾水面舰船标模DTMB 5415计算兴波阻力系数,所得结果在大部分速度区间内误差控制在5%以下。