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

文章基于全非线性势流理论对三维液舱晃荡进行了数值模拟,其控制方程由无奇异边界积分方程法(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样条法来光顺波面。同时,在远方辐射控制面上采用多次透射公式方法(Multitransmitting Formula,MTF)来进行消波,文中得到的结果与理论解进行了比较,结果表明该方法可用来有效模拟全非线性的数值波浪水池。     

基于多次透射公式和无奇异边界元法模拟全非线性数值波浪水池

文章基于势流理论对全非线性的三维数值水池进行了模拟,其控制方程由无奇异边界积分方程法(Desingularized Boundary Integral Equation Method,DBIEM)进行离散求解,在求解全非线性的自由面微分方程时,文中采用混合欧拉—拉格朗日法(Mixed Eulerian-Lagrangian,MEL)和四阶Adams-Bashforth-Moulton(ABM4)预报—修正方法,为了避免结果发散即增强数值稳定性,文中采用B样条法来光顺波面.同时,在远方辐射控制面上采用多次透射公式方法(Multi-transmitting Formula,MTF)来进行消波,文中得到的结果与理论解进行了比较,结果表明该方法可用来有效模拟全非线性的数值波浪水池.     

浅水中下潜物体有航速辐射问题的数值解

开发了一种基于非均匀有理B样条(NURBS)的Rankine源高阶面元法用于求解浅水中下潜物体有航速辐射问题。基于势流理论,流场中的扰动速度势用物体表面和自由面上分布的Rankine源来表达,同时用映像法计及水底的影响;采用NURBS曲面精确表达物面和自由面,在求出边界面上的未知源强后,物体表面的速度势分布用B样条表示;采用配置方法对边界积分方程进行求解,采用高斯—勒让德公式计算方程中的积分,同时开发了一种解决数值积分中奇异性问题的方法。采用文中数值方法求解了浅水中下潜圆球的有航速辐射问题,计算得到了下潜圆球水动力系数和自由面波形,所得水动力系数计算结果和他人的计算结果进行了比较,其吻合程度良好,从而验证了文中开发方法的有效性。     

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

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

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

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

基于N-S方程的航行船舶辐射问题数值模拟

基于N-S(Navier-Stokes)方程,进行了航行船舶辐射问题的数值模拟,包括强制垂荡和强制纵摇的模拟.计算了船舶垂荡和纵摇的附加质量和阻尼.数值模拟中,控制方程--BANS(Reynolds Averaged Navier-Stokes)方程和连续性方程使用有限体积法离散,非线性自由面采用VOF方法处理.文中给出了以不同航速前进的船舶强制垂荡和纵摇的力与力矩,以及船舶垂荡和纵摇的附加质量系数和阻尼系数,并与DUT(Delft University of Technology)的试验数据进行了比较,二者吻合良好.     

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

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

数值波浪水池及顶浪中船舶水动力计算

基于粘流理论建立了三维数值波浪水池,模拟了非线性波浪,并对规则波顶浪中前进的拘束船模的水动力进行了计算.数值模拟中,控制方程-RANS方程和连续性方程使用有限体积法离散,非线性自由面采用VOF方法处理;在入口边界模拟柔性造波板运动产生入射波,使用位于波浪水池尾部的人工阻尼区消波.给出了非线性规则波的模拟结果以及规则波顶浪中前进的拘束船模的水动力计算结果,并与理论解及DUT(Delfi University of Technology)的试验数据进行了比较,二者吻合良好.     

考虑定常兴波影响的辐射问题去奇Rankine源方法数值计算

随着中高速船舶在海洋资源开发中的广泛应用,航速对船舶在波浪中摇荡运动的影响受到重视。本文基于势流理论,采用非均匀有理B样条(NURBS)表达物面,应用去奇Rankine源方法求解流场速度势;推导了兴波问题的二阶近似非线性自由面边界条件,以及考虑二阶近似非线性兴波影响的辐射速度势的线性自由面和物面边界条件。本文编程首先数值计算了Wigley船型的近似非线性兴波问题,将计算的线性和非线性兴波阻力系数与文献结果做了比较;然后数值计算了计及线性定常兴波影响的下潜圆球有航速辐射问题,将计算的附加质量和阻尼系数与文献结果比较。本文计算结果与文献结果吻合程度较好,验证了本文数值计算方法的有效性。文中还给出了兴波波形和辐射波形。     

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

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

离散格式和网格数量对过渡型艇粘性阻力的影响

采用对比性数值分析方法研究高速气泡船减阻效果，针对某优良过渡型艇，为喷气需要预先进行断阶，不计自由面影响，采用有限体积法和SIMPLE算法，湍流模型为k-ε两方程，耦合求解N-S方程和连续性方程，数值分析不同离散格式、不同网格数量对艇体在不喷气状态下遭受的特性阻力影响。计算结果显示：二阶迎风格式比一阶迎风格式数值结果与试验结果符合较好；网格数量达到95万，计算结果与试验结果符合较好。     

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

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

耦合VOF的非线性k-ε模型三维溃坝数值模拟

溃坝是一种伴随有巨大破坏力的灾害性的水流现象。三维溃坝的数值模拟能够为溃坝水流的研究分析提供科学的依据。针对瞬间全溃模型水流的复杂性,考虑水库下游障碍物的滞水及阻水作用,利用开源程序OpenFOAM建立了耦合VOF对流方程的三维非线性(NL)k-ε紊流数学模型求解雷诺时均纳维斯托克斯方程。采用有限体积法进行离散,利用PISO算法进行数值求解。通过模拟下游为三角形障碍物、矩形障碍物以及90°弯道的溃坝实验,并将模拟结果与实验结果进行对比分析,初步证明了该模型的准确性和鲁棒性。     

基于体积力法的全附体KCS型船模PMM运动数值模拟

为高效、精确地求解船舶操纵水动力导数,以全附体KCS型船模为研究对象,基于RANS方程及VOF模型,在star-ccm+平台上采用体积力法模拟螺旋桨作用,计及自由面兴波及航行姿态变化对船模水动力的影响,开展斜航、纯艏摇、漂角和艏摇组合3种平面运动机构(Planar Motion Mechanism,PMM)运动的数值模拟,将横向力Y、艏摇力矩N、横倾力矩K的数值模拟结果与试验结果进行对比,并依照PMM试验的动力学方程,通过最小二乘法拟合和傅立叶积分,将数值模拟结果与试验结果的时历曲线进行分析处理,求得操纵水动力导数。研究结果表明:该方法用于数值模拟PMM运动是可行的;求得的水动力导数中Y′vvv,N′vvv误差略大,其余水动力导数误差均在15%以内。     

高速三体船波浪中运动与增阻CFD计算研究

基于数值波浪水池技术,对波浪中高速三体船运动及增阻进行CFD计算研究。控制方程—RANS方程和连续性方程使用有限体积法离散,非线性自由面采用VOF方法处理;在入口边界模拟柔性造波板运动产生入射波,使用位于波浪水池尾部的人工阻尼区消波。首先对规则波顶浪中单个主船体的运动和增阻进行了计算,并与高速细长体理论计算结果进行了比较;随后进行了三体船运动和增阻计算,分析了侧片体对主船体阻力增加的影响。为高速三体船耐波性研究提供了数值工具。     

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

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

高阶呼吸子式畸形波波形的数值研究

畸形波是一种波高极大的罕见波浪,对海洋平台的结构以及工作人员的人身安全具有极大威胁。由于在开阔海域中对畸形波的观测存在困难,通过数值手段对其波形进行研究具有十分重要的意义。四阶非线性薛定谔方程又被称为Dysthe方程,能够更加准确地描述波浪的非线性演化。该文通过数值求解Dysthe方程,研究了三阶非线性薛定谔方程的高阶呼吸子解(最高至5阶解)的演化过程。提出了四阶离散步长—虚拟波谱法,得到了具有更高守恒精度的数值解。同时,通过坐标变换以及拟合对高阶呼吸子式畸形波的波形进行了深入研究。     

电厂贮灰坝的渗流数值分析

建立了饱和—非饱和渗流模型,采用空间有限元离散、隐式时间差分方法求解描述饱和—非饱和渗流场非线性抛物型方程,提出了计算电厂贮灰坝渗流的数值计算方法。将辽源发电厂贮灰坝进行二维渗流数值计算,计算结果与水电模拟实验结果进行比较,表明结果是可靠的。     

基于OpenFOAM的全回转推进器非定常性能预估

本文选择某一全回转吊舱推进器为研究对象,运用任意网格面元法(AMI)CFD数值计算方法,计算采用的开源数值计算软件OpenFoam对全回转推进器非定常敞水水动力性能进行了研究分析.控制方程离散采用有限体积法,湍流模型选取SSTκ-ω 模型,速度压力耦合采用PIMPLE算法求解.通过与试验结果比较,验证了基于OpenFoam的非定常数值计算方法在全回转推进器敞水性能预报的可靠性和有效性.