1、第一章 绪论1.1 计算流体力学简介（一） 什么是计算流体力学1以计算机作为模拟手段，运用一定的计算技术寻求流体力学各种复杂问题的离散化数值解。 数值解而不是解析解； 计算技术起关键作用； 与计算机的发展紧密相关2计算流体力学、理论流体力学、实验流体力学是流体力学研究工作的三种主要手段既互相独立又相辅相成 理论分析具有普遍性各种影响因素清晰可见、为实验和计算研究提供依据 实验研究仍是研究工作的基石，数值研究的许多方面都密切依赖于实验研究：实验提供数据；计算结果需由实验验证；观察实验现象分析实验数据以建立计算模型等等。 数值模拟是特殊意义下的实验，也称数值实验（二） 计算流体力学研究工作的方向1
2、 与现代计算技术的发展相关联的研究方向：（与计算物理，计算力学发展、图形学、网格技术等）2 与离散数学的理论研究相关连的研究方向；离散化理论、边界条件数值处理的稳定性分析、格式的熵条件等3 在一些相关学科的边缘上寻求新的发展点；4 解决众多相关学科的的科研工作和工程实际提出的与流体力学问题有关的各类复杂的问题机械、航天航空、气象、海洋、石油、环境（包括气动噪音控制） 、建筑、（三） 计算流体力学研究工作的优势、存在的问题和困难1优势：“数值实验”比“物理实验”具有更大的自由度和灵活性，例如“自由”地选取各种参数等“数值实验”可以进行“物理实验”不可能或很难进行的实验；例如：天体内部地温度场数值
3、模拟，可控热核反应地数值模拟“数值实验”的经济效益极为显著，而且将越来越显著；2问题与不足流动机理不明的问题，数值工作无法进行；数值工作自身仍然有许多理论问题有待解决；离散化不仅引起定量的误差，同时也会引起定性的误差，所以数值工作仍然离不开实验的验证；1.2 流体力学微分方程的数学性质当微分方程转化为差分方程并用数值方法求解时，不同类型的微分方程，其数值处理方法各异，其中包括提法的适定性、物理解的性质、差分格式的适用性等；在一些特殊的问题中，甚至通过差分格式的特殊技巧来改变方程的数学性质（一） 一阶拟线性微分方程组的分类对于一阶拟线性微分方程组的向量形式： FxtiUA其中 为 n 阶矩阵nu
4、uU21若： 的特征值为 ,则：A的 根即 0),.21(IAni当 n 个特征值全部为复数时，称方程在 (t , xi)平面上为纯椭圆型；当 n 个特征值全部为互不相等的实数时，称方程在(t , xi)平面上为纯双曲型；而当 n 个特征值全部为实数，但有部分为相等的实数时，称方程在(t , xi)平面上为双曲型；当 n 个特征值全部为零时，称方程在 (t , xi)平面上为纯抛物型；当 n 个特征值部分为复数、部分为实数时，称方程在(t , xi)平面上为双曲椭圆型；二阶拟线性方程组，可以通过降阶法进行类似的分析。(实例见后)（二） 流体力学控制方程数学分类的举例：1 二维定常理想流体 流动
5、的 Euler 方程payvxuayvxupxyuvxyv220)(10)(写成向量形式： 0yxUBA0yxUCpvuUupA01vpvB01 222 22201)(0auvauvvauvauvuBAC1求矩阵 C 的特征值得：24,32,1 4220)()()(auvvvauu如果： 双 曲 椭 圆 型两 个 实 根 ， 两 个 复 根 ，） 四 个 实 根 ， 双 曲 型） 10212Mv2。二维非定常理想流体流动的 Euler 方程yxtUBABACDC110t求 C 的特征值，结论与定常相同：得到在 X-Y 平面的方程性质；求 D 的特征值，得： auauu11,10)()432,
6、为四个实根，即方程在 Y-t 平面为双曲型；所以 Euler 方程可以在时间座标方向推进，而在定常问题中能否推进计算，必须根据流动是否为超音速（M 与 1 的关系）来定。3定常不可压缩 Navier Stokes 方程的数学分类)(1022yvxypyvxu uyvxu 降阶法：令： xvhyugvxfhxvfuyhxfgvfuhyvxypgfu0)(122以下的分析与一阶拟线性方程组的讨论相似，结论为定常 N-S 方程为椭圆型。4非定常不可压缩 Navier Stokes 方程5定常可压缩 Navier Stokes 方程6非定常可压缩 Navier Stokes 方程类7抛物化 N S 方
7、程利用边界层流动的概念，设 X 方向为主流方向，即考虑有： 2yx把流动方向的二阶偏导数略去， （注意与边界层方程不同的是一阶偏导数都将保留！）结论是定常 N-S 方程此时变为抛物型方程。1.3 模型方程以及在计算流体力学中的应用（一） 。模型方程的引入简化对差分格式的性质的讨论及考核必须反映物理问题的最基本的特征，且方便于进行理论分析例如：以流涡函数描述二维流动问题时有方程： 2 22)(yxyvxut 又如： )(122yVxpyVvxtV 模型方程可以提炼为： 一维 Burger 方程。2uxut（二） 。几个典型的模型方程 一维波传播方程： 0xat 一维热传递方程： 2ut 一维对流
8、扩散方程： 2xat Laplace 方程： 022yux Burger 方程： 2xt 无粘 Burger 方程： 0ut其中前 4 个方程为线性方程，可求出解析解，后两个方程为非线性方程，也可以求出解析解。 Burger 方程的解析解：（131） ：粘性系数， 时为无粘方程。2xuut0解： 时，可令未知函数具有如下的形式：0（132）xtxu2),(其中 是待定的二阶可微分函数，将其代入（131）式：,t )(2)(2ttxt xtu)(422xxxu 3222 )(xxx代入（131） ，则得 0xtx不妨设 为满足抛物方程得解 ，即： 0xt（133）则（132）给出了 Burger
9、 方程的解析解的一般形式。若 的初始条件为 ，则由（132）给出的的对应),(txu)()0,(xfu于 的初始条件是：),(t )()(ep),(0xFdfx由（133）给出的 Burger 方程的通解是：deFttx tx4)(2)41),(再代入（132）可得 的解析解。特别指出，粘性 Burger 方程的解),(xu是连续的。无粘 Burger 方程解的间断性： )()0,(xfxut类比于线性的波传播方程，可以设通解为： )(),(utxftxu讨论： tftufuxxx 1)1( 若 （135）0),01()(sstftf即表示在 t-x 座标中的某一个特定的点，其对应的 u 为
10、u(s),使上式为）（ sxt,)(0。此时 且,xu 0(0)()( ssutxft即 ，必定在 点发生解的间断，间断的位置由（135）式),(t)(,sSx确定1.4 计算流体力学的控制方程（一） 。概述 物理问题的解析解通常依赖于一些特殊形式的方程（如采用若干简化假设、选用高阶形式的方程，以避免方程个数过多等） ；而对于计算流体力学的求解，由于数值计算的特点，不惧怕方程的复杂，因此可以基于更一般的方程形式，例如，可以由最基本的物理守恒律导出的普遍形式作为出发方程。 对于求解一个实际问题，往往需要在任意曲线坐标系下来描述，在数值计算中，需要贴体坐标系。因而，要求在任意的非正交曲线坐标系下来
11、描述物理问题是极为平常的情况。 （此点于理论的解析解是极为不同的） 对于理论分析，采用守恒或非守恒变量，守恒方程或非守恒方程，通常没有本质的差别，但在离散的数值计算中，守恒型与非守恒型将可能导致很大的差别，故方程的守恒性是计算流体中，必须特别注意的问题。（二） 。基本守恒律的数学描述1连续方程： 0)(Vt2动量方程 Ft)()(为粘性应力张量*Ip *)()( FIpVt3能量方程：)()( VqEt21Vp)*()()( VTkq )*()()( VTkFVpEt （三） 。直角坐标系下的守恒型方程不计质量力（或质量力有势） ，理想流体、0zGyExFtUtEwvuupwvuFt)(2vp
12、Euvt)(2wpEvuGt)(2若考虑粘性，则： 0111 zyxtUxTkwvuFzxyxy01 yTkwvuEzyxyx1zTkwvuGzzyxzyx01任意非正交曲线坐标系下的守恒形式的方程补充知识：符号规定：直角坐标系： ，单位向量：),(321x3,21ie任意曲线坐标系： 或,),(21关系： )(321xii指标重复，规定为哑指标，进行遍历求和。哑指标一般写成一上一下，以便标识。1.协变基向量：（Covariant base Vectors） 3,21)(lim0irraiiiii 协变基向量 表示的是沿 坐标线的切线方向的单位向量。ii2.逆变基向量：（Contravaria
13、nt base Vectors ）3,21iexajiii逆变基向量 表示的是 的等值面的法方向的单位向量。i consti在正交坐标系中， ，但在任意的非正交坐标系中，协变基不平行于逆变基。ia/iidrdjijias)()(2定义： jijiijgg是二阶对称张量ij微元弧长： iiiidgaidrs指 标 不 累 加 ）此 式 中(同时根据拉梅(Lame)系数的定义: iighds微元面积： 3221)(dad不 累 加kjkjkji ,微元体积： 321321321321 )( dadadV另外，有 wvuwvu321321 321321312321332212132321det dg
14、dg dggaadadVij 向量的散度： 31313132130 )()( )()(limi iikjikji ijikji kjii iiVaAagaAgdgA当 为常向量时，应有：A31310)( 0)(0ikjii kjiaag i ikjAgA梯度： 310)(limi ikji kjiVagd协变基向量与逆变基向量间的关系：)(1)(133kjlnmi lnmililii agagaikja)( likjilkjlil aag )()(1jiija将 式代入散度的表达式 iiiimiii kji gAagAgaAg )(1)(11)(31 )(ii张量的散度 ： )(1 (1( )
15、(1. .3kii ilkliii kjiieTg ageTgaaTgT 利用以上有关张量在任意曲线坐标系中的表示式，并将其代入基本守恒律方程： 标量方程：若令 即为连续方程。0)(1iibgtBt VB向量方程0)(1.ikieTgttA令 此为理想流体（无粘）条件的应力张量。IpV, 0)()( iiii agt 上式是向量表示式，将其写成分量式，理论上可以有许多的投影方式，实际使用中，常将任意的曲线坐标系与直角坐标系结合使用。说明如下：在直角坐标系中， leukwjviu在任意的曲线坐标系中： iaV协变基向量逆变 分量而且： jijii ux,将上式点乘 ：me0)()( jiilii
16、l expgeVgtgmimiuu例如：m=1 上式为： 0)( )3()2()1( xpguVxpguVxpgVtg 有些文献中习惯将曲线坐标系中的 逆变分量 写成： ；由（ ）式:WUV)3()2()1( ,wvuzyxvWzyx0GEFtUtEwvugUUpEwvugFtzyx)(VpEwvugtzyx)(WpEwvugGtzyx)(对于粘性问题：动量方程增加了 一项，*11* ,*kijiii egag jiijkijjilkli xeea ,*所以 N-S 方程： 0)()( * jimjimiimiim xgxpguVgtug 能量方程中的粘性项部分： ljjiji ugV*,1)*(
