24小时热门版块排行榜

>论坛更新日志 (4238)
>文献求助 (437)
>导师招生 (381)
>虫友互识 (381)
>休闲灌水 (259)
>考博 (172)
>硕博家园 (127)
>招聘信息布告栏 (120)
>博后之家 (98)
>论文投稿 (93)
>公派出国 (89)
>基金申请 (73)
>考研 (68)
>找工作 (60)
>论文道贺祈福 (58)
>教师之家 (41)

返回列表

【悬赏金币】回答本帖问题，作者cure727将赠送您 5 个金币

cure727

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 50.2
散金: 5
帖子: 45
在线: 10.3小时
虫号: 27735406
注册: 2021-11-14
性别: MM
专业: 建筑环境与结构工程

[求助] DPM颗粒采样UDF 已有1人参与

大家好，我想请教一下DPM里的有关UDF的一个小问题：
编写DPM采样对应的UDF的时候，结尾会涉及到如何保证不在重复追踪颗粒，我看大多数人最后的代码都写的是如下这种：
#if REMOVE_PARCELS
p->stream_index=-1;
#endif
但我也查阅了相关的UDF的help文档（2020版），里面的两个案例均采用的是如下这种结束语句
#if REMOVE_PARTICLES
MARK_TP(tp, P_FL_REMOVED);
#endif
我想问一下这两个结束语句有何区别呢？希望各位大佬不吝赐教

回复此楼

» 猜你喜欢

每天挣扎在fluent海洋里的科研小白，抑郁ing

1楼 2022-04-25 16:30:01

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

cure727

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 50.2
散金: 5
帖子: 45
在线: 10.3小时
虫号: 27735406
注册: 2021-11-14
性别: MM
专业: 建筑环境与结构工程

完整的UDF如下所示：
#define REMOVE_PARTICLES FALSE

DEFINE_DPM_OUTPUT(Discrete_phase_sample,header,fp,p,t,plane)
{

real flow_time = solver_par.flow_time;
real r, x, y, z;

if(header)
par_fprintf_head(fp," #Time[s] R [m] x-velocity[m/s] y-velocity[m/s] z-velocity[m/s] Drop Diameter[m] Number of Drops Temperature [K] Initial Diam [m] Injection Time [s] \n"

;

if(NULLP(p))
return;

x = p->state.pos[0];
y = p->state.pos[1];
z = p->state.pos[2];
r = sqrt(SQR(x) + SQR(y));

par_fprintf(fp,"%e %f %f %f %f %f %f %e %e %f %e %f \n",
P_TIME(p), x, y, z, p->state.V[0],p->state.V[1],p->state.V[2],P_DIAM(p),p->number_in_parcel,P_T(p), P_INIT_DIAM(p), p->time_of_birth);

#if REMOVE_PARTICLES
MARK_TP(tp, P_FL_REMOVED);
#endif
}