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

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

膨胀效应对激波/湍流边界层干扰的影响

采用直接数值模拟方法对来流马赫数2.9,30°激波角的入射激波与膨胀角湍流边界层干扰问题进行了数值研究。入射激波在壁面上的名义入射点固定在膨胀角角点,膨胀角角度分别取为0°、2°、5°和10°。通过改变膨胀角角度,考察了膨胀效应对干扰区内复杂流动现象的影响规律和作用机制,如分离泡、物面压力脉动特性、膨胀区湍流边界层和物面剪切应力脉动场等。研究发现,膨胀角角度的增大使得分离区流向长度和法向高度急剧降低,尤其是在强膨胀效应下分离泡形态呈现整体往下游偏移的双峰结构。物面压力脉动功率谱结果表明,膨胀角为2°和5°时,分离激波的非定常运动仍表征为大尺度低频振荡,而膨胀角为10°,强膨胀效应极大地抑制了分离激波的低频振荡,加速了下游再附边界层物面压力脉动的恢复过程。膨胀区湍流边界层雷诺剪切应力各象限事件贡献和出现概率呈现逐步恢复到上游湍流边界层的趋势,Görtler-like流向涡结构展向和法向尺度变化剧烈,同时在近壁区将诱导生成大量小尺度流向涡。此外,物面剪切应力脉动场的本征正交分解分析指出,膨胀效应的影响体现在低阶模态能量的急剧降低从而使得高阶模态的总体贡献相对升高。     

展向振荡对激波/湍流边界层干扰的影响

周期振荡作为一种有效的壁面流动控制手段受到广泛关注,而其对激波/湍流边界层干扰的影响目前鲜有研究。本文采用高精度直接数值模拟（DNS）方法对马赫数2.9、12°激波入射角、强振荡下的激波/湍流边界层干扰进行了系统研究。通过与无振荡工况的定量比较,揭示了展向强振荡对干扰区内复杂流动结构的影响规律及作用机制,如分离泡尺度、物面压力脉动非定常特性、物面剪切的非定常特性及统计特征等。研究发现：在展向强振荡作用下,分离点位置提前,间歇区长度增大;同时由于分离泡内强黏性耗散的影响,展向振荡的穿透高度约为分离泡高度的4%,因而对流动结构不会产生实质影响。但展向强振荡会对壁面附近流动造成显著影响,如强振荡诱导的壁面展向速度远大于流向速度,造成流向剪切与展向剪切之间夹角的概率密度函数峰值从0°偏移到80°~90°之间。物面压力及剪切本征正交分解分析表明,展向振荡会导致模态能量从低阶模态向高阶模态转移,降低低频运动的能量占比,增强再附后Görtler涡等壁面附近旋涡结构的强度。     

烧蚀外形方程差分计算方法研究

分析了烧蚀外形方程的性质,研究了人工粘性项对烧蚀外形计算结果的影响,发现人工粘性在抹平非物理波动的同时,会严重影响计算精度.借鉴于CFD中NND格式的思想,构造了求解烧蚀外形方程的无波动、无自由参数的耗散差分格式,得到了满意的计算结果.     

高超声速有攻角锥裙直接数值模拟

采用直接数值模拟方法对有攻角的高超声速7°～34°锥裙开展了数值研究,通过对比0°、90°、180°周向子午面,评估了三维横流效应对激波/边界层干扰的影响规律和作用机制,包括壁面压力、摩阻、热流分布,分离泡非定常运动,再附边界层演化等。研究发现,不同周向站位均出现流动分离,横流区、迎风区内复杂激波结构与边界层相互作用导致壁面压力、摩阻、热流显著升高。热流与压力的比值在干扰区上升后由于再附降低,而热流与摩阻的雷诺比拟关系在分离区则完全失效。分离泡面积脉动的功率谱结果表明,分离泡非定常膨胀/收缩运动呈低频特征,且分离泡呼吸与激波低频振荡在横流区密切相关,在迎风区存在迟滞,而在背风区不相关。速度脉动场的本征正交分解结果表明,分离区的低频特征与低阶模态相应的剪切层附近大尺度结构相关。对下游再附边界层演化分析指出,攻角的存在导致雷诺应力在再附点附近大幅增强,其流向分量的恢复最为迅速,雷诺应力分量的峰值位置在背风区沿流向持续外移,而在迎风区、横流区已迅速向内层恢复。此外,雷诺应力各向异性不变量分布进一步表明干扰下游的近壁区湍流各向异性峰值在背风区弱于迎风区。     

碳/碳化锆复合材料烧蚀机理和计算方法研究

研究了碳/碳化锆复合材料氧化烧蚀机理,发现它们与传统的硅基和碳基材料烧蚀有很大差别。基体ZrC氧化后在表面形成一种膨松状的多孔固态抗氧化膜,氧化膜能有效阻止材料进一步氧化,使烧蚀量大大降低。研究了烧蚀过程抗氧化膜的形成、演化和流失行为,研究了氧化膜中氧气的扩散机制(包括分子扩散和Knudsen扩散),研究了材料原始层表面可能存在的化学反应,建立了分析碳/碳化锆复合材料烧蚀响应的物理数学模型。     

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

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

C/SiC复合材料烧蚀机理和通用计算模型研究

研究了C/SiC复合材料氧化烧蚀机理,发现它们与传统的硅基和碳基材料烧蚀有很大差别。C/SiC的烧蚀取决于氧的分压、表面温度和材料晶态结构及成份,可能出现活性氧化和惰性氧化两种破坏机制。研究了氧化膜的形成和破坏条件,以及氧化膜中氧气的扩散机制,建立了适用于C、Si、SiC和C/SiC(C和SiC可有不同混和比)四种材料烧蚀计算的通用物理数学模型。     

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

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

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

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