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

jiang.yanyuan

请高手帮忙说说冷凝的问题，怎么用UDF模型，我要模拟一个冷凝器，里面是水蒸气冷凝，外面用水冷却。

wanghack

很好，支持一下！！

lw_zhao

请问你们有冷凝的UDF程序吗？

lw_zhao

请问你有冷凝的UDF吗？传一下，谢谢！

lw_zhao

怎么没人有吗？帮帮忙嘛！！！！

lw_zhao

大侠，你用的冷凝的UDF能贴出来看看吗？

lw_zhao

请问热交换模型能模拟板式换热器的冷凝吗?

wajish_zj

好东西，支持

fzhqhd

还是一头雾水

a04031424

冷凝的UDF哪位大哥有阿?目前急需阿....始终解决不了这个问题啊.

Heidiadalheid

一份网上流传比较多的计算冷凝蒸发的udf

1. #include "udf.h"  ‘包括常规宏
   
2. #include "sg_mphase.h"  ‘包括体积分数宏CVOF(C,T)
   
3. #define T_SAT 373    ‘定义蒸发温度100℃
   
4. #define LAT_HT 1.e3   ‘定义蒸发潜热J/Kg
   
5. 
   
6. DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn) ‘液相质量源项UDF
   
7. {
   
8.   Thread *mix_th, *sec_th;  ‘定义计算区域线指针
   
9.   real m_dot_l; ‘定义液相质量转移 kg/(m2.s)
   
10.   mix_th = THREAD_SUPER_THREAD(pri_th);  ‘指向混合区的主相即液相的指针
    
11.   sec_th = THREAD_SUB_THREAD(mix_th, 1);指向单相控制区的气相的指针，气相为第二相
    
12.   if(C_T(cell, mix_th)>=T_SAT)  ‘如果液相单元的温度高于蒸发温度，液相向气相的质量质量转移
    
13. {
    
14.   m_dot_l = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*   
    
15.   fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
    
16.     dS[eqn] = -0.1*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
    
17.               ‘ 定义源项对质量转移偏导，
    
18. }  else {
    
19.     m_dot_l = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
    
20.                 fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
    
21. ‘如果指向混合区液相的单元温度小于蒸发温度，气相向液相的质量转移，液相得
    
22.     dS[eqn] = 0.;  ‘由于是气相向液相转移，所以液相的质量源项对质量转移的偏导为零
    
23. }
    
24. return m_dot_l;
    
25. }
    
26. 
    
27. DEFINE_SOURCE(vap_src, cell, sec_th, dS, eqn) ‘气相质量源项UDF
    
28. {
    
29.   Thread * mix_th, *pri_th;
    
30.   real m_dot_v;
    
31.   mix_th = THREAD_SUPER_THREAD(sec_th); ‘指向混合区的第二相即气相的指针
    
32.   pri_th = THREAD_SUB_THREAD(mix_th, 0); 指向单相控制区的液相的指针，液相为主相
    
33.   if(C_T(cell, mix_th)>=T_SAT) ‘如果混合区单元的温度高于蒸发温度，液相向气相的质量质量转移
    
34. {
    
35.     m_dot_v = 0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
    
36.                 fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
    
37.     dS[eqn] = 0.; ‘由于是液相向气相转移，所以气相的质量源项对来自液相的质量转移的偏导为零
    
38.    } else {
    
39.     m_dot_v = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
    
40.                 fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
    
41.                     ‘如果指向混合区的单元温度小于蒸发温度，气相向液相的质量转移
    
42.     dS[eqn] = -0.1*C_R(cell, sec_th)* fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
    
43.                     由于是气相向液相转移，所以气相的质量源项对自身的质量转移的偏导不为零
    
44. }
    
45.   return m_dot_v;
    
46. }
    
47. 
    
48. DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn) ‘混合模型能量源项UDF
    
49. {
    
50.   Thread *pri_th, *sec_th;
    
51.   real m_dot;
    
52.   pri_th = THREAD_SUB_THREAD(mix_th, 0);‘指向混合区的液相的指针
    
53.   sec_th = THREAD_SUB_THREAD(mix_th, 1); ‘指向混合区的气相的指针
    
54.   if(C_T(cell, mix_th)>=T_SAT) ‘如果混合区的单元温度高于蒸发温度。质量转移由液相向气相转移，吸热，质量转移量前有负号
    
55. {
    
56.     m_dot = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
    
57.                 fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
    
58.     dS[eqn] =  -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)/T_SAT;
    
59.    }  else {                              相反，气相向液相转移则放热
    
60.     m_dot = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
    
61.                 fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
    
62.     dS[eqn] = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)/T_SAT;}
    
63.   return LAT_HT*m_dot;  气化潜热与质量转移率相乘得能量源项W/m3
    
64. }

复制代码

Heidiadalheid

1. #define T_con 373.15 /*液化温度*/
   
2. 
   
3. /*计水蒸气在液化温度下变成水*/
   
4. DEFINE_MASS_TRANSFER(gas_liq_source,cell,thread,from_index,from_species_index, to_index, to_species_index)
   
5. {
   
6. real m_lg;
   
7. Thread *gas = THREAD_SUB_THREAD(thread, from_index);
   
8. Thread *liq = THREAD_SUB_THREAD(thread, to_index);
   
9. m_lg = 0.;
   
10. if (C_T(cell, liq) >= T_con)
    
11. {
    
12. m_lg = -0.1*C_VOF(cell,liq)*C_R(cell,liq)*fabs(C_T(cell,liq)-T_con)/T_con;
    
13. }
    
14. if ((m_lg == 0. ) && (C_T(cell, gas) <= T_con))
    
15. {
    
16. m_lg = 0.1*C_VOF(cell,gas)*C_R(cell,gas)*fabs(T_con-C_T(cell,gas))/T_con;
    
17. }
    
18. return (m_lg);
    
19. }

复制代码

Heidiadalheid

1. #define T_frost 273.15 /*凝固温度*/
   
2. 
   
3. DEFINE_MASS_TRANSFER(liq_ice_source,cell,thread,from_index,from_species_index, to_index, to_species_index)
   
4. {
   
5. real m_lg;
   
6. Thread *liq = THREAD_SUB_THREAD(thread, from_index);
   
7. Thread *soild = THREAD_SUB_THREAD(thread, to_index);
   
8. m_lg = 0.;
   
9. if (C_T(cell, soild) >= T_frost)
   
10. {
    
11. m_lg = -0.1*C_VOF(cell,soild)*C_R(cell,soild)*fabs(C_T(cell,soild)-T_frost)/T_frost;
    
12. }
    
13. if ((m_lg == 0. ) && (C_T(cell, liq) <= T_frost))
    
14. {
    
15. m_lg = 0.1*C_VOF(cell,liq)*C_R(cell,liq)*fabs(T_frost-C_T(cell,liq))/T_frost;
    
16. }
    
17. return (m_lg);
    
18. }

复制代码

zjh432105

回复 34 # lw_zhao 的帖子

你也做冷凝吗，可以交流下吗565489726

zjh432105

回复 42 # Heidiadalheid 的帖子

你好，这个程序只有水的质量变化是吧