基于多块结构网格的并行计算及负载平衡研究

从并行计算流体力学程序的稳定性和效率两大问题入手,针对多块结构网格的通用数据传输方法和基于遗传优化算法的负载平衡方法,并在已有串行多块结构网格程序基础上发展了相应的并行程序.该并行程序以物理区域分割为基础,采用MPI实现消息传递,适用于各种不同的并行机体系结构,具有很好的可移植性.大量数值实验证明,本文发展的并行程序具有良好的稳定性和并行效率,可以进一步应用于大规模实际工程计算.     

复杂外形网格生成技术及应用

以某型双立尾战斗机为例,给出了实际复杂型号飞行器结构网格生成的一种分块方案及实现过程,用一种矢性三次多项式插值的方法生成相邻块公共交界面的网格,用Hilgenstock方法修正源项,求解椭圆型方程生成块内空间网格,使网格线与边界面夹角、与边界面间距控制为期望的值.在所生成的双立尾战斗机分块结构网格中求解Euler方程,获得了合理的全机气动力系数、翼面压力分布和绕流图谱结果.     

非结构三角形网格的一个ENO型有限体积方法

基于双曲型守恒律方程,对非结构三角形网格给出了一种ENO(Essentially Nonoscillatory Scheme)型有限体积格式,方法的主要思想是先对每一个三角形单元构造一个加权的二次插值多项式,而在计算交界面的流通量时采用了两点高斯积分公式以保证格式的整体精度,时间离散采用三阶TVD Runge-Kutta方法.最后给出了该格式收敛的数值阶,并对前台阶问题进行了计算.      

用于复杂流场数值模拟的重叠／嵌套网格方法

介绍了一类用于复杂流场数值模拟的重叠／嵌套网格方法。它是一类多块结构化网格,可用于有或无多重网格加速收敛技术的有限体积、有限差分等计算。其基本思想是把复杂流场的计算区域化分成相对简单的子区域,每个子区域上的网格可以独立生成;网格界面之间的流场信息交换通过三线性插值来实现,通过修改流场解算方法来排除洞中的点。     

真实飞行器网格生成策略及在航空CFD中的应用

从应用于型号气动设计的实用软件研制的需求背景角度,探讨真实飞行器网格生成策略,以几种典型的军机、民机及导弹实用外形为例,研究生成贴体及与边界正交高质量计算网格,并应用于复杂流场计算。经全机亚、跨、超音速、大迎角及发动机进排气复杂流场测试,空气动力和力矩系数以及物面压力分布与实验数据吻合良好,证明本文论述的网格生成方法是实用、有效的。     

利用NURB作曲线和曲面的插值

给出了非均匀有理三次Ｂ－样条插值的方法，利用曲线（面）方程的矩阵表达式导出了反求顶点问题的方程组，使方程组的系数矩阵呈三对角型，易于求解。同时，分别推出几种情况下的端点（边界）条件。利用该结果可以使得曲线（面）易于调整，运算速度快，便于处理。     

发动机进气流动三维瞬态数值模拟研究

采用任意拉格朗日-欧拉法(ALE法)和所提出的动态网格生成新方法,实现了车用发动机进气流动过程的三维瞬态可压缩湍流数值模拟,湍流模型采用经压缩修正的κ-ε双方程模型.计算结果揭示了发动机进气过程的三维瞬态流动特性,为发动机进气系统的优化设计提供了理论依据.     

并行计算中一种非结构网格分割方法

将递归谱对剖分方法应用于流体力学并行计算中的非结构网格分割,以解决负载平衡和最小切割问题。为使用这种方法,计算了网格伴随图的离散Laplacian矩阵的第二特征矢,然后从该特征矢的分量引入网格的对剖分。特征矢计算中应用了Rayleigh商迭代,并进行了一些修正以使收敛强烈地偏向于第二特征矢及考虑逆迭代步中线性方程组的迭代求解。最后,通过非结构自适应网格上Euler方程分区计算的数值结果验证了所发展的网格分割方法。     

基于非线性有限元降阶方法的网格加筋筒壳结构屈曲承载特性研究

非线性有限元数值计算方法是结构屈曲性能评估的主要手段,然而该方法需要基于精细网格迭代求解,降低结构非线性计算规模是其关键技术之一。采用一种结构非线性有限元降阶方法,对均匀轴向压载荷作用下的网格加筋筒壳结构进行了细致的承载稳定性特性分析。该方法将Koiter初始后屈曲理论与Newton弧长法相结合,通过多次使用Koiter摄动展开来跟踪结构非线性平衡路径。计算比较了特征值屈曲载荷和计及几何非线性效应的非线性屈曲载荷值;采用侧向小扰动载荷模拟结构初始几何形状缺陷,分析了缺陷幅度对临界屈曲载荷的影响规律,在验证了非线性有限元降阶方法分析精度的同时,细致考核了该方法在结构非线性屈曲分析以及缺陷敏感性分析中的计算效率优势。     

螺旋桨飞机俯仰力矩特性改进方法

螺旋桨滑流对飞机各部件的气动干扰造成飞机的俯仰力矩特性恶化,危害飞行安全。开展螺旋桨滑流对螺旋桨飞机俯仰力矩特性的影响机理研究和螺旋桨滑流作用下螺旋桨飞机俯仰力矩特性改进方法研究十分重要。采用动态重叠多块结构网格,通过求解非定常雷诺平均可压缩Navier-Stokes方程,数值模拟了某螺旋桨飞机不同拉力系数下降落构型的绕流流场。结果显示螺旋桨滑流对机翼、平尾的气动干扰是导致飞机俯仰力矩特性恶化的主要原因,改进螺旋桨飞机俯仰力矩特性的关键是改进其平尾的升力特性。以拉力系数等于0.4时的降落构型为优化对象,开展了俯仰力矩特性的改进方法研究,包括降低平尾高度、改变平尾上反角、抬高螺旋桨轴线等方法。研究发现在不增大机翼对平尾整体下洗强度的前提下,减小平尾和螺旋桨轴线的垂向距离,可以明显地改善螺旋桨飞机的俯仰力矩特性。