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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » fluent中UDF编译的问题！
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知（长期有效）
51/1返回列表
查看: 2196  |  回复: 6
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
天道酬勤！
1楼2013-07-14 23:55:14
已阅   回复此楼   关注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的回帖
查看全部 7 个回答

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)”之类的。
赞一下(3人) 回复此楼
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的回帖

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是赋值，判断用==
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

leixiaobo713
无止境--环保--近零排放--多相流--晶体生长--数值计算
5楼2013-07-15 12:09:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭