绕船体粘性自由面流动的数值计算

考虑自由面的船舶粘性流动计算由于其复杂性和对计算机硬伯的较高要求，只是到了最近几年才得到较大发展，而在国内尚属起步阶段。本文使用VOF方法计算了绕Wigley船型的粘性自由面流动，流模式选择了SSTκ-ω模型。计算获得了船体上波形以及船波等值线图，同时考察了三个不同横截面上的流动状况。计算得到的阻力还和发表的试验结果作了对比，结果令人鼓舞。     

船波共振对锚泊船的危害和化解方法之探讨

在世界船运史中，锚泊船倾覆的悲剧时有发生。造成倾覆的主要原因除了锚泊船遭遇到风压倾覆力矩大于船舶稳定性力矩这个不可能抗拒的外力之外，作者认为在大风浪中如果陷入船波共振也有导致锚泊船倾覆的可能。一旦遇到这类险情，只要我们能及时采取化解船波共振的操作方法就能迅速使船舶化险为夷。本文试对形成船波共振的条件和化解船波共振的方法进行探讨，供同行们参考。     

不对称流动对航道中航行船舶的影响

本文利用能够处理不对称自由面流动的改进的Dawson's方法来研究当船偏离航道航行时航道岸壁及水深对船舶所受的水动力的影响.计算结果与模型试验进行了比较,结果表明当航道水深不是很浅时,本文提出的方法和计算程序可以给出比较满意的结果,同时还可以发现当船舶偏离中心航道航行时,船两侧湿表面压力以及水面变化也是不对称的,水深以及岸边距离都能影响船体所受水动力及船波.     

基于船波相对运动的砰击压力和上浪载荷预报方法

本文基于船波相对运动理论,对舷侧砰击压力和甲板上浪载荷的数值及试验预报方法分别进行了研究.利用三维势流理论计算船舶与波浪之间的相对运动,可以得到船波相对速度及甲板上浪高度.对于砰击压力通过数值模拟方法得到砰击压力系数后结合船波相对速度来预报;对于上浪载荷则采用考虑船舶航速的溃坝模型结合甲板上浪高度来预报.此外,开展了船舶运动和砰击压力模型试验,船波相对运动由沿着模型横剖面布置在舷侧的浪高仪测量,并且测得了相应位置处的砰击压力.最后分别对船波相对运动和砰击压力的数值结果与试验数据进行了比较分析,同时基于船波相对运动给出了一种砰击持续时间计算方法.     

浅窄航道非线性兴波阻力数值计算

应用Rankine面元法计算船舶在浅窄航道中作匀速直航运动时的非线性兴波阻力．基于势流理论，在船体表面、自由面以及岸壁分布Rankine源来表达流场速度势，利用迭代方法满足非线性自由面边界条件，采用映像法满足水底边界条件，采用网格错位法满足辐射条件．以Wigley数学船型为例进行计算，通过改变航速和航道宽度参数，得到相应的兴波阻力．通过分析不同航速下航道宽度变化对兴波阻力的影响，得到了一些有意义的结论．     

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

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

一种丰满船型的兴波阻力计算方法

提出一种部分计及自由面粘性效应及表面张力影响的完全非线性自由面边界条件,以该条件计算了方形系数Cb为0.8的Todd 60船型的非线性兴波问题,计算结果表明:非线性迭代的收敛性得到改善,兴波阻力被夸大的倾向有所改进。     

非线性兴波数值方法在高速三体船侧体布局方案比较中的应用

为使三体船侧体布局方案阻力性能比较更加可靠,对单体船非线性自由面边界条件的兴波阻力预报程序进行改进,成功开发了考虑升沉和纵倾影响的三体船非线性兴波阻力预报程序,对15种侧体布局方案的Wigley三体船兴波阻力进行计算,并将其与线性薄船理论计算结果、线性自由面条件的Dawson型方法计算结果、模型试验结果进行比较,表明非线性兴波阻力数值方法能够更好地反映三体船兴波阻力特性,用于三体船侧体布局优化设计时较其他2种数值方法更加可靠.     

近水面椭球体非线性兴波的数值模拟

本文开发了一种去奇异面元法用于计算以定常前进速度运动的近水面椭球体引起的三维非线性自由面势流.采用求解非线性深水波的Stokes理论求解该定常非线性问题;物面用NURBS(非均匀有理B样条)表达;孤立的Rankine源分布在物体内部和自由面上方.给出了波型、波阻以及收敛特性随一些参数的变化规律.波阻计算结果和他人的其它数值结果进行了比较,其吻合程度令人满意.     

船舶砰击强度计算方法研究

船舶在航行时发生砰击是一种高度非线性的物理现象。本文基于势流理论Rankine源法,研究用时域非线性方法计算船舶运动响应,获得船波垂向相对运动速度;并选取合适的砰击压力系数,计算船舶砰击载荷。最后用该方法评估了一艘集装箱船艏部砰击强度,验证了该方法的合理性。     

近自由面三维水翼的水动力分析及试验研究

采用面元法对近自由面三维水翼进行势流数值分析并进行了相关试验研究。在数值计算中,将Rankine源和偶极子置于边界面上,用时间步进法模拟水翼的势流场和自由表面波形。在自由面采用非线性自由面边界条件,在尾涡面上采用偶极子布置以满足Kutta条件。文中给出了数值计算模型的参数,对于不同浸深、不同航速和不同攻角下的水翼,计算了水翼表面上的压力分布,水翼的阻力和升力及自由表面波形。数值计算结果与试验结果进行了对比。结果表明,文中方法可用于水翼优化设计、近自由面振动翼运动及水翼船兴波等问题的研究。     

采用Green-Naghdi方程和有限元法计算浅水船波

采用波动方程／有限元法求解Green-Naghdi(G-N)方程计算船舶在有限水深区域的兴波和波浪阻力。把行驶船舶对水面的扰动作为移动压力直接加在Green-Naghdi方程里，以描述运动船体和水面的相互作用，并经此来计算不面波动、船底水动压力和波浪阻力。G－N方程比浅水方程增加一个非线性的频散项，以补充有限水深对浅水船波的影响。采用随船运动网格的有限方法，以Series 60 CB=0.6船作为算例给出浅水船波的计算结果，并与浅水方程的结果进行了比较。计算结果表明，当船速小于临界速度时，由于频散的影响，G－N方程级出的船后尾波波高比浅水方程的结果大，同时波浪阻力也比浅水方程的结果有所提高。当船速大于临界速度时，G－N方程的计算结果与浅水方程基本相同，频率散射无明显影响。     

基于无奇异边界元法模拟三维全非线性液舱晃荡

文章基于全非线性势流理论对三维液舱晃荡进行了数值模拟,其控制方程由无奇异边界积分方程法(Desingular-ized Boundary Integral Equation Method,DBIEM)进行离散求解,在求解全非线性的自由面微分方程时,文中采用混合欧拉—拉格朗日法(Mixed Eulerian-Lagrangian,MEL)和四阶Adams-Bashforth-Moulton(ABM4)预报—修正方法,为了避免结果发散即增强数值稳定性,文中采用B样条法来光顺自由面.在微幅水平激励下,该文中得到的结果与解析解吻合较好.     

基于无奇异边界元法模拟三维全非线性液舱晃荡（英文）

文章基于全非线性势流理论对三维液舱晃荡进行了数值模拟,其控制方程由无奇异边界积分方程法(Desingularized Boundary Integral Equation Method,DBIEM)进行离散求解,在求解全非线性的自由面微分方程时,文中采用混合欧拉—拉格朗日法(Mixed Eulerian-Lagrangian,MEL)和四阶Adams-Bashforth-Moulton(ABM4)预报—修正方法,为了避免结果发散即增强数值稳定性,文中采用B样条法来光顺自由面。在微幅水平激励下,该文中得到的结果与解析解吻合较好。     

MPS方法在细长河道中溃坝流动的应用研究

溃坝问题是我国防灾减灾工作中的一个突出问题.溃坝发生时,会对下游的结构物产生剧烈的砰击作用,并伴随着自由液面的翻卷、破碎等强非线性现象.对这种短时间的自由面强非线性变化,基于网格的计算流体方法往往需要一些特殊的数值处理.文中采用本课题组自主开发的无网格粒子法求解器MLParticle-SJTU对二维带斜坡、三维带90°急弯的溃坝问题进行了模拟.研究了溃坝流的自由面演化,分析了溃坝流遭遇斜坡障碍物的复杂流动现象,并记录下监测点的波高时间历程.通过MPS的数值模拟结果与试验数据进行对比,两者吻合较好.     

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

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

船波自动测试及实时处理系统

本文介绍上海船舶运输科学研究所2号水池以KGY-5A浪高仪和LSI-11/23微机为核心的船波自动测试及实时处理系统。对船波测量、信号传输、激光定位、实时采样及数据处理等部分作了详细说明。测试结果的分析表明,本系统是可靠的,不但大大提高了工作效率,而且提高了测量精度。     

基于船波相对运动的船艏砰击仿真方法

在分析船波相对运动表达式的基础上计算船艏典型剖面的船波相对运动,探讨船艏入水过程中的砰击问题,对比船体某剖面3种入水仿真模型计算所得的砰击载荷,讨论三维外形和航行速度对船艏剖面砰击外载荷的影响。在该分析中,船体在规则波浪中的运动用基于三维势流理论的水动力软件AQWA计算获得,船波相对运动通过理论推导计算获得,用对船艏结构施加强迫运动的方式模拟船波相对运动的真实过程。采用An—sys／Ls—Dyna软件的流固耦合分析进行入水仿真,流体划分为ALE体积网格,船艏视为刚体,划分为Lagrange有限元网格。对比结果表明：在三维模型中,相邻剖面引起剖面最大压力点处的液面变化对该点的砰击压力有增大效果,航速有增大剖面砰击压力的作用,减小船艏底部纵向斜升角有利于降低砰击压力。     

基于船波相对运动的船艏砰击仿真方法

在分析船波相对运动表达式的基础上计算船艏典型剖面的船波相对运动,探讨船艏入水过程中的砰击问题,对比船体某剖面3种入水仿真模型计算所得的砰击载荷,讨论三维外形和航行速度对船艏剖面砰击外载荷的影响。在该分析中,船体在规则波浪中的运动用基于三维势流理论的水动力软件AQWA计算获得,船波相对运动通过理论推导计算获得,用对船艏结构施加强迫运动的方式模拟船波相对运动的真实过程。采用An sys/Ls-Dyna软件的流固耦合分析进行入水仿真,流体划分为ALE体积网格,船艏视为刚体,划分为Lagrange有限元网格。对比结果表明：在三维模型中,相邻剖面引起剖面最大压力点处的液面变化对该点的砰击压力有增大效果,航速有增大剖面砰击压力的作用,减小船艏底部纵向斜升角有利于降低砰击压力。     

解非线性船波问题的一个数值计算方法

本文提出了一种计算完全非线性三维船渡问题的数值计算方法。速度势用考虑了波动的自由面影响的船体湿表面和局部真实的自由面上的简单源(1/r)分布表示,并采用精确的自由面条件。求解过程是一个迭代过程,计算结果与试验值作了比较,吻合较好.