ode15s slows down at reaching steady state SOLVED
조회 수: 2 (최근 30일)
이전 댓글 표시
I'm solving a system of rather stiff equation describing a microkinetic model of methanol formation (30 species, ~70 reactions and vastly different rate constants).
I'm achieving satisfactory speeds of the ode15s solver during simulation since i started using an analitycal jacobian (tremendous increase in speed and stability of the solver). However, as the solution approaches steady state, the solver reduces the speed greatly and gets stuck. When using just the right time interval to get to reduction part, the ratio of sucessful and failed attempts reported increases drastically.
I'm not sure what's causing this behaviour. Is anyone else experiencing similar problems? When using numerical estimation of the jacobian, the solver always speeded up when steady state was reached. Has this something to do with numerical fluctuations in steady state due to floating point number rounding?
Any insight or help would be much appreciated.
EDIT: The problem was fixed by reducing the absolute tolerance of the solver.
댓글 수: 3
Torsten
2017년 1월 12일
I assume you checked the Jacobian
https://de.mathworks.com/help/optim/ug/checking-validity-of-gradients-or-jacobians.html
?
Best wishes
Torsten.
답변 (0개)
참고 항목
카테고리
Help Center 및 File Exchange에서 Ordinary Differential Equations에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!