Numerical analysis of reversible A + B <-> C reaction-diffusion systems

We develop an effective numerical method of studying large-time properties of reversible reaction-diffusion systems of type A + B <-> C with initially separated reactants. Using it we find that there are three types of asymptotic reaction zones. In particular we show that the reaction rate can be locally negative and concentrations of species A and B can be nonmonotonic functions of the space coordinate x, locally significantly exceeding their initial values.