|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
x
改了个蒸发模型,代码是根据那个用烂了的例子改的,但是就是显示出现错误,网上用的办法我也试过没有用,甚至我把TRANSFER后面的内容都加//了调试,结果还是这样的:
1>c:\users\administrator\desktop\1\1\eva_rate.c(3): error C2143: 语法错误 : 缺少“)”(在“常量”的前面)
1>c:\users\administrator\desktop\1\1\eva_rate.c(3): error C2143: 语法错误 : 缺少“{”(在“常量”的前面)
1>c:\users\administrator\desktop\1\1\eva_rate.c(3): error C2059: 语法错误:“<未知>”
1>c:\users\administrator\desktop\1\1\eva_rate.c(3): error C2059: 语法错误:“)”
1>c:\users\administrator\desktop\1\1\eva_rate.c(24): warning C4003: “C_T”宏的实参不足
求问还有没有别的办法改这个,我以为是电脑问题换了一台也是这样。
下面是代码
#include "udf.h" //包括常规宏
DEFINE_MASS_TRANSFER(mlg_r,c,thread,2,-1,3,-1)
{
real m_lg; //蒸发速率
real k1=2.5226e-4; //蒸发速率系数
real k2=8.3117e-8; //指数系数
real a=-0.6044; //初始温差幂指数
real b=3.0260; //平衡温度幂指数
real v_liq=1.6965e-4; //液相体积
real flow_time; //蒸发时间
real T0=343.15; //起始温度70
Thread *liq=THREAD_SUB_THREAD(thread,2); //液相指针
Thread *gas=THREAD_SUB_THREAD(thread,3); //气相指针
m_lg=0;
real flow_time=FLOW_TIME;
//T_sat=312.15; //5000Pa仅限
real T_sat=1657.459/(7.074056-log10(C_P(cell,liq)))+46.13; //由安托尼方程得平衡温度
if (C_T(cell,liq)>=T_sat)
{
m_lg=k1*k2*pow((T0-T_sat),a+1)*pow(T_sat,b)*exp(-k2*pow((T0-T_sat),a)*pow(T_sat,b)* flow_time)*C_VOF(cell,liq)/v_liq;
}
if ((m_lg==0)&&(C_T(cell,gas)<T_sat))
{
m_lg=-0.1*C_VOF(cell,gas)*C_R(cell,gas)*(T_sat-C_T(cell,gas))/T_sat;
}
return (m_lg);
}
|
|