二维喷管的初始流动

基于可压缩Navier-Stokes方程,采用大涡模拟方法与高精度混合WENO/TCD格式,对Ma=1.4的超声速平面射流初始流场进行了数值研究.数值结果清晰地描述了超声速平面射流初始流场的结构特征,包括主涡环与激波结构以及它们演变过程.因主涡环内存在涡导激波对,激波与涡相互作用加速射流剪切层失稳,使剪切层首次卷起形成小涡的位置出现在涡导激波对后,此与亚声速射流情况不同.小涡串卷起成后,相继与涡导激波对相互作用,使激波出现明显的变形并加速主涡环失稳.      

用于可压缩自由剪切流动的湍流混合长度

抓住可压缩流动变密度特性,构造出基于有效涡量的三维von Karman混合长度。湍流模型采用仅依赖湍动能k的单方程KDO(Kinetic Dependent Only)模型,引入新构造的混合长度替换旧尺度得到CKDO模型。为了验证其描述可压缩自由剪切湍流的能力,选择无壁面束缚、密度梯度大和可压缩效应强的自由剪切混合层为算例,其对流马赫数Ma_c=0.8。计算结果表明,KDO模型对混合层的速度分布有着良好的控制和模拟,而经可压缩修正后的CKDO模型与原模型及其他可压缩修正模型相比,所计算的速度分布、主雷诺剪切力和混合层厚度与试验结果更加接近,说明了该混合长度对可压缩混合层这种自由剪切湍流有着良好的刻画能力。     

后向台阶超声速流动的动力混合RANS／LES模拟

针对可压缩湍流流动,构建了1种动力混合RANS/LES模式。运用该模式模拟了超声速后向台阶流动,结果再现了二维后向台阶剪切层的拟序结构,计算得到了统计平均的压力分布,回流区长度与实验结果符合较好。与定系数的一般混合RNAS/LES模拟结果的对比表明,动力混合RANS/LES模式在对流场拟序结构的捕捉、流动时均压力分布方面表现较优。     

LU-AUSMPW混合格式在可压缩无粘和粘性流动数值分析中的应用

本文首先将 AUSMPW格式与三阶 MUSCL格式融合, 给出了其在任意曲线坐标下的三维形式, 并将其与 LU-SGS格式结合, 应用于可压缩 Euler和 Navier-Stokes方程的求解。其次, 构造了一种新的通量限制器。最后, 为了验证 LU -AUSMPW混合格式的性能, 将平面叶栅跨音速无粘流动以及喷管超音速粘性流动作为算例。本文计算结果与文献计算结果和实验数据相符很好, 表明采用 LU-AUSMPW混合格式数值模拟可压缩流场, 具有较高的计算精度、较快的收敛速度和良好的稳定性。      

喷管流场和羽流场的N－S数值模拟

利用时间推进的方法，通过求解可压缩的、雷诺平均的Ｎ－Ｓ方程，并利用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ两层代数湍流模型和Ｐｒａｎｄｔｌ混合长度模型，数值模拟了飞机和火箭发动机喷管的内流和羽流，获得了喷管内流场和羽流场的参数分布，计算结果与试验结果及国外研究成果一致。网格划分采用了自适应网格技术。     

超声速对流推力矢量喷管的数值模拟

利用二次流对流剪切激励主射流方式来产生射流对流推力矢量方法是一种效率很高的矢量化喷流新方法,又称为对流推力矢量法。本文运用三阶精度MUSCL TVD有限体积格式求解二维非定常、可压缩、雷诺平均N-S方程和Spalart-Allmaras湍流模型,计算了对流推力矢量喷管的超声速射流问题,获得几种不同二次流吸气压强下对狭缝和喷管肩部壁面压强分布、剪切层诱导旋涡的结构特征和主射流矢量角等结果,计算结果与实验数据对比取得良好的一致性,并对计算结果进行了详细分析和讨论。     

动量矩方程的剪切积分解析表示式

本文采用两层涡粘性模型和速度剖面的指数律,导出了可压缩动量矩方程的剪切积分解析表示式,从而大大简化了边界层计算。和其它理论结果相比,本结果更接近于实验数据。     

喷管流场和羽流场的N－S数值模拟

利用时间推进的方法，通过求解可压缩的、雷诺平均的Ｎ－Ｓ方程，并利用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ两层代数滴流模型和Ｐｒａｎｄｔｌ混合长度模型，数值模拟了飞机和火箭发动机喷管的内流和羽流，获得了喷管内流场和羽流场的参数分布，计算结果与试验结果及国外研究成果一致。网格划分采和了自适应网格技术。     

守恒型可压缩一维湍流方法及其在超声速标量混合层中的应用

一维湍流（ODT）方法是一种能在一维计算域上遵循湍流基本物理规律的湍流建模方法。通过结合确定性和随机性求解方法,能够在一维计算域上准确捕捉到湍流统计规律,且降维建模可显著减小计算量。ODT方法主要被广泛用于不可压湍流和湍流燃烧研究,若要将其拓展用于模拟高速可压缩湍流,需对建模方法进行深度改进。相比于不可压ODT方法,本文基于欧拉参考框架,针对可压缩湍流的特性,将因变量由原始变量改为有利于减小可压缩湍流模拟误差的守恒通量,并加入了组分求解模块。对确定性和随机性求解模块均进行了相应的深度改进,开发出具有标量混合模拟功能的守恒型可压缩ODT方法。在确定性模块中改为求解以守恒通量为变量的一维截断控制方程,在随机性模块中构造一维涡时,将三联映射的作用对象也相应地由原始变量改为守恒通量,并选用了可保证变密度情况下动量守恒的双核变换。通过模拟空间发展超声速平面湍流混合层并将自相似阶段结果与实验结果比对,验证该方法对可压缩剪切湍流场中标量混合的捕捉精度。守恒型可压缩ODT方法模拟得到的速度场和组分场的平均剖面和脉动强度分布与实验结果准确吻合,精度明显优于传统的耦合梯度扩散亚格子模型的大涡模拟方法（LES-GRAD.DIFF.）以及耦合线性涡（LEM）亚格子模型的大涡模拟方法（LES-LEM）,且该方法的降维处理使其在降低计算成本方面具有显著优势。     

低雷诺数κ-ε封闭模式在三维湍流边界层计算中的应用

本文采用低雷诺数形式的k-ε湍流模型,计算具有大横向流的不可压缩三维湍流边界层。在湍流模型中,应用了涡粘性的各向异性系数。本文所得结果与用混合长模型所得结果以及量测结果的比较表明:对于边界层的不同特征量,湍流模型和涡粘性各向异性对计算结果的影响程度是不同的。总的来说,应用考虑涡粘性各向异性的k-ε模型,较之应用混合长模型,所得计算结果有所改善,但还不是根本性的改善。     

亚跨超音速喷管内流场的整体矢通量分裂算法

从N－S方程出发，采用先进的矢通量分裂算法，对亚跨超音速喷管内流场进行了整体模拟计算。文中对几种不同计算格式的计算效率进行了比较，给出了超音速占优混合流动的直接分割求解方法。计算采用代数湍流模型，跨音速段的计算结果同其它实验和计算结果进行了比较。计算中发现对出口边界条件进行亚音速、超音速和有回流分别处理很有必要；边界层内的M数沿壁面法向变化非常快；不同壁温条件和不同的流动Re数均会改变边界层的厚度。     

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

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

预燃/流向涡掺混超声速燃烧室的稳焰火炬实验研究

 为研究飞行马赫数Ma_flight=4~7的双燃室碳氢燃料超燃冲压发动机燃烧室的原理和工程参数,进行了直连双燃室超声速冷主流和亚燃室稳焰火炬热流的掺混实验和燃烧实验。将进气道输出的超声速气流的10%流量经亚燃进气道导入亚声速预燃室,先低速地与雾化预燃油掺混并建立稳定的预燃。该预燃气流与二次喷入的主燃油掺混而形成富含吸热分解油气的高温射流,再经一组波瓣掺混器与超声速主流在下游流向涡中深入掺混/燃烧,扩大燃区厚度而趋于深入超声流层,以期实现稳定超燃。在总温约为285 K、总压为1.5×10⁶ Pa和1.0×1.0⁶ Pa,燃烧室进口马赫数Ma_inlet=2.5的来流下,对3种不同结构参数的预燃室和一种超燃室,进行了冷态流场和预燃/主燃的喷油/燃烧实验。实验与计算结果表明,冷/热态实验中整个超燃室保持了超声速流动,尽管斜激波系存在一些变化。利用存在的4种旋涡掺混现象,增强超/亚声速流之间的掺混。当采用三波系进气道和较小容积热强度的大体积预燃室和流向涡掺混器,可以形成稳定的高温富油火炬,成为超燃室稳定点火源。在超燃室下层流层的原无预热冷态来流的亚声速和低超声速区域中出现火焰,且其并不破坏超燃室上层的高超声速未燃流动。     

亚超音速机翼最佳弯扭综合设计计算方法

本文介绍一种亚超音速机翼最佳弯扭综合设计的计算方法,它应用了有限基本解方法。分别在亚超音速各选取一个设计点(M数和C_L),进行机翼弯扭设计,其目的是减小与升力相关的阻力。在此基础上,顾及亚音速和超音速这两个设计点的气动力特性,还要兼顾到飞机其它性能和结构上实现的可能性,进行机翼的综合设计。本文分别给出了亚音速最佳弯扭设计,超音速最佳弯扭设计和综合设计的计算结果。经过分析表明,计算结果是合理的。     

Shock train and pseudo-shock phenomena in internal gas flows

The interaction between a normal shock wave and a boundary layer along a wall surface in internal compressible flows causes a very complicated flow. When the shock is strong enough to separate the boundary layer, the shock is bifurcated and one or more shocks appear downstream of the bifurcated shock. A series of shocks thus formed, called “shock train”, is followed by an adverse pressure gradient region, if the duct is long enough. Thus the effect of the interaction extends over a great distance. The flow is decelerated from supersonic to subsonic through the whole interaction region. In this sense, the interaction region including the shock train in it is referred to as “pseudo-shock” in the present paper, as Crocco called it. The shock train and pseudo-shock strongly affect the performance and efficiency of various flow devices. In the present review some fundamental characteristics of the shock train and pseudo-shock are first described. Some simple predictions are made to simulate these very complicated phenomena. Pseudo-shocks appearing in various flow devices are explained. Control methods of the pseudo-shocks are also described. Finally, the current understanding of self-excited oscillation of pseudo-shock is reviewed.     

计算跨音速三维机翼非定常载荷的核函数法

 本文给出了在跨音速流动中,计算三维机翼非定常载荷分布的核函数法和计算结果。此方法是基于线化升力面理论,利用在亚、超音速情况下,联系载荷分布与下洗的积分方程的核函数,采用局部线化法来处理跨音速的非定常绕流问题。在计算中引进了线分布的冲波偶极子,并嵌入适当的正冲波边界条件。 我们从统一的超、亚音速核函数表达式出发,对亚音速区数值积分作了新的简便的处理,对于亚、超音速区相互的诱导下洗给出了自己的数值积分过程。我们算了三个例题,和文献上的数值符合较好。     

亚／超双模态超音速燃烧的实验研究

描述了在超音速燃烧室中如何从亚音速燃烧模态转变为超音速燃烧模态的概念和目前的研究状态，并用实验表明改变氢燃料喷射位置和当量比对燃烧状态的影响，以及燃烧室各截面参数的变化规律。     

煤油在超声速气流中非定常横向喷射的实验观察

本文目的是研究煤油在超声速气流中非定常横向喷射问题,并与静止和亚声速气流中的结果进行对比。改造普通油嘴并采用单次高压燃油喷射系统,得到不同延时的流场阴影照片。结果表明：在超声速气流中,射流柱破碎是由迎风面的表面波引起的,破碎点位于表面波的波谷。雾化过程过分四个区域。射流激波由近射流柱的迎风面,并在人壁反射,非定常喷射的波系演化复杂。煤油在亚声速气流中喷射,不存在射流激波,破碎点沿射流柱上移,穿透深     

C/C复合材料烧蚀形貌测量及烧蚀机理分析

对C／C复合材料试件的表面细观烧蚀进行了常压下的亚、超声速和高压下的亚音速烧蚀形貌的测量；根据测量结果，分析了C／C复合材料在上述情况下的质量损失规律。结果表明：C／C复合材料在亚音速流场的条件下，z向纤维束首先发生剥蚀，当压力升高时，碳布会发生层间剥蚀的现象；而在超音速条件下，碳布更容易发生剥蚀的现象。     

涡轮叶栅超声速流场流动特征与气膜冷却特性

应用shear strain transport(SST) k-ω 两方程湍流模型,对超声速涡轮叶栅通道内气膜冷却特性进行数值研究,得到不同气膜孔倾角和吹风比下叶栅通道内流场流动特征以及气膜冷却效率的变化规律.在激波入射点附近的气膜射流能够向分离区边界层中补充动量,克服逆压力梯度,有效改善由于激波引起的局部过热.亚声速流动状态下的气膜入射角度对冷却效率的影响能够在较大吹风比下得以体现,而超声速主流状态下,气膜冷却效率与入射角度基本无关,说明亚声速的气膜冷却射流对超声速主流的穿透力要弱于对亚声速主流的穿透力;超声速主流条件下,在激波入射位置的气膜冷却效率要高于激波入射位置下游的气膜冷却效率,这与气膜孔出流在当地的湍流度有关.