24小时热门版块排行榜    

查看: 3511  |  回复: 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的回帖

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的回帖

王浩楠

新虫 (初入文坛)

【答案】应助回帖

用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,t,AVG_DIAMETER);
5楼2018-05-20 21:44:49
已阅   回复此楼   关注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, ...

有相似研究者可加微信交流18795866867
6楼2018-05-20 21:46:09
已阅   回复此楼   关注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的回帖

hkj进

新虫 (初入文坛)

楼主能交流下吗?我也在研究碰撞频率这方面q1076449646

发自小木虫Android客户端
8楼2019-09-02 09:52:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 雨夜星猫 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 310求调剂 +4 cxr2024 2024-03-28 5/250 2024-03-28 22:52 by ms629
[考研] 一志愿东北大学材料专硕361分求调剂 +3 许你万丈光芒 2024-03-28 4/200 2024-03-28 21:38 by EBSD
[考研] 323求调剂 +3 橙子味/汽水 2024-03-28 3/150 2024-03-28 20:56 by 20240324
[考研] 315求调剂 +3 Sorappp 2024-03-28 3/150 2024-03-28 20:53 by renwenfeng
[硕博家园] a老师一作学生二作 +8 zgshsjshgd 2024-03-25 16/800 2024-03-28 20:46 by zgshsjshgd
[考研] 287求调剂 +5 重获新的我 2024-03-27 5/250 2024-03-28 20:12 by 麦田里的守望者
[考研] 成绩297化学工程,想调剂,13935055270 +4 努力加油油 2024-03-28 4/200 2024-03-28 20:05 by 麦田里的守望者
[考研] 本科太原理工大学,300分求调剂 +5 藏春山 2024-03-27 5/250 2024-03-28 16:12 by shiyihelulu
[考研] 305求调剂 +3 菠萝菠萝汁 2024-03-28 3/150 2024-03-28 10:11 by 随风而动咦
[有机交流] 溴素溴化反应生成溴苄 100+4 tanjishuang 2024-03-25 6/300 2024-03-28 09:37 by 高树桩
[有机交流] 合成实验紧急求助,各位大神帮忙! 50+4 1595515531 2024-03-26 5/250 2024-03-28 09:33 by 高树桩
[考研] 中北大学化学与化工学院招生 +3 山鬼的信箱 2024-03-25 5/250 2024-03-28 09:06 by 山鬼的信箱
[考研] 昆明理工大学找调剂 +12 法12外狂徒 2024-03-22 12/600 2024-03-28 08:35 by 菠萝菠萝汁
[考研] 化学工程与技术320求调剂 +10 Wjrjxmt 2024-03-26 11/550 2024-03-27 22:39 by weizhongzeng
[考研] 287求调剂 +9 王杨具牢 2024-03-24 11/550 2024-03-27 20:55 by siwuhen
[考研] 一志愿985,总分324,材料与化工专硕求调剂 +10 samuelqyp 2024-03-22 18/900 2024-03-27 20:44 by siwuhen
[考研] 330求调剂 +5 NJFU-LXD 2024-03-26 10/500 2024-03-27 14:28 by 小美啦啦啦啦啦
[考研] 生物学381求调剂,一志愿中山大学 +4 brix贝 2024-03-26 6/300 2024-03-27 12:47 by brix贝
[考研] 304材料专硕求调剂 +4 烟云墨雨飞啊 2024-03-23 4/200 2024-03-26 13:03 by hello呀6
[论文投稿] IEEE Transactions on Industrial Informatics的审稿周期是多久? +3 mollyzhang_2003 2024-03-22 5/250 2024-03-23 07:02 by mollyzhang_2003
信息提示
请填处理意见