第三章数值模拟理论与方法§3.1流体力学的基本方程流体运动所遵循的规律是由物理学三大守恒定律规定的，即质量守恒定律，动量守恒定律和能量守恒定律[44]。（一）连续方程0)(vt（3.1）式中ρ－流体密度u－流体速度分量（二）动量方程（x方向）对于不可压流体（即0v）xpfvuvxux)()()(（3.2）式中γ－运动粘性系数p－压力对于可压缩流体xpfvxuuvxux31（3.3）式中等号后前两项是粘性力y，z方向上的动量方程可类似推出。（三）能量方程vqTkevet（3.4）其中TCev式中等号左边第一项是瞬变项，第二项是对流项，等号右边第一项是扩散项，第二、三项是源项。所以，流体力学基本方程组为：0vtxpfuuvftux)(ypfvvvftvy（3.5）wpfwwvftwwvqeckevfetv§3.2紊流模式理论概况§3.2.1基本方程在自然界中，真实的流体都具有粘性。粘性流体存在两种不同的运动方式和流态，即层流和紊流。而在自然界和工农业生产中所遇见的流体流动大部分都是紊流。三维的N-S方程是目前描述粘性流体运动较为理想的模型，其优点一是应用范围广，在空气、水流、传热等方面均用N-S方程描述；二是对于有分离、旋涡等情况的复杂三维流动更为适用。三维直角坐标下的N-S方程[45]，[46]，即不可压缩粘性流体的动量方程式为：)()()(222222222222222222zwywxwzpFDtDwzvyvxvypFDtDvzuyuxuxpFDtDuzyx（3.6）不可压缩流体的连续性方程为：（3.7）式（3.6）和（3.7）共有四个未知数（u、v、w、p）和四个方程，加上边界条件，从理论上来讲其解是存在的。但是，要直接求解复杂而详细的粘性流体运动是十分复杂和困难的。其原因是：直接求解N-S方程要求求解从反映消散运动的最小涡漩尺度到反映大尺度涡体的所有流动尺度，因而只有对简单情况下才有理论解。0zwyvxu§3.2.2三维N-S方程N-S方程模型的流动计算可分为三种方法：1.直接模拟法（DirectNumericalSimulation，DNS）除稀薄气体等极端条件外，紊流的最小长度尺度远远大于分子运动的长度尺度，故紊流可以作为连续体运动处理。从原理上讲，可以用三维非定常的N-S方程对紊流进行直接计算。这种直接计算不需要紊流模型化，可像层流那样进行数值计算。但是，现实的高雷诺数紊流中，由于其最小尺度很小，若要对最小尺度的紊流进行直接计算，就需要很多的计算时间和庞大的计算机容量。这远远超过现有的计算机能力。当前直接计算法只能用于对低雷诺数紊流进行直接计算，并且用新型巨型向量计算机可取数十万个网格点，但也只能捕捉到较大的紊流涡，网格的网目捕捉不到小涡，从而得到的仅是关于大涡结构的大体结果。将来，即使可能进行精确的直接计算，但为了获得有意义的信息，也必须对大量的计算结果进行统计处理。2.大涡模拟法（LargeEddySimulation，LES）依照紊流的旋涡理论，紊流的脉动与混合主要是有大尺度的涡造成的。大涡从主流中获取能量，分裂后将能量传到较小的涡。大涡的运动为各向异性，随流动情况而不同。小涡主要是耗散能量，几乎各向同性，并且不同流动情况的小涡有许多共性。从而得出大尺度涡模拟的数值方法。即用非定常的（三维且时间相关的）N-S方程确定大涡的特性，不计算小涡。而小涡的效果有近似的模型来处理，即用大涡模拟还可以对那些被直接计算忽略掉的，比如计算网格小的涡，经模型化，进行数值模拟。该方法需要相当大的计算机内存和计算时间。3.雷诺（Reynolds）时均方程法将非定常的N-S方程作时间平均处理。在所得出的时均方程中包含了脉动量乘积的时均值未知数，于是方程个数少于未知数个数，如作出进一步的时均处理将出现更高阶的脉动量乘积的时均值未知数，方程不可能封闭；要是方程封闭，须作出一定的假设。这是工程上普遍采用的方法，因为工程中感兴趣的是时均量。在三维N-S方程计算模型中，雷诺时均方程法是较常使用的一种方法。该方程是在将紊流看成时均运动和脉动运动的基础上建立的。紊流运动的任何变参量都分解为时间平均值和脉动值，例如：'iiiuuu，'ppp等。不可压缩粘性流体的三维N-S方程组作时均处理后的时均方程为：连续性方程：0iixu,0'iixu（3.8）动量方程（雷诺方程）：)(''jiijijjjuuxuxxpFDtuD（3.9）式中：''jiuu为二阶相关项，又称为雷诺应力...
