高超声速飞行器二元进气道试验和计算

设计了一种吸气式面对称高超声速飞行器,针对进气道性能,分别在两座风洞开展通流试验研究。针对第1次风洞试验大攻角状态(α=8°)测量值偏离线性的问题,辅助采用数值模拟手段分析原因,并对试验方案进行改进设计,解决了首次试验出现的问题。结果显示,在典型状态(Ma=5~6)下,进气道起动正常,性能良好,具有一定的抗侧滑能力;随来流马赫数增加,进气道流量系数增大,总压恢复系数减小,计算结果和试验结果一致;试验结果和数值计算的差异主要表现为基本测压方案α>4°后,流量系数和总压恢复系数出现严重的非线性。数值模拟结果表明,主要原因为模型支撑方式及测压方式所引起的偏差,通过改进试验方案,解决了大攻角状态下测量值偏离正常趋势的问题。     

火箭整流罩外气动噪声环境的大涡模拟研究

基于五阶加权本质无振荡（WENO）格式构造隐式大涡模拟方法（ILES）,对跨声速来流条件下（Ma=0.8）火箭整流罩外噪声环境进行数值模拟。通过与风洞试验结果及国外文献进行对比,ILES方法能够在较粗网格下准确预测壁面湍流脉动特性。跨声速流动在壁面折角处出现分离、再附、激波/边界层干扰现象,均方根脉动压力系数出现峰值,同时该区域噪声能量在全频段都较高,易引起结构抖振效应。最后,根据ILES模拟结果,指出工程常用的外噪声经验公式的不足,并提出改进措施。     

洞映射方法的研究和改进

针对传统洞映射方法存储大的缺点,对"广义封闭"的概念进行扩展,提出了最小洞映射方法,该方法允许挖洞曲面结束于网格截断面,有效缩小了洞映射区域,节省了存储开销.对适用于广义封闭的洞映射单元识别方法进行了分析,指出了用物面信息判断映射单元属性可靠性较差.发展了一种新的特别适合广义封闭问题的识别方法:"Inverse mark",使用计算网格结点信息自动识别洞外单元,再作为初始点在网格内部推进.研究表明:"Inverse mark"方法计算效率高,可靠性好,自动化程度高,有效提高了重叠网格方法对缝隙等局部复杂结构的适应性     

火星进入器作强迫震荡运动壁面脉动压力数值模拟

为研究超声速阶段进入器作强迫震荡运动对壁面脉动压力环境的影响规律，本文耦合进入器刚体运动方程与流体力学方程，采用动网格技术，对火星进入器模型开展非定常数值模拟，获取壁面不同位置处的脉动压力信息。研究表明：在超声速阶段，进入器作强迫震荡运动诱导的脉动压力远大于进入器保持相对静止时仅由非定常流动诱导的脉动压力。来流马赫数为1.2时，进入器作强迫震荡运动对脱体激波影响较小，脱体激波强度较弱且形态变化较小，攻角的震荡导致同一测点距离脱体激波的位置发生周期性改变，舱体迎风面及配平翼迎风面的脉动压力环境主要受攻角变化的影响；来流马赫数为3时，进入器作强迫震荡运动对脱体激波的影响较大，脱体激波震荡剧烈，诱导舱体迎风面及配平翼迎风面产生极其恶劣的脉动压力环境，功率谱分析表明激波震荡诱导的脉动压力能量主要集中在30 Hz左右。     

火星进入器壁面脉动压力环境数值模拟

为研究跨超声速阶段来流马赫数、来流攻角及配平翼展开角的变化对进入器壁面脉动压力环境的影响规律，本文采用脱体涡方法对火星进入器模型开展非定常数值模拟，获取壁面不同位置处的脉动压力信息。研究表明：在跨超声速阶段，进入器壁面脉动压力环境随马赫数的增加而趋于减缓。配平翼迎风面分离区受脱体激波影响明显，当来流马赫数较小时，可压缩效应较弱，分离区涡流运动剧烈，诱导的脉动压力环境较强；随着来流马赫数的增加，脱体激波对分离区抑制作用增强，分离区运动受到限制，诱导的脉动压力环境趋于平缓。此外，随着来流攻角增加，配平翼迎风面上再附点的位置向翼根方向转移，从而使翼根处的脉动压力环境趋于恶劣。当配平翼展开180°时，分离区再附点位置基本固定，配平翼迎风面脉动压力环境得到一定程度的减缓。功率谱分析表明，在配平翼迎风面上诱导的脉动压力能量主要集中在中低频区域。     

鲁棒的结构网格自动化重叠方法

网格技术是目前数值模拟中的关键技术之一。重叠网格是一种放宽拓扑要求、减小网格生成难度的网格技术。本文以结构重叠网格为基础,分别针对挖洞、寻点以及洞面优化方法进行了研究和改进,同时完成物面网格重叠,形成了一套鲁棒的、自动化的网格重叠系统。在挖洞方面,结合“最小洞映射”方法,提出“复合式挖洞”方法,节省内存开销;在寻点方面,通过构建格心虚网格,保证搜索空间的连续性,同时结合“有效搜索”思想,排除部分对寻点无贡献的网格点,进而减少ADT叉树节点;在洞面优化上,改变填补判别法则并引入两类受保护洞内点,确保两层插值边界建立,提高鲁棒性;在物面网格重叠上,利用物面投影法完成坐标修正,实现物面附近网格流动变量的准确传递。为验证本文方法,分别对定常翼身组合体DLR-F6绕流和非定常机翼挂载分离过程进行了数值模拟,计算结果与实验结果吻合良好,表明该结构重叠网格系统对多物体间定常、非定常扰流具有较好的数值模拟能力和较高的模拟精度,具有较高的工程应用价值。     

俯仰阻尼导数分量的CFD数值模拟

提出一种直接求解直接阻尼导数的方法,该方法不仅适用于轴对称外形,也适用于非轴对称外形.数值模拟飞行器的非定常强迫沉浮运动和强迫角振动,并在Etkin的非定常气动力模型基础上,辨识得到飞行器的洗流时差导数和俯仰阻尼导数,研究飞行器俯仰通道各阻尼导数的数值计算方法.对弹道外形和基本带翼导弹标模外形及Hyflex升力体外形进行研究分析表明,将阻尼导数分量相加得到的俯仰阻尼导数与直接求解强迫角振动得到的俯仰阻尼导数与试验结果吻合很好,各阻尼导数分量随质心位置的变化趋势也与理论预测相符.对于带翼飞行器,超声速条件下,洗流时差导数在俯仰阻尼导数中占主导作用.      

新型重叠网格洞面优化方法及其应用

高度自动化的洞面优化方法是重叠网格关键技术之一。通过对割补法的分析,针对其洞面切割后洞边界位置不确定性及不可预知性问题,提出将其与物面距离优化准则相结合的新型混合洞面优化方法。该方法针对各物体间体网格,采用物面距离优化准则,将洞边界位置控制在物体中间的位置,使重叠区域更合理、可预测;针对背景网格,将其物面距离设置为较大值,采用物面距离优化准则将其落入物体网格内部的单元一并挖去;针对物面处重叠的网格,采用割补法对其进行切割和填补,保证物面及近壁面网格得以有效的重叠。该方法在保证高适用性的前提下提高了网格重叠质量,自动化程度高,无需人工干预。3个典型复杂流动算例计算结果与实验结果吻合良好,网格重叠区域流场变量传递正确,等值线过渡光滑,流场刻画准确,证明该方法准确可靠。     

一种改进的跨声速旋成体壁面脉动压力经验预测公式

使用大涡模拟（LES）方法对典型整流罩旋成体结构进行非定常数值模拟研究,结果表明引起壁面压力脉动的流场结构随时间和空间变化,不同的流场结构之间相互干扰。同时研究发现,旋成体折角膨胀会造成流动在较短的距离内迅速加减速,同样引起壁面压力脉动。基于上述流动现象,重新划分压力脉动产生区域,通过引入空间修正项,提出一套改进的跨声速旋成体壁面脉动压力经验公式。算例表明,该方法能够显著提高旋成体在跨声速条件下的脉动压力分布的预测精度。     

基于虚网格的格心ADT搜索法

针对重叠网格寻点问题,提出了一种直接在格心网格下操作,基于虚网格的ADT(Alternating Digital Tree)搜索方法.通过对计算网格进行扩展,建立虚边界网格,解决了格心网格因覆盖区域不完整而无法直接建立ADT数据结构的问题.搜索结果即为包含解的格心单元集合,可直接对结果列表遍历以得到合理贡献单元,故完全摒弃了可靠性差的Stencil Walk方法.由虚网格的定义,使寻点在边界附近的处理更为灵活,可以准确给出边界附近贡献单元的有效信息,同时简化了虚网格体系的构建.扩展了重叠边界类型,构建的搜索空间能完整覆盖网格范围,解决了对称面重叠问题.算例研究表明:该方法可靠性好,边界处理能力强,有效提高了重叠网格方法对大型复杂网格的解算能力.