三维可压气体射流在液体中剪切稳定性分析

采用时间线性稳定性分析方法，就无限域不可压有粘液体中的可压缩，无粘气体射流流动的三维扰动稳定性分析，研究较全面分析了气体压缩性，液体粘性，表面张力效应，液气密度比以及滑移速度等重要因素对气体剪切流稳定性的影响，结果表明气体压缩性助长气液剪切流动的不稳定性；液体粘性可在小范围内抑制扰动增长；表面张力作用的提高使扰动最大增长率呈线性衰减，在主导波数区域内，液气密度比愈高扰动愈强，而在最大扰动波数限附近，规律趋于相反，气液间的滑移较显著地助长扰动的不稳定。     

污染效应对超声速燃烧室贫油熄火极限影响实验研究(流场诊断特别专栏)

针对污染效应对超声速燃烧室熄火性能影响问题,采用甲烷燃烧和电阻加热的燃烧室直连式对比实验方法,开展了污染效应对燃烧室贫油熄火极限影响研究。基于超声速燃烧室实验模型,研究了煤油喷嘴结构对贫油熄火极限的影响;采用高速摄像仪和平面激光诱导荧光技术（PLIF）,研究了加热方式对贫油熄火极限的影响,获得了光学测量结果和壁面压力分布,给出了纯净来流和污染来流条件下的燃烧室贫油熄火极限。研究表明：采用甲烷燃烧加热的燃烧室贫油熄火极限对应的当量比高于电阻加热的,分别为0.135和0.105;不同加热方式来流条件下,超声速燃烧室在燃烧不同阶段的稳定模式均为凹槽火焰稳定模式,火焰基底会呈现一定程度的脉动,但稳定在凹槽剪切层内。     

射流与超声速来流剪切层作用的数值模拟研究

本文采用高精度格式求解二维Navier-Stokes方程,分别应用五阶WENO格式、六阶中心差分格式离散对流项和黏性项,时间推进采用三阶Runge-Kutta格式,并应用MPI非阻塞式通信进行并行计算,模拟了超声速射流与同向超声速流动的相互作用,目的在于研究流长产生的两个剪切层的相互作用,发现了有一定夹角的剪切层作用产生周期性次生射流的现象。本文研究了射流/超声速流剪切层相互作用、作用产生的次生射流具有脉冲性等方面,着重对次声射流的周期性等展开分析和讨论。     

超声速平面剪切层多模态混合扰动数值分析

对Ma_c=1.2二维空间发展超声速平面自由剪切层进行扰动模态及流动结构数值分析.使用线性稳定理论和扰动直接数值模拟方法对超声速剪切层中不同流向站位的流动剖面进行扰动参数分析和局部扰动直接数值模拟,获得了剪切层不同站位处的扰动特性.通过特定单频扰动下的直接数值模拟,在剪切层的不同区域内,无相差地促发了慢模态、混合模态和快模态等多种声涡模态,得到了特定单频扰动促发的多模态无相差不稳定波.      

基于Van Leer＋AUSM混合格式求解三维N-S方程

基于多块粘性结构网格,开展了三维N-S方程数值算法的研究。控制方程的空间离散采用有限体积法,在前人工作的基础上,发展了Van Leer＋AUSM混合格式并应用到对流通量的离散中,粘性项采用中心格式离散并利用格林定理计算粘性通量中的导数项,时间推进采用五步R-K法,湍流模型为S-A一方程模型。最后,以M6机翼和某超声速弹丸的粘性流场作为数值算例,计算表明：发展的数值算法对跨声速、超声速流场均具有较高的分辨率,适用于跨、超声速流场的数值模拟。     

嵌边法出流条件在可压缩流直接数值模拟中的应用

以嵌边法为出口边界条件试算了亚、超声速平板边界层扰动波的演化、超声速流中涡的传播、亚声速平板边界层尾缘扰动的演化等,以验证嵌边法作为出流条件在可压缩流中的有效性。然后,探讨了嵌边法有关参数的选取问题。     

基于简化Boltzmann方程的超声速流稀薄效应分析

从近连续流到近自由分子流,计算分析超声速流场的稀薄效应。模型Boltzmann方程先后经简化分布函数和离散速度坐标变换后,采用一个隐式通量修正二阶迎风TVD格式差分求解。壁面取漫反射气体分子模型。在Knudsen数为0.001,0.01,0.1,1.0,10等5种情形下,数值模拟二维圆柱氩气体绕流,观察到了不同的波系及尾迹结构。总阻力系数计算值同实验数据基本吻合,碰撞项Shakov动力学模型反映的流场细节略好于修正BGK模型。稀薄效应趋于弱化激波等超声速流动结构,物体影响域增大。     

超声速空气介质自由扩散混合流涡的演化

在不同谐波扰动促发方式下,观察二维超声速空气/空气自由扩散混合流中涡的形成及演化.考虑含空气组分扩散效应的可压缩Navier-Stokes方程,对流项采用3阶迎风紧致格式离散,输运项采用6阶对称紧致格式离散,非定常时间推进采用3阶紧致存储显式Runge-Kutta方法.在无相差情形下,获得了大尺度基频涡结构的饱和、一次对并、二次对并、三涡对并等现象.在Ma_c=0.3低对流Mach数下,基频涡较饱满,但流向尺度较小.受空气真实气体特性影响,在Ma_c=0.8高对流Mach数下未发现涡/小激波结构.      

扩散抛物化(DP)NS方程组的意义及其在计算流体力学中的应用

本文首先讨论扩散抛物化（DP）NS方程组的早期研究工作：它的提出、数学性质、意义和在CFD中的应用，然后讨论扩散抛物化理论的一些新发展。这些新发展是对NS方程组数值计算进行物理分析的基础上得到的，其中包括NS方程组差分计算时，粘性剪切流对网格间距和格式精度的要求；粘性项只保留剪切项的广义扩散抛物化（GDP）NS方程组，它的性质和应用。由于高Re数流动在NS方程组的差分计算中，网格Re数彼此相差悬殊的特点，产生了计算离散单元守恒方程组的新的算法思路，即离散流体力学（DFD）算法。在DFD算法中需要同时计算三种不同的守恒方程组（Euler,DPNS和NS方程组）。本文讨论了DFD算法中需要同时计算三种不同的守恒方程组（Euler,DPNS和NS方程组）。本文讨论了DFD格式的构造、它的优点和应用。并以超声速绕前后台阶流动为算例，来说明GDPNS方程组的用处和DFD算法的优点。DPNS方程组、GDPNS方程组、DFD算法是高智提出的，对这些问题，他的合作者从理论、算法、数值验证和某些应用又作了系统的研究。     

基于SIMULINK的超声速进气道仿真方法

将进气道一维流动的偏微分方程离散为一系列可积分的线性常微分方程,基于此创建了超声速进气道仿真模型,提出了一种基于SIMULINK的进气道仿真方法,模拟了来流参数和燃烧室反压扰动时流场内各个参数的变化规律。对一个一维变几何超声速进气道进行仿真,获得了结尾正激波的稳态特性及动态特性。该模型定义并输出了激波位置,有效地将一维超声速进气道的分布参数问题转换为零维问题,解决了系统分析与控制系统设计中无法实现集总参数的困难。     

不同马赫数的无粘和粘性流动高阶精度隐式计算方法

 应用隐式时间推进法对不同马赫数的无粘和粘性流动进行数值分析,给出了基于预处理方法的高阶精度隐式求解方法。方程离散采用改进的高阶精度对流迎风分裂格式。该方法通过求解曲线坐标系可压缩Euler和 Navier-Stokes方程,对低速到超音速范围内的无粘和粘性流动的典型问题进行了数值分析。计算结果与文献的数值结果或实验数据对比分析表明,该方法对低速到超音速范围内的无粘和粘性流动问题进行数值分析是可行而有效的。     

自吸气脉冲爆震发动机超声速进气道气动性能的数值研究

采用有限体积法计算了某轴对称超声速进气道,在不同幅度和频率的出口扰动压强下的进气道内结尾正激波的运动情况,得出了进气道内结尾正激波运动特性和扰动压强的频率和幅度的关系。在计算中,本文采用了多块结构化网格,控制体积的界面无粘通量采用三阶迎风格式插值获得,同时采用了minmod通量限制器,以确保在激波处的解的物理特性;扩散通量采用二阶中心差分格式插值获得。定常计算采用当地时间步法,非定常计算采用双时间步法。离散的代数方程采用交替方向迭代法求解。     

Navier—Stokes（NS）方程组差分计算中的物理和网格尺度效应及NS方程组的简化

对于高Re数流动计算，在通常二阶精度NS差分格式和网格数条件下，存在某些粘性项落入修正微分方程截断误差项的问题。这类NS方程组计算实际是计算某种简化NS方程组，而且重复计算误差物理粘性项既浪费机时和内存，误差积累又会对数值解产生不可预测的影响，避免外述缺陷的办法一个提高NS差分格式的精度，另一个是丢掉可能落入截断误差项的物理粘性项，把NS方程组简化为广义NS方程组，广义NS计算避免了误差物理粘性项误差积累对数值解的不可知影响，又可节省内存和机时，对高Re数流体工程计算很有好处。利用广义NS方程组计算超声速绕前向和后向台阶流动的结果表明：广义NS方程组与NS方程组的数值结果很好相符符。     

强粘性层流理论及在粘性流计算中的应用

对高雷诺数流动计算，为了解决部分粘性项甚至全部惯性粘性项落入误差与流动物理尺度的关系问题，本文提出强粘性流动理论。强粘性流中至少有一个粘性项与惯性项同量阶，理论包含了物理尺度各向极限、经典边界层和多层边界层理论为其特例，给出了从经典边界层向物理尺度各向相同极限演化的尺度规律和粘性惯性诸项变化的量阶关系，阐明了粘性与惯性力强相互作用将在剪切层的法向以流向同时“激发”涉尺度结构。对粘性流计算，利用强粘性剪切流尺度律重新标度NS格式的修正微分方程，给出临界网格尺度与流动物理尺度和差分格式精度的关系，得到部分粘性项落入误差和计算结果为非物理数值粘性解的二个判据。并以流场中的边界层、驻点和分离点领域计算为例说明理论的应用，对强粘性剪切流计算、证实部分粘性项甚至全部粘性惯性项落入误差的问题值得重视。     

可压混合层涡并放热效应的涡动力分析

 基于双元单步不可逆燃烧动力学简化模型 ,采用紧致差分算法直接求解二维可压 Navier-Stokes方程 ,研究了时间发展非预混扩散火焰平面反应混合层流动。拟序结构由数值线性稳定性分析方法设定的粘性可压最不稳定小扰动初场促生。随放热强度增大 ,得到了基频涡卷撕裂现象。涡动力分析表明 ,放热主要通过膨胀和斜压这两种强非线性过程抑制大尺度剪切涡的演化     

基于IPDG方法的超声速混合层流动数值模拟

为了满足超声速混合层高精度模拟需求,实现了基于内罚方法的间断伽辽金(IPDG)方法数值模拟.通过将黏性通量作为辅助变量使得Navier-Stokes方程降阶,并利用间断伽辽金方法进行空间离散,最后采用Newton-Krylov隐式方法对空间半离散方程进行时间推进.相对于有限体积法数值精度提高到了 3阶.将该方法应用于对...     

不起动流场对超声速/高超声速进气道自起动性能的影响

对7个典型速域的二维超声速/高超声速进气道加速自起动过程进行了准定常数值仿真,分析了真实临界不起动流场对进气道自起动性能的影响,研究发现:存在一种介于超声速和高超声速临界不起动模式之间的过渡临界不起动模式。当真实不起动流场处于超声速临界不起动模式时,自起动马赫数略大于无黏设计自起动马赫数;处于过渡临界不起动模式时自起动马赫数小于无黏设计自起动马赫数;而该研究中处于高超声速临界不起动模式的进气道,自起动马赫数明显大于无黏设计自起动马赫数。高超声速临界不起动模式下的喉道截面特征气流参数显著偏离无黏临界不起动流场,所以Kantrowitz理论以及基于该理论发展而来的系列方法不适用于预测高超声速进气道自起动性能。     

不同模型堵塞比的超声速风洞二次喉道优化

为了延长采用真空球排气系统的超声速风洞的工作时间,拟采用二次喉道提高真空球临界工作压强。本文针对Ma=4和Ma=7超声速风洞,采用计算流体软件系统地研究了不同发动机堵塞比对二次喉道起动能力的影响以及不同尺度的二次喉道对真空球临界工作压强的影响。研究结果表明,随着发动机堵塞比的增大,二次喉道的临界起动直径相应增大,呈非线性变化,同时,真空球的临界工作压强下降。因此,二次喉道直径应针对不同的模拟状态和发动机堵塞比范围综合优化选取。     

基于超声速边界层转捩的WENO格式对比分析

基于高精度weighted essentially non-oscillation(WENO)格式对马赫数为4.5,雷诺数为10000的超声速来流条件下平板边界层转捩过程进行了大涡模拟.无黏通量分别用5阶、7阶、9阶WENO格式进行离散,黏性通量离散采用4阶中心差分格式,时间推进采用具有总变差递减(TVD)性质的3阶精度Runge-Kutta方法.通过在入口边界叠加一对等幅最不稳定第一模态谐波扰动,分别采用3种WENO格式计算得到了平板层流边界层失稳转捩的演化过程.结果表明:5阶WENO格式的数值耗散明显高于7阶和9阶WENO格式,在捕捉湍流涡和流场脉动特性上存在明显不足.9阶WENO格式的耗散小,能够捕捉到流场更小尺度的涡和高频脉动.研究具有高频脉动特性的问题或者欲捕捉精细涡结构时,建议采用7阶以上的高精度格式.      

滑流层对轴对称外压式进气道稳定性的影响

本文论述了在来流M1.72条件下,一个可变中心锥外压式进气道,强度为5～11％的六种滑流层在各种位置上进入唇口时对进气道稳定性的影响。 试验证明,强度低于10％的滑流层在唇口附近进入该轴对称进气道,并没有引起进气外罩内侧的附面层分离 强度直至11％的滑流层可以在进口处任意位置上进入该进气道,而不产生喘振。文献表明,强度仅为6～7％的滑流层就可以导致二元进气道的喘振,从而显示了轮对称进气道对滑流层进入所引起的喘振有较强的抵抗能力。