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

登錄注冊

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

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

我原是我以為

實(shí)習(xí)版主

應(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é)點(diǎn)數(shù)N y方向節(jié)點(diǎn)數(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;/*計(jì)算精度*/

/*基本參數(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迭代計(jì)算速度*/
  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;
}
}
/*計(jì)算兩次迭代最大誤差*/
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進(jìn)行迭代求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迭代計(jì)算*/
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];
}

/*計(jì)算兩次迭代最大誤差*/
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計(jì)算*/
printf("Re數(shù)為:\n";
Re=2*(De/D)*(De/D)/Wm;
printf("%5.6f\n",Re);

/*Nu計(jì)算*/
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ù)此樓

» 猜你喜歡

26年申博已經(jīng)有4人回復(fù)
材料工程322分已經(jīng)有5人回復(fù)
318求調(diào)劑，計(jì)算材料方向已經(jīng)有7人回復(fù)
292分，材料與化工，申請調(diào)劑已經(jīng)有20人回復(fù)
085600 295分求調(diào)劑已經(jīng)有18人回復(fù)
085801 總分275 本科新能源求調(diào)劑已經(jīng)有11人回復(fù)
材料調(diào)劑已經(jīng)有9人回復(fù)
274求調(diào)劑已經(jīng)有4人回復(fù)
291求調(diào)劑已經(jīng)有24人回復(fù)
282求調(diào)劑已經(jīng)有19人回復(fù)

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

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

[考研] 085600 295分求調(diào)劑 +15	W55j 2026-03-30	18/900	2026-04-01 23:00 by 來看流星雨10
[考研] 282求調(diào)劑 +17	ycy1201 2026-04-01	19/950	2026-04-01 22:54 by peike
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +5	崔wj 2026-03-26	5/250	2026-04-01 22:24 by guanxin1001
[考研] 0703化學(xué)338求調(diào)劑！ +8	Zuhui0306 2026-03-26	9/450	2026-04-01 21:59 by xhai2011
[考研] 材料科學(xué)與工程調(diào)劑 +16	深V宿舍吧 2026-03-30	17/850	2026-04-01 21:00 by 賴春艷
[考研] 307分求調(diào)劑 +14	(o~o) 2026-03-31	15/750	2026-04-01 20:43 by longlotian
[考研] 070300化學(xué)求調(diào)劑 +13	小黃鴨寶 2026-03-30	13/650	2026-04-01 19:31 by 安靖aj
[考研] 一志愿085600中科院寧波所276分求調(diào)劑 +22	材料學(xué)257求調(diào)劑 2026-03-28	23/1150	2026-04-01 19:03 by 客爾美德
[考研] 317求調(diào)劑 +14	蛋黃咸肉粽 2026-03-26	15/750	2026-04-01 15:57 by BruceLiu320
[考研] 08工科275分求調(diào)劑 +12	AaAa7420 2026-03-31	12/600	2026-04-01 15:45 by jp9609
[考研] 求生物學(xué)調(diào)劑 +7	15172915737 2026-04-01	7/350	2026-04-01 15:27 by 閑人終南山
[考研] 085600，材料與化工321分調(diào)劑 +5	大饞小子 2026-03-27	7/350	2026-04-01 15:13 by oooqiao
[考研] 086000生物與醫(yī)藥初試274求調(diào)劑 +5	小叮當(dāng)來了 2026-03-30	6/300	2026-04-01 14:25 by mxj1203
[考研] 348求調(diào)劑 +9	zzzzyk123 2026-04-01	9/450	2026-04-01 14:24 by w蟲蟲123
[考研] 318求調(diào)劑 +8	七憶77 2026-04-01	8/400	2026-04-01 10:37 by Jaylen.
[考研] 318求調(diào)劑 +10	陳晨79 2026-03-30	10/500	2026-03-31 17:37 by 544594351
[考研] 材料求調(diào)劑一志愿哈工大總分298分，前三科223分 +11	dongfang59 2026-03-27	11/550	2026-03-31 16:51 by Wang200018
[考研] 085601 材料工程 313分求調(diào)劑 +6	Ong3 2026-03-27	6/300	2026-03-31 15:18 by yanflower7133
[考研] 一志愿大連理工大學(xué)，機(jī)械工程學(xué)碩，341 +3	西瓜田的守望者 2026-03-30	3/150	2026-03-31 11:08 by asdfzly
[考研] 2026年華南師范大學(xué)歡迎化學(xué)，化工，生物，生醫(yī)工等專業(yè)優(yōu)秀學(xué)子加入！ +3	llss0711 2026-03-28	6/300	2026-03-29 10:26 by llss0711

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

24小時熱門版塊排行榜

我原是我以為

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

» 猜你喜歡

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

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

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

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

查看全部散金貼