python实现数值传热模拟.doc

题目： 流体以均匀速度 u 在半径为 R，长度为 L 圆管中流动（L/R=4），忽略径向温度变化。假定在进口处流体的温度为 TL, 出口处温度为 T0，管内壁温度也为 T0。 其中：无量纲温度φ =(t-T0)/(TL-T0)； 无量纲坐标 X=x/R； 边界条件： 【X=0 时 , φ =0;X=1 时 , φ =1 ?R/λ =1,uR/α =1】 问题：当节点个数为4个时，求X=1，2，3时，各 节 点无量纲 温 度φ值， 并 绘 图 表 示 。 并且，依次求解当节点个数n为8，16，32，64时各节点的φ值； n=8时,求X=1,2，3……6，7，8；时8个节点的φ值； n=16时,求X=1,2，3……14，15，16；时16个节点的φ值； n=32时,求X=1,2,3,……31，32，33；时32个节点的φ值； n=64时,求X=1,2,3……62,63,64；时64个节点的φ值； 然后比较，网格划分粗细程度对数值解的影响。（其中解析解为 ?(X)=eX-e-2Xe4- 平衡关系式： ①控制体内热量的增量 =② 由于热对流进入控制体的净热量+③ 由于导热进入控制体的净热量+④ 流体与管壁间的对流换热进入控制体净热量其中：①稳态：①=0；②热对流② =(??×??×??×????2)[(??i+1+ ??i)/2 - (??i-1+ ??i)/2] = (??×??×??×????2) (??i+1+-??i-1)/2;③导热③ = ????2×λ(??i+1-??i/??? + (??i -1- ??i)/???] = ????2×λ (??i+1-2??i +??i -1)/???;④对流换热④ = 2????×???×? (??0 - ??i); 代入平衡关系式： (??×??×??×????2) (??i+1-??i-1)/2+ ????2×λ (??i+1-2??i +??i -1)/??? +2????×???×(??0 - ??i)=0 同时除以???，得： (??×??×??×????2) (??i+1-??i-1)/2??? + ????2×λ (??i+1-2??i +??i -1)/(???)2 +2????×? (??0 - ??i) =0无量纲化控制方程得： (??×??×??×????) (φi+1-φi-1)/2 X + ??×λ (φi+1-2φi +φi -1)/X2 -2??R×?φi=0又根据?R/λ =1,uR/α =1，得： (φi+1-φi-1)/(2X) + (φi+1-2φi +φi -1)/X2-2φi =0整理得：-(2-X) φi -1+(4+4X2) φi - (X+2)φi+1=0 三节点关系式为： -aiφi-1+biφi-ciφi+1=di (1) 于是有 a1 当节点数 为n时, X =4/n； ai=2-4/n ; bi=4+64/n 2; ci=（4/n)+2; di=0； 递推关系式：φi=Giφ i+1+Qi (2) φi-1=Gi-1φ i+Qi-1 (3) 联立(1)、(3)式，并与得(2)对比得： Gi=ci/（bi-aiGi-1） Qi=（di+aiQi-1）/（bi-aiGi-1） 则有，G1= ci/bi ,GN=0；Q1= d1/b1； 则可依次求得G1 ，Q1；G2 ，Q2；……Gn ，Qn； 由φn=1,代入(3)后可依次求得φn-1，φn-2，φn-3……φ1。 当总节点数为n=4，8，16，32，64时，使用 python 编程软件，写出代码如下:#导入画图库import matplotlib.pyplot as plt#定义函数，以所需要的节点数作为变量，函数返回 Q_S(无量纲温度)数组def Numerical_heat(n):#定义 5 个变量的空列表，以下将公式计算出的结果依次添加到列表中a,b,c,d,G；Q_S=[],[],[],[],[]#构造 a 列表for i in range(1,n):if i==1:aa=0else:aa=-(n/8)+(n*n/16)a.append(aa)#构造 b 列表for i in range(1,n):3if i==1:bb=2else:bb=(n*n/8+2)b.append(bb)#构造 c 列表for i in range(1,n):if i==1:c