计算流体力学一维稳态导热编程作业.docx
TheFiniteVolumeMethodforOne-DimensionalDiffusionProblemsLProblemof1-dimensionalSteady-StateSource-freeHeaiConductionConsidertheproblemofsource-freeheatconductioninaninsulatedrodwhoseendsaremaintainedatconstanttemperaturesof100Cand500Crespectively.Theone-dimensionalproblemsketchedintheFigure1isgovernedbyCalculatethesteadystatetemperaturedistributionintherod.Thermalconductivitykequals1000WmK,cross-sectionalareaAis10"2n2.Fig.1PhysicalModel1网格划分条件L=0.5m,x=£/5=O.lw,女=IoOOW/"ZK,Ta=100,Tb=500A=W2W22方程离散求解域内共有5个节点，节点2、3、4的离散方程：告“+偿f7;=0G+信-“工+,xPE)VxWP)IxPE)I函由于Z(I=kw=k,BXPe=xwp=x,Ae=Aw=A均为常数,方程：a.JT.=a.rT.r+a.T.工aJoCWP/因此对节点2、3、4有离散5”pSWrV'SE'E.k.k.式中。卬=""A>aE""Acip=y+aE，QXQX节点1的离散方程：keAe-kwAw与立QXPExAP袅斗FAj7;=Oz+停XxPE)xP)3&可写为：apTp=awTw+aETE+Sltk了LI°eAe，。卬0，apa”,十SpxPE(2k、S”-也4,S=AvvTxapIAP>同理，节点5的离散方程%Ae生二里-&A,.OXpb(冬Aj&AjO=0Z+aLI瓯PI闲SJ(瓯/可写为：apTp=awTw+aETE+S11k其I4七一0,SAV，ap-aw+aESPxWP=0-ate+(ataI电P)4ZZ2k=0xW'P-AJ%+AlZ)所以得到各节点的“w,ciE,op,Sc,SP的值各节点离散方程系数节点、awaSPSUQp=aw+ESp10Aa×CkA-2Ax29北OXCkA32Aa<£vA002rA3Axx002巨4x4Aax-A加002Ax5Ax0-2AxZR”OX3Ax即:节awaSPSUap=aw+aESp10100-2002007；30021001000020031001000020041001000020051000-200200TB300从而得下述代数方程组3007；=1007；+2007；2007；=1007；+1007；<2007；=1007；+1007；2007；=1007；+1007；3007；=1007；+2007；写成矩阵形式有：300-100000'Tl'-200-100200-10000T200-100200-1000q=000-100200-1000000-100300Z.2007；将7；=1007尸500代入，解得此方程组为:3程序C程序内容如下：#include"stdafx.h,'#include<stdio.h>#defineN5intmain()inti;doubleGL,TA,TB,dx,k,Area,Lenth;doubleAWN+1,AEN+1jAPN+1,SPN+1jSUN+1;doubleaN+l,bN+l,cN+l,fN+ljMN+1jLN+1,UN+1jYN+1>TN+1;Lenth=0.5;TA=100;长度西边界温度TB=500;东边界温度d×=LenthN;k=1000;网格宽度导热系数Area=0.01;截面积第一点计算AW1=0.0;AEl=k*Aread×SPl=-2*k*Areadx;SUl=2*k*Areadx*TA;/最后一点计算AWN=k*Areadx;AEN=0.0;SPN=-2*k*Areadx;SUN=2*k*Areadx*TB;/中间点计算for(i=2ji<=N-lji+)AWi=k*Aread×AEi=k*Aread×SPi=0.0;SUi=0.0;计算Apfor(i=l;i<=N;i+)APi=AWi+AEi-SPi;三对角方程组追赶法求解for(i=lji<=Nji+)ai+l=-AEi;bi=APi;ci=-AWi+l;fi=SUi;)for(i=l;i<=N-l;i+)i=ci;Ll=bl;f0r(i=2;i<=N;i+)Mi=aiLi-l;Li=bi-Mi*Ui-lj)Yl=fljfor(i=2;i<=N;i+)Yi=fi-Mi*Yi-l;TN=YNLN;for(i=N-l;i>=l;i-)Ti=(Yi-Ui*Ti+l)Li;Printf(',方程组ax=b的解为：n");for(i=l;i<=N;i+)printf("T%d=%.3fn"ji,Ti);return0;
