数值传热学第三章扩散方程数值解法.pptVIP

第3章 扩散方程的数值解法及应用 * 第3章 扩散方程的数值解法及应用 3.1 一维导热问题 3.2多维非稳态导热问题 3.3 源项及边界条件处理 3.4离散方程求解方法（LSOR） 3.5管内充分发展流动对流换热数值计算示例 3.5.1充分发展流定义 3.5.2纵向内肋管中的充分发展对流换热 3.5.3矩形通道内的充分发展对流换热 3.1一维导热问题 扩散方程包括了一类由于分子热运动现象，例如导热、扩散和多孔介质中的渗流等现象。 3.1.1控制方程 其中，A为面积因子，λ为导热系数，S为源项，x为空间坐标变量。对于直角坐标系，A=1；对于柱坐标系，A=r；而对于球坐标系A=r2。 3.1.2 方程离散 两边同时乘以A，在控制容积上积分得： 上式假设源项可分成两部分，即 整理得 到此为止，只剩下界面物性λe，λw取值问题了。 3.1.3界面扩导计算 1.算术平均 2.调和平均（热阻串联） 从节点P到节点E的热阻可看成由P到e和由e到E的两部分。 调和平均与算术平均对于非线性非常弱或常扩导问题，二者是一样的。但对于扩导阶跃的情况，算术平均的结果与物理概念相差很远。此时，若把控制容积的界面放到物性阶跃面上，调和平均的物理意义非常明显，且可得到合理的结果。 扩散（导热）方程离散示例 有一锥状杆，x=0处，直径d=d1,温度w为t1;x=L处，d=d2 , t=2,周围流体温度为tf，杆侧面换热系数为h,求杆内温度分布。 首先把区域离散N-1个控制容积，控制容积界面坐标为xu(i),i=1,2,…,N-1,控制容积代表点的坐标为x(i) , i=1,2,…,N. 本问题中源项为： 离散后得 本例中可取 由于源项非线性，需要迭代求解。 从i=2到N-1，扫描所有控制容积，得线性代数方程组： 此方程组共有N-2个未知量，N-2个方程。可用迭代法求解。求解时需要注意此处迭代有两重含义（1）线性代数方程组的迭代问题；（2）非线性迭代（要更新常数项b)问题。 由于两端都是第一类边界条件，故T1＝t1-tf，TN=t2-tf是已知的，只需求内部点的过余温度值。 若导热系数?＝5，h=12 ,tf=20 , t1=150 , t2=30 , L=0.5 , d1=0.025 , d2=0.01,取N=10,50，100分别求解。 1、区域离散，根据给定的L和N，形成控制容积界面数组，xu,和控制容积节点坐标x。要求控制容积非均分，用一个子GRID 程序实现。 3、计算各控制容积的源项Sp,Sc，子程序SOURCE； 4、形成线性代数方程组,aP,aW,aE,b, 子程序COEF; 5、求解线性代数方程组，子程序SOLVE 程序结构成分 2、计算截面积Ap(x), 子程序GEOG 程序框图步骤 1、CALL GRID 2、CALL GEOG 3、赋初始温度场，T=T0=guessed value。 4、CALL SOURCE 5、CALL COEF 6、CALL SOLVE, Get T new values 0、输入已知参数（问题描述参数） 3.2多维非稳态导热问题 有了前面一维稳态导热问题数值解方法的铺垫，多维非稳态导热问题就可迎刃而解了。 一、控制方程 取如图所示的控制容积，离散包括两个空间坐标离散和时间坐标离散。它们都按照控制容积积分法离散。 二、隐式格式 隐式与显式的物理介释：不稳态导热问题，某点温度的变化是因为周围点导热不均衡所致。对于二维问题，周围点简化为四个邻点（NEWS），P点的温度变化是因为NEWS与P导热量不平衡。由于P点和NEWS四点的温度是变化的，在计算传热时，精确的讲应该按时间积分，但由于NEWS及P点的温度对时间的变化关系是未知的，积分不可行，因此必须作出适当的假设。这样就导致了所谓隐格式、显格式和Crank-Nicolson格式。 显然，等号右边的表达式有许多种计算方法，其中常用的有三种。即nt=t，表示传入控制容积的热量温度势差在Δt内不变。nt=t+ Δt，认为传入控制容积的热量取决于t+ Δt时刻的温度差且在Δt内不变。另外一种是取两个时刻温度差的平均值。 显然，当等号右边取t时刻的值时，每个方程中只有一个未知量，可采用递推的方式逐点计算，也就是说各点的温度值不必关联求解，比较简单，但对时间步长有限制。取t+ Δt时刻的值时则不同，每个方程中的未知量除本点温度外，还与邻点的温度有关，需要连立求解。Crank—Nicolson格式也要联立求解。无论从数学还是从物理概念上讲，C-N格式都应该是优秀的，但数值计算表明，它会出现振荡。 其他坐标系也可建立相应的离散方程，总结如下表。 3.3源项及边界条件处理 一、源项线性化（Linearization of source term) For a gi