用混合网格模拟绕X-38粘性流动

采用混合网格，利用DLR所发展的计算软件TAU对X-38高超声速粘性绕流进行了数值计算，网格自适应技术被用来生成有效的计算网格及提高对流场细节的分辨率。计算结果与现有的实验结果和用结构网格得到的结果吻合很好。     

高超声速电离绕流的数值模拟

本文从三维N－S方程出发，采用稳式L－U分解算法和7组分15个反应的化学模型，数值模拟高超声速电离空气绕流。首先采用对称TVD格式、AUSMPW＋格式和Van Leer的矢通量分裂格式计算了高超声速球头绕流，并对它们的计算结果做了对比分析。然后用前两种格式，对RAM－C飞行试验模型三个再入高度（81km、71km、61km）的流场进行了数值模拟，计算的流场电子数密度值和试验测量数据符合较好。     

真实飞行器无粘绕流直角切割网格数值模拟

本文采用直角切割网格，利用MGAERO无粘绕流解算器对三角翼真实飞行器带或不带有外挂物气动特性进行数值计算。计算马赫数为0．6～1．2。研究了网格密度对飞行器气动力特性的影响，同时与风洞实验结果进行了比较，符合程度令人满意。     

三维特体绕流数值模拟中薄层模型的应用

介绍用于近似计算绕三维尖无粘流激波薄层的二层模型。本法以数值确定激波的位置。在计算不同结构形状物体的空气动力特性时能确定气体动力学载荷的近似值。采用这种方法实际上排除了超高速理论中其他已知方法对尖体几何形状所加的约束。     

基于非结构网格的翼身组合体绕流数值模拟

应用投影变换的方法生成三维翼身组合体的表面网格，提出了新的三维空间多层背景网格生成方法，在背景网格建立和三维初始推进面生成的前提下．利用推进阵面法生成了翼身组合体的三维空间网格．应用网格光顺技术和变换方法，改进了初始生成网格的质量。应用格心格式的有限体积法对翼身组合体绕流进行了Euler方程的数值模拟，得到了较为满意的计算结果。     

用非结构动网格方法模拟有相对运动的多体绕流

发展了用于模拟有相对运动的多体非定常绕流，并计算由流场决定的外弹道的数值方法。该方法采用非结构动网格，显式NND有限体 格式，求解非定常Euler方程，并耦合了刚体动力学方程，作为例证，计算了俯仰振动的NACA0012翼型绕流，计算结果与实验及文献结果吻合，最后，整个策略用于计算二维机翼－外挂物分离流场，并确定了外挂物下落轨迹。     

高超声速无粘流场的推进计算

本文将文[1]中用于时间相关法计算的NND格式推广到定常超声速流动的空间推进计算,采用二步的预测、校正方法保证了推进方向的二阶精度,可以证明,这种二阶精度的NND格式具有TVD性质,是MacCormack二步显式格式的推广。本文首先将格式应用于二维平板上斜激波反射流场的推进计算,以检验格式捕捉激波的能力,同时研究了不同的通量分裂方法对格式捕捉激波能力的影响,得到了相当满意的结果。在此基础上,计算了航天飞机简化外形的身部超声速流场,给出了M_∞=10,α=0°,和M_∞=5,α=5°两种状态的部分结果,计算结果清楚地描绘了由于气流在机翼附近受到强烈压缩而产生的内嵌激波与外激波相交的复杂流场结构,与文[7]相比,流场结构更为清晰。     

用非结构网格和欧拉方程计算多段翼型绕流

用Ｄｅｌａｕｎａｙ法划分三角形构造非结构网格时，如果布点不当，会产生边界穿透，使网格生成失败，本文分析了几种布点情形，提出不容易产生边界穿透的布点方式。这一原理已用于生成绕多段翼型的非结构网格。     

粘流与无粘流的相互作用计算

本文总结了粘流/无粘流的各种计算方法和结果。重点在于介绍定常流动中的弱相互作用。首先叙述了弱相互作用的数学模型。给出了不可压流动和跨音速流动中粘流/无粘流相互作用的某些正耦合的计算结果。讨论了在分离区附近边界层正方法失效的原因。然后介绍了边界层反方法和适用于带分离的流动中半反方法耦合的粘流/无粘流的相互作用方法。文中也简单地总结了三维情况的应用和强相互作用。     

机动飞行器中等攻角高超声速无粘绕流数值模拟

本文采用ＮＮＤ差分格式，利用推进－迭代方法，数值模拟了机动飞行器的中等攻角高超声速无粘绕流流场，发展了一种适用于处理中等攻下背风面流场的技术，大量计算表明，结果准确、可靠。气动力力系数和壁面压力分布与实验数据一致，流场波系结构模拟正确。     

基于非结构动网格的非定常激波装配法

从ALE方程出发,通过修正体积的新思路构建离散几何守恒律,基于非结构动网格的信息传递方法得到网格重构以后新网格的流场参数,采用一阶格式捕捉法得到的流场确定装配法的初始激波位置,成功地把激波装配法应用到非定常流动的模拟。模拟超声速飞行环境下整流罩分离的非定常流固耦合过程,在马赫数6条件下装配法和捕捉法的计算结果完全一致,随着马赫数增加,二阶精度格式的捕捉法需要进行调整限制器参数、时间步长等人工干预,但是装配法直到马赫数20计算过程没有出现异常。该基于非结构动网格技术的激波装配方法,与传统的激波装配方法相比,能够用于复杂外形产生的非规则形状激波和非定常流动产生的运动激波,解决二阶精度有限体积法用于高马赫数飞行条件下多体分离过程的流动模拟时遇到的稳定性问题。     

带配平翼钝体高超声速粘性绕流的数值模拟

采用ＮＮＤ差分格式，通过求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，数值模拟了带配平翼钝体的高超声速粘性绕流。文中首先以球头驻点压强和热流为参数，讨论了在差分方程左端采用ＬＵ－ＳＧＳ、ＬＵ－ＡＤＩ和对角化ＡＤＩ三种不同隐式算法收敛效率的异同。然后以带配平翼钝体的高超声速粘性绕流为模型，对所研制的程序进行了计算验证。在计算中采用了代数方法和求解椭圆型方程方法相结合的网格生成技术，针对配平翼外形给出了贴体性、     

一类二维非结构网格DSMC方法的实现策略及其应用

研究了一类二维非结构网格DSMC方法的实现策略。在数据结构方面,设计了局部化的数据组织方式,节约了内存与计算时间。发展了一种跟踪模拟分子迁移的算法,该算法仅需少量的逻辑运算与代数运算,不仅可以快速跟踪模拟分子在网格之间的迁移,而且可以准确判别分子与物面是否相互作用,搜索过程中的附带信息给出了分子与物面碰撞的精确时间与位置,避免了重新计算。引入碰撞距离的思想,既减少存储又保证正确的模拟结果。在程序编制过程中,我们充分展现了Fortran90高级语言的主要特性,引入动态数组、指针、链表以及派生类型数据,编制了计算程序。最后对过渡流域高超声绕流进行了数值试验。     

非结构网格下可压缩黏性流体的隐式算法

采用牛顿线化和Roe近似Jacobians矩阵隐式算法,实现基于非结构网格的三维定常黏性流动问题的求解.分别对平板、跨声涡轮叶片Mark-Ⅱ和CW-22这三种具有代表性的流动进行数值模拟,并详细分析了流动细节,模拟结果与理论以及实验值吻合良好,准确预测了流场内激波、二次流等重要的流动现象.验证了基于非结构网格的可压缩流体隐式计算方法对于黏性流动问题的计算是正确的.      

网格之间流场信息传递的高精度方法和验证

理论分析了插值方法引入误差的原因,提出一种利用时间推进方法和移动网格技术实现网格之间流场信息传递的思想,并构造计算方法。在介绍新的信息传递算法原理的基础上,根据三维情况下传递过程需要满足的时间推进步长稳定性条件,推导出最小推进步数的要求。最后,通过对激波诱导物体运动的非定常流动模拟,比较了信息传递算法和线性插值、加权最小二乘插值方法,计算结果表明,与插值方法相比,信息传递过程没有引入误差,能有效抑制非物理振荡,保持流场计算精度。     

绕凸台的高超声速分离流动研究

本文研究高超声速流动绕三维低凸台的流动特性。凸台高度与边界层厚度之比介于0.5～0.8,凸台周边倾角介于14°～45°。在高超声速风洞中,来流马赫数为5,单位雷诺数为2.6～6.0×10~7/米。实验过程中测量了模型中心线及特殊部位之表面压强分布,根据纹影记录及表面流动显示确定激波与边界层干扰流场特性及分离区的变化。发现分离区在中心线上距凸台最远,而在凸台肩部距周边沿法线方向最近,确定了中心线上最大压强比的位置。     

在非结构动网格上追踪多介质界面的ALE方法

在处理带有介质分界面的多介质流动问题时,把介质边界当作一种特殊的边界,并且在该边界上的点都有着双重的定义,分别对应边界两侧的不同介质状态,通过利用HLLC方法求解任意拉格朗日-欧拉方程(ALE),以及运用界面跟踪和动网格方法分别对一维和二维多介质流进行了数值模拟,并得到了满意的结果.     

基于非结构混合网格的CHN-T1标模气动特性预测

本文使用自行研制的基于非结构混合网格的亚跨超声速流场解算器MFlow,对AeCW-1提供的客机标模CHN-T1进行了数值模拟研究。介绍了非结构混合网格的生成情况,重点分析了网格收敛特性、压力分布、气动特性曲线、流动分离等。计算得到了近似线性的网格收敛特性。随着网格加密,对激波和分离气泡的模拟更精细。尾支撑对气动特性的影响非常明显,特别是对平尾气动特性有很大影响。机翼静气动弹性变形的影响主要是使升力系数和阻力系数减小。湍流模型的QCR修正对大迎角计算结果有较大的影响。计算结果表明,MFlow程序能够准确地预测客机标模的气动特性。