| 查看: 874 | 回复: 15 | |||
[交流]
陶老师数值传热学4.8节长方形通道对流换热相关问题
|
|||
|
fRe结果跟书上差不多，但是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; /*读入参数*/ printf("请输入x方向节点数N y方向节点数M\n"  ;scanf("%d%d",&N,&M); printf("请输入界面长度a 界面宽度b\n" ;scanf("%lf%lf",&a,&b); /*已知参数*/ SD=0; JF=0; dltx=a; dlty=b; De=2*dltx*dlty/(dltx+dlty); D=a;/*特征尺寸*/ N; M; Eps=1.e-6;/*计算精度*/ /*基本参数*/ 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数为:\n" ;Re=2*(De/D)*(De/D)/Wm; printf("%5.6f\n",Re); /*Nu计算*/ printf("Nu数为:\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); } 发自小木虫IOS客户端 |
回复此楼» 猜你喜欢
081700，311，求调剂
已经有15人回复
-
一志愿北京化工085600 310分求调剂
已经有18人回复
-
085600材料与化工301分求调剂院校
已经有4人回复
-
材料与化工371求调剂
已经有14人回复
-
336材料与化工085600求调剂
已经有7人回复
-
材料334求调剂
已经有18人回复
-
331求调剂
已经有8人回复
-
332求调剂
已经有17人回复
-
一志愿南京航空航天大学 材料与化工329分求调剂
已经有4人回复
-
材料专硕322
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- DNS学习记录
已经有16人回复
- 【讨论】做（计算流体力学和计算传热学）相关的请进
已经有37人回复
» 抢金币啦！回帖就可以得到:
-
江西科技师范大学物理化学课题组招生，07,08都可，还有名额
+1/95
-
2026年宁波工程学院杜俊毅副研究员招收化学材料机械类硕士研究生（含调剂）
+1/54
-
一志愿前4位代码0703均可报名-河北大学分析化学专业招收调剂考生
+1/47
-
2026-大连大学～六盘水师范学院联合培养硕士研究生，招收化学、环境等相关专业
+1/38
-
2026年武汉纺织大学材料学院邓波教授团队接收调剂（硕士生4名）
+1/37
-
广西师范大学张凤祥课题组招收5名硕士研究生
+1/36
-
煤炭科学研究总院招收2026年生态学专业硕士调剂生1名
+1/20
-
欧盟ELI项目太赫兹小组招收太赫兹光学博士
+1/19
-
吉林工程技术师范学院 交叉学科研究院 材料与化工 调剂招生
+1/16
-
广东工业大学生态环境与海洋学院张轩教授课题组招收分离膜领域2026年入学博士生1人
+1/9
-
浙江衢州，化材调剂信息（欢迎0854电气、0855机械、0856及各位符合的同学咨询哦）
+1/8
-
福建师范大学环境微生物技术课题组接受2026级硕士调剂。计划招收2-3名。
+1/7
-
湖南理工大学化学化工学院，化学、化工、材料与化工专业研究生调剂名额充足
+1/7
-
英国拉夫堡大学博士招生 电池方向
+1/5
-
齐齐哈尔大学李莉课题组诚招2026级考研调剂生（学硕和专硕）
+1/4
-
湖北师范大学招调剂啦
+1/3
-
厦门大学电磁声学研究院 陈雅鸿副教授课题组 新增两名硕士研究生名额
+1/3
-
哈尔滨工程大学动力学院赵建辉团队招收2026博士研究生
+1/3
-
招聘智慧建造或城乡规划背景硕士生
+1/2
-
北航-机械学院-招26年9月博士
+1/2
简单回复
2017-04-11 22:11
回复
我原是我以为(金币+1): 谢谢参与
这身白衣裙3楼
2017-04-11 22:18
回复
我原是我以为(金币+1): 谢谢参与
n 发自小木虫Android客户端
纳米材料7534楼
2017-04-11 22:21
回复
我原是我以为(金币+1): 谢谢参与
我原是我以为5楼
2017-04-12 19:43
回复
下载电子书6楼
2017-04-12 19:48
回复
我原是我以为(金币+1): 谢谢参与
。。。 发自小木虫Android客户端
tzynew7楼
2017-04-12 19:55
回复
我原是我以为(金币+1): 谢谢参与
2017-04-12 20:00
回复
我原是我以为(金币+1): 谢谢参与
xhmaohan9楼
2017-04-12 20:03
回复
江南的竹10楼
2017-04-12 20:08
回复
我原是我以为(金币+1): 谢谢参与
。 发自小木虫Android客户端
41588116811楼
2017-04-12 20:12
回复
hydzp12楼
2017-04-12 20:42
回复
2017-04-12 21:13
回复
我原是我以为(金币+1): 谢谢参与
skyish14楼
2017-04-12 21:21
回复
我原是我以为(金币+1): 谢谢参与
一 发自小木虫IOS客户端
guanlianwu15楼
2017-04-12 21:32
回复
我原是我以为(金币+1): 谢谢参与
假大空16楼
2017-04-12 21:59
回复
我原是我以为(金币+1): 谢谢参与
