基于CFX的全垫升气垫船船身空气动力模拟

基于Solidworks三维实体建模软件建立了全垫升气垫船三维模型,利用ANSYS CFX计算流体力学软件对全垫升气垫船模型的流场进行数值模拟,计算得到了不同攻角和漂角下相应的空气动力系数;讨论了攻角及漂角的变化对空气动力系数的影响.数值结果表明:随着漂角的增大,横向力也增加,漂角达到90°时横向力最大;随着漂角的增大,纵向阻力降低,漂角为90°时达到0°;漂角继续增加时,纵向阻力为正值,横向力、垂向力增大显著,进而严重影响了全垫升气垫船航行的稳定性.将模拟计算结果与试验数据进行对比,发现计算值与实验值符合较好,验证了CFX软件在全垫升气垫船空气动力学模拟中的适用性和可靠性.     

KVLCC2船-舵系统斜航水动力数值计算

通过求解RANS方程,选取RNGk-ε湍流模型,在模型尺度下对KVLCC2船-舵系统处于不同漂角时的斜航黏性流场进行数值模拟.首先采用3套网格进行网格收敛性分析,然后通过计算结果与试验数据比较,验证数值计算方法的可靠性.在此基础上,计算了舵角为零时不同漂角下以及不同舵角、不同漂角下作用在船-舵系统上的横向力和转首力矩,并对计算结果进行了分析.     

不同头型回转体低速倾斜入水过程流场特性数值模拟

基于有限体积法和VOF多相流模型求解气、水两相流动的RANS方程,结合动网格技术,对回转体低速倾斜入水过程进行数值模拟研究.通过将数值计算结果与试验数据进行比较,验证了数值计算方法的有效性.基于该方法对不同头型回转体低速倾斜入水过程进行分析,得到不同头型条件下回转体入水空泡形态发展规律、回转体运动特性及流体动力特性变化规律.研究结果表明:同一入水深度,入水空泡直径和流场中最小压力值随着回转体头部锥角的增加而增大;不同头型回转体,锥头母线同一位置处压力随着头部锥角的增加而增大;回转体阻力系数与其头部锥角大小直接相关,锥角较大时,阻力系数也较大,速度衰减也较快.      

叶片包角对离心泵流场及脉动特性的影响

针对离心泵非定常流动压力脉动特性,采用滑移网格的大涡模拟技术对叶片包角分别为95°,100°,105°,108°的4副叶轮进行数值模拟.分析了叶片包角对离心泵水力性能、叶轮出口"射流-尾迹"、测点压力脉动频谱特性和叶轮径向力的影响关系.结果表明:随着包角的增大,离心泵的水力性能下降;包角适当增大,会使叶轮射流-尾迹流动结构变弱.在设计工况下,蜗舌附近测点压力脉动最大;在蜗壳螺旋段压力脉动强度沿流动方向逐渐变弱,而在叶轮流道内压力脉动沿流动方向逐渐增强,在叶轮出口处达到最大;而离心泵叶轮所受径向力随着包角的增大而减小,适当地增大包角可以提高离心泵运行的可靠性.     

导管式调距桨水动力性能的数值预报分析

采用全结构化网格,利用稳态多参考系(MRF)方法,对普通调距桨和导管式调距桨变螺距下的水动力性能进行了有效预报.首先对普通调距桨水动力性能进行了计算,计算结果与试验数据符合良好.结果表明:进速系数较小时,较小螺距角才能使螺旋桨的效率达到最大,而进速系数较大时,较大螺距角才能使螺旋桨的效率达到最大;通过分析桨叶表面压力变化,发现螺距角变化较大时,叶表面压力变化不均匀,可能导致噪声性能恶化.随后计算了导管式调距桨的水动力性能,分析得出:随着螺距角减小,导管推力也减小,且大进速小螺距角时,螺旋桨效率下降迅速.最后分析了导管式调距桨的尾流场,得出了尾流场速度随螺距角变化的一般规律.     

CLT桨的尾流场及梢涡特性数值分析

为了探究尾流收缩叶梢负载(CLT)桨的尾流场及梢涡特性,基于商业计算流体力学软件Star-CCM+对CLT桨P1727的尾流场进行了数值模拟.计算采用DDES(延迟分离涡)方法在进速系数J=0.5的工况下对三套网格进行了不确定度分析,确认了数值方法的可行性,然后对螺旋桨的敞水性能及尾流场进行了计算.研究表明:尾流场的网格细化对螺旋桨的推力及扭矩只有微弱影响,对桨后的流场具有较大影响;尾流区域轴向速度分布分为加速流与自由流,进速系数越大,加速流与自由流的分界线向外扩散的趋势越弱;相邻梢涡之间会发生融合,进速系数越大,融合得越晚,梢涡强度越弱;CLT桨相较于常规桨多出一个端板涡,并且会与梢涡融合;梢涡在向下游发展的过程中,会有明显的收缩现象.     

沉深变化对螺旋桨水动力性能影响的数值分析

为了分析桨轴沉深对螺旋桨水动力性能的影响,采用VOF方法对某型桨在静水状态下、不同沉深时的水动力性能进行了研究.采用混合网格技术.根据沉深的不同选用不同的计算域模型来进行网格的划分.通过静水中某沉深下螺旋桨推力系数、转矩系数以及效率的计算结果与实验值的对比,验证了研究结果的正确性.由结果分析得知:随着进速系数的增大,桨轴沉深对螺旋桨推力系数、转矩系数的影响逐渐减小,而对效率的影响逐渐增大,且螺旋桨对自由液面的抽吸位置逐渐后移,扰动越来越小.     

风向角对高速列车驶出隧道过程中气动效应的影响

为研究风向角对驶出隧道过程中高速列车气动效应的影响,以某型高速动车组列车为研究对象,采用数值模拟方法对隧道内气动压力、列车风风速、流场分布及列车气动荷载进行分析。通过与动模型试验结果进出对比,验证数值模拟方法的准确性。研究结果表明：隧道壁面气动压力峰值及变化幅值最大值出现在隧道内部,且出现位置到隧道出口距离与风向角有关;背风侧气动压力受风向角影响更大,气动压力变化幅值随风向角增大呈现先减小后增大再减小的趋势;出口处列车风风速随风向角增大基本呈现先增大后减小的趋势,30°风向角时列车风风速最大,但迎、背风侧列车风风速峰值出现时刻不同;随着风向角增大,流场分布不对称性增强,列车绕流特性由流线型绕流逐渐过渡到钝体绕流,流动分离点到头车鼻尖的距离呈现先增大后减小最后再增大的变化规律,隧道内流动结构愈加复杂;气动横向力、升力变化幅值随风向角增加呈现先增后减趋势,头车横向力系数最大变化幅值分别是中车、尾车的2.4倍和2.6倍,升力系数最大变化幅值分别是中车、尾车的1.1倍和1.5倍,故保证头车安全是控制整车运行安全的关键;侧风下高速列车驶出隧道情形下的最不利风向角为30°,此时头车发生列车事故风险...     

基于CFD的对转桨扇发动机性能确定技术研究

以某对转桨扇发动机三维模型为对象,采用CFD方法对桨扇滑流开展数值模拟。利用分区拼接网格技术对桨扇旋转区域和外流场域进行网格划分及拼接;在此基础上采用雷诺平均N-S方程,雷诺应力项采用RNG 湍流模型,基于滑移网格方法,开展了针对不同来流Ma、桨叶角及桨叶转速的流场的数值计算。结果表明：桨叶角对于拉力影响较大,在桨叶角30°变化到35°,桨扇的拉力变化范围可达37.8%,功率变化可达28.4%。桨叶拉力随着转速的增大而增大,在转速达到1600rpm后,前后排桨叶的拉力系数最大相差33%,当转速继续增大,气流通过前后排桨叶偏转角度大,气流失速严重,桨叶拉力增幅减小。本文通过对该型桨扇发动机流场的计算,定量获取了桨扇的拉力、功率特性,为后续桨扇发动机的推力计算提供参考。     

倒车时螺旋桨横向力计算方法

介绍了关于倒车螺旋桨横向力计算的问题.通过对倒车螺旋桨状态的分析,提出一种螺旋桨横向力计算的理论方法.类比进车时排出流对舵影响的方法,推导出倒车排出流的计算公式,同时对斜流横向力的系数进行了修正.依据理论分析,分别对不同的船型进行仿真研究,最后以实船进行验证.仿真结果表明,计算方法是合理有效的.     

自适应四叉树网格下的N-S方程数值求解模型

提出了一种自适应四叉树网格下的N-S方程数值求解模型.网格能够根据涡度值大小进行自动加密或合并,以达到在不显著增加计算量的前提下,提高重点区域分辨率的目的.模型中采用了无条件稳定的MacCormack格式计算对流项,采用修正的中心差分格式离散压力泊松方程,并提出了在树型网格下黏性项的变通离散格式.通过算例证明,利用新模型所得到的压力泊松方程的数值解具有二阶精度,速度解的精度超过一阶.计算得到的方腔流中轴线上速度分布与Ghia计算结果一致,圆柱绕流中拖曳力系数和升力系数与实测结果一致.方腔流算例还表明,在相同分辨率情况下采用自适应网格计算时间可减少近一半.     

二维物体入水数值模拟

基于势流理论和边界元法模拟了二维剖面以恒定速度入水问题,采用双点模型处理自由面物面交点、4阶龙格-库塔法对自由面进行迭代求解,并利用自由面平滑技术和网格重新划分技术阻止自由面不稳定现象,根据射流区截断模型处理射流区,采用辅助函数法求解物面压强分布以保证数值解的准确性和稳定性.探讨了网格划分对数值结果的影响以及剖面斜升角对最大压强值以及垂向力的影响.预报了不同斜升角的V型剖面的自由面形状和剖面压强分布、最大压强值以及垂向力,并与相似解进行了比较.比较结果验证了本文方法的有效性和准确性.     

凸轮泵内部瞬态流场的动网格数值解析

根据凸轮泵啮合特性建立了渐开线方程,当型线方程的压力角介于40°~50°时,共轭转子相对运动时不产生干涉效应.为了阐明压力角对转子型线方程和几何参数的影响,基于FLUENT动网格模型和局部网格重构的网格变形技术,应用RNG k-ε湍流模型及PISO算法对凸轮泵内部流场进行数值解析.通过定量监测排出端流量波动和吸入端压力脉动,获得不同压力角和不同工况下,流量波动和压力脉动的瞬态特性.结果表明:凸轮泵排出端流量和吸入端压力呈周期性变化,当转子压力角逐渐增大时,二者有减小的趋势;转速是影响流量波动和压力脉动的主要因素,随着转速升高,排出端流量和吸入端压力提升明显,脉动现象明显;当转速升高到一定程度后,排出端速度与吸入端压力的脉动值增加较少.     

复合材料螺旋桨水动力性能与结构响应数值研究

通过CFD计算复合材料螺旋桨水动力载荷,应用有限元法计算复合材料桨结构响应,并基于双向耦合算法实现复合材料螺旋桨在均匀来流下的流固耦合数值模拟.研究了不同进速系数和不同铺层角度工况复合材料螺旋桨的水动力性能及结构响应.结果表明,在进速系数J ≤ 0.8时,复合材料螺旋桨的推进效率明显高于同工况下的金属螺旋桨.随进速系数增加,复合材料螺旋桨推进效率先增加后减小,在进速系数J = 0.8时取得最大值.桨叶总变形和等效应力分布与铺层角度有较大关系.复合材料螺旋桨的螺距角较金属螺旋桨减小,减小的螺距角与桨叶攻角的变化相匹配,自适应地提高了该螺旋桨的推进效率.      

拖式与推式吊舱推进器水动力性能数值仿真及对比

为研究吊舱推进器在直航和回转工况下的推进性能,通过RANS方法结合标准k-ε湍流模型对拖式与推式吊舱推进器在不同工况下的水动力性能进行对比分析.分别计算拖式与推式吊舱推进器直航与回转工况下的推力系数与转矩系数,并与试验结果进行对比以验证数值计算的准确性.直航工况推进特性曲线、回转工况的推力系数与转矩系数曲线对比结果表明,数值计算可以准确预报吊舱推进器的水动力性能.根据该模型计算了斜流工况时两种推进器的推力、转矩系数以及轴向力和侧向力系数随角度的变化规律.结果表明,本文提出的水动力模型可以准确预报不同工况下吊舱推进器的水动力性能.     

不同叶顶间隙对泵喷推进器性能的影响

基于均质多相流的RANS方程与分块网格技术,结合剪切应力SSTk-ω湍流模型,对E779A桨进行了非空化与空化数值模拟,数值计算结果与实验值符合较好,验证了数值计算方法的合理性与有效性.在此基础上,对不同叶顶间隙的泵喷推进器进行了非空化与空化数值模拟计算,重点对叶顶间隙处流场进行了精细后处理.结果表明:未空化时,间隙尺寸越大,泵喷推进器整体效率值越低,最大间隙与最小间隙之间的效率降幅最大为8.62%,叶顶间隙尺寸直接对泵喷推进器的最优进速系数产生影响.当空化发生时,在最优进速系数下的推进器效率下降幅度最小;间隙尺寸越大,空化后的推进器整体效率越低;在较高转速下,较大间隙的空化面积在叶稍处有更为明显的激增量.     

不同攻角下水轮机导叶流道湍流特性研究

采用κ-ε湍流模型对不同导叶攻角下的导叶流道进行三维数值模拟,采用SIMPLEC算法实现速度、压力变量的分离求解。计算结果表明,导叶压力面的压力系数随着攻角增大而明显增大,而吸力面的压力系数随着攻角增大而明显减小。随着攻角增大,导叶边界层提前分离,形成脱落涡,导叶的升力明显增加。在导叶出水边的尾迹涡,其湍动能很大,随着流动进入转轮流道,引起水轮机水力振动。     

导管结构改进后导管桨的水动力性能研究

对19A导管桨的试验模型进行改进,在导管内壁后半部分增加突起结构来提高导管桨的推进效率.应用计算流体力学中的SSTk-ω模型对19A导管Ka4-55螺旋桨的水动力性能进行计算,得到了导管桨在各个进速系数下的推力系数、转矩系数和螺旋桨效率以及导管上的推力系数.通过与实验结果进行对比,验证计算方法的可靠性.对导管改进结构参数进行优选,得到设计进速系数下效率较高的导管模型.对比结果显示:改进导管可以提高导管桨效率,设计进速系数J=0.5处,效率能够提高2.62%.通过对比改进前后导管桨压力分布和导管内速度分布,分析了改进导管桨水动力性能提高的机理.     

转速对贯流泵装置流道水力参数影响的数值分析

为了研究贯流泵装置进出水流道水力性能受转速的影响规律,以贯流泵装置为研究对象进行全流道的数值计算,通过将能量性能数值预测结果与物理模型试验结果对比,验证了数值模拟结果的有效性。在考虑水泵与进出水流道内流相互影响条件下定量分析转速对进出水流道水力性能参数的影响规律。结果表明:(a)不同转速时进水流道水力损失比与流量系数的关系曲线变化趋势基本相同。相同转速时,随着流量系数的增大,进水流道水力效率逐渐减小,流道出口环量先减小、后增大;相同流量系数时,进水流道出口面的平均环量随转速的增加而增加。(b)不同转速时,随着流量系数的增大,出水流道进口入流涡角相对值先减小、后增大。在相同流量系数时,随着转速的增加,出水流道的静圧比逐渐增加,出水流道进口面的偏流角分布基本相同。     

运输机升阻特性仿真与风洞试验相关性方法

以M3飞机配装某型发动机三维模型为研究对象,采用CFD方法对整机流场进行数值计算。对于带动力短舱模型,利用分区拼接网格技术对发动机内流场和飞机外流场进行网格划分和拼接;在此基础上采用雷诺平均N-S方程,基于S-A湍流模型,开展了不同发动机状态、马赫数及攻角的仿真计算。以FNPR=1.61时,试验获得的升阻系数作为基准,在不同攻角下,获得CFD计算结果的修正因子,结果表明：修正后的数值计算结果与风洞试验获取的升力特性曲线,贴合程度好,在攻角小于14°内误差小于3%;修正后的阻力特性曲线整体趋势与风洞试验一致,误差小于10%,阻力系数都是随攻角的增大而增大,且在攻角大于10°后快速增大。