传热学上机报告.docVIP

传热学上机报告 学号： 姓名： 一维稳态导热的数值计算 1.1物理问题 一个等截面直肋，处于温度t∞=80的流体中。肋表面与流体之间的对流换热系数为h=45W/(m2?℃),肋基处温度tw=300℃，肋端绝热。肋片由铝合金制成，其导热系数为λ=110W/(m?℃)，肋片厚度为δ=0.01m，高度为H=0.1m。试计算肋内的温度分布及肋的总换热量。 1.2数学描述及其解析解 引入无量纲过余温度θ=t-t∞tw-t∞，则无量纲温度描述的肋片导热微分方程及其边界条件： x=0,θ=θw=1 x=H, 其中 上述数学模型的解析解为： 1.3数值离散 1.3.1区域离散 计算区域总节点数取N。 1.3.2微分方程的离散 对任一借点i有： 用θ在节点i的二阶差分代替θ在节点i的二阶导数，得： 整理成迭代形式： (i=2,3……,N-1) 1.3.3边界条件离散 补充方程为： 右边界为第二类边界条件，边界节点N的向后差分得：，将此式整理为迭代形式，得： 1.3.4最终离散格式 (i=2,3……,N-1) 1.3.5代数方程组的求解及其程序 假定一个温度场的初始发布，给出各节点的温度初值：，，….，。将这些初值代入离散格式方程组进行迭代计算，直至收敛。假设第K步迭代完成，则K+1次迭代计算式为： (i=2,3……,N-1) 程序: #includestdio.h #includemath.h #define N 11 #define p 8.1818182e-3 void main() { printf(\t一维稳态导热----------周骏\n); printf(已知：h=45，t∞=80, tw=300, λ=110, δ=0.01, H=0.1\n); double t1=300.0,t2=80.0,T[N+1]={1,1,1,1,1,1,1,1,1,1,1,1},k=0,b=0; int i,j; for(i=0;;i++) { k=T[N-1]; for(j=2;j=N-1;j++) T[j]=(T[j+1]+T[j-1])/(2+p); T[N]=T[N-1]/(1+0.5*p); if(fabs(T[N-1]-k)=1e-6) break; } printf(坐标\t\t温度\n); for(i=1;i=N;i++) printf(%.2f\t\t%.2f\n,b+=0.01,T[i]*(t1-t2)+t2); } 截图： 二维稳态导热的数值计算 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-1) 补充四个边界上的第一类边界条件得： (j=1,2,3……,M) (j=1,2,3……,M) (i=1,2,3……,N) (i=1,2,3……,N) 2.4程序及截图 程序 #includemath.h #includestdio.h #define N 5 #define M 5 void main() { printf(\t二维稳态导热----------周骏\n); printf(矩形区域，边长L=W=1，假设区域内无内热源，导热系数为常熟，三个边温度为T1=0，一个边温度为T2=1。\n); double T[N][M],m=0,x