Mixing-Cell Model for the Diffusion-Advection Equation
The unsteady-state advection-diffusion equation is solved with a mixing-cell model. The results obtained with this simple model are in excellent agreement with the analytical solution when the correct choice of time and space steps is made. The equation describing the one-dimensional transport of a reactive component in porous media is:[more]
where is concentration, is time, is distance, is the dispersion coefficient, is the pore velocity, and is the reaction constant.
The initial and boundary conditions are
The analytical solution to this equation is :
The mixing-cell model is an explicit finite-difference approximation to the diffusion-advection equation. The time derivative is approximated with a forward difference, while a backward difference is used to approximate the space derivative. For the cell, this approximation is written as:
. The left-hand side and the first term on the right-hand side describe the transport according to the mixing-cell model. In order to eliminate the effect of the second right-hand side term, the following relationship must hold :
Taking the dispersion to be proportional to , then , thus
, where because cannot be negative.
The results are shown as a three-dimensional plot of the analytic solution (blue curve) and the mixing-cell model (red markers) for several values of the variables , , , and , which are obtained with , where is the column length. The results of the two methods are in very close agreement.[less]
 M. Th. van Genuchten and W. J. Alves, "Analytical Solutions of the One-Dimensional Convective-Dispersive Solute Transport Equation," US Department of Agriculture, Technical Bulletin No. 1661, 1982. www.ars.usda.gov/sp2UserFiles/Place/20360500/pdf_pubs/P0753.pdf.
 H. C. Van Ommen, "The 'Mixing-Cell' Concept Applied to Transport of Non-reactive and Reactive Components in Soils and Groundwater," Journal of Hydrology, 78(3–4), 1985 pp. 201–213. doi:10.1016/0022-1694(85)90101-5.