带壁面射流燃烧室燃烧数值模拟

利用基于非结构化网格有限体积法对三维有壁面射流的燃烧室内两相流动和燃烧进行了数值研究.对气相流动在Euler坐标系下求解,而对液滴相则利用Lagrange方法进行追踪求解.计算区域采用四面体网格进行划分,气相流场用SIMPLEC计算方法,对液滴相采用了欧拉隐式方法.考虑了液滴相与气相的完全双向耦合作用,分别采用了Spalding液滴蒸发模型和涡破碎(EBU)燃烧模型,数值计算结果与文献中实验数据吻合较好.      

多重网格加速的LUSGS算法用于喷管流场数值模拟

采用有限体积法,二阶精度插值的Roe格式对三维流场N-S方程进行空间离散;采用LUSGS隐式时间推进法和多重网格技术加快计算收敛速度;采用独特的显式算法克服了两方程湍流模型刚性对流场计算收敛速度的影响,并给出了该方法的理论依据。通过轴对称喷管内流场、二元喷管内外流场的数值模拟验证了采用以上方法在计算效率上的提高。最后将RNG模型计算结果与SA模型计算结果及实验数据做对比分析。     

基于非结构重叠网格的二维N—S方程求解与应用研究

本文发展了非结构重叠网格技术.通过构造任意单元类型网格重叠的网格间边界确定算法、宿主单元搜索算法以及网格间插值算法,把非结构重叠网格推广应用于求解NS方程,并结合湍流模型模拟二维湍流流场.N-S方程的求解采用有限体积法,通量计算采用Osher格式,湍流模型采用SA模型.为验证本文方法的正确性和有效性,模拟了RAE2822翼型和某三段翼型绕流湍流流场,结果表明采用非结构重叠网格模拟复杂外形湍流流场不但网格生成效率高,而且容易推广应用于具有相对运动的非定常问题,具有较高的工程应用价值.     

关于变形网格“几何守恒律”概念的讨论

在国内外文献中应用包含有变形的动网格进行流动数值模拟时,要满足所谓的几何守恒律条件;通过对这一概念最早提出推导过程和近年来应用发展情况分析以后发现,逻辑上存在不适当之处.从流体力学基本理论出发,对包含有变形的动网格技术的理论基础进行讨论,证实几何守恒律是流体力学控制方程的伴随方程或退化方程;通过简单模型分析有限体积方法离散过程,发现计算过程中不满足几何守恒律所引起的非物理现象本质上是目前有限体积离散处理中微元体积计算方法不符合物理定律所致;因此在变形动网格计算方法中不存在必须要满足的几何守恒律.最后根据理论提出物理背景清晰的微元体积计算方法,数值验证可行.     

运动边界非定常全速流场的数值模拟

将预处理引入到ALE格式的非定常可压缩N-S方程中,发展了一种能有效模拟含运动边界的三维非定常全速流场的数值方法.控制方程在非结构动态网格上用双时间步推进并用隐式、LU-SGS迭代求解;提出了一种特征变量的边界条件以适合预处理后的特征系统.几个典型算例的结果表明,预处理方法将可压缩流方法的计算能力拓展到了低马赫数范围,有效地提高了收敛速度和计算精度.     

低速增压风洞支架干扰的数值研究

对气动院FL-9低速增压风洞尾撑支架干扰进行数值计算.通过合理的拓扑结构,生成了YF-16模型的分区对接网格.计算对模型的典型安装状态进行研究,获得了全机气动特性和支架干扰量,探讨了此安装状态下弯刀对飞机气动力的干扰特性,并重点研究了Re数对干扰量的影响规律.     

带稳定板装置弹射座椅偏航稳定性能研究

针对国内弹射救生领域现状,设计了一种适合于国内某型座椅的弹射救生系统稳定板。运用计算流体力学(CFD)数值计算方法,对改装稳定板前后的人椅系统进行了高速流场下的静态和动态特性分析,并与风洞试验结果进行对比,两者吻合较好。改装稳定板前后的数值计算对比表明:在不对俯仰和滚转性能产生不利影响的情况下,稳定板的改装对人椅系统的偏航性能有比较明显的改善。六自由度动态计算比较真实地模拟了弹射出舱的过程,也进一步验证了稳定板对改善偏航稳定性能所起到的作用。数值计算完全可以应用于弹射座椅的气动研究,同时稳定板的设计对工程上的应用具有一定的参考价值。     

适于黏性计算的自适应笛卡儿网格生成及其应用

复杂外形的黏性网格生成已经成为计算流体力学(CFD)在工程应用中的主要瓶颈,高效而鲁棒的网格生成技术研究显得尤为迫切。使用基于物面外形自适应加密的笛卡儿网格方法生成计算网格,投影方法拟合壁面边界,将投影得到的柱形单元沿法向分层获得适于黏性计算的网格。应用特征恢复技术恢复凹面特征并改善凸面处的网格质量。流场解算中应用最小二乘法重构获得2阶精度,Roe格式计算无黏通量,Holmes-Connell方法计算黏性通量。加入Spalart-Allmaras一方程湍流模型以获得湍流解。上下对称高斯 赛德尔迭代(LU-SGS)格式隐式时间推进,对网格重排序保证格式的稳定性和平衡性。算例显示网格生成方法是快速的和鲁棒的,流场计算精度满足工程需求,因此该方法非常适合应用于实际工程计算。     

预设堵块法检测进气道自起动能力的数值研究

在进气道内预先设置轻质堵块迫使进气道不起动,堵块被气流吹出后流道恢复畅通,为激波风洞提供了一种检测进气道自起动能力的方法。为了深入认识预设堵块方法的检测过程,将堵块简化为一个自由度的刚体运动,采用k-ωSST湍流模型,结合铺层动网格技术,对该检测方法进行了二维非定常数值模拟。在Ma_∞=5.9条件下,采用改变堵块质量的方式,获得了三种典型的检测过程。通过分析堵块在气动力作用下的运动与进气道非定常波系演化相互耦合的过程,揭示了堵塞作用产生和消除的流动机理。结果表明,堵塞时间没有显著地改变进气道的自起动能力,预设堵块方法成功检测到进气道在Ma_∞=4.9自起动。此外,建立了预估堵块运动的匀加速模型,预估的堵塞时间与数值模拟结果较为符合。在Ma_∞=5.9条件下,堵块的临界质量约为3.8g。考虑到激波风洞实验时间短暂,应选择合适的堵块。     

基于阶跃响应的飞行器单独静、动导数数值求解方法

基于阶跃响应方法,结合刚性动网格技术,对飞行器的单独静、动导数的精细化数值计算进行了研究.以纵向为例,通过给物面施加恒定附加攻角,求解得到阶跃响应运动过程的非定常气动力,求导得到静导数.同样给物面施加恒定的俯仰角速度,并同时强迫物面平动以抵消俯仰转动产生的附加攻角影响,可由非定常气动力求导得到动导数值.分别利用NACA0012翼型和三维SACCON飞翼无人机进行了计算验证,各攻角下的静、动导数值与文献、试验结果吻合得很好,最大误差不超过5%.结论表明:基于阶跃响应的单独静、动导数直接模拟方法计算耗时仅为传统强迫振动方法的21%,效率相对较高,且可推广到横航向的动导数计算,为飞行器的稳定性研究提供参考.