24小时热门版块排行榜    

查看: 3587  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 数理学部函评几号结束? +4 科研孤勇者 2024-05-16 4/200 2024-05-17 00:37 by D调的XF
[教师之家] 南京林业大学老师首聘期考核不过,服药自杀 +30 babu2015 2024-05-12 33/1650 2024-05-16 23:14 by Siuhong Yiu
[论文投稿] 投稿成功后又想撤回 +3 otani 2024-05-16 3/150 2024-05-16 23:06 by lizhengke06
[教师之家] 博士去高校就是为了有寒暑假吗? +9 wenwen0825 2024-05-16 14/700 2024-05-16 23:05 by sakuraai
[硕博家园] 耐高温垫片求购 +3 Sexyflea 2024-05-16 4/200 2024-05-16 22:20 by cnofc
[考博] 寻找高分子领域博导 +8 小政雅图 2024-05-14 11/550 2024-05-16 22:10 by 单调的胖子
[硕博家园] 领养小孩 +14 Yyds55 2024-05-14 14/700 2024-05-16 22:06 by 鱼翔浅底1
[基金申请] 12个本子最多让给1a2b +16 地球e村长 2024-05-13 19/950 2024-05-16 17:48 by sars518
[基金申请] 粗大,科研之友提示有一篇10多年前的文章被阅读 +9 babu2015 2024-05-16 11/550 2024-05-16 14:53 by TsolmonZ
[分析] 液相色谱基线上移 5+3 江上鸢飞 2024-05-13 3/150 2024-05-16 11:49 by gwmgyp
[基金申请] 南林海归博士首聘期 没拿到国自然,被降级后自杀 +31 babu2015 2024-05-13 33/1650 2024-05-16 09:57 by Siuhong Yiu
[教师之家] 宋凯 +10 yexuqing 2024-05-14 10/500 2024-05-16 09:21 by shl2112501
[硕博家园] 北京航空航天大学计算机学院罗洪斌课题组招收2024年学术型博士研究生 +3 yanfeienter 2024-05-12 7/350 2024-05-15 22:52 by yanfeienter
[论文投稿] 当前scientific reports还值得投稿么? +4 lizhengke06 2024-05-14 5/250 2024-05-15 20:16 by vmuch520
[硕博家园] 博导推荐,申博,博士招生 +5 大林挺好 2024-05-13 6/300 2024-05-15 14:51 by zxl_1105
[教师之家] 相比过去,现在高校导师水平越来越高(可见招聘条件),研究生越来越差 +14 苏东坡二世 2024-05-11 18/900 2024-05-15 11:04 by 清晨和风一起
[考研] 求考研导师(24年底我才考,材料类,提前联系) +5 笑笑宝公主 2024-05-12 6/300 2024-05-15 10:53 by 13096145
[论文投稿] 投稿竟然倒在了Scientific Reports上。 +12 lizhengke06 2024-05-10 14/700 2024-05-15 09:34 by zqwcr
[基金申请] F06送审了吗? +5 昂首走 2024-05-11 9/450 2024-05-15 06:53 by yp_nupt
[论文投稿] Journal of Electrical Engineering&amp;Technology Reviews Completed 快一周了 10+4 qweasd12345 2024-05-12 5/250 2024-05-13 10:00 by bear2007
信息提示
请填处理意见