传热学上机C程序源答案之一维非稳态导热的数值计算(共4页).docx

精选优质文档-倾情为你奉上二维稳态导热的数值计算2.1物理问题一矩形区域，其边长L=W=1，假设区域内无内热源，导热系数为常数，三个边温度为T1=0，一个边温度为T2=1，求该矩形区域内的温度分布。2.2 数学描述对上述问题的微分方程及其边界条件为： x=0，T=T1=0 x=1，T=T1=0 y=0，T=T1=0 y=1，T=T2=1该问题的解析解：2.3数值离散2.3.1区域离散 区域离散x方向总节点数为N，y方向总节点数为M，区域内任一节点用I,j表示。2.3.2方程的离散对于图中所有的内部节点方程可写为：用I,j节点的二阶中心差分代替上式中的二阶导数，得：上式整理成迭代形式：(i=2,3,N-1)，(j=2,3,M-