版块导航: 正在加载中...

客户端APP下载

调剂小程序

登录注册

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

24小时热门版块排行榜

SNOWANGELLEI

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 5
帖子: 5
在线: 1.4小时
虫号: 3342568
注册: 2014-07-28

[求助] 急需大神指点FLUENT中UDF问题！

在UDF中编写如下程序，想输出每次迭代后的COP值，但是总会出现迭代错误，我是新人，求教各位大神！拜托拜托！
#include"udf.h"
#include"stdio.h"

real Q[8760];
int i=-2;
int k=1;

real cop[8760];
FILE *Tem_t;
FILE *fp;

DEFINE_ON_DEMAND(readfile)
{
  int j;
  float t1,t2;
  Tem_t=fopen("load.txt","rb"

;
  for (j=0;j<8760;j++)
{
  fscanf(Tem_t,"%f",&t1);
  Q[j]=t1;
}

}

DEFINE_PROFILE(Tout1,t, position)
{
face_t inlet,outlet;
Thread *thread;
real temp1,cop;
int ID=16;
Domain *domain;
domain=Get_Domain(1);

begin_f_loop(inlet,t)
{
   thread=Lookup_Thread(domain,ID);

      begin_f_loop(outlet,thread)
   {
         temp1=F_T(outlet,thread);
   }
      end_f_loop(outlet,thread)
  cop[k]=5*(-0.9830*(0.94395*pow((temp1-303)/23,2)+0.8683*(temp1-303)/23+1)+1.9709);

if(Q[i/2]>0) /*如果满足条件*/
{
   F_PROFILE(inlet,t,position)=temp1+Q[i/2]*(1+1/cop[k])/4.2/100/2/0.278;
}
else/*否则*/
{
   F_PROFILE(inlet,t,position)=temp1+Q[i/2]*(1-1/cop[k])/4.2/100/2/0.278;
}


}

end_f_loop(inlet,t)
}

DEFINE_PROFILE(Tout2,t, position)
{
face_t inlet,outlet;
Thread *thread;
real temp1,cop;
int ID=14;
Domain *domain;
domain=Get_Domain(1);

begin_f_loop(inlet,t)
{
   thread=Lookup_Thread(domain,ID);

      begin_f_loop(outlet,thread)
   {
         temp1=F_T(outlet,thread);
   }
      end_f_loop(outlet,thread)
      cop[k]=5*(-0.9830*(0.94395*pow((temp1-303)/23,2)+0.8683*(temp1-303)/23+1)+1.9709);
  fp=fopen("shuchu2.txt","a"

;

fprintf(fp,"%g ",cop[k]);

if(Q[i/2]>0) /*如果满足条件*/
{
   F_PROFILE(inlet,t,position)=temp1+Q[i/2]*(1+1/cop[k])/4.2/100/2/0.278;
}
else/*否则*/
{
   F_PROFILE(inlet,t,position)=temp1+Q[i/2]*(1-1/cop[k])/4.2/100/2/0.278;
}

}

end_f_loop(inlet,t);
k=k+1;
i=i+1;

}

» 猜你喜欢

依托企业入选了国家启明计划青年人才。有无高校可以引进的。已经有7人回复
遇见不省心的家人很难过已经有24人回复
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）已经有6人回复
AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性已经有6人回复
有院领导为了换新车，用横向课题经费买了俩车已经有9人回复
酰胺脱乙酰基已经有13人回复
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗已经有4人回复
有时候真觉得大城市人没有县城人甚至个体户幸福已经有10人回复

1楼 2015-01-09 11:21:08

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

有关fluent UDF的问题已经有15人回复
Ansys Fluent 编译UDF及加载UDF出现错误已经有14人回复
fluent中udf如何编写内部压力变化已经有7人回复
fluent编写UDF求助已经有3人回复
急，在FLUENT中如何同时调用两个UDF？已经有11人回复
Fluent中编写UDF之后无法初始化。。。已经有5人回复
fluent中UDF导入出错本人新手，求大神指点已经有6人回复
FLUENT udf 怎么获取某一时刻参量的值已经有11人回复
fluent中UDF编译的问题！已经有11人回复
求助LINUX下调用Fluent软件及UDF问题已经有5人回复
FLUENT中编写UDF时指针变量是否要初始化已经有4人回复
关于fluentUDF的一个小问题！已经有5人回复
fluent udf请教已经有16人回复
fluent udf的应用已经有20人回复
全部身家求大神帮忙解决fluent压力脉冲怎么编写UDF的问题已经有8人回复
fluent运行UDF时的出现'xcopy' 不是内部或外部命令····· 已经有15人回复
fluent UDF代码调用预定义程序的问题（如fortran程序）已经有15人回复
【求助】FLUENT UDF 气液两相流汽蚀汽化潜热(问题详见内容) 已经有21人回复
udf编程通过，但是在FLUENT中运行出错，麻烦udf高手帮忙看看已经有4人回复
【求助】单机多核fluent并行udf 已经有13人回复
【求助】fluent计算中，udf编译通过，初始化出错已经有11人回复

点击这里搜索更多相关资源

科研从小木虫开始，人人为我，我为人人

相关版块跳转我要订阅楼主 SNOWANGELLEI 的主题更新