学位论文 > 优秀研究生学位论文题录展示

基于非结构网格方法的重叠网格算法研究

作　者: 田书玲
导　师: 伍贻兆
学　校: 南京航空航天大学
专　业: 流体力学
关键词: 计算流体力学非定常流场非结构重叠网格并行算法多体相对运动外挂物投放直升机旋翼
分类号: TP301.6
类　型: 博士论文
年　份: 2008年
下　载: 630次
引　用: 8次
阅　读: 论文下载

内容摘要

多体间具有相对运动的非定常绕流问题是航空工程中一类常见的问题。在这类问题中,物体间存在着较强的气动干扰,而这种气动干扰常常是非线性和非定常的,并对飞行器的气动性能有较大的影响,甚至会对飞行安全构成威胁。近年来,随着数值方法的不断发展以及计算机软硬件的突飞猛进,采用计算流体力学方法模拟多体间具有相对运动的非定常绕流问题已成为可能。Nakahashi提出的基于非结构的重叠网格思想结合了非结构网格方法和结构重叠网格方法的优点,适合应用于多体间具有大幅相对运动的复杂外形非定常绕流问题模拟研究,具有较大的实际工程应用前景。但是,与非结构网格技术和结构重叠网格技术相比,基于非结构的重叠网格方法还处于起步发展阶段,其中的一些关键技术,如网格间插值边界定义方法、适用于非结构重叠网格的网格间插值方法、高效宿主单元搜索算法、向任意网格类型的推广、基于非结构重叠网格方法的分布式并行计算算法以及适用于多体相对运动非定常绕流流场模拟的动态非结构重叠网格算法等方面都需要进行深入研究。因此,本文的主要工作是开展基于非结构网格方法的重叠网格算法研究,其目的是发展出适合于多体相对运动非定常绕流问题模拟的高效计算算法,并开发出相应的非定常计算软件。本文首先讨论了所采用的流动控制方程以及数值求解方法。控制方程采用非定常可压缩Euler/N-S方程,空间离散采用有效体积方法,通量计算采用Jameson中心格式,时间离散采用双时间步长方法,在伪时间上采用多步Runge-Kutta迭代,湍流模拟中采用Spalart-Allmaras一方程湍流模型计算湍流粘性系数,其物面距计算采用了一种基于阵面推进的高效计算方法。其次,开展了非结构重叠网格的网格间边界定义方法研究,并在Nakahashi的网格间边界定义方法的基础上,提出了一种能够处理包含任意多块子网格、任意单元类型子网格的重叠网格系统的优化网格间边界定义算法;为提高非结构重叠网格系统的网格间边界定义效率,在宿主单元搜索中,发展了一种适合于非结构网格数据结构的基于阵面推进与ADT方法相结合的相邻单元搜索算法,大大提高了大型非结构重叠网格系统的网格间边界定义的速度;为适应混合非结构重叠网格的网格间信息交换,在进行非结构重叠网格的网格间插值方法研究中提出了适合于混合类型单元网格的网格间插值方法。然后,为把非结构重叠网格方法推广应用到多体间具有相对的非定常绕流问题模拟中,开展了基于动态非结构重叠网格的非定常Euler/N-S方程求解算法研究,提出了一种适合于多体相对运动模拟的动态非结构重叠网格算法,并且开展数值实验研究,表明该算法在非定常问题模拟中不但具有较高的计算效率和较强的鲁棒性,而且不需任何人工干预,具有很高的自动化程度。最后,开展了基于动态非结构重叠网格方法的非定常Euler/N-S方程分布式并行算法研究。在采用区域分裂算法并行求解Euler/N-S方程时,提出了一种子网格独立分区的非结构重叠网格分区策略,并在子网格分区时根据网格节点的性质(活动节点、非活动节点或者插值点)进行加权分区,使得每个进程上的计算负载达到平衡。进程间通信通过调用MPI标准库函数实现。为测试所发展的基于动态非结构重叠网格的并行计算算法的正确性和并行效率,开展了一系列的数值实验研究,结果表明该并行算法不但能够得到与串行一样的正确结果,而且具有较高的并行计算效率。

全文目录

摘要  4-6
ABSTRACT  6-17
第一章绪论  17-29
  1.1 研究背景  17-19
  1.2 动态网格技术发展  19-24
    1.2.1 动态重叠网格方法  19-22
    1.2.2 动态非结构网格方法  22-23
    1.2.3 动态自适应直角坐标网格方法  23-24
  1.3 非结构重叠网格方法发展及其研究现状  24-26
  1.4 本文主要工作及创新点  26-29
    1.4.1 本文主要工作  26-28
    1.4.2 本文的特色和创新点  28-29
第二章 Euler/N-S 方程数值求解方法  29-59
  2.1 引言  29-30
  2.2 流动控制方程  30-32
  2.3 流动控制的空间离散方法  32-36
    2.3.1 基于非结构网格的有限体积方法  32-34
    2.3.2 Jameson 中心格式  34-35
    2.3.3 粘性通量项的计算方法  35-36
  2.4 流动控制方程的时间离散方法  36-40
    2.4.1 定常方程时间离散方法  36-38
    2.4.2 非定常方程时间离散方法  38-40
  2.5 边界条件  40-45
    2.5.1 物面边界条件  40-42
    2.5.2 远场边界条件  42-45
    2.5.3 对称边界条件  45
  2.6 湍流模型  45-49
    2.6.1 Spalart-Allmaras 湍流模型  45-47
    2.6.2 非结构网格物面距计算方法  47-49
  2.7 预处理 Euler/N-S 方程数值求解方法  49-55
    2.7.1 预处理可压控制方程  50-53
    2.7.2 预处理Euler/N-S 方程数值求解方法  53-55
  2.8 数值算例  55-57
    2.8.1 NACA0012 翼型绕流模拟  55
    2.8.2 DLRF6 翼身组合体绕流模拟  55-57
  2.9 本章小结  57-59
第三章非结构重叠网格方法  59-91
  3.1 引言  59-60
  3.2 非结构重叠网格的网格间边界定义方法  60-65
    3.2.1 网格节点分类方法  60-61
    3.2.2 网格间边界定义  61-63
    3.2.3 网格间边界的优化  63-65
  3.3 宿主单元搜索方法  65-68
  3.4 混合非结构网格重叠方法  68-71
  3.5 任意多块非结构网格重叠方法  71-74
  3.6 网格间信息交换方法  74-80
    3.6.1 网格间插值边界重构  74-78
    3.6.2 非结构重叠网格插值方法  78-80
  3.7 数值算例  80-89
    3.7.1 RAE2822 翼型定常粘性绕流  80-83
    3.7.2 某三段翼型定常粘性绕流  83-84
    3.7.3 M6 机翼定常绕流  84-87
    3.7.4 三球体定常粘性绕流  87-89
  3.8 本章小结  89-91
第四章基于动态非结构重叠网格法的多体相对运动非定常绕流模拟方法研究  91-115
  4.1 引言  91-92
  4.2 刚体六自由度运动方程及其求解方法  92-96
    4.2.1 刚体六自由度运动方程  92-93
    4.2.2 刚体六自由度运动方程数值求解方法  93-94
    4.2.3 刚体位置和姿态的更新方法  94-96
  4.3 基于动态非结构重叠网格法的多体运动模拟算法  96-98
  4.4 三维外挂物投放模拟研究  98-105
  4.5 直升机旋翼流场模拟研究  105-113
    4.5.1 Caradonna 旋翼悬停流场模拟  105-107
    4.5.2 Caradonna 旋翼前飞流场模拟  107-109
    4.5.3 Robin 直升机机身-旋翼组合前飞流场模拟  109-113
  4.6 本章小结  113-115
第五章动态非结构重叠网格法并行计算技术  115-146
  5.1 引言  115-116
  5.2 并行环境及并行计算基本概念介绍  116-118
    5.2.1 并行环境  116-117
    5.2.2 并行计算基本概念  117-118
  5.3 基于非结构重叠网格的 Euler /N-S 并行方法  118-130
    5.3.1 非结构重叠网格的负载平衡分区策略  118-124
    5.3.2 并行非结构重叠网格计算的数据交换方法  124-128
    5.3.3 动态非结构重叠网格计算的并行算法  128-130
  5.4 验证算例  130-140
    5.4.1 M6 机翼定常绕流并行模拟  130-135
    5.4.2 三维外挂物投放过程并行数值模拟  135-140
  5.5 并行动态非结构重叠网格法的应用研究  140-144
  5.6 本章小结  144-146
第六章全文总结与展望  146-149
  6.1 全文总结  146-147
  6.2 进一步工作展望  147-149
参考文献  149-158
致谢  158-159
在学期间的研究成果及发表的学术论文  159

基于非结构网格方法的重叠网格算法研究

内容摘要

全文目录

相似论文