不可压缩翼型绕流数值研究

使用标准k-ω模型及与色散模型相耦合的k-ω模型分别计算了NACA0012翼型和NACA 4412翼型的低速绕流问题.NACA 0012翼型计算了其来流雷诺数为2.88×106,攻角从0°到15°范围内的流动结构、翼型表面压力分布和升力、阻力特性;NACA 4412翼型计算了临界雷诺数为1.52×106,攻角为13.87°时的流动分离和翼型表面压力系数 ,并与实验数据进行对比.结果表明:在同等条件下,使用与色散模型相耦合的k-ω模型计算得到的NACA 0012翼型的升力和阻力系数比标准k-ω模型提高精度约5%,NACA 4412翼型的表面压力系数精度提高了约3%,进一步验证了其可信性,可将其进一步应用到低速飞行器的气动计算中.      

湍流模型对发动机排气系统流场和红外辐射特征计算的影响

计算分析了采用k-ω SST(shear-stress transport)模型、RNG(renormalization group)k-ε模型配合标准壁面函数(RNG k-ε+标准壁面函数)、RNG k-ε模型配合增强型壁面函数法(RNG k-ε+增强壁面函数)3种湍流模拟方法,对典型涡扇发动机轴对称收敛喷管排气系统流场、温度场、组分质量浓度场及红外辐射特征计算结果的影响.研究表明:对于这3种方法,排气系统腔体内壁面压力分布基本一致;但在温度场、组分质量浓度分布以及红外辐射特征方面的差别较大.RNG k-ε+增强壁面函数计算得到的内涵壁面与燃气流的温度、CO₂与H₂O的质量浓度和红外辐射强度值最大,RNG k-ε+标准壁面函数的计算结果其次,而k-ω SST 模型的计算结果最小.通过与试验结果对比,发现k-ω SST模型对排气系统红外辐射特征的模拟更加合理一些.      

使用GAO-YONG湍流方程组对翼型绕流的计算

将GAO-YONG湍流方程组应用到翼型复杂绕流场的计算中,研究了机械能方程中Du/Dt项和各向异性湍流黏性系数分量中的实度系数Cs的影响,给出了控制Du/Dt项作用区域和选择实度系数Cs值的判别函数,对NACA4412翼型最大升力迎角下的低速绕流算例、AS240翼型8°迎角下的转捩算例、AS240翼型19°迎角下的分离流算例进行了模拟,与实验结果比较,得到了比较好的结果.结果表明GAO-YONG湍流方程组能够准确预测翼型前缘层流到湍流的转捩问题,翼型尾迹流动松弛问题和翼型大迎角分离流.     

湍流模型对压气机数值模拟精度的影响

开发了一个求解叶轮机全三维黏性流场的程序,采用了目前叶轮机数值模拟中常用的且模拟性能较好的三种湍流模型:k-ε模型,k-ω模型,Spalart-Allmaras(S-A)模型,计算了跨声压气机NASA(National Aeronautics and Space Administration) Rotor 37的流场,并与实验进行了性能参数、三维流场和角区流动的对比分析.考察在相同的数值计算平台上进行比较这几个湍流模型对压气机的数值模拟性能的影响.最终结果表明:对于具有强剪切、存在分离流的复杂的叶轮机三维流场来说,k-ω模型数值模拟精度更高,相比其他两个模型具有一定的优势.      

用温度修正的k-ε模型研究排气系统红外辐射特性

为研究温度修正的k-ε湍流模型对排气系统红外辐射特性计算精度的改善程度,计算了某轴对称收敛喷管在3~5μm波段的空腔-喷流组合红外辐射特性,并与试验结果对比.喷管流场及温度场采用有限体积法(FVM)求解Navier-Stokes(N-S)方程得到,湍流模型分别采用标准k-ε模型和温度修正的k-ε模型.红外辐射特性的计算采用有限体积法求解吸收-发射性介质条件下的三维辐射传输方程,并考虑大气衰减作用.计算结果表明:温度修正的k-ε模型可有效预测排气系统热喷流的核心区长度,从而获得与试验结果一致的红外辐射强度.与标准k-ε模型相比,采用温度修正的k-ε模型计算得到的红外辐射强度的最大相对误差由64%降至了4.3%,说明温度修正的k-ε模型适用于精确模拟排气系统的红外辐射特性.      

NACA4412翼型低速绕流数值计算中湍流模型对比

使用Spalart-Allmaras(S-A)、SSTk-ω、Gao-Yong 3种湍流模型对NACA4412翼型低速绕流进行了数值计算,研究了尾迹流动松弛效应。对流项采用Roe格式离散,扩散项采用二阶中心格式离散,离散后的控制方程用多步Runge-Kutta显示时间推进法求解。计算中对翼型尾缘流松弛效应进行了分析,比较了翼型表面压力系数、速度剖面、雷诺应力等的分布,3种湍流模型总体上能够较好地模拟NACA4412翼型低速绕流。SSTk-ω模型对流动细节及升力系数计算最好,Gao-Yong模型对翼型平均速度剖面及雷诺剪切应力分布计算最准确。     

翼型大攻角绕流的数值模拟

以求解二维N-S方程数值模拟NACA0012翼型大攻角状态的可压绕流特性;N-S方程是在贴体坐标系中给出的,以代数方法生成C型网格系统。采用LU-ADI格式计算,为提高格式的稳定性在隐式和显式部分分别添加了2阶和4阶人工粘性项。应用BaldwinLomax湍流二层代数模型模拟了大攻角时前缘分离涡的形成,旋涡对流及其非定常现象。在某些Mach数和攻角下NACA0012翼型的湍流解具有周期性。通过比较,本文数值计算结果同实验及国外相应的数值计算结果基本吻合。     

轴承腔油气两相流数值模拟的试验验证

为探索一种适于模拟轴承腔油气两相流流动特性的数值计算模型,建立了一套简化结构的轴承腔试验装置,通过试验研究获得了不同工况对轴承腔通风口出气比例和回油口滑油温升的影响规律.应用Fluent软件,采用Mixture多相流模型和Eulerian多相流模型分别搭配标准k-ε湍流模型和RNG(renormalization group) k-ε湍流模型对选定试验工况进行数值模拟.对比结果表明:当转速为2200r/min时,使用Mixture多相流模型搭配RNG k-ε湍流模型计算的滑油温升与试验值的相对误差仅为2.72%,体现了该计算模型在高转速时的精确性,并在此基础上获得了油气两相流速度场、压力场和温度场的分布特征.      

绕翼型分离流结构的数值研究

利用大涡模拟技术,对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟,详细地给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数。在数值模拟中,采用了弱压缩流的控制方程,用贴体坐标技术进行了网格生成。     

新型等离子体激励器对流动分离点控制

在西北工业大学低湍流度风洞中采用新型等离子激励器对NACA0015翼型进行表面流动分离点的控制实验。实验风速为20m／s和35m／s,迎角为0°～16°。并参照压力分布的实验结果对流动控制的效果进行了对比分析。结果表明：翼型表面的气流分离点只要落在等离子体激励所形成的激励区内,分离点都会被推迟到靠近等离子体激励器的最末端电极处。证明等离子激励器能够对翼型表面的分离点进行有效控制。     

表面凹凸充气机翼的气动特性研究

充气机翼在变体飞机和飞艇中具有潜在的应用前景。充气机翼的结构特征与传统硬质机翼显著不同,其外形与传统机翼相比最大的差异在于表面的片条状鼓包,这种外形带来的气动特性、气弹行为等越来越受到人们的关注。以NACA0015翼型为原形,设计制作充气机翼模型,并利用CFX对传统光滑的NACA0015翼型和凹凸表面的0015F2翼型进行定常和非定常气动行为分析。结果表明:充气机翼的凹凸表面外形增加翼型的失速攻角,但其升力线斜率及升阻比都较光滑翼型要小;0015F2翼型的速度梯度过度区大于NACA0015翼型;0015F2翼型在每一个凹槽区生成驻涡,驻涡的存在使得充气机翼的附面层呈现紊流附面层的特性,驻涡的外移改变了机翼后缘的尾涡形成,推迟了分离,使得失速攻角增大。     

翼型失速的非线性动力学特性及其控制

以NACA0012翼型为研究对象,应用CBS( Characteristic-Based Split)有限元方法对翼型绕流流场进行数值模拟,着重对翼型失速现象的动力学机理和特征进行研究.首先给出了静止翼型在不同迎角下流场的瞬时流线图以及静态翼型的迎角-升力系数关系曲线,并且分析出此曲线具有滞后、突跳等典型的非线性动力学...     

求解平面翼型亚,跨声速绕流的一个新方法

本文提出了求解平面翼型亚、跨声速绕流的一个新方法。引入流函数和Ｖｏｎ Ｍｉｓｅｓ变换后，亚、跨声速绕机翼无旋流动的基本方程组被化为以流线纵坐标ｙ为未知量的单个二阶偏微分方程－流线控制方程。并通过变换将物理平面上的无限域变为计算平面上有限的矩形域，而后在计算平面采用有限差分线松弛迭代法求解。作为算例，计算了对称翼型ＮＡＣＡ００１２－３４和非对称翼型ＮＡＣＡ４４１２的亚、跨声速有攻角绕流，所得数值结果     

基于N-S方程的地面效应数值模拟研究

研究了二维翼型在两种不同地面条件-移动地面和固定地面条件下的地面效应及流动特性问题。其中流场分析采用了不可压缩RANS方程,湍流模型采用Spalart-Allmaras一方程湍流模型。通过对NACA4412翼型在两种不同的地面条件及不同迎角和不同地面高度情况下的绕流计算,结果表明,在两种不同的地面条件下,升力及力矩的差别较小,而固定地面条件的阻力略小于移动地面条件下的阻力。     

不同迎角的翼型气弹特性风洞实验研究

基于可在不同迎角下作沉浮、俯仰两自由度运动的翼段振动装置,在低速风洞中分别针对普通薄翼型NA-CA0012和风力机翼型NREL S809进行气动弹性测试,得到不同实验状态的气动弹性振动时域响应。分别观察到经典颤振和失速颤振现象,并证明了迎角改变对两种翼型颤振特性的影响。     

基于动网格法的翼型启动过程数值模拟

采用动网格法对NACA4421翼型以15°攻角启动过程进行了数值模拟。计算给出了启动过程中尾缘启动涡的生成、脱落与绕翼型环流充分发展的瞬态流场及气动力特性变化曲线,并对升力数据进行了拟合。计算结果表明,启动瞬间,上翼面最大负压和上下翼面最大压强差均出现在翼型后半段,随后逐渐向前缘移动,最终稳定在前缘点附近。下翼面最大正压点和上下翼面压差随弦向位置的最大变化率则始终维持在前缘点附近。加速过程中,整个翼型受到的升力近似于瞬时速度的二次幂指数的规律变化。加速段结束后,翼型转入匀速运动的瞬间出现升力小幅下降的现象,之后逐渐回升至稳定升力。     

任意翼型有小分离气泡跨音速绕流的迭代算法

 本文给出一种带小分离气泡的任意翼型粘性跨音速绕流的计算方法,采用有粘-无粘干扰迭代的概念。无粘流的全速势方程用AF差分格式在保角变换法生成的计算网格中求解,粘流附面层方程用C-S盒式法求解,用逆算法消除分离点处的奇性。本文对Ma_∞=0.8,Re_∞=2×10~6,迎角α=3.5°和4°的NACA64A010翼型粘性绕流进行了计算,结果与实验相比较,吻合良好。     

无轴涵道旋翼气动特性数值研究

采用基于非结构网格的滑移网格技术对无轴涵道旋翼与传统涵道旋翼的气动特性进行了非定常Euler方程数值模拟。分别考察了无轴涵道旋翼与传统涵道旋翼气动特性差异,中心孔径、涵道扩散角、涵道翼型对无轴涵道旋翼气动特性的影响。涵道扩散角变化范围为-6°~10°,涵道翼型选择NACA66、NACA0018及NACA4415三种翼型进行研究。研究发现:涵道翼型类型对无轴涵道旋翼的拉力分配影响较大,对称翼型能减弱旋翼上方低压涡从而涵道能产生更大的拉力；无轴涵道旋翼比传统涵道旋翼具有更优的拉力性能,转速为18 000 r/min时其总拉力是后者的1.185倍；减小无轴涵道旋翼的中心孔径能提高总拉力值,但整体耗能也将随之提高；涵道拉力占比越高的翼型,其最佳涵道扩散角越大,该状态下涵道拉力占比最高。      

用粘流－无粘流相互作用方法计算翼型的分离和失速

本文用势流－边界层相互作用方法计算低速翼型的分离和失速。势流用对称面元法。边界层用改进的滞后掺混法，考虑了高阶项影响，适用于计算分离。文中对粘流－无粘流耦合方法作了改进。改进的半反－局部联立耦合方法，考虑了相邻点之间的作用，收敛性较好。计算了ＮＡＣＡ４４１２翼型在不同迎角下的压力分布和气动力。计算结果与实验符合良好。算例表明，翼型高升力状态计算必须包括尾流的作用，也应当对势流计算压力与实际压力之间的差别进行修正。     

SLAU格式在大涡模拟方法中的应用

采用格心格式的有限体积法求解三维大涡模拟控制方程,成功地将最新改进的SLAU格式推广到大涡数值模拟方法中。SLAU格式的空间重构方式采用五阶WENO格式,时间推进采用LU-SGS隐式时间推进。湍流模式为最新提出的基于随机分析的动态亚格子模型。计算了圆柱绕流、NACA4412翼型绕流、NACA0012翼型绕流三个算例,分别从可压缩较大迎角翼型绕流、可压缩跨音速翼型绕流以及不可压缩低马赫数圆柱绕流三个方面验证了将SLAU格式应用到湍流大涡数值模拟计算中的可行性。结果表明,结合SLAU格式的大涡模拟方法可以达到较好的数值模拟效果,同时间接检验了基于随机分析的动态亚格子模型的数值模拟效果。