传热学编程实习_何鹏举.docVIP

传热学数值解法 学号：0811010210 姓名：何鹏举 1练习题一：一维稳态导热的数值计算 1.1物理问题 图1示出了一个等截面直肋，处于温度=80的流体中。肋表面与流体之间的对流换热系数为,肋基处温度，肋端绝热。肋片由铝合金制成，其导热系数为，肋片厚度为，高度为H=0.1m。试计算肋内的温度分布及肋的总换热量。 1.2数学描述及其解析解 引入无量纲过余温度，则无量纲温度描述的肋片导热微分方程及其边界条件如下： 其中 （其中符号含义与教科书一致） 上述数学模型的解析解为： 按（1-4）计算得到的在肋内的各点温度有表1给出。 表一 等截面直肋内各点温度 坐标x（m） 温度t（） 1.3数值离散 1.3.1区域离散 在对方程（1-1）~（1-3）进行数值离散之前，应首先进行计算区域的离散。计算区域的离散如图所示，总节点数取N。 1.3.2微分方程的离散 由于方程（1-1）在计算区域内部处处成立，因而对图1所示的各离散点亦成立。对任一借点i有： 用θ在节点i的二阶差分代替θ在节点i的二阶导数，得： 整理成迭代形式： 1.3.3边界条件离散 上面得到的离散方程式（1-5），对所有内部节点都成立，因此每个内部节点都可得出一个类似的方程。事实上，式（1-5）表达的是一个代数方程组。但这个方程组的个数少于未知数θi（i=1，2，……，N）的个数。因此，还需要根据边界条件补充两个方程后代数方程组才封闭。左边界（x=0）为第一类边界条件，温度为已知，因此可以根据式（1-2）直接补充一个方程为： 右边界为第二类边界条件，有图1中边界节点N的向后差分来代替式（1-3）中的导数，得： 将此式整理为迭代形式，得： 1.3.4最终离散格式 1.3.5代数方程组的求解及其程序 方法：高斯-赛德尔迭代方法 实施步骤：首先假定一个温度场的初始分布，即给出各点的温度初值，将这些初值带入方程组（1-6）中进行迭代计算，直至收敛。 C程序内容如下： #includestdio.h #includemath.h #define N 11 main() { int i; float cha;/*cha含义下面用到时会提到*/ float t[N],a[N],b[N]; float h,t1,t0,r,D,H,x,m,A,p; /*r代表λ,x代表Δx，D代表δ*/ printf(\t\t\t一维稳态导热问题\t\t); printf(\n\t\t\t\t\t\t----何鹏举\n); printf(\n题目：补充材料练习题一\n); printf(已知：h=45，t1=80, t0=200, r=110, D=0.01, H=0.1 (ISO)\n); /*下面根据题目赋值*/ h=45.0; t1=80.0; t0=300.0; r=110.0; D=0.01; H=0.1; x=H/（N-1）; A=3.1415926*D*D/4; p=3.1415926*D; m=sqrt((h*p)/(r*A)); /*x代表步长，p代表周长，A代表面积*/ printf(\n请首先假定一个温度场的初始分布，即给出各节点的温度初值：\n); for(i=0;iN;i++) { scanf(%f,t[i]); a[i]=(t[i]-t1)/(t0-t1); b[i]=a[i];/*这里b[i]用记录一下a[i]，后面迭代条件及二阶采用温度初场要用到*/ } /*采用一阶精度的向后差分法数值离散*/ cha=1; while(cha0.0001) { a[0]=1; for(i=1;iN-1;i++) a[i]=(a[i+1]+a[i-1])/(2+m*m*x*x); a[N-1]=a[N-2]; cha=0; for(i=0;iN;i++) cha=cha+abs(a[i]-b[i]); cha=cha/N;/*cha代表每次迭代后与上次迭代各点温度差值的平均值*/ } for(i=0;iN;i++) t[i]=a[i]*(t0-t1)+t1; printf(\n\n经数值离散(一阶精度的向后差分法)计算得肋片的温度分布为:\n); for(i=0;iN;i++) printf(%4.2f\t,t[i]); printf(\n\n); getchar(); /*采用二阶精度的元体平衡法数值离散(温度初值还用设定的初场，便于比较)*/ for(i=0;iN;i++) a[i]=b[i]; cha=1; while(cha0.0001) { a[0]=1; for