高温高压喷管化学和热力学非平衡流计算

本文给出了高驻室压力条件下的化学、振动、电离非平衡喷管一维流动的方程组（包括Ｖ－Ｄ耦合和电离激发效应）。为了克服高压条件下喉道前空间推进步长小的困难，采用局部平衡等效γ的计算方法，对两类典型风洞（小尺寸电弧风洞和大尺寸激波风洞）进行了数值计算，并对不同驻室压力、Ｖ－Ｄ耦合、电子激发效应、不同化学模型进行了研究。计算表明，高驻室压力能够有效抑制喷管的非平衡流动。在Ｐ０＝２００ＭＰａ，Ｔ０＝８０００Ｋ     

关于一维非定常气─固二相流特征理论的讨论

本文研究了一维非定常气－团二相流的特征，对于目前常用的、包含一定近似的运动方程组，它有两个虚部很小的复根，研究表明这种很小的虚部出现与方程中引入的近似有关。从十多个算例的结果看，这样的虚部不影响人们在工程计算中用该方程组求解初值问题。针对具有上述特征的气－团二相流方程组，本文发展了一种类特征线方法，它保留了双曲型方程特征线方法的若干优点，但它允许方程的特征根中有微小的虚部。含灰气体激波管流动的计算表明，类特征线法具有精度高、数值耗散和弥散小等优点     

跨音速喷管化学反应流数值模拟

用时间相关的半隐格式有限差分数值方法求解了化学非平衡反应跨音速喷管流场，在喷管收敛段，流动接近化学平衡状态，控制方程的刚性问题严重，数值积分困难。通过对时间差分项隐式离散、对空间差分项显式离散，流场边界采用参考平面上的特征线法计算等，成功地解决了由于化学反应有限速率带来的数值解不稳定问题。该格式简单、需要计算机存贮空间少。本文完成了一维和轴对称非平衡化学反应喷管流动计算，并与化学平衡流和冻结流的计算结果做了比较。     

圆管内高超声速流的发展和激波形态的实验研究

本文的研究目的是了解沿圆管内发展的高超声速流的结构。实验是在加拿大多伦多大学的高超声速炮风洞内完成的，风洞的自由射流马赫数Ｍ＿∞＝８．３０，总温Ｔ＿（ｔ∞）＝１０００Ｋ，总压Ｐ＿（ｔ∞）＝２６．５ＭＰａ单位雷诺数Ｒ＿ｅ＝３．２×ｌ０￣７。壁面静压以及内流中的若干截面内的皮托压力和静压测量结果揭示，管内产生的激波主要是斜激波形态，而且存在着较强的激波与边界层的相互干扰。实验发现，近管中心线的高超声速流动有不稳定现象；壁面的边界层，基本上是湍流边界层，特别是干扰区的下游；圆管出口的周向内壁面的顶壁面静压对管的攻角异常敏感。本文提供的结果，可以指导高超声速流的计算流体力学的方法和进展。     

磁场对高超声速弱电离气体流动的影响

对偶极子磁场作用下的三维钝头体高超声速黏性绕流的化学非平衡流动进行了数值模拟。电导率利用组分公式计算,化学模型为7组元、6反应模型。应用诱导磁场方法将磁场与流场耦合,控制方程的空间离散采用有限体积法,扩散项用中心差分格式,对流项采用二阶迎风格式,时间推进采用3步Runge-Kutta法。计算结果表明,外加偶极子磁场使激波脱体距离增加,壁面摩擦系数和表面热流密度增大。与无磁场作用时相比,在0.153T外加磁场作用下,冻结流中的激波脱体距离增加约3倍,局部壁面摩擦系数最大增加63%,局部表面热流密度最大增加61%;而化学非平衡流中的激波脱体距离增加约0.5倍,局部壁面摩擦系数最大增加47%,局部表面热流密度最大增加31%,且化学非平衡流中激波层内温度的最大值约为冻结流中的64%。     

超音速多圆弧(MCA)叶型叶珊的研究

本文提出了含脱体激波的叶栅有旋超音速流场的特征线起始线计算方法。采用激波嵌入法解完全的欧拉方程,计算了含脱体激波的多圆弧叶栅的超音速进口流场。本文发展了一种将超音速流场计算与亚音速扩散损失计算的尾迹法相关联的计算超音速叶栅损失的方法,可分别确定激波损失和亚音速扩散损失。      

热化学非平衡来流条件下热化学模型影响研究

为研究热化学非平衡来流条件下热化学模型等计算设定对斜激波压缩流动计算结果的影响,针对尖劈构型和相应的前缘钝化构型的高焓激波风洞实验,采用多种计算设定开展详细的数值模拟研究。计算结果表明,计算采用不同热化学模型,以及来流设定为振动冻结/平衡/非平衡状态,会导致斜激波激波角等参数存在一定差别,其中激波角差别可达约2%。当来流速度一定时,过斜激波后分子内能增量在平动转动能和振动能上的分配方式的差别决定了激波角的差别。前缘钝化情形下,采用不同计算设定所得激波角之间的关系和尖前缘构型的规律一致;但是,采用不同计算设定所得斜激波到壁面距离之间的关系和尖前缘构型的规律有差别,这源于钝化前缘的激波脱体距离的影响。对于自由来流下的斜激波压缩流动问题,若考虑了分子振动能激发但未考虑热力学非平衡(例如热完全气体模型、考虑空气反应的单温度模型等),就斜激波激波角等参数而言,计算误差比量热完全气体模型计算误差更大。     

高超声速化学非平衡粘性激波层方程的数值解

本文讨论了高超声速粘性激波层方程数值计算时差分格式引起的物理失真问题。具体分析了全隐格式格式粘性的影响,并作了数值试验。为了验证隐式结果的可靠性。在超声速激波风洞中测量了钝锥的表面压力分布,并与计算结果作了比较,两者基本一致。 本文采用隐式有限差分法数值计算了高超声速化学非平衡粘性激波层绕细长球锥的流动。计算时采用连续方程和法向动量方程耦合求解的方法以解决细长体远后身区计算中的问题。应用网格技术和加强系数矩阵主对角元素优势的方法提高了化学非平衡流计算的雷诺数范围。文中给出了高超声速化学非平衡流的计算结果,并与其它文献的结果作了比较。     

用随机选取法（RCM）计算高温平衡气体的激波管流动

本文用RCM方法计算了高温平衡气体的激波管流动。在求解Riemann问题时,激波关系式采用精确解,而对于稀疏波关系式,本文建立了等效比热比的局部理想气体近似关系式,从而使RCM方法大为简化,并具有足够精度。     

均匀磁场中高超声速弱电离气体流动数值模拟

采用7组元化学模型并应用组分公式计算电导率,通过求解黏性MHD(magneto-hydro-dy-namics)方程组,研究了不同强度均匀磁场对三维钝头体高超声速绕流化学非平衡流动的影响.结果表明,随着磁场的增强,激波脱体距离逐渐增加;总阻力系数和壁面温度逐渐减小.在B_y=0.03 T磁场作用下,与无磁场的结果相比,化学非平衡流中的激波脱体距离增加约7%,总阻力系数减小约5%,局部壁面温度最大降低74%;而冻结流中的激波脱体距离增加约43%,总阻力系数减小约6.9%,局部壁面温度最大降低18%.      

燃气射流起始冲击波形成机理的实验研究

利用实验研究的方法对燃气射流起始冲击波的形成机理进行研究。目的在于研究其发生、发展的规律,为寻求降低其危害性的研究提供参考。实验中使用了激光Ｍｏｉｒｅ偏拍仪对起始冲击波场的形成初期进行了光学显示,同时还使用了压阻式压力传感器对起始冲击波超压进行了测量。通过对测量结果的分析,认识了燃气射流起始冲击波的物理本质且总结出其形成和发展的一般规律,它在近场的高峰值超压对火箭武器的承载设备具有破坏作用。     

Precessing motion in stratified radial swirl flow

Vortex/flame interaction is an important mechanism for unsteady combustion in a swirl combustion system. Technology of low emission stirred swirl(Te LESS), which is characterized with stratified swirl flow, has been developed in Beihang University to reduce NOXemission. However,large-scale flow structure would be induced in strong swirl flow. Experiments and computational fluid dynamics(CFD) simulation were carried out to investigate the unsteady flow feature and its mechanism in Te LESS combustor. Hotwire was firstly applied to testing the unsteady flow feature and a distinct mode with 2244 Hz oscillation frequency occurred at the pilot swirl outlet.The flow mode amplitude decayed convectively. Large eddy simulation(LES) was then applied to predicting this flow mode and know about its mechanism. The deviation of mode prediction compared with hotwire test was 0.8%. The spiral isobaric structure in pilot flow passage indicates that precessing vortex core(PVC) existed. The velocity spectrum and phase lag analysis suggest that the periodic movement at the pilot outlet was dominated by precessing movement. Negative tangential momentum gradient reflects that the swirl flow was unstable. Another phenomenon was found out that the PVC movement was intermittently rotated along the symmetric axis.     

二维涡轮级的定常/非定常数值模拟

为了解定常与非定常计算差别对工程应用参考价值的影响，本文通过求争二维Navier-Stokes方程，针对一级涡轮叶栅分别进行了定常／非常数值模拟。结果发现，对于叶片表面压力分布，第一排导叶定常计算与非定常计算得到一致的结果，但是在动叶上差别较大，尤其在叶排小间距时。当叶排间距拉大后，定常计算结果趋近于非定常计算结果时间平均值。     

单级环境静叶等离子体流动控制机理仿真研究

在轴流压气机等离子体扩稳研究中,针对单转子压气机流动控制的研究较多,而针对单级环境下静叶流动控制的研究却很少.采用静叶轮毂轴向等离子体激励方式,通过数值模拟方法研究单级环境下静叶流场特性,揭示轴流压气机静叶等离子体流动控制扩稳机理.结果表明:等离子体激励器的轴向位置对单级轴流压气机的扩稳效果影响显著,越靠近叶片前缘,扩稳效果越好;布置在静叶通道后半部的等离子体激励器无法提高压气机的稳定性,而在静叶前缘施加轴向等离子体激励时,近轮毂区气流被诱导加速,主流的轴向速度提高,有效抑制了静叶近轮毂区吸力面的流动分离,静叶近轮毂区的堵塞减小,使得单级轴流压气机的稳定性提高.     

吸力面附面层抽吸在三维高负荷扩压叶栅中的作用机制

对某矩形高负荷扩压叶栅在不同弦向位置开设全叶高抽吸槽的5组方案进行了数值研究,分析了抽吸槽弦向位置等参数对抽吸量分布规律的影响;通过叶栅实验探究了局部展向抽吸方案的效果.数值仿真的计算域包含吸附叶片内部的真空腔,边界条件按照实验条件设置.研究发现:全叶高抽吸方案的抽吸量沿展向大致呈C型分布;叶高中部和端部的主要抽吸效果都体现在叶高中部流场,端部的抽吸量对叶栅角区的回流有一定的抑制效果.抽吸量沿展向的分布规律受叶栅流道和叶片内腔流场的共同作用,因此应根据三维高负荷扩压叶栅流场的具体特性对吸力面抽吸槽/孔进行细化设计.      

高超声速主流中横向喷流干扰非定常特性研究

脉冲发动机产生的横向喷流往往伴随着非定常效应,由于其持续时间非常短,给流场模拟和测量带来很大的困难.目前在地面设备中难以进行精确的试验,本文采用数值模拟方法对横向喷流的非定常效应进行了研究.重点对发动机启动和关车过程中典型时刻的流场进行了细致的刻画,分析了流动参数和模型气动力系数的动态变化过程,并给出了喷流压比和攻角变化所带来的影响.研究表明,喷流启动和关闭过程中,飞行器所呈现出的动态气动特性与稳态喷流时有着明显的差异,尤其是压心位置和法向力系数将发生显著的变化.在飞行器设计中,必须从安全控制的角度出发,考虑这一非定常效应的影响.     

多级压气机非定常流动的数值模拟

进行了多级压气机动静叶片排干扰的三维非定常无粘流动数值模拟研究，采用高分辨率格式求解非定常Euler方程组。对某以级压气机的非定常无粘流场的计算结果表明，即使在不考虑粘性影响的情况下，压气机动静叶排干扰的非定常效应仍然是比较显著的，各叶片排所受到的非定常影响是不相同的，计算得到的叶片表面压力周期性波动的结果可为叶片强度计算提供较为准确的数据。     

弯曲管道内湍流流动的数值模拟

为了对弯曲管道内的湍流流动特性进行有效的预测，用显示代数应力模式对－典型的弯曲管道内流动进行了数值模拟和分析，计算结果表明：（1）显式代数应力模式能对弯曲管道内流动的平均量（如速度、压力、摩阻等）进行较为准确的预测；（2）对模式中的ε方程做适当修正后，该模式能较好地模拟出流动曲率效尖对湍流特性的影响，是对雷诺应力模式的一个比较好的近似。此外，该模式在形式上与标准k－ε两方程模式类似，使用方便，具有较强的工程实用价值。     

气液两相分层流剪切应力的不确定度分析

基于水平气液两相分层流的压力梯度、波状液层和气壁剪切应力实验数据,对两流体剪切应力进行不确定度计算,分析各不确定度分量对液壁和气液界面剪切应力结果的影响。由剪切应力可以直接计算分层流动量平衡方程的本构参数,对本构参数的关联作了合理的解释。气液界面剪切应力受液层高度和压力梯度的准确性影响较大,而液壁剪切应力却相反,和气液界面剪切应力最大程度相关。通过和实验数据对比发现,液壁摩擦因子无法像气壁摩擦因子那样可以用单相管流的指数型关系式来描述,而是应该结合气液界面摩擦因子建立一个基于两相各流体特征参数的有效关联式。     

圆环管道中曲率对流动稳定性的影响

为了研究圆环管道中曲率对流动稳定性的影响,选取了展向带有曲率的轴向流动和轴向带有曲率的周向流动两种圆环管道模型进行分析.采用线性稳定性理论分别计算了两种模型下不同曲率的稳定性和扰动演化情况,并和槽道流进行比较,从而得到了曲率对基本流动、特征值和特征函数的影响.结果表明：曲率越大,基本流向内侧偏转越大,扰动波流向波数越大,增长率越小,中性曲线的不稳定区域越小,特征函数的对称性越差.因此曲率对边界层内的扰动起着稳定的作用,并且在一定范围内曲率越大,稳定性越好.