NS方程激波计算的摄动有限差分方法

摄动有限差分(PFD)方法从一阶迎风差分格式出发,将差分系数展开为网格步长的幂级数,通过提高修正微分方程的逼近精度来获得更高精度的差分格式。由于格式基于一阶迎风格式,因此具有迎风效应、网格节点少等特点。本文首先通过对Burgers方程的摄动差分格式的推导,将摄动有限差分格式引入时间相关法的计算,并构造了守恒形式的摄动有限差分格式,然后推广到一维Navier-Stokes方程组的计算。数值比较研究表明:本文构造的NS方程摄动有限差分格式具有比一阶迎风较高的精度和分辨率,而且保持了一阶迎风格式的无振荡性质。     

涵道尾桨气动特性及翼型CFD分析

以Jameson/TVD混合格式对低速时的OAF和NACA63a312翼型的升阻特性进行了计算。该方法的特点是采用有限体积法进行离散,并用龙格-库塔时间推进法进行求解。通过将计算的结果与NASA公布的试验曲线进行对比,验证了该混合格式用于旋翼气动力计算的可行性。为了进一步研究旋翼的拉力和功率,本文建立了以微圆环动量理论和叶素理论为基础的计算模型,并用该模型计算了涵道尾桨在不同条件下的拉力和功率情况,计算的结果具有一定的代表性。      

高阶精度格式WCNS-E-5在亚跨声速流动中的应用研究

采用高阶精度非线性紧致加权格式WCNS-E-5和Baldwin-Lomax模型,求解雷诺平均Navier-Stokes方程,开展了典型翼型与机翼的湍流流动数值模拟研究.对方程中粘性项采用的四阶精度差分近似以及网格导数求解与边界格式的四阶精度,保证了高精度算法的实现.计算结果表明:本文算法能够准确地模拟这些翼型与机翼的亚跨声速流场,得到与实验测量十分吻合的壁面压力分布,计算结果对网格的依赖性小.     

超声速底部喷流干扰流场数值模拟

数值模拟了不同马赫数,不同喷流压比下的轴对称超声速底部喷流干扰流场,采用LU隐式算法进行数值求解并引入了Baldwin-Lomax代数湍流模型.采用分区网格将弹身与底部区域合为一个整体进行计算,得到了清晰的流场结构和弹体表面及底部的压力分布,并与试验结果进行了比较.数值模拟结果表明超声速底部喷流干扰流场结构复杂,有、无喷流时底部流场有很大不同, 喷流对底压分布有明显影响,进而对轴向力系数影响显著.     

对我国高等教育扩招的正负效应的思考

我国高等教育扩大招生规模,并进入到马丁.特罗教授理论中的“大众化“阶段,这是中国高等教育的进步。一方面,扩招带来的机遇和动力显而易见,这样可以促进经济发展,缩小我国同发达国家差距,使更多的人可以接受到良好的教育;另一方面,扩招带来的负效应不容忽视,这主要表现在高等教育质量下滑,学费高昂和就业压力增大等问题上。为此提出一些建设性意见,希望从根本上减少由于扩招后的负效应给我国高等教育带来的种种冲击。     

使用高阶逆风通量差分裂格式的悬停旋翼流场数值模拟

为减少Jameson二阶中心差分有限体积法导致的旋翼尾迹的数值耗散,将三阶逆风格式(MUSCL)与通量差分裂方法相结合,建立了一个三维雷诺平均N-S方程数值模拟悬停状态旋翼流场的方法。为了充分考虑旋翼尾迹对流场的影响,采用周期性边界条件和由动量理论导出的远场边界条件。为进一步减少尾迹数值耗散和便于添加上述边界条件,采用了嵌套网格方法。然后,进行了算例计算,给出了桨叶表面的压强分布,与可得到的试验数据及二阶中心差分方法的计算结果进行了对比,并针对几种不同桨尖形状的旋翼悬停流场进行计算,数值结果显示:后掠桨尖可减弱超临界流动。此外,还计算和分析了旋翼下方不同位置上的涡量分布。计算结果表明,本文方法能够有效地减少旋翼尾迹的数值耗散。      

高分辨率格式在非定常无粘可压缩流动中的应用研究

本文发展了一种计算非定常无粘可压缩流动的高分辨率格式。通量计算采用AUSMPW+Rie-mann求解器,其中使用了基于单调性的压力加权函数f,并与Roe Riemann求解器进行比较。为提高精度,使用五阶WENO格式进行左右状态插值,且时间推进采用三阶TVD Runge-Kutta方法。用该方法对移动接触间断问题、Sod激波管问题、二维激波反射问题和双马赫反射问题进行了数值计算,数值结果表明基于五阶WENO插值的AUSMPW+格式有很高的激波及接触间断分辨能力,并比基于五阶WENO插值的Roe格式有更高的计算效率。     

详细反应机理的航空发动机二维燃烧流场计算

对某航空发动机主燃烧室简化模型进行了基于正庚烷(n-C₇H₁₆)详细反应机理的燃烧流场数值模拟.在二维贴体坐标系下,计算过程采用了k-ε双方程模型来预估湍流特性,用EDC(eddy dissipation con-cept)湍流燃烧模型预估反应速率,用PSR(perfectly stirred reactor)模型的方法处理复杂化学反应项,得出了基于详细反应机理的温度场、组分浓度场,验证了机理文件中反应的某些组分与特定自由基之间的关系,对比了总包反应EBU-SRK(eddy break up-simplified reaction kinetics)的温度场,结果基本相符.      

基于非结构动态网格的翼型积冰过程数值模拟

讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法.外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格-库塔方法.采用四阶龙格-库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率.使用非结构动态网格技术对积冰后的翼型网格进行快速调整.在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性.     

预处理AUSMDV格式在微型扑翼绕流数值模拟中的应用

由于微型扑翼的低速低雷诺特性,对其粘性绕流的数值模拟存在较大困难。传统的中心空间离散格式不具有格式耗散而需要加入带有自由度的人工耗散,给计算带来不确定因素;此外,在低速低雷诺条件下,RANS方程组的条件数太大,导致计算收敛缓慢甚至不收敛,需要引入预处理方法加速收敛。本文使用具有低耗散特性的AUSMDV空间离散格式结合预处理方法对扑翼的非定常粘性绕流进行了数值模拟,非定常推进采用了含双时间迭代的LU-SGS隐式时间推进,湍流模型采用了BL模型。计算结果表明,加入预处理能够克服AUSMDV格式难以收敛的缺点,并且与中心格式相比,基于预处理的AUSMDV格式能够更准确地模拟扑翼的非定常粘性绕流,通过将本文计算结果与文献结果对比验证了本文方法的有效性。在此基础上研究了扑翼相关参数对气动特性的影响,对进一步了解扑翼气动力产生机理,指导扑翼的气动设计具有参考价值。