The model was developed from numerical and analytical solutions of the diffusion forms of the mass continuity equation and the first law of thermodynamics. Analytical solutions provided a functional framework to evaluate nonlinear interactions obtained in the numerical solutions. Dimensional analysis was used to define the relationships between the soil properties used in the model.