24小时热门版块排行榜    

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

hkj进

新虫 (初入文坛)

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

发自小木虫Android客户端
8楼2019-09-02 09:52:47
已阅   回复此楼   关注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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 博后特别资助状态变化 +24 随梦而飞2017 2024-05-30 35/1750 2024-06-01 10:10 by 青岛阳仔
[硕博家园] 每天学术时间不能保证,能保证的只有: +8 hahamyid 2024-05-27 8/400 2024-06-01 00:33 by 鱼翔浅底1
[考博] 广东以理材料系碳点与功能材料课题组 — 2博士名额 / 科研助理 +4 小城夜很美 2024-05-27 11/550 2024-05-31 21:26 by 小城夜很美
[考博] 申请2024或2025年博士研究生 +5 嘟噜嘟1 2024-05-29 11/550 2024-05-31 19:27 by 嘟噜嘟1
[基金申请] 工材01送了吗? +11 xiaopang8958 2024-05-25 19/950 2024-05-31 18:31 by WORLD0256
[基金申请] 离职后国自然项目不能变更单位,在新单位还能申请新的国自然项目吗 5+3 宜兰 2024-05-31 6/300 2024-05-31 12:05 by felicity6056
[考博] 24or25材料专业申博 +3 农夫三拳有点痛 2024-05-30 10/500 2024-05-31 11:17 by 安塔瓦拉多
[基金申请] 入职高校3年发表10+SCI,尽人事听天命 +29 kaoyan250 2024-05-27 40/2000 2024-05-31 08:44 by Xiaolin81
[文学芳草园] 物是人非 +4 myrtle 2024-05-30 4/200 2024-05-30 15:05 by mapenggao
[电化学] 有锂金属软包电池代工公司代做软包嘛 +5 yyywj995 2024-05-26 6/300 2024-05-30 08:33 by 那年的旧城
[硕博家园] 又想换工作 +15 brightmj 2024-05-27 20/1000 2024-05-29 23:25 by zyqchem
[博后之家] 2024公派博后申请 +4 326lhpqk 2024-05-27 5/250 2024-05-29 20:03 by @古月胡
[基金申请] 信息学部函评结束了吗? +6 ducan21 2024-05-28 7/350 2024-05-29 12:10 by WORLD0256
[论文投稿] 真急着毕业,CPB主编终审17天了,邮件催稿了两次,就是一点动静没有 5+3 kkkk夏 2024-05-28 6/300 2024-05-29 11:18 by hitsdu
[基金申请] E10开始送了,希望有好运 +5 sail 2024-05-27 5/250 2024-05-28 18:36 by 芝小芝
[论文投稿] EI期刊审稿人邮箱问题 5+3 shier妈妈 2024-05-27 4/200 2024-05-28 14:53 by topedit
[有机交流] 机理求助 120+4 15147165026 2024-05-26 10/500 2024-05-28 14:42 by 江东闲人
[基金申请] 面上基金会评专家,有回避机制吗? +4 huang1991js 2024-05-27 4/200 2024-05-27 19:08 by 星火12
[硕博家园] 研0 +6 控制调剂yl 2024-05-25 11/550 2024-05-26 23:37 by sakuraai
[硕博家园] 2024博士招生 +3 big 混子 2024-05-26 3/150 2024-05-26 20:47 by 宁多缺毋滥
信息提示
请填处理意见