矩形/非结构混合网格技术及在二维/三维复杂无粘流场数值模拟中的应用

建立了一套模拟复杂无粘流场的矩形／非结构混合网格技术,其中非结构网格仅限于物面附近,发挥非结构网格的几何灵活性,以少量的网格模拟复杂外形;同时在以外的区域采用矩形结构网格,发挥矩形网格计算简单快速的优势,有效地克服全非结构网格计算方法需要较大内存量和较长ＣＰＵ时间的不足。混合网格系统由修正的四分树（２Ｄ）／八分树（３Ｄ）方法生成。本文将ＮＮＤ有限差分与ＮＮＤ有限体积格式有机地融合于混合网格计算,消除了全矩形网格模拟曲面边界的台阶效应,同时保证了网格间的通量守恒。该算法被成功地运用于ＮＡＣＡ００１２翼型的跨、低超声速绕流,二维弯管内运动激波的非定常流,以及三维捆绑火箭的超声速绕流的数值模拟。计算结果表明本方法在模拟复杂无粘流场方面的灵活性和高效性。     

烧蚀外形转捩流动的数值模拟

本文采用雷诺平均Navier-Stokes方程和修正Baldwin-Lomax代数湍流模型,以Roe二阶流通量差分分裂格式进行离散,数值模拟了典型烧蚀外形的高超声速流场,重点讨论了流动转捩的判定及其与复杂流态的关联.结果表明:数值模拟出的球钝锥外形转捩点位置与工程估算结果符合较好;对工程估算中存在困难的奶头状凹陷烧蚀外形,数值模拟也可给出符合流场结构的转捩点位置.     

专用工具材料定额管理系统的开发

本文简要介绍了计算机专用工具材料定额管理系统的设计原则和主要功能,及其计算机辅助管理产生的效果。     

三角翼大攻角绕流场数值模拟研究

本文基于N-S方程采用Jameson的显式中心有限体积格式和Roc的上风格式对三角翼大攻角绕流场进行初步数值模拟，空间网格为H-0型网格。计算结果表明两种计算格式均可较好的模拟大攻角绕流场。本文计算结果为今后开展全机大攻角绕流场数值模拟提供了计算经验。     

多重网格法求解离心压气机内部流场

本文提出了一种新的多重网格法,并将之与显式时间推进法、有限体积差分格式、粘性体积力法配合使用求解离心压气机内部流场。数值计算表明,作者多重网格法的计算结果与试验结果吻合比较好,同时该多重网格法的收敛速度比不用多重网格法加快了将近两倍。     

基于Steger-Warming矢通分裂的NND格式

将Steger-Warming矢通分裂方法和NND格式相结合构成一种基于Steger-Warming矢通分裂且能自动捕捉到激波的差分算法.通过对激波管和对可燃气体DDT过程的数值分析,证明该算法解决了Steger-Warming矢通分裂方法在零特征点附近存在数值解震荡的问题,有着良好的计算精度和计算效率.文中采用的计算方法对各种有关矢通分裂计算方法的实现有一定的参考价值.     

六阶精度的群速度直接控制紧致格式及其应用

差分数值解产生非物理振荡的直接原因是在于非均一的波群的群速度,因此本文利用直接群速度控制的方法重新构造具有六阶精度的紧致型差分格式,以达到改善激波数值解的目的。对所构造的格式的精度及数值解随控制参数变化的行为进行了分析。最后文中给出了一些典型算例,证明了文中所构造的格式具有方法简便,物理含义清楚,精度高,网格基架点小,和捕捉激波能力较强的优点。      

UNIX下报表技术研究及报表工具开发

提出并讨论了UNIX下报表工具需要解决的三个核心问题:报表数据,报表格式和报表输出。构造了一个报表工具模型。在该模型中,引入了报表数据定义语言和报表格式描述语言,用相应的解释器解释它们。最后,以文中的报表工具模型为基础,实现了一个UNIX下的报表工具。     

湍流模型离散精度对数值模拟影响的计算分析

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用五阶空间离散精度的加权紧致非线性格式(WCNS),通过改变物面法向第一层网格间距,开展了剪切应力输运(SST)两方程模型不同离散精度的数值分析。主要目的是为高阶精度格式在复杂外形上的应用提供技术支撑。计算模型包含了低速NLR 7301两段翼型和高速RAE2822翼型,研究内容主要包括湍流模型的二阶精度离散和五阶精度离散两种方式对收敛历程、边界层湍流黏性系数分布、边界层速度分布、压力系数分布以及气动特性的影响。在与试验数据对比的基础上,计算结果表明:对于不同的第一层物面法向网格间距,湍流模型离散精度对低速绕流计算结果有比较明显的影响,对于高速小迎角附着流动计算结果影响不明显;相对于湍流模型二阶精度离散,湍流模型高阶精度离散网格敏感性较弱,具有更高的数值模拟精度,但收敛性略差。     

基于超声速边界层转捩的WENO格式对比分析

基于高精度weighted essentially non-oscillation(WENO)格式对马赫数为4.5,雷诺数为10000的超声速来流条件下平板边界层转捩过程进行了大涡模拟.无黏通量分别用5阶、7阶、9阶WENO格式进行离散,黏性通量离散采用4阶中心差分格式,时间推进采用具有总变差递减(TVD)性质的3阶精度Runge-Kutta方法.通过在入口边界叠加一对等幅最不稳定第一模态谐波扰动,分别采用3种WENO格式计算得到了平板层流边界层失稳转捩的演化过程.结果表明:5阶WENO格式的数值耗散明显高于7阶和9阶WENO格式,在捕捉湍流涡和流场脉动特性上存在明显不足.9阶WENO格式的耗散小,能够捕捉到流场更小尺度的涡和高频脉动.研究具有高频脉动特性的问题或者欲捕捉精细涡结构时,建议采用7阶以上的高精度格式.