激波边界层的相互作用对扰动波传播的影响

利用线性稳定性理论和直接数值模拟研究了带有入射斜激波的、来流马赫数Ma=4．5条件下的平板边界层的失稳特性。重点考察了在由于激波边界层相互干扰，平板边界层上形成分离区，又进而产生激波、膨胀波和旋涡等复杂流动现象的流场上游，引入小扰动的TS波后，扰动波传播通过带有这些复杂流动现象的流场时，扰动波的发展变化特点。通过对流场中扰动波（包括基本波和衍生波）演化特征的分析，研究分离和激波等复杂流动现象对平板边界层稳定性的影响特点。数值模拟发现，激波的出现不同于一般的压缩波，在亚声速区与超声速区对扰动波演化的影响是不同的，此外，分离对扰动有稳定作用。     

超声速来流与侧向射流作用下的三维流场结构

研究超声速流动与侧向射流作用下的三维流场结构对飞行器表面流体分离控制以及超燃发动机燃料混合等具有重要作用。基于大涡模拟(LES)方法和高阶TCD/WENO混合格式,对来流马赫数Ma=4,平板上侧向射流的流场结构特性进行了数值模拟。模拟结果表明,流场内的主要波系由弓形激波、分离诱导的激波以及桶形激波组成,这与相关实验结果相符。同时,计算结果还清晰显示了射流与主流相互作用下,射流上游的分离区流动呈现高度瞬态特性,分离区内旋涡结构不断变化,数目由2个到4个再到6个不断发展,平板表面流线复杂。射流尾涡摆动,涡管上扬,呈现非定常特性。     

跨声速压气机转子的二次流旋涡结构

为了明确跨声速轴流压气机内部流场结构,数值模拟了NASA Rotor37转子,结合λ2准则分析流场参数,探索流动的规律和旋涡结构。研究发现,压气机转子的旋涡模型主要由马蹄涡、壁角涡、径向涡、脱落涡、泄漏涡、诱导涡和分离涡等7个旋涡组成。马蹄涡吸力面分支耗散,压力面分支向相邻的吸力面发展。壁角涡与脱落涡位于叶根角区,引起流动损失和角区失速。径向涡位于激波后吸力面的分离区内,它扩大吸力面分离、引起低能流体向叶顶堆积。激波与叶尖泄漏在叶顶通道中形成3涡:泄漏涡、诱导涡和分离涡,而叶栅通道出口存在分离涡和由泄漏涡与诱导涡合成的叶顶通道涡。泄漏涡与诱导涡破碎在流道中间产生的堵塞区,分离涡造成吸力面尾缘的低速区,共同触发跨声速压气机的失稳。     

连续双扫掠激波/湍流边界层干扰流动特性研究

为了探究连续双扫掠激波/湍流边界层干扰的流动特性,采用仿真方法对一双尖鳍/平板物理模型进行研究。结果表明:双扫掠激波/湍流边界层干扰形成的两个干扰区存在明显的相干现象,虽然第一道扫掠激波/边界层干扰流动仍具有典型的准锥形相似特性,但受其干扰所形成的非均匀流的影响,第二道扫掠激波/边界层干扰却不再具有准锥形相似特性,同时第二个干扰区将影响其上游临近气流的运动甚至影响第一个干扰区的再附线和分离线等结构。两个干扰区形成各自的λ波结构,并且沿着流向两个干扰区内的激波结构相互汇聚,最终合并为单个更强的λ波结构;不仅如此,两个干扰区内还形成了复杂的旋涡结构,包括一级主旋涡和二级主旋涡,这些旋涡向下游运动,最终融合成一个尺度更大的锥形主旋涡。     

隔离段内超声速流动摩擦阻力分析

以超燃冲压发动机等截面隔离段内部流动阻力特性为研究背景,以数值模拟和实验测量为研究手段,研究马赫数2.5来流条件下,平板流动和带激波反射流动条件下的壁面摩擦阻力.实验测量得到了超声速流动下的壁面摩擦阻力.研究发现,当激波反射形成局部微小分离时,激波和边界层的相互作用使得分离区下游的壁面剪切应力恢复为比分离区上游稍高水平.存在激波反射时,壁面总的摩擦阻力略大于无激波时的平板流动.     

基于DES方法的三角翼激波-涡干扰流场数值模拟

采用基于Spalart-Allmaras湍流模型的脱体涡模拟（DES）方法,数值求解Navier-Stokes方程,模拟绕尖前缘三角翼的跨音速流动,并对三角翼上翼面的复杂激波-旋涡干扰流场进行了分析。与NASA兰利研究中心的NTF风洞实验结果对比分析表明,DES方法能很好地模拟跨音速三角翼上的旋涡流动。随着攻角由中度攻角增加到大攻角,支架附近的激波越来越强,对主分离涡的干扰作用越来越大,直至出现激波干扰导致的涡破裂。激波的形状、位置及涡破裂位置均与实验结果吻合良好。     

高背压进气道中内外流耦合作用的大涡模拟

采用大涡模拟研究了出口堵塞比为50.8%的轴对称进气道流动,重点考察了内外流耦合作用下流动的非定常特性。采用国家数值风洞（NNW）工程仿真软件进行数值模拟,得到的壁面平均压力、瞬时压力分布与试验数据符合良好。分析表明：为匹配出口背压,进气道在喉道区域形成激波串结构,使内流道流场分为上游超声速区、中部激波串区以及下游亚声速区;在激波串区,剧烈的逆压梯度产生了分离激波、激波串、分离区及分离剪切层等复杂结构;伴随着激波串运动和边界层大尺度分离,进气道壁面压力出现宽频脉动特征。脉动压力的时空分布表明：内流道脉动压力以扰动波的形式传播,为此建立的声反馈模型能较好地预测亚声速区的主导频率。相关性分析表明：激波串运动受上下游流动耦合作用,其中,频率为St=0.7的运动主要受上游流动影响,频率为St=0.9的运动主要受下游压力扰动波影响。     

真实气体效应对升力体舵面局部流动分离的影响

数值模拟分析了高马赫数低雷诺数条件下激波边界层干扰、激波与激波相互作用、流动分离再附等流动现象的特点以及高温真实气体效应的影响。分别采用量热完全气体、平衡气体、化学非平衡气体模型对升力体由于舵面偏转引起的局部流动分离情形进行了数值模拟。研究了飞行高度、壁面温度及来流马赫数对流动分离的影响。计算结果表明:真实气体效应使空气在边界层内发生离解反应,边界层内温度降低,粘性减小,动能损失减小,克服逆压梯度的能力更强,从而使分离区明显减小。分离区的减小改变了分离/再附激波的位置和强度,进而对局部压力及热流分布产生重要影响;随高度增加,平衡气体较完全气体分离区相对减小量增大,平衡气体效应对流动分离/再附现象的影响越大;壁温对分离区影响较大,随壁温升高,分离区增大;随马赫数增大,分离区减小,真实气体和完全气体的差异增大,真实气体效应的影响更加显著。     

激波-边界层-分离流相互干扰三维湍流的数值模拟

本文采用数值方法求解时间相关三维可压缩雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组,模拟激波—边界层—分离流相互干扰三维湍流流动。湍流模型为Ｂａｄｗｉｎ－Ｌｏｍａｘ两层代数模型,改进后用于三维内流问题。采用单元中心有限体积法离散流场控制方程,ＶａｎＬｅｅｒ矢通量格式计算无粘通量,中心差分法计算粘性通量,ＬＵＳＧＳ时间推进格式计算定常流场。本文以二元跨音速扩压器内三流动为算例,数值模拟较强激波—边界层—分离流相互干扰维湍流流动,并与实验结果进行了比较。数值模拟结果,在激波强度、分离点位置和再附点位置等方面,与实验结果吻合较好。     

TSTO并联分离激波/边界层干扰流动特性分析

针对两级入轨飞行器的缩比模型，通过试验与数值模拟相结合的方式，在马赫数6条件下开展典型级间距状态的激波/边界层干扰流场研究，详细分析干扰区壁面及空间的流动结构与特性。结果表明：试验中模型壁面边界层在激波入射之前为层流状态，在强激波干扰后迅速转捩为湍流状态，因此试验结果在第一道激波作用结束之前与层流计算结果吻合，而在第一道激波作用结束之后与湍流计算结果一致；激波/边界层干扰呈现复杂的三维流动特征和明显的开放结构，强激波在壁面形成的高压区呈弧状向下游展开，轨道级头部产生的入射激波在级间来回反射，强度依次递减；同时，干扰区内存在展向弯曲的主分离线与再附线、沿流向排列的二次分离线与再附线、流动剪切形成的旋涡结构以及包括鞍点、结点、焦点在内的临界点；层流边界层受到激波作用形成的分离区明显大于湍流边界层，同时开放特征更为显著。     

Plasma waves at Venus

Many significant wave phenomena have been discovered at Venus with the plasma wave instrument flow on the Pioneer Venus Orbiter. It has been shown that whistler-mode waves in the magnetosheath of the planet may be an important source of energy for the topside ionosphere. Plasma waves are also associated with thickening of the ionopause current layer. Current-generated waves in plasma clouds may provide anomalous resistance resulting in electron acceleration, possibly producing aurora. Ion-acoustic waves are observed in the bow shock, and appear to be a feature of the magnetotail boundary. Lastly plasma waves have been cited as evidence for lightning on Venus.     

扰动波对轴流压气机旋转失速的影响

提出了考虑轴流压气机进出口气流角和非定常损失变化的二维旋转失速稳定性模型,进行了数值分析和讨论。计算结果表明压气机前方的扰动波对压气机旋转失速稳定性有一定的影响,压气机本身的结构形式以及压气机上、下游结构对扰动波的传播速率有重要影响     

燃烧室宽度对煤油旋转爆震波传播模态的影响—“爆震专刊”

为研究基于煤油的旋转爆震波的传播特性,以煤油和含氧量40%的富氧空气作为燃料和氧化剂,基于燃烧室外径均为100 mm的无内柱燃烧室和燃烧室宽度分别为32 mm、26 mm和20 mm的环形燃烧室开展了对比实验。不同氧化剂流量下,共观察到四种燃烧波模态,分别为爆燃模态、准稳定爆震模态、双波对撞模态和稳定旋转爆震模态。无内柱燃烧室中,氧化剂流量较低时无法维持旋转爆震波的稳定传播,出现爆燃模态和准稳定爆震模态;当氧化剂流量超过120 g/s时,可以得到旋转爆震模态,旋转爆震波峰值压力超过0.7 MPa,平均传播速度为1750 m/s。对于环形燃烧室,旋转爆震波的传播速度仅为1245~1465 m/s,明显低于无内柱燃烧室中的传播速度。随环形燃烧室宽度减小,对应旋转爆震波模态的工况范围更窄,传播速度更慢。在本研究对应的工况范围内,增大燃烧室宽度,更有利于基于煤油的旋转爆震波的稳定传播。     

二维超声速剪切层的非线性失稳过程分析

利用直接数值模拟方法，模拟了二维超声速剪切流动中的小扰动波随时间发展在空间上的传播过程，由此分析其稳定性，通过相图分析，李雅普诺夫指数分析和分维分析，本文得到了二维超声速自由剪切层流动由线性到非线性失稳过程的数学特征，结果分析表明，在本文的计算条件下，二维超声速混合层流动在失稳过程中出现了间歇现象，并进入了混沌运动的初期阶段。     

A self-similar solution of a curved shock wave and its time-dependent force variation for a starting flat plate airfoil in supersonic flow

The problem of aeroelasticity and maneuvering of command surface and gust wing interaction involves a starting flow period which can be seen as the flow of an airfoil attaining suddenly an angle of attack. In the linear or nonlinear case, compressive Mach or shock waves are generated on the windward side and expansive Mach or rarefaction waves are generated on the leeward side. On each side, these waves are composed of an oblique steady state wave, a vertically-moving one-dimensional unsteady wave, and a secondary wave resulting from the interaction between the steady and unsteady ones. An analytical solution in the secondary wave has been obtained by Heaslet and Lomax in the linear case, and this linear solution has been borrowed to give an approximate solution by Bai and Wu for the nonlinear case. The structure of the secondary shock wave and the appearance of various force stages are two issues not yet considered in previous studies and has been studied in the present paper. A self-similar solution is obtained for the secondary shock wave, and the reason to have an initial force plateau as observed numerically is identified. Moreover, six theoretical characteristic time scales for pressure load variation are determined which explain the slope changes of the time-dependent force curve.     

小区域交通流协调控制技术及应用

提出了一种在线式小区域交通信号协调控制策略,其可行性已在实际应用得到证实。通过设置在交叉口处的车辆检测器不断采集车辆信息,实时调整关联路口间的相位差,并通过合理的交通配时,实现主干道双向绿波控制,从而有效地提升道路的通行能力。     

Subsonic impulsively starting flow at a high angle of attack with shock wave and vortex interaction

Impulsively starting flow, by a sudden attainment of a large angle of attack, has been well studied for incompressible and supersonic flows, but less studied for subsonic flow. Recently, a preliminary numerical study for subsonic starting flow at a high angle of attack displays an advance of stall around a Mach number of 0.5, when compared to other Mach numbers. To see what happens in this special case, we conduct here in this paper a further study for this case, to display and analyze the full flow structures. We find that for a Mach number around 0.5, a local supersonic flow region repeatedly splits and merges, and a pair of left-going and right-going unsteady shock waves are embedded inside the leading edge vortex once it is sufficiently grown up and detached from the leading edge. The flow evolution during the formation of shock waves is displayed in detail. The reason for the formation of these shock waves is explained here using the Laval nozzle flow theory. The existence of this shock pair inside the vortex, for a Mach number only close to 0.5, may help the growing of the trailing edge vortex responsible for the advance of stall observed previously.     

Transitional wave configurations between Type III and Type IV oblique-shock/bow-shock interactions

The interactions of oblique/bow shock waves are the key flow phenomena restricting the design and aerothermodynamic performance of high-speed vehicles. Type III and Type IV Shock/Shock Interactions(SSIs) have been extensively investigated, as such interactions can induce abnormal aerodynamic heating problems in hypersonic flows of vehicles. The transition process between these two distinct types of shock/shock interactions remains unclear. In the present study, a subclass of shock/shock interaction configuration is revealed and defined as Type IIIa. Type IIIa interaction can induce much more severe aerodynamic heating than a Type IV interaction which was ever reported to be the most serious in literature. The intense aerodynamic heating observed in this configuration highlights a new design point for the thermal protection system of hypersonic vehicles. A secondary Mach interaction between shock waves in the supersonic flow path of a Type III configuration is demonstrated to be the primary mechanism for such a subclass of shock/shock interaction configuration.     

脉冲爆震环境下涡轮性能及气热负荷的时序演变规律

针对脉冲爆震发动机涡轮部件对剧烈时变来流条件和爆震波的非定常响应问题，对脉冲爆震发动机典型工况下某单级 涡轮开展了3维非定常数值仿真，详细讨论了脉冲爆震环境对涡轮流通能力、作功能力、流动损失、温度分布及受力等关键气热特 征的影响。结果表明：爆震波在导叶上游的传播和反射会显著影响涡轮的瞬态流通特征，导致涡轮进口流动在正向和逆向间反复 变化；爆震波压缩作功使进口温度大幅升高，而导叶反射波则会使流体温度进一步升高，甚至超过来流温度的峰值；转静子叶片轴 向间隙内呈现复杂的爆震波干涉与反射结构，在此影响下所研究的涡轮转子来流攻角变化范围超过100°，从而引起涡轮流通能 力、流动结构及损失的剧烈时序变化；在爆震波的冲击下，涡轮导叶排瞬态轴向力超过涡轮稳态设计点的120倍，周向负荷超过稳 态设计点的40倍，涡轮动叶轴向力和周向负荷则可达到稳态设计点的6~7倍，给涡轮结构强度造成极大的影响；爆震波和反射波 对工质的压缩作功可使导叶表面流体的最高瞬态温度达到导叶表面流体周期平均温度的3.5倍以上，动叶的也可达到2.8倍以上， 使叶片冷却的难度增大，且可能引发严重的烧蚀问题。