面向化学非平衡流的CFD并行计算技术和大规模并行计算平台研究
发布时间:2020-05-01 21:49
【摘要】:吸气式高超声速技术是指飞行马赫数大于5、以吸气式发动机及其组合发动机为动力、在大气层和跨大气层中实现高超声速远程飞行的飞行器技术。吸气式高超声速技术的核心是超燃冲压发动机技术和机体/推进一体化飞行器技术。吸气式高超声速飞行包含许多复杂的物理、化学现象,大幅度增加了地面试验的费用、周期和难度。随着计算技术的发展,计算流体力学(CFD)可以在相对较短的时间内完成流动参数的分析研究,提供流场的详细特征,弥补风洞试验的局限性,因此在超燃冲压发动机的设计和性能研究方面得到了广泛的应用。与此同时,超燃冲压发动机流场模拟对计算流体力学也提出了更高的要求,超燃冲压发动机流场模拟涉及到各种组分方程、化学非平衡效应、湍流方程、输运系数的计算等,其计算量相当大,一般都需要在高性能并行计算机系统上进行大规模并行计算。 实际的CFD并行计算必须考虑串并行结果是否一致、大规模并行时的并行效率以及并行可扩展性等问题,在功能扩充方面要求CFD计算支持多种物理化学模型、差分格式和隐式求解方法并能适用于不同的网格体系。本文主要从以上几个方面对超燃冲压发动机流场模拟的CFD并行计算进行了研究,针对每个方面均提出了创新的见解,开发了主要应用于超燃冲压发动机流场模拟的面向化学非平衡流的三维大规模并行软件平台AHL3D。 在大规模并行技术方面,本文提出了一种静态负载平衡策略,采用一维优先的规则分块算法和基于自动重分块的不规则分块算法分别处理大规模和中小规模并行的问题。为保证复杂计算区域问题的通信效率和正确性,将通信分为三步:面通信、物理边界面棱通信、通信边界面棱通信,保证了通信时点对的对应关系,避免了复杂拓扑结构时某些特殊网格点的重复通信和通信的奇异性。 将流水线并行技术与CFD相结合,提出了一种基于流水关系有向图的多维流水线并行实现算法,根据各求解子域的邻接几何关系和变量依赖关系分析,自动判断流水线并行是否可行,自动确定流水方向和流水维数,实现了复杂计算区域问题不同并行规模串并行结果的完全一致。 实现了对接网格、错位拼接网格和重叠网格的多块并行计算。提出了一种“迷路算法”应用于重叠网格“挖洞”,解决了其它算法难以处理凹型包络面的问题,算法的健壮性较好,实现简单,,而且算法计算复杂度低,计算效率高。在多层次嵌套重叠情况下的通信时序控制方面,提出了重叠关系有向图避免通信等待和重复插值。 在大规模线性代数方程组的迭代求解方法上,提供了三种Krylov子空间迭代
【图文】:
只要重叠包络面封闭,不论该包络面外形如何复杂,迷路算法都可以准确地得到包络面内外点。因三维算例不易直观显示,所以采用一个较复杂的二维外形作为示例。如图4.7所示,该正方形网格为主网格,有两个外形较为复杂的从网格,图4.8的阴影部分即为得到的洞内点。 1015图4.7二维算例重叠网格示意图图4.8二维算例对应“挖洞”结果4.7.2迷路法效率测试在“挖洞”算法的处理效率方面,我们在主频为3G的 hitelXeon计算节点上实测了不同算例迷路算法和整个“挖洞”过程的墙上时间,具体结果如下。可以看出,该“挖洞”算法对于大网格数据量、多层次嵌套重叠和复杂组合体问题仍能得到较高的效率。表4.1部分典型算例“挖洞”效率示意典型算例主网格数(万)从网格数(万)迷路时间(秒)挖洞时间(秒)三维凹型圆柱三维单机翼飞行器整机模型13,1518.60104.3521.63 5.9278,01 0.0030 0.0048 0.123 0.65 0.5812.054.7.3找重方法测试我们在主频为3G的 hitelXcon计算节点上进行了正确性验证和时间效率测试,正确性验证以最终的找重结果文件为准。因kd树搜索时可以判断一个点落在两个第73页
该正方形网格为主网格,有两个外形较为复杂的从网格,图4.8的阴影部分即为得到的洞内点。 1015图4.7二维算例重叠网格示意图图4.8二维算例对应“挖洞”结果4.7.2迷路法效率测试在“挖洞”算法的处理效率方面,我们在主频为3G的 hitelXeon计算节点上实测了不同算例迷路算法和整个“挖洞”过程的墙上时间,具体结果如下。可以看出,该“挖洞”算法对于大网格数据量、多层次嵌套重叠和复杂组合体问题仍能得到较高的效率。表4.1部分典型算例“挖洞”效率示意典型算例主网格数(万)从网格数(万)迷路时间(秒)挖洞时间(秒)三维凹型圆柱三维单机翼飞行器整机模型13,1518.60104.3521.63 5.9278,01 0.0030 0.0048 0.123 0.65 0.5812.054.7.3找重方法测试我们在主频为3G的 hitelXcon计算节点上进行了正确性验证和时间效率测试,正确性验证以最终的找重结果文件为准。因kd树搜索时可以判断一个点落在两个第73页
【学位授予单位】:解放军信息工程大学
【学位级别】:博士
【学位授予年份】:2006
【分类号】:TP338.6
本文编号:2647096
【图文】:
只要重叠包络面封闭,不论该包络面外形如何复杂,迷路算法都可以准确地得到包络面内外点。因三维算例不易直观显示,所以采用一个较复杂的二维外形作为示例。如图4.7所示,该正方形网格为主网格,有两个外形较为复杂的从网格,图4.8的阴影部分即为得到的洞内点。 1015图4.7二维算例重叠网格示意图图4.8二维算例对应“挖洞”结果4.7.2迷路法效率测试在“挖洞”算法的处理效率方面,我们在主频为3G的 hitelXeon计算节点上实测了不同算例迷路算法和整个“挖洞”过程的墙上时间,具体结果如下。可以看出,该“挖洞”算法对于大网格数据量、多层次嵌套重叠和复杂组合体问题仍能得到较高的效率。表4.1部分典型算例“挖洞”效率示意典型算例主网格数(万)从网格数(万)迷路时间(秒)挖洞时间(秒)三维凹型圆柱三维单机翼飞行器整机模型13,1518.60104.3521.63 5.9278,01 0.0030 0.0048 0.123 0.65 0.5812.054.7.3找重方法测试我们在主频为3G的 hitelXcon计算节点上进行了正确性验证和时间效率测试,正确性验证以最终的找重结果文件为准。因kd树搜索时可以判断一个点落在两个第73页
该正方形网格为主网格,有两个外形较为复杂的从网格,图4.8的阴影部分即为得到的洞内点。 1015图4.7二维算例重叠网格示意图图4.8二维算例对应“挖洞”结果4.7.2迷路法效率测试在“挖洞”算法的处理效率方面,我们在主频为3G的 hitelXeon计算节点上实测了不同算例迷路算法和整个“挖洞”过程的墙上时间,具体结果如下。可以看出,该“挖洞”算法对于大网格数据量、多层次嵌套重叠和复杂组合体问题仍能得到较高的效率。表4.1部分典型算例“挖洞”效率示意典型算例主网格数(万)从网格数(万)迷路时间(秒)挖洞时间(秒)三维凹型圆柱三维单机翼飞行器整机模型13,1518.60104.3521.63 5.9278,01 0.0030 0.0048 0.123 0.65 0.5812.054.7.3找重方法测试我们在主频为3G的 hitelXcon计算节点上进行了正确性验证和时间效率测试,正确性验证以最终的找重结果文件为准。因kd树搜索时可以判断一个点落在两个第73页
【学位授予单位】:解放军信息工程大学
【学位级别】:博士
【学位授予年份】:2006
【分类号】:TP338.6
【引证文献】
相关期刊论文 前2条
1 刘鑫;陆林生;陈德训;;非结构网格并行计算预处理方法研究[J];计算机科学;2012年03期
2 刘鑫;陆林生;;拼接网格通量守恒插值算法研究[J];计算机应用与软件;2012年02期
相关博士学位论文 前1条
1 李义进;燃气轮机全通流数值模拟研究[D];哈尔滨工程大学;2012年
本文编号:2647096
本文链接:https://www.wllwen.com/kejilunwen/jisuanjikexuelunwen/2647096.html