计算流体力学 全套讲义.DOC

第一章 绪论 第一节 计算流体力学：概念与意义 计算流体力学概述 任何流体运动的规律都是由以下3个基本定律为基础的：1）质量守恒定律；2）牛顿第二定律（力＝质量×加速度），或者与之等价的动量定理；3）能量守恒定律。这些基本定律可由积分或者微分形式的数学方程（组）来描述。把这些方程中的积分或者（偏）微分用离散的代数形式代替，使得积分或微分形式的方程变为代数方程（组）；然后，通过电子计算机求解这些代数方程，从而得到流场在离散的时间/空间点上的数值解。这样的学科称为计算流体（动）力学（Computational Fluid Dynamics，以下简称CFD）。CFD有时也称流场的数值模拟，数值计算，或数值仿真。 在流体力学基本方程中的微分和积分项中包括时间/空间变量以及物理变量。要把这些积分或者微分项用离散的代数形式代替，必须把时空变量和物理变量离散化。空间变量的离散对应着把求解域划分为一系列的格子，称为单元体或控制体（mesh，cell，control volume）。格子边界对应的曲线称为网格（grid），网格的交叉点称为网格点（grid point）。对于微分型方程，离散的物理变量经常定义在网格点上。某一个网格点上的微分运算可以近似表示为这个网格点和相邻的几个网格点上物理量和网格点坐标的代数关系（这时的数值方法称为有限差分方法）。对于积分型方程，离散物理量可以定义在单元体的中心、边或者顶点上。单元体上的积分运算通常表示为单元体的几何参数、物理变量以及相邻单元体中物理变量的代数关系（这时的数值方法称为有限体积方法和有限元方法）。所谓数值解就是在这些离散点或控制体中流动物理变量的某种分布，他们对应着的流体力学方程的用数值表示的近似解。由此可见，CFD得到的不是传统意义上的解析解，而是大量的离散数据。这些数据对应着流体力学基本方程的近似的数值解。对于给定的问题，CFD研究的目的在于通过对这些数据的分析，得到问题的定量描述。在这一点上，CFD与实验研究有类似之处。另一方面，CFD直接处理的是描述流动的数学模型：微分或积分形式的方程组及其边界条件。在这一点上，CFD与理论流体力学又是相同的。 CFD可以应用于所有与流体运动相关的领域。无论在那个领域中，为了获得问题满意的答案，CFD的研究通常应该遵循下面的步骤： 第一，问题的界定和流动区域的几何描述。应明确要解决的问题中流场的几何形状，流动条件和对于数值模拟的要求。几何形状通常来源于对于已知流动区域的测量。如果处于设计阶段，流场的几何形状可能不是完全确定的，在这种情况下，必须知道对于流场的几何形状有哪些限制条件，并根据这些限制条件或其他初步设计手段确定流场的假定形状，然后根据模拟的结果对几何形状进行不断调整，在多次模拟的过程中，逐步确定最终几何形状。流动条件可以包括流动的雷诺数、马赫数、边界处的速度、压力等等。对于数值模拟的要求包括：数值模拟的精度和所花费的时间，所感兴趣的流动参数等。 第二，选择主控方程和边界条件。主控方程指在数值计算过程中要求解的方程。在问题确定后，必须选择流动的主控方程和边界条件。一般认为，在牛顿流体范围内，所有的重要流动现象都可以用Navier－Stokes方程来描述。但是，为了提高计算的效率，有时可以选择经过简化的数学模型（如果这种简化仍能保留流动的物理本质，满足对于数值模拟的要求的话）。简化模型包括势流方程，Euler方程，边界层方程，薄层近似的N－S方程等等。根据问题的特点，可以考虑定常或非定常、可压或不可压的流动模型。边界条件可以有固体壁面条件，来流、出流条件，周期条件，对称条件等。边界条件通常依赖于主控方程，如在固体壁面，Euler方程要求采用不可渗透条件，而N－S方程要求无滑移条件。如果有必要，我们还需要采用一些附加的物理模型，最典型的例子就是湍流模型。虽然N－S方程可以描述湍流流动，但是直接采用原始的N－S方程计算湍流流动（称为直接数值模拟）要求网格点的数量非常多，因而计算量非常大，这是目前的计算机所不能承受的。所以人们通常采用经过Reynolds平均的N－S方程，为了封闭这个方程就必须采用某种湍流模式。 第三，确定网格划分策略和数值方法。在CFD中，网格划分可以有各种不同的策略，如结构网格、非结构网格、组合网格、重叠网格等。网格可以是静止的，也可以是运动的（动网格），还可能根据数值解动态调整（自适应网格）。CFD中的数值方法有有限差分、有限体积、有限元、谱方法等。数值方法和网格划分策略是相互关联的。例如，如果采用有限差分方法，通常要选用结构化网格；而有限体积方法和有限元方法则可以适应于结构和非结构网格。根据网格划分策略和数值方法，最终应该形成数值求解基本方程和边界条件的计算机程序或软件。这些程序，可以是针对某一问题自行编制的，也可