|
我编了这样一段程序。解释的时候出现:line 15: invalid type for binary expression: float * array 2 of float.的错误。line 15
以下是程序:
#include "udf.h"
#include "sg_mphase.h"
#define T_SAT 319
#define LAT_HT 2.392e6
DEFINE_SOURCE(vap_src, cell, pri_th, dS, eqn)
{
Thread * mix_th, *sec_th;
real m_dot_v,q;
mix_th = THREAD_SUPER_THREAD(pri_th);
sec_th = THREAD_SUB_THREAD(mix_th, 1);
if(C_T(cell, mix_th)<=T_SAT)
{
q=-2*C_K_L(cell,pri_th)*C_T_G(cell,pri_th)*C_VOF(cell,pri_th); /*此行提示错误*/
m_dot_v=q/LAT_HT;
dS[eqn] = -2*C_K_L(cell,pri_th)*C_T_G(cell,pri_th)/LAT_HT;
}
else
{
q=C_VOF_G(cell,sec_th);
m_dot_v=q/LAT_HT;
dS[eqn] = 0;
}
return m_dot_v;
}
DEFINE_SOURCE(liq_src, cell, sec_th, dS, eqn)
{
Thread *mix_th, *pri_th;
real m_dot_l, q;
mix_th = THREAD_SUPER_THREAD(sec_th);
pri_th = THREAD_SUB_THREAD(mix_th,0);
if(C_T(cell, mix_th)<=T_SAT)
{
q=2*C_K_L(cell,sec_th)*C_T_G(cell,sec_th)*C_VOF_G(cell,sec_th);
m_dot_l=q/LAT_HT;;
dS[eqn] =2*C_K_L(cell,pri_th)*C_T_G(cell,pri_th)/LAT_HT;
}
else
{
q=-2*C_K_L(cell,pri_th)*C_T_G(cell,pri_th)*C_VOF_G(cell,pri_th);
m_dot_l=q/LAT_HT;
dS[eqn]=0.;
}
return m_dot_l;
}
DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn)
{
Thread *pri_th, *sec_th;
real m_dot, q;
pri_th = THREAD_SUB_THREAD(mix_th, 0);
sec_th = THREAD_SUB_THREAD(mix_th, 1);
if(C_T(cell, mix_th)<=T_SAT)
{
q=2*C_K_L(cell,pri_th)*C_T_G(cell,pri_th)*C_VOF_G(cell,pri_th);;
dS[eqn] = 0.;
}
else
{
q=-2*C_K_L(cell,sec_th)*C_T_G(cell,sec_th)*C_VOF_G(cell,sec_th);
dS[eqn]=0.;
}
return q;
}
我的udf编的不咋样。你帮我看看。我定义的是气液相变的源项。多谢啦 |
|