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

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

水面舰船附近气泡运动数值模拟

当气泡在自由面和水面舰船附近运动时,两者的存在均会对气泡脉动产生影响。基于不可压缩势流理论,采用边界元方法对边界积分方程进行求解。针对建立自由面需要大量的网格,且处理自由面与水面结构交界面时数值不稳定,采用考虑自由面效应的格林函数取代基本格林函数,通过与自由面附近气泡轴对称模型的计算结果进行对比,验证了该方法的有效性。通过对舰船结构附近的气泡运动射流特性进行模拟,发现当药包在舷侧和自由面附近爆炸时,气泡射流可能不会完全作用在舰船上,自由面效应的存在削弱了气泡的打击能力。     

波浪中航行两船模池壁干扰效应的数值研究

对规则波中顶浪并行航行的两船模,基于频域移动脉动源格林函数建立了有效表达池壁效应的三维数学模型,该模型采用分布源积分方法求解近距航行两船模的波浪作用力和运动响应,并利用镜像原理将水池格林函数表示为开阔水域的格林函数和由镜像点源组成的无穷级数.分析移动脉动源格林函数的远场传播模态及其特性可知,当τ>0.272(τ为斯特劳哈尔数)时各传播波模态均不在点源的上游出现,故此时点源波系的影响范围可由环形—扇形波的半楔角确定,在此基础上计及船模辐射—绕射波的池壁反射作用对自身和另一船模的干扰,依据船型要素、水池宽度和两船相对位置,采用六个几何临界半楔角表征船体的辐射—绕射波,提出了并行航行两船模的池壁效应判别方法.数值研究表明:不同镜像次数下水池格林函数的数值差异主要体现在距船艏较远处,在兼顾计算精度与效率的原则下可选取有限项级数的水池格林函数作为积分内核,而池壁效应对近距航行两船运动响应的影响主要体现在长波区.     

波浪中航行两船模池壁干扰效应的数值研究（英文）

对规则波中顶浪并行航行的两船模,基于频域移动脉动源格林函数建立了有效表达池壁效应的三维数学模型,该模型采用分布源积分方法求解近距航行两船模的波浪作用力和运动响应,并利用镜像原理将水池格林函数表示为开阔水域的格林函数和由镜像点源组成的无穷级数。分析移动脉动源格林函数的远场传播模态及其特性可知,当τ0.272(τ为斯特劳哈尔数)时各传播波模态均不在点源的上游出现,故此时点源波系的影响范围可由环形—扇形波的半楔角确定,在此基础上计及船模辐射—绕射波的池壁反射作用对自身和另一船模的干扰,依据船型要素、水池宽度和两船相对位置,采用六个几何临界半楔角表征船体的辐射—绕射波,提出了并行航行两船模的池壁效应判别方法。数值研究表明:不同镜像次数下水池格林函数的数值差异主要体现在距船艏较远处,在兼顾计算精度与效率的原则下可选取有限项级数的水池格林函数作为积分内核,而池壁效应对近距航行两船运动响应的影响主要体现在长波区。     

基于三维时域格林函数高阶面元法的船体运动计算

针对波浪中航行船体的运动问题,基于三维时域自由面格林函数并采用高阶面元离散,编程计算了浮体的运动响应.采用自由面格林函数,只需在船体湿表面布源,减少系数矩阵维度,避开自由面网格划分;同时高阶面元法(HOBEM)在计算速度势空间分布时能够提高精度.本文通过计算Wigley-Ⅲ型船在波浪中有航速前行时的运动,与实验进行对比,发现结果符合较好,证明这种方法的有效性.     

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

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

应用混合格林函数法计算波浪中浮体运动及离散参数取值的算例分析

通过假想的直壁控制面将流场分割成内域和外域,分别引入Rankine源和自由面Green函数并结合控制面上的连续条件进行耦合求解,应用混合格林函数法对浮体在波浪中的运动进行了时域计算。以wigley型船为例,针对不同的时间步长、自由面和控制面的不同区域范围、面元划分的不同密度等进行了一系列数值计算和比较,讨论了离散区域大小、网格密度等参数对计算结果的影响,并给出其合理的取值范围。在此基础上,对一艘LNG船在不同浪向入射波中的运动响应进行计算,并与试验数据对比。结果证明,该方法能有效解决外飘船型时域计算的发散问题,可以对浮体在斜浪中的运动进行时域模拟。     

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

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

基于多重透射改进的二维半理论及水动力系数预报

水动力系数是船舶耐波性能预报、操纵性能预报,滑行艇在静水中的运动稳定性及海豚运动预报的重要参数.探讨高效、稳定的水动力系数预报方法对准确预报上述问题具有重要的意义.该研究将地震波领域的多重透射公式应用于简单格林函数二维半理论,探讨了远方边界条件,用四阶龙阁库塔法处理自由面条件,进行详细的多重透射阶数、横剖面数、附加因子、自由面网格数量以及计算域的大小对计算结果的影响的探讨,并以Wigley船型为例计算了静水航行并做强迫振荡时的水动力系数.数值计算结果与切片方法、二维半理论复杂格林函数方法、三维频域方法及和试验结果进行了比较,比较表明多重透射公式能够有效地模拟二维半远方辐射条件,提高计算效率和稳定性,对二维半理论简单格林函数方法是一种有效的改进.     

基于三维非线性区域分割法的船舶运动时域计算研究

基于三维势流理论,构建船舶在波浪中航行的非线性区域分割法。用假想的控制面将计算域分为内域和外域,其中外域采用自由面格林函数法,内域考虑自由面非线性,采用满足移动兴波线性化自由面条件的Rankine面元法,并通过过渡函数使内外域速度势及其导数在控制面上连续。内外域耦合求解得到扰动速度势及流体动压力。其中自由面运动学边界条件、动力学边界条件及物面条件m项中的移动兴波速度势及波面抬升是基于非线性兴波计算得到的,较线性化自由面条件更为精确。以数学船型Wigley-3为对象,讨论了不同因素对计算收敛性的影响,验证了该方法的可靠性和准确性;以S60和S175两种具有外飘的实际船型在波浪中的运动响应为例,将计算结果与试验及其他方法进行了对比,验证了三维时域非线性区域分割法在实船运动预报中的准确性。     

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

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

直接边界元法时域全非线性兴波计算

文章在时域中采用直接边界元的方法计算模拟了以稳定航速运动的近水面潜体全非线性兴波现象.从格林定理出发,在边界面上布置Rankine源和偶极子,采用时间步进的方法计算模拟潜体近自由面运动的兴波现象.每个时步自由面上的势和节点位置采用全非线性拉格朗日自由面边界条件求解,重新生成自由面网格准备下一个时步计算的初始条件.文中给出了数值计算的结果并与试验和线性计算的结果相比,证实了该方法是可靠的,并可以利用该方法建立数值波浪水池.     

格林函数剩余区间高阶导数数值计算研究

点-偶混合分布法求解三维频域有航速船舶水动力时会涉及到自由面格林函数高阶导数的计算。为提高计算精度,将区间划分为普通解析计算区间和缓变剩余区间,并对剩余区间进行双重Chebyshev逼近,得到若干子区域的Chebyshev系数。通过数值计算算例,表明采用双重Chebyshev级数逼近剩余区间内的格林函数及高阶偏导数,至少达6D精度。     

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

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

基于多层面改进多域法的航行船舶时域水动力分析

针对航行船舶三维时域水动力分析问题，本文采用多域高阶边界元法（MDHOBEM, Multi-Domain Higher Order Boundary Element Method）开展研究，通过假设控制面，将计算域分为内域和外域两部分，内域采用Rankine面元法，外域使用自由面格林函数法，内外域联立求解。从内域、外域、波浪激励三个层面改进多域法的精度、计算效率，以适应不同工况需求。内域重新推导非线性边界条件，考虑入射波对绕辐射问题计算的影响，提高大幅波浪中运动计算精度；外域采用垂向积分形式自由面格林函数，提高格林函数积分计算的收敛性、稳定性及计算效率；波浪激励采用高阶谱（HOS, High Order Spectral）方法生成演化的波场，考虑波浪自身的非线性，提高对非线性现象的捕捉和模拟能力。以C11船为对象，研究内域非线性改进对船舶大幅运动计算精度的提高；以S175集装箱船为对象，验证外域垂向积分对运动及波阻增加计算效率和稳定性的提高；对比S175在线性波场和非线性波场中的运动，论证考虑非线性波场对运动计算的影响；以某内倾船为对象，研究不同层次改进对参数横摇现象捕捉的影响，验证非...     

船舶有航速时三维时域格林函数计算方法研究

船舶工业CAE软件势流求解器的开发对势流计算方法的可靠性提出要求,本文针对船舶有航速时深水三维时域格林函数计算方法开展研究。首先引入三维时域自由面格林函数,其瞬时项的求解采用Hess&Smith方法,自由面记忆项的求解采用美国密歇根大学Beck团队的方法。其次,分区域推导出三维时域自由面格林函数记忆项及其导数的编程所用的数学展开表达式。然后,采用间接法模型求解分布源强和速度势,并采用压力直接积分法求出时域辐射力。最后,针对Wigley I船型开展了辐射力计算,并与公开的试验及计算结果进行对比,验证本文方法的可靠性。     

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

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

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

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

规则迎浪中船舶参数横摇的三维时域预报方法研究

参数横摇是船舶因复原特性改变而引起的典型非线性现象。文章采用三维时域方法预报规则迎浪中船舶的参数横摇运动。该方法引入匹配面将流域分为内域和外域,内域中采用Rankine源来满足物面条件和线性自由面条件,而外域中应用时域格林函数来满足线性自由面条件和远场辐射条件。数值方法中,Froude-Krylov力和恢复力是通过对船舶瞬时湿表面积积分获得,同时考虑了横摇、垂荡和纵摇三自由度之间的耦合作用,以及非线性横摇阻尼的影响。数值结果与试验结果吻合很好,说明该方法可以有效地预报参数横摇。     

规则迎浪中船舶参数横摇的三维时域预报方法研究（英文）

参数横摇是船舶因复原特性改变而引起的典型非线性现象。文章采用三维时域方法预报规则迎浪中船舶的参数横摇运动。该方法引入匹配面将流域分为内域和外域,内域中采用Rankine源来满足物面条件和线性自由面条件,而外域中应用时域格林函数来满足线性自由面条件和远场辐射条件。数值方法中,Froude-Krylov力和恢复力是通过对船舶瞬时湿表面积积分获得,同时考虑了横摇、垂荡和纵摇三自由度之间的耦合作用,以及非线性横摇阻尼的影响。数值结果与试验结果吻合很好,说明该方法可以有效地预报参数横摇。