| 查看: 2857 | 回复: 5 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
用UDF定义电池的热源，以下是电池发热表达式以及程序，结果发散了，问题在哪？
|
||
|
#include "udf.h" DEFINE_SOURCE(heat_gen,c,t,dS,eqn) { real i=15.; real c1=15.; real F=96485.; real source; real t1=RP_Get_Real("flow-time" ); real T=C_T(c,t); real soc=1-i*t1/c1; if(0<=soc<=0.77) { if(293<=T<303) source=2.258e-6*pow(soc,-0.3952)*i*i-i*T*(99.88*soc-76.67)/F; else if(303<=T<313) source=1.857e-6*pow(soc,-0.2787)*i*i-i*T*(99.88*soc-76.67)/F; else source=1.659e-6*pow(soc,-0.1692)*i*i-i*T*(99.88*soc-76.67)/F; } else if(0.77<soc<=0.87) { if(293<=T<303) source=2.258e-6*pow(soc,-0.3952)*i*i+i*T*30/F; else if(303<=T<=313) source=1.857e-6*pow(soc,-0.2787)*i*i+i*T*30/F; else source=1.659e-6*pow(soc,-0.1682)*i*i+i*T*30/F; } else { if(293<=T<303) source=2.258e-6*pow(soc,-0.3952)*i*i+i*T*20/F; else if(303<=T<313) source=1.857e-6*pow(soc,-0.2787)*i*i+i*T*20/F; else source=1.659e-6*pow(soc,-0.1692)*i*i+i*T*20/F; } dS[eqn]=0; return source; }  ~]M@IE8{WRB79@99SPKIW_T.png  5O]L2TN5DYKXJFH733W}LPD.png  )}6FJ5$8DTIZSDUCX9GP({R.png  {8([`H`A{52%JV5I~K]0Q3G.png |
回复此楼» 猜你喜欢
博士申请都是内定的吗？
已经有12人回复
-
谈谈两天一夜的“延安行”
已经有5人回复
-
氨基封端PDMS和HDI反应快速固化
已经有11人回复
-
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有11人回复
-
论文投稿求助
已经有4人回复
-
Applied Surface Science 这个期刊。有哪位虫友投过的能把word模板发给我参考一下嘛
已经有3人回复
-
投稿精细化工
已经有6人回复
MoonLee777
新虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 2274.6
- 红花: 17
- 帖子: 839
- 在线: 78小时
- 虫号: 5947368
- 注册: 2017-03-12
- 性别: GG
- 专业: 流体力学
5楼2018-05-27 11:37:15
2楼2018-05-23 17:17:59
MoonLee777
新虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 2274.6
- 红花: 17
- 帖子: 839
- 在线: 78小时
- 虫号: 5947368
- 注册: 2017-03-12
- 性别: GG
- 专业: 流体力学
3楼2018-05-26 00:36:31
4楼2018-05-27 10:27:14