多孔表面抑制第二模态失稳的最优开孔率和孔半径分析

多孔介质、微槽道和超声波吸声材料等可用于抑制高超声速边界层第二模态扰动波。通过声波在无穷长小管道中的传播模型给定多孔介质表面的扰动边界条件,采用时间模式的线性稳定性理论分析了多孔介质表面对边界层稳定性的影响。对马赫数6热壁平板边界层的考察表明,多孔介质表面不仅可以大大推迟第二模态扰动波的中性点,还可以大大抑制最不稳定第二模态扰动波的幅值增长率。为了找到最优控制参数,在较广的参数范围内考察了多孔介质表面的开孔率和孔半径对第二模态扰动波幅值增长率的影响,得到了最优开孔率和孔半径沿流向的分布。最后,还考察了基本流的当地非平行性(法向速度)对控制效果和最优控制参数的影响。     

不同迎角下脊形前体绕流数值模拟研究

脊形前体有较强的背风涡流场,不同的前体形状对前体涡流场和气动力有很大的影响。本文针对脊形前体飞行器大迎角湍流大分离流动计算的困难,采用IDDES混合湍流模型,以及与之匹配的非定常算法,研究了不同来流迎角下脊形前体的气动特性,以及背风涡非定常演化、破裂的细致流动结构。选取了不同脊形角,以及不同上、下高宽比的脊形前体进行计算。计算结果表明,在迎角较小时,随着迎角的增大,前体主涡会逐渐增强,在迎角较大时,前体主涡破裂;在相同迎角下,脊形角较小时,前体涡较强,涡升力也更大;对于相同脊形角的前体,当上半截面高宽比较小时,前体主涡强度较大,前体涡破裂临界迎角较小,即会提前破裂。     

高超声速6°迎角圆锥边界层背风流向涡稳定性分析

针对马赫数6、迎角6°圆锥边界层背风面流向涡开展了直接数值模拟和稳定性分析,其中,稳定性分析工具包括空间BiGlobal方法和面推进的PSE3D方法。研究发现,有迎角圆锥边界层背风中心线附近会自发形成一对对转的强流向涡。这对流向涡或因低频反对称模态和高频对称模态的共同作用而在下游破碎。     

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

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

基于RANS/LES混合方法的分离流动模拟

飞行器在大迎角、快速俯仰机动时,流场中含有大尺度、非定常的涡结构,传统雷诺平均Navier-Stokes （RANS）模型不能准确模拟流场结构,根据国际上相关研究的发展趋势,需要采用混合RANS/大涡模拟（LES）模型来对复杂分离流动进行准确模拟。本文对基于分区混合与湍流尺度混合的双重RANS/LES混合计算模型进行发展与应用。通过典型简化模型的静、动态湍流大分离流动,测试和验证所采用的脱体涡模拟（DES）类方法,重点研究改进的延迟DES （IDDES）模型在动态问题应用中的正确性和有效性,并对所采用的数值模拟方法和相应的计算软件的可靠性、鲁棒性以及精度进行了考核验证。典型算例包括超声速圆柱底部流动、跨声速方腔流动、NACA0015机翼深失速分离涡模拟等。计算表明：发展的IDDES类混合计算模型可有效解决对数层不匹配的问题;对于定态非定常分离流动,DES、DDES、IDDES等模型计算结果差别不大,随着流动的非定常特性增强,IDDES模型的优势逐渐显现;对于动态非定常分离流动,则需要采用IDDES类模型。     

MF-1弹道设计与蒙特卡罗飞行仿真

MF-1是我国的首次高超声速空气动力学基础问题研究模型飞行试验。本文研究建立了MF-1试验飞行器弹道设计与拉偏仿真数学模型,开展了MF-1弹道设计和拉偏仿真计算。对结果分析表明,对于MF-1这类无控飞行器,某些条件下可能在飞行过程中发生滚转共振现象,导致飞行总迎角迅速发散,气动过载显著增加,飞行弹道严重偏离设计弹道,甚至可能出现解体风险。通过严格控制飞行器加工和安装偏差,有助于降低滚转共振发生概率。     

方形截面导弹摇滚特性数值研究

采用隐式耦合求解RANS方程与滚转运动方程的方法研究了方形截面导弹的摇滚运动特性。耦合求解过程采用隐式处理提高了求解精度,采用强迫摇滚得到了自由滚转静导数和动导数;采用与摇滚方程的耦合求解得到了方形截面导弹模型在不同迎角下的自由摇滚历程。两种研究结果均显示随着迎角的增加,方形截面导弹将由滚转稳定变为不稳定,并形成极限环振荡,即摇滚;随着迎角进一步增加将形成滚转发散。两种分析方法所取得的结果与试验结果吻合较好。     

前体涡非对称分离机理及前缘吹气控制研究

通过设计对称性算法,求解层流Navier-Stokes方程,数值模拟了细长体在低超声速情况下前体背风涡随攻角演化的规律,在此基础上,进一步研究了前缘吹气对前体涡演化和侧向力特性的控制机理.根据数值模拟结果及分析,倾向于支持在层流框架内,前体涡的非对称失稳是一种对流不稳定机制,要想根据需要产生对称或不对称的前体涡,就必须外加持续的扰动.在约16°～48°攻角区间内,前缘吹气可产生规律性较好的侧向力,有可能直接利用前体涡进行横侧向控制.为工程实用化,需提高前缘吹气的激励收益效费比.     

马赫6柱-裙构型激波/湍流边界层干扰摩阻统计特性

为探究激波/边界层干扰表面摩阻统计特性,对马赫6柱-裙构型的激波/湍流边界层流场进行了直接数值模拟,推导了柱坐标形式的时均摩阻分解公式,将其与对流效应、流向不均匀性效应、分子黏性耗散效应、曲率效应和湍动能耗散关联起来,在此基础上对摩阻脉动和时均摩阻的统计特性进行了研究。概率密度分布结果表明干扰区内摩阻脉动概率密度明显偏离正态分布,摩阻脉动在该区域的间歇性较强,这与湍流脉动在激波干扰后存在显著内在可压缩性效应有关;谱功率密度结果表明摩阻脉动的能量在干扰前中频部分的能量占据主导,峰值对应频率约为0.14,干扰区后附近区域峰值位置转移到高频区域。采用摩阻分解公式对激波干扰前后的时均摩阻进行分解,结果表明时均摩阻的主要贡献来自分子黏性耗散效应和湍动能耗散效应;在激波干扰前分子黏性应力效应占据主导位置,干扰后湍动能耗散效应占据主导位置;激波干扰后对流效应明显增强,曲率效应的贡献在激波干扰前后基本保持不变;流向不均匀项贡献由正贡献变为负贡献,且所占比例有所上升,这主要是由压力梯度项的变化导致的。对分解后各项积分内函数沿法向分布的分析表明,分解后各项在边界层内分布呈现较大的差异,特别地,湍动能耗散项...     

火箭发动机尾罩热分离流场的数值模拟

潜射导弹的发动机一般用尾罩加以保护以防止海水侵蚀,发动机点火后尾罩存在热分离过程,本文应用非定常数值模拟方法和动网格技术数值模拟了发动机尾喷管超声速流场的建立过程,分析了高温燃气与尾罩碰撞后反射流动对发动机防护罩的非定常载荷,计算结果表明本文计算方法可以用来数值研究广泛存在的气流瞬态冲击效应。