超声速绕钝板熵层不稳定性的研究

为了研究熵层的不稳定性,首先选择马赫数为4.5、前缘为圆头的0°钝板模型,利用直接数值模拟得到基本流场,然后在此基础上进行线性稳定性分析,同时通过经验方法确定熵层和边界层区域。研究发现,对于某一流向站位的基本流剖面,存在两个广义拐点,分别位于边界层和熵层。针对熵层的广义拐点,进行稳定性分析研究了熵层中不稳定模态特性,得出该模态的不稳定区域集中在前缘附近,且扰动增长率很小。通过直接数值模拟熵层扰动的演化,发现在前缘附近扰动幅值也是缓慢上升,之后缓慢下降。另外,距离前缘较远处的数值模拟得到的扰动幅值与LST结果比较,二者吻合得很好,说明熵层在下游具有很好的平行性。     

内埋式武器舱的流动及气动特性分析

探讨了战斗机内埋式武器舱的一些研究进展,对武器舱形状、来流马赫数、边界层厚度、噪声等对武器安全投放的影响进行综合分析,重点讨论了武器舱几何形状对流动特性的影响;基于N-S方程,采用数值模拟方法,分析了几种不同长深比二维武器舱的流动特性,以及三维武器舱锯齿舱门对气动特性的影响;探讨了噪声特性及噪声抑制技术,并对武器投放中...     

真实气体效应对升力体舵面局部流动分离的影响

数值模拟分析了高马赫数低雷诺数条件下激波边界层干扰、激波与激波相互作用、流动分离再附等流动现象的特点以及高温真实气体效应的影响。分别采用量热完全气体、平衡气体、化学非平衡气体模型对升力体由于舵面偏转引起的局部流动分离情形进行了数值模拟。研究了飞行高度、壁面温度及来流马赫数对流动分离的影响。计算结果表明:真实气体效应使空气在边界层内发生离解反应,边界层内温度降低,粘性减小,动能损失减小,克服逆压梯度的能力更强,从而使分离区明显减小。分离区的减小改变了分离/再附激波的位置和强度,进而对局部压力及热流分布产生重要影响;随高度增加,平衡气体较完全气体分离区相对减小量增大,平衡气体效应对流动分离/再附现象的影响越大;壁温对分离区影响较大,随壁温升高,分离区增大;随马赫数增大,分离区减小,真实气体和完全气体的差异增大,真实气体效应的影响更加显著。     

可压缩剪切层与边界层流动相互作用的稳定性研究

为了探究可压缩剪切层与边界层组合流动稳定性规律,采用线性稳定性理论对这种组合流动的线性稳定性问题进行研究。平均流动模型取双曲正切函数描述的可压缩剪切层与可压缩相似边界层速度剖面的直接叠加。采用整体数值方法 (Global)求解组合流动线性稳定性方程,获得全部不稳定模态。结果表明,在这样的组合流动中,存在两个不稳定模态,一个与剪切层相关的不稳定模态,另一个是与边界层相关的不稳定模态。在此基础上,研究了剪切层与壁面之间的距离、扰动频率、雷诺数、马赫数、速度梯度、温度梯度对这种组合流动稳定性的影响。马赫数增大对不稳定模态起抑制作用,剪切层对边界层的不稳定模态的促进、抑制作用与马赫数、扰动频率有关,边界层对剪切层不稳定模态基本上起抑制作用。     

熵层对高超声速二维钝楔气动参数的影响

采用薄激波层理论和流管质量守恒相结合的方法,分析了高超声速二维钝楔边界层外缘熵的分布规律,研究了熵层对边界层外缘密度、马赫数以及壁面气动加热等气动参数的影响.结果表明:边界层外缘熵在倒圆-肩部区下降最为剧烈,熵层对气动参数的影响在高超声速下不可忽略,特别是使转捩区和湍流区的气动加热增加约53.6%.因此,将表面流态控制在层流模式对高超声速飞行器热防护具有重要意义.      

Influence of entropy layer for hypersonic on aerodynamics parameters 2-D blunt wedge

采用薄激波层理论和流管质量守恒相结合的方法,分析了高超声速二维钝楔边界层外缘熵的分布规律,研究了熵层对边界层外缘密度、马赫数以及壁面气动加热等气动参数的影响.结果表明：边界层外缘熵在倒圆-肩部区下降最为剧烈,熵层对气动参数的影响在高超声速下不可忽略,特别是使转捩区和湍流区的气动加热增加约53.6%.因此,将表面流态控制在层流模式对高超声速飞行器热防护具有重要意义.     

不同波系配置的鼓包压缩面引起的激波/边界层干扰特性

为了对不同波系配置的鼓包压缩面引起的锥形激波和来流边界层干扰的流场特性进行对比分析,分别以总偏转角相同的单锥、双锥和等熵锥等3种外锥流场为基准流场,在来流马赫数为2.0的条件下,采用流线追踪法设计了3种鼓包压缩面。在此基础上,采用数值仿真的方法对这3种鼓包压缩面在不同马赫数下的三维激波/边界层干扰流场进行了对比分析,结果表明:相比于单锥鼓包,双锥和等熵鼓包流场中的涡量集中区明显减小,且涡核处的涡量值也有所降低;对于双锥和等熵鼓包,由于受到后续波系的干扰,其表面流动的拓扑结构不再呈现出锥形相似性;双锥与等熵鼓包具有和单锥鼓包相同的边界层排移能力,且这一能力不随来流马赫数的变化而变化。此外,采用多波系配置能提高鼓包外压缩系统的总压恢复系数,波系数目越多,提升的程度就越大,且这一优势在高马赫数下更加明显。      

飞行参数对射流矢量喷管内流场影响的数值模拟

采用三维雷诺平均N-S方程和k-ε湍流模型对固体火箭发动机燃气二次喷射矢量喷管内外流场进行数值模拟,研究飞行参数对矢量喷管内主/次流耦合流场的影响。结果表明,自由来流造成喷管出口附近压力变化,通过喷管内壁面附面层作用影响喷管内部流动特性。飞行马赫数为超声速,或者环境压强较低时,侧向力大小与飞行参数无关。在低空亚声速流中,侧向力随飞行马赫数增加迅速减小,随环境压强降低而减小。     

高超声速层流气动热预测混合算法研究

提出了兼顾激波和边界层模拟的混合算法--在激波区逼近NND格式,在流动光滑区逼近三阶偏置格式,采用了Steger-Warming、KFVS和Vanleer三种典型的矢通量分裂方法,对来流马赫数为10、基于头部半径的雷诺数为1×105和1×106的球头绕流的壁面热流进行了计算.计算结果表明:混合算法明显降低了热流对网格的依赖性.在此基础上完成了钝双锥和双椭球壁面热流的预测,得到了与实验测量一致性比较好的结果.     

等离子体控制边界层流动仿真研究

使用泊松方程、漂移-扩散模型和N-S方程对介质阻隔面放电等离子体控制边界层流动进行了一体化数值模拟。仿真结果表明介质阻隔面放电等离子体通过推-拉空气可以对低速气流进行控制,正向推力和反向拉力均可以在距壁面一定高度处形成射流,正向射流自身增加边界层动量,更重要的是通过引射主流将主流中的能量输运到边界层中,从而减小边界层厚度;反向射流与来流相撞后形成回流区将边界层向外推,能够增加边界层厚度。     

激波边界层的相互作用对扰动波传播的影响

利用线性稳定性理论和直接数值模拟研究了带有入射斜激波的、来流马赫数Ma=4．5条件下的平板边界层的失稳特性。重点考察了在由于激波边界层相互干扰,平板边界层上形成分离区,又进而产生激波、膨胀波和旋涡等复杂流动现象的流场上游,引入小扰动的TS波后,扰动波传播通过带有这些复杂流动现象的流场时,扰动波的发展变化特点。通过对流场中扰动波（包括基本波和衍生波）演化特征的分析,研究分离和激波等复杂流动现象对平板边界层稳定性的影响特点。数值模拟发现,激波的出现不同于一般的压缩波,在亚声速区与超声速区对扰动波演化的影响是不同的,此外,分离对扰动有稳定作用。     

γ-Re_θ模式应用于高速边界层转捩的研究

应用γ-Reθ转捩模式对超声速、高超声速边界层转捩进行数值模拟。γ-Reθ模式通过求解关于当地雷诺数和间歇因子两个输运方程给出转捩起始位置和转捩区长度等信息。本文对马赫数3.5至7范围内的四种算例进行计算,研究了来流雷诺数、攻角变化和头部钝化半径等关键参数的变化对γ-Reθ模式预测流动转捩性能的影响。给出了壁面摩擦阻力系数、热流值与实验值的对比以及壁面附近间歇因子等值线分布等计算结果。γ-Reθ模式能够正确预测攻角和钝化半径变化时转捩位置和流动参数的变化趋势;在较大雷诺数时计算结果与实验值吻合很好。γ-Reθ模式对于超声速或高超声速边界层转捩的模拟仍需修正和改进。     

等离子体控制边界层法流动仿真研究

使用泊松方程、漂移-扩散模型和N-S方程对介质阻隔面放电等离子体控制边界层流动进行了一体化数值模拟.仿真结果表明介质阻隔面放电等离子体通过推-拉空气可以对低速气流进行控制,正向推力和反向拉力均可以在距壁面一定高度处形成射流,正向射流自身增加边界层动量,更重要的是通过引射主流将主流中的能量输运到边界层中,从而减小边界层厚度;反向射流与来流相撞后形成回流区将边界层向外推,能够增加边界层厚度.     

等离子体合成射流对钝头激波的控制与减阻

为了验证等离子体合成射流对超声速流动的流动控制和减阻效果,在考虑热完全气体效应的情况下,工程拟合等离子体热力学属性和输运系数,利用能量源项模型对超声速平板和球头等典型流场结构进行了数值模拟,并取得了与实验较为一致的结果。研究结果表明：对于马赫数为2的超声速平板流动,等离子体合成射流能有效干扰边界层的发展,并诱导产生一系列大尺度结构;对于马赫数为3的球头流场,等离子体合成射流能显著改变激波脱体距离与球头的阻力特性;在放电后第1个周期内,合成射流能使球头平均阻力降低6.3%,而在射流峰值情况下使阻力降低32.0%,激波脱体距离增加2倍,实现了激波控制和流动减阻的预期效果。     

进口边界层对采用弯叶片的平面扩压静叶栅流场性能的影响

给定不同型式的进口边界层,在两种不同亚音速条件下对一平面扩压静叶栅的弯叶片流场进行了数值模拟。结果表明弯叶片对扩压叶栅的改善的能力受进口边界层的特征影响。这种影响分为两个方面:(1)边界层厚度的影响和(2)边界层动量损失厚度的影响。边界层越厚或动量损失厚度越大,在低马赫数条件下弯叶片对吸力面角区密流增加越明显,从而更大程度地提高了端区的流动性能,降低了叶栅损失。在高马赫数条件下,若边界层越厚或动量损失厚度越大,角区密流虽变化不大,但因端区损失较大,其性能的提高会给叶栅总性能的改善带来较大的收益。      

基于边界层转捩的高超声速进气道特性研究

为了探索边界层非强迫转捩对进气道性能的影响,采用数值计算的方法开展了边界层转捩对轴对称混压式高超声速进气道流场特性的研究。研究表明:随着进气道中心锥锥尖钝化半径增大,边界层转捩先推迟。当锥尖钝度大到一定程度时,边界层转捩位置前移。随着钝化半径进一步增大,边界层转捩再次推迟,转捩位置逐渐后移。来流湍流度越大,边界层越不稳定,边界层转捩越易发生。与湍流边界层相比,考虑边界层转捩时进气道的总压恢复系数及流量系数较高、热载荷及阻力系数较小,Ma=6.5时喉道处总压恢复系数最高上升17.3%,进气道阻力最大下降17.4%。边界层转捩对壁面热流密度分布影响较大,但对壁面压力分布影响较小。钝化影响进气道的自起动性能,随着钝化半径增大,自起动马赫数升高,而边界层转捩对进气道自起动性能影响较小。     

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

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

高超声速钝头体边界层转捩试验

为研究高超声钝头体边界层转捩以及头部钝度对转捩的影响,在FD-14和FD-14A两座激波风洞中开展了热流、压力扰动和高速纹影显示等综合测量。试验结果表明,转捩雷诺数关于钝度雷诺数的变化显示出转捩反转的趋势。压力扰动的功率谱密度（PSD）分析结果以流向离散分布云图形式显示,边界层高速纹影图像显示了第二模态波的发展、湍流的生成和熵层对边界层结构的显著影响。大头部钝度带来的强熵梯度熵层流动对边界层压力扰动频谱特性和流动结构影响显著,在转捩反转机理中起到重要作用。此外,马赫数对转捩的影响不容忽视。     

典型构型空腔模型设计与流动/噪声特性研究

空腔高速流动包含着复杂波系、旋涡、剪切层等空气动力学典型流动特征,空腔研究中的模型设计、来流马赫数和边界层等参数模拟,以及典型空腔流动/噪声计算与试验结果的差异、缘由等问题亟需解决。综合参考国内外空腔模型参数选取和设计方案,在分析前期前缘平板尖劈构型和角度对来流边界层形态和厚度的影响规律的基础上,通过改变空腔前缘平板尖劈构型和角度来解决上述参数模拟不准和引起计算与试验结果差异等问题,提出了一种典型结构平板-空腔标模(C201)的参数选取和设计方案,给定前缘平板尖劈角度为5°,并利用理论分析、数值计算和风洞试验对该模型设计方案进行了考核验证。综合利用表面流谱、空间流场结构和静/动压测量方法等,获得了亚跨超声速条件下(Ma:0.6~2.0)C201空腔流动/噪声的基本试验数据,并深入分析了该模型的流动/噪声特性及来流马赫数、攻角等参数的影响规律,可为数值计算方法的验证、空腔流动/噪声机理分析和控制方法的构建等提供基本验证数据和借鉴依据。     

预处理方法在混合笛卡尔网格中的应用研究

使用基于混合笛卡尔网格的预处理方法,对低马赫数下的定常/非定常流动问题进行了数值模拟研究。网格生成中,在物面附近生成贴体结构网格,其余计算区域使用笛卡尔网格进行填充,两套网格之间通过查找"贡献单元"实现流场数据间的传递。同时,使用基于密度的预处理方法,发展了一套可求解从低速(极低马赫数)到常规马赫数的N-S方程求解器,时间离散使用隐式双时间步LU-SGS格式,空间离散使用具有二阶精度的格心格式有限体积方法。非定常流动问题的数值模拟过程中,使用逆距离插值来进行新现网格单元上流场参数确定。通过对NACA0012翼型在低马赫数下的定常绕流和动态失速问题的数值模拟研究表明:使用预处理方法能够加快低马赫数下数值计算的收敛速度和提高计算的准确性,基于混合笛卡尔网格的N-S方程求解器能够模拟包含运动边界的低速不可压非定常流动问题。耦合预处理技术的混合笛卡尔网格方法给包含低速和常规马赫数的复杂运动边界问题的求解提供了新的思路。