Stiffness in chemical reaction systems is a frequently encountered computational problem, arising when different reactions in the system take place at different time-scales. Computational savings can be obtained under time-scale separation. Assuming that the system can be partitioned into slow- and fast- equilibrating subsystems, it is then possible to efficiently simulate the slow subsystem only, provided that the corresponding kinetic laws have been modified so that they reflect their dependency on the fast system. We show that the rate expectation with respect to the fast subsystem's steady-state is a continuous function of the state of the slow system. We exploit this result to construct an analytic representation of the modified rate functions via statistical modelling, which can be used to simulate the slow system in isolation. The computational savings of our approach are demonstrated in a number of non-trivial examples of stiff systems.
Efficient stochastic simulation of systems with multiple time scales via statistical abstraction
Bortolussi L;
2015
Abstract
Stiffness in chemical reaction systems is a frequently encountered computational problem, arising when different reactions in the system take place at different time-scales. Computational savings can be obtained under time-scale separation. Assuming that the system can be partitioned into slow- and fast- equilibrating subsystems, it is then possible to efficiently simulate the slow subsystem only, provided that the corresponding kinetic laws have been modified so that they reflect their dependency on the fast system. We show that the rate expectation with respect to the fast subsystem's steady-state is a continuous function of the state of the slow system. We exploit this result to construct an analytic representation of the modified rate functions via statistical modelling, which can be used to simulate the slow system in isolation. The computational savings of our approach are demonstrated in a number of non-trivial examples of stiff systems.File | Dimensione | Formato | |
---|---|---|---|
prod_424141-doc_151154.pdf
non disponibili
Descrizione: Efficient stochastic simulation of systems with multiple time scales via statistical abstraction
Tipologia:
Versione Editoriale (PDF)
Dimensione
237.75 kB
Formato
Adobe PDF
|
237.75 kB | Adobe PDF | Visualizza/Apri Richiedi una copia |
prod_424141-doc_157709.pdf
accesso aperto
Descrizione: postprint
Tipologia:
Versione Editoriale (PDF)
Dimensione
326.57 kB
Formato
Adobe PDF
|
326.57 kB | Adobe PDF | Visualizza/Apri |
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.