随着CFD技术的不断发展以及计算机和计算能力的迅速提高，CFD已经成为飞行器研究和设计中有力的工具，已经广泛运用于飞机的气动设计和优化。升力和阻力是影响飞行器飞行的两个最重要的因素。几乎所有的气动分析都是为了尽最大可能的在给定阻力的情况下最大化升力，或者在给定升力的情况下最小化阻力。对飞机外形的升力和阻力的计算成为对其气动分析的基础。精确地计算飞行器的气动性能，特别是飞机的阻力，是多年来空气动力学家所追求的目标。由于飞行器的阻力比升力要低一个数量级，所以无论是实验还是CFD计算，都很难获得阻力的准确值，而阻力的微小差异会严重影响飞行器的性能和有效载荷，阻力的准确计算成为气动力计算的一个最重要目标。

对于阻力的计算，无论是实验方法还是数值模拟方法都面临着巨大挑战。实验方法不能完全的描述真实流场，比如对于雷诺数，风洞结果还不能达到真实情况的量级，还有风洞璧面对流场的影响，对于其修正方法还依赖于经验，因此不可避免存在误差问题。从数值模拟的角度考虑，即使采用最精确的方法也很难保证压力分布和切向力分布的计算结果完全准确，因此由此积分得到的阻力很难表征其真实结果。从国际第一届和第二届阻力预测会议上可以看出，阻力计算结果的分散度较大，精确的阻力预测仍然是一项挑战。对典型的跨声速飞机而言，1个阻力单位（0.0001）也就是约占总阻力0.4%的预测水平还没达到。通过数值仿真来精确预测气动阻力，已经成为CFD研究的一个重要目标。虽然提高阻力计算精确度可以通过提高计算格式的精度，但是二阶以上的高阶精度计算格式作为工程应用是无法接受的，因此如何在已有的二阶精度格式下去努力提高阻力计算的准确度是我们面临的重要任务。

近些年来随着电子器材及机械加工小型化技术方面的进步，人们越来越关心微型飞行器MAVs(Micro Air Vehicles)的发展。MAVs体积比以前所研究的飞行器小，工作于包括低雷诺数环境在内的独特的气动领域。由于其特殊性能，MAVs能够应用于如监视、生物化学药剂探测、战场勘测、交通监测及城市信息收集等方面。而微型旋翼飞行器在这些方面相比固定翼飞行器有着明显的优势，特别是当飞行器要求保持不动(悬停)或在相当紧凑的环境中作机动飞行时旋翼飞行器的优势更为突出。涉及微型旋翼飞行器发展的主要问题之一便是旋翼的空气动力计算。它对于飞行安全和旋翼性能的提高均起着重要作用。

对于微型旋翼飞行器的气动特性计算人们大多采用叶素理论方法简单地把二维翼型计算结果积分得到三维旋翼的升力或扭矩等参数，该方法适用于工程上的估算，但存在着计算精度不高的问题。而本文通过直接求解三维Eider方程可以得到更为详细的流场细节，如速度场分布，压力场分布等，能为微型旋翼的设计提供更准确的设计依据。

微型旋翼一般工作在低速、低雷诺数条件下，此时空气接近于不可压。选取适用于大型直升机旋翼的可压缩流动计算程序求解低速流动容易出现收敛性困难的问题。通常求解不可压流动方程可采取的方法有:(1)MAC方法，(2)SIMPLE方法，(3)拟压缩性方法； (4)预处理矩阵方法。蔡伟明通过对不可压连续性方程进行拟压缩修正得到新的控制方程，该方程选用与桨叶固连的旋转坐标系，方程的求解采用有限体积空间离散和龙格库塔时间推进格式。求解出的物理变量值是旋转坐标系下的绝对物理量。由于悬停流场的准定常特性，故整个计算过程只对一个桨叶进行，其他桨叶的影响由周期性条件引入。远场边界条件采用Riemann不变量方法来处理。运用当地时间步长、压强阻尼和隐式残值光顺等技术使计算加速收敛到定常状态。 2100433B

物体与空气作相对运动时作用在物体上的力，简称气动力。它由两个分布力系组成：一是沿物体表面面元法线方向的法向分布力系，另一是在表面面元切平面上的切向分布力系。空气动力通常就是指这两个力系的合力。以飞行器(如飞机)为倒，为便于对飞行器的运动规律进行分析，通常取一个原点位于飞行器重心的气流坐标系，将空气动力分解为三个方向上的分量。设坐标系的x轴平行于气流方向且正向与气流方向相反，y轴在飞行器对称面内与x轴垂直且正向指向飞行器上方，z轴垂直于xy平面，指向右翼，则合力在x、y、z三个轴上的分量分别称为阻力、举力和侧向力。若空气动力作用点与飞行器重心不重合，则飞行器还受到一个合力矩的作用，它在x、y、z三个轴上的分量分别称为滚转力矩、偏航力矩和俯仰力矩。飞行器所受的空气动力与它的飞行速度、高度和飞行姿态有关。空气动力的分布和大小是飞行器结构和强度设计的依据，而且关系到飞行器的飞行性能、操纵性能和稳定性。空气动力学的一个主要任务就是确定飞行器的空气动力。确定空气动力需要知道空气的性质和运动规律。相应于低速流动、亚声速流动、跨声速流动、超声速流动、高超声速流动、稀薄气体流动和高温气体流动等不同情况，空气动力的分析有不同的理论和实验方法。

气动力计算近场法

数值计算外形气动力最常用的方法是对外形的表面压强和应力的积分，这也称为近场法。它的准确度取决于积分点的数量，表面曲率的变化和积分方法的精度。在某种意义下，近场法是计算阻力最自然的方式，但进一步地分析发现，即使流场计算结果比较理想，这种方法也不容易得到精度较高的阻力，其原因主要在于数值误差对阻力的计算有重要的影响，在近场法中，数量级较小的压差阻力的积分是由数量级比它大的吸力和阻力之差确定的，在这种情况下，很容易产生数值误差，除非能够给出足够精度的压力分布，但实际情况往往并非如此，对摩擦阻力的计算也会出现类似的困难，在转棙点附近，摩擦力系数的变化非常显著，计算出来的摩擦阻力也会出现很大的数值误差。

另一方面，采用近场法计算出来的阻力是最后的积分值，不能给出阻力各分量的值。

气动力计算远场法

基于控制体的途径去测量阻力，最简单方法就是通过测量模型中尾流上平行于来流方向的动量损失，然而这种方法要求测量整个尾流区域，同时又要注意风洞璧面的影响，这其中是困难重重。Betz 改进了这种方法，考虑到风洞璧面的影响，修正了积分公式，同时缩小了模型后的积分区域。但是 Betz 并没有解决由涡产生这部分阻力，忽视了由于有限翼展后面存在自由涡而产生的阻力，而这个阻力是有限翼展机翼产生升力所必须付出的阻力代价。机翼后方自由涡面上的流体微团旋转所需的能量，必须由机翼提供一个附加的推力来克服诱导阻力才能维持有升力的飞行。后来 Maskell 进一步修正了 Betz 方法中的一些问题，重新推导出一个积分公式，在飞行器模型尾流的一个修正的区域积分，可以同时计算出型阻和波阻。自此以后在 Betz 和Maskell 研究阻力方法基础上，对此进行各种改进研究，在此基础上产生传统的远场积分法和尾迹面积分法。

传统远场法计算阻力基于积分形式的动量定理，将总阻力表示为下游无穷远处假设平面上的积分，从中分解出各种阻力成分，再对每一个阻力分量得出一阶精度（相对于阻力大小）的计算公式。这种方法的优点是并不需要知道飞行器外形的详细几何信息，也可以很容易地将阻力分解（型阻和诱导阻力），有利于设计人员根据给出的结果给出好的减阻方法，而且这种方法产生的数值误差是各阻力分量的高阶小量，因而计算出来的总阻力精度要比采用近场法计算结果的精度高。但是由于远场处的网格比较稀疏，再加上流场解的耗散，传统的远场法并不能满足计算精度要求。

尾迹面法是对传统远场法的改进，假设将控制体的进口面和旁侧曲面移至无限远，传统远场积分的升力和阻力就可以写成在出口面（尾迹面）上的积分，通过在尾迹面上对相对流动参数进行积分得到升力和阻力，一般尾迹面是选择在接近飞行器后缘不远处，保证积分区域含有涡量，在同样网格精度下精确度更高，流场解的耗散也很小，同时也具有传统远场积分法的优点。

《计算空气动力学并行编程基础》由国防工业出版社出版。

根据锅炉本体各部受热面的结构计算锅炉通风过程的全压降,以便选择合理的通风装置。

第1章基本方程形式

1.1引言

1.1.1计算空气动力学的作用

1.1.2计算空气动力学的特点

1.1.3计算空气动力学的步骤

1.2热力学基础

1.2.1气体状态方程

1.2.2比热、内能和焓

1.2.3黏性系数

1.2.4传热系数

1.3直角坐标系下微分形式N—S方程

1.3.1笔续介质假设

1.3.2方程通用形式

1.3.3质量方程

1.3.4动量方程

1.3.5能量方程

1.3.6N—S方程组

1.4无量纲化

1.4.1特征量的选取方法

1.4.2方程组无量纲化因子

1.4.3方程组无量纲化

1.4.4特征速度不同引起的表达式差异

1.5下线坐标系下微分形式N—S方程

1.5.1度量系数

1.5.2坐标变换

1.5.3有限差分法

1.6积分形式N—S方程

1.6.1直角坐标系下积分形式N—S方程

1.6.2曲线坐标系下积分形式N—S方程

1.6.3有限体积法

1.6.4有限体积法度量系数

1.6.5有限体积法单元体积

1.7有限差分法与有限体积法

1.7.1有限差分法与有限体积法的异同

1.7.2有限差分一有限体积混合方法

第2章时间离散格式

2.1显式Runge—Kutta法

2.2隐式方法残差项线性化

2.2.1对流通量Jacobian矩阵与分裂

2.2.2黏性通量Jacobian矩阵与谱半径

2.2.3有限体积法七对角矩阵

2.2.4有限差分法十三对角矩阵

2.3隐式线性化定常流动求解方法

2.3.1当地时间步长

2.3.2定常流方程组形式

2.3.3点松弛SGS解法

2.3.4线松弛SGS解法

2.3.5LU—SGS解法

2.3.6LU—ADI解法

2.3.7Jacobi迭代法

2.4隐式线性化非定常流动求解方法

2.4.1单时间步法

2.4.2双时间步法

2.5初始条件

2.5.1给定原始变量

2.5.2给定飞行高度、马赫数、攻角与侧滑角

2.5.3给定雷诺数、马赫数、温度、攻角与侧滑角

2.6收敛判据

2.6.1定常流动迭代收敛判据

2.6.2非定常流动亚迭代收敛判据

第3章空间离散格式

3.1曲线坐标系下的导数计算

3.2对流项半点重构

3.2.1二阶MUSCL重构

3.2.2五阶显式WCNS重构

3.3对流通量导数

3.3.1矢通量分裂格式（FVS）

3.3.2通量差分裂格式（FDS）

3.3.3混合格式（AUSM）

3.4黏性通量导数

3.4.1完全N—S方程

3.4.2薄层近似方程

3.5再论通量Jacobian矩阵

第4章边界条件处理

4.1虚网格概念

4.2特征分析

4.2.1边界信息传播方向

4.2.2特征线相容关系式

4.3远场边界条件

4.3.1远场特征变量边界条件

4.3.2远场Riemann边界条件

4.4壁面边界条件

4.4.1无黏壁面

4.4.2黏性壁面

4.5对称边界条件

4.6内边界条件

第5章加速收敛技术

5.1多重网格方法

5.1.1基本循环FAS格式

5.1.2多重网格循环策略

5.1.3传递算子

5.2预条件Newton—Krylov算法

5.2.1非线性函数

5.2.2非精确Newton法

5.2.3Krylov子空间方法

5.2.4预条件JFNK算法

第6章串流模型方程

6.1雷诺时均控制方程

6.1.1雷诺时均方法

6.1.2时均方程无量纲化

6.2零方程湍流模型

6.2.1Baldwin—Lomax湍流模型

6.2.2Degani修正的B—L湍流模型

6.3一方程湍流模型

6.3.1Baldwin—Barth湍流模型

6.3.2Spalart—Allmaras湍流模型

6.4两方程湍流模型

6.4.1Wilcox k—ω模型

6.4.2Menter’s k—ω模型

6.4.3EASM Gatski—Speziale k—ω模型

6.4.4EASM Girimaji k—ε模型

6.4.5Abid k—ε模型

6.4.6EASM Gatski—Speziale k—ε模型

6.5初始条件与边界条件

6.6应用湍流模型的注意事项

6.7两方程湍流模型通用形式

第7章高温气体效应

7.1气体模型

7.1.1质量比焓与内能

7.1.2质量比熵

7.1.3混合气体自由焓

7.2高温气体化学平衡流

7.2.1Gibbs最小自由焓方法

7.2.2能量守恒方程计算温度

7.3高温气体热化学非平衡流

7.3.1高温气体非平衡方程组

7.3.2对流项Jacobian矩阵

7.3.3组分生成源项及其Jacobian矩阵

7.3.4振动源项及其Jacobian矩阵

第8章前后处理方法

8.1数据文件格式

8.1.1网格与流场文件格式PLOT3D

8.1.2边界文件格式Generic.inp

8.1.3Tecplot格式的可视化文件

8.1.4CGNS格式转换

8.2网格奇点与搜索

8.2.1奇点重构的问题描述

8.2.2重构算法及时间复杂性分析

8.2.3数值实验结果与讨论

8.3网格二次剖分技术

8.3.1网格块分组问题

8.3.2二次剖分方法

8.3.3剖分法的选择策略

8.4气动特性计算

8.4.1气动力

8.4.2气动力矩

8.4.3升/阻力与压心

8.4.4气动热

第9章并行编程基础

9.1并行计算与并行编程

9.1.1并行计算相关概念

9.1.2并行计算机系统结构与分类

9.1.3并行编程模型与实现技术

9.1.4并行编程的常用模式

9.2并行程序设计

9.2.1并行程序流程设计

9.2.2MPI进程并行

9.2.3OpenMP线程并行

9.3并行程序性能优化技术

9.3.1提升并行程序性能的技术

9.3.2数值实验结果与讨论

9.4并行软件框架PETSc

9.4.1PETSc软件体系结构

9.4.2PETSc软件功能

9.4.3基于PETSc编程

9.4.4多DMDA创建实例

参考文献

主要符号表2100433B