冷却塔 UDF

xinxinwang

请教达人，在冷却塔模拟中，填料区里面发生的传质现象怎么表述，查阅了相关文献，有两个源项UDF需要编译，一个是质量源项，一个是能量源项，但我不会，请教达人们，提供相关的UDF参考一下~~谢谢了

xinxinwang

帮忙啊

auryyuan

我在做空化的模拟，编写的UDF如下
DEFINE_CAVITATION_RATE(user_cav_rate, c, t, p, rhov, rhol, mafv, p_v, h_gas, f_gas, m_dot)                             {
                        
                        
Thread **pt=THREAD_SUB_THREADS(t);                     
real vofv=C_VOF(c,pt[1]);                                                      
real vofl=1-vofv;                                                              
real hl=C_H(c,pt[0]);                                                           
real hv=C_H(c,pt[1]);                                                         
real p_vapor=*p_v;                                                            
real T=C_T(c,t);                                                               
real Cp=C_CP(c,t);                                                        
real P=C_P(c,t);
real rho_L= *rhol;                                                           
real rho_V= *rhov;                                                        
real p_operating, dp;                                                         

  p_operating=RP_Get_Real("operating-pressure");                                   
           
  dp=P+p_operating;                                                            
                                         
                        
   

  if(dp>p_vapor)                                                                     
                 
*m_dot = 1000*3*0.597*vofv*sqrt(2*(dp-p_vapor)/(3*958.345))*(2-(sqrt(3/PI)*(rho_L*Cp*fabs(373-T)/((hv-hl)*rho_V)*sqrt(C_UDSI(c,t,0)/0.009)))/(sqrt(2*(dp-p_vapor)/(3*958.345))));
                 

                     
          else
                  
                    
*m_dot = -5.0*3*5.0*100*0.597*vofl*sqrt(2*(p_vapor-dp)/(3*958.345))*(2+(sqrt(3/PI)*(rho_L*Cp*fabs(T-373)/((hv-hl)*rho_V)*sqrt(C_UDSI(c,t,0)/0.009)))/(sqrt(2*(p_vapor-dp)/(3*958.345))));
                 
}

该函数便已没有错误，但是计算的话就会出现错误提示
Error:wrong return type:float udf fucntion expected
如果我直接赋值
*m_dot = 0.1也还是同样的错误


另外，函数名中已经包括很多变量，例如液体密度
如果我想取出液体密度 是应该怎么取啊？

real rho_L= *rhol;
还是
real rho_L= C_R(c,t);

求指教啊