24小时热门版块排行榜     石溪大学接受考研调剂申请>

【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 3544  |  回复: 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的回帖

moplye

禁虫 (初入文坛)

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

3楼2018-05-19 14:20:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

qiuchenwang

新虫 (小有名气)

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

发自小木虫IOS客户端
2楼2018-05-19 09:46:08
已阅   回复此楼   关注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的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 没学上 +7 季向阳 2024-04-26 13/650 2024-04-27 20:38 by 一条咸鱼.
[高分子] 请问UV灯是365nm的,那么选光引发剂的波长选多少的?要完全一致吗? +5 engledd2004 2024-04-21 5/250 2024-04-27 18:32 by quirinus
[考博] 博士招生-211农业院校 +5 NYC917 2024-04-26 7/350 2024-04-27 17:43 by 小梅菌
[考博] 材料方向24博士申请/一作SCI三篇 +3 白天不碰 2024-04-24 5/250 2024-04-27 17:41 by 安塔瓦拉多
[找工作] 普通院校药学硕士,做合成的,感觉找不到工作 +11 pom戴墨镜 2024-04-24 19/950 2024-04-27 15:15 by 奔赴远方
[基金申请] 两类问题算是白选了~ +8 jurkat.1640 2024-04-23 13/650 2024-04-27 12:03 by 淀粉搬运工
[基金申请] &quot;颜宁:基础研究应顶天立地&quot;能做到基础研究同时顶天立地的才是牛人 +6 zju2000 2024-04-24 6/300 2024-04-27 10:17 by kcmn1000
[论文投稿] 投稿RSC旗下杂志突然间看不到投稿状态,里面投稿的文章信息不见了,有遇到过的吗? 50+3 dlying 2024-04-26 3/150 2024-04-27 09:58 by bobvan
[论文投稿] 求论文投稿期刊推荐 15+3 mihudawang 2024-04-26 4/200 2024-04-27 09:56 by bobvan
[硕博家园] 聊天 +12 暮色恋伊人 2024-04-22 13/650 2024-04-27 08:40 by WASM
[考研] 学硕专硕 +5 小蜗牛* 2024-04-26 5/250 2024-04-26 16:43 by 鱼翔浅底1
[考研] 0854-0855调剂 +8 shangannum1 2024-04-21 12/600 2024-04-26 16:42 by yz仔
[硕博家园] 考研,求职还是考编? +15 xizj 2024-04-21 24/1200 2024-04-26 11:49 by Kan客
[教师之家] 某种做法不行。说过几遍了。同学还那样做。再那样做就给低分 +4 河西夜郎 2024-04-24 4/200 2024-04-26 08:51 by Quakerbird
[考博] 求博导 +6 好okjh 2024-04-21 10/500 2024-04-25 14:04 by 好okjh
[基金申请] 国社科项目,你们学校都限额申报吗? +7 屡战屡败 2024-04-21 10/500 2024-04-25 12:10 by 屡战屡败
[考博] 24年 申博 化学/材料 一作6篇sci +9 wangyp123 2024-04-23 11/550 2024-04-24 19:01 by bangbangbiu
[考博] 博士招生 +4 zx179 2024-04-24 7/350 2024-04-24 15:01 by H考研成功
[论文投稿] 期刊推荐 20+4 木颜尘ip 2024-04-22 7/350 2024-04-24 10:06 by bobvan
[考博] 研二光催化6月底4篇2区 +7 wjtab 2024-04-22 11/550 2024-04-23 06:59 by byron2012
信息提示
请填处理意见