版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

我原是我以為

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 343.8
帖子: 163
在線: 84.1小時
蟲號: 4254249

[交流] 陶老師數(shù)值傳熱學(xué)4.8節(jié)長方形通道對流換熱相關(guān)問題

fRe結(jié)果跟書上差不多，但是Nu差別好大，有大神幫我看看嘛

#include <math.h>
#include <stdio.h>
#define NT 70
void main()
{
/*定義變量*/
int N,M,i,j,Iter;
double W[NT][NT],W0[NT][NT],sita[NT][NT],sita0[NT][NT],T[NT][NT],De,Re,Nu;
double dltx,dlty,dx,dy,Wm,D,x,y,a,b;
double Eps,DT,DTmax,SD;
double ap,aw,ae,as,an,JF,b1,b2,JA,ap1;
FILE*fp;

/*讀入?yún)?shù)*/
printf("請輸入x方向節(jié)點數(shù)N y方向節(jié)點數(shù)M\n"

;
scanf("%d%d",&N,&M);
printf("請輸入界面長度a 界面寬度b\n"

;
scanf("%lf%lf",&a,&b);

/*已知參數(shù)*/
SD=0;
JF=0;
dltx=a;
dlty=b;
De=2*dltx*dlty/(dltx+dlty);
D=a;/*特征尺寸*/
N;
M;
Eps=1.e-6;/*計算精度*/

/*基本參數(shù)*/
dx=dltx/(N-1);
dy=dlty/(M-1);
ae=aw=dy/dx;
an=as=dx/dy;

/*速度邊界條件處理*/
/*固體表面速度為0值*/
for(j=0;j<M;j++)
  {
W[0][j]=0;
W[N-1][j]=0;
}
for(i=1;i<N-1;i++)
{
W[0]=0;
W[M-1]=0;
}
/*速度場賦初值*/
for(j=1;j<M-1;j++)
  for(i=1;i<N-1;i++)
{
W[j]=0.00008;
W0[j]=W[j];/*對上一次迭代值賦值*/
}
Iter=0;
/*Gauss-Seidel迭代計算速度*/
  loop1:
for(j=1;j<M-1;j++)
{
for(i=1;i<N-1;i++)
  {
  ap=ae+aw+an+as;
b1=(dx*dy)/(D*D);
W[j]=ae/ap*W[i+1][j]+aw/ap*W[i-1][j]+an/ap*W[j+1]+as/ap*W[j-1]+b1/ap;
}
}
/*計算兩次迭代最大誤差*/
Iter=Iter+1;
DTmax=0.0;
for(j=1;j<M-1;j++)
  for(i=1;i<N-1;i++)
{
DT=fabs(W[j]-W0[j]);
W0[j]=W[j];
if(DT>DTmax)
DTmax=DT;
}
if(DTmax>Eps)
goto loop1;

for(j=1;j<M-1;j++)
for(i=1;i<N-1;i++)
{
SD+=W[j]*dx*dy;
}
Wm=SD/(a*b);
printf("無量綱速度平均值:\n";
printf("Wm=%5.6f\n",Wm);

/*邊界條件處理*/
/*固體表面無量綱為0,取T=JA*sita進行迭代求JA*/
for(j=0;j<M;j++)
  {
sita[0][j]=0;
sita[N-1][j]=0;
}
for(i=1;i<N-1;i++)
{
sita[0]=0;
sita[M-1]=0;
}
/*sita場賦初值*/
for(j=1;j<M-1;j++)
  for(i=1;i<N-1;i++)
{
sita[j]=0.003;
sita0[j]=sita[j];/*對上一次迭代值賦值*/
}

Iter=0;
/*Gauss-Seidel迭代計算*/
loop2:
for(j=1;j<M-1;j++)
{
for(i=1;i<N-1;i++)
{
JF+=W[j]*sita[j];
}
}
JA=a*b/JF/Wm;
for(j=1;j<M-1;j++)
for(i=1;i<N-1;i++)
{
b2=dx*dy*JA*W[j]/(D*D*Wm);
ap1=ap-b2;
sita[j]=ae/ap1*sita[i+1][j]+aw/ap1*sita[i-1][j]+an/ap1*sita[j+1]+as/ap1*sita[j-1];
}

/*計算兩次迭代最大誤差*/
Iter=Iter+1;
DTmax=0.0;
for(j=1;j<M-1;j++)
  for(i=1;i<N-1;i++)
{
DT=fabs(sita[j]-sita0[j]);
sita0[j]=sita[j];
if(DT>DTmax)
DTmax=DT;
}
if(DTmax>Eps)
goto loop2;

/*Re計算*/
printf("Re數(shù)為:\n";
Re=2*(De/D)*(De/D)/Wm;
printf("%5.6f\n",Re);

/*Nu計算*/
printf("Nu數(shù)為:\n";
Nu=0.25*(De/D)*(De/D)*JA;
printf("%5.6f\n",Nu);

for(j=1;j<M-1;j++)
  for(i=1;i<N-1;i++)
  {
   T[j]=JA*sita[j];
  }

  fp=fopen("result.dat","w";
  printf("無量綱溫度:\n";
  for(j=M-1;j>=0;j--)
  {
   for(i=0;i<N;i++)
   {
   fprintf(fp,"%5.6f",T[j]);
   }
   printf("\n";
  }
  fclose(fp);
}

發(fā)自小木蟲IOS客戶端

回復(fù)此樓

» 猜你喜歡

求調(diào)劑已經(jīng)有3人回復(fù)
265求調(diào)劑已經(jīng)有4人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有6人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有12人回復(fù)
286分人工智能專業(yè)請求調(diào)劑愿意跨考！已經(jīng)有3人回復(fù)
329求調(diào)劑已經(jīng)有5人回復(fù)
申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復(fù) 已經(jīng)有4人回復(fù)
材料學(xué)碩318求調(diào)劑已經(jīng)有5人回復(fù)
一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑已經(jīng)有6人回復(fù)
081700化工學(xué)碩調(diào)劑已經(jīng)有3人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

DNS學(xué)習(xí)記錄已經(jīng)有16人回復(fù)
【討論】做（計算流體力學(xué)和計算傳熱學(xué)）相關(guān)的請進已經(jīng)有37人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2017-04-11 21:30:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

簡單回復(fù)

syhorchid2樓

2017-04-11 22:11 回復(fù)

我原是我以為(金幣+1): 謝謝參與

這身白衣裙3樓

2017-04-11 22:18 回復(fù)

我原是我以為(金幣+1): 謝謝參與

n 發(fā)自小木蟲Android客戶端

納米材料7534樓

2017-04-11 22:21 回復(fù)

我原是我以為(金幣+1): 謝謝參與

我原是我以為5樓

2017-04-12 19:43 回復(fù)

發(fā)自小木蟲IOS客戶端

下載電子書6樓

2017-04-12 19:48 回復(fù)

我原是我以為(金幣+1): 謝謝參與

。。。發(fā)自小木蟲Android客戶端

tzynew7樓

2017-04-12 19:55 回復(fù)

我原是我以為(金幣+1): 謝謝參與

slytjiaofei8樓

2017-04-12 20:00 回復(fù)

我原是我以為(金幣+1): 謝謝參與

xhmaohan9樓

2017-04-12 20:03 回復(fù)

我原是我以為(金幣+1): 謝謝參與

發(fā)自小木蟲Android客戶端

江南的竹10樓

2017-04-12 20:08 回復(fù)

我原是我以為(金幣+1): 謝謝參與

。發(fā)自小木蟲Android客戶端

41588116811樓

2017-04-12 20:12 回復(fù)

我原是我以為(金幣+1): 謝謝參與

發(fā)自小木蟲IOS客戶端

hydzp12樓

2017-04-12 20:42 回復(fù)

我原是我以為(金幣+1): 謝謝參與

發(fā)自小木蟲IOS客戶端

lingyun2008101813樓

2017-04-12 21:13 回復(fù)

我原是我以為(金幣+1): 謝謝參與

skyish14樓

2017-04-12 21:21 回復(fù)

我原是我以為(金幣+1): 謝謝參與

一發(fā)自小木蟲IOS客戶端

guanlianwu15樓

2017-04-12 21:32 回復(fù)

我原是我以為(金幣+1): 謝謝參與

假大空16樓

2017-04-12 21:59 回復(fù)

我原是我以為(金幣+1): 謝謝參與

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主我原是我以為的主題更新

返回列表

普通表情龍兔虎貓高級回復(fù) (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 294求調(diào)劑材料與化工專碩 +14	陌の森林 2026-03-18	14/700	2026-03-19 22:38 by 學(xué)員8dgXkO
[考研] 0703化學(xué)調(diào)劑，六級已過，有科研經(jīng)歷 +12	曦熙兮 2026-03-15	12/600	2026-03-19 19:42 by maocaozhuxi
[考研] 317求調(diào)劑 +3	申子申申 2026-03-19	6/300	2026-03-19 14:16 by 申子申申
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7	小材化本科 2026-03-18	7/350	2026-03-19 10:46 by Linda Hu
[考研] 274求調(diào)劑 +6	S.H1 2026-03-18	6/300	2026-03-19 09:34 by 花店相見
[考研] 354求調(diào)劑 +4	Tyoumou 2026-03-18	7/350	2026-03-18 21:45 by Tyoumou
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +7	步川酷紫123 2026-03-13	7/350	2026-03-18 17:12 by 盡舜堯1
[考研] 08工科 320總分求調(diào)劑 +5	梨花珞晚風(fēng) 2026-03-17	5/250	2026-03-18 14:49 by haxia
[考研] 311求調(diào)劑 +11	冬十三 2026-03-15	12/600	2026-03-18 14:36 by 星空星月
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 303求調(diào)劑 +4	睿08 2026-03-17	6/300	2026-03-18 11:01 by Iveryant
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.
[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 22408總分284求調(diào)劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425
[考研] 中科大材料與化工319求調(diào)劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3	麥田耕者 2026-03-14	3/150	2026-03-14 20:04 by 外星文明
[考研] 304求調(diào)劑 +7	7712b 2026-03-13	7/350	2026-03-13 21:42 by peike

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

我原是我以為

[交流] 陶老師數(shù)值傳熱學(xué)4.8節(jié)長方形通道對流換熱相關(guān)問題

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

» 搶金幣啦！回帖就可以得到: 查看全部散金貼

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

» 搶金幣啦！回帖就可以得到:

查看全部散金貼