24小时热门版块排行榜    

查看: 4246  |  回复: 9

ddup885

银虫 (正式写手)

[求助] 利用DEFINE_EXCHANGE_PROPERTY编的蒸发和冷凝的udf,大侠们帮我看看对不对

#include "udf.h"                                                /*包括常规宏*/
#define  pi 3.1415                                            /*圆周率常数*/
#define  R 8.314                                             /*gas constant 8.314j/mol.K*/
#define  moleweight 0.018                                   /*SI unit water mole weight 0.018kg/mol */
#define A 7.19621
#define B 1730.63
#define C 233.426
                                    
DEFINE_EXCHANGE_PROPERTY(liq_gas_source, cell, thread, liq_i, gas_i)
{
  real m_lg;
  real P_SAT;
  Thread *gas = THREAD_SUB_THREAD(thread, gas_i);
  Thread *liq = THREAD_SUB_THREAD(thread, liq_i);

  P_SAT=exp(A-B/(C_T(cell,gas)+C-273.15))*1e3;     /*定义饱和蒸汽压与液相温度的关系式,即antoine公式*/
  m_lg = 0.;
if(C_P(cell,liq)   m_lg=(P_SAT-C_P(cell,liq))*sqrt(moleweight/(2*pi*R*C_T(cell,gas)));  /*蒸发的质量*/
if(C_P(cell,gas)>P_SAT)
  m_lg=(P_SAT-C_P(cell,gas))*sqrt(moleweight/(2*pi*R*C_T(cell,gas)));  /*凝结的质量*/

  return (m_lg);
}

上述程序是说关于水的蒸发和冷凝的,当水的压强小于饱和蒸汽压时,水分蒸发;实际压强大于饱和蒸汽压时,水汽凝结。
编完了之后不是很确定,有几个疑问
(1)DEFINE_EXCHANGE_PROPERTY宏是否将正的返回值传递给液相源相,而将负的返回值传递给气相源相?
(2)用了上述宏之后是否就不需要再编写两相源相的宏了呢?
(3)我引用的蒸发率的单位是kg/(s.m2),和连续性方程中的kg/(s.m3)不一样。可是蒸发率表示的是单位时间单位面积的蒸发质量,如何能和连续性方程中的单位变成一样的呢?
(4)在上述宏之外,是否还需要给能量加一个源相?

请和我研究方向相似的虫友多多帮忙,谢谢
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

ddup885

银虫 (正式写手)


ddup885: 回帖置顶 2012-08-29 14:56:55
xiegangmai: 金币+1, 谢谢参与 2012-10-10 21:20:39
原来可以直接用空穴模型 cavitation解决这种非沸腾蒸发。
4楼2012-08-29 14:56:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

fangdian

木虫 (正式写手)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
ddup885: 金币+2, 有帮助 2012-08-27 16:15:52
xiegangmai: 金币+1, 谢谢参与 2012-11-12 19:25:40
没有质量转变的源项?DEFINE_MASS_TRANSFER   
还有,能量转移呢?
2楼2012-08-24 17:05:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ddup885

银虫 (正式写手)

引用回帖:
2楼: Originally posted by fangdian at 2012-08-24 17:05:17
没有质量转变的源项?DEFINE_MASS_TRANSFER   
还有,能量转移呢?

我是利用相间的质量传递率定义的(phase interaction),看fluent帮助好像这个就是两相间质量的传递,所以我想是不是就可以替代质量源相。能量转移正在编,还没加进去。
3楼2012-08-24 17:10:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guolei8619

铁虫 (初入文坛)

【答案】应助回帖


xiegangmai: 金币-1, 谢谢参与,非应助请勿选择“应助回帖” 2012-10-10 21:20:54
你好,我是做平板热管模拟的,也在写蒸发冷凝的UDF。想请教一个问题。m_lg=(P_SAT-C_P(cell,liq))*sqrt(moleweight/(2*pi*R*C_T(cell,gas)));  /*蒸发的质量*/这个质量的蒸发量计算是采用的经验公式吗?
5楼2012-10-10 16:44:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuxia8495

银虫 (小有名气)

【答案】应助回帖

xiegangmai: 谢谢参与,非应助请勿选择“应助回帖” 2012-11-12 19:25:56
引用回帖:
5楼: Originally posted by guolei8619 at 2012-10-10 16:44:45
你好,我是做平板热管模拟的,也在写蒸发冷凝的UDF。想请教一个问题。m_lg=(P_SAT-C_P(cell,liq))*sqrt(moleweight/(2*pi*R*C_T(cell,gas)));  /*蒸发的质量*/这个质量的蒸发量计算是采用的经验公式吗?

类似于Langmuir方程形式,可又不是,我也想知道楼主为何如此选择??谢谢
fluent交流群:247890687
6楼2012-11-09 08:56:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuxia8495

银虫 (小有名气)

【答案】应助回帖

★ ★
xiegangmai: 金币+2, 鼓励讨论交流 2012-11-12 19:26:14
内容已删除
fluent交流群:247890687
7楼2012-11-09 09:13:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huizheng

铁虫 (初入文坛)

【答案】应助回帖

★ ★
xiegangmai: 金币+2, 谢谢参与 2013-01-08 21:06:20
lz的这个UDF用的是define_exchange_property,但是在uDF的帮助中定义相变传质用的是define_mass_transfer,并且我也试过了,在phase下面定义相变传质的时候define_exchange_property是不对的
8楼2013-01-07 16:50:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sogoodaaa

银虫 (正式写手)

【答案】应助回帖

★ ★
xiegangmai: 金币+2, 鼓励讨论交流! 2013-06-07 20:49:03
我针对楼主提出的第三个问题提点看法哈,这个问题也是我疑惑好久的,就是一般计算冷凝速率,其量纲是质量/(时间。面积),就是计算在一定面积上的冷凝量,但是有时冷凝的面积不容易测量,这时就会采用容积传质系数来计算通量,办法是再乘以一个比表面积,即质量/(时间*面积)*面积/体积,这样量纲就一致了。
开心的过每一天
9楼2013-06-07 12:39:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ailianyi99

新虫 (初入文坛)

10楼2015-05-20 15:32:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ddup885 的主题更新
信息提示
请填处理意见