求大神看看能量方程源项udf问题出在哪里（源项由空化过程中物质汽化产生）
#include "udf.h"DEFINE_SOURCE(energy1, cell, t, dS, eqn)
{
Thread **pt;
real re=199176;
real MDv=1.138;
real MDl=806.08;
real Ps=C_P(cell, t);
real av=C_VOF(cell, pt);
real MDm=av*MDv+(1-av)*MDl;
real RB=pow(3*av/((1-av)*4*3.14159*pow(10, 8)), 1/3);
real Pv=0.0583*pow(C_T(cell,t), 4)-8.9018*pow(C_T(cell,t), 3)+517*pow(C_T(cell,t), 2)-16031*C_T(cell,t)+278654;
real source;
if(Ps<Pv)
{
source=-(3*av*(1-av)*MDv*MDl*pow(2*(Pv-Ps)/3/MDl, 1/2)/RB/MDm)*re;
dS=-(3*av*(1-av)*MDv*MDl*re/RB/MDm)*(1/(3*MDl*pow(2*(Pv-Ps)/3/MDl, 1/2)));
}
else
{
source=(3*av*(1-av)*MDv*MDl*pow(2*(Ps-Pv)/3/MDl, 1/2)/RB/MDm)*re;
dS=-(3*av*(1-av)*MDv*MDl*re/RB/MDm)*(1/(3*MDl*pow(2*(Ps-Pv)/3/MDl, 1/2)));
}
return source;
}
目前已经发现问题出在第九行（如何获得网格中气相体积/质量分数）。以上为错误的udf，通过添加能量方程源项后一步都算不了。但是把第九行气相分数改成常数可以计算，有没有大神教教小弟第九步到底怎么获得网格内的气相体积分数。能指明个大致方向就万分感谢了。
页:
[1]