基于混合网格的三维Navier-Stokes方程并行算法

提出了一种基于混合网格的三维Navier-Stokes方程的并行计算方法。Navier-Stokes的求解采用了基于面的有限体积方法,该方法适用于任何网格类型。采用一方程Spalart-Allmaras模型来计算紊流黏性。并行计算采用区域分裂的方法,利用METIS网格分区系统实现了各节点的加载平衡。节点间的数据交换通过调用MPI库函数来实现,采用非阻断通讯的方式来减少数据交换时间。充分利用FORTRAN90的动态存储特性来减少对内存的需求。最后,通过对绕DLR-F6外形(翼身组合体+挂架+发动机短舱)黏性流动的数值模拟,验证了该并行程序的准确性,高性能并行计算以及处理复杂几何外形的能力。     

航天飞行器脉动压力数值计算方法综述

航天飞行器以高马赫数飞行时,由于边界层的复杂流动将产生强烈的脉动压力环境。由此会引发仪器设备强烈的振动与噪声环境,从而大大降低了系统的可靠性,甚至导致飞行失败。长期以来,国内外在飞行器脉动压力预示方面开展了大量的理论分析、数值计算与试验研究,取得了一定的成果。随着计算流体力学和脉动压力计算方法的日趋成熟,数值方法正日益成为脉动压力环境预测的主要工具。针对飞行器脉动压力环境的基本特征,对目前常用的脉动压力数值计算方法和现有研究成果进行了较为全面的介绍,对新一代航天飞行器的研制提供了一定的参考依据。     

考虑半径补偿的三坐标测量数据精确配准方法

叶片模型的配准定位是叶片形状检测分析的一个核心环节,针对涡轮叶片的三坐标测量数据在基于迭代最近点的配准定位中与CAD模型表面点的配准关系不对应的问题,提出一种考虑半径补偿的精确配准方法。在求取配准点集的过程当中考虑到半径补偿,使得测量点集与待配准点集之间的关系相对应,进而获取更高的配准精度。进行基于仿真数据的配准试验及对比,给出了基于实测数据的应用与分析,验证了该方法的精确性与实用性。结果表明,考虑半径补偿的三坐标测量数据精确配准方法是可行的。     

Apram模型上双曲型方程初边值问题数值方法的并行性研究

针对Ａｐｒａｍ模型［１］，研究了双曲型方程初边值问题差分格式的并行实现方法。文中不仅讨论了差分格式在Ａｐｒａｍ模型上的各种并行处理技术，而且对并行算法的可扩放性进行了讨论，最后设计了适合于Ａｐｒａｍ模型结构特征的并行计算程序，对二维Ｅｕｌｅｒ方程正规激波反射问题的计算实践，结果令人满意     

三维弹性边界元分析中面力不连续和1/r积分奇异性问题的处理

本文将文[2]中的方法推广应用于解决三维弹性面积分中1/r的奇异性问题,从数值计算角度,较详细地讨论了消除1/r奇异性的原理。为了取得较高的数值精度及较少的运算时间,本文对高斯点数目的选择原则也作了分析与讨论,提出了一个实用的经验公式。针对三维弹性问题中常见的边界面力不连续问题,提出了一种简单的处理方法,即直接从边界离散化后的边界元方程出发,按力不连续点所处的积分单元分别进行处理。两个典型算例的数值试验结果表明,文中所用的方法是行之有效的。     

矩形喷管内外喷流流场的数值模拟

基于 Ｆｒａｖｅ 平均的 Ｎ Ｓ方程和 Ｂ／ Ｌ 湍流模型,采用了 Ｊａｍ ｅｓｏｎ 格式和 Ｍ Ｕ Ｓ Ｃ Ｌ 格式,对外流在亚声、超声和零 Ｍ ａｃｈ 数三种状态时矩形喷管的喷流流场进行了数值模拟,探讨了内外喷流及喷管内的流动特征,并与实验数据进行了比较,计算结果与实验结果吻合较好,验证了计算所使用的 Ｎ Ａ Ｐ Ａ 软件的可靠性,从而为用计算流体力学的方法模拟具有两股来流的三维湍流流场提供了一种有效的手段。     

三维Euler方程的自适应八叉树结构直角网格算法

给出了三维自适应八叉树结构直角网格的生成方法,采用了以几何外形为基础的网格加细方法,并解决了网格生成过程中的若干难点,使网格生成更简便迅速,产生的网格更适合于Ｅｕｌｅｒ方程计算。在计算中采用以中心差分为基础的Ｊａｍｅｓｏｎ的有限体积法和以面为基础的通量计算方法,减少了工作量。在物面边界处采用通量分解的方法实现边界条件,简单易行。本文对前机身、单个及多个外挂等问题进行了数值实验。结果表明,计算与实验符合较好。     

基于Nash均衡与进化计算的协调航迹规划

在对策论框架下研究多无人飞行器的协调航迹规划问题,提出了一种基于Nash平衡和进化计算的航迹规划算法.协调航迹规划是一个多目标合作优化问题,也是一个非合作对策问题.飞行器编队的分布结构提供了一种隐性偏好信息,有助于找到飞行器编队的最优协调航迹和最优飞行代价分配.本文算法使用特定的航迹个体表示方法和进化算子,按照层次方法处理各种航迹约束.不同飞行器间的协同和竞争关系定义在航迹评价函数里.仿真结果验证了基于博弈论航迹规划算法的可行性和优越性.     

低马赫数条件下几何折转角对平面扩压叶栅弯叶片流场性能的影响

将低马赫数大转角弯叶片数值计算结果与实验结果进行了比较,结果表明计算与实验结果吻合得较好,在计算中弯叶片降低能量损失系数为11%,在实验中弯叶片降低能量损失系数为13%.为了深入讨论压气机中弯叶片的适用条件,在0°迎角下,只改变几何折转角,共做了4套平面叶栅的数值分析,其结果表明大转角弯叶片的效果较为明显,比如在59.5°折转角叶栅中,最佳弯角为30°的正弯叶片降低能量损失系数的17.9%;而在39.5°折转角叶栅中,最佳弯角15°正弯叶片只减小了损失系数的0.7%,当进一步增加弯角时,损失反而增加.可以证明在低马赫数条件下叶型折转角的大小是在压气机中是否使用弯叶片的一个重要参考因素.     

基于N-S方程串并行计算的机翼优化设计

对典型的大展弦比和中等展弦比机翼在0.7~0.9马赫数区间进行了串行和并行计算的数值校验,证实了以N-S (Navier-Stokes)方程为主控方程的串行和并行流场求解器的正确性,并讨论了并行效率和加速比.结合Powell算法,讨论了在确定的机翼平面形状和翼型的条件下,以升阻比最大为目标的三维机翼截面翼型最大厚度与扭角的优化设计.算例结果表明,厚度的非线性分布和负的扭角会改善机翼流场的流动状态,使机翼的升阻比得到提高,优化设计方法是可行的.