相变传热与流体流动数值分析资料报告作业1.docx

《相变传热与流体流动数值分析资料报告作业1.docx》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、相变传热与流体流动数值分析作业1学院（系）：能源与动力学院专业：能源与环境工程学生姓名：学号：指导教师：完成日期：理工大学DalianUniversityofTechnologyTheFiniteVolumeMethodforDiffusionProblemsSubjects：I.Considertheproblemofsource-freeheatconductioninaninsulatedrod.Theequationgoverningone-dimensionalsteadystateconductiveheattr

2、ansferisddxkdTdx=0,wherekequals1000W/m/K.Theendsaremaintainedatconstanttemperaturesof100℃and500℃,thecross-sectionareaAis0.01㎡.II.AlargeplateofthicknessL=2cm,thethermalconductivityk=0.5W/m/K,uniformheatgenerationq=1000kW/m3.ThefaceAandBareattemperaturesof100℃and200℃

3、respectively.ThegoverningequationisddxkdTdx+q=0.III.Thereisacylindricalfinwithuniformcross-sectionareaA.Thebaseisatatemperatureof100℃andtheendisinsulated.Thefinisexposedtoanambienttemperatureof20℃.ThegoverningequationisddxkAdTdx-hp(T-T∞)=0.Andm=hp/(kA)=25m-2,L=1m.S

4、olution://王佳琪-作业1.cpp:定义控制台应用程序的入口点。//#include#include#include#includeusingnamespacestd;#defineN5voidmain(){doublea[N-1],b[N],c[N-1],Sp[N],l[N-1],u[N],f[N],y[N],T[N];//追赶法使用的数组inti,j;//j题号doublek,L,A,TA,TB,q,m,x;//边界变量cout<<"j=

5、";cin>>j;if(j==1j==2){cout<<"k=,L=,A=,TA=,TB=,q="<>k>>L>>A>>TA>>TB>>q;x=L/N;for(i=0;i

6、c[i];elseif(i==N-1)b[i]=-Sp[i]-a[i-1];elseb[i]=-Sp[i]-c[i]-a[i-1];}for(i=0;i>m>>L>>TA>>TB;x=L/N;for(i=0;i

7、=aea[i]=-1/x,c[i]=-1/x;for(i=0;i

8、TA;}}u[0]=b[0];for(i=1;i