Codice16/2010
TitoloPositivity and conservation properties of some integration schemes for mass action kinetics
Data2010-06-04
Autore/iFormaggia, Luca; Scotti, Anna
LinkDownload full text
AbstractThe numerical schemes approximating chemical reactions according to the mass action law should reproduce at least two properties of the corresponding physical system: mass conservation and nonnegativity of the concentrations. This paper analyzes the equations of mass action kinetics providing a proof of the existence, uniqueness, and positivity of the solution under mild hypothesis on the reaction rate and the stoichiometric coefficients. We then consider some classic integration schemes in terms of conservation, positivity and accuracy compared to schemes tailored for production-destruction systems, and propose an original scheme which guarantees conservation, nonnegativity of the solution and has order of convergence between two and three.