24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

xiaoguai3

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 41
帖子: 23
在线: 8.4小时
虫号: 1543104
注册: 2011-12-19
专业: 可再生与替代能源利用中的

[求助] 新手求教UDF问题~

能解释，但一添加到边界或初始化就报错。
Error:
FLUENT received fatal signal (ACCESS_VIOLATION)
1. Note exact events leading to error.
2. Save case/data under new name.
3. Exit program and restart to continue.
4. Report error to your distributor.
Error Object: ()

我的UDF如下，二维模型，简单的固体传热问题，想在模拟过程的同时累加整个固体区域的所有得热量，达到一定值后改变热流边界值。由于刚接触UDF不久，不知道问题出在何处。请指教！d= Get_Domain(1);这个语句有问题吗？
#include "udf.h"
real Q_tot=0;
DEFINE_PROFILE (unsteady_heatflux, thread, position)
{
Domain *d;
real cp=880;
real density=2180;
real tempn, templ, volume, Q;
Thread *t;
cell_t c;
d= Get_Domain(1); /*Get the domain using Fluent utility *//*Loop over all cell threads in the domain*/
thread_loop_c(t,d) /*Compute Q */ /*Loop over all cells */
begin_c_loop(c,t)
{
real time= RP_Get_Real("flow-time"

;
real b=(int)(time/3600)+1;
int i=(int)((b-1)/24);
real Heat[3]={278.5, 280.7, 278.5};
volume=C_VOLUME(c,t); /* get cell volume */
templ=C_T_M1(c,t); /*Get cell tempertuare of previous step*/
tempn=C_T(c,t); /*Get cell tempertuare*/
Q=cp*density*volume*(tempn-templ);
Q_tot +=Q;
if (Q_tot<= Heat )
F_PROFILE(c,t,position)=60;
else F_PROFILE(c,t,position)=0;
}
end_c_loop(c,t)
}

回复此楼

1楼2012-01-04 20:57:23

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

gyctju

金虫 (正式写手)

应助: 314 (大学生)
金币: 938.9
散金: 227
红花: 48
沙发: 2
帖子: 992
在线: 427.5小时
虫号: 1569491
注册: 2012-01-08
专业: 分离过程

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
xiegangmai(金币+2): 谢谢应助！ 2012-01-08 22:58:47

可能你想法错了。边界条件部分不应该跟累加传热量的循环放在一起，这样的实际效果跟你想的是不一样的。
还有边界条件设定时应该用begin-f-loop
建议修改成如下形式试试
#include "udf.h"
real Q_tot=0;
DEFINE_PROFILE (unsteady_heatflux, thread, position)
{
Domain *d;
real cp=880;
real density=2180;
real tempn, templ, volume, Q;
Thread *t;
cell_t c;
d= Get_Domain(1); /*Get the domain using Fluent utility *//*Loop over all cell threads in the domain*/
real time= RP_Get_Real("flow-time";
real b=(int)(time/3600)+1;
int i=(int)((b-1)/24);
real Heat[3]={278.5, 280.7, 278.5};
thread_loop_c(t,d) /*Compute Q */ /*Loop over all cells */
begin_c_loop(c,t)
{

volume=C_VOLUME(c,t); /* get cell volume */
templ=C_T_M1(c,t); /*Get cell tempertuare of previous step*/
tempn=C_T(c,t); /*Get cell tempertuare*/
Q=cp*density*volume*(tempn-templ);
Q_tot +=Q;
}
end_c_loop(c,t)
}
begine_f_loop(c,t,)
{
if (Q_tot<= Heat )
F_PROFILE(c,t,position)=60;
else F_PROFILE(c,t,position)=0;
}
end_f_loop