清华大学 任玉新---高等计算流体力学

主要内容任玉新1.Basics2.Methods for compressible flows1) The mathematical properties of Euler equations2) Shock wave and entropy conditions3) Riemann problem and the Godunov scheme4) Approximate Riemann solvers: HLL solver and Roe solver5) TVD scheme6) ENO/WENO scheme7) The compact scheme13.Methods for incompressible flows1) The staggered and the colocated grids2) The MAC method3) The SIMPLE method4) The projection method5) Other methods6) Solution of N-S equations on the nonstaggered gridReferences[1] E. F. Toro, Riemann solvers and numerical methods for fluid dynamics, Springer, 1997 (First edition)[2] J.D. Anderson, Computational fluid dynamics: basics with applications, Springer (清华大学出版社影印版)2[3] Barth and Deconinck (eds.) High order method for computational physics, Lecture Notes in Computational Science and Engineering, 9. Springer, 1999[4] Ferziger and Peric, Computational method for fluid dynamics, Springer, 1996[5] T. J. Chung, Computational fluid dynamics, Cambridge University Press, 2002[6] J. W. Thomas, Numerical partial differential equations: conservation laws and elliptic equations. Texts in applied mathematics 33, Springer, 1999[7] 吴子牛，计算流体力学基本原理，科学出版社， 2002.[8] Sherrie L. Krist, Robert T. Biedron, Christopher L. Rumsey，CFL3D User＇s Manual, The NASA Langley Research Center,Hampton, VA[9] S. K. Lele, J. Comput. Phys. 103, 16 (1992)[10] S. Pirozzoli, J. Comput. Phys. 178 (2002)[11]Yu-Xin Ren, Miao＇er Liu, Hanxin Zhang, J. Comput. Phys. 192 (2003)3FTP: 166.111.37.201Usr:cfdPasswd:cfd2005Email：ryx@高等计算流体力学讲义(1)第一章计算流体力学基本原理第1节流体力学基本方程一、非定常可压缩Navier－Stokes方程45不计品质力的情况下，在直角坐标系中，守恒型N －S 方程可以写为下列向量形式：()()()0v v v t x y z∂∂-∂-∂-+++=∂∂∂∂U F F G G H H ， (1) 其中u v w E ρρρρρ⎛⎫ ⎪ ⎪ ⎪= ⎪ ⎪ ⎪⎝⎭U 2()u u p uv uw E p u ρρρρρ⎛⎫ ⎪+ ⎪ ⎪= ⎪ ⎪ ⎪+⎝⎭F 2()v vu v p vw E p v ρρρρρ⎛⎫ ⎪ ⎪ ⎪=+ ⎪ ⎪ ⎪+⎝⎭G 2()w uw vw w p E p w ρρρρρ⎛⎫ ⎪ ⎪⎪= ⎪+ ⎪ ⎪+⎝⎭H ，0xx xyv xzxx xy xz T u v w kx ττττττ⎛⎫⎪ ⎪⎪=⎪ ⎪⎪∂+++ ⎪∂⎝⎭F 0xy yy v yzxy yy yz T u v w k y ττττττ⎛⎫⎪ ⎪ ⎪= ⎪ ⎪ ⎪∂+++ ⎪∂⎝⎭G ，60xz zyv zzxz zy zz T u v w k z ττττττ⎛⎫⎪ ⎪ ⎪= ⎪ ⎪ ⎪∂+++ ⎪∂⎝⎭H 。

如果忽略N －S 方程中的粘性和热传导，得到的简化方程为Euler 方程：0t x y z∂∂∂∂+++=∂∂∂∂U F G H 。

方程、称为向量守恒型方程。

其重要特点是：连续、动量和能量方程被写为统一形式。

其中，,,,,,,v v v U F G H F G H 均為列向量，U 是方程的解向量，稱為守恆變數；,,,,,v v v F G H F G H 稱為通量（flux ），具體說,,F G H 為無粘通量，,,v v v F G H 為粘性通量。

所谓守恒型方程，是空间导数项为散度的形式的方程。

，式所示的向量型守恒方程，实际上仍然是散度形式。

显然，，式的另一种等价形式为：0t ∂+∇=∂UE ， (2) 其中()()()v v v =-+-+-E F F i G G j H H k，或=++E Fi Gj Hk，通量張量，,,i j k 為直角坐標系三個坐標軸方向的單位基向量。

把式在任意固定的控制体上积分，并利用Gauss 公式，有70Sd dS t Ω∂Ω+=∂⎰⎰⎰⎰⎰U E n 。

这就是守恒积分型方程。

可见，守恒的微分、积分型方程之间有直接的联系。

式是我们以后将要讲到的有限体积方法的出发方程，而、或是则是有限差分方法的出发方程。

二、流体力学方程的简化形式根据具体流动状态，N －S 方程可以进行各种简化。

简化的形式及其适用条件是理论流体力学的重要研究内容之一。

这里我们对于各种简化方程作一归纳，见下图：89图1.N －S 方程的简化形式三、 曲线坐标系中的基本方程当求解域的形状比较复杂时，计算流体力学方法通常在曲线坐标系中实施。

因此，有必要得到曲线坐标系中流体力学基本方程的形式。

在曲线座标中，向量可以采用在直角坐标中的分量形式，也可以采用协变或逆变分量，基本方程也将因此呈现出不同的形式。

最简单，应用也最普遍的形式是向量分量为直角坐标系中的分量。

下面，我们讨论这种情况下的流体力学基本方程。

直角坐标到曲线座标的变换及其逆变换关系为：(,)(,)x y x y ξξηη== (3)(,)(,)x x y y ξηξη== (4)1、导数的变换对于一阶偏导数，根据链式求导法则，有x x x φφφξηξη∂∂∂=+∂∂∂ 。

同理可得10y y y φφφξηξη∂∂∂=+∂∂∂。

对于二阶偏导数，有222222222222222()()()[][]()2()x x xx xx x x xx xx x x x x x x xx xx x x x x x x x x φφφξηξηφφφφξηξηξηξηφφφφφφξηξξηηξηξηξξηξηηφφφφφξηξξηηξηξξηη∂∂∂∂=+∂∂∂∂∂∂∂∂∂∂=+++∂∂∂∂∂∂∂∂∂∂∂∂=+++++∂∂∂∂∂∂∂∂∂∂∂∂∂=++++∂∂∂∂∂∂。

同理可得222222222()2()yy yy y y y y y φφφφφφξηξξηηξηξξηη∂∂∂∂∂∂=++++∂∂∂∂∂∂∂， 222222()xy xy x y x y y x x y x y φφφφφφξηξξξηξηηηξηξξηη∂∂∂∂∂∂=+++++∂∂∂∂∂∂∂∂。

把导数的变换关系代入微分方程，就可以得到微分方程在计算平面中的形式。

以直角坐标系中的Laplace 方程22220x yφφ∂∂+=∂∂ 为例，把上述二阶导数的变换关系代入上述Laplace 方程，得11222222222[()()]2[][()()]()()0x y x x y y x y xx yy xx yy φφφξξξηξηηηξξηηφφξξηηξη∂∂∂+++++∂∂∂∂∂∂++++=∂∂。

2、度量系数及其计算方法在導數的座標變換公式中涉及到下列座標變換係數：,,,x y x y ξξηη。

这些系数称为座标变换公式对应的度量系数（metrics ）。

我們看到，為了求解計算平面中的偏微分方程，如错误！未找到引用源。

式，必頇確定度量係數（有時還包括,,,,,xx xy yy xx xy yy ξξξηηη等）的離散值。

那么，这些度量系数如何计算呢？由于一般情况下，我们只知道座标变换关系、的离散运算式，度量系数一般也要通过有限差分方法近似计算。

但是，直接構造,,,x y x y ξξηη的差分近似是不容易的。

以x ξ為例，根據偏導數的意義，x ξ為y 保持不變時ξ隨x 的變化，如圖2所示，網格點P 處的x ξ的計算公式應為：12)Q Px P Q Px x ξξξ-=-。

由於Q 一般不是網格點，因此,Q Q x ξ是未知的，只能通過插值方法確定。

另一方面，我們可以定義逆變換(4)式的度量係數,,,x x y y ξηξη。

在贴体坐标系中，这些度量系数的有限差分离散非常简单。

如果采用中心差分离散，有131,1,,,1,1,1,1,,,1,1,)2)2)2)2i j i j i j i j i j i j i j i j i j i j i j i j x x x x x x y y y y y y ξηξηξηξη+-+-+-+--=∆-=∆-=∆-=∆ 。

這就提示我們，如果能夠找到,,,x y x y ξξηη和,,,x x y y ξηξη之間的關係，我們就可以得到,,,x y x y ξξηη等的計算方法。

下面，我们推导二者之间的关系。

由变换关系式，有x y x y d dx dy d dx dyξξξηηη=+=+，写成矩阵形式x y x y d dx d dy ξξξηηη⎡⎤⎡⎤⎡⎤=⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦ 。

由逆变换式，有x x dx d y y dy d ξηξηξη⎡⎤⎡⎤⎡⎤=⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦。

、式中的矩阵称为正变换和逆变换的Jacobi 矩阵。