Ion exchange on glass substrates : numerical determination of concentration profiles

A numerical method of calculating ion exchanged concentration profiles is presented. Two programmes based on this method have been implemented. The first takes into account the Poisson equation (no simplifying assumption) while the second is based on the 'local electroneutrality' simplification. Numerical analysis of the results obtained through the two types of simulations shows no significant differences corroborating the usual quasi-neutral modelling.