二维结构/非结构混合网格的生成及其流场模拟

采用分区求解方法在二维结构/非结构混合网格中对粘性绕流作了数值模拟。利用层面推进法生成绕物体的O形结构网格,利用推进阵面法在结构网格外生成非结构网格。在结构网格区内,采用B-L湍流模型,求解雷诺平均N-S方程。在非结构网格区内求解Euler方程。两区的解通过在交界面上满足通量守恒条件进行耦合。数值算例表明了方法的有效性和实用性。     

共轭梯度方法在求解不可压N-S方程中的应用

采用非结构网格上的SIMPLE(Semi-implicit Method for Pressure-linked Equations)算法、k-ε湍流模型求解了三维不可压Navier-stokes方程,对低速不可压粘性流场进行数值模拟,在求解压强方程时采用共轭梯度方法,求解其它变量方程时,采用预处理的BICG(biconjugate gradients)算法.对NACA0012翼型绕流流场和飞艇绕流流场进行数值模拟并对结果进行分析,取得较好的结果.     

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

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

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

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

用N－S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二级时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型￣［１］。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

用N—S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二维时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

三维非结构混合网格生成与N-S方程求解

 提出了一套较为通用的, 完全自动化的非结构混合网格生成方法, 在物面粘性作用区, 运用一种改进的推进层方法生成三棱柱形和金字塔形网格; 在其他流动区域采用阵面推进方法生成四面体网格。采用一种改进精度的格心有限体积法和Spalart-Allmaras 一方程湍流模型对三维N-S 方程进行了求解, 并以M6 机翼、DLR-F4 翼身组合体和某运输机外形的粘性流场作为数值算例, 验证了上述网格生成和流场求解方法的正确性和实用性。     

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

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

飞机增升装置气动力特性计算方法研究

本文给出了飞机增升装置气动力特笥的数值计算方法和计算结果，为了计算增升装置的缝隙效应和剪刀口效应，求解了粘性三维Navier-Stokes方程。为了正确 模拟大舵偏和阻力板后的大面积分离，采用了反映涡流特笥的新类型的两方程湍流模型。     

NASA TrapWing高升力标模数值模拟研究

针对NASA TrapWing高升力全展襟翼构型，采用计算流体力学（ CFD）方法进行三维复杂流场仿真模拟，考察网格尺度和湍流模型对高升力模型气动特性的影响。采用“超立方体”概念，生成绕TrapWing模型的不同网格密度高质量多块结构网格，通过求解雷诺平均Navier-Stokes方程，研究网格尺度对高升力模型气动特性的影响。在此基础上选取中等规模计算网格，考察Spalart Allmaras和Menter k-ωSST湍流模型对高升力全展构型流场模拟能力，分析湍流模型对气动特性的影响。研究结果表明：网格策略具有较好的网格收敛性；湍流模型对翼稍附近位置上翼面压力系数的预测稍有影响，SA湍流模型预测的压力系数较SST更接近实验值。确认研究工作为大型飞机增升装置数值模拟提供了一定的参考。