UDF提示未聲明和“=”未能把“double”轉(zhuǎn)變?yōu)椤皉eal *”
#include "udf.h"
#include "sg_mphase.h"
DEFINE_SOURCE(liq_m,c,l_t,dS,eqn) /*液相質(zhì)量輸運(yùn)源項(xiàng)*/
{
real m_dot_l; /*液相質(zhì)量輸運(yùn)量*/
real stc; /*液氧表面張力系數(shù)*/
real fv; /*氣相質(zhì)量分?jǐn)?shù)*/
real pres_vap,pres_sat,pres_turb; /*液氧飽和蒸汽壓(考慮湍流影響)*/
Thread *m_t,*v_t; /*計(jì)算區(qū)域指針*/
m_t=THREAD_SUPER_THREAD(l_t); /*混合區(qū)液相指針*/
v_t=THREAD_SUB_THREAD(m_t,1); /*單相區(qū)氣相指針*/
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”轉(zhuǎn)變?yōu)椤皉eal *”]()
這個該如何修改呢?
 沒金幣了,抱歉啊,還望大家?guī)蛶兔? |