| 查看: 3551 | 回复: 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; }  1.png |
回复此楼» 猜你喜欢
氧化硅薄膜玻璃 求购
已经有4人回复
-
vasp极化强度
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有238人回复
-
vasp计算的一些问题
已经有1人回复
-
国科温州研究院 项目聘用有意“二战”的优秀本科毕业生
已经有31人回复
-
SCAPS-1D学习交流
已经有2人回复
-
深圳大学跟哈工大联合博士
已经有11人回复
-
Journal of alloys and compounds期刊
已经有2人回复
-
请教一个概率问题
已经有0人回复
感谢参与，应助指数 +1
|
本帖内容被屏蔽 |
3楼2018-05-19 14:20:06
5