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

24小时热门版块排行榜     石溪大学接受考研调剂申请>

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » fluent沸腾两相流里，关于仅发生在气液相界面传质问题UDF求助！
【调剂】北京石油化工学院2024年16个专业接受调剂
41/1返回列表
查看: 2552  |  回复: 3
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

blueberryzy

新虫 (初入文坛)

[交流] fluent沸腾两相流里，关于仅发生在气液相界面传质问题UDF求助！已有3人参与

各位老师，求助一个关于fluent udf的问题，现在我想根据截图里的这个公式，通过udf添加一个描述气液相界面传热传质的模型，下面是我编写的UDF：
/***********************************/
#include "udf.h"
#include "sg.h"
#include "sg_mphase.h"
#include "math.h"
#include "stdio.h"
#include  "mem.h"/*包含梯度的头文件*/
#include "metric.h"
#include "flow.h"
#include "prop.h"

#define T_SAT 20.3
#define LAT_HT 447

cell_t c;

DEFINE_ADJUST(store_VOF_gradient, domain)/*将体积分数梯度储存在UDMI中*/
{
Thread *t;
Thread *ppt;
Thread **pt;

int phase_domain_index=1;
Domain *pDomain = DOMAIN_SUB_DOMAIN(domain,phase_domain_index);
Alloc_Storage_Vars(pDomain,SV_VOF_RG,SV_VOF_G,SV_NULL);
Scalar_Reconstruction(pDomain, SV_VOF,-1,SV_VOF_RG,NULL);
Scalar_Derivatives(pDomain,SV_VOF,-1,SV_VOF_G,SV_VOF_RG, Vof_Deriv_Accumulate);

mp_thread_loop_c (t,domain,pt)
{
if (FLUID_THREAD_P(t))
{
ppt = pt[phase_domain_index];
begin_c_loop (c,t)
{
C_UDMI(c,t,0) = NV_MAG(C_VOF_G(c,ppt));
C_UDMI(c,t,1) = C_VOF_G(c,ppt)[0];
C_UDMI(c,t,2) = C_VOF_G(c,ppt)[1];
//C_UDMI(c,t,3) = C_VOF_G(c,ppt)[2];
}
end_c_loop (c,t)
}
}
Free_Storage_Vars(pDomain,SV_VOF_RG,SV_VOF_G,SV_NULL);
}

/*定义质量与能量的源项*/
DEFINE_SOURCE(liq_src,cell,pri_th,dS,eqn)/*液相质量源项*/
{
Thread *mix_th, *sec_th;
real m_dot_l;
real T_X, T_Y;
real keff;

mix_th = THREAD_SUPER_THREAD(pri_th);

sec_th = THREAD_SUB_THREAD(mix_th, 1);

keff = C_VOF(cell,pri_th)*C_K_L(cell,pri_th)+C_VOF(cell,sec_th)*C_K_L(cell,sec_th);

if (NULL != THREAD_STORAGE(pri_th,SV_T_G))
{
T_X = C_T_G(cell,mix_th)[0];
}
if (NULL != THREAD_STORAGE(pri_th,SV_T_G))
{
T_Y = C_T_G(cell,mix_th)[1];
}

if(C_VOF(cell,sec_th) >0.05&&C_VOF(cell,sec_th) <1&&C_T(cell,mix_th) >= T_SAT)
{
m_dot_l = -keff*(T_X*C_UDMI(c,pri_th,1)+T_Y*C_UDMI(c,pri_th,2))/LAT_HT;
dS[eqn] = -C_K_L(cell,pri_th)*(T_X*C_UDMI(c,pri_th,1)+T_Y*C_UDMI(c,pri_th,2))/LAT_HT;
}
else
{
m_dot_l = keff*(T_X*C_UDMI(c,pri_th,1)+T_Y*C_UDMI(c,pri_th,2))/LAT_HT;
dS[eqn] = C_K_L(cell,sec_th)*(T_X*C_UDMI(c,pri_th,1)+T_Y*C_UDMI(c,pri_th,2))/LAT_HT;
}

return m_dot_l;

}

DEFINE_SOURCE(vap_src,cell,sec_th,dS,eqn)
{
Thread *mix_th, *pri_th;
real m_dot_v;
real T_X, T_Y;
real keff;

mix_th = THREAD_SUPER_THREAD(sec_th);

pri_th = THREAD_SUB_THREAD(mix_th,0);

keff = C_VOF(cell,pri_th)*C_K_L(cell,pri_th)+C_VOF(cell,sec_th)*C_K_L(cell,sec_th);

if (NULL != THREAD_STORAGE(pri_th,SV_T_G))
{
T_X = C_T_G(cell,mix_th)[0];
}
if (NULL != THREAD_STORAGE(pri_th,SV_T_G))
{
T_Y = C_T_G(cell,mix_th)[1];
}

if(C_VOF(cell,sec_th) >0.05&&C_VOF(cell,sec_th) <1 &&C_T(cell,mix_th) >= T_SAT)
{
m_dot_v = keff*(T_X*C_UDMI(c,pri_th,1)+T_Y*C_UDMI(c,pri_th,2))/LAT_HT;
dS[eqn] = C_K_L(cell,sec_th)*(T_X*C_UDMI(c,pri_th,1)+T_Y*C_UDMI(c,pri_th,2))/LAT_HT;
}
else
{
m_dot_v = -keff*(T_X*C_UDMI(c,pri_th,1)+T_Y*C_UDMI(c,pri_th,2))/LAT_HT;
dS[eqn] = -C_K_L(cell,pri_th)*(T_X*C_UDMI(c,pri_th,1)+T_Y*C_UDMI(c,pri_th,2))/LAT_HT;
}

return m_dot_v;

}
DEFINE_SOURCE(enrg_src,cell,mix_th,dS,eqn)
{
Thread *pri_th, *sec_th;
real m_dot;
real T_X, T_Y;
real keff;

pri_th = THREAD_SUB_THREAD(mix_th, 0);

sec_th = THREAD_SUB_THREAD(mix_th, 1);

keff = C_VOF(cell,pri_th)*C_K_L(cell,pri_th)+C_VOF(cell,sec_th)*C_K_L(cell,sec_th);

if (NULL != THREAD_STORAGE(pri_th,SV_T_G))
{
T_X = C_T_G(cell,mix_th)[0];
}
if (NULL != THREAD_STORAGE(pri_th,SV_T_G))
{
T_Y = C_T_G(cell,mix_th)[1];
}

if(C_VOF(cell,sec_th) >0.05&&C_VOF(cell,sec_th) <1 &&C_T(cell,mix_th) >= T_SAT)
{
m_dot = -keff*(T_X*C_UDMI(c,pri_th,1)+T_Y*C_UDMI(c,pri_th,2))/LAT_HT;
dS[eqn]= 0;
}
else
{
m_dot =keff*(T_X*C_UDMI(c,pri_th,1)+T_Y*C_UDMI(c,pri_th,2))/LAT_HT;
dS[eqn]=0;
}

return LAT_HT*m_dot*1000;
/***********************************/
在fluent里模拟后，UDF可以编译，但是计算过程发现没有相变发生，有没有老师可以帮忙看一下UDF?谢谢！

fluent沸腾两相流里，关于仅发生在气液相界面传质问题UDF求助！
截图1.jpg
回复此楼

» 猜你喜欢
1楼2018-10-23 09:54:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

★
小木虫: 金币+0.5, 给个红包，谢谢回帖
本帖仅楼主可见
赞一下
2楼2019-04-29 16:43:16
 已阅   申请仿真EPI   回复此楼   编辑   查看我的主页

344127102

铁虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
2楼: Originally posted by lijicheng009 at 2019-04-29 16:43:16
您好！

金虫
我也利用fluent做相变模拟，可是我利用fluent自带的蒸发-冷凝模型计算的时候，我想要的是自由液面蒸发的结果，可算一段时间后，液相里面也产生了气泡，这显然是不合理的，不知道这个模型能不能算纯 ...

我们也遇到了同样的问题，不知道你们解决了吗
赞一下 回复此楼
3楼2020-02-17 23:58:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1757761207

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
3楼: Originally posted by 344127102 at 2020-02-17 23:58:50
我们也遇到了同样的问题，不知道你们解决了吗...

请问你的问题解决的怎么样了？

发自小木虫Android客户端
赞一下 回复此楼
4楼2020-12-27 15:06:27
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 blueberryzy 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件) 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 真的好想读博！ +7 wangzhe_bs 2024-04-22 9/450 2024-04-23 22:20 by wangzhe_bs
[找工作] 工作问题求助 +3 向前?@ 2024-04-18 5/250 2024-04-23 21:25 by 向前?@
[硕博家园] 博士白读了 +33 Da_Meng_Zi 2024-04-21 36/1800 2024-04-23 17:22 by 李_q_y
[有机交流] 紧急求助，有谁用过三甲基硅醇钾TMSOK吗？ 1000+3 genius2008 2024-04-19 8/400 2024-04-23 16:57 by czyzsu
[硕博家园] 考研，求职还是考编? +10 xizj 2024-04-21 18/900 2024-04-23 16:37 by 别吵i
[考博] 博士申请 +4 zhangxm0525 2024-04-17 12/600 2024-04-23 16:35 by 莹亮的星空
[考研] 本科211 277分求调剂 +3 0854通信工程 2024-04-18 6/300 2024-04-23 15:39 by H.W.Y
[考博] 申博成果界定是根据Jcr分区还是中科院分区 +4 我属驴核动力驴 2024-04-22 4/200 2024-04-23 12:18 by lyfbangong
[基金申请] 颜宁:基础研究应顶天立地,立地即从需求端作为问题的出发点,立项背景就是面向真实需求 +10 zju2000 2024-04-18 10/500 2024-04-23 11:46 by 稻花香
[教师之家] 国内科研环境差导致国内科研水平差的根源 +8 zju2000 2024-04-18 8/400 2024-04-23 10:40 by Quakerbird
[基金申请] 特别资助审核状态 +12 pantray 2024-04-17 18/900 2024-04-23 09:50 by y20150913
[论文投稿] Chemical Engineering Journal投稿3周了，一直显示With editor状态。这是送审了吗？ 10+3 yifeng11 2024-04-20 7/350 2024-04-23 09:29 by bobvan
[考博] 研二光催化6月底4篇2区 +7 wjtab 2024-04-22 11/550 2024-04-23 06:59 by byron2012
[考博] 25申博记录贴 +6 我属驴核动力驴 2024-04-18 12/600 2024-04-22 15:50 by 诩雾
[考博] 青岛科技大学顾玉超教授课题组招收博士生 +3 yanqingwang 2024-04-19 3/150 2024-04-22 15:21 by 理科生yes
[考研] 312求调剂 +3 Leroic000 2024-04-19 3/150 2024-04-20 18:19 by xinaishan
[基金申请] 估计今年青基又没戏 +11 忆念7 2024-04-18 11/550 2024-04-20 17:29 by mygod_nju
[有机交流] 紫外光谱 50+3 54胡 2024-04-17 3/150 2024-04-19 10:48 by Nanamiwww
[论文投稿] CCS Chemistry投稿求助 30+3 wfqtriumph 2024-04-17 4/200 2024-04-19 08:56 by Bletilla
[论文投稿] 发个文章，结果被某学报请的外审PUA了 +3 ziniu168 2024-04-18 8/400 2024-04-18 21:40 by xli1984
信息提示
关闭
请填处理意见
关闭