找回密码
 注册
查看: 9606|回复: 15

小妹对水蒸发udf 的疑问,想请教各位大虾!

[复制链接]
发表于 2009-3-16 15:34:14 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?注册

x
各位大虾,我做的是水遇高温气体蒸发成水蒸汽的问题,用的是mixture模型,要把水和水蒸汽之间的传质加入到质量源项,我在网上找到了一段有关水蒸发的udf,其中的温度项是用C_T(cell,mix_thread),但fluent 的udf 帮助中一个相似的例子中对应的温度项却是用C_T(cell,liq_thread)来表示的,请问各位路过的大虾,这个温度到底应该是用液相温度C_T(cell,liq_thread)表示,还是用混合相温度C_T(cell,mix_thread)?网上的udf如下:(请留意笑脸处)
DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn) ‘液相质量源项UDF

{

Thread *mix_th, *sec_th;  ‘定义计算区线指针

  real m_dot_l; ‘定义液相质量转移 kg/(m2.s)

  mix_th = THREAD_SUPER_THREAD(pri_th);  ‘指向混合区的主相即液相的指针

  sec_th = THREAD_SUB_THREAD(mix_th, 1);指向单相控制区的气相的指针,气相为第二相

  if(C_T(cell, mix_th) >=T_SAT)  ‘如果液相单元的温度高于蒸发温度,液相向气相的质量质量转移

{

  m_dot_l = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*   

  fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;

    dS[eqn] = -0.1*C_R(cell, pri_th)*fabs(C_T(cell, mix_th)  - T_SAT)/T_SAT;

              ‘ 定义源项对质量转移偏导 ,



                               }

  else {

    m_dot_l = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*

                fabs(T_SAT-C_T(cell,mix_th):) )/T_SAT;

‘如果指向混合区液相的单元温度小于蒸发温度,气相向液相的质量转移,液相得



    dS[eqn] = 0.;  ‘由于是气相向液相转移,所以液相的质量源项对质量转移的偏导为零



       }

   

  return m_dot_l;

}
fluent帮助中的udf:
if (C_T(cell, liq):) )  >= T_SAT)
     {
       m_lg = -0.1*C_VOF(cell,liq)*C_R(cell,liq)*
               fabs(C_T(cell,liq):)  -T_SAT)/T_SAT;

[ 本帖最后由 h04581803 于 2009-3-16 07:36 编辑 ]
 楼主| 发表于 2009-3-16 15:37:28 | 显示全部楼层
终于发送成功了,自己先顶一个,期待各位高手的指点!!!!
 楼主| 发表于 2009-3-17 09:50:08 | 显示全部楼层
gallongallon版主帮我看下拉,很着急,谢谢了,版主帮忙!
 楼主| 发表于 2009-3-18 13:47:52 | 显示全部楼层
怎么问的问题总没人回答啊,
发表于 2009-3-24 09:00:04 | 显示全部楼层
我觉得这个问题用温度差做推动力很不对,难道没到饱和温度就不蒸发了?又不是沸腾。
我在想是不是按浓度差 的 差来推动,表现出来 就是按照饱和蒸汽压和气体压力的差值来做推动?
 楼主| 发表于 2009-3-27 09:55:09 | 显示全部楼层
我只想算最后有多少水变成水蒸汽,不考虑力的因素,而且计算中水的压力是不会降到饱和蒸汽压以下的,楼上的可否解释下,我不是太明白你的意思。。。。。
发表于 2009-3-31 09:29:41 | 显示全部楼层

回复 6# h04581803 的帖子

我也是做一个类似的问题, 但是还没有把模型建出来。你的弄出来了?
先说你的这个东西,压力不降到饱和压以下,怎么蒸发?
蒸发的推动力不应该是(饱和蒸汽压-气相压力),比如水常温蒸汽压3000左右,大气里气相水压力。。。很小吧,设为a,然后推动力就是(3000-a),然后水沸腾了就是(100000-a)了,就更有利于蒸发,当然沸腾还有别的扰动啥的。

反正目前为止我的课题我的是没有弄出来。
9974055,加我q,讨论讨论。
 楼主| 发表于 2009-4-2 16:09:30 | 显示全部楼层
我想我们做的不是一个东西,我做的是喷射,喷射点压力是很高的,不可能会降到水的饱和蒸汽压,也没必要考虑,谢谢楼上的指点!
发表于 2010-11-21 16:07:02 | 显示全部楼层

回复 1# h04581803 的帖子

你在网上找的那个程序的第5行
mix_th = THREAD_SUPER_THREAD(pri_th);  ‘指向混合区的主相即液相的指针
即已经定义过mix_th是指向液相的,所以这个和帮助中的UDF其实是一样的哦。
发表于 2010-11-23 08:22:52 | 显示全部楼层

回复 6# h04581803 的帖子

你只想算多少水蒸发了的话,那你用这个udf肯定是不对的。
发表于 2012-9-18 13:23:33 | 显示全部楼层

回复 10# 浪子轻狂 的帖子

那要怎样算水的蒸发量,就是可以写出来能表示出水蒸发量的UDF吗?
发表于 2012-9-18 17:04:45 | 显示全部楼层
网上流传的这个水蒸发的UDF真是害人不浅。首先,收敛性就估计有问题。
发表于 2013-3-1 10:29:19 | 显示全部楼层

回复 12# gearboy78 的帖子

那水遇高温气体蒸发成水蒸汽的udf要怎样编啊?
发表于 2013-3-1 10:30:55 | 显示全部楼层

回复 11# shmily默默21 的帖子

水遇高温气体蒸发成水蒸汽的udf要怎样编?还有的就是mixture模型应该怎样设置,要设为三相吗?
发表于 2013-5-19 20:59:54 | 显示全部楼层

回复 1# h04581803 的帖子

你的问题解决了吗?我遇到了同样的问题,用的也是这个UDF
您需要登录后才可以回帖 登录 | 注册

本版积分规则

快速回复 返回顶部 返回列表