第五章--导热问题的数值解.pptxVIP

高等传热学内容;第五章导热问题的数值解;5-1导数的有限差分近似表达式;5-1导数的有限差分近似表达式;5-1导数的有限差分近似表达式;这里以二维稳态导热为例讨论有限差分法的应用。为了进行数值求解，首先要把求解的区域离散化。有限差分法要求来用正交的网格，在直角坐标系中就是矩形网格，如图5-2所示。网线的交点称为节点，在区域内的节点称为内节点，落在边界上的节点为边界节点。网线之间的距离称为步长。x方向的步长记为△x，y方向的步长记为△y。为了便于计算，需要对节点编号。若节点P的坐标为(x,y)，x＝i△x，y＝i△y，则用(i,j)表示节点的位置，节点的温度t(x，y)则记为ti,j。

数值求解的第二步是建立节点方程，即对每一个节点写出一个代数方程。建立节点方程的方法之一是从微分形式出发，也就是在微分方程或边界条件中用差商近似微商，可以得到以节点温度为未知量的代数方程。;5-2稳态导热的数值分析;有内热源的二维稳态问题，在常物性条件下可由泊松方程描述：

(5-2-1)

区域内的所有点，包括内节点(i,j)都应满足以上的方程。把节点(i,j)处的二阶偏导数用对应的差商来近似，有

(5-2-2)

(5-2-3)

把以上两式舍去截断误差并代入式(5-2-1)，就得到内节点的差分方程

(5-2-4);如果采用正方形的网格，即△x＝△y，且无内热源(qV＝0)，则式(5-2-4)简化为

(5-2-5)

注意到该差分方程的截断误差是O[(△x)2+(△y)2]。利用边界条件也可导得边界节点的差分方程。若边界节点(1,j)和(n,j)分别满足绝热边界条件和第三类边界条件：

(5-2-6)

(5-2-7)

把式(5-l-7)的朝前差分格式代入式(5-2-6)，可得边界节点(1，j)的差分方程

(5-2-8);注意到，该节点方程的截断误差为△x量级，与内节点差分方程的截断误差相比，其精度低于一个量级。为了使各节点方程的精度能够均衡，可以利用“虚节点”的概念对边界节点进行适当的处理。

假设在边界节点(l,j)的外面还有个节点(0,j)，且。注意该处的边界还是维持绝热，而节点(1,j)可以按内节点处理，得到

(5-2-9)

很显然，该节点方程的截断误差是(△x)2量级的。

对于第三类边界条件的节点(n,j)，出于同样的考虑，也可以假设有一个虚节点(n+1,j)。这样，一方面边界节点(