24小时热门版块排行榜    

查看: 1077  |  回复: 2
【悬赏金币】回答本帖问题,作者无敌小张!将赠送您 5 个金币

无敌小张!

新虫 (初入文坛)

[求助] fluent的udf

大神们,有段fluent的udf,但是看不懂,有没有大神能给解答一下是什么意思
研一萌新没啥金币,直接分享,大家有需要就拿去用
#include "udf.h"
#define U 18
#define I 200
#define PI 3.1415928
#define R 0.003
#define n 0.8
#define v 0.003
#define L 0.005
#define u0 1.26e-6

DEFINE_SOURCE(S_X,cell,thread,dS,eqn)
{
         real xc[ND_ND];   
         real t;
         real source;
         real x,y,z;
     real r;

         C_CENTROID(xc,cell,thread);
         x=xc[0];   
         y=xc[1];
         z=xc[2];

         t=RP_Get_Real("flow-time";
         r=sqrt((x-v*t)*(x-v*t)+y*y);
         if(r<R)
         {
                 source=-1*u0*I*I/(4*PI*PI*R*R*r)*exp(-1*r*r/(2*R*R))
                         *(1-exp(-1*r*r/(2*R*R)))*(1-z/L)*(1-z/L)*x/L;
                 dS[eqn]=0;
         
         }
         else
                 source=dS[eqn]=0;
             return source;

}
DEFINE_SOURCE(S_Y,cell,thread,dS,eqn)
{
         real xc[ND_ND];
         real t;
         real source;
         real x,y,z;
     real r;

         C_CENTROID(xc,cell,thread);  
         x=xc[0];
         y=xc[1];
         z=xc[2];

         t=RP_Get_Real("flow-time";
         r=sqrt((x-v*t)*(x-v*t)+y*y);
         if(r<R)
         {
                 source=-1*u0*I*I/(4*PI*PI*R*R*r)*exp(-1*r*r/(2*R*R))
                         *(1-exp(-1*r*r/(2*R*R)))*(1-z/L)*(1-z/L)*y/L;
                 dS[eqn]=0;
         
         }
         else
                 source=dS[eqn]=0;
             return source;

}
  DEFINE_SOURCE(S_Z,cell,thread,dS,eqn)
{
         real xc[ND_ND];
         real t;
         real source;
         real x,y,z;
     real r;

         C_CENTROID(xc,cell,thread);
         x=xc[0];
         y=xc[1];
         z=xc[2];

         t=RP_Get_Real("flow-time";
         r=sqrt((x-v*t)*(x-v*t)+y*y);
         if(r<R)
         {
                 source=-1*u0*I*I/(4*PI*PI*L*r*r)
                         *(1-exp(-1*r*r/(2*R*R)))*(1-z/L)*(1-z/L);
                 dS[eqn]=0;
         
         }
         else
                 source=dS[eqn]=0;
             return source;


  }
DEFINE_PROFILE(heat_flux,thread,index)
{

         real xc[ND_ND];
         face_t f;
         real t;
         real heat,T0=300;
         real x,y,z;
     real r;
         real hc = 80, sigmab = 5.67e-8, epsilon = 0.06;
     t=RP_Get_Real("flow-time";
         begin_f_loop(f,thread)
         {
           F_CENTROID(xc,f,thread);
             x=xc[0];
             y=xc[1];
             z=xc[2];
         r=sqrt((x-v*t)*(x-v*t)+y*y);
                 if(r<sqrt(6)*R)
           heat=n*I*U/(2*PI*R*R)*exp(-1*r*r/(2*R*R));
                 else
                   heat=0;
                 
                 F_PROFILE(f, thread, index) = heat - hc*(F_T(f,thread), - T0) - sigmab*epsilon*(pow(F_T(f,thread), 4.) - pow(T0, 4.));

         }
     end_f_loop(f,thread)
     
  }
回复此楼

» 猜你喜欢

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

heartbute

新虫 (初入文坛)

2楼2021-09-12 00:03:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

朴勋医生

新虫 (小有名气)

最后一个面热源,前面都是体热源,你也做激光加工吗

发自小木虫Android客户端
3楼2021-10-25 09:00:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 无敌小张! 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 自然基金交了,近段时间尽量别看标书 +12 hdxnhjtx 2024-03-27 12/600 2024-03-28 22:07 by dxcharlary
[考研] 323求调剂 +3 橙子味/汽水 2024-03-28 3/150 2024-03-28 20:56 by 20240324
[考研] 315求调剂 +3 Sorappp 2024-03-28 3/150 2024-03-28 20:53 by renwenfeng
[考研] 双一流高校南林,化学/化工/材料 调剂联系 +10 wang333666 2024-03-24 15/750 2024-03-28 20:50 by kanghw0101
[考研] 总分323,一志愿华东理工 +10 努_力66 2024-03-23 10/500 2024-03-28 20:10 by 麦田里的守望者
[考研] 318化工求调剂 +5 蜜糖不甜Z 2024-03-23 5/250 2024-03-28 15:56 by qfGexu
[考研] 985收调剂 +23 A_Laity! 2024-03-27 32/1600 2024-03-28 13:38 by kawhichen
[考研] 320化学工程与技术求调剂 +8 Wjrjxmt 2024-03-27 8/400 2024-03-28 10:46 by 随风而动咦
[论文投稿] 找人指导sci +4 知欧猪- 2024-03-27 4/200 2024-03-28 10:38 by SPE翻译
[考研] 一志愿985,427求调剂,求求 +11 echo586 2024-03-25 21/1050 2024-03-28 08:41 by echo586
[考研] 279求调剂一志愿华东理工材料与化工085602专硕 +4 yt0312 2024-03-24 7/350 2024-03-27 16:21 by barlinike
[考研] 315求调剂 +4 Sorappp 2024-03-27 5/250 2024-03-27 13:48 by 瀚海星云
[有机交流] 硅胶板保存问题及活化 10+3 2651230415 2024-03-23 4/200 2024-03-27 10:27 by mrzhl1986
[考研] 化学工程与技术320求调剂 +12 Wjrjxmt 2024-03-23 14/700 2024-03-26 13:01 by hello呀6
[考研] 308求调剂 +5 catherinetyx 2024-03-24 5/250 2024-03-26 13:01 by hello呀6
[论文投稿] ACS Catalysis的投稿状态 60+4 myrice112 2024-03-22 11/550 2024-03-26 08:20 by qq6986253
[考研] 317求调剂 +4 (?????)????? 2024-03-25 5/250 2024-03-25 23:58 by DJ35
[考研] 电子信息380 求调剂 +3 Arthur2 2024-03-25 8/400 2024-03-25 20:46 by 计算机专硕345调
[考研] 一志愿北京工业大学,材料工程,总分301求调剂 +9 hss0418 2024-03-24 14/700 2024-03-25 20:28 by EBSD
[考研] 327求调剂 +4 忘却七秒情 2024-03-22 4/200 2024-03-24 14:01 by 小雨哥
信息提示
请填处理意见