Modeling the kinetics of a precipitation dissolution reaction occurring in a porous medium where diffusion alsotakes place leads to a system of two parabolic equations and one ordinary differential equation coupled with a stiff reaction term. This system is discretized by a finite volume scheme which is suitable for the approximation of thediscontinuous reaction term of unknown sign. Discrete solutions are shown to exist and converge towards a weak solution of the continuous problem. Uniqueness is proved under a Lipschitz condition on the equilibrium gap function. Numerical tests are shown which prove the efficiency of the scheme.

We present the numerical analysis on the Poisson problem of two mixed Petrov-Galerkin finite volume schemes for equations in divergence form $\mathop{\rm div}\nolimits\varphi(u,\nabla u)=f$ . The first scheme, which has been introduced in [CITE], is a generalization in two dimensions of Keller's box-scheme. The second scheme is the dual of the first one, and is a cell-centeredscheme for u and the flux φ. For the first scheme, the two trial finite element spaces are the nonconforming space of Crouzeix-Raviartfor the primal unknown u and the div-conformingspace of Raviart-Thomas for the flux φ. The two test spaces are the functions constant per cell both for the conservative and for the fluxequations.We prove an optimal second order error estimate for the box scheme and we emphasize the link between this scheme and the post-processing of Arnold and Brezzi of the classical mixed method.