| 查看: 613 | 回复: 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;  这个该如何修改呢？  没金币了，抱歉啊，还望大家帮帮忙 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有11人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有11人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有10人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有12人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有10人回复
