注册 登录
流体中文网 返回首页

gearboy的个人空间 http://cfluid.com/?15372 [收藏] [复制] [分享] [RSS]

日志

UDF求DPM撞击壁面的角度

已有 995 次阅读2013-4-17 14:10 |个人分类:UDF使用技巧

使用DEFINE_DPM_BC宏来实现。记得hook。

#include "udf.h"
DEFINE_DPM_BC(usr_bc_reflect,p,t,f,f_normal,dim)
{
        real cos_angle,angle;
        cos_angle=NV_DOT(f_normal,P_VEL(p))/MAX(NV_MAG(P_VEL(p)),DPM_SMALL);
        cos_angle=MIN(1.0,MAX(-1.0,cos_angle));
        angle=90.0-acos(cos_angle)/M_PI*180.0;
        Message0("particle id=%d,angle between particle velocity and wall face=%f degree\n",p->part_id,angle);
        Reflect_Particle(p,f_normal,dim,f,t,t,f);  //if "trap" condition, then use "Trap_Particle(p);"
        return PATH_ACTIVE; // if "trap" condition, then use "return PATH_END;"
}

发表评论 评论 (1 个评论)

回复 passflame 2013-6-12 23:22
楼主你好!麻烦帮忙看看这个问题好吗?dpm模型中如果离散相和连续性发生传热传质,如果离散相经过某个区域发生强化传热传质现象,那应该怎么修改该区域离散相的传热传质系数?

facelist doodle 涂鸦板

您需要登录后才可以评论 登录 | 注册

返回顶部