24小时热门版块排行榜    

查看: 3605  |  回复: 7
【悬赏金币】回答本帖问题,作者雨夜星猫将赠送您 50 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

雨夜星猫

新虫 (初入文坛)

[求助] fluent如何输出颗粒对壁面的碰撞速度和角度。已有2人参与

我用fluent模拟气固管流冲蚀,得到了壁面的冲蚀速率:
1、问题:fluent导出壁面数据后,缺乏颗粒对壁面的平均碰撞速度(轴向、径向、切向)和平均角度(入射角、反射角)这五个变量。 一个壁面位置受多次颗粒撞击,速度和角度有多个,需要统计一下平均。
2、UDF刚接触,按照模板修改了一点点,不知道怎么继续写下去,求高人指路!

#include "udf.h"
DEFINE_DPM_BC(bc_reflect,p,t,f,f_normal,dim)
{
real alpha; /* angle of particle path with face normal */
real vn=0.;
real nor_coeff = 0.;
real tan_coeff = 0;
real normal[3];
int i, idim = dim;
real NV_VEC(x);

for (i=0; i<idim; i++)
normal = f_normal;

if(p->type==DPM_TYPE_INERT)
{
alpha = M_PI/2. - acos(MAX(-1.,MIN(1.,NV_DOT(normal,P_VEL§)/
MAX(NV_MAG(P_VEL§),DPM_SMALL))));
if ((NNULLP(t)) && (THREAD_TYPE(t) == THREAD_F_WALL))
F_CENTROID(x,f,t);

   /* 定义壁面反射系数 */
  nor_coeff=0.993-0.0307*alpha+4.75e-4*alpha*alpha-2.61e-6*alpha*alpha*alpha
  tan_coeff=0.998-0.029*alpha+6.43e-4*alpha*alpha-3.56e-6*alpha*alpha*alpha

  /* calculate the normal component, rescale its magnitude by
     the coefficient of restitution and subtract the change */

  /* Compute normal velocity. */
  for(i=0; i<idim; i++)
    vn += P_VEL(p)*normal;

  /* Subtract off normal velocity. */
  for(i=0; i<idim; i++)
    P_VEL(p) -= vn*normal;

  /* Apply tangential coefficient of restitution. */
  for(i=0; i<idim; i++)
    P_VEL(p) *= tan_coeff;

  /* Add reflected normal velocity. */
  for(i=0; i<idim; i++)
    P_VEL(p) -= nor_coeff*vn*normal;  

  /* Store new velocity in P_VEL0 of particle */
  for(i=0; i<idim; i++)
    P_VEL0(p) = P_VEL(p);
  
  return PATH_ACTIVE;
}
return PATH_ABORT;
}

fluent如何输出颗粒对壁面的碰撞速度和角度。
1.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雨夜星猫

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by 王浩楠 at 2018-05-20 21:44:49
用UDM编写,参考UDF manual中的UDM
这是输出平均粒径的UDM;
F_UDMI(f,t,AVG_DIAMETER) = (P_DIAM(p) + num_in_data * F_UDMI(f,t,AVG_DIAMETER)) / (num_in_data + 1);
C_UDMI(c0,t0,AVG_DIAMETER) = F_UDMI(f, ...

UDMI需要初始化吗?
写在DEFINE_DPM_BC里面还是DEFINE_DPM_EROSION里面?
7楼2018-05-21 19:41:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

qiuchenwang

新虫 (小有名气)

楼主知道答案后 能不能拜托贴yixia

发自小木虫IOS客户端
2楼2018-05-19 09:46:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

moplye

禁虫 (初入文坛)

感谢参与,应助指数 +1
本帖内容被屏蔽

3楼2018-05-19 14:20:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qiuchenwang

新虫 (小有名气)


独孤神宇: 金币+1 2018-05-19 18:35:21
引用回帖:
3楼: Originally posted by moplye at 2018-05-19 14:20:06
通过udf编写获得...

谢谢啊 是因为要求平均 所以写udf 还是因为没有可能通过particle history导出来?

发自小木虫IOS客户端
4楼2018-05-19 18:00:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 博后面上今天有bug可以看到是否资助? +15 lyfbangong 2024-06-12 22/1100 2024-06-13 13:45 by Breezy恋
[有机交流] 原料反应完了,怎么知道是产物还是中间体 +6 小胡在努力 2024-06-11 8/400 2024-06-13 13:33 by 091602
[硕博家园] 硕博巨婴,也许才刚刚开始 +19 SNaiL1995 2024-06-12 45/2250 2024-06-13 13:02 by SNaiL1995
[有机交流] 车间生产,真空度很高,温度很高,但减压蒸馏速度很慢。 40+8 召唤鬼泣lL 2024-06-13 16/800 2024-06-13 12:43 by lzhy001
[基金申请] 为什么我的博后基金还在流动站审核中?不会是学院给我卡了吧? +6 王凯12 2024-06-13 11/550 2024-06-13 12:23 by MarkDeng
[找工作] 江西双非一本和四川双一流高校如何选择? 5+3 寒山敲钟 2024-06-12 11/550 2024-06-13 11:52 by wangtaotall
[找工作] 今年的就业形势真的很差吗 +13 215747018 2024-06-07 14/700 2024-06-13 11:26 by wenwen0825
[论文投稿] ACS Nano投稿后分配到副编辑手里12天了,能确定送审了吗? +5 潇洒怡惜 2024-06-12 9/450 2024-06-13 11:20 by 投必得科研顾问
[基金申请] 75批博后基金 +4 kyukitu 2024-06-13 5/250 2024-06-13 11:00 by muyunatti
[基金申请] 连续两年医学口青年项目初审体会 +9 进击的荣耀 2024-06-09 15/750 2024-06-13 08:54 by 进击的荣耀
[论文投稿] 摩擦磨损论文投稿 +3 jmysan 2024-06-12 3/150 2024-06-13 08:36 by 莱茵润色
[食品] 食品博士导师 +4 小李醒yy 2024-06-11 5/250 2024-06-12 19:15 by 投必得科研顾问
[硕博家园] 迷茫的毕业季 +7 独孤老狗 2024-06-08 7/350 2024-06-12 19:13 by 投必得科研顾问
[基金申请] 看到大家一直在说代表作NC,不明觉厉 +27 翰海2022 2024-06-07 29/1450 2024-06-12 15:30 by 3115321
[教师之家] 难上的课安排给我,杂事安排给我。评优时我就没份儿了,即便我做得不比绝大多数人差 +15 河西夜郎 2024-06-08 15/750 2024-06-12 12:51 by oooooo?o
[硕博家园] 求助 +6 LYWwrz 2024-06-09 9/450 2024-06-11 13:12 by powerhours
[公派出国] 公派出去要准备多少钱? +7 天空之光 2024-06-06 9/450 2024-06-11 08:49 by tfang
[论文投稿] 审稿交流 +5 huyugan123 2024-06-08 7/350 2024-06-10 10:51 by huyugan123
[基金申请] 什么时候会评呢? +7 qq632458 2024-06-07 9/450 2024-06-10 10:07 by Howard28
[基金申请] B01化学口的,近几年要什么样的评价才能上会呀? +4 salmon95 2024-06-06 5/250 2024-06-09 09:45 by chem-cai
信息提示
请填处理意见