A Modified Effective Binary Diffusion Model

Uphill diffusion, the diffusion of a component up its own concentration gradient, has often been observed in experimental studies of diffusion i synthetic multicomponent systems and natural silicate melts. To apply experimental data to natural processes, it is necessary toquantify this effect. A modified effective binary diffusion model is developed in this paper. The diffusive flux of a component is assumed tobe proportional to its activity gradient (instead of concentration gradient). The activity coefficient is related to its concentration at quasi-equilibrium distribution i a compositional continuum with the help of the two-liquid partition concept. The concentration at quasi-equilibrium distribution is approximately xpressed as a simple function of SiO 2 or SiO2+A1203 concentration. The diffusion of SiO 2 or SiO2+A1203 isassumed tobe effectively binary. With these assumptions, diffusion of a major or trace component can be described by two adjustable parameters, •D (the intrinsic effective binary diffusivity) and Cf (a parameter r lated to two-liquid partition coefficients). The simple model fits well the experimental uphill diffusion profiles produced by diffusive crystal dissolution i to an andesitic melt using constant diffusivities when compositional v riation across the whole profile is small. The parameters fit to the profiles are generally physically meaningful. The model is able to predict he direction of diffusion (occurrence orabsence of uphill diffusion) in previous experiments and can be used to simulate the development of chemical concentration a d isotopic ratio profiles during diffusion. When compositional dependence of diffusivities is allowed, the model may be applied to systems with greater compositional v riation. This model can be tested in several ways and should be useful in studying diffusion in a variety of natural systems.