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

phil

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

linda

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

1. #include "udf.h"
   
2. #include "sg_mphase.h"
   
3. #define T_SAT 373 /＊定义饱和温度*/
   
4. #define LAT_HT 1.e3 /＊定义水蒸汽潜热*/
   
5. 
   
6. DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn) /＊定义液相源项*/
   
7. {
   
8. 
   
9.   Thread *mix_th, *sec_th; /＊混合相、第二相*/
   
10.   real m_dot_l; /＊蒸发冷凝速率*/
    
11. 
    
12.   mix_th = THREAD_SUPER_THREAD(pri_th);
    
13.   sec_th = THREAD_SUB_THREAD(mix_th, 1);
    
14. 
    
15.   if(C_T(cell, mix_th)>=T_SAT){
    
16.     m_dot_l = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
    
17.                 fabs(C_T(cell, pri_th) - T_SAT)/T_SAT;
    
18. 
    
19.    dS[eqn] = -0.1*C_R(cell, pri_th)*
    
20.       fabs(C_T(cell, pri_th) - T_SAT)/T_SAT;
    
21.                                }
    
22.   else {
    
23.     m_dot_l = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
    
24.                 fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
    
25. 
    
26.     dS[eqn] = 0.;
    
27. 
    
28.        }
    
29.    
    
30.   return m_dot_l;
    
31. 
    
32. }
    
33. 
    
34. DEFINE_SOURCE(vap_src, cell, sec_th, dS, eqn) /＊定义气相源项*/
    
35. 
    
36. {
    
37. 
    
38.   Thread * mix_th, *pri_th; /＊混合相、第一相*/
    
39.   real m_dot_v;
    
40. 
    
41.   mix_th = THREAD_SUPER_THREAD(sec_th);
    
42.   pri_th = THREAD_SUB_THREAD(mix_th, 0);
    
43. 
    
44.   if(C_T(cell, mix_th)>=T_SAT){
    
45.     m_dot_v = 0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
    
46.                 fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
    
47.     dS[eqn] = 0.;
    
48.                                }
    
49.   else {
    
50.     m_dot_v = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
    
51.                 fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
    
52. 
    
53.     dS[eqn] = -0.1*C_R(cell, sec_th)*
    
54.       fabs(C_T(cell, sec_th) - T_SAT)/T_SAT;
    
55. 
    
56.        }
    
57.    
    
58.   return m_dot_v;
    
59. 
    
60. }
    
61. 
    
62. DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn)
    
63. {
    
64.   Thread *pri_th, *sec_th;
    
65.   real m_dot;
    
66.   pri_th = THREAD_SUB_THREAD(mix_th, 0);
    
67.   sec_th = THREAD_SUB_THREAD(mix_th, 1);
    
68. 
    
69.   if(C_T(cell, mix_th)>=T_SAT){
    
70.     m_dot = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
    
71.                 fabs(C_T(cell, pri_th) - T_SAT)/T_SAT;
    
72. 
    
73.     dS[eqn] = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)/T_SAT;
    
74.                                }
    
75.    
    
76.   else {
    
77.     m_dot = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
    
78.                 fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
    
79. 
    
80.     dS[eqn] = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)/T_SAT;}
    
81. 
    
82.   return LAT_HT*m_dot; /＊能量大小*/
    
83. }

复制代码

chbo76

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

crazyfish

学习中，期待高手指点

hyacinth

FLUENT 6.3 Tutorial Guide=〉Using the VOF Model

crazyfish

fluent6.2有这个吗？

cora

原帖由 crazyfish 于 2007-6-24 09:00 发表

                               
登录/注册后可看大图

fluent6.2有这个吗？

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

16 Using the VOF Model

phil

原帖由 hyacinth 于 2007-6-24 07:15 发表

                               
登录/注册后可看大图

FLUENT 6.3 Tutorial Guide=〉Using the VOF Model

这个我看了，它模拟的是墨水掉瓶子里，没涉及到相变什么的啊....

chbo76

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

crazyfish

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

phil

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

crazyfish

期待大侠指点中:@(

cora

原帖由 crazyfish 于 2007-6-26 09:51 发表

                               
登录/注册后可看大图

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

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

zjj7070702004

jimee_lzj

想问一下14楼的，你计算的初场如何给定的？！