弹丸旋转空气动力效应非定常数值模拟

准确计算马格努斯力和力矩对旋转弹箭设计、弹道计算和稳定性研究都至关重要。采用非定常雷诺平均Navier-Stokes(RANS)方程对高速旋转的SOCBT弹丸进行了数值模拟,纵向气动特性在整个攻角范围内都较好地与试验值保持一致,而侧向气动特性马格努斯力和力矩在小攻角范围与试验数据吻合较好,在大攻角范围却存在一定的差异。采用延迟分离涡模拟(DDES)方法的计算结果有较为明显的改善,对比研究表明分离点位置对马格努斯效应有着显著影响。表明DDES方法对于提高旋转弹箭马格努斯效应的数值模拟精度有较大的潜力。     

壁面函数在激波诱导分离流动中的应用

以数值模拟激波-附面层干扰引起的流动分离问题为研究背景,发展了基于有限体积方法的雷诺平均Navier—Stokes（RANS）方程的流场数值模拟方法。利用壁面函数模型得到壁面剪切应力,通过修正壁面粘性通量,构造了一种新的湍流边界处理方法,并将其耦合到RANS方程和SSTk-ω湍流模型的数值求解中;同时,针对激波诱导引起的附面层流动分离问题,提出一种附面层网格加密技术,能够自适应加密分离区内附面层网格,使得在流动分离区域也能够使用壁面函数模型。数值算例表明,壁面函数模型能够降低数值模拟结果对网格的依赖性;同时也验证了壁面函数耦合附面层网格自适应方法,在处理激波诱导引起的附面层流动分离问题时的有效性和准确性。     

正攻角状态下尾迹对低压涡轮附面层的渗透及转捩的影响

为研究正攻角状态下尾迹对低压涡轮附面层的渗透以及转捩过程的影响,分别在设计攻角和+10°攻角下,对高负荷低压涡轮叶型的吸力面附面层流动进行了数值模拟与实验。数值模拟使用CFX软件,采用大涡模拟模型。结果表明+10°攻角下,尾迹对附面层转捩过程的促进作用比设计攻角下更为显著,这是由于正攻角下的尾迹中心射流与吸力面切向的夹角更大,使得尾迹扰动渗透入附面层的深度更深,尾迹放大Klebanoff条纹强度更强,尾迹诱导转捩起始位置更靠上游。     

非设计状态下尾迹输运对高负荷低压涡轮附面层的影响

为了研究非设计状态下的上游尾迹在叶栅通道内的形态演化与发展,分析尾迹与叶片吸力面附面层的相互作用,基于高负荷低压涡轮(LPT)叶型Packb对非设计状态下尾迹输运进行了研究。研究主要通过数值模拟的方法进行,使用CFX软件,利用LES模型耦合Smagorinsky亚格子模型。讨论了非定常来流0°和+10°攻角工况下,尾迹与吸力面附面层相互作用的差异。分析发现,+10°攻角工况时,尾迹对附面层转捩的促进作用较0°攻角工况时更为显著;+10°攻角工况时,尾迹与附面层的相互作用时间更长,尾迹诱导转捩的起始位置更靠上游。     

混压式进气道与弹体一体化流场数值模拟

基于Favre平均的N－S方程和B／L代数湍流模型,采用Jameson格式和矩阵人工粘性,对X布局的混压式超声速进气道与弹体的一体化流场进行了数值模拟,研究了解决了数值计算中的附面层抽吸及混压式进气 启动问题,得出了进气道内外流场的马赫数分布和速度分布,讨论了绕弹体非均匀来流条件下,进气道的位置,附面层抽吸及攻角对进气道性能的影响和进气道外型形状对弹性气动力的影响,结果表明,对混压式进气道必须进行附面层抽吸,进气道的位置对进气道的性能有很大的影响,进气道的外型型面对弹体的气动性能有一定的影响。     

高压燃气涡轮径向内冷叶片气热耦合的数值分析

采用气热耦合方法对采用径向内冷方式的MarkⅡ型跨声速高压燃气涡轮金属导叶进行数值模拟,通过分析叶片通道内的传热和流动过程发现叶片表面附面层内流动非常复杂,包含层流流动、转捩和湍流流动状态,所以只有使用转捩模型计算的叶片附面层内流动与实际情况相符,叶片壁面温度和换热系数分布与实验结果吻合的较好,使用其他湍流模型由于不能准确描述附面层内流动而使得计算结果误差相对较大,但是所有的湍流模型都能很好的模拟附面层以外流动.      

高阶精度格式WCNS在三角翼大攻角模拟中的应用研究

采用5阶精度的加权紧致非线性格式(WCNS-E-5)数值模拟了65°三角翼的大攻角绕流流场,主要目的是考核高阶精度格式WCNS在大攻角旋涡流动方面以及跨声速流场的激波附面层干扰、涡破裂位置的模拟能力,重点研究不同网格规模和湍流模型对尖前缘三角翼涡系之间的相互作用的影响。通过求解任意坐标系下的雷诺平均N-S方程,采用5阶精度的加权紧致非线性格式(WCNS-E-5)和多块对接结构网格技术,两种湍流模型分别是一方程SA和两方程SST湍流模型,在与相应试验结果对比的基础上,详细研究了WCNS-E-5格式在跨声速大攻角旋涡流动中的表现,以及不同网格规模、两种湍流模型对主涡二次涡相互作用、涡破裂位置和表面压力分布的影响。本文的研究结果表明,高阶精度格式WCNS-E-5能成功应用于三角翼的跨声速大攻角流动,网格规模的增加进一步提高流场分辨率,SST湍流模型相对SA湍流模型在三角翼大攻角流动中具有更好的适用性。     

用N－S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二级时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型￣［１］。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

H型叶栅通道流计算网格的改进

提出一种壁面正交、尾迹区网格节点分布逐渐均匀化叶栅通道网格生成方法。为了研究壁面处网格正交性对数值计算结果的影响，推导出壁面附近流动参数一阶导数离散误差与网格线交角的数学关系。采用N-S方程分别对平板紊流附面层流和叶栅通道流进行计算分析。结果表明：壁面处网格正交性对附面层内流动计算预测精度影响较大，而对势流区流动影响很小；尾迹区切向网格节点逐渐均匀化分布能更准确计算尾迹区流动。     

用N—S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二维时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

船尾形状对旋成体马格努斯效应的影响

尖头旋成体和船尾形状是子弹、炮弹及火箭弹等抛射体上常用的布局形式。研究表明船尾布局具有减小底部阻力、增大射程的作用,但此时旋成体的马格努斯效应增大,对运动稳定性产生不利影响。为了解释这种流动现象,对三维旋转弹流场进行了数值模拟,对从亚声速到超声速下的旋成体马格努斯力和力矩进行了分析,重点对标准形状和船尾形状两种底部进行了比较。结果表明,相对于标准形状,在所有来流下船尾形状都起到了增大马格努斯效应的作用,并且马格努斯力和力矩与船尾角成正比。为了揭示其流动机理,选择代表性计算状态对两种布局马格努斯力矩系数分布、边界层厚度分布和边界层位移厚度分布进行了对比分析,结果表明,在亚跨声速下船尾马格努斯效应由绕拐角的加速流动引起,使当地压力系数幅值增大;在超声速下船尾马格努斯效应由船尾段的气流膨胀引起,使旋成体左右两侧的边界层位移厚度畸变增大。上述两种效应都使马格努斯力矩增加,对于亚声速流动来说,该效应发生在柱段与船尾段连接位置;对于超声速流动来说,该效应发生在连接点以后的船尾段上。当来流速度在声速点附近时,上述两种效应都可能发挥作用,使船尾形状的旋成体马格努斯效应大幅增加。     

Simulation of the flow around a rotating circular cylinder taking into account laminar-turbulent transition

We consider a problem of simulating a viscous incompressible fluid flow past a rotating cylinder in the region of critical Reynolds numbers as well as the simulation of separated bubbles in the boundary layer and the reverse Magnus effect. The Fluent 12 software is used for numerical simulation. The transition k-kl-ω model of turbulence makes it possible to consider the laminarturbulent flow regime transition.     

转动园盘上可压缩流动换热的计算研究

本文将边界层方程无因次化,用“BOX”算法求解。研究计算了转动园盘上层流区、过渡区和紊流区的扭矩系数和局部换热系数以及变物性的影响,并研究了适合计算转盘问题的双层紊流模型。     

Effects of Mach numbers on Magnus induced surface pressure

Experimental and numerical methods were used to investigate the Magnus phenomena over a spinning projectile. The pressure force acting on the surface of a spinning projectile was measured for various cases by employing a relatively novel experimental technique. A set of miniature pressure sensors along with a data acquisition board, battery and storage memory were placed inside a spinning model and the surface pressure were obtained through a remotely controlled system. Circumferential pressures of the model for both rotational and static conditions were obtained at two different free stream Mach numbers of 0.4 and 0.8 and at different angles of attack. The results showed the ability of this new test method to measure the very small Magnus force via surface pressures over the projectile. The results provide a deep insight into the flow structure and illustrate changes in the cross-flow separation locations as a result of rotation. Similar results were obtained by the numerical simulations and were compared with the experimental data.     

旋转坐标系下分区计算的LU隐式方法

将隐式时间迭代方法应用于并行计算是CFD研究中的热点问题,研究发现LU-SGS格式在旋转坐标系下对分区边界比较敏感,对旋翼悬停状态进行分区并行计算时碰到了计算发散的问题。针对此问题,对基于LU分解的3种隐式时间迭代格式（LU-SGS、DP-LUR和HLU-SGS）进行了对比研究,设计了静止流场下的旋转网格算例对计算方法进行测试。结果表明,LU-SGS格式在网格边界采用简化处理方法,当边界的逆变速度增大时会引起数值误差放大,误差的不断积累导致计算发散。DP-LUR和HLU-SGS格式通过在边界单元采用雅克比迭代算法,能有效消除分区边界影响,使计算格式在大CFL数条件下保持稳定,其中HLU-SGS继承了LU-SGS迭代效率高的特点。在此基础上,采用并行化的LU隐式方法对Caradonna-Tung旋翼进行了并行计算,针对两套稀密程度不同的背景网格,在产生复杂分区边界的计算网格条件下,均获得了气动力和旋翼尾迹都充分收敛的流场。通过计算,分析了背景网格密度对预测气动力性能和捕捉尾迹流场的影响,计算结果与试验值吻合良好,验证了当前隐式计算方法适用于旋转坐标系并行计算,适合于推广至其他大规模并行分区的流场计算。     

用热线风速仪研究轴流压气机的非定常流场及端壁附面层

根据轴流压气机特点和热线风速仪特性,对测试技术作了一些新尝试。本文总结这些方法,并提供若干典型结果。大量试验表明这些方法有效。     

二维旋转湍流边界层内气固两相流动的数值模拟

本文用修正的ＤＯＲＯＤ－ＦＥＭ方法对不可压流体二维旋转湍流边界层进行了数值模拟，得到了与实验数据较符合的边界层形状因子和壁面摩擦系数的计算结果。采用单向耦合（Ｏｎｅ－ｗａｙ Ｃｏｕｐｌｉｎｇ）方法，计算了固粒在湍流边界层和无粘流中的运动轨迹，给出了壁面被固粒撞击后的冲蚀磨损分布，分析了近壁流场性态与壁面冲蚀磨损的关系。     

可压紊流边界层掺混方法力学模型的改进

 本文给出了二维可压紊流边界层的一种改进的掺混算法,方法简单,计算机时短,可用于计算激波后的紊流边界层的发展。六个算例的计算结果表明,本方法具有良好的数值精度。同P.D.史密斯的方法相比,结果有明显的改进。     

计算压气机特性的基元叶片法

本文介绍一种轴流压气机非设计性能的预测方法。该法特点是 ,在流线曲率法解径向平衡方程基础上 ,引入环壁附面层损失的计算 ,考虑叶型存在对气流熵增和总焓变化的影响 ,并计入二次流 /径向间隙流引起的损失 ,激波、尾迹损失等。结果与试验数据对比符合工程要求。