尖拱弹身横向喷流数值模拟

本文通过数值求解层流Navier—Stakes方程来模拟超声速流场中横向超声速喷流的干扰流场，采用分块对接网格和“O”型网格技术，精确模拟喷口截面形状并生成高质量的贴体计算网格，对尖拱弹身的横向超声速喷流的干扰流场进行了数值模拟，研究了喷口附近的流场结构、涡系结构、波系结构及流动分离等流场特性，得到喷流流场的干扰效应，计算结果与实验值吻合较好。     

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

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

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

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

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

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

无陀螺捷联惯导系统捷联方案研究

提出了无陀螺捷联惯导系统的惯性测量单元沿单轴相对载体旋转的捷联方案,对该方案进行了理论分析和数值仿真。结果表明,合理安装９ 个加速度计并适当选取捷联矩阵,则捷联方案无需高精度测量惯性测量单元相对载体旋转角速度也能大大提高载体角速度的精度。为满足相同的导航精度要求,本方案可降低对加速度计精度的要求约１０００ 倍。     

高精度差分格式在机翼分离流计算中的应用

用一个模型方程分析了边界层方程的数值稳定性 ,指出稳定性问题会随计算雷诺数减小而变得严重 ,因而对较低的Re数 ,计算分离边界层流动采用高阶精度差分格式十分必要。本文同时给出用四阶精度差分格式求解在流动分离情况下边界层方程的技术方法 ;对一个大展弦比后掠翼在攻角等于 1 4°时计算了机翼上表面的分离线 ,并与文献上公布的用别的测算方法测出结果做了比较 ,两者总体上符合良好     

稀疏非结构网格上的亚声速流高精度数值模拟

采用高阶间断有限元法在非结构网格上数值求解二维亚声速Euler方程。数值结果表明,尽管采用的非结构网格非常稀疏,但通过采用真实曲线物面边界和高阶的基函数仍可得到高精度的数值解。另外,对于超低速情况,方程无需经过任何特殊处理就可以得到收敛的数值解。由于采用牛顿一般最小余量法(Newton-GMRES)时通常需要较好的初始值,本文设计了一种阶谱循环过程来提高数值求解的鲁棒性。     

可压缩流动的高分辨率低耗散混合格式研究(英文)

针对可压缩流场的计算,提出了一种混合型激波捕捉格式来获得对于光滑区域和流场间断都具有高精度的解。该混合方法采用有限体积形式的五阶WENO格式来计算流场的光滑区域,而采用MUSCL格式来捕捉流场中的间断。该混合格式结合WENO格式的低耗散特性和MUSCL格式的高分辨率特性。两种格式之间通过一个间断感受因子进行切换。文中选择了三种典型的间断感受因子进行对比研究。WENO格式和MUSCL格式均为激波捕捉格式,因此混合开关中的常数的选择范围相对较宽。选取了若干典型算例对于混合格式的性能进行了研究,并对不同格式的计算精度和效率进行了对比。数值结果表明,与WENO和MUSCL相比,本文提出的混合方法在较广的范围内均表现出良好的特性。     

基于超声速边界层转捩的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阶以上的高精度格式.      

面向混合网格高精度阻力预测的熵修正方法

Roe通量差分分裂格式具有耗散小、分辨率高等特点,在亚跨超声速流场模拟中得到广泛应用,但空间离散时需要对Roe平均矩阵的特征值进行熵修正,使用熵修正会增大耗散,影响阻力的求解精度。本文针对非结构混合网格中Roe格式熵修正的特点,通过改进传统的Roe格式Harten-Yee熵修正方法,提出了一种可提高非结构混合网格黏性计算精度的Harten-Yee熵修正改进方法。利用改进后的方法,完成了DLR-F4翼身组合体算例的计算和对比分析,改进后的熵修正方法残差收敛特性与原始Harten-Yee熵修正一致,计算精度提高,计算结果和无熵修正时基本一致,说明改进后的熵修正方法既保留了使用熵修正带来的程序鲁棒性等优点,同时把熵修正对阻力预测精度的影响降到了最低,验证了方法的有效性。采用改进后方法对第3届AIAA阻力会议的计算模型DLR-F6翼身组合体进行了详细的模拟,分析了网格收敛性和雷诺数的影响。结果表明：改进的Harten-Yee熵修正,更加适用于非结构混合网格的黏性计算,计算精准度达到国际同类CFD软件水平,进一步验证了改进方法的可靠性。