液晶显示方法研究管内马赫数4拟似冲击波流动

在一个压力-真空型超声速风洞中, 剪切应力敏感的液晶表面流动显示技术被用来研究方管内马赫数4拟似冲击波的超声速流动.实验是在日本室兰工业大学的超声速风洞中进行.它提供了关于测试表面流动结构的定性信息,诸如湍流边界层分离和再附, 以及边界层流动的维数等.液晶显示结果与纹影实验照片和油膜方法进行比较,证实了它是一种非常有效的表面流动显示工具.     

马赫数2.5超/亚声流流向涡掺混超燃室冷态实验

结合对应的数值计算,对一种带亚声速预燃室和流向涡掺混器的超声速燃烧模型燃烧室实验台,在其进口马赫数为2.5的来流条件下,进行了冷态流场的实验研究.实验测得其壁面静压分布和激波系结构与流场的CFD计算结果基本一致.实验结果表明,模型燃烧室全流场超声速,达到设计状态.马赫数2.5下的冷态实验数据和CFD计算数据为进行点火实验提供了依据.      

液晶流动可视化方法研究拟似冲击波的内部超声速流

在一个压力一真空超声速风洞中,剪切应力敏感液晶流动可视化技术被应用来研究方管内马赫数2拟似冲击波（pseudo-shockwave）的超声速流动。它主要提供关于整个流动的定性信息,诸如湍流边界层分离、再附着位置以及流动的维数等。而且液晶也反映了表面流线,分离区内的涡流和管道流动的角效应。使用两种不同黏度的液晶分别进行实验,分析黏度对结果的影响。液晶实验的结果与纹影照片所得结果比较吻合,说明了液晶是一个非常有效的流动可视化工具。     

导流格栅对发射箱内冲击波强度影响仿真

针对发射箱内添加导流格栅导致箱内冲击波强度变化的问题,构建了燃气流与后盖开启过程相互耦合的数值仿真模型,并基于动网格技术进行了仿真计算。研究结果表明:添加导流格栅后,箱内及前盖冲击波强度均有所增加,前盖压强峰值增加约 5%,导流格栅安装位置越靠近发动机喷管,箱内冲击波强度增幅越大。研究结果可为发 射箱内导流格栅设计提供理论参考。     

CFX-5．5对GBM-01标模气动力计算应用研究

本文介绍了采用通用CFD数值模拟软件CFX-5．5对GBM-01标模气动力特性的数值计算过程和计算结果。计算马赫数为0．6-2．024，计算结果与实验结果具有相同的趋势．文中对计算结果与实验结果的差异进行了分析。本文研究工作可为CFX-5．5在工程中的实际应用提供参考依据。     

升力体式密闭救生舱无控制再入研究

将救生舱选定为升力体模型,利用数值模拟的方法动态研究救生舱再入过程中的稳定性和减速性,分析无主动控制的救生舱返回的可行性。针对高超声速和马赫数小于4两种工况,分别采用牛顿理论方法和CFD方法进行气动力计算;利用气动特性计算结果求解六自由度轨道仿真动力学方程进行轨迹计算。研究结果显示,救生舱初始再入条件合适时,其再入过程具有稳定性;仅依靠目前救生舱外形的结构减速不能满足减速性要求,需辅助其他减速措施。     

液体煤油超声速燃烧数值分析的湍流模型

利用计算流体力学软件对煤油在所设计的双模态燃烧室内的喷雾燃烧进行了数值模拟。采用k-ε三种形式、k-ω两种形式和雷诺应力湍流模型,计算了在飞行马赫数为5,煤油与空气的当量比为0.551时的情况。通过将数值模拟结果与实验数据进行对比,证明这些湍流模型都可用于煤油超声速燃烧的数值计算。并且,k-εRNG模型和k-ωSST模型的计算结果与实验数据最为接近,可以认为这两个模型在模拟煤油超声速燃烧流场中是非常有效的;而雷诺应力模型效果较差。     

基于定常N-S方程的导弹滚转数值模拟研究

基于商业CFD软件FLUENT定常N-S方程,在不同计算资源下数值模拟了某标准导弹在不同飞行马赫数、滚转速度和迎角时的滚转特性,并与AEDC报告实验结果进行了对比。计算结果验证了定常算法数值模拟滚转阻尼系数的可行性及其在不同计算资源下计算精度的可优化性,而且定常算法与非定常计算相比消耗更少的计算资源。     

汽车风洞支撑干扰扣除方法研究

借鉴航空风洞镜像法原理,提出一种针对汽车风洞支撑气动干扰扣除的方法。在汽车风洞中完成汽车模型和天平支撑连接或分离的两次风洞试验,得到了汽车模型的气动力和支撑对汽车模型的干扰力之和。为了扣除模型支撑对汽车模型气动力测量的影响,在进行风洞试验的同时,应用CFD软件进行工况完全一致的数值仿真,计算模型支撑对不同车型气动阻力的影响。最后,通过归纳的修正公式将数值仿真获得的影响数值转化为汽车风洞试验的影响数值,获得汽车风洞试验的最终阻力系数结果。结合风洞试验和数值仿真的研究结果表明两种方法互相验证,互相补充,可以解决汽车风洞试验时支撑干扰扣除的实际工程问题。     

超声速横喷干扰湍流场数值模拟

本文探讨了超声速平板横向喷流干扰湍流流场的数值模拟方法。分别采用Roe、Harten-Yee及AUSM 格式对三维雷诺平均Navier-Stokes方程进行离散,湍流模拟应用Baldwin-Lomax代数模型,讨论了B-L模型中Fy最大值的不同选取方式。分析了数值格式、网格划分对计算结果的影响。将来流马赫数为5、喷流马赫数为3的实验结果与计算结果进行比较,选择了与实验符合较好的数值模拟方法,为进一步利用数值计算方法研究喷流干扰流场特性奠定基础。     

波转子非定常泄漏流动机理

针对制约波转子性能的泄漏问题,通过提取波转子中与泄漏相关的主要流动现象并建立波转子单通道泄漏模型,对非定常泄漏流动机制进行了详细的数值研究,对本文给出的非定常泄漏损失预测模型进行了数值验证。结果表明:转子通道中存在不同程度的压力波动,波动幅值与间隙宽度有关;连续反射膨胀波、周期性出现的弓形激波及其反射激波是压力波动的根本原因;间隙内部泄漏过程存在3个主要的流动阶段;泄漏过程中通道激波传播速度不变、波后时均压力不变;在一定间隙宽度范围内,激波马赫数、激波静增压比与无量纲间隙宽度均呈线性关系,当间隙宽度从0增大到0.08时,激波马赫数衰减7.3%,激波静增压比衰减10.1%;泄漏流动通过泄漏产生的主膨胀波对激波传播过程施加影响,通道激波衰减本质上是理想激波与主膨胀波叠加效应的结果;泄漏损失预测模型与数值结果吻合良好。     

密切曲面内锥乘波前体进气道设计和试验研究

介绍了密切曲面内锥乘波前体进气道（Osculating Inward turning Cone Waverider Inlet,OICWI）的一体化设计方法,对该型乘波前体进气道的性能进行了数值分析,针对该型一体化乘波前体进气道完成了风洞试验研究。理论设计结果和设计状态无粘模拟结果一致,设计状态下的计算结果表明,前体进气道具有较高的总压恢复、较好出口流场均匀度及较高的流量捕获率。试验研究结果表明,改型一体化前体进气道在马赫数5～7条件下顺利启动,流场波系及压力分布同数值分析结果吻合。     

基于LES方法的平板非定常激波/湍流边界层干扰研究

以高超声速发动机进气道湍流分离控制为应用背景,采用大涡模拟(LES)方法进行马赫数为3.0(唇口附近马赫数约为3.0)的激波/湍流边界层干扰(SWTBLI)流场机理研究.利用扰动循环引入的方法,先得到充分发展湍流场,然后根据斜激波关系式引入激波的方法进行激波/湍流干扰模拟.研究结果显示:充分发展湍流场在激波作用下产生逆...     

一种二元高超声速进气道起动特性的尺度效应研究

对一种二元混压式进气道三维流场进行了数值和实验考察,研究了不同尺度进气道模型自起动性能的变化。结果表明,在相同来流单位雷诺数条件下,随着模型尺度的减小,进气道自起动马赫数有所提高,起动性能有所降低。同时对不同尺度模型进行雷诺数匹配,发现在相同雷诺数下,不同尺度模型的起动性能相近,表明雷诺数是影响不同缩尺模型起动性能不同的主要原因。在可获得的实验结果范围内,数值模拟所得到的自起动结果基本与之相符。此外,对实验中发现在低雷诺数下进气道反而呈现出自起动特征的异常现象进行了初步分析,通过数值模拟比较指出了低雷诺数下来流偏向层流流态,可能会导致进气道呈现一种“起动”状态。     

用CSCM格式模拟二维收缩管道内的流动

 采用Lombard等人提出的CSCM(Conservative Supra-Characteristics Method)差分方法数值模拟二维收缩管道内波系干扰及激波边界层干扰等问题。所得结果同已往计算结果进行了比较。通过压力场和速度场的计算结果可看到Ma数对流场的影响。在Ma数较小时激波角较大,可看到激波经两次相交一次反射逐渐变弱的情况。Re数改变对流场也有明显影响,且给出由于激波干扰引起边界层分离的计算结果。最后还给出收缩管道内激波由正激波逐步形成斜激波的非定常过程。     

欠膨胀椭圆射流的流动结构

为了探讨超声速欠膨胀椭圆射流的流动特征,采用大涡模拟（LES）方法与高精度混合格式对出口压力比N分别为14、24、40的欠膨胀椭圆射流流场结构进行数值模拟。结果清晰描述了欠膨胀椭圆射流的三维(3 D)结构特征与发展规律,并分析了因喷嘴方位曲率不一致而导致长轴与短轴平面上激波结构出现差异的原因。另外,结果还发现:当N为24时,短轴平面上射流域内的激波结构已由正规反射转变为马赫反射,但长轴平面上仍维持正规反射,而当N为40时,长轴与短轴平面上的激波结构均为马赫反射结构,由此可知喷嘴方位曲率变化越平缓,马赫反射形成所需的出口压力比越大。      

斜激波极值规律的边界层影响

采用边界层理论与斜激波/膨胀波精确算法，建立一种结合Eckert参考温度法和Illingworth-Stewartson变换法优势的边界层权重算法，用于研究超声速黏性楔面边界层位移厚度对斜激波极值规律的影响。分别应用层流Navier-Stokes方程和湍流Navier-Stokes方程的CFD解算器对边界层新模型进行了算例精度评估。在来流马赫数为1.2~2.4和楔面角为3°~20°的范围内，压强比的相对误差小于0.1%。计入层流与湍流边界层影响的理论模型研究表明，边界层影响使得最优马赫数增加；对于层流边界层，最优马赫数增量约为0.001 5~0.003 3；对于湍流边界层，最优马赫数增量约为0.002 8~0.006 1。     

氢燃料双模态燃烧室模态转换

为研究双模态燃烧室的模态转换规律,采用计算流体动力学(CFD)方法,对氢燃料模型燃烧室在不同当量油气比、飞行马赫数、燃料喷射方式下的流场进行了数值模拟,并与试验结果进行对比,两者相互吻合.研究结果表明:当量油气比提高、飞行马赫数降低及双面喷射燃料均使燃烧室更趋于亚燃工作模态.     

CFD/CSD-based flutter prediction method for experimental models in a transonic wind tunnel with porous wall

To predict the flutter dynamic pressure of a wind tunnel model before flutter test, an accurate Computational Fluid Dynamics/Computational Structural Dynamics (CFD/CSD)-based flutter prediction method is proposed under the conditions of a 2.4 m × 2.4 m transonic wind tunnel with porous wall. From the CFD simulations of the flows through an inclined hole of this wind tunnel, the Nambu’s linear porous wall model between the flow rate and the differential pressure is extended to the porous wall with inclined holes, so that the porous wall can be conveniently modeled as a boundary condition. According to the flutter testing approach for the current wind tunnel, the steady CFD calculation is conducted to achieve the required inlet Mach number. A time-domain CFD/CSD method is then employed to evaluate the structural response of the experimental model, and the critical flutter point is obtained by increasing the dynamic pressure step by step at a fixed Mach number. The present method is applied to the flutter calculations for a vertical tail model and an aircraft model tested in the current transonic wind tunnel. For both models, the computed flutter characteristics agree well with the experimental results.     

矩形转圆形高超声速内收缩进气道数值及试验研究

采用压力梯度先增大后减小压升规律轴对称基准流场,结合流线追踪及截面渐变技术设计了矩形转圆形内收缩进气道模型,并采用4°斜楔模拟飞行器前体,对前体、进气道一体化模型进行了数值模拟和风洞试验,初步得到了该进气道的流场结构及总体性能.设计点和接力点的数值模拟结果表明该进气道可在马赫数Ma=4～6状态下正常工作,且具有良好的总...