24小时热门版块排行榜    

查看: 2310  |  回复: 9

luckyhuman1

铜虫 (小有名气)

[求助] fluent 怎样模拟壁面吸收组分

我的模型采用组分运输,计算的是液体和液体,没有反应,在壁面要消耗掉组分,这个过程要怎样模拟啊?是在壁面上定义消耗掉的源项吗?还是采用其他方法?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kylafree

至尊木虫 (知名作家)

【答案】应助回帖


感谢参与,应助指数 +1
xiegangmai(金币+1): 谢谢应助! 2012-03-11 16:20:42
可以通过在壁面附近定义消耗掉的组分源项来实现
2楼2012-03-11 15:44:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luckyhuman1

铜虫 (小有名气)

引用回帖:
: Originally posted by kylafree at 2012-03-11 15:44:30:
可以通过在壁面附近定义消耗掉的组分源项来实现

您好,我是这样写的一段UDF,但是加进去之后的计算速度太慢了,我一个晚上没管才迭代了几次
您能不能给我看看这个有什么问题,
#include "udf.h"
#define Coating_weight 0.06
#define Coat_Al    0.004
#define Strip_wide  0.5
#define Strip_velocity  2.3
DEFINE_SOURCE(Al_source,c,t,dS,eqn)
{
real source;
real CC[ND_ND];
int zone_ID = 5;
real volume=0.0;
real y;
real area;
Thread *inletstrip;
Thread *t0;

y=0.9-0.2*Strip_velocity*cos(33.5);
area=0.2*Strip_velocity*Strip_wide;

inletstrip = Lookup_Thread(Get_Domain(1),zone_ID);
t0=THREAD_T0(inletstrip);

begin_c_loop(c,t0)
{
  C_CENTROID(CC,c,t0);
C_UDMI(c,t0,0)=CC[1];

if (CC[1]>=y&&CC[1]<=0.9)
volume+=C_VOLUME(c,t0);
}
end_c_loop(c,t0)
C_UDMI(c,t0,1)=volume;

if (C_UDMI(c,t0,0)>=y&&C_UDMI(c,t0,0)<=0.9)
{source=-5*Coating_weight*Coat_Al*area/volume;
C_UDMI(c,t0,2)=source;
}
else source=0;
return source;
}

这个要在壁面的上的一段消耗,其他地方不消耗,所以有个判断位置的语句,然后基本的消耗源是个常数,麻烦您给看看啦,都改好久也不知道怎么改了,不要体积的求和的时候计算还比较快,加了那一段就没法算了
3楼2012-03-14 10:30:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lbn045

金虫 (小有名气)

【答案】应助回帖

DEFINE_SOURCE默认就是在cell中循环 我觉得你把循环的语???去掉会??会好些
4楼2013-05-30 02:11:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lbn045

金虫 (小有名气)

【答案】应助回帖


xiegangmai: 金币+1, 谢谢参与! 2013-06-02 10:39:03
DEFINE_SOURCE 默认在cellz中循环 所以我认为你把循环语句去掉 会不会好些
5楼2013-05-30 02:13:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lbn045

金虫 (小有名气)

请问Coating_weight,Coat_Al,area的物理意义都是什么?我也想编一个类似的源项,源项的单位是kg/m3.s,要是非稳态我还能弄懂,但是稳态的就不会弄了。所以想请教你一下,呵呵
6楼2013-06-28 01:58:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (著名写手)

本帖仅楼主可见
7楼2015-04-07 17:14:58
已阅   申请仿真EPI   回复此楼   编辑   查看我的主页

星星在我家门

木虫 (著名写手)

热爱互联网

有没有人可以指导一下,我写的define_profile改变边界处组分分数?可以实现同样的功能吗?
于无声听惊雷
8楼2018-04-27 20:50:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

low_w

新虫 (小有名气)

引用回帖:
2楼: Originally posted by kylafree at 2012-03-11 15:44:30
可以通过在壁面附近定义消耗掉的组分源项来实现

您好,请问,可否具体说一下如何确定其是壁面附近的单元,进而进行源项的定义,小白求指导,必重谢!
9楼2018-09-28 22:43:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

low_w

新虫 (小有名气)

引用回帖:
3楼: Originally posted by luckyhuman1 at 2012-03-14 10:30:02
您好,我是这样写的一段UDF,但是加进去之后的计算速度太慢了,我一个晚上没管才迭代了几次
您能不能给我看看这个有什么问题,
#include "udf.h"
#define Coating_weight 0.06
#define Coat_Al    0 ...

您好,请问怎么确定该网格是近壁面处的网格呀?万分感谢!
10楼2018-09-28 22:45:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 luckyhuman1 的主题更新
信息提示
请填处理意见