应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 陶老师数值传热学4.8节长方形通道对流换热相关问题
161/1返回列表
查看: 930  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

我原是我以为

新虫 (小有名气)

[交流] 陶老师数值传热学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客户端
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2017-04-11 21:30:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
syhorchid2楼
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): 谢谢参与
slytjiaofei8楼
2017-04-12 20:00   回复  
我原是我以为(金币+1): 谢谢参与
xhmaohan9楼
2017-04-12 20:03   回复  
我原是我以为(金币+1): 谢谢参与
发自小木虫Android客户端
江南的竹10楼
2017-04-12 20:08   回复  
我原是我以为(金币+1): 谢谢参与
。 发自小木虫Android客户端
41588116811楼
2017-04-12 20:12   回复  
我原是我以为(金币+1): 谢谢参与
发自小木虫IOS客户端
hydzp12楼
2017-04-12 20:42   回复  
我原是我以为(金币+1): 谢谢参与
发自小木虫IOS客户端
lingyun2008101813楼
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): 谢谢参与
相关版块跳转 我要订阅楼主 我原是我以为 的主题更新
161/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6 Aj1rhIDL5ixY 2026-09-14 6/300 2026-09-17 14:17 by ObVzqQQh1rrE
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6 40ms4Wlo4umh 2026-09-14 6/300 2026-09-17 13:53 by ObVzqQQh1rrE
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 yaCh4X3wd045 2026-09-14 4/200 2026-09-17 13:17 by ObVzqQQh1rrE
[找工作] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-17 12:39 by ObVzqQQh1rrE
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +5 s3fFTmArrBt6 2026-09-13 5/250 2026-09-17 12:05 by ObVzqQQh1rrE
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 vZfe6xYu34yj 2026-09-14 4/200 2026-09-17 05:50 by k06BKNblrGRH
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 Aj1rhIDL5ixY 2026-09-14 5/250 2026-09-17 05:14 by k06BKNblrGRH
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-17 04:28 by k06BKNblrGRH
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-17 04:26 by k06BKNblrGRH
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-17 04:16 by k06BKNblrGRH
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-17 03:16 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-17 02:38 by Tql5LQhh5rLK
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-16 22:26 by Tql5LQhh5rLK
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-16 19:15 by rRvxEj3zorCy
[考博] 申博发邮件 +4 Lmengk 2026-09-15 4/200 2026-09-16 00:05 by 木心想继续深造
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:23 by CgyNCDVNhGVg
[教师之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 07:02 by 5BDX0d0WFp7t
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
信息提示
关闭
请填处理意见
关闭