24小时热门版块排行榜    

查看: 336  |  回复: 0

SNOWANGELLEI

新虫 (初入文坛)

[求助] 急需大神指点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;

}
回复此楼

» 猜你喜欢

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

智能机器人

Robot (super robot)

我们都爱小木虫

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

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 SNOWANGELLEI 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 此成果不能导入原因:元数据必填信息不完整,可 进行补充。 +3 Kittylucky 2026-03-02 4/200 2026-03-02 09:48 by kingkocxr
[考研] 材料调剂 +5 爱擦汗的可乐冰 2026-02-28 6/300 2026-03-02 09:36 by 破风α
[考研] 材料工程专硕283求调剂 5+3 ,!? 2026-03-02 3/150 2026-03-02 09:04 by 材子momo
[考研] 哈工大计算机刘劼团队招生 +3 hit_aiot 2026-03-01 5/250 2026-03-02 07:48 by 得鹿梦鱼111
[考研] 求调剂 +5 yunziaaaaa 2026-03-01 6/300 2026-03-01 23:57 by ccp273206157
[考研] 材料化工调剂 +12 今夏不夏 2026-03-01 13/650 2026-03-01 23:32 by L135790
[考研] 292求调剂 +6 yhk_819 2026-02-28 6/300 2026-03-01 23:23 by 向上的胖东
[考研] 0856调剂 +5 刘梦微 2026-02-28 5/250 2026-03-01 22:30 by wang_dand
[考研] 26考研报考西工大材料308分求调剂 +3 weizhong123 2026-03-01 3/150 2026-03-01 21:42 by 公瑾逍遥
[考研] 299求调剂 +3 Y墨明棋妙Y 2026-02-28 5/250 2026-03-01 21:01 by tangxiaotian
[考研] 306分材料调剂 +4 chuanzhu川烛 2026-03-01 5/250 2026-03-01 19:48 by 无际的草原
[考研] 一志愿中南大学理学化学 +4 15779376950 2026-03-01 5/250 2026-03-01 19:00 by Fff-1
[考研] 272求调剂 +6 材紫有化 2026-02-28 6/300 2026-03-01 18:58 by 18137688336
[考研] 321求调剂一志愿东北林业大学材料与化工英二数二 +4 虫虫虫虫虫7 2026-03-01 7/350 2026-03-01 16:52 by caszguilin
[考研] 311求调剂 +6 亭亭亭01 2026-03-01 6/300 2026-03-01 15:41 by 324616
[考研] 304求调剂 +6 曼殊2266 2026-02-28 7/350 2026-03-01 15:14 by wjLi2017
[考研] 课题组接收材料类调剂研究生 +3 gaoxiaoniuma 2026-02-28 4/200 2026-03-01 14:30 by jjj三跨
[考研] 材料284求调剂,一志愿郑州大学英一数二专硕 +10 想上岸的土拨鼠 2026-02-28 10/500 2026-03-01 14:12 by yc258
[考研] 304求调剂 +3 52hz~~ 2026-02-28 5/250 2026-03-01 00:00 by 52hz~~
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
信息提示
请填处理意见