钝体驻点线附近粘性激波层的数值方法

一、引言 在研究高超音速钝头体粘性气体绕流时,驻点区域的气动力参数和热力学参数的预报也是十分重要的。当粘性层很薄时,一直采用古典边界层方程计算;但是,考虑粘性气体与无粘流的相互干扰时,边界层模型已不适用,可用粘性激波层方程来描述极其广泛的来流马赫数及雷诺数的流场。 近年来,不少人致力于驻点线附近流场的计算,有的跨激波,有的分离激波。但     

用N-S方程求解柔壁翼型绕流

用 N- S方程数值模拟带有自适应柔壁被动控制的翼型粘性绕流。用时间推进法求解这种准定常绕流。在时间历程中 ,柔壁自适应变化 ,最后与整个流场一起趋于定常。比较了柔壁和刚壁两种翼型的压力场和等马赫线及等密度线等     

直升机旋翼悬停流场的粘性数值模拟

在旋转坐标系下求解可压缩NS方程,对Caradona-Tung旋翼的跨声速悬停流场进行了数值模拟,计算中采用嵌套网格方法,并且应用了多重网格技术加速流场收敛。计算结果表明,多重网格技术显著改善了流场的收敛效率,计算量仅为单重网格的四分之一。文中同时还对无粘和粘性流场进行了对比,NS方程考虑了物面的粘性作用,能正确模拟出激波和边界层的干扰现象,计算得到的桨叶表面压力分布与实验值更加吻合。     

跨音叶栅非定流场数值模拟

本文将二维平面振荡叶栅的非定常流场分解为定常平均流场和以简谐振荡为运动规律的非定常小扰动的叠加 ,分别求解定常的 N- S方程和非定常线化 N- S方程 ,能对不同分离形式进行数值模拟 ,有效地处理分离流动。计算结果表明 ,非定常分离、粘性效应和跨音激波对叶栅气动弹性稳定性有极其重要的影响 ,并且说明分离区内的脉动压力对叶栅振荡的响应只比分离区外低一个数量级 ,因而在处理非定常分离流动时不能简单地将分离区内的压力脉动视为零。     

超声速主流中横向喷流干扰流场数值模拟

本文用二阶迎风 TVD格式求解层流 N- S方程 ,数值模拟了超声速主流与横向喷流相互作用复杂干扰流场 ,得到了合理的涡系结构 ,计算波系结构及沿物面压力分布也与实验结果一致。计算还表明 ,TVD格式不但具有较强的捕捉激波的能力 ,而且用来模拟复杂涡系结构也是有效的     

细长体大攻角绕流的数值模拟

用数值模拟方法研究了细长体大攻角的绕流特性（求解在周向３６０°范围内进行）。数值模拟的出发方程是三维可压流的全Ｎ－Ｓ方程,格式为二阶迎风ＴＶＤ格式。通过数值模拟研究了附近,攻角１９°、３２．５°和５０°状态下,细长体绕流场中的激波、分离流动及表面压力分布等重要特征。     

跨音速飞机模型试验洞壁干扰数值模拟的初步研究

分别以固壁条件和洞壁附近实测的压力分布,模拟各类实壁和透气壁试验段的洞壁边界条件,利用N-S方程数值求解模型在风洞中的绕流场,得出洞壁干扰对跨音速模型绕流的影响和气动力修正结果。初步研究结果表明,该方法能较有效地模拟模型在跨音速风洞中的绕流场,GBM04A模型在0.6m风洞中的试验结果经洞壁干扰修正后与无干扰参考结果吻合较好。     

三角翼过失速非定常洞壁干扰修正

分别以固壁条件和洞壁附近实测的压力分布,模拟各类实壁和透气壁试验段的洞壁边界条件,利用N-S方程数值求解模型在风洞中的绕流场,得出洞壁干扰对跨音速模型绕流的影响和气动力修正结果。初步研究结果表明,该方法能较有效地模拟模型在跨音速风洞中的绕流场,GBM04A模型在0.6m风洞中的试验结果经洞壁干扰修正后与无干扰参考结果吻合较好     

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

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

跨声速三维非线性洞壁干扰的数值计算

分别以固壁条件和洞壁附近的压力分布模拟各类实壁和透气壁试验段的洞壁边界条件，利用Ｅｕｌｅｒ方程和Ｎ－Ｓ方程数值求解模型在风洞中的绕流场，得出洞壁干扰对跨声速模型绕流和气动力的影响，初步的研究结果表明，该方法能较有效地模拟模型在跨声速风洞中的绕流场，经洞壁干扰修正后的ＧＢＭ－０４Ａ模型在０．６ｍ风洞中的试验结果与无干扰参考结果吻合较好     

基于LES方法的平板非定常激波/湍流边界层干扰研究

以高超声速发动机进气道湍流分离控制为应用背景,采用大涡模拟(LES)方法进行马赫数为3.0(唇口附近马赫数约为3.0)的激波/湍流边界层干扰(SWTBLI)流场机理研究.利用扰动循环引入的方法,先得到充分发展湍流场,然后根据斜激波关系式引入激波的方法进行激波/湍流干扰模拟.研究结果显示:充分发展湍流场在激波作用下产生逆...     

可压缩粘性流绕翼型的N─S方程自适应网格解

从Ｎ－Ｓ方程出发,采用ＬＵ－ＡＤＩ隐式分解方法求解绕翼型的可压缩粘性流动。湍流模式采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模式和ｑ－ω二方程微分模式以研究湍流边界层的非平衡效应;为更好地捕捉激波,对网格进行了自适应调整。通过对ＮＡＣＡ００１２和ＲＡＥ２８２２两种翼型在亚音速和跨音速不同马赫数和雷诺数下进行大量计算,表明该方法对改进湍流计算、提高激波分辨率有较好的效果     

用CSCM格式模拟二维收缩管道内的流动

 采用Lombard等人提出的CSCM(Conservative Supra-Characteristics Method)差分方法数值模拟二维收缩管道内波系干扰及激波边界层干扰等问题。所得结果同已往计算结果进行了比较。通过压力场和速度场的计算结果可看到Ma数对流场的影响。在Ma数较小时激波角较大,可看到激波经两次相交一次反射逐渐变弱的情况。Re数改变对流场也有明显影响,且给出由于激波干扰引起边界层分离的计算结果。最后还给出收缩管道内激波由正激波逐步形成斜激波的非定常过程。     

跨音速粘流的计算

 概述了一种采用粘流/无粘流相互作用原理计算跨音速粘流的方法。采用熵修正激波算子计及跨越激波时的熵增,形成的非等熵位势方法可比传统的位势方法更准确地计算无粘流动。提出了一种流向速度型,结合其它辅助关系式导出了三维湍流边界层积分方程反方法。用此方法可求得粘流解。利用半反耦合方式耦合了无粘流和粘流解。数值算例表明,计算结果与实验结果吻合;且对计算机的要求较低。     

高超声速黏性干扰相关性参数的修正

为了更精确地将黏性干扰相关性理论应用到实际工程上,基于高超声速平板黏性干扰相关性理论和斜激波理论,针对当前广泛采用的黏性干扰相关性参数(∨)∞进行了修正,建立了气动力关联性更好的参数(∨)∞.采用计算流体力学(CFD)数值方法模拟了航天飞机轨道器外形OV-102在典型状态下的黏性干扰效应,考察了修正后黏性干扰相关性参数(∨)∞对于气动力数据的关联特性.验证结果表明:修正后的黏性干扰相关性参数(∨)∞对于黏性诱导的各气动力数据均有良好的关联特性,并且相比未修正的黏性干扰相关性参数(∨)∞更加有效.     

一种修正的Osher通量在高阶WCNS中的特性

将一种基于多维修正的Osher通量运用于高阶加权紧致非线性格式(WCNS)中,该修正通量主要在垂直于激波面的界面上增加耗散,能够改善Osher通量的激波捕捉稳定性,同时对边界层和接触间断的分辨率影响非常小。对修正的Osher通量在高阶WCNS中的特性进行研究,通过数值模拟测试了基于Osher通量的WCNS的激波稳定性、热流预测精度、边界层模拟能力、激波边界层干扰模拟能力,并与Steger-Warming通量和Roe通量进行了对比。结果表明修正后的Osher通量比Harten修正的Roe通量具有更好的激波捕捉鲁棒性,而边界层、驻点热流值和激波边界层干扰的模拟则明显优于Steger-Warming通量。上述结果说明了基于修正的Osher通量的高阶WCNS具有较好的激波捕捉特性、热流预测精度和边界层计算能力。     

Viscosity effects on weak irregular reflection of shock waves in steady flow

The viscosity effects on weak shock wave reflection are investigated with the Navier–Stokes and DSMC flow solvers. It is shown that the viscosity plays a crucial role in the vicinity of three-shock intersection. Instead of a singular triple point, in viscous flow there is a smooth three shock transition zone, where one-dimensional shock jump relations cannot be applied. At the flow parameters corresponding to the von Neumann reflection, when no inviscid three-shock solution exists, the computations predict an irregular shock-wave configuration similar to that observed previously in experiments. The existence of a viscous zone in the region of shock-wave interaction allows a continuous transition from the parameters behind the Mach stem to the parameters behind the reflected shock, which is impossible in the inviscid three-shock theory. Results of numerical simulations agree qualitatively with the conclusions of the approximate viscous theory, which describes the flow in the vicinity of the three-shock intersection.     

考虑动叶顶部间隙影响的跨音涡轮级的数值模拟

本文给出了一种基于区域分解算法的可考虑动叶顶部间隙和机匣刮削作用的涡轮级计算方法。并对NASA实验涡轮级进行了计算, 计算结果和实验结果吻合得很好。通过分析, 可以发现动叶顶部间隙流动将引起动叶顶部气流的欠偏转, 从而降低顶部的做功能力, 并且这一影响区域会达到一个相当大的范围。在动叶叶片顶部尾缘存在激波的情况下, 尾缘激波和顶部泄漏流动相互作用使激波结构变得十分复杂和无序, 这种无序的结构必然带来较大的损失。      

A Structured Mesh Euler and Interactive Boundary Layer Method for Wing/Body Configurations

To compute transonic flows over a complex 3D aircraft configuration, a viscous/inviscid interaction method is developed by coupling an integral boundary-layer solver with an Eluer solver in a ＂semi-inverse＂ manner. For the turbulent boundary-layer, an integral method using Green＇s lag equation is coupled with the outer inviscid flow. A blowing velocity approach is used to simulate the displacement effects of the boundary layer. To predict the aerodynamic drag, it is developed a numerical technique called far-field method that is based on the momentum theorem, in which the total drag is divided into three component drags, i.e. viscous, induced and wave-formed. Consequently, it can provide more physical insight into the drag sources than the often-used surface integral technique. The drag decomposition can be achieved with help of the second law of thermodynamics, which implies that entropy increases and total pressure decreases only across shock wave along a streamline of an inviscid non-isentropic flow. This method has been applied to the DLR-F4 wing/body configuration showing results in good agreement with the wind tunnel data.     

跨声速粘性流中的翼型设计与分析

本文采用一种数值方法设计和分析跨声速翼型。几个算例表明,对于激波/边界层弱干扰和强干扰情形,该方法的结果与风洞实验值吻合良好。但是,当翼型表面出现分离时,两者的差别较大。