| 查看: 2888 | 回复: 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 |
回复此楼» 猜你喜欢
交叉科学部支持青年基金，对三无青椒是个机会吗？
已经有4人回复
-
招博士
已经有6人回复
-
限项规定
已经有8人回复
-
国家基金申请书模板内插入图片不可调整大小?
已经有5人回复
-
国家级人才课题组招收2026年入学博士
已经有5人回复
-
Fe3O4@SiO2合成
已经有6人回复
-
青年基金C终止
已经有4人回复
-
青椒八年已不青，大家都被折磨成啥样了？
已经有7人回复
-
为什么nbs上溴 没有产物点出现呢
已经有10人回复
-
救命帖
已经有11人回复
6楼2018-05-28 09:04:24
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