第三章 微分方程数值计算原理.pptVIP

(3-40) 由式（3-27），对于六点差分格式（即C-N格式）有： (3-41) 将（3-38）式、（3-39）式、(3-40)式和(3-41)式描绘成函数图形，如图3-3所示。可以看出，当f值很小时，各种差分格式很接近于精确解；但当f值逐渐增大时，不同差分格式的情况就大相径庭。在f较大，如f=1/6时，C-N格式与精确解最为接近；当f继续增大时，C-N格式与完全隐式格式比较接近精确解。 需要指出的是，上述分析只能说明由于对 项采取不同的差分近似而导致的不同差分方程的相对准确性，并不能说明一个具体差分格式的绝对精度。事实上，差分格式的准确性除取决于差分方程本身外，还可能受到边界条件差分化的影响。 图3-3 不同差分格式的精度和稳定性 §3.4直接差分法与边界条件的差分化 §3.4.1直接差分法 §3.4.2边界条件的差分化 §3.4.1直接差分法 除了基于导热微分方程用差分代替微分来建立差分方程外，还可将热传导定律直接应用于经网格剖分后的单元体，通过能量平衡列写差分方程。这种方法就是所谓的直接差分法，有时又被称为单元热平衡法。如图3-4所示的二维节点（i，j），其能量平衡为： 当λx=λy，△x=△y时，整理得： 这是二维显式差分方程，与用差分代替微分由微分方程得到的结果相同，同样用这种方法也可得到各种隐式差分格式。 直接差分法的优点有三： 首先它不涉及偏微分方程等较为复杂的数学概念及处理，物理意义直观清晰，容易理解。 其次，它能够用于不规则网格剖分，在处理复杂几何形状方面具有独特优势。 此外，它充许不同单元采用不同的热物性值，因而可方便地处理不同材料组成的系统问题。 §3.4.2 边界条件的差分化 ①给定温度的边界条件 这种边界条件的差分方程可直接写出： （n=0, 1, 2, 3, …） （3-42） 式中Tw为物体表面或边界处的温度，Tb为边界节点温度值。Tw可以是一个恒定值，也可以是一个随时间变化的函数，但无论如何，它都是已知的。 ②对流换热的边界条件 如图3-5所示，Tl，m，p为对流换热的边界节点温度，Tl-1，m，p为与其相邻的内节点温度，Tl,m,p-1，Tl,m,p+1，Tl,m+1,p，Tl,m-1,p为其相邻边界节点的温度，环境温度为T∞，对流换热系数为h。单元热平衡式如下： * 第三章 微分方程的数值计算原理 §3.1 概述 §3.2 非稳态温度场的差分格式 §3.2.1显式差分格式——时间前向差分 §3.2.2完全隐式差分格式——时间向后差分 §3.2.3六点差分格式（C-N格式）——时间中心差分 §3.2.4带权差分格式 §3.3不同差分格式的特性 §3.3.1差分格式的稳定性 §3.3.2差分格式的精度 §3.4直接差分法与边界条件的差分化 §3.4.1直接差分法 §3.4.2边界条件的差分化 §3.5 有限差分与有限元网格剖分 §3.1概述 数值计算法是一种对微分方程的近似求解 基本思想是把本来求解变量随空间、时间连续分布和变化的问题，转化为在空间领域与时间领域的有限个离散点上求变量值的问题。由于解是由一系列的数值表达出来的（解析解则是由函数式表达），故称之为数值计算法（或数值解法）。 目前主要的数值计算法有： 有限差分法（FDM） 有限元法（FEM） 一般以前者应用较多 （3-1） 有限差分法（Finite Difference Method）的物理基础仍然是各种守恒定律，其数学基础则是用差商代替微商 对于连续函数T（x），其微商（导数）为 上式中 即为差商，其中△T与△x均不为零，而 是当△x→0时差商的极限值，即微商。这里通过△x→0将 与 联系了起来。在△x未到达零之前， 始终只是 的近似值。可以把 趋于 的过程看作是由近似向精确的过渡；反之，用 代替 就意味着由精确向近似过渡。二者的差 即为用差商代替微商后所产生的误差。 将函数值