| 查看: 793 | 回复: 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; } } } |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有8人回复
-
最失望的一年
已经有17人回复
-
为什么nbs上溴 没有产物点出现呢
已经有6人回复
-
求推荐博导
已经有4人回复
-
存款400万可以在学校里躺平吗
已经有34人回复
-
求助一下有机合成大神
已经有4人回复
-
求推荐英文EI期刊
已经有5人回复
-
26申博
已经有3人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
疑惑？
已经有5人回复
想你睡不着
新虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1209.6
- 散金: 40
- 红花: 3
- 帖子: 871
- 在线: 277.8小时
- 虫号: 1814174
- 注册: 2012-05-13
- 专业: 海洋工程
3楼2017-03-12 10:14:58