CFD中文网

其实我觉得展示UDF 特别是数值型的 有一个matlab的预测验证图会比较直观

组分还是DPM都可以 也都有缺点 一般性的11%以下的 采用dpm好些 我是做除尘的 略知一二 但是焊烟工艺没接触过

GCC下你的注释可能还是需要使用/**/这种 然后你最新的控制台输出 好像你的UDF已经编译成功了 你可以试试看行不行 GCC下其实UDF的调试更为方便

GCC和windows的MS C complier的语言要求性是不同的 有时候GCC的要求更加严格一些 你的代码的第18行能否贴一下

GCC只要内核没问题 一般Linux下Fluent可以UDF直接编译 不需要额外的环境变量配置 否则你需要查看 对应的Fluent makefile文件
刚开始可以使用Ubuntu这样的桌面Linux系统先熟悉一下 之后再使用Journal file 因为一般性的UDF还是需要不少调试环节的 有一个图形显示可能会好点 个人不成熟的意见 仅当参考

你的方法是建立ISO面 评价标准是体积分数50%的网格 这未必是准确的判定
第一你需要有一个准确的判定标准 比如网格体积分数大致介于什么之间的是属于边界属于接触
之后就是根据这个判断 利用udf和解析几何的办法判断下就行
udf主要使用c_vof就够了

@红豆沙 普通DEM 但是算力方面是一个问题 我最近在PT上发表了相关方面 能够在宏观上稍微解决这个问题
trap方法是碰到了直接脱离计算域 而且dpm本身不考虑体积 ddpm如果是沉积分散得较散 不适合 所有dem体系都是瞬态的 因为要考虑碰撞

后处理CFDPOST直接导出面数据 或者使用udf直接导出到.csv 然后直接进入matlab随便你怎么处理了

@深蓝 你的链接发的并没有写关键代码 这只是部分的 这部分只是排错 计算距离代码并没列出 实际上这部分的实现如果你对UDF有一定的了解并不是很困难 关键是执行效率的问题 如果不是每步执行的 是一步执行的话 是比较简单的 你的方法就可以满足 如果是每步执行的 先排除哪个也就是哪个if写在头 是执行效率的关键 你需要把最大概率不满足的条件写在判断的最开始 其他的都是解析几何的问题 可以找node 也可以找F_centroid

因为也看到了你在小木虫上的提问
我这里还是说一下
你这个示例 可能不是fluent做的 他更像是DEM的resolved加浸入边界法做的 这是属于CFDDEM的长项
国内北大力学王博士搞这个的 你可以跟他联系一下
你如果在fluent里面做这种 因为个数比较少 实际上还是比较适合6DOF的
那么你用6DOF一些反应就可能需要完全自主编程了
你用纯的DEM 实际上DEM是一个拉格朗日系 他中间对流体的阻碍只是来自于曳力等插值力 他不会真正有一个这样明显的流场轨迹 所以你可以考虑6DOF加上UDF的办法来做这个案例
但是UDF不会很简单 因为6DOF就不当做DPM系处理了 那么很多自带模型就用不了了
如果你说的颗粒一动不动 那么固定网格也可以做 只是这种方法可能没有太多的学术价值了

利用DPM_SCALAR_UPDATA

SV_locate_point (real [], CX_Cell_Id *);
你可以用这个函数来获取这个点的所在网格的信息 要自己去看下头文件
如果这是你的要求的话
而如果不是的话 那么只是一个纯解析几何的问题

@veen 你的尺度很大 这种属于一个地域级别的研究 在这种方式下 传统采用简单滑移模型 来处理颗粒项 你可以了解下 对应到fluent 根据我的经验 采用RSM+Mixture可能好点

理论确定方法 可以采用 Kozeny-Carman 和 egrun 公式 但是实际操作一般通过试验测得 理论公式的误差很大

直接的点输出 更好的办法采用cfd-post比较直接
如果是UDF内部 用SV_location 还是SV_locate 你自己再查下 我现在不太记得 用这个函数来找到点所处位置网格 然后输出数据就行 但是这个SV函数是一个比较蠢的算法 他比较慢 在并行下的效率未必高 这点需要注意

使用adjust宏 如果是用步数 以N_ITER做判断 否则用CURRENT_TIME或者N_TIME做判断
被模为0即可 在一些道路的污染物扩散方面 这种要求比较常见

@百善孝为先0626 liukan-12@163.com 真的有非常疑难的问题 确实影响到了你的一些重大进程 你可以联系我

@东岳 感谢东岳老师 目前我还在完善过程中 这个模型目前还存在相当多的问题 这个模型主要是集中在工程应用上 他的大部分数据是利用单管试验来获取的 所以是一个靶向非常集中的宏观模型 对其他工艺的适用性可能非常一般 再次感谢李老师的提议 我想在更进一步完善这个模型之后 对一些研究多孔过滤的人才能真正起到一定的意义

@wwzhao 再次感谢