计算流体力学基原理及其发展现状之我见.docVIP

计算流体力学基本原理及其发展现状之我见 一、计算流体力学基本原理 计算流体力学或计算流体动力学，英文Computational Fluid Dynamics，简称CFD，是用电子计算机和离散化的数值方法对流体力学问题进行数值模拟和分析的一个分支。 计算流体力学是目前国际上一个强有力的研究领域，是进行传热、传质、动量传递及燃烧、多相流和化学反应研究的核心和重要技术，广泛应用于航天设计、汽车设计、生物医学工业、化工处理工业、涡轮机设计、半导体设计、HAVCR 等诸多工程领域，板翅式换热器设计是CFD 技术应用的重要领域之一。 为了说明计算流体力学主要方法，需先了解流体力学运动的基本方程的性质和分类。流体力学的基本方程是在19 世纪上半叶由C.-L.-M.-H.纳维和G.G.斯托克斯等人建立的，称为纳维－斯托克斯方程，简称N-S方程。 1、非定常可压缩Navier－Stokes方程 不计质量力的情况下，在直角坐标系中，守恒型N－S方程可以写为下列向量形式： ， 其中 ， ， 。 如果忽略N－S方程中的粘性和热传导，得到的简化方程为Euler方程： 。 方程、称为向量守恒型方程。其重要特点是：连续、动量和能量方程被写为统一形式。其中，均为列向量，是方程的解向量，称为守恒变量；称为通量（flux），具体说为无粘通量，为粘性通量。 所谓守恒型方程，是空间导数项为散度的形式的方程。（3）式所示的向量型守恒方程，实际上仍然是散度形式。显然，（3）式的另一种等价形式为： ， 其中 ， 或 ， 通量张量，为直角坐标系三个坐标轴方向的单位基矢量。把式在任意固定的控制体上积分，并利用Gauss公式，有 。 这就是守恒积分型方程。可见，守恒的微分、积分型方程之间有直接的联系。式是我们以后将要讲到的有限体积方法的出发方程，而、或是则是有限差分方法的出发方程。 2、流体力学方程的简化形式 根据具体流动状态，N－S方程可以进行各种简化。简化的形式及其适用条件是理论流体力学的重要研究内容之一。这里我们对于各种简化方程作一归纳，见下图： 图1.N－S方程的简化形式 3、曲线坐标系中的基本方程 当求解域的形状比较复杂时，计算流体力学方法通常在曲线坐标系中实施。因此，有必要得到曲线坐标系中流体力学基本方程的形式。在曲线坐标中，矢量可以采用在直角坐标中的分量形式，也可以采用协变或逆变分量，基本方程也将因此呈现出不同的形式。最简单，应用也最普遍的形式是矢量分量为直角坐标系中的分量。下面，我们讨论这种情况下的流体力学基本方程。 直角坐标到曲线坐标的变换及其逆变换关系为： （1）、导数的变换 对于一阶偏导数，根据链式求导法则，有 。 同理可得 。 对于二阶偏导数，有 。 同理可得 ， 。 把导数的变换关系代入微分方程，就可以得到微分方程在计算平面中的形式。以直角坐标系中的Laplace方程 为例，把上述二阶导数的变换关系代入上述Laplace方程，得 。 （2）、度量系数及其计算方法 在导数的坐标变换公式中涉及到下列坐标变换系数：。这些系数称为坐标变换公式对应的度量系数（metrics）。我们看到，为了求解计算平面中的偏微分方程，如式，必须确定度量系数（有时还包括等）的离散值。那么，这些度量系数如何计算呢？由于一般情况下，我们只知道坐标变换关系、的离散表达式，度量系数一般也要通过有限差分方法近似计算。但是，直接构造的差分近似是不容易的。以为例，根据偏导数的意义，为保持不变时随的变化，如图2所示，网格点处的的计算公式应为： 。 由于一般不是网格点，因此是未知的，只能通过插值方法确定。 另一方面，我们可以定义逆变换式的度量系数。在贴体坐标系中，这些度量系数的有限差分离散非常简单。如果采用中心差分离散，有 。 这就提示我们，如果能够找到和之间的关系，我们就可以得到等的计算方法。下面，我们推导二者之间的关系。 由变换关系式，有 ， 写成矩阵形式