| 查看: 861 | 回复: 2 | ||
[求助]
求大神帮忙看看这个动量源项为什么加载不上？哪错了，
|
|
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; } } } |
回复此楼» 猜你喜欢
应该是93bebmhtak前后十一个字符比较关键
已经有36人回复
-
关于Filecode分析方法
已经有9人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有10人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
FileCode能看出啥？
已经有21人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
1061692107
新虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 103.3
- 红花: 1
- 帖子: 102
- 在线: 27.6小时
- 虫号: 1697095
- 注册: 2012-03-17
- 性别: GG
- 专业: 煤炭地下开采
2楼2017-03-09 16:15:51
想你睡不着
新虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1209.6
- 散金: 40
- 红花: 3
- 帖子: 871
- 在线: 277.8小时
- 虫号: 1814174
- 注册: 2012-05-13
- 专业: 海洋工程
3楼2017-03-12 10:14:58
