Until 1991, the simulation of coupled systems, in which electrochemical reactions are coupled with homogeneous chemical reactions in such a way that at least some of the relevant partial differential transport equations (pdes) contain terms for more than one species, was considered hard. The reason is that many such reaction mechanisms lead to thin reaction layers [1], necessitating the use of ...