30P-30N多段翼型复杂流场数值模拟技术研究

采用"亚跨超CFD软件平台"(TRIP2.0)数值模拟了30P-30N多段翼型的复杂流场,主要目的是考核湍流模型、转捩位置对多段翼型压力分布和典型站位速度型的影响。本文通过求解任意坐标系下的雷诺平均的N-S方程,采用多块对接结构网格技术,在与相应试验结果对比的基础上,详细研究了SA一方程湍流模型、SST两方程湍流模型、不同的转捩位置对该翼型压力分布和典型站位速度型的影响。本文的研究结果表明,采用全湍流模拟方式可以较好地模拟该多段翼型的压力分布,但对速度型的模拟精度较差;模拟试验的转捩位置可以改善主翼附面层与前缘缝翼边界层尾迹区的模拟精度;采用微吸气技术推迟前缘缝翼的转捩位置,可以进一步提高缝翼尾迹区的数值模拟精度。     

一种剪切修正的涡黏性亚格子模型及其在槽道湍流中的应用

提出了一种利用壁面附近强剪切对涡黏系数进行修正的亚格子模型,称为剪切修正的涡黏模型.该模型能使涡黏系数在壁面附近自动满足壁面修正,在脉动为零的地方,涡黏系数自动为零,保证了壁面附近涡黏系数趋于零的特性.用该模型对槽道湍流壁面无吹吸、有吹吸,较低和较高雷诺数的情况进行了数值模拟,与直接数值模拟的结果进行了对比,结果吻合较好.      

适于黏性计算的自适应笛卡儿网格生成及其应用

复杂外形的黏性网格生成已经成为计算流体力学(CFD)在工程应用中的主要瓶颈,高效而鲁棒的网格生成技术研究显得尤为迫切。使用基于物面外形自适应加密的笛卡儿网格方法生成计算网格,投影方法拟合壁面边界,将投影得到的柱形单元沿法向分层获得适于黏性计算的网格。应用特征恢复技术恢复凹面特征并改善凸面处的网格质量。流场解算中应用最小二乘法重构获得2阶精度,Roe格式计算无黏通量,Holmes-Connell方法计算黏性通量。加入Spalart-Allmaras一方程湍流模型以获得湍流解。上下对称高斯 赛德尔迭代(LU-SGS)格式隐式时间推进,对网格重排序保证格式的稳定性和平衡性。算例显示网格生成方法是快速的和鲁棒的,流场计算精度满足工程需求,因此该方法非常适合应用于实际工程计算。     

湍流边界层厚度对三维空腔流动的影响

采用脱体涡模拟(DES)方法开展了不同湍流边界层厚度(TTBL)下的三维空腔非定常流动数值计算。空腔长、宽、深比例为5:1:1,来流马赫数为0.85,雷诺数为13.47×10⁶ m^-1,各工况湍流边界层厚度比值为1:2:4:8。研究结果表明,湍流边界层厚度对自由剪切层的发展、空腔底部静态压力分布、脉动压力及空腔流动类型均有重要影响,且随着边界层厚度的增大,下游剪切层覆盖的范围会增大,但是剪切层增长率降低;空腔前后静态压力压差减小、压力梯度下降;腔内局部测点的脉动压力声压级下降,各阶声压峰值频率向低频方向偏移;空腔流动类型往开式流动方向转换。     

气氧/煤油富燃燃气发生器燃烧试验研究

为了研究航天煤油的富燃燃烧特性,设计了一个气氧/煤油富燃燃气发生器,并进行了室压范围1~4MPa、混合比范围0.25~0.52的燃烧试验,获得了混合比和室压对特征速度、燃气温度及燃烧效率的影响规律。试验结果表明,采用火花塞点火方式可以实现该燃气发生器在高度富燃条件下的可靠点火;在本文研究范围内,试验测得的特征速度和燃气温度明显低于化学平衡理论预测值,燃烧效率介于63.9%~88.3%之间。混合比是特征速度、燃气温度及燃烧效率的最大影响因素,室压的影响远小于混合比。燃烧效率与实际燃气温度与理论燃气温度比值的平方根成正比。     

基于涡流发生器的风洞侧壁附面层控制

为降低风洞侧壁附面层对半模型数据的影响,在前期数值模拟的基础上,研制了一种适用于2.4m跨声速风洞半模型试验段侧壁的梯形涡流发生器,并进行了试验验证.结果表明:加装涡流发生器效果明显,亚声速范围内能够使附面层厚度降低20%~30%,对主气流均匀性影响可忽略;加装后半模型零升阻力系数降低,升力线斜率增大,压力中心向机身移动,体现了明显的附面层减薄效果,证明所研制的涡流发生器可应用于亚声速半模型试验中.      

整体涡对跨声速单转子性能和稳定性影响

为研究整体涡影响轴流压气机性能和稳定性的机理,设计了一种叶片式旋流发生器,并对旋流发生器和跨声速单转子进行了联合数值模拟研究.计算结果表明:反向整体涡导致堵塞点流量、进口攻角和总压比增加,效率和稳定工作范围降低;同向整体涡导致堵塞点流量、进口攻角、总压比和稳定工作范围减小,效率增加;反向整体涡导致50%叶高以上吸力面附面层分离以及叶尖泄漏流堵塞加剧是压气机失稳的主要原因.      

转捩模型在翼型气动性能预测 中的应用研究

边界层转捩位置的准确预测对于提高飞行器气动性能的预测精度具有重要意义.选取与k-ω SST湍流模型相耦合的γ-Reθt模型,以零压力梯度平板为研究对象,通过求解基于有限体积法的雷诺平均N-S方程验证该模型自动捕捉流动转捩的准确性;将该模型应用于传统有压力梯度的NACA 0012翼型的流场特性和气动性能的研究中,并与原始k-ω SST模型的计算结果及全湍流试验数据进行比较.结果表明:远场边界距离对翼型阻力系数有较大的影响;与无转捩模型相比,γ-Reθt转捩模型对翼型阻力系数的预测精度有一定程度的提高;对于二维模型,γ-Reθt转捩模型难以捕捉翼型表面的三维效应和非定常分离特性.     

吸气式高超声速飞行器机体推进一体化技术研究进展

吸气式高超声速一体化飞行器最显著的特点是子系统之间的耦合较其他类型飞行器更加强烈,这使得其设计具有挑战性。所有的子系统之间部件相互干涉,包括:气动、推进、控制、结构、装载和热防护等,特别是机体与超燃冲压发动机之间的耦合最为突出。飞行器的前体和后体下壁面既是主要的气动型面,又是超燃冲压发动机进气道外压缩型面和尾喷管的膨胀型面,在产生推力的同时也产生升力和俯仰力矩。机体与发动机的强耦合作用对飞行器的推力、升力、阻力、俯仰力矩、气动加热、机身冷却、稳定性和控制特性有直接的影响。本文介绍了国内外机体推进一体化技术的研究进展,重点介绍了中国空气动力研究与发展中心(CARDC)的相关研究工作,包括:密切曲锥曲面乘波进气道和基于双激波轴对称基准流场内转式进气道设计方法、独创的大尺度脉冲式燃烧加热风洞一体化飞行器带动力试验技术和高超声速内外流耦合数值模拟技术等。对高速飞行中激波边界层相互干扰、流动分离机理、可压缩湍流转捩及其控制、超燃冲压发动机燃烧流动机理等相关基础问题也进行了研究,强调了对高效高精度计算方法的迫切需求。     

无网格法耦合RNG k-ε湍流模型在亚、跨声速翼型黏性绕流中的数值模拟

基于移动最小二乘无网格方法,耦合RNG(Re-Normalisation Group)k-ε湍流模型求解雷诺平均Navier-Stokes方程。采用AUSM(Advection Upstream Splitting Method)+-up迎风格式求解数值通量,应用在高度各向异性点云结构中取得良好结果的点云重构技术结合移动最小二乘法拟合空间导数,并用三阶SSP(Strong Stability Preserving)型Runge-Kutta显式时间推进格式求解离散后的控制方程。在此基础之上,实现了对NACA0012、RAE2822翼型亚、跨声速黏性绕流的数值模拟,给出了翼型表面压力系数分布曲线、不同位置处的平均速度剖面、马赫数等值线等计算结果,并与实验值及相关文献数值模拟结果进行比较,结果吻合较好。表明所发展的结合点云重构技术的无网格方法耦合RNGk-ε湍流模型能够成功模拟翼型亚、跨声速黏性绕流,验证了所提算法的有效性,并拓展了无网格方法求解湍流流动的途径。