有限差分法基本原理较好.pptVIP

时间导数用一阶向前差商近似代替： 空间导数用二阶中心差商近似代替： 取 ，则最终的差分方程： 显式有限差分模板： 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 100 100 0 0 0 0 0 0 0 0 0 100 100 100 100 100 100 100 100 100 100 100 100 50 50 62.5 62.5 68.8 68.8 0 25 25 37.5 37.5 45.3 0 0 12.5 12.5 21.9 21.9 0 0 0 6.25 6.25 14.1 0 0 0 0 6.25 6.25 0 0 0 6.25 6.25 14.1 0 0 12.5 12.5 21.9 21.9 0 25 25 37.5 37.5 45.3 50 50 62.5 62.5 68.8 68.8 如仍取 而为缩短计算时间，时间步长 取 ，则最终的差分方程： 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 0.0 0.5 1.0 1.5 100 100 0 0 0 0 0 0 0 0 0 100 100 100 100 100 100 100 0 200 0 100 -100 0 0 100 0 0 0 0 0 0 0 0 0 0 0 100 0 100 -100 100 0 200 有限差分法基本原理 流体的控制方程 流体的控制方程 数值离散概述 有限差分法求解流动控制方程的基本过程是：首先将求解区域划分为差分网格，用有限个网格点代替连续的求解域，将待求解的流动变量（如密度、速度等）存储在各网格点上，并将偏微分方程中的微分项用相应的差商代替，从而将偏微分方程转化为代数形式的差分方程，得到含有离散点上的有限个未知变量的差分方程组。求出该差分方程组的解，也就得到了网格点上流动变量的数值解。 离散网格点 差分和逼近误差 差分概念： 设有 的解析函数 ，函数 对 的导数为： 、 分别是函数及自变量的微分， 是函数对自变量的导数，又称微商。上式中的 、 分别称为函数及其自变量的差分， 为函数对自变量的差商。 差分的三种形式（一阶）： 向前差分 向后差分 中心差分 与其对应的差商的三种形式（一阶）： 向前差商 向后差商 中心差商 差分和逼近误差 由导数（微商）和差商的定义可知，当自变量的差分（增量）趋近于零时，就可以由差商得到导数。因此在数值计算中常用差商近似代替导数。 差分和逼近误差 差分和逼近误差 用泰勒级数展开可以推导出导数的有限差分形式。 差分和逼近误差 差分和逼近误差 逼近误差：差商与导数之间的误差，表明差商逼近导数的程度。 由函数的 Taylor 级数展开，可以得到逼近误差相对于自变量差分的量级，称为用差商代替导数的精度。 差分和逼近误差 差分和逼近误差 差分和逼近误差 差分和逼近误差 差分和逼近误差 二阶中心差分： 二阶中心差分： 差分和逼近误差 差分方程的建立过程 差分相应于微分，差商相应于导数。只不过差分和差商是用有限形式表示的，而微分和导数是以极限形式表示的。如果将微分方程中的导数用相应的差商近似代替，就可以得到有限形式的差分方程。 模型方程 为了抓住问题的实质，同时又不使讨论的问题过于复杂，常用一些简单的方程来模拟流体力学方程进行讨论分析，以阐明关于一些离散方法的概念。这些方程就叫做模型方程。常用的模型方程： 对流方程： 对流－扩散方程： 热传导方程： Poisson方程： Laplace方程： 差分方程的建立过程 以对流方程说明差分方程的建立过程。 1.划分网格 选定步长 和 ，然后在坐标平面用平行于坐标轴的两族直线划分网格： 2.针对某一点，用差商近似代替导数 对流方程在 点为 差分方程的建立过程 时间导数用一阶