非结构网格及混合网格复杂无粘流场并行计算方法研究

在分布式存储网络异构编程PVM环境下,本文将原串行计算程序改造为并行计算程序,提高了计算效率.利用非结构网格,对球体及航天飞机外形绕流流场进行了计算,取得了较好的性能结果.最后将该方法应用于混合网格的计算,求解了多种复杂外形的无粘绕流,如翼身组合体和DDF1C导弹,得到了满意的结果.     

直排波瓣喷管引射器流场计算κ-ε模型的选择

应用N-S方程,对一种实验引射性能较好的直排波瓣喷管引射器内复杂流场进行了数值模拟.采用三维贴体坐标,同位网格,研究了标准κ-ε、Chen κ-ε和RNG κ-ε湍流模型应用于波瓣喷管引射器流场计算的可行性和计算精度.结果表明,上述三种湍流模型均能揭示混合段中流向涡的发展规律.总的看来,Chen κ-ε模型和RNG κ-ε模型计算结果与测量值差别最小,远远好于标准κ-ε模型.     

三维非结构网格的欧拉方程聚合多重网格法

介绍一种基于三维非结构网格的聚合多重网格技术,用于求解Euler方程定常解,加速收敛效果十分明显。求解Euler方程的基本算法采用有限体积中心差分格点格式。聚合多重网格技术按一定的规则将细网格的若干控制体积聚合在一起,组成一个较粗的网格。数值实验表明本文多重网格法十分有效。     

翼－身－尾组合体绕流的Euler方程计算

本文采用重迭网格技术和ＥＪｅｒ方程计算翼－身－尾组合体绕流。对算身与尾。身部分采用各自的Ｈ。Ｏ型网格，Ｅ＊ｅ／方程求解采用１ａｍｅｓｏｎ的有限体积法，即中心差分近似和显式Ｒｕｎｇｅ·Ｋｕｔｔａ时间推进。采用前后区交替迭代使前后两区通过重迭区交换信息。本文用ＮＡＣＡＴＮ４０４１翼－身，尾模型为例，计算的空气动力特性与实验符合较好。     

低速大迎角三角翼旋涡流动的有粘与无粘数值模拟

通过人工压缩性方法求解了三维定常不可压Euler方程和雷诺平均Nayierstokes方程，应用Beam－Warming近似因子分解格式及其对角化形式，采用BaldwinLomax代数湍流模型。计算模型为70°三角翼，分析了有粘与无粘流动旋涡特性的影响，尤其是对Euler方程模拟大迎角分离涡、涡破裂的能力及其局限性提供了数值验证，计算结果与实验结果相一致。     

应用混合网格分析复杂积冰翼型气动性能

通过发展一种点-点对接的结构/非结构混合网格生成方法,避免了复杂积冰外形难以划分结构化网格的问题,并对复杂积冰翼型的气动性能进行了分析计算.在靠近积冰边界的内层采用了非结构网格以拟合复杂的边界,在非结构网格外采用结构化网格以节省存储空间和计算时间.为了检验网格对选用湍流模型的影响,整个流动区域分别采用SST和SA湍流模型,求解了雷诺平均N-S方程.数值计算结果表明SST模型更适于模拟复杂分离流动,积冰对翼型的气动性能造成了严重的影响.     

双欧法与四元数法的应用比较

对解决欧拉方程奇异性的双欧拉和四元数法进行了对比研究，四元数法从理论上讲比较完美，但实际应用中存在较大的累积计算误差，从而影响计算精度；双欧法利用正、反欧拉方程间精华区倒挂关系进行分区交替运算，把精华区扩展到全域，不仅根除了奇异性，而且计算误差小。因此，对于解决欧拉方程奇异性来讲，双欧法要优于四元数法。     

用Weibull模数表征碳纤维的力学性能

碳纤维属于脆性材料,随机分布的缺陷控制其抗拉强度。用Weibull最弱连接理论及经验方程式可阐明和计算抗拉强度与CV值以及与直径的关系。本文提出用拉格朗日乘子多元回归计算的双参数双模态韦氏模型方法,并用直接拟合法计算三参数双模态韦氏模型,得到了很好的结果。     

航天器对接中接近至首次接触阶段的数值仿真

确定接触点问题是航天器对接过程中接近与接触阶段的仿真的关键问题。以内导向瓣异体同构周边式对接机构为例，给出了确定接触点的方法，建立了确定接触点的数学模型、接近与接触阶段两航天器的刚体运动方程，并对接近与接触阶段进行了数值仿真。仿真结果表明：所建立的数学模型和动力学方程是正确的。