Positivity and conservation properties of some integration schemes for mass action kinetics

Keywords

Code:
16/2010
Title:
Positivity and conservation properties of some integration schemes for mass action kinetics
Date:
Friday 4th June 2010
Author(s):
Formaggia, Luca; Scotti, Anna
Download link:
Abstract:
The 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.