24小时热门版块排行榜    

CyRhmU.jpeg
查看: 544  |  回复: 0

hangs999111

新虫 (小有名气)

[求助] UDF提示未声明和“=”未能把“double”转变为“real *”

UDF提示未声明和“=”未能把“double”转变为“real *”
#include "udf.h"
#include "sg_mphase.h"
DEFINE_SOURCE(liq_m,c,l_t,dS,eqn)    /*液相质量输运源项*/
{
real m_dot_l;                     /*液相质量输运量*/
real stc;                           /*液氧表面张力系数*/
real fv;                            /*气相质量分数*/
real pres_vap,pres_sat,pres_turb;   /*液氧饱和蒸汽压(考虑湍流影响)*/


Thread *m_t,*v_t;                   /*计算区域指针*/

  m_t=THREAD_SUPER_THREAD(l_t);      /*混合区液相指针*/
  v_t=THREAD_SUB_THREAD(m_t,1);      /*单相区气相指针*/

  fv=C_VOF(c,v_t)*C_R(c,v_t)/(C_VOF(c,v_t)*C_R(c,v_t)+C_VOF(c,l_t)*C_R(c,l_t));

  if(C_T(c,m_t)>2503)/*饱和温度*/
    {
     stc=0;
     pres_sat=2.561e7;
    }
  else if(C_T(c,m_t)>1156)
    {
     stc=109.7;
    /* pres_sat=133.322*pow(10,20.6695-526.97/C_T(c,l_t)-6.7062*log10(C_T(c,l_t))+0.012926*C_T(c,l_t)-C_T(c,l_t)*C_T(c,l_t)*(9.8832e-13));*/
     
     pres_sat=1000000*pow(2.71828,11.9463-12633/C_T(c,l_t)-0.4672*log(C_T(c,l_t)));
  
  }
  else
    {
     stc=109.7;
     pres_sat=23546.4;
    }
  pres_turb=0.39*C_R(c,l_t)*C_K(c,m_t);
  pres_vap=pres_sat+0.5*pres_turb;
if(C_P(c,m_t)<=pres_vap)             /*液相汽化*/
   {
    m_dot_l=-0.02*sqrt(C_K(c,m_t))/stc*C_R(c,l_t)*C_R(c,v_t)*sqrt((2*(pres_vap-C_P(c,m_t)))/(3*C_R(c,l_t)))*(1-fv);
    dS=-0.02*sqrt(C_K(c,m_t))/stc*C_R(c,l_t)*C_R(c,v_t)*sqrt((2*(pres_vap-C_P(c,m_t)))/(3*C_R(c,l_t)));
   }
else                                 /*气相液化*/
   {
    m_dot_l=0.01*sqrt(C_K(c,m_t))/stc*C_R(c,l_t)*C_R(c,l_t)*sqrt((2*(C_P(c,m_t)-pres_vap))/(3*C_R(c,l_t)))*fv;
    dS=0;
   }
return m_dot_l;
UDF提示未声明和“=”未能把“double”转变为“real *”

这个该如何修改呢?
没金币了,抱歉啊,还望大家帮帮忙
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hangs999111 的主题更新
信息提示
请填处理意见