各向异性网格自适应计算在超燃模拟中的应用

采用各向异性网格自适应求解技术,将其应用于DLR超燃冲压发动机燃烧室中的超燃模拟.开展了3个算例,包括采用滑移壁面条件的燃烧室冷流场模拟、采用无滑移壁面条件的冷流场模拟及采用无滑移壁面条件的反应流模拟.模拟中,各向异性网格自适应计算捕捉到了如激波、射流、边界层、火焰面等具有各向异性特征的大梯度区域,并利用各向异性网格进行了很好的加密.对比利用各向同性网格的初始流场计算,各向异性网格自适应计算使基于滑移条件、无滑移条件的冷流计算及反应流计算的网格单元数量分别下降了36.2%,36.4%和36.8%,有效降低了计算规模,而且流场大梯度区域的计算结果更准确,辨析度更好.结果表明:对于像超燃这类具有各向异性特征的问题,各向异性网格系统比各向同性网格系统有更好的计算效率及准确性,同时也表明基于Mach数场构造的各向异性网格系统可以有效应用于超燃计算.      

绕圆头回转体通气空化流型的实验研究

为了深入了解通气空化流动现象,利用高速全流场显示技术,对绕圆头回转体通气空化流型进行实验研究。结果表明,重力效应和通气量对通气空化的多相流流型起主要作用。定义了弗洛德数和通气率两个无量纲数,将绕圆头回转体通气空化分为5种多相流流型,即透明空泡、透明气弹、透明分层、水气混合以及半透明水气混合。流动参数对流型的影响分为2个阶段,即重力起主要作用阶段和重力效应不明显阶段。在重力起主要作用阶段,通气率一定时,随着弗洛德数的增大,附着弹体的空泡倾斜程度变小,弹体上表面的断裂空泡转变为贴着弹体壁面的稳定空泡;弗洛德数一定时,随着通气率的增大弹体上表面断裂空泡的尺度不断增大。在重力效应不明显阶段,通气率一定时,随着弗洛德数的增大,雷诺数变大,流场的湍流强度增大,空泡尾流区域水气交换的程度加剧;弗洛德数一定时,随着通气率的增大,通气空化数减小,绕弹体的云雾状空泡逐渐转变为透明空泡。最后,进一步分析了重力影响下透明空泡脱落的非定常过程,以及反向射流作用下云雾状空泡交替脱落的非定常过程。     

(高)超声速流动试验技术及研究进展

近年来,与高速飞行器相关的(高)超声速流动受到了极大的关注。这类流动所具有的非定常性、强梯度和可压缩性对试验方法和风洞设计技术提出了挑战。超声速纳米示踪平面激光散射(NPLS)技术是由作者所在团队研发的非接触光学测试技术。它能够以较高的空间分辨率来揭示超声速三维流场的一个瞬态剖面的时间解析的流动结构。介绍了NPLS技术以及基于NPLS开发的密度场测量、雷诺应力测量和气动光学波前测量等方法,并回顾了这些技术在超声速边界层、超声速混合层、超声速压缩拐角、激波/边界层相互作用和光学头罩绕流等流动中的应用,清晰地再现了边界层、混合层、激波等典型流场结构及其时空演化特性。另外,为了模拟和研究高空大气条件下边界层自然转捩和超声速混合层的转捩特性,介绍了高超声速静风洞、超-超混合层风洞的设计技术以及层流化喷管的设计方法。     

间隙高度对超声速膨胀器流场及性能的影响

采用三维雷诺平均N-S方程和标准k-ε湍流模型,对不同间隙高度超声速膨胀器的流场和性能进行了数值研究,结果表明:间隙高度显著影响三维流道内的局部流动特性和超声速膨胀器的整体性能,随间隙高度增加,气流最高相对马赫数降低,高速区范围逐步缩小;泄漏涡增强,尺度变大,横向和径向运动明显,泄漏损失增加,但激波及激波附面层相互作用的损失降低;超声速膨胀器的膨胀比先增大后减小,等熵绝热效率持续降低.下端壁、吸力面附近低能流体之间以及与壁面的摩擦损失和间隙泄漏损失是有间隙超声速膨胀器三维流道内损失的主要来源,超声速膨胀器的间隙高度宜在0.9%h₀~1.5%h₀之间选取.      

变循环发动机组合变几何调节方案

基于面向对象的设计思想,设计出了一个双涵道变循环发动机(VCE)的性能计算模型,并选取了双涵道VCE具有代表性的典型工作点,进行了双涵道VCE的组合变几何调节的相关研究.结果表明:在双涵道VCE的亚声速巡航和超声速巡航工作点,合理调节风扇导流叶片角度可以更大程度发挥双涵道VCE性能优势;在亚声速巡航工作点时选取的方案4和超声速巡航工作点时选取的方案4均比双涵道VCE设计点的变几何调节方案性能更优.     

一类新的速度入流条件在管道流模拟的应用

基于经典边界层理论,发展了一类简化的、按边界层厚度描述的速度入流边界条件。基于入流条件,对低速后向台阶流动和超声速压缩拐角流动进行了数值模拟。通过细致比较空间流场、边界层内速度型、表面摩擦阻力和压力分布等特性,对该类型边界条件进行了计算确认。结果表明,所发展的简化边界层厚度速度入流条件提法正确、合理,具有描述方法简单、鲁棒性好、适合工程计算的优点,可显著简化不同流速的管道流动的数值模拟。     

气动热环境下多翼导弹温度场数值模拟

导弹在飞行过程中与大气摩擦产生的热量,会降低材料的强度极限和飞行器结构的承载能力,使结构产生热变形,破坏部件的气动外形并影响飞行器的安全飞行,带来安全隐患。通过气动加热工程算法确定表面热流密度,利用有限元分析软件MSC Nastran和MSC Patran,对四翼导弹进行有限元计算,给出温度场分布情况,为导弹的设计和制造提供依据。     

亚/超声速楔状流层流边界层速度与温度相似解及拟合解

利用相似变换获得了楔状流层流边界层无量纲流函数的3阶非线性常微分方程,用Runge-Kutta法求解微分方程获得了不同楔形角楔状流层流边界层无量纲速度随相似变量的变化曲线;推导了亚声速和超声速楔状流层流边界层无量纲温度关于相似变量的2阶线性齐次和非齐次微分方程,获得了温度分布的通解,恒壁温条件下亚声速楔状流和绝热壁面条件下超声速楔状流层流边界层无量纲温度解析解及指数函数形式的拟合解.以楔形角为0为例利用相似变换研究了超声速条件下气体压缩性及黏度随温度变化等因素对层流边界层速度与温度的影响,得出不可压缩常物性与可压缩变物性条件下无量纲速度相对误差绝对值小于9.8%的结论.研究表明:Pr越大贴近壁面处无量纲温度变化越剧烈;超声速条件下壁温低于绝热壁温时黏性耗散作用可以使层流边界层气体温度从壁面到主流间出现先升高后降低的变化.      

宽裕度超声叶型气动优化设计

通过与德国航天局设计的超声预压缩叶型PAV-15试验数据比对,确定高精度超声叶栅流场计算方法,研究表明:根据激波位置分段调整流管厚度可提高计算与试验结果的一致性。为提高超声叶栅稳定工作裕度并保证设计点性能,建立根据目标裕度估算喘点反压方法和优化设计方法。对两个超声叶型进行多目标优化,优化结果表明：优化叶栅可减小设计工况槽道激波入射角、减小激波及激波附面层干扰损失;气动喉道前移、结尾正激波后移,提高叶栅耐反压能力。两个优化叶型在保持总静压比不变的前提下,稳定裕度均达到设计目标,设计点损失也有所下降。     

超声速进气道压力估算方法及验证

进气道载荷的预示和限制是超声速飞行器设计中的关键问题。以典型颌下进气超声速飞行器为研究对象,对其进气道流场进行数值仿真,研究不同马赫数、攻角、侧滑角及余气系数条件下的进气道压力特性;针 对进气道压力工程估算及设计需求,使用无量纲和解耦的方法,对进气道压力经验公式进行拟合;反算飞行试 验中的进气道压力,并与测量结果进行对比。结果表明:进气道压力随马赫数增大而增大,随余气系数增大而 减小;正常工作包线内,较小的攻角、侧滑角对进气道压力影响不明显;进气道压力经验公式计算值与飞行试验 测量值符合较好,具有较高的精度。