24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

死天

专家顾问 (著名写手)

专家经验: +60
应助: 399 (硕士)
金币: 3648.9
红花: 89
帖子: 1336
在线: 443.8小时
虫号: 2643057
注册: 2013-09-09
专业: 流体力学
管辖: 仿真模拟

[交流] fluent Degassing Boundary（脱气边界UDF）已有7人参与

RT。想要的就自己下吧。具体说明见附件

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : UDF.pptx

2015-03-31 13:27:19, 267.78 K

» 本帖已获得的红花（最新10朵）

天堂鱼520

88zhanchi88

» 猜你喜欢

基金申报已经有5人回复
基金委咋了？2026年的指南还没有出来？已经有7人回复
国自然申请面上模板最新2026版出了吗？已经有17人回复
纳米粒子粒径的测量已经有8人回复
疑惑？已经有5人回复
计算机、0854电子信息（085401-058412）调剂已经有5人回复
Materials Today Chemistry审稿周期已经有5人回复
溴的反应液脱色已经有7人回复
推荐一本书已经有12人回复
常年博士招收(双一流，工科) 已经有4人回复

1楼 2015-03-31 13:28:07

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

88zhanchi88

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 149.9
散金: 217
红花: 3
帖子: 65
在线: 118.9小时
虫号: 2569661
注册: 2013-07-28
性别: GG
专业: 工程热物理相关交叉领域

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

送红花一朵

您好，我在做旋转填料床的模拟，强化气液相的接触，边界条件需要定义为气体进口，液滴出口。
因此我想用 Degassing Boundary的udf来在出口边界上定义一个源项。
我想问两个问题：
1.由于源项只能加载到cell zone 里，因此通过什么语句找到出口边界的位置？
  以下help文件中的例子可以直接用吗？需要修改么？
2.以下几个源项分别加载到哪些相里呢？混合相还是液相还是气相？

#include "udf.h"
#include "sg.h"
#include "sg_mphase.h"
#include "flow.h"
#include "mem.h"
#include "metric.h"

DEFINE_SOURCE(degassing_source, cell, thread, dS, eqn)
{
real source;
Thread *tm = THREAD_SUPER_THREAD(thread);
source = -C_R(cell,thread)*C_VOF(cell,thread)/CURRENT_TIMESTEP;
C_UDMI(cell,tm,0) = source;
dS[eqn] = -C_R(cell,thread)/CURRENT_TIMESTEP;
return source;
}

DEFINE_SOURCE(x_prim_recoil, cell, tp, dS, eqn)
{
real source;
Thread *tm = THREAD_SUPER_THREAD(tp);
Thread *ts;
ts = THREAD_SUB_THREAD(tm,1);
source = -C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP*C_U(cell,tp);
dS[eqn] =-C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP;
return source;
}

DEFINE_SOURCE(x_sec_recoil, cell, ts, dS, eqn)
{
real source;
Thread *tm = THREAD_SUPER_THREAD(ts);
source = -C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP*C_U(cell,ts);
dS[eqn] = -C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP;
return source;
}

DEFINE_SOURCE(y_prim_recoil, cell, tp, dS, eqn)
{
real source;
Thread *tm = THREAD_SUPER_THREAD(tp);
Thread *ts;
ts = THREAD_SUB_THREAD(tm,1);
source = -C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP*C_V(cell,tp);
dS[eqn] =-C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP;
return source;
}

  DEFINE_SOURCE(y_sec_recoil, cell, ts, dS, eqn)
{
real source;  Thread *tm = THREAD_SUPER_THREAD(ts);
source = -C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP*C_V(cell,ts);
dS[eqn] = -C_R(cell,ts)*C_VOF(cell,ts)/CURRENT_TIMESTEP;
return source;
}