空域空间网格标识原理及应用

在对空域空间有限元的相关研究中,构建空域空间的网格标识系统,是将信息技术应用于空域管理、推进空域数字化的重要基础。在美军全球区域参考系统和通用地理位置参考系统基础上,提出了一种同时具备标识平面和立体空间位置能力的参考系统。通过在地球表面以上特定高度之下的空间内,建立一种递归剖分策略方法,将空域空间划分为若干空间网格,不同尺度下的网格既是一种地理空间位置参考系统,也是一种组成空域的基本空域体,以此实现空间位置参考与空间位置标定的混合表达,为开展空域数字化管理提供底层模型。      

AUSM+-up格式在无网格算法中的推广

将AUSM -up格式推广运用到了无网格算法中.在点云离散的基础上,运用曲面拟合和AUSM -up格式求得数值通量;采用四步龙格-库塔方法进行时间推进,并引入当地时间步长和残值光顺等加速收敛措施.为了验证本文方法的计算精度和鲁棒性,对NACA0012翼型跨声速流动、某三段翼型低速绕流、平行错位NACA0012双翼超声速、高超声速流动进行了数值模拟.     

飞机结冰中水滴撞击特性欧拉法的网格影响分析（英文）

过冷水滴撞击特性计算研究是飞机结冰预测与防除冰系统性能分析的基础,欧拉法是过冷水滴运动及撞击特性计算的常用方法之一,分析欧拉法在计算复杂表面水滴运动撞击时的准确性具有重要意义。以NACA0012翼型、冰风洞风道、S形进气道及三段翼为研究对象,利用欧拉法采用不同网格计算获得了水滴运动及局部水收集系数,结果表明：当水滴运动没有受到上游部件的导流或者遮挡时,欧拉法计算结果具有网格无关性;当水滴运动受到上游部件影响发生轨迹偏转时,欧拉法计算结果具有网格相关性,不同网格计算得到的液态水含量及局部水收集系数不一致,需要进一步考虑水滴运动受到上游效应影响时网格对欧拉法计算结果的影响。     

随机场TPDF在非结构网格平台的开发及检验

为实现对燃烧室内部燃烧过程的高精度数值模拟,增强航空发动机燃烧室研究手段,采用开源的非结构网格平台Saturne,开发可以耦合详细化学反应机理,高精度解析湍流燃烧过程的概率密度函数输运方程（TPDF）湍流燃烧模型。采用随机场方法求解TPDF方程,在原有程序基础上研发了TPDF程序模块及配合的加速算法等模拟单元,发展了针对燃烧室燃烧过程模拟的软件功能。采用射流火焰Flame D和旋流火焰TECFLAM对新软件进行了测试,结果表明：模拟结果与试验数据趋势一致,精度尚可。对某型燃烧室性能模拟验证中,计算的出口温度分布与试验结果一致,平均温度误差小于3.8%,新的软件可进一步应用于燃烧室性能的研究。     

超音速分离线喷管大摆角状态下化学烧蚀动态仿真

针对超音速分离线喷管大摆角状态下化学烧蚀导致的壁面退移,基于动网格技术建立了相应的动态仿真模型,实现了对不同燃烧室条件下喷管化学烧蚀率的预示。初步稳态计算得到喉部烧蚀率为0.048 6 mm/s,高出试验结果5.67%,验证了仿真设置的合理性。以此状态结果为瞬态计算的初场,进行相应的化学烧蚀动态仿真计算。该喷管的矢量角放大系数在0.5 s仿真时间内因壁面退移减小了0.42%,对称面下侧分离线结构附近因燃气流动受阻成为喷管烧蚀最严重的位置,烧蚀率为0.074 5 mm/s。增大燃烧室压强或温度,会导致同周向位置的分离线后侧与前侧壁面烧蚀率比值减小。对于分离线附近型面变化较小处,压强5.5 MPa增加到7.5 MPa,该比值减小了10%,温度3200 K增加到3600 K减小了15%。     

高速进气道抛罩过程内流场演化特性研究

为探究进气道整流罩打开过程中内流场的演化规律,针对下颌式进气道及其配套整流罩,基于重叠网格法求解N-S方程,详细分析了抛罩过程中内流场演化过程与特征流场结构。研究表明,堵罩时,整流罩前端发展出V形激波,整流罩尾部形成后向台阶且产生一对涡矢量方向相反的流向涡,开罩过程中,整流罩前端的V形激波与尾部的后向台阶均会对内流场流动产生影响,捕获流量在中间某时刻达到峰值且大于稳定状态。对于开罩过程流动准定常且内流道无大分离区的流场,进气道在V形激波打在唇口时完成起动且整流罩在此后不影响内流场流动。假设的匀加速转动与重构出的真实转动吻合较好,对于更精确的运动过程模拟可采用三次函数的角速度变化律。     

多维对流通量重构方法的改进

为了研究多维对流通量重构格式中不同形式主方向对激波稳定性、压力场计算方面的影响,分别采用以笛卡尔方向、密度梯度及压力梯度作为主方向进行计算,对结果分析后发现,单纯依赖物理量梯度的主方向会造成压力等值线曲折的非物理解。为了消除这种非物理现象,采用一种基于压力函数的混合方式来确定主方向。单层及展向多层高超声速圆柱绕流问题模拟分析的结果表明:混合方式所确定的主方向不仅可以消除笛卡尔坐标系确定的主方向计算带来的红玉现象,同时还可以消除单纯的密度或压力梯度作为主方向所带来的压力场"zigzag"的非物理现象;通过与结构网格结果的比较表明:网格相关格式驻点处热流的相对误差为29.2%,而混合主方向确定方法所计算驻点热流相对误差只为4%,表明此改进方法计算得到的壁面热流是合理的。     

基于超网格的重叠网格守恒插值方法

保证重叠网格边界数据插值的守恒性是计算流体力学面临的一大难题。基于超网格方法与格心格式有限体积法,发展了一种新型混合重叠网格守恒插值方法。详细研究了用网格切割技术在重叠网格边界构造局部超网格,用网格求交算法合理地扩大贡献单元模板,建立了一种适用于任意网格类型的隐式并行重叠网格守恒插值方法,以超网格为媒介可实现重叠网格二阶精度的守恒插值。数值结果表明,本文方法对二阶分布的流场变量具有严格的守恒性,相比三线性插值方法和逆向距离权插值方法,本文方法减小了数值误差,提高了重叠网格边界的插值精度,加快了计算收敛速度,改善了重叠区域网格尺度相差较大时流场的光滑性和连续性。     

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

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

高效非结构网格变形与流场插值方法

非结构动网格在含动边界的非定常流动模拟中应用十分广泛。为了提高网格变形的计算效率和对复杂外形的适用性,发展了高效的K近邻-径向基函数（KNN-RBF）动网格方法。并采用高效的交界面节点一致性方法实现了网格变形的并行化。由于边界的运动和变形形式多种多样,在自主开发的计算流体力学（CFD）程序中实现了动网格的用户自定义函数（UDF）功能模块,可以灵活地设置此类问题中的边界运动。在模拟大变形问题时,采用了全局网格重构的策略对网格质量较差的区域进行处理。提出了快速查找新旧网格最近点的两级KD（K-Dimensional）树方法,并根据最近点在新旧网格间进行流场插值。从而建立了高效健壮的并行流场插值方法。通过对若干算例进行模拟,验证了所发展方法的正确性和健壮性。