傳熱大作業(yè)兩種邊界條件匯總2404

YReportof'HeatTransfer'ComputerPractice二維導(dǎo)熱物體溫度場(chǎng)的數(shù)值模擬 IV11二維導(dǎo)熱物體溫度場(chǎng)的數(shù)值模擬假設(shè)在垂直紙面方向上冷空氣及磚墻的溫度變化很小，可以近似地予以忽略。在下列兩種情況下試計(jì)算：(1)磚墻橫截面上的溫度分布；(2)垂直于紙面方向的每米長(zhǎng)1122tt，ttt2h2tIV·上述問(wèn)題為二維矩形域內(nèi)的穩(wěn)態(tài)、無(wú)內(nèi)熱源、常物性的導(dǎo)熱問(wèn)題，對(duì)于這樣的物理??x2?y2toCt=0oC1122IVfa((m，n)xyemn的離散以網(wǎng)格線的交點(diǎn)作為需要確定溫度值的空間位置，即節(jié)點(diǎn)，節(jié)點(diǎn)的位置已該點(diǎn)在兩個(gè)t=(t+t+t+t)①對(duì)于(m，n)為邊界節(jié)點(diǎn)時(shí)：節(jié)點(diǎn)：tt=(t+2t+t)tt=(t+t)tt=(t+2t+2t+t)<m<17)時(shí)的節(jié)點(diǎn)的溫度為t=0C?！A(i,)j=0計(jì)算新的內(nèi)節(jié)點(diǎn)和邊界點(diǎn)溫否是結(jié)束IVIV·等溫邊界程序運(yùn)行結(jié)果：對(duì)流邊界程序運(yùn)行結(jié)果：·溫度分布圖·對(duì)流邊界節(jié)點(diǎn)溫度分布圖：IV·programsuiyiREAL::t(13,16),ta(13,16)REAL::dt,dtm,Q1,Q2,QREAL::LAMD=0.53REAL::epsilon=1.e-6!偏差t!迭代初值ta=0.0t(:,1)=30.0t(6,i)=0.0odoi13t(i,6)=0.0odoi!二到五行內(nèi)部節(jié)點(diǎn)溫度doj5t(i,j)=0.25*(t(i-1,j)+t(i+1,j)+t(i,j-1)+t(i,j+1))oodoi!六到十二行內(nèi)部節(jié)點(diǎn)溫度·tijtijtijtijtij+1))doi!絕熱邊界節(jié)點(diǎn)溫度計(jì)算tititit12,i))doi!絕熱邊界節(jié)點(diǎn)溫度計(jì)算titititi,15))mhenprint'writef.1))',advance='no')t(i,:)writef.1))',advance='yes')t(i,1:6)·goto1fQQtiti*LAMDQQtiti*LAMDQQtt+0.5*(t(13,1)-t(13,2)))*LAMDQ!初始內(nèi)表面導(dǎo)熱量QtitiLAMDQQtiti*LAMDQQtt))+0.5*(t(13,5)-t(13,6)))*LAMDQ=(Q1+Q2)/2.0IV·IV·programsuiyiREAL::t(13,16),ta(13,16)REAL::dt,dtm,Q1,Q2,Q,eREAL::LAMD=0.53integerhntegerhintegertoutflow=30integer::t_innerflow=10REAL::epsilon=1.e-6t=20.0ta=0.0doj5t(i,j)=0.25*(t(i-1,j)+t(i+1,j)+t(i,j-1)+t(i,j+1))oodoi=6,12!六到十二行內(nèi)部節(jié)點(diǎn)溫度ojt(i,j)=0.25*(t(i-1,j)+t(i+1,j)+t(i,j-1)+t(i,j+1))oodoi!絕熱邊界節(jié)點(diǎn)溫度計(jì)算t(13,i)=0.25*(t(13,i-1)+t(13,i+1)+2*t(12,i))odoi!絕熱邊界節(jié)點(diǎn)溫度計(jì)算·titititi,15))tihtoutflowLAMDtititi(2+0.1*h1/lamd)tjhtoutflowLAMDtjtjtj(2+0.1*h1/lamd)tjhtinnerflowLAMDtjtjtj/(2+0.1*tihtinnerflowLAMDtititi/(2+0.1*thtoutflowlamd(t(2,1)+t(1,2))/2)/(lamd+h1*0.1)thtoutflowlamdtt(lamd+h1*0.1/2)thtinnerflowlamdttlamd+h2*0.1/2)thtoutflowlamdtt(lamd+h1*0.1/2)thtinnerflowlamdtt12,6))/2)/(lamd+h2*0.1/2)tlamdtttth2*0.1*t_innerflow)/(3*lamd+h2*0.1)m·henprint'writef.1))',advance='no')t(i,:)writef.1))',advance='yes')t(i,1:6)goto1fQQtjtj*LAM