24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 3397  |  回复: 5

jyl432

新虫 (初入文坛)

[交流] fluent计算多相流问题时,UDF能否提取单元内分相速度,以及单元平均速度 已有4人参与

请问在vof两相流计算中,因为编写udf是要用到两相之间的速度差(Ug-Ul)。UDF中能否使用命令提取出cell中各分相的速度,或是两相间的滑移速度如何提取。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

znon

金虫 (小有名气)

★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
xiegangmai: 金币+2, 谢谢参与! 2013-03-24 14:32:09
曳力系数的例子中就有,你看一下就知道该如何处理了。
DEFINE_EXCHANGE_PROPERTY(custom_drag,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;

/* find the threads for the gas (primary) */
/* and solids (secondary phases)          */

thread_g = THREAD_SUB_THREAD(mix_thread, s_col);/* gas phase  */
thread_s = THREAD_SUB_THREAD(mix_thread, f_col);/* solid phase*/

/* find phase velocities and properties*/

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 Reynold's 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.281632*pow(void_g, 1.28);
else
  bfac = pow(void_g, 9.076960);

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楼2013-03-21 17:47:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

malu430

金虫 (小有名气)

★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
xiegangmai: 金币+1, 谢谢参与! 2013-03-24 14:32:20
vof模型cell中各相好像是共用一个速度场的吧 速度应该是一样的啊
3楼2013-03-22 14:46:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

陈wen

新虫 (初入文坛)

★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
xiegangmai: 金币+1, 谢谢参与! 2013-03-24 14:32:29
没必要用到两相之间的速度差,直接提取就可以了
4楼2013-03-22 20:08:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsqtcc

金虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
3楼: Originally posted by malu430 at 2013-03-22 14:46:19
vof模型cell中各相好像是共用一个速度场的吧 速度应该是一样的啊

用VOF模型也可以提取两相的速度差,用NV_DD(,,,)。
闪电侠
5楼2015-03-17 23:19:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsqtcc

金虫 (正式写手)

你好,请问问题解决了吗?
闪电侠
6楼2015-03-17 23:20:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jyl432 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 362求调剂 +3 西南交材料专硕3 2026-03-31 3/150 2026-03-31 10:39 by jp9609
[考研] 347求调剂 +8 山顶见α 2026-03-25 8/400 2026-03-31 10:14 by 不吃魚的貓
[考研] 看看我!301求调剂 081700 +3 一个小明同学 2026-03-24 3/150 2026-03-31 09:40 by longlotian
[考研] 308求调剂 +8 墨墨漠 2026-03-25 8/400 2026-03-31 09:15 by dgd2013
[考研] 哈尔滨工业大学材料与化工专硕378求调剂 +3 塔比乌斯 2026-03-30 3/150 2026-03-30 22:55 by 无际的草原
[考研] 085701环境工程求调剂 +11 多久上课 2026-03-27 12/600 2026-03-30 21:21 by 研究僧导导
[考研] 288资源与环境专硕求调剂,不限专业,有学上就行 +11 lllllos 2026-03-30 11/550 2026-03-30 21:19 by 研究僧导导
[考研] 一志愿华东师范大学有机化学专业,初试351分,复试被刷求调剂! +7 真名有冰 2026-03-29 8/400 2026-03-30 21:08 by maaj02
[考研] 277跪求调剂 +8 1915668 2026-03-27 12/600 2026-03-30 21:01 by dophin1985
[考研] 317求调剂 +8 十闲wx 2026-03-24 8/400 2026-03-30 21:00 by peike
[考研] 一志愿北京化工大学材料与化工(085600)296求调剂 +25 稻妻小编 2026-03-26 25/1250 2026-03-30 20:11 by 滴滴上岸呀
[考研] 0703化学321分求调剂 +10 三dd. 2026-03-30 11/550 2026-03-30 19:24 by markhwc
[考研] 332求调剂 +6 @MZB382400 2026-03-28 6/300 2026-03-30 16:57 by 无际的草原
[考研] 一志愿211,335分,0856,求调剂院校和导师 +7 倾____萧 2026-03-27 8/400 2026-03-30 09:37 by longlotian
[考研] 求调剂 +10 张zz111 2026-03-27 11/550 2026-03-30 09:17 by 无际的草原
[考研] 085600,材料与化工321分求调剂 +10 大馋小子 2026-03-28 10/500 2026-03-29 23:35 by 飞行日记西
[考研] 一志愿双一流机械285分求调剂 +4 幸运的三木 2026-03-29 5/250 2026-03-29 14:49 by Miko19
[考研] 356求调剂 +3 gysy?s?a 2026-03-28 3/150 2026-03-29 00:33 by 544594351
[考研] 315分求调剂 +7 26考研上岸版26 2026-03-26 7/350 2026-03-28 04:05 by fmesaito
[考研] 机械学硕310分,数一英一,一志愿211本科双非找调剂信息 +3 @357 2026-03-25 3/150 2026-03-26 16:34 by by.MENG
信息提示
请填处理意见