| 查看: 1283 | 回复: 10 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
duwei789551新虫 (小有名气)
|
[求助]
求助，两个uds变量相除就显示发散是什么问题？ 已有1人参与
|
|
|
如题，需要计算保存为两个标量方程源项的比值。当只用这一段代码时 r=C_UDSI(c,c_thread,1)/C_UDSI(c,c_thread,0); 计算一开始就显示uds-1发散； 为避免发散出现，对分母做出限制 if (C_UDSI(c,c_thread,0)<=1.0) { r=0.00000001; } else { r=C_UDSI(c,c_thread,1)/C_UDSI(c,c_thread,0); } 这样的话算到途中，uds迭代曲线将要出现的时候显示uds-1发散； 不加r=C_UDSI(c,c_thread,1)/C_UDSI(c,c_thread,0);这句话就没有问题，不知是什么原因 谢谢解答！ |
回复此楼» 猜你喜欢
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有237人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
duwei789551
新虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 709.4
- 帖子: 65
- 在线: 39.5小时
- 虫号: 3854710
- 注册: 2015-05-07
- 性别: GG
- 专业: 多相流热物理学
6楼2016-10-13 14:42:26
qc12345
至尊木虫 (职业作家)
- 应助: 146 (高中生)
- 金币: 13721.3
- 红花: 26
- 帖子: 3303
- 在线: 706小时
- 虫号: 3071487
- 注册: 2014-03-19
- 专业: 园艺学与植物营养学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
duwei789551: 金币+100, ★★★★★最佳答案, 解决问题 2016-10-13 20:05:11
感谢参与，应助指数 +1
duwei789551: 金币+100, ★★★★★最佳答案, 解决问题 2016-10-13 20:05:11
| 設s=C_UDSI(c,c_thread,0)再相除 |
» 本帖已获得的红花（最新10朵）
2楼2016-10-12 17:59:34
duwei789551
新虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 709.4
- 帖子: 65
- 在线: 39.5小时
- 虫号: 3854710
- 注册: 2015-05-07
- 性别: GG
- 专业: 多相流热物理学
送红花一朵 |
我的程序是这样的 标量方程1源相： real s； ... s=C_UDSI(c,c_thread,0； if (C_UDSI(c,c_thread,0)<=1.0) { r=0.00000001; } else { r=C_UDSI(c,c_thread,1)/s; } J=pow(rou,2.0)/roul； C_UDSI(c,c_thread,0)=J; C_UDMI(c,c_thread,0)=J; dS[eqn]=0.0; return J; 源相方程2 real s； ... s=C_UDSI(c,c_thread,0； if (C_UDSI(c,c_thread,0)<=1.0) { r=0.00000001; } else { r=C_UDSI(c,c_thread,1)/s; } ... yihao=C_UDMI(c,c_thread,0)*r+C_UDSI(c,c_thread,0)*rou； C_UDSI(c,c_thread,1)=yihao; C_UDMI(c,c_thread,5)=yihao; dS[eqn]=0.0; return yihao; 改写成这样还是算到迭代曲线出现时发散，不知道有什么问题 |
3楼2016-10-13 11:28:41
qc12345
至尊木虫 (职业作家)
- 应助: 146 (高中生)
- 金币: 13721.3
- 红花: 26
- 帖子: 3303
- 在线: 706小时
- 虫号: 3071487
- 注册: 2014-03-19
- 专业: 园艺学与植物营养学
4楼2016-10-13 14:00:53
