In this paper, we study the numerical approximation of a system of partial dif-ferential equations describing the corrosion of an iron based alloy in a nuclear waste repository. In particular, we are interested in the convergence of a numerical schem
e consisting in an implicit Euler scheme in time and a Scharfetter-Gummel finite volume scheme in space.
In this paper, we consider a drift-diffusion system describing the corrosion of an iron based alloy in a nuclear waste repository. In comparison with the classical drift-diffusion system arising in the modeling of semiconductor devices, the originali
ty of the corrosion model lies in the boundary conditions which are of Robin type and induce an additional coupling between the equations. We prove the existence of a weak solution by passing to the limit on a sequence of approximate solutions given by a semi-discretization in time.
