24小时热门版块排行榜

返回列表

黄涛神

捐助贵宾 (小有名气)

应助: 2 (幼儿园)
金币: 182.1
红花: 1
帖子: 141
在线: 289小时
虫号: 4743146
注册: 2016-06-01
专业: 金属材料的凝固与结晶学

[求助] 求大神帮忙看看这个动量源项为什么加载不上？哪错了，

DEFINE_SOURCE(s_l_interaction_resistance_z,c,t,dS,eqn)
{
real source;
real fs,fl,rho,rhos,rhol,u,v,w,a,b,cc;
u=C_U(c,t),v=C_V(c,t),w=C_W(c,t);
rhos=7200.,rhol=6900.;
rho=(1-C_VOF(c,t,1))*rhos+C_VOF(c,t,1)*rhol;/*as to rho=C_R(c,t);*/
fs=rhos*(1-C_VOF(c,t,1))/rho;
fl=1-fs;
a=C_UDMI(c,t,5)-2.*C_VOF(c,t)*C_UDMI(c,t,5);
b=C_UDMI(c,t,3)-2.*C_VOF(c,t)*C_UDMI(c,t,3);
cc=C_UDMI(c,t,4)-2.*C_VOF(c,t)*C_UDMI(c,t,4);

if ((C_VOF(c,t,1)>0&&C_VOF(c,t,1)<1.)&&(C_VOF(c,t,0)=0))
{
if((THREAD_STORAGE(t,SV_UDS_I(10))!=NULL)&&(T_STORAGE_R_NV(t,SV_UDSI_G(10))!=NULL))
{
source=-ND_DOT(C_UDSI_G(c,t,10)[0],C_UDSI_G(c,t,10)[1],C_UDSI_G(c,t,10)[2],u,v,w)-\
C_UDSI(c,t,10)*(C_DUDX(c,t)+C_DVDY(c,t)+C_DWDZ(c,t));

dS[eqn]=-(C_DUDX(c,t)+C_DVDY(c,t)+C_DWDZ(c,t))*rho*fs*fl-2*rhos*w*rhol*(rho*a-C_VOF(c,t)*(1-C_VOF(c,t))*(rhol-rhos)*C_UDMI(c,t,5))/(rho*rho)-\
rhos*rhol*C_VOF(c,t)*(1-C_VOF(c,t))*C_W_G(c,t)[2]/rho-rhos*rhol*C_VOF(c,t)*(1-C_VOF(c,t))*C_W_G(c,t)[0]/rho-\
rhos*u*rhol*(rho*b-C_VOF(c,t)*(1-C_VOF(c,t))*(rhol-rhos)*C_UDMI(c,t,3))/(rho*rho)-rhos*rhol*C_VOF(c,t)*(1-C_VOF(c,t))*C_W_G(c,t)[1]/rho-\
rhos*v*rhol*(rho*cc-C_VOF(c,t)*(1-C_VOF(c,t))*(rhol-rhos)*C_UDMI(c,t,4))/(rho*rho);
return source;
}
}
}

回复此楼