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

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

激波/平板湍流边界层干扰的数值模拟

采用理性GAO-YONG可压缩湍流模型,模拟了激波/平板湍流边界层干扰现象,结果表明:计算所得到的壁面压力分布、摩阻系数分布和速度型分布均与实验值吻合很好,并且比较准确地预报出了入射斜激波与平板湍流边界层干扰所引起的边界层分离点和再附点等流动特性.由此表明GAO-YONG可压缩湍流模型能够准确地用来模拟激波/平板湍流边界层的干扰流动.      

超燃进气道激波/湍流边界层干扰

针对超燃进气道湍流边界层/激波干扰引起的分离问题,采用基于5阶WENO数值格式的大涡模拟(LES)方法开展流场湍流非定常预测,旨在分析进气道湍流化技术实现进气道起动的可行性。研究表明,平板激波/湍流边界层干扰(STBLI)问题,LES方法能够清晰、可靠预测反射、分离激波形成过程及激波与充分发展湍流边界层的相互干扰,定量结果与试验一致;进气道研究方面,层流状态下,激波干扰产生强分离,导致进气道堵塞,而采用湍流化控制后试验和计算均表明流场分离明显减小,流场稳定且无明显堵塞现象,进气道可以起动,总压恢复系数达到要求,该结果表明,利用强湍流化减弱分离,实现进气道起动思想是可行的。      

Sajben跨声速扩压器分离流动中湍流模式数值研究

为研究湍流模式对激波/湍流边界层干扰内流流动的影响,提高数值计算准确度,使用SA,SST k-ω,非线性EASM k-ω,Gao-Yong四个湍流模式对Sajben扩压器内激波/湍流边界层干扰流动进行了数值计算。对流项采用Roe格式离散,扩散项采用二阶中心格式离散,离散后的控制方程用多步Runge-Kutta显示时间推进法求解。文中展示了四个湍流模式计算得到的壁面压力、速度剖面、摩阻系数等分布。计算值与实验值符合很好,四个湍流模式总体上能够较好地模拟扩压器内激波/湍流边界层干扰复杂分离流动。Gao-Yong湍流模式对分离区内的压力、速度型的模拟更加准确,而非线性EASM k-ω模式对分离再附点位置计算最理想。     

运用GAO-YONG湍流模型对扩压器内跨声速流动的数值模拟

运用GAO-YONG可压缩湍流方程组,采用同位网格SIMPLE算法,对扩压器跨声速流动中的二维激波/湍流边界层干扰现象进行了数值模拟。将计算得到的流场的时均参数与实验值进行比较,数值模拟结果在激波强度、壁面压力分布以及分离点和再附点位置等方面,与实验值吻合较好,表明GAO-YONG可压缩湍流方程组能够比较准确的模拟较强激波/湍流边界层干扰流动,从而进一步为GAO-YONG湍流模型的正确性及其在可压缩流场模拟方面的适用性提供了佐证。     

管内正激波／湍流边界干扰的数值模拟

本文采用时均Ｎ－Ｓ方程和Ｂａｌｄｗｉｎ／Ｌｏｍａｘ代数湍流模型计算了典型拉伐尔喷管内正激波与湍流边界层干扰流场。     

高马赫数下激波湍流边界层干扰数值研究

应用GAO—YONG可压缩湍流方程组数值模拟了入射斜激波／平板湍流边界层相互干扰现象,计算了来流马赫数为5．0,激波入射角度分别为15．876°、23．287°两种不同激波干扰强度下的流场。计算程序中的对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge—Kutta显式时间推进法求解空间离散后的控制方程。计算较好地模拟了高马赫数下的激波／湍流边界层干扰的流场结构,位移边界层厚度,动量损失厚度等,也比较准确地预测了平板壁面压力、摩阻系数等气动力参数的分布。     

弯曲后掠压缩拐角激波/湍流边界层干扰特性研究

为研究内转式进气道前体压缩激波与机体边界层之间的弯曲后掠压缩拐角激波/湍流边界层干扰现象,对矩形捕获型线、直母线圆锥基准流场生成的内转式进气道压缩型面进行简化,并利用数值仿真方法对简化模型进行计算,分析并对比了非耦合和耦合情况下弯曲后掠压缩拐角激波/湍流边界层干扰特性.结果表明:非耦合模型所形成的分离区呈弯刀形,分离区...     

激波/湍流边界层干扰物面剪切应力统计特性

为了揭示激波/湍流边界层干扰区内物面剪切应力统计特性的演化规律，采用直接数值模拟方法对来流马赫数2.9、12°激波角的入射激波与平板湍流边界层相互作用问题进行了研究。通过与风洞试验数据的比较分析，验证了计算结果的可靠性。系统地探究了干扰区内物面剪切应力的典型特征，如预乘谱、概率密度分布和相干结构等。研究结果表明，分离激波的低频振荡运动对流向及展向分量的预乘谱均没有实质影响，其脉动仍以高频特征为主，低频能量变化较小。干扰区内流向剪切应力概率密度函数分布变化剧烈，分离泡内对数正态分布规律不再适用，而展向剪切应力在干扰区内与正态分布较为接近。相较于上游湍流边界层，分离泡内物面剪切应力矢量夹角与幅值的联合概率密度变化显著，峰值概率降低，峰值范围增大。此外，流向剪切应力脉动场的本征正交分解分析指出，主能量模态与分离激波的低频振荡以及下游再附边界层内的Görtler-like流向涡结构密切相关。     

高超声速边界层感受性问题高精度数值模拟

通过五阶WENO格式和六阶对称紧致格式以及三阶TVD R-K结合的方法,对存在强激波和小扰动相互干扰的高超声速边界层感受性问题进行研究.结果表明:此方法能够模拟边界层内不稳定波的产生和发展以及小扰动和激波相互干扰等现象,因此能广泛应用于含激波的感受性等问题的可压缩湍流直接数值模拟以及具有激波和边界层干扰等复杂流场的计算.     

基于大涡模拟的压缩拐角激波/边界层干扰研究

为了提高对激波/边界层相互干扰的基本理解,采用大涡模拟(LES)对来流马赫数Ma=2.9,转角为24°的压缩拐角激波与湍流边界层相互干扰进行了研究。采用回收/调节方法作为入口湍流生成技术,并在超声速平板湍流边界层上进行了验证。采用涡识别方法和数值纹影图等流动显示方法,研究了干扰区内激波与边界层相互干扰的结构变化特征。通过对比分析湍动能和雷诺正应力在不同流向位置的分布规律,研究表明:经过激波干扰后湍动能主要集中在边界层的外层,并在拐角附近形成一个低湍动能区;雷诺正应力流向分量和法向分量在边界层内的分布呈现为单峰结构,而展向分量呈现为双峰结构。运用间歇因子对分离激波的大尺度流向运动进行研究,发现激波围绕着平均分离点作前后运动,运动的尺度等于进口湍流边界层厚度的72%。证实了拐角下游G?rtler流向涡对的存在,并对其展向分布和空间演化特性进行了详细研究。     

求解可压缩流动的同位网格SIMPLE方法研究

在Rhie-Chou动量插值的基础上,推导了同位网格可压缩SIMPLE算法.经过无粘流超音速凸包算例和激波/湍流边界层干扰算例计算发现,如果对流项采用高阶有界HLPA格式,密度插值采用一阶迎风和中心差分的混合格式,这种算法能够很好地模拟凸包超音流的流动现象,在采用了新型GAO-YONG湍流模型后也能够较好地模拟激波/湍流边界层干扰.      

跨声速翼型绕流计算边界层反方法的研究

本文用有限差分方法,通过有粘/无粘迭代计算了二元翼型的跨声速绕流问题。在边界层粘性区域内考虑了层流、转捩及湍流流动,当边界层内出现分离时,使用边界层反方法,采用代数湍流模型。算例表明,对激波/边界层弱干扰和强干扰情况,该方法的结果与风洞实验结果吻合良好,对于求解边界层小分离流场是一种好的近似方法。     

压缩拐角激波与旁路转捩边界层干扰数值研究

为了研究激波与旁路转捩边界层的干扰机理,采用直接数值模拟(DNS)方法对来流马赫数Ma∞=2.9,24°压缩拐角内激波与转捩边界层的相互作用进行了系统的研究。考察了旁路转捩干扰下压缩拐角内分离区形态和激波波系结构的典型特征。比较了转捩干扰与湍流干扰流动结构的差异,并分析了造成差异的原因。研究了拐角内转捩边界层的演化特性,探讨了转捩干扰下脉动峰值压力和峰值摩阻的分布规律及形成机制。研究结果表明:相较于湍流干扰,两侧发卡涡串的展向挤压使得分离区起始点以V字型分布,且分离激波沿展向以破碎状态为主,激波脚呈现多层结构;拐角内的干扰作用急剧加速了边界层的转捩过程;转捩干扰下的拐角内峰值脉动压力以单峰结构出现在分离区的下游,同时干扰区内的强湍动能和高雷诺剪切应力使得其局部峰值摩阻系数要高于湍流干扰。     

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

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

使用GAO—YONG湍流模型数值研究管道凸起流动

用GAO—YONG可压缩湍流方程数值模拟了Delery管道凸起跨音流场中的激波湍流边界层干扰现象。分析了GAO-YONG可压缩湍流方程组对湍流的非平衡、多尺度、各向异性等特性的描述能力。计算中对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge-Kutta显式时间推进法求解了空间离散后的控制方程。计算很好地模拟到了压力平台区、“入”波结构等典型激波湍流边界层干扰的流动现象,也得到了壁面压力分布、平均速度剖面以及雷诺应力分布等,并与相应的实验数据进行了对比分析,两者符合很好。     

超声速膨胀角入射激波/湍流边界层干扰直接数值模拟

为了揭示膨胀效应对激波/湍流边界层干扰区内复杂流动现象的影响规律,采用直接数值模拟方法对来流马赫数2.9、30°激波角的入射激波与10°膨胀角湍流边界层相互作用问题进行了数值研究。系统地探讨了激波入射点分别位于膨胀角上游、膨胀角角点和膨胀角下游3种工况下膨胀角干扰区内若干基本流动现象,如分离泡、物面压力脉动及激波非定常运动、湍流边界层统计特性和相干结构动力学过程等。结果表明,激波入射点流向位置改变对分离区流向和法向尺度的影响显著,尤其是当激波入射点位于角点及其下游区域。研究发现,膨胀角干扰区内物面压力脉动强度急剧减小,分离区内压力波向下游传播速度将降低而在膨胀区内将升高,膨胀效应极大地抑制了分离激波的低频振荡运动。相较于入射激波与平板湍流边界层干扰,入射激波流向位置改变对膨胀角再附区速度剖面对数区及尾迹区影响显著,将导致其内层结构参数升高而外层降低,近壁区内将呈现远离一组元湍流状态的趋势。此外,流向速度脉动场本征正交分解分析指出,主模态空间结构集中在分离激波及剪切层根部附近而高阶模态以边界层内小尺度正负交替脉动结构为主。低阶重构流场结果表明,前者对应为分离泡低频膨胀/收缩过程而后者表征为分离泡高频脉动。     

应用GAO-YONG湍流模式数值模拟三维激波/湍流边界层干扰

应用GAO-YONG可压缩湍流模式数值模拟了三维激波/湍流边界层干扰算例之一——单鳍流动.攻角20°,来流马赫数2.93,雷诺数9.8×105.对流项和扩散项分别采用Roe格式和二阶中心差分格式计算.Runge-Kutta显示时间推进方法求解了半离散的控制方程.包括壁面压力分布,边界层内流动偏移角等在内的计算值与试验数据进行了比较.准确地预测出了三维激波/湍流边界层干扰流场的主要流动特性——λ波结构,主分离涡核,膨胀区,滑移线等.计算与Alvi等提出的单鳍流动的理论模型符合很好,得到了平板表面压力以及分离线、再附线等在单鳍流动中所独有的半圆锥特性.      

边界层转捩与压缩拐角分离流动的非定常作用

为了认识激波与转捩边界层之间的相互作用,选取压缩拐角模型为研究对象,在高超声速风洞中开展了激波/转捩边界层干扰试验研究。试验在FD-20炮风洞中进行,试验马赫数为8,雷诺数0.5×107~2×107/m。试验采用了薄膜电阻温度计和常规压力传感器,分别测量压缩拐角周围的热流和压力分布。根据干扰区上下游的边界层流态,将试验分为两部分:层流/湍流干扰和转捩/湍流干扰。对比分析了边界层转捩发生在干扰区时的热流分布、压力分布以及脉动热流的非定常特性。研究结果表明,激波/转捩边界层干扰的热流和压力分布特征,不同于常规的层流干扰和湍流干扰,其介于二者之间。层流/湍流干扰的热流和压力分布特征类似于层流干扰;转捩/湍流干扰的热流和压力分布特征类似于湍流干扰。互双谱分析结果表明,当边界层转捩发生在分离区时,转捩脉动与分离泡脉动同时出现并增长。当两者幅值足够大时,转捩脉动会与分离泡脉动发生非线性耦合作用。这种耦合作用会诱导出新的频率特征的脉动结构,从而使分离区内的脉动量显著增加。     

压缩拐角激波/边界层干扰的可压缩湍流模型研究

为研究压缩拐角激波/边界层干扰问题,抓住可压缩流动中的密度变化特性,利用构造的可压缩Von Karman尺度,基于KDO(Kinetic Dependent Only)湍流模型,发展出可压缩湍流模型CKDO(Compressible Kinetic Dependent Only)。通过对8°,16°,20°和24°压缩拐角算例的数值模拟,测试了其对可压缩、激波/边界层干扰这一湍流难题的求解能力。计算结果表明,总体上CKDO模型对壁面压力和壁面摩擦阻力系数的捕捉能力优于其它模型,并且随着压缩拐角角度的增大,其描述该流动的能力更加突出。CKDO模型在24°压缩拐角处计算的分离区大小仅比实验大10%左右,明显比其它模型结果好。这表明CKDO模型在模拟激波/边界层干扰这一类流动中有较好的适用性。