基于修正RNG k-ε模型的叶片泵非设计工况数值模拟

对于非设计工况下流体机械的流场计算,由于其偏离设计工况,分离流动严重等特点,使得用普通的湍流模型数值模拟精度不高。本文采用一种将涡粘性系数Cμ与湍流脉动动能和湍流耗散率的变化相关联的方法,对RNG k-ε湍流模型进行修正。为验证其对非设计工况下叶片泵的预测精度,分别采用修正模型和原始RNG模型,数值模拟了一比转速为439的叶片泵的内部流场,得到了其扬程、效率曲线及在典型工况下的内部流场流线分布。通过比较改进前后RNG湍流模型的数值计算结果,并与实验结果对比,得出修正RNG模型更能准确预测出非设计工况下叶片泵的流场特性。     

弱可压缩移动粒子半隐式方法（WC-MPS）与湍流大涡模型（LES）模拟溃坝流动

应用弱可压缩移动粒子半隐式法（WC-MPS）结合湍流大涡模型(LES)模拟了自由液面流动问题。该方法通过流体状态方程实现流场压力的直接计算,并根据大涡模拟Smagorinsky模型来计算湍流项。最后应用该方法,针对两个溃坝自由液面流动问题,特别考虑不同粒子直径以及湍流Smagorinsky系数,进行数值对比研究以及实验对比分析,得到了采用较小粒子直径和较大湍流Smagorinsky系数模拟大变形自由液面具有较好结果的结论。     

浮体吃水比对窄缝流体运动特性的影响

本文结合粘性流体理论和VOF模型,采用Waves2Foam开展了浮体吃水比对浮体间窄缝流体运动特性影响的研究。结果表明,同吃水浮体间窄缝流体振幅分布与实验结果一致,数值水池的水深对数值结果有影响;同时发现,浮体吃水比对多浮体系统窄缝流体共振的波浪幅度和共振频率有非线性影响。浮体吃水增加,共振频率下降,迎浪侧浮体吃水增加,导致反射波系数增加,流体振荡波幅显著下降,背浪侧浮体吃水增加,导致反射波系数下降,流体振荡波幅上升。同时,由流体粘性引起的粘性耗散随着窄缝流体振荡幅度的增加而同步增加。最后,对涡量演化的分析说明了窄缝流体振荡幅度的变化。     

基于SPS-SPH方法的自由表面流动数值模拟

SPS-SPH方法是为了克服经典SPH方法中人工粘性过于分散的缺点,以及克服改进的SPH方法中层流粘性不能准确捕捉湍流流体运动的不足。该方法借鉴大涡模拟方法的思想提出了一个亚粒子尺度(Sub-Particle Scale,简称SPS)涡粘性项并将它附加在Morris的层流粘性项之后。分别采用SPS-SPH方法和经典SPH方法对2D溃坝问题进行数值模拟并与实验数据进行了比较,结果表明,SPS-SPH方法的计算结果比经典SPH方法更加接近实验结果,该方法合理可靠,能够应用于数值模拟自由表面流动的相关问题。将SPS-SPH方法进一步运用于求解二维和三维典型自由表面流动问题,计算结果表明,SPS-SPH方法对于自由表面流动的数值模拟问题具有较高的精确度,对解决此类问题有较大的参考价值。     

均匀流中潜艇水下运动表面尾迹的数值模拟

本文通过基于RANS方程和VOF方法的非定常粘性数值方法对潜艇在均匀流中水下运动的表面尾迹特征进行了探索性研究,较好地解决了用粘性流方法来模拟远场波浪中的数值耗散问题,数值计算得到的潜艇运动尾迹特征与诸多文献和理论分析结果基本一致.本文均匀流中的数值模拟将成为潜艇在分层流运动数值模拟的基础.     

用基于SST模型的DES方法数值模拟圆柱绕流

分离涡数值模拟方法(DES)在物面附近求解雷诺平均Navier-Stokes方程,在其他区域采用大涡模拟方法.兼具前者计算量小的优点和后者能模拟大分离湍流流动、计算精度高的优势.采用基于剪切应力传输(SST)两方程湍流模型的DES方法对粘性不可压缩流体的圆柱绕流问题进行数值模拟.通过对所得的速度场、压力场、阻力系数Cd、升力系数Cl、斯特罗哈数St等结果的分析及与文献上的实验和计算数据的比较,说明DES方法对于层流及高亚临界雷诺数的圆柱绕流模拟是合理的.     

模拟斜坡堤上越浪量一种新的数值模式

利用FLUENT软件,基于粘性不可压缩流体的雷诺时均N-S方程和k-ε湍流模型,采用有限体积法对控制方程进行离散,建立数学模型。基于质量源造波和动量源消波的方法,在连续方程和动量方程中加入造波源项和消波源项,采用VOF方法捕捉自由表面,开发了无反射造波的二维数值波浪水槽,通过数值模拟值和理论值的对比,对该数值波浪水槽进行了验证。进而,利用该数值波浪水槽模拟斜坡堤的越浪问题,分析了平台宽度、斜坡坡度对越浪量的影响,并将数值模拟结果和前人物模试验结果相对比,二者趋势一致,吻合良好,表明了该模型在处理斜坡堤越浪问题方面的有效性和精确性,从而建立了利用FLUENT模拟斜坡堤越浪问题一种新的数值模式。     

基于RANS方法的导管推进器梢隙流动数值模拟北大核心CSCD

以导管推进器为研究对象,采用雷诺时均纳维斯托克斯(RANS)方法对其梢隙流动进行数值模拟研究。通过对网格类型、湍流模型的适用性研究以及对梢部流场的研究探讨,初步建立了基于RANS梢隙流动的数值模拟方法;数值计算结果与实验结果吻合较好。对比分析发现,结构化网格与非结构化网格相比能捕捉到梢隙流动中更加细节的流场信息,如壁面边界层流动等,更适合于梢隙流动的数值模拟。3种湍流模型SST?k-ω,RNG?k-ε及RSM的计算结果基本一致,都能有效模拟梢隙流动。通过间隙区域流场分析发现,梢隙流动的驱动力主要是叶面与叶背之间的压差,受壁面边界层流动的影响。流体进入间隙时流动分离形成间隙分离涡,间隙泄漏流穿过间隙与吸力面侧流体相互作用,卷起形成梢隙涡,在约37.5%弦长位置处形成并附着在桨叶壁面发展,大约在75%弦长位置与桨叶分离进入尾流场中。研究获得了梢隙涡的起始、发展、脱落的变化过程。     

平板上直立圆柱三维绕流的CFD计算研究

利用计算流体动力学(CFD)方法,对粘性不可压缩流体流过平板上直立圆柱的绕流进行了三维数值模拟,采用了有限体积法和SIMPLE计算程式,利用不可压缩的Navier-Stokes方程,模拟了雷诺数在200的层流绕流流动。在建立贴体的结构化精细网格之后,成功模拟了马蹄涡,与经典的理论对比具有较好的吻合性。同时得到了交接部区域的剪切力云图。计算结果表明,主/附体结构交接部区域流动机理十分复杂,在这个区域有马蹄涡的存在,是产生交接部冲蚀的根本原因。计算结果表明本文所用方法在模拟含有复杂流动机理的圆柱绕流中是可靠性的。     

基于RANS方法的导管推进器梢隙流动数值模拟

以导管推进器为研究对象,采用雷诺时均纳维斯托克斯(RANS)方法对其梢隙流动进行数值模拟研究。通过对网格类型、湍流模型的适用性研究以及对梢部流场的研究探讨,初步建立了基于RANS梢隙流动的数值模拟方法;数值计算结果与实验结果吻合较好。对比分析发现,结构化网格与非结构化网格相比能捕捉到梢隙流动中更加细节的流场信息,如壁面边界层流动等,更适合于梢隙流动的数值模拟。3种湍流模型SST?k-ω,RNG?k-ε及RSM的计算结果基本一致,都能有效模拟梢隙流动。通过间隙区域流场分析发现,梢隙流动的驱动力主要是叶面与叶背之间的压差,受壁面边界层流动的影响。流体进入间隙时流动分离形成间隙分离涡,间隙泄漏流穿过间隙与吸力面侧流体相互作用,卷起形成梢隙涡,在约37.5%弦长位置处形成并附着在桨叶壁面发展,大约在75%弦长位置与桨叶分离进入尾流场中。研究获得了梢隙涡的起始、发展、脱落的变化过程。