This article presents an approach to solving problems related to galvanic corrosion that involve moving boundaries (due to preferential corrosion of one of the metals in the system). The method incorporates an adaptive (node based, finite difference) grid technique for the treatment of boundary-related singularities that arise in the calculation of the electric potential. Simulation of the time evolution of the damage done by the corroding interface is performed using of a level set formulation. An analysis of the convergence of the method and a comparison with experimental data from the literature are included.