24小时热门版块排行榜    

查看: 1734  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 标题:捡漏预警|08工科/09农学调剂!英语要求低,过线即有机会! +7 马超放烟花 2026-03-07 12/600 2026-03-10 23:16 by Equinoxhua
[考研] 288求调剂 +12 王晓阳- 2026-03-09 17/850 2026-03-10 22:16 by peike
[考研] 材料工程307,求调剂 +6 我要燃烧你的梦 2026-03-08 6/300 2026-03-10 22:08 by peike
[考研] 材料与化工(0856)304求B区调剂 +4 邱gl 2026-03-10 7/350 2026-03-10 21:22 by 英语大词典
[考博] 26申博求助 +3 跳跃饼干 2026-03-10 4/200 2026-03-10 21:15 by Tntcnn
[考研] 一志愿山东大学105500药学专硕,总分302求调剂 +5 五维天空 2026-03-04 13/650 2026-03-10 19:02 by 清风月
[考研] 材料与化工,291,求调剂 +12 咕噜咕噜123123 2026-03-05 13/650 2026-03-10 13:50 by allen-yin
[考研] 化工284求调剂 +12 小黑想变白 2026-03-07 12/600 2026-03-10 09:17 by Wahxp
[考研] 求调剂,不管什么专业,我是可塑造的人才一枚,希望遇到知己老师捞捞我 +4 13102137290 2026-03-06 5/250 2026-03-07 21:21 by lissomchan
[考研] 070300化学求调剂292分 +3 打烊eee 2026-03-07 3/150 2026-03-07 20:04 by houyaoxu
[考研] 材料专硕323求调剂 +5 李白26 2026-03-07 5/250 2026-03-07 16:46 by njzyff
[考研] 一志愿哈尔滨工业大学0856材料与化工,前三科206,总分283,求调剂 +7 26考研求调剂 2026-03-06 7/350 2026-03-07 14:49 by Miko19
[考研] 一志愿中石油,化学工程,301求调剂 +4 豚馒枳枳 2026-03-04 4/200 2026-03-06 22:15 by 星空星月
[考研] 复试调剂 +7 呼呼?~+123456 2026-03-05 10/500 2026-03-06 16:07 by @飒飒飒飒
[考研] 08工科求调剂 +3 隆LLL 2026-03-06 4/200 2026-03-06 15:26 by 隆LLL
[考研] 274求调剂 +9 一个学习者 2026-03-04 9/450 2026-03-05 20:56 by 黑衣馒头人
[考研] 376求调剂 +3 王浩然sam 2026-03-04 3/150 2026-03-05 07:48 by bxbo
[考研] 293求调剂 +3 是乐渝哇 2026-03-04 3/150 2026-03-04 23:12 by wutongshun
[考研] 295求调剂 +6 等春来, 2026-03-04 6/300 2026-03-04 16:11 by zhukairuo
[考研] 322,求调剂 +3 菜菜爱玩 2026-03-04 3/150 2026-03-04 12:15 by xiongyaxuan
信息提示
请填处理意见