隐式紧耦合SST和TNT湍流模型的高速流动数值模拟

将SST(shear stress transport)和TNT(turbulent/non-turbulent)湍流模型输运方程与平均流场控制方程进行隐式紧耦合求解,结合当地时间步长方法和湍流源项隐式处理确保求解过程的快速和稳定.采用AUSMPW+(AUSM by pressure-based weight functions)格式和LU-SGS(lower-upper symmetric Gauss-Seidel)隐式紧耦合方法对高超声速压缩拐角流动、锥柱裙流动和超声速非对称激波/边界层干扰问题进行了数值模拟.计算结果与实验值的对比表明:SST模型和TNT湍流模型可以很好地预测15°压缩拐角流动的壁面压力和热流密度;随着压缩拐角的增大,计算结果与实验值偏差增大;可压缩性修正对压缩拐角流动的压力和热流密度分布有很大影响,对超声速非对称激波/边界层干扰流动影响很小;隐式紧耦合方法比显式紧耦合方法具有更好的收敛特性.      

传感器局部温度差异对压缩拐角热流测量的影响

针对高超声速飞行器在进行热流测量时存在的传感器与其周围热防护材料之间的热匹配性问题,即冷点效应问题,采用自主研发的计算程序研究了冷点效应对压缩拐角表面热流分布的影响,得到了在不同压缩拐角特征流动区域内冷点效应对其表面热流分布的影响机制。计算结果分析表明：冷点效应的存在会显著地提高冷点区域内的表面热流值;在压缩拐角的分离域内,冷点效应造成的热流增长幅度较大,在再附点附近,冷点效应造成的热流增长幅度较小。     

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

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

壁温对压缩拐角流动影响的数值模拟研究

压缩拐角是高超声速飞行器上的典型非连续区域,其分离/再附结构对局部热环境有较大影响。本文采用基于Navier-Stokes方程的自研程序,开展了典型压缩拐角外形的气动热环境数值模拟研究,分析了不同壁温条件下压缩拐角的热环境分布规律。计算结果表明:随着壁温升高,流场整体结构变化不大,但由于近壁面流体物性参数变化较大,导致拐角分离涡分离点向前移动、再附点向后移动,拐角处干扰区扩大;压缩拐角大部分区域的热流随着壁温的升高而减小,且干扰区内热流减小的幅度比无干扰区更大,但热流并不完全遵循随壁温升高而减小的规律。另外,通过对比变壁温计算热流和热壁修正公式修正热流,发现热壁修正公式在干扰区存在精度降低、适用性不足的问题。本文的研究进一步加深了壁温对压缩拐角流动影响的认识,缩小了热壁修正公式的适用范围。     

压缩拐角区的压力,热流工程计算方法

 <正> 引言 从压缩拐角分离流动的试验及数值计算出发,建立起一个实用范围较广精度较高的工程估算程序,方便快速地估算出压缩拐角分离区和再附区的峰值压力、峰值热流;分离区和再附点位置;分离区长度等。这是气动设计部门关注的课题。这些工作是在文献(3,4)中的经验公式的基础上,按照分离流动流谱,分区进行迭代修正计算,对二维情况给出了完整     

侧壁约束下压缩激波/湍流边界层相互作用低频不稳定性实验研究

侧壁面约束条件下压缩激波/湍流边界层相互作用（Shock-Wave/Turbulent Boundary-Layer Interactions, SWTBLIs）呈现出区别于经典二维SWTBLIs的强三维特性，尤其是低频不稳定性问题自发现以来就受到了广泛关注，而其与下游压缩环境之间的关联性尚不明确。本文基于带侧壁约束的超声速压缩拐角构建三维受限的SWTBLIs流动，并通过自由射流风洞实验进行了系统研究，其中来流Ma∞=2.5、湍流边界层发展厚度对矩形流道宽度的占比约0.08、压缩拐角在12°~24°的宽范围内变化。通过不同压缩拐角下流场结构演化及壁面压力脉动信号谱特性的对比，揭示了侧壁面约束条件下压缩拐角对低频大尺度振荡时间-空间演化行为的影响规律和作用机制。研究发现：压缩拐角角度较小时，侧壁面约束导致的“拖尾效应”有助于延缓压缩拐角附近流动分离模式的过渡，并抑制压缩激波的不稳定性；压缩拐角达到20°及以上时，压缩拐角前沿逐渐发展为大尺度流动分离模式，侧壁面约束的角区将首先发展出能量集中的低频脉动，并逐步演化为两种峰值频率分别为约50Hz、200Hz量级的不同大尺度间歇性低频振荡叠加行为，且伴随着压缩激波形态的显著改变；两种低频振荡行为都属于SWTBLIs流场的全局性不稳定振荡，但沿着展向振荡强度差异显著，角区内的振荡强度远远高于中心主流区域。     

超声速二维湍流分离流动的数值解法

本文采用层流中引导简化N-S方程的方法,建立了简化的湍流运动方程组,分别用零方程、一方程模型计算了绕二维压缩拐角的湍流分离流动。数值求解用显隐式的反扩散格式。     

多尺度湍流模式在压缩拐角绕流中的应用

 应用4 个多尺度涡粘性湍流模式对一个二维压缩拐角流动进行了数值模拟,通过与实验结果以及两个基准模式(一个单尺度k-ε两方程模式和一个零方程模式) 的计算结果的比较,结合压缩性对湍能输运的作用、湍流的非平衡性对多尺度湍流模式在拐角非平衡湍流中的应用进行了评估。     

高超声速飞行器体襟翼局部分离流动数值研究

采用全Navier-Stokes方程的计算流体力学模拟技术研究了类HTV-2高超声速飞行器体襟翼局部分离流动特性,分析了不同飞行高度、壁面温度、飞行迎角等对流动分离特性的影响。研究表明:随着高度增加,壁面附近压力分布发生改变,沿着流向所形成的逆压梯度不断减小,使得体襟翼与飞行器表面附近的分离区减小;壁面温度增加导致压缩拐角(体襟翼与飞行器表面间)上游的流向速度梯度减小,即壁面粘性力减小,从而使分离区增大;随着迎角增加,压缩拐角内的逆压梯度增大,但上游流向速度梯度增加带来的黏性力增大更为明显,使流动更不容易发生分离,即分离区减小。     

高能电弧等离子体激励控制双压缩拐角激波/边界层干扰实验研究

超声速/高超声速飞行器进气道入口处多采用多级压缩构型,其诱导的激波/边界层干扰严重影响进气道效率和飞行器的气动性能,因此对双压缩拐角激波/边界层干扰进行流动控制具有较强的应用背景。在来流速度为Ma=2.0的风洞内,针对三种典型的双压缩拐角构型,开展了高能流向脉冲电弧放电阵列调控双压缩拐角激波/边界层干扰的实验研究,并对激励流场的高速纹影图像进行了空间梯度阈值处理和均方根处理。结果显示,在激励的作用下两道分离激波的强度均减弱,验证了利用高能流向脉冲电弧放电阵列控制双压缩拐角激波/边界层干扰的可行性。在分析控制效果的基础上,获得了在不同构型拐角的流场中前驱冲击波列和控制气泡的演化规律,结合控制效果的时序特征,最终揭示了高能流向脉冲电弧放电阵列作用于双压缩拐角激波/边界层干扰的前驱控制和接续控制的接力控制机理。     

一类新的速度入流条件在管道流模拟的应用

基于经典边界层理论,发展了一类简化的、按边界层厚度描述的速度入流边界条件。基于入流条件,对低速后向台阶流动和超声速压缩拐角流动进行了数值模拟。通过细致比较空间流场、边界层内速度型、表面摩擦阻力和压力分布等特性,对该类型边界条件进行了计算确认。结果表明,所发展的简化边界层厚度速度入流条件提法正确、合理,具有描述方法简单、鲁棒性好、适合工程计算的优点,可显著简化不同流速的管道流动的数值模拟。     

高超声速马赫反射

通过实验和数值计算模拟了激波绕４０度压缩拐角的流动。实验是在气动中心激波管中完成的，数值模拟采用了层流的Ｎ－Ｓ方程，对无粘项，粘 英和化学源项分别采用迎风ＴＶＤ格式，中心差分格式和点隐式的方法进行差分离散，计算了平衡与非平衡条件下的马赫反射流动，其数值计算结果与实验结果符合甚好。     

转捩对压缩拐角激波/边界层干扰分离泡的影响

为了研究转捩对压缩拐角内分离泡结构的影响,进行了来流马赫数2.9,24°压缩拐角激波/转捩边界层干扰的直接数值模拟（DNS）。通过在拐角上游平板的不同流向位置处添加周期性吹吸扰动激发流动转捩,使得转捩不同阶段进入拐角入口,从而在拐角内产生激波/转捩边界层的相互干扰。计算得到的平均速度剖面、壁面压力分布以及分离泡大小与风洞试验及以往直接数值模拟结果吻合较好,验证了计算结果的可靠性。研究了转捩过程对角部干扰区内分离泡结构的影响规律,分析比较了不同转捩阶段下角部分离区内湍动能的生成、耗散和分配机制。研究结果表明：转捩初期的拟序涡结构对分离泡尺度及形状影响最大,发卡涡包在角部拐点附近发生展向融合,并在角部区域形成湍流斑,此时分离泡尺度最小,形状呈现中间高两边低的山峰型。随着转捩的发展,分离区内湍动能生成和近壁区的耗散逐步降低,此时输运项起到了主要的平衡作用。     

一种抑制激波-边界层相互作用的新型无源被动控制

为了探索一种抑制激波-边界层相互作用的新型无源、被动控制方法,采用有限体积法计算流体力学程序,分别对压缩拐角和激波入射平板两种典型流动的附面层分离进行了无源被动控制的数值研究.数值结果表明,新型无源被动控制方法消除了20°压缩拐角流动所产生的分离泡;对激波入射平板所诱导的附面层分离,采用的自适应无源控制方法可将分离区长度减小为无控制时的58％、总压恢复可比无控制时提高7.2％、x=0.85m截面最大回流速度较无控制时减小了69.3％.从而验证了无源被动抑制控制激波-边界层相互作用的可行性.     

超声速压缩拐角激波/边界层干扰动力学模态分解

压缩拐角激波与边界层干扰问题广泛存在于高速飞行器的外部和内部流动中,其非定常复杂流场结构对飞行器气动性能影响显著。动力学模态分析将有助于进一步加深理解激波与边界层干扰流场不同特征频率对应的流动结构及动力学特性,为揭示其复杂流动机理提供参考。本文采用动态模态分解(DMD)方法对来流马赫数为2.9、24°压缩拐角内激波与超声速边界层干扰下的非定常流动进行了模态分析。评估了稀疏改进动态模态分解方法在压缩拐角流动中的适用性,研究了湍流干扰和转捩干扰下典型特征频率对应的动力学模态空间结构差异及其原因,分析了转捩边界层展向非均匀性对低频/高频模态动力学机制的影响规律。研究发现,湍流干扰与转捩干扰下拐角干扰区内均存在两类截然不同的动力学模态:低频模态和高频模态。低频模态结构集中在分离激波及分离泡剪切层的根部,表征为分离泡的大尺度膨胀和收缩运动;高频模态空间分布则以平均声速线附近正负交替结构为主,对应为边界层内不稳定波沿剪切层往下游的传播。转捩边界层的展向结构对低频模态运动特性影响明显,而对高频模态的影响则相对较小。     

存在烧蚀的高超声速压缩拐角流动数值研究

耦合求解热化学非平衡流控制方程和烧蚀壁面边界条件,进行存在石墨烧蚀的压缩拐角流场数值模拟。流场化学反应采用16组元(N2,O 2,NO,N,O,NO +,N+2,O +,N+,CO,CO 2,C,C2,C3,CN,e-)29个反应的非平衡模型,热力非平衡的双温度模型下,不同反应采用不同控制温度。石墨材料表面反应包括碳的氧化反应、碳催化的O 原子复合反应和碳的升华反应。对15°、18°、24°压缩拐角模型,在自由流 Ma =10~30,总焓值6~55 MJ/kg 范围,分别进行无烧蚀的壁面催化与非催化条件和石墨烧蚀条件下的流场计算,分析各类条件下的流场结构、流动分离特性以及流场热化学参数分布特点,研究壁面条件对流动特性的影响。结果表明：流动分离可能性和分离区范围随着压缩拐角斜面倾角增大而增大,随来流马赫数增大而减小;相对于低壁温条件,无烧蚀的辐射平衡壁温和壁面烧蚀条件下流动分离区增大,斜面上压力、摩阻和热流峰值点也有所后移。     

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

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

SSG模型在压缩拐角激波/湍流边界层干扰模拟中的应用

采用非线性方案模化雷诺应力再分配项的二阶矩雷诺应力输运模型——SSG模型和两方程线性涡粘模型SST模型,对24°高雷诺数二维压缩拐角激波/湍流边界层相互干扰流动进行数值模拟。通过在壁面压力分布、壁面摩擦阻力系数分布、截面速度分布及分离区大小等方面与实验数据比较,综合评估了上述两种模型模拟该问题的能力。结果表明:两种湍流模型在上述各方面都取得了与实验较为一致的结果,但相比较而言,除SST模型在流动再附后对壁面摩擦阻力系数分布的模拟更接近实验值外,SSG模型结果在上述其他方面均更接近实验值。本文工作为压缩拐角激波/湍流边界层干扰这类流动数值模拟的湍流模型选择提供了参考。     

高超声速风洞压敏漆试验技术

提出了适用于高超声速风洞开展压敏漆(PSP)试验研究的关键技术及解决办法。采用自主设计的PSP校准系统及测试系统,考核了代号为EC-PSP的压力敏感涂料在高超声速条件下的适用性、图像处理软件功能以及高温条件影响下数据处理方法的可行性。以压缩拐角模型为例开展了马赫数为5的高超声速PSP技术验证性风洞试验研究,辅以红外测温方法获得模型表面连续温度分布。试验结果表明在高超声速风洞开展的PSP试验技术研究清晰地捕获了基于压力变化的压缩拐角模型表面流动特征,实现了连续压力分布的测量。     

非线性EASM在激波与湍流边界层相互作用中的数值应用

采用基于k-ω湍流模型的非线性显式代数应力模型(EASM)对超燃冲压发动机常用的超声速凹槽、压缩拐角和侧壁压缩进气道简化模型的激波与湍流边界层的相互作用进行了计算,主要研究了EASM 模型对壁面压强、摩阻、Stanton数和壁面摩擦力线的计算精度,计算结果与SST(shear stress transport)模型进行...