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

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

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

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

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

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

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

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

带运动顶盖的封闭方腔内湍流流场的计算

本文运用κ-ε模型求解了带运动顶盖的封闭方腔内的二维湍流流场,其流动雷诺数为10000,对乘方定律格式(PLDS)和二次上风差分格式(QUICK)进行了比较,采用了PISO求解方法、截面连续性修正并对k、ε方程的源项作特殊处理,使计算的收敛性和稳定性明显提高。计算出的主涡、二次涡位置及流场速度分布值与实验吻合良好。计算所用网格为12×12,在Honeywell DPS 8/49计算机上运行的CPU时间对PLDS与QUICK分别为1.9分和1.8分。     

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

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

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

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

结合FLUENT软件的叶型粘性逆命题气动设计方法

为了有效解决粘性条件下的反问题,发展了一种满足粘性假设的逆命题气动设计方法,并将其通过UDF与FLUENT软件结合,构建了逆命题设计平台。该方法给定叶型的目标压力分布,采用目标和当前压力分布的差计算出网格变形量修正叶型,直到获得满足目标压力分布的型线。采用本平台对一组抛物线叶栅进行计算验证,之后完成了涡轮叶栅的重新设计。抛物线和涡轮叶栅的逆命题计算结果的压力分布均与目标吻合良好,证明了本平台的有效性。逆命题计算后压力残差下降约85%~95%,计算时间约为分析计算的4~7倍。     

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

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