武器分离及舱门开启过程数值模拟研究

基于非结构重叠网格技术,发展了一套能较好模拟空腔流动和武器分离问题的数值计算方法.通过模拟典型空腔流动和WPFS模型分离轨迹,并与试验结果比较,说明本文方法可较好模拟空腔复杂流动和武器分离轨迹.采用发展的数值方法对某内埋弹舱舱门开启过程和内埋武器分离轨迹进行了数值模拟,研究了舱门运动的非定常效应对舱门气动特性的影响.研究表明:一定条件下,舱门运动产生的非定常效应对舱门气动特性影响较大,舱门和舱门控制系统设计时必须考虑；在本文给定条件下,弹舱内的内、外侧导弹均能安全分离.通过研究,为内埋弹舱的舱门设计和内埋武器分离安全评估提供了技术手段和依据.     

面向非结构混合网格高精度阻力预测的梯度求解方法

针对非结构混合网格的特点,通过改进传统的Green-Gauss梯度求解方法,提出了一种可提高非结构混合网格黏性计算精度的节点型Green-Gauss梯度求解方法。利用改进后的方法,完成了DLR-F4翼身组合体算例的计算和对比分析。改进后的梯度求解方法残差收敛更好,下降量级更多,阻力系数和试验吻合更好,激波区域压力分布和分离区域流场细节的模拟更精确,说明改进后的梯度求解方法有效提高了程序的鲁棒性和阻力预测精度,验证了方法的有效性。采用改进后的方法对第5届AIAA阻力预测研讨会的通用研究模型（CRM）进行了详细的模拟分析,结果表明：改进的梯度求解方法更加适用于非结构混合网格的黏性计算,计算精准度达到国际同类CFD软件水平,进一步验证了改进方法的可靠性。     

超声速空腔流动数值模拟研究

采用Roe格式、S-A一方程湍流模型和LU-SGS隐式算法对三个不同类型的空腔进行了模拟,并把数值计算结果和试验结果进行了比较,二者一致.通过对流过不同L/D空腔的可压缩流动进行模拟,得到了不同L/D空腔的底部压力分布、空腔内流线、等压线和等马赫数云图,分析了空腔流动类型随L/D增大转变的机理.     

分离涡模拟方法在内埋弹舱机弹分离模拟中的应用

为了研究内埋弹舱非定常流场对武器投放分离特性的影响,基于分离涡模拟（DES）方法和非结构重叠网格技术,建立了内埋武器机弹分离高精度数值模拟方法。在此基础上,首先通过典型空腔流动和机弹分离算例验证了数值方法的有效性,然后采用平均雷诺应力模拟（RANS）和DES两种方法探究了机弹分离过程中弹舱腔内涡结构的演化历程,分析了两种方法在机弹分离流场结构模拟中的差异,研究了内埋弹舱流场结构对导弹分离特性的影响。计算结果表明：机弹分离过程中,高强度涡结构逐渐向弹舱腔体后缘聚集,同时弹舱剪切层流动结构被导弹破坏,导弹出弹舱阶段俯仰方向所受气动力矩急剧变化。相比RANS方法,DES方法捕捉到的高强度小尺度涡结构对导弹压力分布影响较大,使得分离后期导弹俯仰角、俯仰角速度产生明显差异,DES方法更适合内埋武器舱分离特性精细化研究。该研究结果可为内埋武器机弹分离研究的计算方法选择和安全分离策略提供一定的参考依据。     

基于SA湍流模型的DES方法应用

使用基于SA(Spalart-Allmaras)湍流模型的分离涡模拟(DES,Detached-Eddy Simulation)方法,对圆柱绕流和超声速圆柱底部流动等高雷诺数大分离流动进行了数值模拟.无粘通量的计算采用二阶迎风格式,时间推进使用二阶隐式双时间步法.计算表明相对于雷诺平均N-S方程(RANS,Reynolds-Averaged Navier-Stokes)模拟,时空均为二阶精度的DES模拟能够分辨更为精细的分离区三维旋涡结构,对分离区主要湍流结构的捕捉和平均压力分布的预测取得了满意的结果.     

脉冲爆震燃烧室内缓燃向爆震转捩数值模拟

为了研究脉冲爆震发动机燃烧室内火焰加速及缓燃向爆震转捩过程，利用7组分8反应的氢气详细化学反应机理进行 了2维数值模拟。结果表明：在火焰传播的初始阶段，障碍物、旋涡与火焰的相互作用是主导火焰加速的主要因素。在燃烧区域 产生的膨胀波向前传播并诱导未燃混气流动；障碍物后的回流区及障碍物顶端脱落的旋涡使火焰面拉伸、卷曲，增加火焰面面积， 同时提高燃烧放热强度，使火焰传播速度加快。在火焰发展的后期阶段，激波与障碍物的相互作用会使脱落在已燃混气中的未燃 混气微团发生局部爆炸，产生新的激波进一步推动火焰加速，缩短火焰锋面与前导激波的距离。由于交错型障碍物增加了火焰传 播的距离，其缓燃向爆震的转捩时间比采用对称型障碍物时的更长。     

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

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

火箭子级垂直回收布局气动特性及发动机喷管影响

垂直回收可重复使用运载火箭是运载火箭发展的一个重要方向,大长细比的火箭子级垂直再入过程属于典型的非规则钝头体绕流,与传统低阻力流线体飞行器气动特性差异较大。采用风洞试验辅以数值仿真分析的方法,对基于栅格舵的火箭子级垂直回收构型基本气动特性和非规则钝头体绕流情况进行了研究,获得了发动机外露喷管和栅格舵对火箭子级垂直回收气动特性的影响规律,给出了火箭子级垂直回收布局设计建议。结果表明：火箭子级倒飞状态下肩部区域会在小迎角下产生大分离流动,外露发动机喷管左右两侧诱导出较强的分离涡结构,与火箭尾翼、肩部大分离流动相互作用;垂直回收构型在超声速下阻力会一直处于较高的水平,不同马赫数下压心移动量较大,倒飞时发动机外露喷管会产生较大的干扰静不稳定力矩,其量值与栅格舵提供的静稳定控制力矩基本相当,在火箭子级垂直回收方案设计时需要引起注意。     

雷诺数对飞机气动特性的影响研究

采用CARDC自主开发的基于非结构混合网格技术的大规模并行解算器MFlow,对运输类飞机和战斗机气动特性的雷诺数影响规律进行了研究。通过对不同翼型气动特性的雷诺数影响进行数值模拟研究,探讨了雷诺数对运输类飞机和战斗机气动特性影响不同的原因,得到了翼型厚度和弯度是雷诺数影响气动特性的关键参数的结论。     

面向混合网格高精度阻力预测的熵修正方法

Roe通量差分分裂格式具有耗散小、分辨率高等特点,在亚跨超声速流场模拟中得到广泛应用,但空间离散时需要对Roe平均矩阵的特征值进行熵修正,使用熵修正会增大耗散,影响阻力的求解精度。本文针对非结构混合网格中Roe格式熵修正的特点,通过改进传统的Roe格式Harten-Yee熵修正方法,提出了一种可提高非结构混合网格黏性计算精度的Harten-Yee熵修正改进方法。利用改进后的方法,完成了DLR-F4翼身组合体算例的计算和对比分析,改进后的熵修正方法残差收敛特性与原始Harten-Yee熵修正一致,计算精度提高,计算结果和无熵修正时基本一致,说明改进后的熵修正方法既保留了使用熵修正带来的程序鲁棒性等优点,同时把熵修正对阻力预测精度的影响降到了最低,验证了方法的有效性。采用改进后方法对第3届AIAA阻力会议的计算模型DLR-F6翼身组合体进行了详细的模拟,分析了网格收敛性和雷诺数的影响。结果表明：改进的Harten-Yee熵修正,更加适用于非结构混合网格的黏性计算,计算精准度达到国际同类CFD软件水平,进一步验证了改进方法的可靠性。