计算流体力学清华大学完整版课件

计算流体力学引论TheElementsofComputationalFluidDynamics计算流体力学引论TheElementsofComput1教材：任玉新,陈海昕.《计算流体力学基础》，

清华大学出版社,北京,2006。预修课程：流体力学、

偏微分方程数值解法、

计算机语言和编程基础。

参考书目：J.D.Anderson,Jr.ComputationalFluidDynamics-TheBasiswithApplications,McGraw-Hill,NewYork,1995.J.H.Ferziger,M.Peric.ComputationalMethodforFluidDynamics,Springer—Verlag,Berlin,2002.计算流体力学引论教材：任玉新,陈海昕.《计算流体力学基础》，预修课2教学目的、要求

本课程是流体力学及相关学科(地球科学、环境流体力学、化学、石油工程等)研究生的专业基础课，主要讲述计算流体力学基础理论及其应用。本课程重点介绍有限差分和有限体积方法的基本概念、基本理论和部分典型数值方法，阐释计算流体力学求解问题的思路，使学生能够掌握计算流体力学的基本概念，具备初步解决模型问题的能力。计算流体力学引论教学目的、要求

本课程是流体力学及相关学科(地球科学、环境3课程考核：作业（30%）

期末考试（70%）

期末考试：闭卷笔试

计算流体力学引论课程答疑：周二，13:30~15:10，N606

课程考核：作业（30%）期末考试：闭卷笔试4计算流体力学引论TheElementsofComputationalFluidDynamics计算流体力学引论TheElementsofComput5第一章绪论§1.1计算流体力学的概念与意义§1.2流体力学的基本方程§1.3流体力学方程组的类型判别

第一章绪论§1.1计算流体力学的概念与意义6§1.1计算流体力学的概念与意义1、流体运动遵循3个基本定律：1)质量守恒定律；2)动量守恒定律；3)能量守恒定律2、流体的本构模型和状态方程

控制方程（Governingequations）偏微分方程（方程组）或积分形式的方程（方程组）

流体运动的复杂性主要表现为控制方程的高度非线性和流动区域几何形状的复杂性等，导致对绝大多数流动问题无法得到解析解。高速计算机的发展，使得计算流体动力学（ComputationalFluidDynamics,CFD）逐渐成为一门独立学科。§1.1计算流体力学的概念与意义7计算流体力学(CFD)：通过数值方法求解流体力学控制方程，得到流场的离散的定量描述，并以此预测流体运动规律的学科。在CFD中，首先，把控制方程中的积分、微分项近似地表示为离散的代数形式，把积分、微分形式的控制方程转化为一组代数方程，这个过程称为控制方程的离散化(discretization)；所采用的离散化方法称为数值方法或数值格式。然后，通过电子计算机求解这些代数方程组，得到流场在离散的时间/空间点上的数值解(numericalsolution)。CFD也被称作流场的数值模拟、数值计算、数值仿真等。计算流体力学(CFD)：通过数值方法求解流体力学控制方程，得8计算流体力学的研究步骤第一，问题的界定和流动区域的几何描述。

流场的几何形状：源于对已有流动区域的测量或者新的产品和

工程的设计结果。

流动条件：雷诺数、马赫数、边界处的速度及压力等

对数值模拟的要求：精度、所花费的时间。第二，选择控制方程和边界条件。

在牛顿流体范围内，用Navier-Stokes方程描述。

根据问题的特点，可以考虑定常或非定常，可压或不可压的流动模型。简化的数学模型：势流方程，Euler方程，边界层方程，

薄层近似的Navier-Stokes方程等。

边界条件通常依赖于控制方程。

固体壁面条件，来流、出流条件，周期性条件，对称条件等

附加的物理模型：湍流模型，化学反应等。计算流体力学的研究步骤第一，问题的界定和流动区域的几何描述。9第三，确定网格划分策略和数值方法。

网格划分：结构网格、非结构网格、组合网格、重叠网格。

网格可以是静止的，也可以是运动的，还可以根据数值解动态调整（自

适应网格）。数值方法：有限差分、有限体积、有限元、谱方法等。

数值方法和网格划分策略是相互关联的。第四，程序设计和调试。

在网格划分策略和数值方法的基础上，编制、调试数值求解流体运动方程

的计算机程序或软件。第五，程序验证和确认。

验证(Verification)：Theprocessofdeterminingthatamodelimplementationaccuratelyrepresentsthedeveloper’sconceptualdescriptionofthemodelandthesolutiontothemodel.

确认(Validation)：Theprocessofdeterminingthedegreetowhichamodelisanaccuraterepresentationoftherealworldfromtheperspectiveoftheintendedusesofthemodel.第三，确定网格划分策略和数值方法。10第六，数值解的显示和评估计算感兴趣的力、力矩等；应用流场可视化软件对流场进行显示、分析；对数值方法和物理模型的误差进行评估等。第六，数值解的显示和评估11计算流体力学典型流程物理模型数学模型网格生成离散方法选择时、空离散边界条件离散解代数方程组验证与确认流场显示结果分析计算流体力学典型流程物理模型数学模型网格生成离散方法选择时、12举例：自然循环回路内的流动与传热特性举例：自然循环回路内的流动与传热特性13物理模型：

(1)空间维数：1D、2D、3D(2)时间特性：定常、非定常

(3)流动性质：无粘/粘性、可压缩/不可压缩、层流/湍流

(4)流体物性：常物性、变物性Geometricparameter:HeightHWidthWLengthofheatsink(source)L

TubediameterdRayleighnumberRaHeatsourcetemperatureThHeatsinktemperatureTcOperationpressureP物理模型：Geometricparameter:14数学模型：控制方程

定解条件

初始条件：

边界条件：固体壁面上无滑移；

恒温热源、恒温热沉，

其余为绝热壁面。数学模型：15网格划分：网格划分：16数值算法：离散方法：

FDM、FVM、FEM……空间离散：对流项，粘性项，源项……时间离散：显式、隐式边界离散：来流、出流、固壁、远场、周期性……求解代数方程组数值算法：离散方法：17数值解的验证与确认：数值解的验证与确认：18流场显示及结果分析：流场显示及结果分析：19计算流体力学的特点及意义实验研究理论研究计算流体力学优点：借助各种先进仪器，给出多种复杂流动的准确、可靠的观测结果，这些结果对于流动机理的研究和与流体运动有关的机械和飞行器的设计具有不可替代的作用。缺点：费用高昂，周期很长，有些流动条件难以通过实验手段来模拟。优点：可以给出具有一定适用范围的简洁明了的解析解或近似解析解，这些解析解对于分析流动的机理和预测流动随参数的变化非常有用。缺点：只能研究简单流动问题，能够得到解析解的流动问题为数不多，远远不能满足工程设计的需要。发展CFD的主要动机：利用高速电子计算机，克服理论研究和实验研究的缺点，深化对于流体运动规律的认识并提高解决工程实际问题的能力。优点：原则上可以研究流体在任何条件下的运动，使得我们研究流体运动的范围和能力都有本质的扩大和提高。费用低，周期短。计算流体力学的特点及意义实验研究优点：借助各种先进仪器，给出20§1.2流体力学基本方程守恒型积分方程牛顿流体本构关系Stokes流体假设§1.2流体力学基本方程守恒型积分方程牛顿流体本构关系S21守恒型微分方程积分型方程和微分型方程在意义上有微妙差别：积分型方程允许在控制体内部流动参数有间断；微分型方程假定流动参数是可微的，因而是连续的。因此，积分型方程是比微分型方程更为基本的方程，尤其是流场中确实存在间断时。状态方程守恒型微分方程积分型方程和微分型方程在意义上有微妙差别：状态22直角坐标系下的守恒型方程Navier-Stokes

方程直角坐标系下的守恒型方程Navier-Stokes方程23计算流体力学清华大学完整版ppt课件24计算流体力学清华大学完整版ppt课件25Euler方程等价形式Navier-Stokes

方程中，Euler

方程中，CFD中，守恒型方程是使用最频繁的一种形式。Euler方程等价形式Navier-Stokes方程中，26边界条件黏性流动的适定边界条件：

在固体壁面上速度满足无滑移条件：

温度条件可以是下面三种之一：无黏流动的适定边界条件在固体壁面上速度满足不可穿透条件边界条件无黏流动的适定边界条件在固体壁面上速度满足不可穿透27§1.3偏微分方程的分类及数学性质一阶拟线性方程组Euler方程：一阶非线性偏微分方程组

Navier-Stokes方程：二阶非线性偏微分方程组流体力学的基本方程都可以写成一阶拟线性方程组的形式。对一阶导数项而言，是线性方程组；如果B,A是U的函数，则整个方程组是非线性的，称之为“拟线性方程组”。§1.3偏微分方程的分类及数学性质一阶拟线性方程组对一阶28考虑一维守恒型Euler方程（一阶）考虑一维守恒型Euler方程（一阶）29令令30考虑Laplace方程（二阶）作业一：根据类似的方法，将Navier-Stokes方程写成一阶拟线性方程组的形式考虑Laplace方程（二阶）作业一：根据类似的方法，将Na31特征线理论分析拟线性方程的特征线和相容关系具有重要意义。通过引入特征线和相容关系，可以把偏微分方程的某种线性组合化为常微分方程。在有些情况下，还可以由此得到解析解。考虑一般形式的有两个自变量的拟线性方程组，它的分量形式特征线理论32计算流体力学清华大学完整版ppt课件33计算流体力学清华大学完整版ppt课件34双曲型方程的定义双曲型方程的定义35计算流体力学清华大学完整版ppt课件36计算流体力学清华大学完整版ppt课件37双曲、抛物和椭圆型方程的数学性质不同类型的方程，如双曲、抛物、椭圆型方程具有不同的数学行为，对应着不同的物理过程；因而，也应采用不同的数值方法求解。双曲、抛物和椭圆型方程的数学性质不同类型的方程，如双曲、抛物38计算流体力学清华大学完整版ppt课件39计算流体力学清华大学完整版ppt课件40计算流体力学清华大学完整版ppt课件41计算流体力学清华大学完整版ppt课件42计算流体力学清华大学完整版ppt课件43计算流体力学清华大学完整版ppt课件44计算流体力学清华大学完整版ppt课件45计算流体力学清华大学完整版ppt课件46计算流体力学清华大学完整版ppt课件47计算流体力学清华大学完整版ppt课件48流体力学方程组的其它类型流体力学方程组的其它类型49计算流体力学清华大学完整版ppt课件50计算流体力学清华大学完整版ppt课件51计算流体力学清华大学完整版ppt课件52计算流体力学清华大学完整版ppt课件53计算流体力学引论TheElementsofComputationalFluidDynamics计算流体力学引论TheElementsofComput54第二章有限差分方法基础§2.1有限差分方法概述§2.2导数的数值逼近方法§2.3差分格式的性质§2.4发展方程的稳定性分析第二章有限差分方法基础§2.1有限差分方法概述55§2.1有限差分方法概述以一维非定常热传导方程为例，介绍有限差分方法的概念、简单构造方法和求解过程。2.1.1基本方程和定解问题方程(2.1.1)和初边条件(2.1.2)构成了一个适定的定解问题。有限差分方法：对于一个偏微分方程，如果把方程中的所有偏导数近似地用代数差商(AlgebraicDifferenceQuotient)代替，则可以用一组代数方程近似地替代这个偏微分方程，进而得到数值解，这种方法称为有限差分方法(FiniteDifferenceMethod)。§2.1有限差分方法概述以一维非定常热传导方程562.1.2求解域及偏导数的离散化为了用有限差分方法求解式(2.1.1)，需要把其中的偏导数表示为代数形式，为此，首先要把自变量从连续的分布变为离散形式。这个过程称为求解域的离散化。

1.空间求解域的离散化把空间求解域分为M段（均匀剖分）

2.时间变量的离散化把感兴趣的时间段(t=T之前)分为N段（均匀剖分），则时间方向的求解域可以划分为2.1.2求解域及偏导数的离散化为了用有限差分57

求解域被划分为一系列离散的时空网格点图2.1求解域的离散化

3.解的离散表示目标：求出所有网格点上物理量u的近似解。求解域被划分为一系列离散的时空网格点图2.158

4.导数的数值逼近把方程中的偏导数项近似表示为代数形式。4.导数的数值逼近把方程中的偏导数项近似表示为代59计算流体力学清华大学完整版ppt课件60计算流体力学清华大学完整版ppt课件61计算流体力学清华大学完整版ppt课件622.1.3差分格式同一偏导数可以有不同的近似方法，不同的导数近似方法导致方程的不同的有限差分近似。FTCS(ForwarddifferenceinTime,CentraldifferenceinSpace)格式时间方向用前差近似，空间二阶导数用中心差分近似。对初始条件和边界条件的离散化式(2.1.9)~(2.1.12)称为方程(2.1.1)的一个有限差分方程或有限差分格式(finitedifferencescheme)。2.1.3差分格式同一偏导数可以有不同的近似方法，不同的632.BTCS(BackwarddifferenceinTime,CentraldifferenceinSpace)格式时间方向用后差近似，空间二阶导数用中心差分近似。在研究数值方法时，通常把tn时刻的物理量视为已知量，而把tn+1时刻的物理量作为待求的未知量。因此，式(2.1.13)可以改写成2.BTCS(Backwarddifference642.1.4差分方程的求解FTCS格式可以改写为可见，在FTCS格式中，某一点的数值解只依赖于前一时间步的三个点，如图2.2所示。图2.2：FTCS格式的模板点2.1.4差分方程的求解FTCS格式可以改写为可见，在65FTCS格式的求解过程FTCS格式的求解过程662.BTCS格式可以改写为跟FTCS格式不同，BTCS格式中同时涉及到n+1时刻的多个未知量，不能递推求解，称为隐式格式(implicitscheme)。图2.3：BTCS格式的模板点BTCS格式的求解过程2.BTCS格式可以改写为跟FTCS格式不同，BTC67计算流体力学清华大学完整版ppt课件68计算流体力学清华大学完整版ppt课件692.1.5用时间相关方法求解定常问题考虑非定常热传导方程和定解条件2.1.5用时间相关方法求解定常问题考虑非定常热传导方程70计算流体力学清华大学完整版ppt课件71BTCS格式的求解过程FTCS格式的求解过程BTCS格式的求解过程FTCS格式的求解过程72§2.2导数的数值逼近方法2.2.1精度分析在上一节，我们得到了一阶偏导数的前差、后差和中心差分近似，以及二阶导数的中心差分近似。这些近似方法逼近偏导数的程度如何呢？可以用Taylor展开式进行分析。§2.2导数的数值逼近方法2.2.1精度分析73计算流体力学清华大学完整版ppt课件74一般来讲，对偏导数的近似精度越高，差分格式的精度越高。一般来讲，对偏导数的近似精度越高，差分格式的精度75例：一维非定常热传导方程的FTCS格式中涉及的导数差分近似的精度。例：一维非定常热传导方程的FTCS格式中涉及的导数差分近似的762.2.2导数差分近似的待定系数法2.2.2导数差分近似的待定系数法77计算流体力学清华大学完整版ppt课件78计算流体力学清华大学完整版ppt课件79计算流体力学清华大学完整版ppt课件802.2.3导数差分近似方法的差分算子法1.差分算子的定义算子，一种前置运算符。算子和它后面的作用量一起代表一种确定的运算过程。在引入差分算子的定义之前，先介绍一种特殊的算子——移位算子。移位算子的运算规则为移位算子的下标表示移位的方向，上标表示移位的步数。2.2.3导数差分近似方法的差分算子法1.差分算子的定81差分算子：移位算子和可以表示为移位算子函数的算子。差分方法中常用的算子：差分算子：移位算子和可以表示为移位算子函数的算子。差分方法中822.差分算子之间的关系2.差分算子之间的关系83所有的差分算子均可用Taylor展开式来估算截断误差项的量级。所有的差分算子均可用Taylor展开式来估算截断误差项的量级843.微分算子与差分算子的关系3.微分算子与差分算子的关系854.导数的近似根据差分算子之间的转化关系，可以建立微分算子与其它差分算子之间的联系，从而得到导数的数值近似公式。即：4.导数的近似根据差分算子之间的转化关系，可以86即：与待定系数法得到的结果一致。即：与待定系数法得到的结果一致。87即：即：885.紧致格式从上面的推导可以看出，导数的有限差分近似精度越高，所需要的模板点越多。对于一阶导数，一般需要5个点才能得到四阶精度的差分近似。模板点数太多不仅使数值方法变得复杂，也给边界附近的处理带来一定困难。紧致格式：用较少的模板点构造导数的高阶近似。5.紧致格式从上面的推导可以看出，导数的有限差89计算流体力学清华大学完整版ppt课件90基于Pade近似的导数近似方法，称为紧致格式(compactscheme)。基于Pade近似的导数近似方法，称为紧致格式(compac91计算流体力学清华大学完整版ppt课件92计算流体力学清华大学完整版ppt课件93§2.3差分格式的性质2.3.1范数的定义及性质1.向量范数§2.3差分格式的性质2.3.1范数的定义及性质1.942.算子范数2.算子范数95计算流体力学清华大学完整版ppt课件962.3.2差分格式的精度差分格式是微分方程的近似，通常用局部截断误差(localtruncationerror)衡量差分格式逼近微分方程的程度。2.3.2差分格式的精度差分格式是微分方程的近似，通常用97计算流体力学清华大学完整版ppt课件98计算流体力学清华大学完整版ppt课件99如果时间步长和空间步长之间满足一定的关系，FTCS格式时间方向可达到二阶精度，空间方向可达到四阶精度。根据差分格式精度的定义，按照上面的分析，FTCS格式时间方向是一阶精度，空间方向是二阶精度。如果时间步长和空间步长之间满足一定的关系，FTCS格式时间方1002.3.3差分格式的相容性截断误差是在网格点上逐点定义的。定义中每个网格点上的数值解构成一个解向量，每一个网格点上的截断误差也构成一个向量。因此，可以用向量范数来刻画差分格式的局部截断误差。2.3.3差分格式的相容性截断误差是在网格点上逐点定义的1012.3.4差分格式的收敛性和稳定性1.差分方程的矩阵形式考虑线性的发展方程（双曲型方程和抛物型方程）的差分格式。发展型方程的一般形式：以非定常热传导方程的FTCS格式为例，将差分格式写成矩阵形式：FTCS格式：解向量记为：考虑到边界条件，则差分格式可以写为：2.3.4差分格式的收敛性和稳定性1.差分方程的矩阵形1022.整体截断误差局部截断误差：差分方程逼近微分方程的程度整体截断误差：差分方程的解逼近微分方程的精确解的程度2.整体截断误差局部截断误差：差分方程逼近微分方程的程度整103计算流体力学清华大学完整版ppt课件1043.差分格式的收敛性和稳定性差分格式的收敛性对于保证数值解的有效性是非常重要的。如果差分格式是收敛的，那么，当计算网格足够密时，数值解将相当接近精确解。差分格式的稳定性等价于差分方程数值解的一致有界性。3.差分格式的收敛性和稳定性差分格式的收敛性对于保证数值解105计算流体力学清华大学完整版ppt课件106上述定理建立了算子范数的一致有界性与稳定性之间的关系。当差分格式稳定时，整体截断误差和局部截断误差量级相同。上述定理建立了算子范数的一致有界性与稳定性之间的关系。当差分107计算流体力学清华大学完整版ppt课件108Lax等价性定理是计算流体力学中的一个重要定理。直接分析差分格式的收敛性比较困难，而稳定性分析则比较简单。Lax定理告诉我们，在一定条件下，收敛性和稳定性是等价的；通过稳定性分析，即可确定差分格式的收敛条件。4.稳定性的意义Lax等价性定理是计算流体力学中的一个重要定理。直接分析差分109§2.4发展方程的稳定性分析2.4.1矩阵方法§2.4发展方程的稳定性分析2.4.1矩阵方法110计算流体力学清华大学完整版ppt课件111计算流体力学清华大学完整版ppt课件112计算流体力学清华大学完整版ppt课件113计算流体力学清华大学完整版ppt课件1142.4.2VonNeumann稳定性理论2.4.2VonNeumann稳定性理论115计算流体力学清华大学完整版ppt课件116计算流体力学清华大学完整版ppt课件117计算流体力学清华大学完整版ppt课件118计算流体力学清华大学完整版ppt课件119计算流体力学清华大学完整版ppt课件120计算流体力学清华大学完整版ppt课件121计算流体力学清华大学完整版ppt课件122计算流体力学清华大学完整版ppt课件123计算流体力学清华大学完整版ppt课件1242.4.3稳定性分析实例2.4.3稳定性分析实例125计算流体力学清华大学完整版ppt课件126计算流体力学清华大学完整版ppt课件127计算流体力学清华大学完整版ppt课件128计算流体力学清华大学完整版ppt课件129计算流体力学清华大学完整版ppt课件130计算流体力学清华大学完整版ppt课件131计算流体力学清华大学完整版ppt课件132计算流体力学清华大学完整版ppt课件133计算流体力学清华大学完整版ppt课件134计算流体力学引论TheElementsofComputationalFluidDynamics计算流体力学引论TheElementsofComput135第三章发展型模型方程的有限差分

和有限体积方法§3.1一阶线性对流方程的差分格式§3.2抛物型模型方程——对流扩散方程的

差分格式§3.3有限体积方法§3.4差分格式数值解的性质第三章发展型模型方程的有限差分§3.1一阶线性对136§3.1一阶线性对流方程的差分格式讨论双曲型模型方程：一阶线性对流方程线性对流方程的差分格式和流体力学中Euler方程的差分格式以及Navier-Stokes方程中对流项的差分格式有密切的关系，因此，掌握其差分格式的构造方法具有非常重要的意义。本节中，介绍的差分格式构造方法包括：基于导数逼近基于特征理论基于时间展开基于算子分裂§3.1一阶线性对流方程的差分格式讨论双曲型模型方程：一1373.1.1基于导数逼近的差分格式构造差分格式的最简单的方法。采用前差、后差和中心差等离散方法，直接近似微分方程中的导数项。

1.

Euler显式格式时间方向：前差。空间方向：中心差。3.1.1基于导数逼近的差分格式构造差分格式的最简单的方138

2.

Euler隐式格式时间方向：后差。空间方向：中心差。2.Euler隐式格式时间方向：后差。空间方向：139

3.蛙跳(Leap-Frog)格式时间方向：中心差分。空间方向：中心差分。3.蛙跳(Leap-Frog)格式时间方向：中心140在满足稳定性的条件时，放大因子等于1，格式具有零耗散，称为中性稳定的。

4.一阶迎风(upwind)和顺风(downwind)格式时间方向：前差。空间方向：前差或后差。Courant–Friedrichs–Lewy在满足稳定性的条件时，放大因子等于1，格式具有零耗散，称为中141计算流体力学清华大学完整版ppt课件1423.1.2基于特征线理论的差分格式，CFL条件特征性质是双曲型方程的重要特点。在构造差分格式时，考虑微分方程的数学物理性质，有助于得到性态较好的差分格式。3.1.2基于特征线理论的差分格式，CFL条件特征性质是143计算流体力学清华大学完整版ppt课件144计算流体力学清华大学完整版ppt课件145计算流体力学清华大学完整版ppt课件1463.1.3基于时间展开的差分格式3.1.3基于时间展开的差分格式147计算流体力学清华大学完整版ppt课件148计算流体力学清华大学完整版ppt课件1493.1.4基于算子分裂方法的格式3.1.4基于算子分裂方法的格式150计算流体力学清华大学完整版ppt课件151计算流体力学清华大学完整版ppt课件152计算流体力学清华大学完整版ppt课件153计算流体力学清华大学完整版ppt课件154计算流体力学清华大学完整版ppt课件155计算流体力学清华大学完整版ppt课件1563.1.5边界条件的数值处理3.1.5边界条件的数值处理157计算流体力学清华大学完整版ppt课件158§3.2抛物型模型方程—对流扩散方程的差分格式§3.2抛物型模型方程—对流扩散方程的差分格式1593.2.1求解域的离散和边界条件的处理3.2.1求解域的离散和边界条件的处理160计算流体力学清华大学完整版ppt课件1613.2.2差分格式3.2.2差分格式162计算流体力学清华大学完整版ppt课件163计算流体力学清华大学完整版ppt课件164计算流体力学清华大学完整版ppt课件1653.2.3近似因式分解方法3.2.3近似因式分解方法166计算流体力学清华大学完整版ppt课件167计算流体力学清华大学完整版ppt课件168计算流体力学清华大学完整版ppt课件169计算流体力学清华大学完整版ppt课件1703.2.4多维问题差分格式的稳定性分析3.2.4多维问题差分格式的稳定性分析171计算流体力学清华大学完整版ppt课件172§3.3有限体积方法3.3.1积分型守恒方程§3.3有限体积方法3.3.1积分型守恒方程1733.3.2空间控制体3.3.2空间控制体1743.3.3有限体积方法的全离散形式3.3.3有限体积方法的全离散形式175计算流体力学清华大学完整版ppt课件176计算流体力学清华大学完整版ppt课件177计算流体力学清华大学完整版ppt课件178计算流体力学清华大学完整版ppt课件1793.3.4有限体积方法的半离散形式3.3.4有限体积方法的半离散形式180计算流体力学清华大学完整版ppt课件181计算流体力学清华大学完整版ppt课件182计算流体力学清华大学完整版ppt课件183§3.4差分格式数值解的性质3.4.1修正方程§3.4差分格式数值解的性质3.4.1修正方程184计算流体力学清华大学完整版ppt课件185计算流体力学清华大学完整版ppt课件1863.4.2差分格式的耗散和频散3.4.2差分格式的耗散和频散187计算流体力学清华大学完整版ppt课件188计算流体力学清华大学完整版ppt课件189计算流体力学清华大学完整版ppt课件190计算流体力学清华大学完整版ppt课件191计算流体力学清华大学完整版ppt课件192计算流体力学清华大学完整版ppt课件193计算流体力学清华大学完整版ppt课件194计算流体力学清华大学完整版ppt课件195计算流体力学引论TheElementsofComputationalFluidDynamics计算流体力学引论TheElementsofComput196第五章可压缩流动数值模拟概述§5.1控制方程§5.2激波间断和广义解§5.3激波捕捉方法§5.4有限差分和有限体积方法§5.5Navier-Stokes方程中黏性项的离散§5.6时间步长的计算§5.7边界条件的处理第五章可压缩流动数值模拟概述§5.1控制方程197§5.1控制方程§5.1控制方程198计算流体力学清华大学完整版ppt课件199计算流体力学清华大学完整版ppt课件200计算流体力学清华大学完整版ppt课件201计算流体力学清华大学完整版ppt课件202计算流体力学清华大学完整版ppt课件2035.1.1守恒型的Navier-Stokes方程5.1.1守恒型的Navier-Stokes方程204计算流体力学清华大学完整版ppt课件2055.1.2守恒型Euler方程5.1.2守恒型Euler方程206§5.2激波间断和广义解5.2.1激波的形成§5.2激波间断和广义解5.2.1激波的形成207计算流体力学清华大学完整版ppt课件2085.2.2广义解5.2.2广义解209计算流体力学清华大学完整版ppt课件2105.2.3熵条件5.2.3熵条件211计算流体力学清华大学完整版ppt课件212计算流体力学清华大学完整版ppt课件213计算流体力学清华大学完整版ppt课件214§5.3激波捕捉方法§5.3激波捕捉方法2155.3.1守恒格式和Lax-Wendroff定理5.3.1守恒格式和Lax-Wendroff定理216计算流体力学清华大学完整版ppt课件217计算流体力学清华大学完整版ppt课件218计算流体力学清华大学完整版ppt课件2195.3.2人工黏性和格式黏性5.3.2人工黏性和格式黏性220计算流体力学清华大学完整版ppt课件221计算流体力学清华大学完整版ppt课件222计算流体力学清华大学完整版ppt课件223计算流体力学清华大学完整版ppt课件224§5.4有限差分方法和有限体积方法§5.4有限差分方法和有限体积方法2255.4.1有限体积方法–方案A5.4.1有限体积方法–方案A226计算流体力学清华大学完整版ppt课件227计算流体力学清华大学完整版ppt课件228计算流体力学清华大学完整版ppt课件2295.4.2有限体积方法–方案B5.4.2有限体积方法–方案B230计算流体力学清华大学完整版ppt课件231计算流体力学清华大学完整版ppt课件2325.4.3有限差分方法5.4.3有限差分方法2335.4.4有限差分方法与有限体积方法的异同5.4.4有限差分方法与有限体积方法的异同234计算流体力学清华大学完整版ppt课件235§5.5Navier-Stokes方程中黏性项的离散5.5.1Navier-Stokes方程的有限体积和有限差分格式§5.5Navier-Stokes方程中黏性项的离散5.2365.5.2黏性通量的计算方法5.5.2黏性通量的计算方法237计算流体力学清华大学完整版ppt课件238计算流体力学清华大学完整版ppt课件239计算流体力学清华大学完整版ppt课件240§5.6时间步长的计算§5.6时间步长的计算241计算流体力学清华大学完整版ppt课件242计算流体力学清华大学完整版ppt课件243计算流体力学清华大学完整版ppt课件244§5.7边界条件的处理§5.7边界条件的处理2455.7.1特征分析5.7.1特征分析246计算流体力学清华大学完整版ppt课件247计算流体力学清华大学完整版ppt课件248计算流体力学清华大学完整版ppt课件2495.7.2固壁边界5.7.2固壁边界250计算流体力学清华大学完整版ppt课件251计算流体力学清华大学完整版ppt课件2525.7.3远场边界5.7.3远场边界253计算流体力学清华大学完整版ppt课件254计算流体力学清华大学完整版ppt课件255计算流体力学清华大学完整版ppt课件2565.7.4Navier-Stokes方程的边界处理5.7.4Navier-Stokes方程的边界处理2575.7.5虚拟网格和虚拟控制体5.7.5虚拟网格和虚拟控制体258计算流体力学引论TheElementsofComputationalFluidDynamics计算流体力学引论TheElementsofComput259第六章可压缩流动的数值计算方法§6.1中心型格式§6.2迎风型格式§6.3高分辨率格式第六章可压缩流动的数值计算方法§6.1中心型格式260§6.1中心型格式§6.1中心型格式2616.1.1Lax-Wendroff格式6.1.1Lax-Wendroff格式262计算流体力学清华大学完整版ppt课件263计算流体力学清华大学完整版ppt课件264计算流体力学清华大学完整版ppt课件265计算流体力学清华大学完整版ppt课件2666.1.2MacCormack格式6.1.2MacCormack格式267计算流体力学清华大学完整版ppt课件268计算流体力学清华大学完整版ppt课件269计算流体力学清华大学完整版ppt课件270计算流体力学清华大学完整版ppt课件2716.1.3Jameson的中心型有限体积格式6.1.3Jameson的中心型有限体积格式272计算流体力学清华大学完整版ppt课件273计算流体力学清华大学完整版ppt课件274§6.2迎风型格式6.2.1一维线性波动方程组的迎风格式§6.2迎风型格式6.2.1一维线性波动方程组的迎风275计算流体力学清华大学完整版ppt课件276计算流体力学清华大学完整版ppt课件2776.2.2Euler方程的迎风型有限差分格式6.2.2Euler方程的迎风型有限差分格式278计算流体力学清华大学完整版ppt课件279计算流体力学清华大学完整版ppt课件280计算流体力学清华大学完整版ppt课件281计算流体力学清华大学完整版ppt课件282计算流体力学清华大学完整版ppt课件283计算流体力学清华大学完整版ppt课件284计算流体力学清华大学完整版ppt课件285计算流体力学清华大学完整版ppt课件2866.2.3Euler方程的迎风型有限体积格式6.2.3Euler方程的迎风型有限体积格式287计算流体力学清华大学完整版ppt课件288计算流体力学清华大学完整版ppt课件289计算流体力学清华大学完整