带自由表面水流的分离涡模拟

目前水动力学问题的数值求解仍然广泛采用RANS模型,特别对于大尺度的地表水流的数值模拟,该类模型计算效率较高,但其难以给出湍流场的小尺度涡结构.深入研究复杂地形条件下的水流运动,需要发展诸如LES等的高级数值模型.天然地表水流,地形、边界复杂,雷诺数通常较高,高级数值模拟受限于计算能力的限制,还难以工程应用.RANS和LES混合模型是目前具有工程应用的一种混合数值模型,其已经在CFD领域获得了长足发展.然而对于地表水流运动的数值模拟,相关工作还较少.本文将分离涡模型(DES),即一种RANS和LES的混合模型,应用于带自由表面的地表水流运动,建立了一套数值仿真模型.模型基于有限体积法,水平面内采用非结构计算网格,垂向为结构化网格,对流项离散格式采用二阶TVD格式,并行基于OpenMP语言库.算例表明DES模型有助于揭示复杂地形条件下带自由表面水流的大涡拟序结构.     

基于非结构网格的不可压N-S方程多矩有限体积法

提出了一种基于三角形及四面体非结构网格的有限体积法(FVM),用以鲁棒且精确地求解不可压粘性流动问题.与传统的FVM方法仅将体积分平均值(VIA)作为计算变量的做法不同,本文提出的方法将VIA及点值(PV)同时作为计算变量并在每个迭代步进行计算更新.VIA以通量形式进行计算以确保数值守恒,PV可以通过控制方程的不同形式进行求解更新,无需守恒,因此可以采用非常高效的方法进行求解.将PV作为增加的变量使得紧致网格模板得以实现更高阶精度的重构,而且由此获得的数值模型对于非结构网格变得更鲁棒.本文针对二维/三维的三角形/四面体非结构网格提出了数值格式,给出了几个基准测试算例,验证了本文提出的数值方法在采用非结构网格求解不可压粘性流动问题时的精确性和鲁棒性.     

跨声速串列转子失速机制的数值研究

为了解决跨声速串列转子的低裕度问题,就必须了解跨声速串列转子的流场结构与失速机制.设计了叶尖切线速度为450m/s,负荷系数为0.56的高负荷跨声速串列转子.基于数值模拟的结果,分析了该串列转子在0.5mm叶尖间隙下的叶尖流场结构与失速机制,并在此基础上分别探讨了叶尖间隙和前、后排叶片周向位置对串列转子特性的影响和失速机制的变化.结果表明:前排叶片的叶尖区域是 影响串列转子稳定性的关键;随着叶尖间隙的增加,串列转子的失速机制也发生变化,从前排叶片叶尖区域的尾迹与径向潜流堵塞后排叶片通道转变为前排叶片叶尖泄漏流堵塞;在较大周向相对位置(后排叶片压力面周向远离前排吸力面)的情况下,串列转子获得最好的效果,随着周向相对位置(PP)的增加,失速部位从后排叶片转移至前排叶片.      

蜂窝式轴心通风器油气分离性能计算

为对航空发动机蜂窝式轴心通风器油气分离效率进行研究,建立考虑油气双向耦合的流场计算方法及油滴/壁面相互作用模型,在验证通风阻力及油气分离效率可靠性的基础上,对不同转速、通风流量和环境温度下蜂窝式轴心通风器的油气分离效果进行计算和分析.结果表明:转速的增加会使油气分离效率得到提升,而通风流量和环境温度的增加则导致油气分离效率的降低.蜂窝孔结构的加入对通风阻力影响不大,却对通风器的滑油分离过程起主要作用,计算表明其对滑油分离贡献率在80%以上.      

New constitutive relation for eddy viscosity models

A new turbulent constitutive relation was directly derived from Boussinesq's hypothesis and mixing length theory, and then implemented in the standard k-ε model. The performance of this constitutive relation was validated in zero pressure gradient flat-plate boundary layer flow, fully-developed turbulent channel flow and separated flow in a plane asymmetric diffuser. The investigation demonstrated that, this new constitutive relation gave very accurate results in the former two basic cases and provided significant improvement in prediction of separated and reattachment points in the plane asymmetric diffuser. Separation and reattachment points at x/H=7.5 and 29 were calculated accurately in comparison to experimental results, and the static pressure coefficient of 0.82 was very close to large eddy simulation calculation. These results are very encouraging but further verification and extensive application of the new constitutive relation to other two-equation eddy viscosity model are needed.      

温度对飞机燃油射流泵工作特性影响

建立了射流泵数学模型,采用数值计算方法对其工作特性进行了模拟,搭建了实验台,采用RP-3燃油,通过实验数据验证了模型的正确性.在此基础上,分别选择了RP-3和JET-A两种燃油,计算了从-40℃~40℃温度下,不同操作压力对射流泵流量比的影响.结果表明:燃油温度高于0℃时,温度对射流泵工作特性影响并不明显,随着温度的降低,射流泵性能与标准温度下（20℃）偏差越来越大,且压比越高此偏差越明显.当温度高于0℃时,采用RP-3和JET-A燃油性能相差不大,RP-3燃油性能略优于JET-A燃油,但是在低温下,RP-3燃油性能远高于JET-A燃油.因此当燃油温度较低且压比较高时,必须充分考虑温度对射流泵特性的影响.      

平箔片楔形高度对气体止推箔片轴承特性影响

针对波箔型气体止推箔片轴承,建立了箔片结构二维薄板模型,并通过有限差分法和有限元法耦合求解可压缩气体Reynolds方程和气膜厚度方程,获得了给定轴承载荷条件下轴承气膜压力分布、气膜厚度分布、平箔片变形量和功率损耗等轴承特性.通过对比楔形高度分别为25,70,100,200,300μm时的轴承特性仿真结果研究了楔形高度对轴承性能的影响.结果表明:降低楔形高度使轴承气膜压力分布更均匀,并降低了平箔片的局部集中载荷.但楔形高度存在一个最佳值,使得达到相同轴承载荷所需的最小气膜厚度最大,并且轴承具有最小功率损耗,提高了轴承的工作效率.该结果为气体止推箔片轴承的结构设计提供的理论参考.      

水下超声速射流对上浮水雷受力特性影响研究

固体火箭发动机具有功率密度大、推力大等优势,常被用于上浮水雷的推进器。然而发动机在水下工作时燃气射流使流场压力发生剧烈脉动,进而影响发动机的推力性能及上浮水雷的受力特性。基于VOF多相流模型和理想气体模型,建立了上浮水雷在不同工作状态(欠膨胀、完全膨胀、过膨胀)下的数值模型,研究了水下燃气喷射流对上浮水雷的受力特性影响。结果显示,欠膨胀工况时,发动机推力平缓,大小为12.2KN,上浮水雷受力未出现负值;完全膨胀时,射流发生颈缩、胀鼓现象位置距离雷体较远,导致发动机推力及水雷壳体受力振荡不剧烈;过膨胀工况时,射流发生颈缩、胀鼓现象距离喷管较近,发动机推力发生剧烈脉动,产生21.37%的振荡幅度,胀鼓现象发生时,流场压力显著降低使得上浮水雷后体受力减小,壳体阻力增大,上浮水雷最大产生27KN的负推力。     

基于自适应网格的超声速液体射流破碎过程研究

为研究超声速气流中液体横向射流的破碎过程,采用脉冲背景光方法和VOF方法开展了实验和数值研究。为提高液体横向射流中气液界面和气流场特征捕捉的精确性,采用自适应网格技术对于气液界面、激波出现位置进行网格细化,计算得到了较为精细的气液界面、激波特征及涡系结构。研究结果表明:在低成本仿真模拟条件下,利用自适应网格计算得到的射流轨迹和轮廓与实验吻合较好,射流轨迹的最大误差为10%;射流初始段在超声速气流条件下,仍然存在一段高度约为1.9倍喷孔直径且圆柱形态保持较好的连续光滑液柱。随着喷注压降的升高,液柱的长度逐渐增大;主流气体流经液柱发生三维绕流,在射流附近和近壁面区域形成不断演化的反转涡对,反转涡对的形成加速了液体射流一次破碎过程。     

Al/Ni含能多层膜在Al2O3陶瓷/不锈钢焊接中的应用

采用磁控溅射设备,生长AuSn合金做焊料层、Al/Ni含能多层膜做热量提供层,实现了不锈钢和Al_2O_3间的异质材料自蔓延高温扩散焊。利用SEM、XRD和DSC等测试手段表征AuSn合金和Al/Ni含能多层膜的微观形貌、相成分和放热量;用万能试验机测试焊接接头的力学性能。结果表明,AuSn合金的质量比基本达到80∶20,而多层膜的层状结构清晰,反应热达到1 239 J/g。焊接实验结果表明,仅使用AuSn焊料时,剪切强度仅为46 MPa,在增加Al/Ni含能多层膜后,其剪切强度可达90 MPa,强度提高了约一倍。焊接接头的界面显微形貌和相结构研究表明,剪切强度的增强主要是Al/Ni多层膜提供了额外能量使得界面处的反应剧烈,陶瓷金属化层与中间层的反应加剧,形成了新的反应生成物。