On a Variational Method for Stiff Differential Equations Arising from Chemistry Kinetics

For the approximation of stiff systems of ODEs arising from chemistry kinetics, implicit integrators emerge as good candidates. This paper proposes a variational approach for this type of systems. In addition to introducing the technique, we present its most basic properties and test its numerical performance through some experiments. The main advantage with respect to other implicit methods is that our approach has a global convergence. The other approaches need to ensure convergence of the iterative scheme used to approximate the associated nonlinear equations that appear for the implicitness. Notice that these iterative methods, for these nonlinear equations, have bounded basins of attraction