轴流压气机尾流撞击效应机理的数值模拟

通过精确给定转子进口边界条件,以模拟上游静子尾迹对转子流场的影响,这样利用单排的计算结果模拟转/静干涉引起的非定常流动对流场气动参数的影响,因为相邻叶片排的叶片数不同,从而避免传统方法中常用的叶片数约化存在的难以正确模拟非定常激励频率的问题,对计算时间和计算能力的要求也可减少.针对一个低速转子的数值模拟结果显示,如果能合理地组织来流尾迹的扫过频率和尾迹亏损的幅度,对转子气动性能的改善是很显著的,从而为风扇/压气机中控制和利用非定常流动提供了一个可行的方向.     

喷水推进泵临界空化工况空化流态试验

在舰艇空间尺寸的限制下,提高推进器的综合性能是舰船设计的关键。空化是喷水推进泵的主要噪声源和重要设计参数。为降低舰艇的噪声,提高喷水推进泵的空化性能,十分有必要研究空化流动结构演变规律及其对喷水推进泵性能的影响。以某型喷水推进泵为研究对象,采用高速摄像技术开展了喷水推进泵临界空化工况空化流动结构演变规律试验研究。在空化发生和发展过程中,喷水推进泵空化流动结构包括片状空化、云状空化、叶顶间隙空化、叶顶泄漏涡空化和垂直空化涡。试验捕捉了空化演变的物理过程,分析了各空化流动结构对喷水推进泵性能下降的影响。结合数值模拟和前人研究,阐述了空化流动结构形成机理及其对喷水推进泵性能的影响。研究结果为喷水推进泵内的空化现象提供了新的认识,也为空化性能预测方法的研究给予了借鉴和指导。      

用GAO-YONG湍流模式对分离流动的数值模拟

采用GAO-YONG湍流模式对二维管道内三角形钝体绕流问题以及槽道内后台阶分离流动进行了数值模拟.求解程序基于开源数值计算平台OpenFOAM (Open Field Operation and Manipulation),数值模拟结果很好地预测了钝体绕流问题中漩涡脱落的尾流的流动趋势以及后台阶流动中分离再附的流动结构,同时分析了速度分布以及摩擦系数等参数并与实验值进行了对比,结果符合很好.这表明GAO-YONG湍流模式对大分离流动有较好的预测能力,对工程实践具有指导性作用.     

用GAO-YONG湍流模式对分离流动的数值模拟

采用GAO-YONG湍流模式对二维管道内三角形钝体绕流问题以及槽道内后台阶分离流动进行了数值模拟.求解程序基于开源数值计算平台OpenFOAM (Open Field Operation and Manipulation),数值模拟结果很好地预测了钝体绕流问题中漩涡脱落的尾流的流动趋势以及后台阶流动中分离再附的流动结构,同时分析了速度分布以及摩擦系数等参数并与实验值进行了对比,结果符合很好.这表明GAO-YONG湍流模式对大分离流动有较好的预测能力,对工程实践具有指导性作用.     

超空泡水下航行体的结构动力响应特性

通过综合有关文献的理论,得到作用于超空泡水下航行体尾部的冲击载荷与航行速度的关系,然后利用有限元法研究了在不同速度条件下航行体的结构动力响应和变形特性.得到了加速度响应的主频率及其取值区间、航行体内的最大应变值及其分布位置;并给出了加速度响应主频率以及航行体内最大应变值与运动速度的关系.计算结果表明:存在2个加速度响应主频率,当航行体航速相对较低时,应采用第1个主频率作为结构设计时的参考频率,速度较高时应采用第2个主频率;计算还表明:航行体内最大应变值随着航速的增加而显著增加.结果对超空泡水下航行体的结构强度分析和结构优化设计有着指导意义.     

内孔结构对旋转空化发生器水力特性的影响

空化普遍存在于工业生产、舰船推进、航空航天等领域.使用高速相机对闭式循环试验台的旋转空化发生器内部流动进行试验观测，然后基于RNG k ε湍流模型和Zwart Gerber Belamri空化模型开展了旋转空化发生器流场的数值模拟.结果表明数值模拟结果与试验观测数据吻合较好，验证了本文数值模拟方法的准确性和可靠性.受离心力作用，转子孔内的压力从内径到外径逐渐增大.内孔底部由于轴向旋涡的作用，产生较大旋涡并耗散能量.内孔顶部流体与腔体区主流相互作用，在内孔顶侧产生碰撞并形成较小漩涡.不同内孔结构对旋转空化发生器的空化效果具有重要影响，在相同工况下椭圆形内孔的空化率小于圆锥形和圆柱形，说明型线光滑的内孔结构局部损失小，产生的空化较弱.不同内孔结构下，旋转空化发生器内部压力脉动的主频为转频fi或24fi(对应圆周方向开孔排数)，压力脉动最大幅值出现在进口管侧和出口管侧，主要原因是受动静干涉作用影响.     

超声速场中的反向喷流数值模拟

采用了高分辨率的N-S(Navier-Stockes) 方程数值模拟方法,对钝头体头部反向喷流进行研究.计算并分析了头部反向喷流现象中喷流马赫数、来流马赫数、攻角等因素对流场细致结构的影响,并对反向喷流减少阻力和减少气动加热的原理进行了深入分析和探讨.研究分析表明这些因素通过决定弓形激波、马赫盘的强度和位置,来影响回流区的大小和位置,从而使有反向喷流时阻力系数所受到的影响比无喷流时更大.      

超空泡鱼雷外形优化

为了优化超空泡鱼雷的外形参数,获得同一空化数下的最大航行体体积,基于经验公式,考虑到空泡与超空泡雷体的关系,通过对圆盘空化器产生的超空泡形态近似处理和对超空泡鱼雷外形分段分析,推导了航行体最优外形下的半径计算公式.计算了一定空化数下航行体优化后的体积,与已有的某超空泡航行体原型做了对比,验证了优化结果的有效性.结果表明,外形优化能够增加航行体的容积和航程,提高超空泡鱼雷的性能.这为提高超空泡鱼雷的航程提供了一个方向,对今后超空泡航行体外形的优化设计提供参考.      

再入飞行体远场二维湍流尾迹的工程计算方法

建立了再入体再入大气层形成的尾迹的工程模型,并运用边界层理论求解二维等压湍流尾迹流动。根据Ｐｒａｎｄｔｌ的湍流切应力的混合长度理论和Ｔａｙ－ｌｏｒ的涡量交换理论分别获得速度分布和温度分布,根据Ｌｙｋｏｕｄｉｓ的远场尾迹理论计算了轴心参数分布,并且获得了尾迹长度、自由电子数密度等尾迹参数。计算结果与同类文献一致。文章的目的是为计算尾迹的雷达散射截面提供初始参数。     

基于流动控制的水下航行体流体动力特性分析

利用CFX对水下航行体流体动力特性进行了数值模拟研究.求解了混合介质的剪切压力输运湍流模型控制方程、雷诺-平均奈维尔-斯托克斯方程以及各相间的质量输运方程,采用三维数值模拟方法对比分析了不同流动控制方案的水下航行体空泡形态特征、表面压力分布和阻力系数变化特性,讨论了不同参数对减阻效果的影响.结果表明,对于有横流作用的水下垂直发射航行体,多相流动控制不仅可以降低通气空泡的不对称性和航行体阻力,同时可以均匀迎流面和背流面压力,从而实现航行体多相流空泡形态及表面压力特性等的控制.此外,通气口的位置对减阻效果具有显著影响.      

多块网格计算跨音压气机带叶尖间隙流场

开发了三维数值模拟程序研究轴流跨音转子叶尖间隙流动,应用高雷诺数k-ε湍流模型加壁面函数的方法,计算了轴流跨音转子NASA Rotor37在设计转速下的流场.叶尖间隙采用分区的H型网格和主流区连续对接耦合计算,没有用间隙模型,也没有考虑Vena收缩效应而减小间隙量.在用有限体积法对Navier-Stokes方程和湍流方程进行空间离散的过程中采用了交错网格的方法将N-S方程与湍流方程紧密地耦合在一起,从而提高了计算精度.计算结果和实验数据进行了详细的比较和分析.结果表明,中部叶展具有与实验结果非常一致的流场特征,根、尖区流场则因涡粘假设和激波问题的存在而使流动细节与实验结果略有偏差.     

肼电弧喷射推力器工作方式可替换性数值研究

在正确理解电弧喷射推力器工作过程物理机理的基础上建立了其简化的物理模型,并基于此物理模型利用耦合电磁源项的N-S方程组来表达其数学模型,对相同结构尺寸的低功率肼电弧喷射推力器加电弧和不加电弧时的工作性能进行了比较计算.采用二阶精度无波动、无自由参数的耗散差分格式(NND格式)和显式时间推进法求出数值解.计算结果给出了有电弧和无电弧时推力器的流场结构,预计了不同工作方式及不同工作条件下推力器的推力和比冲,证实了推力器工作方式的可替换性.      

缝道参数对多段翼型气动性能的影响

数值模拟研究缝道参数对着陆构型多段翼型气动性能的影响,采用有限体积法求解雷诺平均N-S(Navier-Stokes)方程,湍流模型采用S-A(Spalart-Allmaras)模型,给出了多段翼型升力系数与表面压力系数随各缝道参数的变化规律.计算结果表明:采用的计算模型可用于二维多段翼型的计算;缝道参数对多段翼型气动性能的影响显著,在各个缝道参数的变化范围内升力系数均存在最优值,适当的缝道参数可以抑制后翼上表面逆压梯度,消除后翼分离,并保持前翼尾迹不与后翼边界层相掺混,提高多段翼型的气动性能.      

用于叶轮机械流场计算的交错网格有限体积法

在对各种交错网格离散格式的奇偶不耦合数值误差进行了详细的傅里叶分析和比较后,采用一种紧致的交错网格有限体积法将纳维-斯托克斯方程和k-ε湍流模型方程紧密的耦合在一起求解,发展了一套三维定常流程序.通过对NASA Rotor 67设计和非设计工况内部流场的数值模拟结果与实验的比较、分析,初步验证了程序的准确性和可靠性.      

某高空螺旋桨气动特性数值模拟与风洞试验

以某高空螺旋桨为研究对象,利用k-ε湍流模型,采用三维Navier-Stokes方程,并通过数值计算模拟该螺旋桨的复杂绕流,分析其气动特性.按照等雷诺数等前进比准则,进行了该螺旋桨地面缩比风洞试验.经过比较发现,计算结果与试验结果吻合良好.验证了数值模拟计算高空螺旋桨气动特性的正确性,为高空螺旋桨的设计和试验提供参考依据.     

PANS方法在双圆柱绕流数值模拟中的性能分析

PANS方法以传统的雷诺平均Navier-Stokes(RANS)方程为模板,引入控制参数对原方程进行修改,是一类能够有效模拟分离湍流流动的混合方法.采用基于Menter 剪切应力输运(SST)湍流模型的PANS方法,选取双圆柱绕流作为研究对象, 分析双圆柱流动结构,验证方法的准确性与可靠性,考察PANS方法中模型参数选取对于计算结果的影响,评估不同类型PANS方法的性能.研究表明:模化湍动能比例这一模型参数对计算结果有显著影响,且计算域内统一的模型参数取值难以处理复杂流动问题.在此基础上使方法中模化湍动能比例取值随流场物理信息和网格尺度变化,所得结果与SST DES计算结果及实验结果吻合较好,表明该方法适用于复杂湍流数值模拟.      

跨音速风扇全环叶片颤振特性的流固耦合分析

发展了求解叶片颤振问题的流固耦合计算方法和全环叶片振动的气动弹性模型,在每一时间步同步求解流体运动方程和叶片振动方程并交换边界信息;流体域求解了非定常雷诺平均N-S方程,得到每一步由于叶片变形而引起的流场变化;叶片变形则由积分叶片表面受到的气动力并求解结构动力学方程得到.颤振分析是在全环叶片模型上进行的,并解除了预先设定叶片间相位角的限制.此方法的显著特征是在一次气动弹性计算过程中,可同时分析叶片多个固有模态、多个节径下的气动弹性稳定性,大大提高了使用时域法进行叶片排气弹分析的计算效率.考察了NASA rotor 67风扇全环模型在堵塞点、最高效率点和近喘点3个气动工况下,节径变化对叶片气动弹性稳定性的影响,给出了不同模态下气弹最不稳定状态对应的叶片振动节径形式.结果表明,振动形式对于叶片气动弹性稳定性的影响很大.     

前飞旋翼三维湍流场的数值模拟

采用计算流体力学CFD(Computational Fluid Dynamics)方法,分析正常前飞和近地前飞时直升机的流场与性能.在直角坐标系中,采用有限体积法和压力耦合方程的半隐式法SIMPLER(Semi-Implicit Method for Pressure-Linked Equations Revised)算法求解三维定常不可压的湍流N-S(Navier-Stokes)方程.在分析中,旋转的螺桨被描绘成沿螺桨桨叶展向分布、与本地流动参数相关以及时间平均的动量源项.计算方法还包括K-ε湍流模型和壁面函数法等措施.流场分析和性能预测同实验测量数据的良好的一致性表明,这种CFD方法可以有效的分析直升机的具体设计问题.      

固体火箭发动机燃烧室三维流动数值计算

采用k-ε双方程湍流模型,以SIMPLE计算程式求解翼柱型装药固体火箭发动机燃烧室内三维非定常不可压流N-S方程.在建立翼柱型装药简化模型的基础上采用边界标志法来表达燃面推移,采用从二维到三维的初场给定方法结合多重网格法求出了多个时间步的非稳态流场结构.计算结果表明燃烧室内旋涡运动呈现一定的空间与时间分布,周向加质导致燃气通道横断面上轴向速度分布的极不均衡.      

基于弱耦合的翼伞气动变形数值模拟

对不带气室翼伞和带气室翼伞的气动特性和结构变形进行三维数值模拟,分析气室对翼伞气动特性和结构变形的影响.流场方面,采用有限体积法求解Navier-Stokes控制方程,选用剪切应力输运(SST,Shear-Stress Transport)k-ω湍流模型;结构方面,假定翼伞有初始形状,使用插值方法映射传递流固交界面的压力数据,通过ANSYS计算翼伞伞衣的气动变形.结果表明:气室对阻力影响不大,对升力的影响主要表现在大迎角情况下;翼伞柔性和尺寸大小对其有利迎角的范围影响不大;翼伞的气动变形和最大主应力主要集中在气室中前部分,带气室翼伞由于有肋片加强的缘故,伞衣气动变形较小,一定程度上保证了翼伞的气动特性.