第一章数值计算基本概念

1、.第1章 数值计算根本概念1.1 概述关键词：CFD、微分方程à离散方程、连续解à离散点上的解1.1.1 CFD数值流体力学一般称为CFDComputational Fluid Dynamics, 为流体力学的一个重要支柱。CFD即利用离散方法discretization method, 将微分方程简化成代数方程式，通过计算机近似求解流体微分方程的方法。它的解是一些小的空间和时间上的区域上的解，称为离散点。CFD 同理论、实验并列。被人注目的理由之一是，它为计算机利用的力学计算力学的一面，特别是它为超级计算机的重要利用领域之一。此外，利用高度的图形处理，可将其结果表示非常美

2、丽的图象，对年轻人非常有魅力。因此，流体的数值模拟，在许许多多的领域内得到了利用。有许多人是在对数值计算方法理解的根底上，自己编程进展模拟。也有相当一部分人是用商用程序进展模拟。CFD包括面很广泛，从采用良好的工程设计方法，到详细求解Navier-Stokes方程；从简单流动到非常复杂的流动。简单的可能在几秒时间内就能完成，复杂的需要在最大的超级计算机上用几百个小时才能完成。完美的CFD应满足以下条件：· 适用任何问题· 计算速度快· 能得到精度高且可信度高的结果· 程序简单，谁都能简单使用· 记忆容量少其实不然，以上的要求互相矛盾，至今无一程

3、序能满足。提醒：连续介质àNavier-Stokes, 非连续介质àBoltzman1.1.2 微分方程的求解方法将连续的数据用离散的数据来记录，称为离散化discretization。在离散的点之间用光滑曲线通过内插来连接。这样，即使对于假的离散数据，只要在头脑内想象成连续的函数即可认为在对微分方程进展求解。这样，只要如今的时间和空间，就可根据这些离散数据对想象进展预测。数值流体力学的问题一般是要理解每时每刻流场的变化过程。即对支配方程式进展积分求解。实际上是求空间离散点网格上的压力、速度等物理量。图示：离散化、控制方程、压力，速度，温度光滑曲线1.2 数值求解方法的根本

4、组成关键词：数学模型、方程离散化方法、坐标、空间离散、网络、求解方法、收敛准那么1.2.1 数学模型1. 控制方程类型提醒：主导方程、支配方程根本偏微分方程的形式: 2D提醒：微分形式和积分形式 11 提醒：空间:x, y 时间空间：t,x; t,x,y对于求解域内的任一点xo, yo· 双曲型方程: , 过该点有两条实的特征线如当ac<0 异号，波动方程· 抛物型方程: , 过该点有一条实的特征线如当ac=0 ，非稳态导热· 椭圆型方程: , 过该点无实的特征线如当ac>0同号，稳态导热i. 椭圆型方程 相当于平衡问题或稳态问题。影响区域是椭圆的。与

5、时间无关。空间的闭区域。又称为边值问题。例如：稳态导热问题。稳态扩散问题。闭区域（xo,yo）求解特征：所有点联立求解。用直接法或迭代法。提示：稳态、边值、互相影响 边界ii. 抛物型方程时间步进性问题或相当于时间的步进性问题。又称为初值问题。影响区域以特征线为分界限，与主流方向垂直。例如：1D非稳态导热时间步进；2D稳态边界层型的流动和换热问题扩散忽略，主流方向步进求解特征：从的初值开场，逐步推进，依存获得适宜定边界的解。求解代数方程的量可为一维的，可节约容量。（xo,yo）前一时刻 tt+Dt物理意义：分布与瞬时以前的情况和边界条件相关。时间步进t, x推进下游的分布仅与上游的变化相关主流

6、步进xiii. 双曲型方程xo,yo也是步进问题。但依赖区域仅在两条特征区域之间。t, x例如：无粘性流体的非稳态问题；无粘性流体的稳态超音速流动。2. 流动类型偏微分方程组或积分方程组及边界条件。必须选择应用的目的：· 不可压缩ßà可压缩· 非粘性的ßà粘性· 湍流ßà层流· 2维或3维· 单相ßà多相· 。由此可以选择不同的简化守恒方程。1.2.2 控制方程的离散化方法discretization methodi. 有限差分法finite differ

7、ence method FDM 微分方程使用网络节点，选择微分的近似方法。· 将区域离散成有限个网格，通常为构造化网格；· 选择方程各项的差分形式Taylor展开；· 对每个节点建立差分方程；· 整理出关于节点上未知数的非线性代数方程式。提示：网格、微分方程、差分形式、差分方程、代数方程 ii. 有限体积法finite volume method FVM 积分方程使用控制体积，选择外表和体积积分的近似方法。· 将区域离散成有限个控制体积，适用任何形状的网格；· 选择未知函数对时间和空间的局部分布曲线线性或曲线分布；· 对每个

8、CV进展空间外表、体积和时间的积分；· 整理出关于节点上未知数的代数方程式。特点：适用任何形状的网格，可用复杂几何形状与坐标类型无关提示：网格、积分方程、分布曲线、外表和体积分、代数方程 iii. 有限单元法finite element method FE选择函数和权重函数。· 将区域离散成有限个体积或单元element，2D时通常为三角型或多边型；· 选择每个单元解的近似函数形式例如：线性形状函数，与单元角上的值相关；积分权重· 选择积分方程的权重函数；· 对每个节点值的积分残差为零，求出离散方程；· 整理出关于节点上未知数的非线性

9、代数方程式刚度矩阵。特点：有限单元法通常适用于不规那么的求解区域。提示：网格、积分方程、分布函数、权重函数、积分残差为零、刚度矩阵 iv. 频谱法spectral schemesv. 边界元法boundary element methodsvi. 分区自动化cellular automata不同的方法影响精度，求解问题的难度，编程和调试的难度，计算的速度。精度越高，涉及的网点就越多，系数矩阵就越大，需要的内存就越高，由此不得不使用粗网格，结果反而影响精度。目前一般二阶精度为最正确选择。1.2.3 坐标和根本矢量系统· Cartesian coordinate system 直角坐标系

10、统· Cylindrical coordinate system 柱坐标系统· Spherical coordinate system 球坐标系统· Curvilinear orthogonal coordinate system 曲线正交坐标系统· Non-orthogonal coordinate system 非正交坐标系统· 挪动的或静止的选择的方法依赖与目的流动。可能会影响离散方法和网格类型的选择。也可以根据矢量或张量表达的需要，选择坐标系。1.2.4 空间区域的离散化i. 计算区域domainii. 网格gridiii. 网格线gr

11、id lineiv. 格子cellv. 节点grid pointer，node, center node· 计算节点computational node, FDM· 节点FVMvi. 控制容积control volume，CVvii. 界面face1.2.5 数值网格numerical gridi. 构造化网格structured grid或称规那么网格regular grid· 网格线：自己不交，以其它线只交一次。· 节点可用一组坐标下标唯一表示， 例i,j,k· 相邻节点坐标用 ±1 表示 优点：使用广泛缺点：只适宜几何简单的计算

12、区域ii. 块构造化网格block-structured grid· 在同一个计算区域上有两种或以上不同标准的网格划分。通常使用的有粗网格、精细网格的· 粗区域可以是不规那么，可以重叠· 细网格为构造化网格细网格：123456789101112131415粗区域： I II III IViii. 非构造化网格unstructured grid· 主要用于有限体积法和有限单元法内· 格子控制体积或单元形状任意· 相邻节点数无限制常用格式形状有：2D：三角形、多边型；3D：蜂窝等，通常格子的生成有专门的格子生成方式grid generat

13、ion1.2.6 离散方程的求解方法离散化产生一个大的非线性代数方程系统。求解方法取决于问题本身。· 非稳态流动问题：使用求解初值问题的方法marching in time时间步进,在每一个时间点上，求解一个椭圆问题。· 稳态流动问题：- 准时间步进pseudo-time marching法- 等效迭代方法由于方程是非线性的，通常需要迭代。这些方法对方程使用逐次线性化，产生的线性系统几乎都是采用迭代技术来求解。提示：非线性方程系统、求解方法、问题本身、非稳态、稳态1.2.7 收敛准那么- 内迭代：求解线性方程- 外迭代：处理非线性项，和使方程耦合。何时停顿某个迭代从精度和效

14、率来说都是非常重要的。1.2.8 数值求解方法的特性提示：有解、解有界、计算收敛1. 相容性consistency 当网格跨度趋近于零时，离散差分方程接近微分方程。截断误差逐渐为零。2. 稳定性stability à不稳定: instability, 任何误差不会放大。- 暂态问题：只要真正解有解，数值解也有界。- 迭代方法：计算不发散稳定性很难判断，最常用的方法为 von Neumann方法。但是，在求解复杂的非线性的耦合方程的，并具有复杂边界条件的方法往往是很难得到稳定的结果，而需要经历和本能。许多求解方法需要限制时间步长，和采用低松弛。3. 收敛性convergence

15、4;divergence ; 收敛性ß 相容性+稳定性当网格跨度趋近于零时，离散差分方程的解接近微分方程的解。收敛与稳定同样很难判断，往往采用数值实验：逐步精化网格，如方法是稳定而且收敛的，那么结果将收敛到一个与网格大小无关的解。4. 守恒性conservationànon-conservation由于求解的方程都是守恒方程，数值结果也应是守恒的。这不仅要保证部分的守恒，也要保证总体的守恒。使用有限体积法或基于严格的守恒形式进展的离散，那么可保证每个控制体的守恒。其它离散方法那么要充分注意守恒问题。守恒问题在求解方法中是非常重要的特性。非守恒方法会导致人工源阱的产生。但非守

16、恒方法有时如采用极小的网格能保证相容性和稳定性，而产生正确的解。但一般因采用粗的网格，故建议使用守恒形式。5. 真实性realizibility对于特别复杂的情况，如湍流、燃烧、多相流动，要考虑能保证求得物理上现实的解。这不一定是数值的问题，可能是模型的问题，是否能真正的描绘物理现象。模型的问题也可能导致非物理的解或是数值方法发散。6. 精度accuracy误差分为：· 模型误差· 离散误差· 收敛误差1.2.9 流动根本方程式1. 流动守恒定理Conservation Lawsi. 连续性方程质量守恒方程, continuity equation 12 单位体积

17、的质量流束的散度净量为内部密度的增加速度。用全导数表示： 13 对于不可压缩流体， 。ii. 动量守恒方程momentum equation守恒形式： 14 左边：运动量的增加速度右边：单位体积内对流引起的动量变化、压力差作的功、粘性引起的动量变化、外力应用连续性方程，生成：非守恒形式： 15 为控制体积受到的加速度。对于Newton流体，x方向速度在y方向的变化 16 17 18 19 110 111 iii. 能量守恒方程energy equation流体总能量=内能+动能+势能 112 该形式用起来不方便，通常消去动能。将动量方程式乘于速度的矢量，为： 113 改变上式成为： 114 用

18、此式消去能量方程的动能项，使能量方程成为 115 此式表示内能的变化。可以写成温度的形式，由： 116 使用Maxwell函数，使能量方程成为： 117 q用Fourier公式： 118 对牛顿流体而言为粘性系数和耗散系数的积，耗散系数为： 119 高速气体流动时膨胀效果变的重要时、产生很大的速度梯度情况下，需考虑此0项。一般不考虑。能量方程的特殊形式：· 理想气体 120 压力一定，密度一定 连续性方程 121 · 焓的表示 à 122 · 热传导方程 123 2. Euler 方程式，Navier-Stokes 方程式通常有将包含连续方程和能量方程的

19、守恒方程式的矢量形式全体统称为Euler无粘性或 Navier-Stokes Equation的。i. 三维流动守恒方程的一般表达式 124 125 126 127 l:第二粘性系数, m:粘性系数 c:体积粘性系数 128 Stokes假定: 129 上式为可压缩流体的NS方程式的积分系形式。 不计粘性时为Euler方程式。ii. 守恒方程的微分形式 130 131 132 133 134 135 136 137 对于压缩性流体，理想气体的音速a为 138 iii. 守恒方程的积分形式 139 140 iv. 势方程Potential Equation对于超音速或低速流动流动，常使用势方程。

20、对于不存在涡流的非冲击波流动，涡度为零。可定义一个势函数F，满足： 141 v. 初值问题的特征initial value problem开展型方程：evolution equationvi. 冲击波Shock wave1.2.10 守恒方程式的守恒形式和非守恒形式1. 守恒形式 142 2. 非守恒形式 143 1.2.11 典型偏微分方程形式1. 波动方程双曲型方程双曲型方程有初值问题和初边值问题。一维非线性波动方程为最典型的2阶双曲型方程 144 线性波动方程为1阶线性双曲型方程 145 2. 扩散方程抛物线型方程一维形式： 146 3. 稳态方程椭圆型方程，Laplance方程，抛物线

21、性的稳态问题边值问题： 147 椭圆型方程为满足适宜性条件，必须给出所有边界的边界条件。1.2.12 数值流体力学的根本分类1. 不可压缩流体压力差主导的流动问题。涡流法、MAC法、拟压缩法、SIMPLE方法，有限元、边界元法。2. 可压缩流体密度波主导的流动问题。速度主导。势方程，微小扰动的势方程，波动问题，超音速流动问题，冲击波问题TVD方法， 有限元方法3. 湍流流动湍流主导问题。湍流模型：0方程1方程模型，k-e方程，壁面函数法4. 多相流流动多界面和部分平均化问题。燃烧、稀薄流、多相流、电磁流体问题。空泡率，滑速比，均匀流、2流体模型，构造方程，相间输运，状态方程、多相流动，拉格朗日

22、坐标系5. 挪动界面的流动具有明显连续挪动界面问题。自由外表，外表张力，网格系统和挪动的耦合，适体网格，波，数值波Eular法、VOF法Volume of Fluid，有限元法6. 微观模型解析详细解释微小尺度上的流动问题。粒子模型，分子模型，大涡法，格子模型7. 格式形成法复杂几何形状的问题。适体网格、坐标变换、构造网格，非构造网格的形成，适宜解形式的网格，复杂形状的网格建立1.2.13 坐标系1. 根本坐标系- 直角坐标- 圆筒坐标- 球坐标坐标系连 续 性 方 程直角坐标 圆筒坐标球坐标2. BFC适体坐标Boundary-Fitted Coordinate. 为沿着边界形状的曲线坐标系。它是通过坐标变换，将实空间上的复杂形状变换成其它空间上的简单直交坐标。一般，生成曲线坐标格子的方法有：· 利用复数变换的传统方法zàgz, z=gz, z =x+iy, z=x+ih, z平面的图形在z平面上的等角投影的坐标变换。· 利用内插函数的代数方法利用例如Lagrange等的内插多项式，求出通过所涉及的格子的曲线或与曲面垂直的曲线，将解析领域内部的格子点进展补正。这种方法对于三维问题很难。复杂领域内的，不实用。