| 查看: 878 | 回复: 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; } } } |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
国社科又开始会评了，不知道这次命运如何
已经有13人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有13人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
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
30