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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 流化床曳力系数UDF
101/1返回列表
查看: 2748  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

mokalin

新虫 (初入文坛)

[求助] 流化床曳力系数UDF

简单的液固两相流化床，用欧拉双流体模型，两相间的曳力系数UDF怎么编写，曳力系数公式已经有了，UDF不会编，有没有高手能给个类似的.c文件
回复此楼

» 猜你喜欢

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

47425804

木虫 (正式写手)

【答案】应助回帖

★
xiegangmai(金币+1): 谢谢参与！ 2011-10-18 17:55:21
UDF手册里有例子可以参考。
赞一下(1人) 回复此楼
2楼2011-10-18 15:10:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mokalin

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by 47425804 at 2011-10-18 15:10:47:
UDF手册里有例子可以参考。

能给个附件吗，链接也行，谢谢
赞一下 回复此楼
3楼2011-10-18 15:13:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

NaOHWTO

金虫 (小有名气)
★
xiegangmai(金币+1): 谢谢参与！ 2011-10-18 21:12:27
我有，自己编的，已经通过各种曳力模型计算，不过不会给出来的。按照HLEP里面的例子编吧，仔细看看可以搞定的。
赞一下(1人) 回复此楼
4楼2011-10-18 19:52:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

47425804

木虫 (正式写手)

【答案】应助回帖

★
mokalin(金币+5): 谢谢 2011-10-19 13:47:24
spc08(金币+1): 感谢应助 2011-10-19 15:18:29
引用回帖:
3楼: Originally posted by mokalin at 2011-10-18 15:13:21:
能给个附件吗，链接也行，谢谢

就是fluent自带的UDF手册，讲DEFINE_EXCHANGE_PROPERTY宏的那一节，有自定义曳力模型的例子。
我也给你贴上来吧
/******************************************************************
   UDF for customizing the default Syamlal drag law in ANSYS FLUENT  
*******************************************************************/
   
#include "udf.h"

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

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人) 回复此楼
5楼2011-10-19 13:21:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mokalin

新虫 (初入文坛)
引用回帖:
5楼: Originally posted by 47425804 at 2011-10-19 13:21:49:
就是fluent自带的UDF手册，讲DEFINE_EXCHANGE_PROPERTY宏的那一节，有自定义曳力模型的例子。
我也给你贴上来吧
/******************************************************************
   UDF for customizi ...

谢谢
回复此楼
6楼2011-10-19 13:47:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tacaro

至尊木虫 (知名作家)

小木虫大能猫骑士勋章
赞5楼~~~~~
回复此楼
7楼2011-10-20 09:36:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangxu86

铜虫 (小有名气)
引用回帖:
5楼: Originally posted by 47425804 at 2011-10-19 13:21:49:
就是fluent自带的UDF手册，讲DEFINE_EXCHANGE_PROPERTY宏的那一节，有自定义曳力模型的例子。
我也给你贴上来吧
/******************************************************************
   UDF for customizi ...

请问楼主，知道里面的参数0.281632和9.076960是怎么得到的吗？？
赞一下 回复此楼
8楼2011-10-23 15:33:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

47425804

木虫 (正式写手)
引用回帖:
8楼: Originally posted by zhangxu86 at 2011-10-23 15:33:04:
请问楼主，知道里面的参数0.281632和9.076960是怎么得到的吗？？

这个我也不知道，要看文献。
赞一下 回复此楼
9楼2011-10-24 08:50:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

一月的河

木虫 (小有名气)
引用回帖:
8楼: Originally posted by zhangxu86 at 2011-10-23 15:33:04
请问楼主，知道里面的参数0.281632和9.076960是怎么得到的吗？？...

我也只是在其他关于流化床的例子中看到，默认的常数是0.8和2.65，但是通过数学分析和处理，这些常数分别为0.281 632和9.076 96
赞一下 回复此楼
10楼2013-08-22 16:50:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mokalin 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 302求调剂 +4 呼呼呼。。。。 2026-03-17 4/200 2026-03-17 22:20 by 77Qi
[考研] 26考研求调剂 +6 丶宏Sir 2026-03-13 6/300 2026-03-17 16:13 by 醉在风里
[考研] 材料与化工专硕调剂 +5 heming3743 2026-03-16 5/250 2026-03-17 14:03 by 勇敢太监王公公
[考研] 0854可跨调剂，一作一项核心论文五项专利，省、国级证书40+数一英一287 +3 小李0854 2026-03-16 3/150 2026-03-17 13:40 by 热情沙漠
[考研] 一志愿南京大学，080500材料科学与工程，调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 283求调剂 +3 听风就是雨； 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 环境工程调剂 +6 大可digkids 2026-03-16 6/300 2026-03-16 17:16 by barlinike
[考研] 0703化学调剂，求各位老师收留 +8 秋有木北 2026-03-14 8/400 2026-03-16 15:21 by 哦哦123
[考研] 机械专硕调剂 +3 笨笨兔子 2026-03-12 3/150 2026-03-15 20:02 by 栗子粥?
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 331求调剂（0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 341求调剂 +3 番茄头--- 2026-03-10 3/150 2026-03-13 23:07 by JourneyLucky
[考研] 求调剂（材料与化工327） +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] （081700）化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
信息提示
关闭
请填处理意见
关闭