应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 曳力函数的的UDF
161/2返回列表12下一页
查看: 3496  |  回复: 15
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者zhangxu86将赠送您 5 个金币

zhangxu86

铜虫 (小有名气)

[求助] 曳力函数的的UDF已有2人参与

请问各位大侠，我的模型使用的是欧拉模型，物质是沙和水，请问一下，编写沙水两者之间的曳力时，是使用液液交换系数还是使用液体-固体交换系数呢？下面的是我看到一篇流化床立体上的程序，不知道为什么计算两次K和f呢？上面一段程序使用的就是Syamlal-O‘Brien，下面一段就不明白了。。。。。
/* This is a routine for customizing default Syamlal drag law in Fluent 6.
   The default drag law uses 0.8 (for void<=0.85) and 2.65 (void>0.85) for
   bfac. This is for a min fluid vel of 25 cm/s. The current drag law has been
   tuned for a min fluid vel of 8 cm/s and uses 0.28 and 9.07 for these
   parameters. */

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

# define pi 4.*atan(1.)
#define diam2 3.e-4

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;

/* find the threads for the gas (primary) and solids (secondary phases).
These phases appear in columns 2 and 1 in the Interphase panel respectively*/

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 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.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;

}

DEFINE_EXCHANGE_PROPERTY(custom_drag_ihme, 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, cd, eg,
      void_g, k_g_s;

/* find the threads for the gas (primary) and solids (secondary phases).
These phases appear in columns 2 and 1 in the Interphase panel respectively*/

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 reynolds number*/

reyp = rho_g*abs_v*diam2/mu_g;

cd = (24./(reyp+SMALL)) + 5.48*pow((reyp+SMALL),-0.573) + 0.36;

void_g = C_VOF(cell, thread_g);/* gas vol frac*/

eg = pow(void_g,-2.65);

k_g_s = (3./4.)*(cd*(1.-void_g)*abs_v*rho_g*eg)/diam2;

return k_g_s;

}
回复此楼

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

我的收藏

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2011-10-25 21:29:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

47425804

木虫 (正式写手)
★
臭水沟(金币+1): 谢谢交流~~ 2011-10-26 10:07:31
计算两次K和f？
不明白什么意思。
哪儿计算两次K和f了？
赞一下(1人) 回复此楼
2楼2011-10-25 22:52:01
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangxu86

铜虫 (小有名气)
引用回帖:
2楼: Originally posted by 47425804 at 2011-10-25 22:52:01:
计算两次K和f？
不明白什么意思。
哪儿计算两次K和f了？

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

k_g_s = (3./4.)*(cd*(1.-void_g)*abs_v*rho_g*eg)/diam2;
是两次计算交换系数K值，不知道下面的K值是干什么用的呢？
赞一下 回复此楼
3楼2011-10-26 09:12:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

47425804

木虫 (正式写手)

【答案】应助回帖

★
臭水沟(金币+1): 谢谢交流~~ 2011-10-27 16:10:28
引用回帖:
3楼: Originally posted by zhangxu86 at 2011-10-26 09:12:50:
k_g_s = (1.-void_g)*rho_s*fdrgs/taup;

k_g_s = (3./4.)*(cd*(1.-void_g)*abs_v*rho_g*eg)/diam2;
是两次计算交换系数K值，不知道下面的K值是干什么用的呢？

没有计算两遍，这是两个函数，也就是两个曳力模型，互相没有关系的。
赞一下(1人) 回复此楼
4楼2011-10-26 10:13:08
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangxu86

铜虫 (小有名气)
引用回帖:
4楼: Originally posted by 47425804 at 2011-10-26 10:13:08:
没有计算两遍，这是两个函数，也就是两个曳力模型，互相没有关系的。

不是一个曳力就可以了吗，一个模型怎么可以选择两个曳力模型呢？这个一个流化床的案例，我觉得只要一个就可以了啊
赞一下 回复此楼
5楼2011-10-26 19:30:10
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

47425804

木虫 (正式写手)

【答案】应助回帖

★ ★
臭水沟(金币+2): 谢谢应助~~ 2011-10-27 16:10:41
引用回帖:
5楼: Originally posted by zhangxu86 at 2011-10-26 19:30:10:
不是一个曳力就可以了吗，一个模型怎么可以选择两个曳力模型呢？这个一个流化床的案例，我觉得只要一个就可以了啊

那要看具体的cas里是怎么用的了
说不定有好几种颗粒，分别用了各自的曳力模型；
也说不定就是写了两个放在里面，图方便而已。
赞一下(1人) 回复此楼
6楼2011-10-26 22:06:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangshuaibyx

金虫 (小有名气)
引用回帖:
6楼: Originally posted by 47425804 at 2011-10-26 22:06:31
那要看具体的cas里是怎么用的了
说不定有好几种颗粒，分别用了各自的曳力模型；
也说不定就是写了两个放在里面，图方便而已。...

就是气固两相，能用两个曳力模型吗？是不是只能用一个，上面程序中有两个但是可能有一个用不到。
赞一下 回复此楼
7楼2014-06-21 21:45:14
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ajige123

禁虫 (初入文坛)
本帖内容被屏蔽
8楼2014-11-27 15:14:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liubao2012

木虫 (小有名气)
引用回帖:
8楼: Originally posted by ajige123 at 2014-11-27 15:14:36
这是流化床中最简单的UFD了，气固流化选择曳力的时候就是两个，custom_drag_syam是气体的曳力系数，custom_drag_ihme是固体的曳力系数，楼主还没做过两维流化床的最简单的CASE吧，要不然怎么没有选择过曳力系数呢？...

为什么气相和固相都会有曳力系数啊？曳力系数本身不就是气固之间的唯一一个作用系数吗？
那是不是所有的曳力模型里面都要有DEFINE_EXCHANGE_PROPERTY(custom_drag_ihme, cell, mix_thread, s_col, f_col)
和DEFINE_EXCHANGE_PROPERTY(custom_drag_ihme, cell, mix_thread, s_col, f_col)啊？
赞一下 回复此楼
9楼2015-12-28 22:14:11
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ajige123

禁虫 (初入文坛)
本帖内容被屏蔽
10楼2015-12-29 10:18:27
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhangxu86 的主题更新
161/2返回列表12下一页
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国家资助计划 +5 bearbo 2024-07-01 6/300 2024-07-02 08:28 by litleprince
[基金申请] 史上最严的基金限项 +8 MathGLA 2024-07-01 8/400 2024-07-02 08:16 by Howard28
[找工作] 博士有找工作搭子吗？ +11 如花美卷18 2024-06-30 14/700 2024-07-02 01:41 by 专业控
[基金申请] 咨询 +5 李长云 2024-06-30 6/300 2024-07-02 00:22 by caiqing
[论文投稿] 论文润色 +6 自强不息a?a 2024-06-26 7/350 2024-07-01 21:57 by 自强不息a?a
[基金申请] 2024上会名单已出 +8 dream200606 2024-06-28 13/650 2024-07-01 13:51 by dream200606
[访问学者] 放榜是先短信还是先系统？ +8 wyjecho666 2024-06-29 10/500 2024-07-01 10:54 by maolC
[论文投稿] 两个月二审还没完成 +7 建议Accept 2024-06-29 11/550 2024-07-01 10:20 by 建议Accept
[论文投稿] 投稿求助 54+3 s090604054 2024-06-30 4/200 2024-07-01 10:11 by 恕在下直言
[基金申请] 中了博后面上系统能查到分数和排名吗？ +5 老黑不黑 2024-06-29 6/300 2024-07-01 06:19 by yqliu130
[基金申请] 是否上会都不知道的三无人员，翻论坛想看会评消息 +7 sparknow 2024-06-27 8/400 2024-07-01 00:19 by 临川树荫
[访问学者] no news is good news +8 wyjecho666 2024-06-29 12/600 2024-06-30 13:49 by chang_jiang
[基金申请] 博后面上如何提高创新分？ +5 yuyiang 2024-06-29 10/500 2024-06-29 18:24 by yuyiang
[有机交流] 大佬们，打的核磁氢谱与chemdraw预测的有偏差 +8 来了个 2024-06-26 9/450 2024-06-29 17:58 by if今年的吧
[硕博家园] 考博英语 +8 198新用户 2024-06-25 15/750 2024-06-29 15:02 by Lucaswu
[硕博家园] 有机合成硕士毕业8年，想读博士，但又无法全职 +5 青灯孤案 2024-06-26 7/350 2024-06-29 12:16 by wod河畔
[论文投稿] Scitific Reports +5 智晨君 2024-06-27 5/250 2024-06-29 10:37 by zhulei513
[基金申请] 2024年度国家资助博士后研究人员计划（B、C 档）也不知道什么时候出 +11 SCU船船船 2024-06-28 15/750 2024-06-29 00:21 by 666中
[考博] 申博找导师 +4 疏影横斜水清浅3 2024-06-26 5/250 2024-06-28 11:15 by 安塔瓦拉多
[有机交流] 酚羟基甲基化 10+4 A好运来啦啦啦 2024-06-25 6/300 2024-06-26 10:08 by 88817753
信息提示
关闭
请填处理意见
关闭