24小时热门版块排行榜    

Znn3bq.jpeg
北京石油化工学院2026年研究生招生接收调剂公告
查看: 984  |  回复: 3

xujinchao

铜虫 (小有名气)

[求助] 请问fluent两相流拽力模型,加载到三相流中,就出错,请问需要怎么设置或修改 已有1人参与

谢谢大家

请问fluent两相流拽力模型,加载到三相流中,就出错,请问需要怎么设置或修改
11.png
回复此楼

» 收录本帖的淘帖专辑推荐

模拟相关

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xujinchao

铜虫 (小有名气)

/***************************************************************
   UDF for customizing the Gidaspow drag law in Fluent
****************************************************************/

#include "udf.h"
#include "sg_mphase.h"

#define diam2 1.7e-5

DEFINE_EXCHANGE_PROPERTY(custom_drag_syam, cell, mix_thread, s_col, f_col)
{
Thread *thread_g, *thread_s;
real x_vel_g, x_vel_s, y_vel_g, y_vel_s, abs_v, slip_x, slip_y,
      rho_g, rho_s, mu_g, reyp, afac,
      bfac, void_g, vfac, fdrgs, taup, k_g_s;

thread_g = THREAD_SUB_THREAD(mix_thread, s_col);
thread_s = THREAD_SUB_THREAD(mix_thread, f_col);

x_vel_g = C_U(cell, thread_g);
y_vel_g = C_V(cell, thread_g);
x_vel_s = C_U(cell, thread_s);
y_vel_s = C_V(cell, thread_s);
slip_x = x_vel_g - x_vel_s;
slip_y = y_vel_g - y_vel_s;
rho_g = C_R(cell, thread_g);
rho_s = C_R(cell, thread_s);

mu_g = C_MU_L(cell, thread_g);
/*compute slip*/
abs_v = sqrt(slip_x*slip_x + slip_y*slip_y);

/*compute reynolds number*/
reyp = rho_g*abs_v*diam2/mu_g;

/* compute particle relaxation time */
taup = rho_s*diam2*diam2/18./mu_g;
void_g = C_VOF(cell, thread_g);/* gas vol frac*/

/*compute drag and return drag coeff, k_g_s*/
afac = pow(void_g,4.14);

if(void_g<=0.85)
  bfac = 0.8*pow(void_g, 1.28);
else
  bfac = pow(void_g, 2.65);

vfac = 0.5*(afac-0.06*reyp+sqrt(0.0036*reyp*reyp+0.12*reyp*(2.*bfac-
              afac)+afac*afac));
fdrgs = void_g*(pow((0.63*sqrt(reyp)/vfac+4.8*sqrt(vfac)/vfac),2.))/24.0;

k_g_s = (1.-void_g)*rho_s*fdrgs/taup;

return k_g_s;
}
2楼2015-08-12 23:22:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xujinchao

铜虫 (小有名气)

没有人啊,只能自己顶一下了
3楼2015-08-13 14:48:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aaaalu

铜虫 (初入文坛)

【答案】应助回帖

你的意思是说用你编的这个UDF来计算三相流中气固之间的作用力吧?
你先检查一下物性。
然后调一下初始条件,比如设气含率0.00001
至少先让这个Case跑起来。
4楼2015-12-07 22:09:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xujinchao 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 22408 调剂材料 +7 我叫ez 2026-04-06 8/400 2026-04-07 17:12 by 蓝云思雨
[考研] 293求调剂 +3 勇远库爱314 2026-04-06 3/150 2026-04-07 11:15 by hugr
[考研] 考研调剂 +7 15615482637 2026-04-04 7/350 2026-04-06 22:56 by chenzhimin
[考研] 308求调剂 +11 墨墨漠 2026-04-06 11/550 2026-04-06 22:56 by 875465
[考研] 287分求调剂 有专利国奖一志愿哈工大085406 +6 白易辰 2026-04-06 7/350 2026-04-06 22:46 by 875465
[考研] 一志愿生物与医药,296分,求调剂 +11 66鹿 2026-04-03 13/650 2026-04-06 21:45 by 52305043001
[考研] 287求调剂 +3 通信学硕081000 2026-04-03 4/200 2026-04-06 21:03 by going home
[考研] 一志愿 江南大学 085602 化工专硕 338分求调剂 +15 路痴小琪 2026-04-05 15/750 2026-04-06 09:27 by cql1109
[考研] 372分,材料与化工,一志愿湖南大学,求调剂 +3 蓝笺片 2026-04-01 3/150 2026-04-06 09:04 by 无际的草原
[考研] 085600,320分求调剂 +7 大馋小子 2026-04-01 8/400 2026-04-05 21:19 by 学员8dgXkO
[考研] 工科08专硕机械275求调剂 +3 AaAa7420 2026-04-02 3/150 2026-04-05 13:26 by jp9609
[考研] 282电子信息0854专硕调剂 +4 202451007219 2026-04-02 6/300 2026-04-04 21:55 by laoshidan
[考研] 316求调剂 +9 墨辰_Orion926 2026-04-04 9/450 2026-04-04 21:35 by lbsjt
[考研] 一志愿华北电力大学(北京),材料科学与工程学硕265,求调剂 +11 yelck 2026-04-03 12/600 2026-04-04 19:52 by dongzh2009
[考研] 建环,能源,土木老师路过看一看!!! +5 嘿嘿uu 2026-04-01 5/250 2026-04-03 11:47 by znian
[考研] 085801 总分275 本科新能源 求调剂 +19 bradoner 2026-04-01 23/1150 2026-04-03 10:07 by linyelide
[考研] 专硕 351 086100 也是考的材科基 本科也是材料 +8 202451007219 2026-04-02 8/400 2026-04-03 09:50 by 蓝云思雨
[考研] 材料调剂 +4 一样YWY 2026-04-03 4/200 2026-04-03 09:48 by 蓝云思雨
[考研] 22408 266求调剂 +3 masss11222 2026-04-02 3/150 2026-04-02 18:11 by 笔落锦州
[考研] 材料求调剂 +10 呢呢妮妮 2026-04-01 13/650 2026-04-02 09:17 by olim
信息提示
请填处理意见