24小时热门版块排行榜    

查看: 3074  |  回复: 6

leixiaobo713

银虫 (正式写手)

[求助] fluent中UDF编译的问题!

本人现在做fluent模拟的问题,需要用UDF对热源进行编译。其中热源是随时间变化的,当时间t=0时,热源为一个关于温度的函数,当t不等于0时,热源为0(也就是材料仅有一个初始热源)。下面为所编程序,但是老是出错,请高手指教,感激不尽!
#include "udf.h"
DEFINE_SOURCE(source_terms,c,t,dS,eqn)
{
real A=0.01;
real B=0.02;
real C=0.001;
real D=0.001;
real T=C_T(c,t);
real t=RP_Get_Real("flow-time" );
real source;
if(t=0)
  source=A*T+B*pow(T,2)+C*pow(T,3)+D*pow(T,4);
else
  source=0;
dS[eqn]=0;
return source;
}
结果interpreted时出现下面提示:
t definition shadows previous definition
Error: C:\Documents and Settings\Administrator\\327\300\303\346\Source terms.c: line 14: invalid expression type for if: float.
不知如何解决啊,里面既有时间t, 也有real T=C_T(c,t);

[ Last edited by leixiaobo713 on 2013-7-14 at 23:57 ]
回复此楼

» 收录本帖的淘帖专辑推荐

Fluent UDF

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

天道酬勤!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pointh

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
leixiaobo713: 金币+10, ★★★★★最佳答案 2013-07-15 09:03:30
xiegangmai: 金币+2, 谢谢参与! 2013-07-24 20:58:09
有两个问题:
1.你想用时间变量t,但这个t在宏DEFINE_SOURCE里面已经定义了;
2.“if (t=0)”这种写法不太好,因为计算机里面的数和我们想的不太一样,我们认为是0,但计算机里面它可能是0.00000000001这种,因此你判断等于时一定要注意一下。像这种你最好写成“if (t<1e-10)”之类的。
2楼2013-07-15 09:01:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leixiaobo713

银虫 (正式写手)

引用回帖:
2楼: Originally posted by pointh at 2013-07-15 09:01:11
有两个问题:
1.你想用时间变量t,但这个t在宏DEFINE_SOURCE里面已经定义了;
2.“if (t=0)”这种写法不太好,因为计算机里面的数和我们想的不太一样,我们认为是0,但计算机里面它可能是0.00000000001这种,因此 ...

多谢你的回复与指点,对于时间在DEFINE中已经定义,我是注意到的,因为是初学者,我不能确定DEFINE中定义的t是否就是时间t,因为C_T(c,t)为返回网格的温度。我在其他的一些程序中发现,对于一些参数,如热源,流速等随时间的变化,对t的定义都是通过real t=RP_Get_Real("flow-time" );来进行定义的。
天道酬勤!
3楼2013-07-15 09:15:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leixiaobo713

银虫 (正式写手)

引用回帖:
2楼: Originally posted by pointh at 2013-07-15 09:01:11
有两个问题:
1.你想用时间变量t,但这个t在宏DEFINE_SOURCE里面已经定义了;
2.“if (t=0)”这种写法不太好,因为计算机里面的数和我们想的不太一样,我们认为是0,但计算机里面它可能是0.00000000001这种,因此 ...

我按照你的说法把程度改了下,如下:
#include "udf.h"
DEFINE_SOURCE(source_terms,c,t,dS,eqn)
{
real A=0.01;
real B=0.02;
real C=0.001;
real D=0.001;
real T=C_T(c,t);
real source;
if(t<1e-10)
  source=A*T+B*pow(T,2)+C*pow(T,3)+D*pow(T,4);
else
  source=0;
dS[eqn]=0;
return source;
}
但是interpreted后出现下面的提示:
line 10: invalid type for binary expression: pointer to structure < float.
天道酬勤!
4楼2013-07-15 09:32:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wudeshou

铁杆木虫 (正式写手)

★
xiegangmai: 金币+1, 谢谢参与! 2013-07-24 20:58:19
if(t=0)
   source=A*T+B*pow(T,2)+C*pow(T,3)+D*pow(T,4);
else
   source=0;

t=0是赋值,判断用==

» 本帖已获得的红花(最新10朵)

无止境--环保--近零排放--多相流--晶体生长--数值计算
5楼2013-07-15 12:09:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leixiaobo713

银虫 (正式写手)

送红花一朵
引用回帖:
5楼: Originally posted by wudeshou at 2013-07-15 12:09:40
if(t=0)
   source=A*T+B*pow(T,2)+C*pow(T,3)+D*pow(T,4);
else
   source=0;

t=0是赋值,判断用==

真是多谢指点!!!
天道酬勤!
6楼2013-07-15 13:48:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

二亚

铜虫 (初入文坛)

太厉害了,呵呵
信念的力量是伟大的
7楼2013-09-28 17:22:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 leixiaobo713 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +10 ncTVBM6AOJHF 2026-11-04 22/1100 2026-11-06 17:18 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 7/350 2026-11-06 16:43 by 5sJFDKp2kwuI
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 16:43 by 5sJFDKp2kwuI
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ELBZ6WwEaRj6 2026-11-04 14/700 2026-11-06 16:32 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 16:23 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 15/750 2026-11-06 13:49 by BkA9Bi5cMhkD
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +7 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 13:47 by BkA9Bi5cMhkD
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ncTVBM6AOJHF 2026-11-04 17/850 2026-11-06 13:47 by BkA9Bi5cMhkD
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +8 ncTVBM6AOJHF 2026-11-04 18/900 2026-11-06 13:36 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 13:36 by BkA9Bi5cMhkD
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 05:32 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 12/600 2026-11-06 05:32 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 11/550 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +7 ncTVBM6AOJHF 2026-11-04 22/1100 2026-11-06 04:47 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 9/450 2026-11-06 02:15 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
[基金申请] 各位基础老师们,看完中央《从严治党》的决定后,对学术界更加有信心了吧 +5 iaeyuan 2026-11-01 5/250 2026-11-04 14:43 by jinchengzi
信息提示
请填处理意见