中国CFD史

本文简要回顾了中国的CFD史,早在20世纪70年代,在钱学森教授的建议下,一些中国学者对计算流体动力学(CFD)这一新型学科产生了浓厚兴趣并开始研究.从那时起,中国的CFD逐渐发展繁荣起来,特别是从世纪之交,中国CFD非常活跃,为航空航天和其他民用领域做出了巨大贡献.在中国,特别注重CFD的研究与物理分析的结合.本文首先介绍了M5A的概念,它概括了中国CFD的主要研究领域;随后,介绍了设计格式的原则,以此原则为指导发展了一系列格式,比如NND格式、ENN格式和紧致格式.最近,发展了一系列高阶精度数值格式并应用于实际外形的复杂流动计算.再后,讨论了求解气动方程的网格尺寸准则.最后,同顾了如流动分离、流动拓扑结构和结构稳定性理论、旋涡沿其轴向演化和Hopf分叉理论、动导数和飞行稳定性等代表性工作.     

混合型三阶格式及关于两种DES算法的比较计算研究

为开展RANS/LES混合模拟,在传统三阶迎风偏置格式的基础上,提出一类混合型的三阶计算格式,格式在中高波数范围具有可调的耗散水平.在此基础上,文中通过计算比较研究两种RANS/LES混合算法:一种是基于Spalart-Allmaras一方程模型的DES模型(DES-SA),另一种是基于混合长度模型的DES模型(DES-ML),使用 发展的三阶格式对圆柱绕流进行模拟,并将得到的结果与文献比较,对计算格式、DES-ML算法进行了初步探讨.     

一类中心型三阶格式及应用

通过分析格式的耗散关系和重线性化技术,发展相应格式构造技术,提出了一类中心型的三阶非线性格式,格式具有与NND格式相当的网格模板.由于格式的线性形式具有很好的耗散特性,同时重线性化技术的使用使格式在光滑区尽量采用线性形式,从而理论上提升了计算对旋涡等结构的分辨和保持能力;另一方面,格式的非线性机制使格式在非光滑区退化为类NND格式,实现耗散和稳定机制.验证计算表明,格式能够开展粘性和旋涡及脉动流场的计算应用,算法形式简单,便于工程应用,容易进行有限体积推广.     

考虑平衡气体效应的乘波体特征线设计方法及气体模型影响

乘波体在高空高超声速飞行时,流场温度急剧升高,气体热力学性质发生改变,其气动特性与采用完全气体模型得到的结果存在差异,因此在飞行器设计中开展真实气体效应影响研究具有重要的工程价值。本文针对锥导乘波体设计提出了一种考虑平衡气体效应的特征线方法（Method of characteristic, MOC）,基于不同半锥角和不同流动捕获曲线生成了不同外形的乘波体,并使用快速估算方法估算平衡气体条件和完全气体条件下各乘波体的气动特性和几何特征参数,研究了不同气体模型对乘波体设计结果的影响。研究结果表明,高马赫数下基于平衡气体模型设计得到的乘波体升阻比、容积率和俯仰力矩系数等特性参数相比完全气体结果存在差异,并且平衡气体效应对不同外形乘波体的影响程度和影响规律不同。     

三维可压缩非定常流的壁面分离判据及其分离线附近的流动形态

通过引入依赖于密度的物面法向速度变换wr=-1/ρh1h2∫ozh1h2(e)ρ/(e)tdz,描述物理速度空间(ua,va,wa=u,v,w+wr)具有无滑移壁面条件的三维可压缩非定常连续方程可转换成变换速度空间(u,v,w)内具有无滑移条件定常连续方程.因此,采用定常壁面分离的分析方法和结论,再通过变换和研究wr的贡献,给出了三维可压缩非定常壁面分离的判则以及分离线附近的流动形态.研究指出,二维和三维情况下,都出现伴有壁外附着的壁面分离情况.数值模拟证实了理论和结论.     

基于旋转爆震三维流场结构分析的计算模型对比研究

为了研究不同计算模型对三维旋转爆震数值模拟的影响,分别基于Euler方程、N-S方程、RANS方法和IDDES方法并结合滑移、无滑移壁面边界条件,耦合氢气/空气的有限化学反应速率模型（7组分8基元反应）,采用高分辨率的五阶有限差分格式WENO-PPM5离散对流项,对圆环筒型燃烧室内的三维旋转爆震波进行数值模拟,得到了不同模型下旋转爆震波的流场结构、传播特性和推力性能。对比了不同计算模型对流场结构的影响。当采用滑移壁面边界条件时,Euler方程和N-S方程的计算结果较为一致。当使用无滑移壁面边界条件时,边界层的存在会导致可燃混气与燃烧产物之间接触面上的爆燃燃烧区域沿壁面向上游渗透,增大爆燃区域的范围;且爆震波锋面非受限侧变窄、严重变形,不同计算方法计算的变形程度有所不同。     

一种混合型四阶格式、基于特征的边界条件及其应用

本文描述并应用差分格式的构造原则,采用了时空转换的格式构造方法构造了一个四阶精度的差分格式,文中还进一步提出将四阶格式和二阶ＮＮＤ格式混合的方法,即在光滑区格式保持四阶精度、在激波附近降为二阶ＮＮＤ格式,来增强格式捕捉激波的能力。同时,本文基于特征传播的思想,建立了一种开边界的数值计算方法,并采用将边界条件嵌入控制方程的方法以及新的简化形式的边界处理方法,完整给出了无粘、粘性物面和开边界的的边界计     

三维方柱不可压缩绕流的大涡模拟计算

用非定常不可压缩流动求解方法对三维方柱绕流问题作了大涡模拟计算研究.求解程序采用虚拟压缩方法和双时间步的时间推进方案,对流项的差分采用高阶精度的数值差分格式.通过与实验和文献结果比较分析,验证了基于虚拟压缩方法的湍流大涡模拟计算方案对于复杂的三维钝体绕流计算是可行的.     

圆锥激波与平板边界层干扰的数值模拟研究

为进一步研究参数变化以及三维效应对圆锥激波/平板边界层之间相互干扰的影响,使用两方程Menter-SST模型,针对来流马赫数为2时的三维圆锥激波与平板边界层的相互干扰现象进行了数值模拟与定性、定量分析。分别研究了圆锥激波发生器半锥角和来流单位雷诺数变化对干扰区流动的影响,总结了参数变化引起的流动分离变化规律;此外,还计算了与三维计算的中心对称面上的入射激波等效的二维情形,并将三维结果与二维情形进行对比,对比结果显示中心对称面上的壁面压力系数、分离涡尺寸、涡量分布等与相应的二维情形存在明显差异。