phil 发表于 2007-6-19 19:17

VOF怎么模拟蒸发和冷凝啊？

:'( 我模拟的液气有明显界面，而且表面张力起主要作用，但是涉及到相变的饱和蒸汽压等数据在vof里没有，换其他的模型表面张力又不能涉及了，请问该怎么办呢？

linda 发表于 2007-6-20 00:06

没做过类似的问题，不过前一段时间见过一个沸腾方面的UDF，看看对你是否有帮助


#include "udf.h"
#include "sg_mphase.h"
#define T_SAT 373 /＊定义饱和温度*/
#define LAT_HT 1.e3 /＊定义水蒸汽潜热*/

DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn) /＊定义液相源项*/
{

Thread *mix_th, *sec_th; /＊混合相、第二相*/
real m_dot_l; /＊蒸发冷凝速率*/

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, pri_th) - T_SAT)/T_SAT;

   dS = -0.1*C_R(cell, pri_th)*
      fabs(C_T(cell, pri_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 = 0.;

       }
   
return m_dot_l;

}

DEFINE_SOURCE(vap_src, cell, sec_th, dS, eqn) /＊定义气相源项*/

{

Thread * mix_th, *pri_th; /＊混合相、第一相*/
real m_dot_v;

mix_th = THREAD_SUPER_THREAD(sec_th);
pri_th = THREAD_SUB_THREAD(mix_th, 0);

if(C_T(cell, mix_th)>=T_SAT){
    m_dot_v = 0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
                fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
    dS = 0.;
                               }
else {
    m_dot_v = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
                fabs(T_SAT-C_T(cell,mix_th))/T_SAT;

    dS = -0.1*C_R(cell, sec_th)*
      fabs(C_T(cell, sec_th) - T_SAT)/T_SAT;

       }
   
return m_dot_v;

}

DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn)
{
Thread *pri_th, *sec_th;
real m_dot;
pri_th = THREAD_SUB_THREAD(mix_th, 0);
sec_th = THREAD_SUB_THREAD(mix_th, 1);

if(C_T(cell, mix_th)>=T_SAT){
    m_dot = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
                fabs(C_T(cell, pri_th) - T_SAT)/T_SAT;

    dS = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)/T_SAT;
                               }
   
else {
    m_dot = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
                fabs(T_SAT-C_T(cell,mix_th))/T_SAT;

    dS = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)/T_SAT;}

return LAT_HT*m_dot; /＊能量大小*/
}

chbo76 发表于 2007-6-21 17:36

好是好，可有谁能把上面程序中所实现的公式给介绍一下啊？

crazyfish 发表于 2007-6-23 15:36

学习中，期待高手指点

hyacinth 发表于 2007-6-24 07:15

FLUENT 6.3 Tutorial Guide=〉Using the VOF Model

crazyfish 发表于 2007-6-24 09:00

fluent6.2有这个吗？

cora 发表于 2007-6-24 21:10

原帖由 crazyfish 于 2007-6-24 09:00 发表 http://www.chinavib.com/forum/images/common/back.gif
fluent6.2有这个吗？

有的，http://202.41.85.84/doc/fluent6.2/help/html/tg/main_pre.htm

16 Using the VOF Model

phil 发表于 2007-6-25 13:21

原帖由 hyacinth 于 2007-6-24 07:15 发表 http://www.chinavib.com/forum/images/common/back.gif
FLUENT 6.3 Tutorial Guide=〉Using the VOF Model
这个我看了，它模拟的是墨水掉瓶子里，没涉及到相变什么的啊....

chbo76 发表于 2007-6-25 18:02

模拟蒸发和冷凝要涉及汽化压强和潜热的问题。
比较复杂。

crazyfish 发表于 2007-6-26 09:51

那就是说需要编一个复杂的UDF了？头疼:@L

phil 发表于 2007-6-26 15:38

就是啊，我就不知该怎么把饱和蒸汽压还有潜热添进去...
斑竹帮帮忙啊...:'( :'( :'( :'(

crazyfish 发表于 2007-7-1 10:35

期待大侠指点中:@(

cora 发表于 2007-7-4 11:01

原帖由 crazyfish 于 2007-6-26 09:51 发表 http://www.chinavib.com/forum/images/common/back.gif
那就是说需要编一个复杂的UDF了？头疼:@L

是这样mass transfer和heat transfer都需要考虑

zjj7070702004 发表于 2008-1-17 19:31

condensation on the cold wal

jimee_lzj 发表于 2008-1-17 22:50

想问一下14楼的，你计算的初场如何给定的？！
页: [1] 2 3 4
查看完整版本: VOF怎么模拟蒸发和冷凝啊？