Yes, it is possible to do so. See the work below. One can use the finite element method to describe the concentrations of the various chemical agents as a function of position and posit some reaction kinetics at the surface and in solution (e.g. metal dissolution may be first order with respect to the concentration of H+). This can all be done with the flexible multiphysics in COMSOL.