流体力学第十章.ppt

第 10 章 计算流体力学简述 2. 计算流体力学发展及常用数值方法 2.1 计算流体力学的发展 一般认为计算流体力学是从20世纪60年代中后期逐步形成和发展起来的。 计算流体力学大致经历了四个发展阶段：无粘性线性、无粘性非线性、雷诺平均的N-S方程以及完全的N-S方程。 计算流体力学已经发展成为以数值手段求解流体力学物理模型、分析其流动机理为主线，包括计算机技术、计算方法、网格技术和可视化后处理技术等多种技术的综合体。 2.2 计算流体力学常用数值方法 计算流体力学不同数值方法的共同点： 离散化和代数化 计算流体力学数值方法的基本思想： 将原来连续的求解区域划分成网格或单元子区域，在其中设置有限个离散点（称为节点），将求解区域中的连续函数离散为这些节点上的函数值；通过某种数学原理，将作为控制方程的偏微分方程转化为联系节点上待求函数值之间关系的代数方程，求解所建立起来的代数方程以获得求解函数的节点值。 在计算流体力学中应用比较广泛的数值方法： 有限差分法、有限元法、边界元法、有限体积法、 有限分析法 2.3 计算流体力学主要商业软件简介 计算流体力学商业软件的三个主要部分： 前处理模块、解算模块和后处理模块 几个知名度较高的大型通用商业软件： CFX FLUENT PHOENICS ANSYS STAR-CD 3.2 差分格式及其基本构造方法 1. 差商、逼近误差 设有x的解析函数y = f (x) ，则 称为函数y = f（x）在点x的差商，而 、 分别称为函及自变量的差分 称 表示的差商称为一阶向前差商， 类似地，称 一阶向后差商； 称 为一阶中心差商。 而称 二阶中心差商。 从以上的差商表达式可知，所取差商不同，逼近误差也不同。对于一阶导数来说，向前差商和向后差商的逼近误差相同，都是 量级，具一阶精度；而一阶和二阶中心差商的逼近误差同是（ ）2量级，都具有二阶精度。对于一阶导数的差商逼近，中心差商的精度比向前差商和向后差商都要高。 2. 差分网格、差分方程 【例11-1】一维对流方程初值问题的差分格式构建。微分方程和初值条件为 图10-1 网格剖分 在上述的网格剖分下，可写出对流方程定解问题的差分格式如下： 像上式这样的时间导数采用向前差商近似，空间导数采用中心差商近似的格式称为FTCS（Forward Time Central Space）格式。 上述一维对流方程初值问题也可采用时间和空间都向前差分的格式来近似，这时差分方程成为: 像上式这样的时间导数和空间导数都采用向前差商近似的格式称为FTFS（Forward Time Forward Space）格式。 若采用时间向前差分、空间向后差分，则可得到一维对流方程的FTBS（Forward Time Backward Space）差分格式： 3.3 差分方程的相容性、收敛性和稳定性 采用有限差分法求解由微分方程所描述的流体流动问题时，在确定差分离散格式是否可用之前必须回答三个问题： 当差分网格的时间与空间步长趋向于零时，差分方程是否充分逼近原微分方程； 差分格式的真解是否充分逼近原微分方程的精确解； 差分格式的近似解与真解之间的误差是否有界。 这三个问题在有限差分法理论中分别称为相容性，收敛性和稳定性。 相容性 所谓相容性是指当差分网格步长（时间步长和空间步长）都趋向于零时，差分方程的截断误差也趋向于零，表明差分方程在形式上逼近原微分方程。 若在每一网格节点上差分方程与原偏微分方程都是等同的，则称差分方程与原微分方程相容。 若网格步长趋向于零时，差分方程的截断误差不趋向于零，导致差分方程在形式上不逼近于原微分方程，则称差分方程与原微分方程不相容。 2. 收敛性 对于差分网格上的任意计算节点，它对应于原微分方程定解域上的点 ，设差分格式在此点的精确解为 ，而原微分方程的精确解为u ，则二者之差： 称为差分格式在此网格节点上的离散化误差。 如果在任一网格节点上，当 、 时，都有差分格式在该点的离散化误差 ，则此差分格式的解将收敛于原微分程的解。我们称此差分格式是收敛的。 3. 稳定性 利用有限差分格式数值求解偏微分方程是按时间逐层推进计算的。在这样逐层步进的数值计算过程中，误差不可避免。在求解的逐层推进过程中，必须分析这些误差传播的情况——如果这些误差