网格节点型有限体积法在滑动轴承润滑特性求解中的应用

本文采用网格节点型有限体积法求解二维瞬态雷诺方程,采用矢量化方法离散雷诺方程,从而构建并得到新的数值模型,该模型既有有限体积法守恒性的特点,又能适用于结构复杂的求解域.通过与相关算例的对比验证,证明了该模型可适用于滑动轴承润滑特性的研究,同时精度也可以得到保证;根据计算得到的压力分布特点,采用不同的网格类型和网格数量重新划分求解域,并对计算结果和时间进行了对比,结果表明可以采用不同的稀疏网格来提高计算效率;通过计算圆锥形织构轴承的压力分布,证实了该模型对复杂求解域的良好适用性.     

基于双向势流—黏流耦合的数值消波模拟

提出一种基于双向势流-黏流耦合的数值消波方法,将计算域分为靠近浮体的黏流计算域和远离浮体的势流计算域,并通过重叠区域进行流场信息匹配。在黏流计算域中求解N-S方程,采用有限体积法对控制方程离散,采用几何VOF方法捕捉自由液面。在势流计算域中求解Laplace方程,引入线性自由液面假设,采用Rankine源法求解。在重叠区域内通过WheelerStretching方法和松弛区域法进行速度场双向信息交换。基于该方法,在OpenFOAM框架下开发了求解器,并在模型尺度下,对不同参数的规则波进行了模拟,还与waves2Foam求解器的计算结果进行对比。计算结果表明,在同等消波效果下,基于论文方法的数值模拟所需的黏流计算域消波区长度为waves2Foam求解器的1/5。     

同步入水双圆柱相互作用的计算研究

为了研究多圆柱入水的相互影响,采用二维CFD方法数值模拟双圆柱同步常速入水物理过程.利用自由液面捕捉法来计算自由液面,将自由液面当作密度域的接触间断来处理,以非均匀密度的不可压缩Euler方程作为控制方程.采用直角切割网格作为离散网格系统和有限体积法进行数值离散,以人工压缩法来实现速度和压力的耦合求解.采用Roe的近似Riemann解算子计算流甬数和一阶隐式的双时间推进法进行时间步进.以三个不同参数的双圆柱入水作为算例,将数值模拟得到的自由液面形状、圆柱压强分布和受力时间历程与单圆柱入水的结果进行了比较.     

潜艇PMM实验的CFD仿真技术研究

操纵性是现代潜艇最重要的总体性能之一.为预报潜艇的操纵性能,需求得所有的水动力导数.借助CFD技术和动网格技术数值模拟小振幅平面运动机构试验,求解了SUBOFF主艇体的水动力导数,并与混合分布法计算的附加质量项作了比较,二者符合良好,从而验证了方法的有效性.流场的求解基于RANS方程和RNG k-ε湍流模型,动网格的处理采用Hrvoje Jasak和Zeljko Tukovic的方法,压力和速度的耦合采用PIMPLE方法解耦.流场的控制方程采用有限体积法离散且分离式求解,网格运动的控制方程采用有限元方法离散,网格的分裂采用体和面同时分裂法.     

自适应网格下的近岸波浪缓坡方程数值模型

在近岸变水深缓坡海域,充分考虑到波浪传播过程中的折射、绕射、反射以及波浪破碎等效应,采用引入破碎项的椭圆型缓坡方程,基于自适应四叉树网格建立了近岸波浪缓坡方程的数值模型,数值模型采用有限体积法求解。该网格可根据计算节点处波长与网格尺寸间的关系自动判定是否进行网格加密划分,使之适应椭圆型缓坡方程每个波长范围内至少布置8个计算节点的特征,提高计算精度的同时又能提高计算效率。采用已有的物理模型试验和实测资料对数值模型进行了验证。结果表明,该模型可以有效地模拟近岸变水深缓坡海域波浪传播。     

二维泥沙数学模型的改进——模型的建立

本文针对河道弯曲、狭长的特点,引进贴体坐标.将复杂的河道计算域转换成规则矩形域,使问题在规则域上求解.拟合坐标系下模型改进表现在以下两方面,一是引进紊动动能和紊动耗散率闭合雷诺方程,使模型能够回答整治工程局部水流结构,如回流现象,二是建立了一种新的悬沙输移模式,该泥沙模型能体现悬沙与床沙交换机理,既适用于水库上游淤积计算,也适用于水库下游冲刷计算.     

基于体积力方法的船舶波浪回转运动数值仿真

研究船舶在波浪环境中的自由操纵运动对船舶操纵性具有重要意义.本文对ONRT全附体模型在规则波中的回转运动进行仿真研究.采用基于结构化动态重叠网格的自研CFD软件HUST-Ship求解RANS方程与船体六自由度运动方程的耦合问题.RANS方程通过有限差法离散,并使用PISO算法进行求解.波浪中操纵运动使用了移动计算域方法...     

深槽水域垂向二维水沙运动数值模拟

采用网格堵塞技术实现实际地形的模拟,采用有限体积法离散方程,采用SIMPLER算法求解耦合的速度场和压力场,并建立了沿宽度积分的垂向二维水流泥沙数学模型。将计算结果和试验资料进行对比,结果吻合良好。     

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

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

回转体三维绕流场数值计算

提出一种满足主尺度和均衡条件的水滴形潜艇回转体数学描述方法，借助采用k-ε湍流模型和有限体积法来求解RANS方程的方法，对用该数学描述方法生成的一种回转体模型进行了数值模拟计算；采用两层壁面模型，以便在近壁区域采用精细网格来模拟回转体边界层内的流场流动特征。同时，对文献[1]给出的一种回转体模型也进行了数值模拟计算。给出了回转体表面压力系数纵向分布，计算了多种来流速度时不同横剖面处的速度以及表面摩擦阻力系数分布曲线，并与有关文献给出的实测值和计算值进行了比较，证明其分布趋势基本一致，而且与实测值更加接近，从而为水滴形潜艇流场计算和阻力预报提供了一种有效的计算方法。     

巡飞器气动特性研究

巡飞器是无人机技术和弹药技术有机结合的产物。采用气动估算和数值计算相结合的方法,研究巡飞器的气动特性。应用Digital Datcom获得气动估算解。采用分区网格技术和Euler方程计算巡飞器绕流。分析了适用于巡飞器复杂外形的网格拓扑,给出了巡飞器分区网格。利用有限体积法求解三维Euler方程。结果表明气动估算与数值计算满足实际工程需要。     

无轴轮缘推进器中可拆卸螺旋桨的性能分析

应用FLUENT数值模拟软件对无轴轮缘推进器螺旋桨进行水动力分析,根据实际运行工况需求确定数值模拟的基本方程、湍流模型、初始条件、边界条件、流体域以及网格划分等参数,建立相关数值模拟模型,计算出不同桨叶数目的螺旋桨在一定进速条件下的桨叶压力分布状况及螺旋桨推力与扭矩,并与常用经验公式计算结果进行比较,验证了无轴轮缘推进器螺旋桨有限体积元模型的正确性。研究结果可为该类型螺旋桨的优化设计提供理论支持。     

带附体潜艇绕流场数值计算与验证

采用κ-ε湍流模型和有限体积法用求解RANS方程的方法,对全附体的SUBOFF AFF-8潜艇模型三维粘性绕流场进行了数值模拟计算.一种结合增强型壁面函数的两层模型被采用,以便在近壁区域采用精细网格来模拟边界层内的流场流动特征.给出了雷诺数Re=1.2×107下模型表面纵中剖面线上半部分压力系数与摩擦阻力系数的分布,以及指挥台围壳与尾翼表面不同高度处的压力系数分布,同时计算了桨盘面处的尾流速度分布,并与有关文献给出的实测值进行了比较和验证,分布趋势基本一致,较相关文献给出的计算结果有相当大的改进,此外从指挥台围壳与尾翼后的速度矢量图观察到了与有关文献相同的涡旋流场结构,从而验证了该数值方法的有效性,为潜艇概念设计中评价设计方案优劣提供了重要的技术手段.     

基于两种计算方法的二维水翼性能预报

对拉普拉斯方程采用Galerkin方法得到积分表达式,建立以流场速度势为未知量有限元总体方程.并采用2种不同的计算模式求解有限元方程,即控制流线法和速度势分解法,计算了绕NACA4412二维翼的表面压力分布,不同攻角的升力系数以及流线分布.讨论了2种方法各自的特点.计算结果与试验数据进行了对比,取得了较为一致的结果,说明了本文方法的有效性和可靠性.     

基于CFD的螺旋桨定常水动力性能预报精度研究

根据螺旋桨的型值参数,在Fluent前处理器ICEM中建立2套不同的网格模型对螺旋桨的敞水性能进行预报,采用计算流体力学(CFD)理论,结合雷诺平均纳维—斯托克斯(RANS)方程和3种不同湍流模型对螺旋桨的敞水性能进行研究,模拟在不同进速系数下的螺旋桨的推力、转矩系数、桨叶表面压力分布以及桨后尾流场的情况等。通过与试验数据的比较表明:采用RANS方法结合不同网格及湍流模型都能准确预报螺旋桨的敞水性能,其中六面体网格模型计算结果及收敛速度更优,且结合雷诺应力模型(RSM)精确度最高,更适合螺旋桨粘性流场的数值计算。     

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

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

船舶螺旋桨周围粘流场的数值模拟方法

本文叙述了通过直接求解雷诺平均力方程亚获得螺旋桨粘流场数值解的方法与数值求解步骤，该方法采用非交错网格系统，利用幂函数格式离散动量方程，预报一校正方法求解速度一压力藕合问题，应用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍 模式求解雷诺应力项使方程组支此来获得螺旋桨粘流场。     

航槽三维流动的数学模型

本文提出了一种模拟挖槽内三维流动的数学模型.该模型将压力分解为动水压力与静水压力之和.由控制体积法导出三维偏微分方程的离散格式.采用Patankar和Spalding提出的压力校正法求解动水压力场.然后求解水位控制方程.用最新求出的水位替代压力校正法中的水位.如此循环迭代直至收敛.该模型无须采用“静水压力假定”和“刚盖假定”.能够模拟不可压缩流体含自由表面流动问题.利用该模型计算得到的挖槽内流速分布与水槽试验值吻合.此外,本文还计算了相当于天然情况各种水深、各种挖深比、各种挖槽底宽下挖槽内的流场.     

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

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

应用分区有限体积法数值计算舵粘性水动力

通过数值求解不可压缩粘性流的RANS方程计算均匀来流中的船舵水动力。采用标准的κ-ε湍流模式封闭方程，并用二阶有限体积法在非交错分区贴体网格上对守恒方程进行离散；采用Rhie和Chow的插值方法确定单元面的质量通量，用SIMPLE法求解压力分布，用壁函数法处理壁面边界条件。对一种NACA冀型舵在两种不同雷诺数下大攻角范围的流动进行了计算，得到的舵水动力与试验结果和他人的计算结果进行比较，在舵未达到失速角前吻合良好，而所确定的舵失速角与最大升力值则均偏小。失速角后舵背流面的流动分离点也能捕捉到。最后分析了雷诺数对舵水动力的影响。