數(shù)值傳熱Fortran程序_第1頁
數(shù)值傳熱Fortran程序_第2頁
數(shù)值傳熱Fortran程序_第3頁
全文預(yù)覽已結(jié)束

下載本文檔

版權(quán)說明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡介

班級：2014級暖通學(xué)碩班學(xué)號：2014214045姓名：李昊民題目：流體以均勻速度在半徑為R，長度為L的圓管內(nèi)流動(dòng)，忽略其徑向溫度變化，假定在進(jìn)口處流體的溫度為TL,出口處流體溫度為TO，管壁溫度也是TO，無量綱溫度Φ=t-TLTL-To,無量綱坐標(biāo)X=xR,X=0時(shí)，ΦO=0，解：選控制體←μ建立熱平衡方程由于對流進(jìn)入控制體的熱量為：（1）由于導(dǎo)熱進(jìn)入控制體的熱量為：（2）由于管壁對流換熱進(jìn)入控制體的熱量為：（3）由(1)+(2)+(3)=0總的平衡方程為：無量綱化及解方程劃分控制體的個(gè)數(shù)設(shè)為n，則無量綱化解得結(jié)果為：（1+2n）Φi+1-(2+32n2)Φi+1程序編制采用fortran編寫程序，變成如下:programshuzhichuanreimplicitnoneintegern,i,krealmreal,allocatable::a(:),b(:),c(:),g(:),o(:),q(:)write(*,*)'請輸入'read(*,*)nallocate(a(n),b(n),c(n),g(n),q(n),o(n))m=ng(1)=(1+2/m)/(2+32/(m**2))o(n)=1doi=1,na(i)=1-2/mb(i)=2+32/(m**2)c(i)=1+2/menddodoi=1,n-1g(i+1)=c(i)/(b(i)-a(i)*g(i))enddowrite(*,*)g(3),a(n),b(n),c(n)k=n100o(k-1)=g(k-1)*o(k)k=k-1 if(k>1)then

溫馨提示

  • 1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽，若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間，僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理，對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論