24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者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
我想问一下这两个结束语句有何区别呢？希望各位大佬不吝赐教

回复此楼

» 猜你喜欢

2026年国自然面上资助率已经有21人回复
基础研究怎么拉横向，学校到款任务越来越多，难以完成拉横向，都有哪些途径啊已经有10人回复
2027年申博已经有3人回复
一个有机合成实验室都需要哪些设备？已经有7人回复
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！已经有22人回复
微信指数没变化，科研之友没阅读已经有18人回复
HXDI做水性聚氨酯乳液，是不是特别容易出渣已经有3人回复
系统今天又提示维护了，估计离放榜不远了已经有15人回复
你们的时间戳变了吗已经有4人回复
产物和副产物价值比较已经有5人回复

每天挣扎在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
}