Adapted numerical modelling of the Belousov–Zhabotinsky reaction

Adapted numerical schemes for the integration of differential equations generating periodic wavefronts have reported benefits in terms of accuracy and stability. This work is focused on differential equations modelling chemical phenomena which are characterized by an oscillatory dynamics. The adaptation is carried out through the exponential fitting technique, which is specially suitable to follow the apriori known qualitative behavior of the solution. In particular, we have merged this strategy with the information coming from existing theoretical studies and especially the observation of time series. Numerical tests will be provided to show the effectiveness of this problem-oriented approach.