有差异基础白色.pptVIP

有限差分法主要步骤： 利用网格将区域离散处理； 构造差分格式 用差分、差商来代替微分、微商，将微分方程离散化为差分方程，并将定解条件离散化； 解线性代数方程组 建立差分格式后，微分方程的求解就可归结为求解一个线性代数方程组，通过解线性代数方程组，得到的是数值解。;差分法的基础是用差商代替微商。 若y=f(x)是连续函数，则它的导数为;2. 差分公式;向前差公式（导数在点xi计算，而差商取fi及向前一点fi+1） ;一般地说，当差分公式的截断误差E=O(△x p)时，则称其具有p阶精度。 ;2 有限差分的基本原理;图 1定解区域网格线;(2) 差分方程的建立（差分格式的构造）;改写成便于计算的形式: ;传热分析用到的物理参数及其单位: ;3 热传导问题 ;初始条件 指某一时刻导热物体的温度分布。 对于稳定导热问题，温度场不随时间变化，时间条件自然消失。 温度随时间变化时，给出某一瞬时物体内部各点温度。 t=0时物体内部的温度分布规律通常为 T|t=0=T0(x,y,z);第二类：已知物体边界上各点的热流密度，即 ;4 稳态传热问题的有限差分方程;首先只考察内部节点。;P;上式可简化为(三维) ;二维稳态问题;(边界条件的差分形式);热流边界条件 ;每一个边界节点只应属于一种边界条件。在两种边界条件交接的节点，可人为规定属于哪一种的边界条件。;5 非稳态的有限差分方程;几何区域离散化。假定区域离散化后，距离步长Δx=xi+1-xi, Δy=yj+1-yj，且Δx=Δy。显然，xi=iΔx；yj=jΔy，i,j=0,1,2,…。 时间域离散化。用n（n=0,1,2,…）将时间区域t≥0离散化，两个时刻的间隔（时间步长）Δt=tn+1-tn，tn=nΔt。;显示差分格式 将导热微分方程应用于时刻n的节点（i, j），可写成;将三式代入式（2），得相应的差分方程为 ;tn+1时刻(i, j)节点的温度Ti,jn+1 ，可以根据自身及其相邻节点在tn时刻的温度来计算，而tn时刻的温度是已知的。因此，结合初始条件和边界条件，就可以计算区域内各节点随时间t增长的温度值Ti,jn。;或写成;等号右端用n+1时刻的一阶向后差商来近似，而等号左端温度对距离的二阶偏微商则对应tn+1时刻，故相应的差分方程为: ;上式包含邻点tn+1时刻的温度值。因此，从tn时刻的值不能简单地计算出（i, j）点tn+1时刻的温度，必须在每一个时间步长内求解一组联立方程才能求得Ti,jn+1（这组方程的数目等于待求温度的节点总数）。故称这种差分格式为隐式差分格式。 隐式差分格式多种多样，式（12）的差分形式称为完全隐式差分格式。其优点是它不受边界条件、步长的影响，是无条件稳定的格式。;对于一维热流公式为;将对应节点（i, j）的微分方程写成如下形式 ;简化为;对于二、三类边界条件，在边界外设立虚节点，使边界节点变换为内节点： 用中心差商近似一阶微商； 边界节点取内节点差分方程。;为提高整个差分格式的计算精度，常用中心差商来代替边界上的一阶微商。为此，在边界外设虚假节点。;节点1变为内节点，其显式差分方程为;同样可得换热边界条件的隐式差分格式：;内节点0显式差分格式为 ;设在x=0处为给定热流q的边界条件，且保持不变，则边界条件可写成;隐式格式;3. 求解的精确性和稳定性 ;指计算误差随步数的增加是否会积累到超出所允许的范围；或者说，最后计算结果对初始条件和边界条件的数据误差及计算中的舍入误差是否敏感的问题。 ;显式有限差分方程稳定性条件 ;在计算起始阶段，容易产生振荡，出现不稳定现象。 缩短时间步长可以提高精度。但步长如果太小，精度虽然得到保证，但计算工作量加大。 比较好的方法是采用变步长计算。在开始阶段选小步长，经短时间后，逐步加大步长。既能保证精度又可节约计算时间。;6 有限差分方程的应用;将求解区域进行网格划分。x, y方向的步长分别为Δx, Δy，节点坐标为(i, j)，i, j为整数。节点温度为Ti,j。 ;给定热流边界条件 ;;把线性代数方程组写成矩阵形式;例（2）二维非稳态问题求解 ;将的求解区域划分成正方形网格。Δx=Δy =0.03m ，节点(i, j)的温度用Ti,j表示。 ;边界条件 ;7 有限差分方程的计算机解法 ;[A]—矩阵，通常是方阵；{T}—节点温度列向量；{B}—右端列向量 ;迭代法 一般比直接法好，系数矩阵是稀疏矩阵时，迭代解法更快些，并且没有舍入误差的积累。;1. 高斯消去法;方法是：第二个方程减去第一个方程乘以a21(1)/a11(1)，第三个方程减去第一个方程乘以a31(1)/a11(1) ，直到最后一个方程。这样方程组变为等价方程组 ;其次，从等价方程的第3个方程开始，消去后面n-2