计算流体力学基础 ppt课件

计算流体力学基础★课时安排：总学时３２小时，２４小时讲课；8+8小时上机练习。★主要相关前修课程　　计算机语言、工程流体力学、高等数学★主要内容　　介绍流场计算的基本概念、基本方法和简单算例第一章概述1.1计算流体力学的发展及特点简述　　流体力学研究三种方法:实验研究、理论 分析 定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析 和数值计算。☆　实验研究　●真实可靠、是发现流动规律、检验理论和为流体机　　械 设计 领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计 提供数据的基本手段。　●实验要受测量技术限制，实验周期长、费用高。☆　理论研究　●在研究流体流动规律的基础上，建立了流体流动基　　本方程。　●对于一些简单流动，通过简化求出研究问 题 快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题 的解析解。　●对于实际流动问题，通常需运用流体力学基本方程，借助于计算机求数值解（计算机数值模拟）—　　计算流体力学CFD。☆计算机数值模拟●数值模拟耗费小、时间短、省人力，并能对实验难以测量的流动进行模拟，如燃烧室、转子通道内。●在航空航天、核工业、天气预报、海浪和风暴潮预报等方面有极广泛应用。●在航空航天方面，可用于计算飞行器飞行过程中周围流场（计算出升力、阻力）。计算航空发动机各部件内部流场，以及整台发动机三维流场。目前国内有一些使用较多的商用软件，如fluent、Star-CD、numeca等。●美国自上二十世纪八十年代后期，由于CFD方法应用，使一台发动机设计时间从10-15年降到5-8年，试验样机数从40-50台降到10台左右。美国NASA主持建立了推进系统数值仿真系统。☆　数值模拟与实验研究、理论分析关系●三者相互依赖、相互促进●数值模拟所占份额会越来越大（计算机技术迅速发展、计算方法的不断改进）。1.2流场数值模拟概念☆流场数值模拟概念●也称为流场计算机模拟，是以计算机为手段，通过数值计算以数据和图像显示，再现研究对象及其内在规律。　　●数值模拟可理解为用计算机做实验。比如一个机翼绕流问题，通过计算可得到其升力、阻力数值；绕流流线、激波位置、流动分离、涡的生成和传播☆流场数值模拟几个步骤　●建立数学模型：根据流动特点建立适当的数学模型（控制方程）　●确定计算方法　　１）控制方程的离散方法：将流体力学基本方程转化成可用计算机语言描述的形式，称为离散方程，有限元、有限差分、有限体积等。２）边界条件的处理：有／无滑移、壁面等温／绝热等。●编制计算机程序或运用已有程序进行计算　　１）网格生成：在流场中按一定规律分布一些点，称为网格节点。此过程通常称为前处理。无限信息空间用若干个点近似表示　　２）流场计算：运用离散方程求出每一网格节点上气动热力参数值，如温度、压力、速度。　　３）计算结果后处理：根据网格节点上参数值，进一步处理出需要得到的信息，如流动阻力、升力、流量、流线等。　　涡轮叶片通道内三维流计算实例压气机转静子表面压力分布涡轮通道内速度分布航天飞机表面网格航天飞机表面流速矢量图航天飞机表面温度分布气流绕圆柱体流动压力分布气流流过汽车风扇流动直升机旋翼运动NACA0012翼型绕流流线图翼型绕流流线图风力机表面极限流线图轴流叶轮计算与实验叶片表面极限流线轴流叶轮计算与实验性能比较轴流叶轮计算与实验流场结构比较第二章流体力学数值计算数学模型及定解条件　☆本章所涉及的基本方程有两类：　　　●流体力学基本方程，基本出发点：质量守恒、动量守恒和能量守恒　　　●简化模型方程：具有流体力学基本方程的某些特性，用于对所对应的流体力学方程理论分析2．1   可压缩非定常粘性流数学模型连续方程：运动方程：能量方程：上述基本方程构成了Navier-Stokes(简称NS)方程。在三维直角坐标系下Navier-Stokes方程为：上述方程组不封闭，还需要补充数学关系式：１）状态方程：２）物性系数与状态参数关系：2．2   不可压缩非定常粘性流数学模型　　　　当来流Ｍ数小于0.2时，为不可压流动，以下为二种不可压粘性流动控制方程。1）不可压流Navier-Stokes方程连续方程：运动方程：能量方程：2）流函数-涡量方程:对于平面流动:平面流动速度与流函数涡量关系:2.3无粘流数学模型1）欧拉方程：2）全位势方程：上式中:α为音速;3）不可压流全位势方程：2.4常用的模型方程　●流体力学基本方程大都为复杂、非线性方程（组），从数值计算角度分析研究比较困难。并且迄今为止还没有形成成熟的理论。　●为了认识基本方程的数学性质，常用一些简单的线性数学方程作为替代进行研究。　●这些方程具有基本方程的某些特征，称之为模型方程1）对流方程：　　★此方程是双曲型方程，形式类同于一维欧拉方程。2）伯格斯（Burgers）方程：　★是一个非线性方程，具有NS方程类似的性态，式中系数β相当于流体的粘性系数。　　3）对流—扩散方程：　★这个方程和伯格斯方程同属双曲—抛物型方程，但它是线性的，比较简单。　★当β=0时，退化成双曲型方程，当α=0时，则变成抛物型方程4）抛物型方程：5）椭园型方程：　★称为泊松方程，其右端函数项f为已知；　★若f=0,则成为拉普拉斯方程。2．5偏微分方程的数学性质及其与流体运动的关系流体力学基本方程及模型方程属偏微分方程（组），由于方程的复杂性通常无法采用积分方法求精确解，但可将其离散进行数值求解。流体力学方程（组）的数值求解需符合流动的物理规律，同时边界条件的给定也要遵循流动的物理规律，因此首先需了解方程的数学性质。2．5．1拟线性偏微分方程组的分类◇拟（准）线性方程组对于流体力学控制方程，所有最高阶偏导数项都是线性的（这些项前仅有一个系数项，系数项是变量的函数、没有最高阶偏导数与偏导数项的乘积）◇拟线性方程（组）的数学性质以下列拟线性方程组为例式中，系数项是x,y,u,v的函数。u,v是因变量为独立变量x,y的函数，并且u,v是x,y的连续函数。将下式：与以上四式组合在一起并写成矩阵形式可得（2.21）（2.22）令矩阵[A]为上式的系数矩阵，即：并将[A]矩阵的第一列用(2.23)式右侧矢量替代构成矩阵[B]（2.23）根据Gramer法则，有同理可求出du,dv,dx,dy计算：在xy平面内任一点P，过P点作一曲线ab，如果点2无限接近于P点，则：☆ab曲线是任意选定的，其选择不影响计算结果。☆但如果选择的方向使得则无法采用(2.24)计算值ef 称为通过P点的特征线(2.24)所谓特征线即为通过xy平面内某点P的曲线，沿此曲线方向无法确定u和v的偏导数值。因此可通过求解：确定特征线。由展开得：进一步可得由上式可确定xy平面内每一点的特征线斜率，从而确定特征线。如果令：则上式可写成：即：令：，如果在 xy平面内某一点有：1），则偏微分方程组(2.21)有两条各不相同特征线，称方程为双曲型；2），则偏微分方程组(2.21)只有一条特征线，称方程为抛物型；3），偏微分方程组(2.21)没有特征线，称方程为椭圆型。双曲型、抛物型和椭圆型实际上是直接借用以下二次曲线性质2．5．2偏微分方程组分类的通用方法以上根据Gramer法则给出了拟线性方程组类型的确定方法。下面介绍另一种方程组类型通用确定方法。为简单起见，假设方程组(2.21)右端项为0，即：定义矢量：这样式(2.29)可写成矢量形式：(2.30)或者：(2.31)（2.29）上式可变成：上式中矩阵的特征值决定偏微分方程组类型。如果特征值全是实数，方程组为双曲型；如果特征值全为复数，方程组为椭圆型。[例]二维无旋、无粘定常可压缩流，流场中有一细长体，如机翼翼型。如果在上游有一小扰动，扰动速度分量为：。根据连续方程、运动方程和能量方程可推得：为自由来流马赫数。确定以上流动的类型。方法一：采用Gramer法则。对照式（2.21）有：而：于是： 因此当流动超音时，方程组为双曲型；当流动亚音时，方程组为椭圆型方法二：采用特征值方法以下方程：可写成以下矢量形式：所以：由：求出特征值因此采用方法二计算结果与方法一相同。由两个结果比较可看出：上式中的矩阵特征值即为特征线在某一点的斜率。2．5．3流体力学控制方程类型及其对流场数值计算的影响根据具体流动特点，采用的流体力学控制方程组可分为:双曲型、抛物型和椭圆型。一、双曲型方程在二维空间坐标(x,y)下有一点P，对于双曲型方程组有二条特征线通过该点，分别称为左特征和右特征。P点的影响区域仅局限于二条特征线之间下游区域，也就是说，P点产生的扰动影响在此区域可感受到，同时也只有此区域可感受到。影响P点区域仅限于二条特征线之间的上游区域，就是说此区域并且也只有此区域的扰动会影响P点。对于控制方程为双曲型方程的流动问题可采用空间推进方法进行求解。如上图，可给定y轴上的流动参数作为初始条件，然后沿着x轴方向一步一步推进即可求得整个流场。可推得以下几种流体力学控制方程组属于双曲型控制方程组。例1：定常无粘超音速流超音速气流流过一双圆弧机翼，在翼型前缘产生弓形激波，激波后气流仍为超音速。可以证明这种流动控制方程组为双曲型（流动可近似采用小扰动方程描述）。对于此流动可在翼型上游设初始边界ab，边界上流动参数取自由流参数，沿x方向向下游推进即可求得整个流场。例2：非定常无粘流对于非定常的欧拉方程组，无论流动是否超音都是双曲型（关于时间是双曲型的）。对于一维非定常流，在xt坐标系中，阴影部份为P点的影响区域；P点解由在x轴上（即初始时刻，t=0），区间ab数值确定。管道内一维波运动为曲型的一维非定常无粘流例子。通常在流场数值计算中更多采用非定常方程时间推进求定常解。只要边界条件不随时间变化，当计算推进时间足够长时，流动趋于定常、流动参数不再随时间变化，这时得到的解即为定常解。采用非定常方法求定常解的求解过程似乎绕了弯道。实际上对于工程中的有些定常流动问题，采用定常流控制方程无法求解。比如：超音速流绕钝头体的流动，属于超音和亚音混合流动问题。超音区域流动属双曲型；亚音区域流动属椭圆型。在流场计算出来以前无法确定超音区和亚音区的分界线，同时目前还没有对于不同类型的流动都适用的求解方法。将此类定常流动控制方程加入非定常项变成非定常流控制方程，而非定常流动方程无论在亚音区还是超音区都属于双曲型方程（关于时间），因此解决了此类流动不能求解的困难。二、抛物型方程根据前面分析，对于抛物型方程通过任一点只有一条特征线。如图2.5，假设过P点有一条垂直于x轴方向的特征线，则P点的扰动将影响特征线右边的阴影区域。抛物型方程与双曲型方程一样可采用空间推进方法求解。首先给定初始边界ac上数据，沿x方向推进即可求得边界ab和cd间的解。例1：附面层流动对于附面层流动，通过对NS方程进行简化处理得到适用于附面层内流动的简化方程组为抛物型。给定附面层进口边界ab和ef上数值，采用沿壁面方向空间推进即可求出整个附面层内流动。壁面采用无滑移边界条件，bc和fg两个外边界采用无粘流计算结果。采用附面层方程计算附面层内流动，需先给定附面层外边界流动参数。附面层外边界流动参数决定附面层厚度发展，附面层厚度又影响附面层外势流区流动。因此附面层与势流区流动相互影响，需采用迭代方法进行流场计算，计算方法复杂目前已少有人采用。附面层流动分析三、椭圆型方程对于椭圆型方程，无特征线或特征方程是虚根。流场中任意一点P的扰动会向周边任意方向传递。因此边界点的数值同样影响流场中任意一点的解。在所有边界上都要给出边界条件。通常边界条件有以下三种类型：1）给定变量u,v数值，此类边界条件称为Dirichlet边界条件；2）给定变量u,v的导数值，此类边界条件称为Neumann边界条件；3)部份边界给定变量u,v数值，部份边界给定变量u,v的导数值，称为混合边界条件。椭圆型方程影响区域例：定常无粘亚音流动控制方程该方程属椭圆型方程。在此关键是流动亚音，因为对于亚音流，流场中一点的扰动理论上可向各方向传递到无限远处。如下图亚音翼型绕流，翼型上游的流线向上折转，翼型下游的流线向下折转。翼型产生的扰动引起整个流场的变化（理论上直至无穷远处）。亚音速翼型绕流2．6流体力学问题的定解条件　★数学方程建立后，为确定解必须给出定解条件　★定解条件包括初始条件和边界条件一、初始条件　★初始条件就是在某一起始时刻给出流场中速度、压力、密度和温度等参数分布。　　●对于定常问题并不需要初始条件　●实际计算，对于非线性方程（组）要进行迭代求解，需要初始条件作为迭代的初值。　　●初始条件给定不影响最后结果，但初始条件的合理性会影响迭代计算收敛速度，甚至于影响收敛性。二、边界条件　★关于各种流动边界上要给多少个边界条件、给出哪些边界条件，目前还没有一个完善的理论。　　★对于绝大多数工程实际中的流动问题，研究人员根据理论分析结合 经验 班主任工作经验交流宣传工作经验交流材料优秀班主任经验交流小学课改经验典型材料房地产总经理管理经验 都能给出合适的边界条件。　下面介绍一些常见的流动边界及边界条件。1）来流边界（进口边界）　●对于外流流动前方边界称为来流边界；对于内流流动，如进气道和叶轮机内流动进口截面称为进口边界。　　●来流边界理论上应在物面上游无穷远处，在那里流动未受扰动易于给出边界条件　　●在此边界上一般给出：总压、总温、气流角等参数　2）下游边界（出口边界）　　●下游边界（外流流动）和出口边界（内流流动）要设定在绕流体的远下游，在那里流动通过充分掺混已比较均匀，这样有利于边界条件的给定。　　●对于亚音速流，通常给出出口边界上静压（又叫出口反压）；　　●对超音速流，由于下游扰动对上游流动没有影响，因而不能给定出口反压。　　●其他未确定参数如速度、密度、温度以及超音速流的静压等，则采用计算区域内部的数值外插求得。3）壁面边界　●速度的给定　　a)粘性流，流体在壁面边界上的速度等于壁面的运动速度，如果壁面静止，则流体速度为零，即无滑移边界条件b)无粘流，流体在边界处的法向速度为零，而切向速度则由计算求得不再为零，即滑移边界条件●温度的给定a)等温壁，给出壁面温度，并假设壁面处流体的温度与壁面温度相同　b)绝热壁，壁面热流量为零，即：　●压力的给定：壁面法向压力梯度为零，即：第三章有限差分近似及其数学性质●计算流体力学任务是将描述流体运动的偏微分方程转化成离散形式，然后在计算机上求出这些方程的解。　●方程的离散方法有：有限差分法、有限元法、有限体积法等　●有限差分法用差商代替微商，将微分方程转化成差分方程。实现偏微分方程的离散化，以适合于计算机编程计算。3.1差分格式基本概念　●对于一个二维定常问题，如果求解域如图示　●在直角坐标系下，变量可表示成：U(x,y)　●流场中任一网格节点表示为（i,j）,i=1，m；j=1，n　●网格点（i,j）上差分计算值表示为　，它是对函数值　　　　的近似。●●空间步长●时间步长●流体力学方程是偏微分方程，主要由一阶和二阶偏导数项组成●差商代替微商3.2常用偏导数项的差分格式和精度分析3．2．1一阶偏导数差分格式　●一阶偏导数通常有　　　　　　　　　　　　。　●常用的有中心差分和向前、向后差分格式。1.向前差分格式　由泰勒级数：由于所以—截断误差　●　　是步长的一次方，称此差分格式为一阶精度，记作：　在数值计算过程中，时间和空间步长取值都很小，因而截断误差R数值也很小，这样确保用差商替代微商有足够的精度。　●忽略掉截断误差项上式为　　向前差分格式，其精度为一阶。２.向后差分格式精度为一阶３.中心差分格式精度为二阶　☆中心差分比向前和向后差分离散精度高。　☆差分格式选择：a)考虑差分格式的稳定性　　b)在边界上适用性3．2．2二阶偏导数差分格式1.普通中心差分2.普通一侧差分格式3.二阶混合偏导数项通常采用中心差分3.3差分方程和相容性●差分方程：偏微分方程中的偏导数项用差商代替得到的差分形式方程1)　的差分方程a)时间向前、空间中心的差分格式(FTCS)　●差分方程的截断误差：差分方程与微分方程之间存在一个误差，对于上方程：　这个差分方程具有一阶时间精度二阶空间精度　●定解条件离散化a)初始条件　　　　　离散形式b)边界条件　　　　离散形式　由差分方程和定解条件采用时间向前推进可求出n=2、3、4…各时间层上内部节点上的函数值。●格式图：表示差分方程相邻网格节间关联性的图形。图中●表示方程在该网格节点上离散，○表示差分方程所涉及的网格节点。(a)FTCS格式　　　　　(b)FTFS格式　　　(a)FTBS格式b)时间向前、空间向前的差分格式(FTFS)c)时间向前、空间向后的差分格式(FTBS)●差分方程相容性分析微分方程:对应的差分方程:截断误差为:如果有:微分方程与差分方程相容微分方程的定解条件为:对应的差分问题的定解条件:截断误差为:如果有:微分方程与差分方程定解条件相容　●有限差分方法求解流体力学问题举例　　⊙在两固定平板间，流体由于平板二端压差驱动，作层状流动。这种流动称为库特（Couette）剪流。　　⊙不考虑端部效应，在每个等x截面流体速度分布完全相同。　　⊙由NS方程可推得关于速度u(y)的微分方程　　⊙边界条件⊙求出解析解为：★有限差分方法进行求解步骤　１）建立基本方程和适当的定解条件　２）网格划分：沿y方向将线段6等分，则空间步长为　　３）偏微分方程及边界条件差分离散离散方程中心差分离散边界条件节点２节点３节点４节点５节点６４）编制计算机程序进行数值计算，得结果：　　☆这个例子虽然非常简单，但反映了流体力学问题有限差分数值计算的全过程。　　☆通常对于工程实际问题，大部份工作量是花在第２步（网格生成）和第４步（编制计算机程序进行数值计算）。3.4差分方程的收敛性　●差分方程收敛性定义：当步长趋于零时，差分格式的解能否趋近于微分问题的解称为差分格式的收敛性　●对差分网格上的任一网格节点（i,n），设差分格式在此点的解为　　，相应的微分问题解为　　　　　●二者之差为　　称为离散误差。如果有：　　此差分格式是收敛的，即差分方程的解收敛于相对应微分问题的解。否则不收敛。　●条件收敛—差分格式在一定的约束条件范围内是收敛的。　●无条件收敛—指差分格式在任何条件下都收敛。　●关于差分方程收敛性举例。　例１：一维问题如下：精确解为：差分方程:或：因此有:……即：进一步：而：因此：即：所以上述差分方程收敛例2:对流模型方程FTBS差分格式:离散误差为:精确解满足:于是有:若:即有:于是:…….将上面n-1个不等式相加得:　●以上例１为无条件收敛；例２为在　　　　　条件下收敛的差分格式。　●对于流体力学问题由于是多个方程组成的非线性方程组，差分格式收敛性的证明目前还是比较棘手的数学问题。　●目前应用比较广泛的是采用冯.纽曼（VonNeumann）方法—通过差分格式稳定性分析来证明收敛性。3．5差分方程的稳定性及稳定性分析●差分格式的依赖区间、决定区域和影响区域初边值问题:　　1）采用FTCS格式离散2）采用FTFS格式离散3）采用FTBS格式离散FTCS格式　FTFS格式　　FTBS格式差分解的依赖区间和决定区域FTCS格式FTFS格式FTBS格式差分解的影响区域同一微分问题，当采用不同的差分格式时，其依赖区间、决定区域和影响区域是不一样的。进而影响差分方程收敛性。　　●考察步长比对误差传递的影响方程:其解为零,即:采用任何差分格式,若计算中不产生误差,有:假设产生计算误差:采用FTBS格式有:　●上例子显示了该格式的影响区域　●同时显示了　　数值不同时，计算误差对后面时间层上节点的计算值所产生的影响。　　１）　　　，所产生的影响在数值上不会扩大。　　２）　　　，所产生的影响在数值上会越来越大，这样随着计算向前推进，误差会将真解湮没，并最终导致数值趋于无穷，计算失败（发散）。　●在采用有限差分法的运算过程中，计算误差总是不可避免的（如计算机舍入误差）　●有些情况下，误差在传播过程中逐渐衰减；而另一些情况下，误差在传播过程中会逐渐递增、积累。　●若计算中产生的误差，在一定条件下逐渐衰减，那么就称这个差分格式在给定的条件下稳定，这个条件就是它的稳定准则。反之则称差分格式不稳定。☆以对流—扩散方程，说明稳定性概念采用FTCS格式离散得：　　　假设在　　时刻以前的运算中不产生任何计算误差，在　　时刻以后的运算中也不产生新的计算误差，只在　　时刻产生了误差。考察误差传播。这样有：　这时有：　于是：两式相减得：上式称为这就是误差传递方程　　★误差传递方程与原差分方程形式相同上式可改写成：　★此式右边第一项是由于对流项而产生的误差增长，第二项为由于扩散项而产生的误差增长。　★讨论这二项引起的误差增长情况。　　▲误差沿节点分布情况可以是各种各样的。无论误差分布呈何种形态，随着计算由n时间层向前推进，稳定格式误差应逐渐减小，而不稳定格式将逐渐增大。　　▲假设某时刻产生的误差沿节点是振荡的，并且振幅沿节点增加１.不考虑粘性，即无扩散项，。　a)☆结果使误差振幅随时间n单调增加,因而不稳定b)　☆误差随着n的增大逐步减小。格式稳定　☆如果过大,校正会过头，称为过冲——稳定条件2.不考虑对流项,当当☆误差随着n的增大逐步减小。格式稳定☆如果过大,校正也会过头★对流项和扩散项同时存在时，它们各自所产生的误差在传递过程中将会相互影响，在一定的约束条件下，差分格式稳定　●稳定性的数学定义将差分解的误差　　扩展成连续函数Z(x,t)，如果则对应差分格式稳定。K是有限常数。　考虑边界和初始值（即定解条件），稳定性定义式写成：　　▲差分问题在初始时刻或某任一时刻引入的误差为小量，此后的解与差分问题的精确解的误差也一定为小量，所以差分格式为稳定格式。取：则：　　●VonNeumann的稳定性分析方法，又称傅氏级数（Fourier）方法考察对流方程：ＦＴＢＳ格式离散：误差传递方程：　☆误差是节点上的离散量，现将其扩展成空间连续量。这样节点误差的傅氏级数为:代入离散方程得：所以对于任意k值：取：Ｇ—误差放大因子如果：即差分方程稳定●考察对流方程FTBS格式稳定条件要使：必须即为差分方程稳定条件●对流—扩散方程FTCS格式差分方程稳定性分析误差传递方程为：稳定性条件误差的傅氏级数简写成3.6差分方程的相容性、收敛性和稳定性的关系　●前面讨论了差分问题的相容性、收敛性和稳定性　●已经知道相容性是收敛性的必要条件　●发现稳定性与收敛性之间有一定的联系　●Lax等价定理：对一个适定的线性微分问题及一个与其相容的差分格式，如果该格式稳定则必收敛，不稳定则必不收敛。　　　　▲适定指适当的定解条件　　　▲若线性微分问题适定，差分格式相容，则稳定性是收敛性的充分必要条件。　　稳定性线性微分问题适定差分格式相容收敛性　　☆由于收敛性的证明通常比稳定性证明要难，故借助于Lax定理，可将收敛性证明转化成稳定性的证明。第四章模型方程的常用差分格式4．1对流方程的差分格式1.逆风差分格式稳定条件截断误差　　●格式分析　　　为ＦＴＢＳ格式１）　☆a相当于动量方程中的速度u　☆　　在节点i采用向后差商，因此称为逆风差分。　☆根据流动的物理规律，流场中某点的流动参数受上游流动影响比下游大，并且速度越大差别越大。当流动超音时，就不再受下游影响。　☆因而从流动机理分析，采用向后差商比中心差商和向前差商稳定性好（实际上这时采用后两种差分格式是不稳定的）。２）为ＦＴＦＳ格式2.Lax-Wendroff格式由泰勒展开式　　☆采用VonNeumann方法分析可得，这种格式的稳定条件为　　　。　　☆由泰勒级数分析可得截断误差为利用微分方程因此有：稳定条件截断误差3．全隐格式　●上述两种差分格式都有较严格的稳定性条件。为了扩大差分格式的稳定范围，还可以构造隐式格式。　●隐式格式就是在差分方程中，n+1时间层上有多个节点函数值出现。　●显式格式就是在差分方程中，n+1时间层上有只有一个节点函数值出现。　　☆采用中心差商逼近式中一阶和二阶空间导数，并略去高阶小项，得到：格式恒稳，即无条件稳定如果：则得差分方程：　●显式格式时间推进求各时间层的节点函数值过程直截了当。　●隐式格式求n+1时间层节点i的函数值，涉及到相邻节点i-1和i+1的函数值。对流方程：4．2扩散方程差分格式1.古典格式稳定条件为：2.三层全隐式格式4．3对流扩散方程差分格式1.中心显式差分格式稳定条件：2.逆风差分格式稳定条件：3.全隐差分格式4.4计算实例●假设在两相距1m的无限大平板间充满水，平板原来都处于静止状态，在某一时刻t=0，上平板突然以恒定速度平动，求在任意时刻t两板间水的速度分布。●由于两板无限大，可忽略端部效应，这样每一个等x截面速度分布相同。因此只需求x=0截面的速度分布。1）建立控制方程和定解条件●运用非定常二维不可压粘性流的基本方程，根据此流动具体情况对方程进行简化，最后可得到流动的控制方程：●边界条件●初始条件2）网格划分　　●采用均匀网格，网格点数为m，边界节点分别为1和m，那么网格空间步长为3）控制方程和定解条件的离散化　　　●采用古典格式进行差分离散（时间向前、空间中心），得4）编制和调试计算机程序稳定条件为●边界条件的离散形式●初始条件的离散形式　●编制程序要注意程序的条理性、可读性、以及通用性等。　●要有好的条理性，程序则要基于结构化思想进行设计，即一个程序要分成若干块，每一块赋予各自功能，块与块间逻辑关系清晰。　●可读性与条理性是紧密相关的，为了增加程序可读性，还要在程序中适当加入说明语句；以及程序的外观布局、变量所采用的符号都要有所考虑。　●通用性表现在两方面，一是避免程序的局限性，比如对于现在的例子，程序中网格节点数m不要是一个确定的数,而作为变量，这样通过赋值语句对m赋值，可方便地给出不同网格节点数的差分解。其二是在保证程序对确定流场计算模拟功能的前提下，最好也能兼顾计算其他相近流场。START网格划分初始给件n=1边界条件内点值计算结果输出否n=n+1程序流程图程序调试要求１．上机前阅读参考源代码２．完成程序调试３．测试稳定性条件４．修改程序，加入图形显示功能4．5多维问题的几种常用差分格式　●实际问题多为空间二维和三维的问题，上面所述的针对空间一维差分方法向二维和三维的推广并没有太大的理论难度。　●本节以扩散方程为例，将一维差分方法向二维推广，其他情况处理方法类似。1.加权平均差分格式二维扩散方程为相应的定解条件为采用均匀网格：1.加权平均差分格式其中　　☆当　　　时格式是显式的；　　　当时格式是隐式的。　　☆对于二维问题，隐式格式直接求解需要解一个大型稀疏代数矩阵，比较麻烦且耗费机时，后面要进一步介绍对此类问题的处理。●稳定性分析差分方程放大因子：令：则：１）　　　　　：　　格式无条件稳定　　　　　　　２）　　　　　：稳定条件：３）　　　　：或2.交替方向隐式格式（ADI方法）　●对于多维问题，采用隐式格式要求解大型稀疏矩阵，而采用显式格式稳定性限制又较严格。　●交替方向隐式格式综合了显式和隐式格式特点，它的基本思想是将差分计算分成两步。　　★第一步，在一个方向上（比如x方向）是隐式的，而另一个方向是显式的；　　★第二步则两个方向交换，即在第一个方向是显式的，而另一方向为隐式。　●在二步计算中，由于只有一个方向是隐式，这样每一步求解的方程组都是三对角方程组，所以求解过程大为简化。同时格式的稳定性条件比之于显式格式也会大为放宽。　●因为计算在二个方向上交替进行，所以叫做交替方向隐式格式（AlternatingDirectionImplicitmethod,ADI）对扩散方程：采用交替方向隐式格式，得到差分表达式为：　　★在进行n+1/2时间层计算时，涉及到x方向　　　　　　　　　三个节点函数值要同时求得，因而与一维隐式格式相同；同样对于n+1时间层也涉及到y方向三个相邻节点函数值。　　★利用VonNeumann方法进行稳定性分析，格式无条件稳定。●交替方向法另一形式　上两式合并可得　　★因此相当于全隐式格式中加入一高阶小量，此格式是无条件稳定●向三维推广　　☆上式表明：对于三维问题，先在x方向采用隐式格式进行第一步计算；再在y方向采用隐式格式进行第二步计算；最后在z方向采用隐式格式进行第三步计算。　　☆这种格式是无条件稳定的。3.时间分裂格式改写成忽去高阶小量分解成由泰勒级数☆这种格式构造的基本思想是将多维问题化为几个一维问题。相当于求解：稳定性条件：　　　●时间分裂差分格式由于其构造和计算方法简单而有较广泛应用*4．6数值效应第五章不可压流场的数值计算●流动分类：◇定常和非定常流；◇有粘和无粘流；◇可压和不可压流。●从数值计算角度考虑，将流场计算分为不可压流计算和可压缩流计算。◇对于可压缩流，无论其为有粘、无粘、定常或非定常都可以采用同一类方法进行计算。◇但可压缩流的计算方法一般不能用于求解不可压流。5．1不可压无粘流场计算的流函数涡量法5．1．1基本方程推导二维不可压粘性流对于无粘流，且忽略体积力流场无旋上式在二维直角坐标系下展开形式泊松方程5．1．2泊松方程的差分求解中心差分若假设若5．1．2泊松方程的差分求解●差分方程求解方法◇要求解(m-2)*(n-2)个代数方程组成的方程组1)采用常规的代数方程组求解方法来求解，比如:高斯消去法◇这种直接方法求解，如果节点数很多，也就是代数方程个数很多时，就会很耗费计算机内存和机时。2)迭代计算方法◇方法简单，求解速度快,且节省计算机内存。　　一、黎曼方法☆迭代方法相当于利用非定常方程求定常解。☆泊松方程是一个定常方程,要求解这个方程，可以建立一非定常方程构造方程的FTCS差分格式若:黎曼迭代 公式 小学单位换算公式大全免费下载公式下载行测公式大全下载excel公式下载逻辑回归公式下载 黎曼迭代的计算步骤1）将连续函数表示的边界条件转化为离散形式左边界：右边界：下边界：上边界：2）给定初场，即给定初始时刻(n=1)假想的内部节点上的函数值（可采用边界值线性插值）3）采用黎曼迭代公式由左到右、由下而上逐点，或者由下而上、由左到右逐点进行迭代。4）收敛程度判断和输出结果。二、点松驰法松驰因子　　　　☆当　　<1时称为低松驰；当　　>1时称为超松驰；而=1即为黎曼迭代。　　　　☆如果迭代计算过程不收敛，则需减小。但越小收敛速度越慢。超松驰的收敛速度快。三、对于非正方形网格（）1)黎曼迭代公式为2)点松驰迭代公式四、线松驰法●黎曼迭代和点松驰迭代法都是点迭代法，即计算是逐点进行的。在计算(i,j)节点函数值时，周围四个节点(i-1,j)、(i+1,j)、(i,j-1)、(i,j+1)上函数值都作为已知。●线松驰法把与(i,j)节点的同一行或同一列上的相邻两个节点函数值不再当成已知，从而需要与(i,j)节点函数值同时求出。即计算是逐行或逐列进行的。四、线松驰法例如在同一列中进行线松驰迭代，黎曼迭代计算变成1)列迭代2)行迭代　　　☆线松驰法达到规定的误差所需要的迭代次数大约是黎曼法的一半。　　　☆但每一次的迭代计算，线松驰法要解一个三对角矩阵。5．1．3追赶法求解三对角方程组求解方程组i=1:i=2:i=3:i=N-1:i=N:……….设写成矩阵形式为方程组求解1)消去过程（追的过程）…….2）回代过程。SUBROUTINEBANFAC(N)DIMENSIONA（100），B(100)，C（100），F（100），XM（100）C求解三对角矩阵C追赶过程NP=N-1DO1J=1,NPJP=J+1A（JP)=A（JP)/B（J)B(JP)=B(JP)-A（JP)*C（J)F(JP)=F(JP)-A(JP)*F(J)CONTINUEC回代过程XM(N)=F(J)/B(N)DO2J=1,NPJA=N-JXM(JA)=(F(JA)-C(JA)*XM(JA+1))/B(JA)2CONTINUERETURNEND附加:将此程序改成C语言程序并进行测试5.1.4计算举例（内置方形体的突然扩张通道流）　●如下图，流体由一小孔流入突然扩张通道，通道内放置一方形体，假设流动无粘、无旋、不可压且定常，求流场中速度分布。　●YL1=0.1,YL2=2.0,XL1=2.0,XL2=2.0,XL3=3.0。●根据流动特点可得流动控制方程为：●边界条件：　　１）沿ABCDEF是一条流线，给定　　　；　　２）沿IHG也是一条流线，给定　　　。　　３）由于出口截面离绕流体较远，因而可设　　　　　●网格划分　　　　如图中所示采用矩形网格，在此要注意划分的网格要尽可能使流动边界与某一条网格线重合，以便于离散化边界条件的给定。　　　　本算例采用x方向70等分，y方向20等分即可以满足这个要求，这时　●方程离散方程采用中心差分离散，并构造点松驰迭代计算公式　●边界条件离散AB段：BC段：ＣＤ段：ＤＥ段：ＩＨ段：ＨＧ段：ＥＦ段：ＦＧ段：因为要求:1)阅读源程序，并对各模块加说明文字；2)完成程序调试，运用Tecplot画出速度矢量图和流函数等值线图；3）改变松驰因子，测试其对收敛性影响。4）改成小孔流出5．2不可压粘性流场计算不可压粘性流的基本方程为在二维直角坐标系下展开得5.2.1不可压粘性流求解的流函数涡量法1．基本方程●引入流函数和涡量消去上方程中的压力项，可得：●如果需要求得压力场采用时间向前、空间中心差分得：2．涡量边界条件　　设壁面以　　速度运动，由于粘性无滑移条件，壁面上沿x方向速度分量为　　又根据无渗透条件：沿壁面v恒为0，因而有　　由此得上面的涡量边界值计算公式是一阶精度。下面构造更高精度的涡量边界值计算公式。如果设为待定系数在边界上：由泰勒展开式得：比较得：开始n=1输入内部节点和边界节点初始值n=n+1由(5.32)式计算内部节点初始值由(5.35)或(5.37)式计算边界节点值由(5.31)式计算内部节点值由(5.35)或(5.37)式计算边界节点初始值由(5.32)式计算内部节点值输出结果收敛是否由(5.33)式计算流函数涡量法计算流程●平板驱动方腔内流场计算　　　☆如图所示的水槽，水槽横截面为方形，边长1m。底部为一活动平板，以恒定速度U0沿x方向平移。槽内水在平板的带动下一起运动。　　　☆设流动雷诺数Re极小，这时动量方程(5.31)中惯性力项比粘性力项小得多可忽略，则涡量控制方程变为　　　☆这种流动称为Stokes流，控制方程为二维抛物型方程。　　　☆采用FTCS格式差分离散得：☆如果：　　　　　并令　　　☆流函数控制方程：☆流场速度分布计算☆流函数和涡量边界条件　１）底部边界：　２）顶部边界：　３）左边界：　４）右边界：要求1）考察稳定性条件;2）不忽略对流项进行计算并进行结果比较;3)将二阶精度涡量边界条件改成一阶精度.5.2.2不可压粘性流求解的原始变量法1．基本方程对于不可压流，控制方程为◇流函数涡量法推广到三维流动较困难。◇原始变量法——直接以速度和压力为变量采用ＮＳ方程进行数值计算方法。☆二维和三维流动均可采用☆计算过程较流函数涡量法复杂◇原始变量法分为：人工压缩方法、投影法和ＭＡＣ法等。下面介绍人工压缩方法。◇不可压非定常流连续方程中没有关于时间的偏导数项，与动量方程性质不同，方程组不易求解。◇对于求定常解，可考虑在连续方程中加一非定常项。◇连续方程中加入：最合适。☆待解的变量：☆动量方程中有；方程组中尚缺项，◇连续方程变为c为音速，伪压缩方程在二维直角坐标系下展开得2．基本方程离散　　虽然压力分布不均匀，但若根据上式计算压力梯度，则都等于零。　　●为了解决这一问题，选用交错网格对上述方程组进行差分离散。　　◇压力梯度项离散上方程中有x方向和y方向的压力梯度，如采用中心差分，即将y方向运动方程在□所示一类点上展开:将连续方程在用○表示节点上展开将x方向运动方程在△所示一类点上展开：◇边界条件　　●对于粘性流，在固壁边界上速度等于壁面的运动速度。　　1)计算　　　　　　时，涉及到点　　　的函数值　　　，此点在求解域外不存在。为此通常采用线性外插法确定一虚假值，即　　　2)计算　　　　　时，涉及到点　　　　的函数值　　　　，此点也在求解域外不存在。采用类似方法可得：第六章可压缩流场的数值计算6．1可压缩无粘流的差分计算6．1．1一维欧拉方程的显示格式一维欧拉方程为◇欧拉方程的差分格式可直接用于NS方程。虽然粘性流的控制方程（NS方程）比无粘流控制方程（欧拉方程）要复杂得多，其实质是比无粘流多了粘性项，通常不影响数值格式的构造，并且由于粘性项存在会增加格式的稳定性。◇着重介绍无粘流的差分计算方法，然后再进一步推广到粘性流的计算。1.MacCormack格式●二步显式差分格式，第一步称为预报步，第二步称为修正步，第一步采用的是一阶精度向后差分，第二步采用一阶精度向前差分，差分格式总精度为空间二阶精度。●该格式广泛用于求解欧拉方程和NS方程，并在航空航天工业的流场计算中得到大量应用。●MacCormack格式的不足之处是在激波附近数值解中有数值振荡，因而在实际使用时要加入人工粘性项。关于人工粘性在下面的多步龙格—库塔格式中将作介绍。2.多步龙格—库塔格式加人工粘性项欧拉方程可采用三步、四步、五步这样一些多步龙格—库塔格式，在此介绍四步龙格—库塔格式。方程可离散成其中　　　　和　　　　分别为二阶和四阶人工粘性项。　　则取决于节点i的压力二阶导数值，因为：在激波和滞止点附近压力变化较骤烈，　　　　也相应较大，引入的二阶人工粘性就增大。6．1．2多维欧拉方程的显示格式1.MacCormack格式◇一维问题的上述差分计算格式可直截了当地向多维推广。在此以二维为例，其欧拉方程为：2.多步龙格—库塔格式6．2可压缩粘性流的差分计算◇当流动不产生分离或不产生较大分离时，采用欧拉方程进行数值模拟是可行的。欧拉方程由于忽略了流体的粘性，无法计算流动的损失和阻力。◇但当流动存在较大的分离区或绕流的相互干挠区时，则需要考虑流体的粘性，即需要采用NS方程进行数值模拟。比如：激波与附面层的相互干挠区◇非定常NS方程是具有双曲—抛物性质的方程。方程由无粘项（对流项）和粘性项组成。在远离固体壁面区域（势流区）方程呈双曲型；而在固体壁面附近（附面层区），粘性占主导地位，这时方程呈抛物型。◇粘性项描述了流动物理量的耗散特性☆它使流场中的流动参数分布趋于平滑。☆粘性项的差分离散一般采用中心差分即可。☆粘性项的引入，从差分格式构造角度考虑，不会带来什么困难。◇附面层区流动参数变化较剧烈☆在此区域网格节点分布要很密集。☆空间步长较小，时间步长也相应受到很大限制，☆大大增加计算工作量。一维可压缩粘性流的差分计算MacCormack格式的差分方程　　●二维和三维的差分格式构造类同与一维。●在第七章网格生成基础上，第八章要专门介绍任意曲线坐标系下三维粘性流计算。一、贴体坐标　　　◇园柱体绕流问题如果流动无粘且无旋，则有：节点A上的函数值处理方法1）2）第七章流场网格生成　　　如果上述问题采用柱坐标，这时有坐标变换关系式：边界条件控制方程◇坐标转换的意义　　●在xy坐标系下曲线边界的问题通过坐标变换转化成了柱坐标下的矩形边界．　●在此称直角坐标系下真实流动区域为求解域；而通过坐标转换所得到的对应变形区域为计算域。　●实现在矩形计算域内数值计算。采用下式表示求解域与计算域之间的坐标变换关系：—为任意曲线坐标◇坐标转换关系　☆一维坐标转换设有一维坐标转换问题　反函数关系所以有☆二维和三维坐标转换二维坐标转换问题反函数因此有而所以☆采用同样方法可得三维坐标转换关系式其中Ｊ——坐标转换矩阵雅可比（Jacobin）。在此以二维流动为例直角坐标系下的二维欧拉方程◇任意曲线坐标系下的基本方程坐标转换关系——任意曲线坐标系下二维守恒型欧拉方程　◇坐标转换关系差分表示　　如果采用流函数等值线构成　　网格线，而势函数等值线构成　网格线，则有：　☆微分方程生成方法有：双曲型方程、抛物型方程以及椭圆型方程生成方法。☆基本思想假设流体在由内边界HABCD和外边界GFE构成的“C”型流道内流动，流体不可压无粘无旋。◇椭圆方程方法网格生成　　在直角坐标系下展开得：　边界条件设定为：　　◇通过数值求解方程组：◇采用坐标变换　　◇得到坐标变换数值关系：◇将　　　　　　　　等值线在物理域中画出，即得到相应的网格。其中　　如果加上一些辅助项能使坐标拉伸或压缩，进而对网格线分布密度进行调整。得：◇网格节点分布控制对于一维流动：它的解是——建立了坐标转换关系。显然这是一个线性转换，坐标上的等间距节点分布变换到x坐标上也是等间距。◇如果方程改成它的解是：为了方便起见，设端点条件为：则有　　——　　坐标上进行等间距节点分布，对应x坐标上已不是等节点间距．　　☆当c＞0时，x=0处的节点加密；　　☆当c<0时，x=0处的节点变稀疏。　　◇在确定的所对应的附近改变节点分布的疏密程度，可以采用下列方程　　☆a>0在　　附近网格加密；a<0在　　附近网格变疏。　　☆c越大，指数函数衰减越快，影响区域越小。☆如果有n个节点附近网格疏密程度要调整，则可采用下列关系式：将上述方法推广到二维问题中，得到：将方程转换得：方程离散形式（取　　　　　）：采用黎曼点迭代◇网格生成实例以下为圆柱体绕流，其物理域流场外边界为一圆，内边界为圆柱体。初始网格生成的网格计算域网格◇翼型绕流第八章三维紊流平均流的有限差分计算　　●由于计算机技术的迅速发展，计算机内存和计算速度都有了很大提高。早先由于计算机水平限制很多流动问题数值计算采用无粘假设。现在几乎所有二维问题以及大部份三维问题都可以采用粘性计算，即采用NS方程。　　●紊流流动是非定常不规则流，目前的计算机水平对此种流动进行数值计算仍存在较大困难。采用紊流平均流NS方程结合紊流模型是当前普遍采用的方法。　　●工程上所涉及的流体力学问题绝大都数属于不规则边界，采用有限差分方法进行计算则必须将直角坐标系下的基本方程转化到任意曲线坐标系下。　　●本章将简述任意曲线坐标系下三维紊流流动四步龙格—库塔计算方法。8．1三维紊流平均流N-S方程　　在任意曲线坐标系下，非定常可压缩三维紊流平均流NS方程可写成：其中：U，V，W为逆变速度分量和分别为层流和紊流粘性系数8．2Baldwin-Lomax双层代数紊流模型　　●紊流模型可分为二类：基于涡粘性概念的紊流模型和雷诺应力及代数应力模型。代数紊流模型假设局部紊流生成与耗散平衡，将涡粘性直接用附面层特征量的代数关系式表示。　　●如Baldwin-Lomax模型(1978)和Michel模型(1969),此类紊流模型计算过程简单，计算速度快，无需求解偏微分方程（组），因而得到广泛应用。8．3控制方程的空间离散对三维NS方程的空间偏导数项采用中心差分进行离散，得：8．4人工粘性项●附面层外势流区由于不存在物理粘性，中心差分会造成奇偶点数值不耦合。为了计算过程的稳定性，防止激波和滞止点附近数值振荡，通常计算过程中都要引入人工粘性。●将方程改写成C（D）、D（Q）分别为对流项和物理粘性项，AD（Q）为人工粘性项。8．5控制方程的时间离散采用四步龙格一库塔对守恒方程进行时间积分，为此将方程改写成：第九章流场数值计算的新方法简介●为了能采用数值方法模拟复杂流动和提高数值计算结果的精度，需要有计算速度更快、内存更大的计算机。此外还有一个可行的途径是改进数值计算方法。9.1多重网格法●一个三维粘性流场的数值计算要花费数小时或数天时间。因此寻求一些加速求解代数方程组的计算方法是很有意义的。●近十多年来发展起来的多重网格法（MultipleGridMethod）是一种加速收敛的比较好的方法。它首先用来求解椭圆型方程，后来又被推广到对时间相关流动问题的加速求解。●求解一个流场首先要进行网格划分，为了使数值解有符合要求的精度，网格划分得要足够细、节点数要足够多。这样计算迭代时间会比较长。（点迭代计算时间为：　　，　　　　　　）●先将网格划得稀疏一些，在此稀疏网格上进行计算得到一个初步结果。然后再加密网格，其初值采用在稀疏网格上计算结果插值得到。●这样得到的初场比较合理，因而可以使迭代计算过程加快。另外稀疏网格上计算不需要有很高的精度，因此为了减少计算时间，稀疏网格上迭代计算次数也不需太多。●多重网格法就是基于上述思想。此外多重网格法需要在粗细网格上交替进行计算，其目的是将误差通过滤波方法减少。9.1.1迭代法的误差衰减黎曼迭代公式以泊松方程为例来讨论，其一般形式为为常系数，为已知函数。边界条件为：误差传递方程将误差项用傅氏级数展开，如：得到：其中：因此放大因子为：设　　　　　　　　　　　　　　上式可变为对于不同的　　　　有不同的误差衰减率，如：　　●对于同一种网格而言，在迭代计算过程中，波数比较大时，误差衰减率就大（放大因子小），而波数小时衰减慢。也就是说误差的高频分量衰减的比较快，低频分量衰减得比较慢。　　●也有些高频分量衰减得比较慢，比如：　　●傅氏级数是将误差分解成不同频率的误差分量之和。由上可看出：在迭代计算过程中，不同频率的误差分量衰减的速度不同，而衰减速度缓慢的误差分量将制约整个迭代计算的收敛速度。　　●以上是点迭代，如果改用线迭代，比如采用列迭代，差分方程为：放大因子表达式为：因此有：　　●与点迭代比较，线迭代可以加速一些高频和中频分量的衰减，但对于低频和某些特别的高频误差分量仍然不能加速衰减。　　●上面讨论的是同一种网格不同波数的误差分量衰减情况。●对于同一个波分量，如果改变了，它的衰减速率也会相应改变。１）则如果　　●由上分析，如果在流场数值计算过程中，采用多种网格尺寸，并且反复在不同的网格上进行迭代计算，会产生互补性，使大部份误差分量较快衰减，迭代过程得以加速。　　●不同尺度的网格称为不同层次，由多个层次的网格构成多重网格。在多层网格上反复交替进行计算以加速收敛的迭代过程就叫做多重网格法。则衰减速率就提高２）9.1.2多重网格法的计算过程设在流场计算中将网格划分为M+1层，网格由粗到细标记为0，1，2，…k…，M；对应k网格层的步长为:记原方程(9.1)为在k层离散得到的差分方程为最后要求解的差分方程设是其近似解为由于解近似引起的亏损，若记则得亏损方程●线性多重网格法是在M层网格上求精确解和在k层网格上求解亏损方程的组合。　　●首先考察只有两层网格的计算过程☆转移计算：　　在计算过程中两层网格节点上的数值要相互交流，　　由细网格到粗网格的转移，表示为：　☆网格构成　　粗网格取细网格的间隔节点构成。粗网格空间步长为细网格两倍；粗网格节点数为细网格的１／４。　由粗网格到细网格的转移，表示为：　限制算子　插值算子　1）限制算子　a)粗网格上的节点函数值直接采用细网格上与之重合节点函数值　b)采用相邻九个节点值进行加权平均，即2）插值算子双重网格的计算过程1)给定细网格上初值，对差分方程作次迭代（如松驰迭代则称为次松驰）。一般=1-2，得到近似解2)粗网格修正a)计算细网格上亏损量b)由细网格向粗网格转移亏损量c)在粗网格上求差分方程的准确解d)由粗网格向细网格转移修正量e)计算细网格上修正后的量然后再以为初值重复1）和2）步，直至收敛为止●多重网格法的计算过程☆网格步长记作:　　对于均匀网格一般取：☆h为0网格层（即最粗网格层）上的网格间距。☆每个网格层上相应的算子为　　　　　　　　　　　　　　　　　　或简记成　转移算子：由粗网格到细网格为插值算子　　（k=1，2，…，M）由细网格到粗网格为限制算子　　（k=1，2，…，M）☆求解方程为：☆计算步骤如下：1）   给定初值　作　次松驰，即：2）   粗网格修正a)     计算M层亏损量---------------------b)     将M层的亏损量向M-1层转移---c)     采用几步松驰近似求解-----------d)     再计算M-1层亏损量--------------e)     将M-1层亏损量向M-2层转移---f)     采用松驰法近似求解--------------g)求最粗网格层上的亏损量计算方程高精度解（准确解）3）由　通过插值算子　得到------●举例，求解：　　　用超松驰法（SOR）、交换方向隐式格式（ADI）及多重网格法（MG）所需时间如下表：边界条件：方法运算次数CPU时间（IBM/370-158）SOR~1200SADI~NlgN127SMG~N7.7S9.2自适应网格　　●在进行网格生成时要求在流动参数变化较骤烈的区域网格点分布密集，比如激波附近、尾流区、附面层区等。　　●在生成网格前要先估计在流场中哪些区域流动参数变化较快。　　●对于粘性绕流可以认定在物面附近的附面层区流动参数沿物面法向变化较快　　●对于有些流动问题，比如超音速流，流场中有无激波、激波位置往往很难预估。　　●所以在计