高压捕获翼构型跨声速流动特性初步研究

高压捕获翼是一种可以在高超声速条件下同时获得高升阻比、高容积率和高升力系数的新型布局概念。为初步分析该类新型布局的宽速域气动特性,以一种圆锥-圆台体组合高压捕获翼原理性构型为计算模型,以典型跨声速条件(马赫数0.92、0°攻角)为计算工况,进行了计算和分析。流场分析结果表明,在跨声速流动条件下,机体与高压捕获翼之间存在比较强烈的气动干扰,且干扰的强烈程度与高压捕获翼-机体间的垂向距离大小直接相关。与不带捕获翼的参考构型相比,增加捕获翼会导致机体尾部分离区范围增大,并在机体与捕获翼之间的开放通道内形成类似于变截面收缩管的流动,致使沿流向方向出现了明显的激波串,进而导致捕获翼下表面壁面压力出现较为明显的波动。同时,由于机体和捕获翼间的垂直距离沿展向方向逐渐增加,导致该波动在对称面附近最为剧烈,然后随展向位置逐渐增加,压力波动逐渐减弱。     

火星大气高温光谱建模与非平衡辐射热流预测

火星科学实验室成功着陆后的热环境重构数据表明,气动辐射加热在火星进入防热设计中具有不同于以往认识的重要影响,未知机制和模型不确定性等问题有待进一步研究。探测器高速进入火星大气产生极高温非平衡气动环境,造成火星气动辐射与常规CO₂红外辐射研究显著不同。针对火星大气高温光谱和辐射热流预测,首先,建立适用于火星大气的高温非平衡光谱辐射模型,获得典型高温条件CO₂光谱结构和辐射强度,与NASA和JAXA试验结果对比,结果显示符合较好。其次,依靠激波管和发射光谱测量技术,开展典型进入条件的辐射强度测量试验,数值结果、试验值和NASA试验结果相互符合较好,验证了光谱模型。最后,对探路者号进行气动辐射加热分析,完成典型进入条件下的非平衡流动和辐射特性计算,基于光线法得到光谱辐射强度沿驻点线的变化,表明高速与低速条件的气动辐射机制存在显著差异;基于有限体积法获得进入器表面辐射热流分布,结果显示辐射热流的分布及变化规律与地球再入显著不同,进入速度6 km/s以下时辐射热流随进入速度增加而减小,同时进入器锥身及肩部的辐射热流高于驻点区域。     

同轴旋转倒置圆台环隙间流动特性分析

针对同轴旋转倒置圆台环隙间流体复杂流动问题, 对其环隙间流动特性进行了实验研究。重点进行染色液流动显示实验和PIV流场测速实验, 对实验结果做定性及定量分析, 研究内筒转速和环隙宽度对环隙间流动特性的影响。染色液流动显示实验和PIV流场测速实验分别定性和定量地展现了环隙间螺旋涡的产生及变化过程。对不同内筒转速和环隙宽度下的螺旋涡涡心运动周期进行分析, 结果表明, 内筒转速升高, 周期减小;环隙宽度增大, 周期增大。运用瞬时流动和时均流场解析了环隙间螺旋涡运动产生机制, 探究内筒转速和环隙宽度对3种雷诺应力大小的影响与分布情况。内筒转速变化, 雷诺径向正应力始终最大;环隙宽度变化, 雷诺切应力始终最小。      

高稳定性电弧等离子体超声速喷管流动装置

笔者在已研制的壁稳电弧等离子体发生器基础上,改进了阳极结构、加置Laval喷管,建立了一套可长期连续运行的高稳定性电弧等离子体超声速喷管流动装置,超声速电弧射流的不稳定度为1％(纯氩)和3％(氮氩),电弧射流具有很好的稳定性、再现性和对称性。测量表明：超声速电弧射流的电子温度沿流动方向呈现出波浪形周期升降的变化规律。     

梯形截面微管道内流场的PIV测量

利用micro-PIV系统测量了去离子水在水力直径为237μm,长度为31mm的梯形截面微管道内的流场结构,得到了不同雷诺数下沿流动方向不同位置微管道中间截面上的速度分布,并利用数值模拟方法计算了相应雷诺数下三维梯形截面微管道内的速度分布。实验结果和数值模拟结果显示实验中的梯形截面微管道内从层流到湍流的转捩发生在Re=1500～1800左右,实验中充分发展段长度Le。微管道水力直径Dk和雷诺数之间的关系可由公式Le／Dk=（0．08～0．085）船表示。     

二维多介质流动问题的界面处理方法

就二维可压缩多介质流动问题的数值模拟,给出了一种新的界面处理方法。通过在界面处构造Riemann问题,利用Riemann问题的解分别定义界面两边流体的边界条件,由于Riemann问题的解准确地描述了界面处流体的流动状态,因此得到了更加准确的界面边界条件。本文将由Riemann问题的解得到的界面速度外推到整个流场,重新定义速度场,避免了由于速度的大梯度变化而导致的Level—Set等值线相互交错,因而得到了更加精确的界面位置。利用该方法对水下激波与柱型气泡相互作用问题进行数值模拟,结果表明该方法能准确地捕捉各种物理现象。     

二维翼型非定常运动的涡流场显示--2俯仰运动和沉浮/俯仰联合运动

介绍了在水洞中进行的二维机翼纯俯仰运动和沉浮/俯仰联合运动时的流动显示实验.实验结果揭示了机翼做纯俯仰运动时尾迹涡街的特性以及与Strouhal数(Ste)的关系.在Ste=0.105附近,涡街排列成一条线;大于该值,为逆Karman涡街;小于该值,为正Karman涡街.机翼做沉浮/俯仰联合运动时,尾迹中总是出现逆Karman涡街.结合数值计算结果说明,机翼做沉浮/俯仰联合运动时,推力及最大有效迎角与沉浮和俯仰运动之间的相位差有关.并讨论了前缘涡与后缘涡相互作用对推力的影响.     

离心泵叶轮内含颗粒的幂律流体的密相两相湍流研究

阐述了幂律流体控制方程的特点，首先对单相的幂律流体与密相液固两相流进行了计算，计算结果与实验数据吻合较好。其次对泵离心叶轮内的含颗粒的幂律流体的两相流动进行了计算，从进口到出口，幂律流体的速度是减小的，压力随半径增大逐渐升高，而且压力面上的压力大于吸力面上的压力。吸力面上的颗粒拟温度小于压力面上的颗粒拟温度，而且靠近壁面处的大于叶道中的。最后比较了液固两相流与含颗粒的幂律流体的两相流的流动。流体相的紊动能有一定的差别，而且颗粒拟温度也有较大差别。     

过渡型凹腔超声速燃烧流动振荡抑制方法

针对长深比为10.0的过渡型凹腔在隔离段入口马赫数为3.0条件下存在的冷流自激振荡现象,提出了一种凹腔内增加肋条抑制振荡的方案。通过试验和数值计算,对该方案抑制振荡的效果进行了检验,并分析了肋条增加前后燃烧室流场结构和燃烧性能的差别。研究发现：通过在凹腔内增加肋条能够消除过渡型凹腔冷流工况下存在的175.8 Hz的自激振荡,燃烧流场也更加稳定;增加肋条后凹腔的稳焰能力有所降低,部分在凹腔未完全燃烧的煤油进入扩张段后继续发生反应,从而使燃烧区向下游延伸、增大,发动机的燃烧效率和净推力分别降低5.4%和8.9%,但推力更加平稳;燃烧室一维平均热流密度峰值由2.9 MW/m²降低至1.8 MW/m²,燃烧室的热环境大幅改善。     

固体火箭发动机燃烧室三维流动数值计算

采用ｋε双方程湍流模型,以ＳＩＭＰＬＥ 计算程式求解翼柱型装药固体火箭发动机燃烧室内三维非定常不可压流ＮＳ 方程．在建立翼柱型装药简化模型的基础上采用边界标志法来表达燃面推移,采用从二维到三维的初场给定方法结合多重网格法求出了多个时间步的非稳态流场结构．计算结果表明燃烧室内旋涡运动呈现一定的空间与时间分布,周向加质导致燃气通道横断面上轴向速度分布的极不均衡．