| 查看: 811 | 回复: 0 | ||||
[交流]
用Matlab的ODE求解器计算反应动力学，如何在反应物量趋近0时停止反应
|
|
求教群里各位大侠，用Matlab的ODE求解器计算反应动力学时，由于反应物的反应动力学数量级比较大，当反应物量趋近于0时，容易积分得到负值；由于ode求解器是对反应物量和产物量变化速率方程组进行求解，用nonnegative会导致修正后反应物消耗量与不需要修正的反应物及产物间不符合化学计量关系。 如燃烧反应： CO+0.5O2=CO2 function dy=fun(...) r=r_co(...); dco=-r; do2=-0.5*r; dco2=r; dy=[dco;do2;dco2]; end 当CO趋近于0，O2>>0时，若对CO采用非负约束，则消耗的CO与O2无法满足计量关系，也就是说结果原子会不守恒。请问大家有没有遇到过这类问题，是如何解决以保证积分不为负呢？ |
回复此楼» 收录本帖的淘帖专辑推荐
 matlab |
» 猜你喜欢
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
UJN物理学专业调剂
已经有12人回复
-
物理学I论文润色/翻译怎么收费?
已经有94人回复
-
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组招收2026硕士研究生1名
已经有0人回复
-
山西大同大学物理学硕士研究生招收调剂生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
一志愿物理所凝聚态物理352求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有1人回复
