12/3/2019
Ali M. Sahlodin
Department of Chemical Engineering
AmirKabir University of Technology
1397 S.H
Methods of event detection and location
Mode transition
Copyright © Ali M. Sahlodin, Dept. of Chemical Engineering, AmirKabir Univ. of Tech. 2
1
12/3/2019
Change in differential index
Example
DAE solver fails at the point of event
ISTATE = -7
MSG = the corrector could not converge (t = 1.38629)
Copyright © Ali M. Sahlodin, Dept. of Chemical Engineering, AmirKabir Univ. of Tech. 3
Transition function must comply with DOF in
the new mode.
Solve re-initialization system
Galan et al (1999) Applied Num. Math
Copyright © Ali M. Sahlodin, Dept. of Chemical Engineering, AmirKabir Univ. of Tech. 4
2
12/3/2019
Simple example revisited.
See your notes for the solution
Copyright © Ali M. Sahlodin, Dept. of Chemical Engineering, AmirKabir Univ. of Tech. 5
No guarantee of consistency
between mesh points
Inefficient integration
Integration failure
ACM Transactions on Modeling and Computer Simulation, Vol. 6, No. 2,
April 1996, Pages 137–165
Copyright © Ali M. Sahlodin, Dept. of Chemical Engineering, AmirKabir Univ. of Tech. 6