FLUENT中水－水蒸汽蒸发过程UDF程序段及其解释 分享 (转)

glise

1. #include "udf.h"  ‘包括常规宏
   
2. 
   
3. #include "sg_mphase.h"  ‘包括体积分数宏CVOF(C,T)
   
4. 
   
5. #define T_SAT 373    ‘定义蒸发温度100℃
   
6. 
   
7. #define LAT_HT 1.e3   ‘定义蒸发潜热J/Kg
   
8. 
   
9. DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn) ‘液相质量源项UDF
   
10. 
    
11. {
    
12. 
    
13. Thread *mix_th, *sec_th;  ‘定义计算区线指针
    
14. 
    
15.   real m_dot_l; ‘定义液相质量转移 kg/(m2.s)
    
16. 
    
17.   mix_th = THREAD_SUPER_THREAD(pri_th);  ‘指向混合区的主相即液相的指针
    
18. 
    
19.   sec_th = THREAD_SUB_THREAD(mix_th, 1);指向单相控制区的气相的指针，气相为第二相
    
20. 
    
21.   if(C_T(cell, mix_th)>=T_SAT)  ‘如果液相单元的温度高于蒸发温度，液相向气相的质量质量转移
    
22. 
    
23. {
    
24. 
    
25.   m_dot_l = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*   
    
26. 
    
27.   fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
    
28. 
    
29.     dS[eqn] = -0.1*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
    
30. 
    
31.               ‘ 定义源项对质量转移偏导 ，
    
32. 
    
33. 
    
34. 
    
35.                                }
    
36. 
    
37.   else {
    
38. 
    
39.     m_dot_l = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
    
40. 
    
41.                 fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
    
42. 
    
43. ‘如果指向混合区液相的单元温度小于蒸发温度，气相向液相的质量转移，液相得
    
44. 
    
45. 
    
46. 
    
47.     dS[eqn] = 0.;  ‘由于是气相向液相转移，所以液相的质量源项对质量转移的偏导为零
    
48. 
    
49. 
    
50. 
    
51.        }
    
52. 
    
53.    
    
54. 
    
55.   return m_dot_l;
    
56. 
    
57. }
    
58. 
    
59. 
    
60. 
    
61. DEFINE_SOURCE(vap_src, cell, sec_th, dS, eqn) ‘气相质量源项UDF
    
62. 
    
63. {
    
64. 
    
65.   Thread * mix_th, *pri_th;
    
66. 
    
67.   real m_dot_v;
    
68. 
    
69.   mix_th = THREAD_SUPER_THREAD(sec_th); ‘指向混合区的第二相即气相的指针
    
70. 
    
71.   pri_th = THREAD_SUB_THREAD(mix_th, 0); 指向单相控制区的液相的指针，液相为主相
    
72. 
    
73.   if(C_T(cell, mix_th)>=T_SAT) ‘如果混合区单元的温度高于蒸发温度，液相向气相的质量质量转移
    
74. 
    
75. {
    
76. 
    
77.     m_dot_v = 0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
    
78. 
    
79.                 fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
    
80. 
    
81.     dS[eqn] = 0.; ‘由于是液相向气相转移，所以气相的质量源项对来自液相的质量转移的偏导为零
    
82. 
    
83.                               }
    
84. 
    
85.   else {
    
86. 
    
87.     m_dot_v = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
    
88. 
    
89.                 fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
    
90. 
    
91. ‘如果指向混合区的单元温度小于蒸发温度，气相向液相的质量转移，气相失
    
92. 
    
93. 
    
94. 
    
95.     dS[eqn] = -0.1*C_R(cell, sec_th)* fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
    
96. 
    
97. 由于是气相向液相转移，所以气相的质量源项对自身的质量转移的偏导不为零
    
98. 
    
99. 
    
100. 
     
101.        }
     
102. 
     
103.   return m_dot_v;
     
104. 
     
105. }
     
106. 
     
107. DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn) ‘混合模型能量源项UDF
     
108. 
     
109. {
     
110. 
     
111.   Thread *pri_th, *sec_th;
     
112. 
     
113.   real m_dot;
     
114. 
     
115.   pri_th = THREAD_SUB_THREAD(mix_th, 0);‘指向混合区的液相的指针
     
116. 
     
117.   sec_th = THREAD_SUB_THREAD(mix_th, 1); ‘指向混合区的气相的指针
     
118. 
     
119.   if(C_T(cell, mix_th)>=T_SAT) ‘如果混合区的单元温度高于蒸发温度。质量转移由液相向气相转移，吸热，质量转移量前有负号
     
120. 
     
121. {
     
122. 
     
123.     m_dot = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
     
124. 
     
125.                 fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
     
126. 
     
127. 
     
128. 
     
129.     dS[eqn] =  -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)/T_SAT;
     
130. 
     
131.                                }
     
132. 
     
133.   else {                                               相反，气相向液相转移则放热
     
134. 
     
135.     m_dot = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
     
136. 
     
137.                 fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
     
138. 
     
139.     dS[eqn] = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)/T_SAT;}
     
140. 
     
141.   return LAT_HT*m_dot;  气化潜热与质量转移率相乘得能量源项W/m3
     
142. 
     
143. }

复制代码

chbo76

好，很详细。
关于里面涉及的公式也给出就更好了。:victory:

杨杨和脸脸

杨杨请问高手，0.1是怎么来的，单位是？源相的单位是否有限制