Numerical modeling of bronze coating formation using concentration dependent diffusivity