ChemFlux™ is able to model fate and transport through the saturated and unsaturated zones. The process of advection, diffusion, adsorption and decay may be modeled. ChemFlux is also able to refine the solution mesh based on any problem variable. Saturation levels may be imported from SVFlux™. The adaptive mesh refinement abilities of ChemFlux allow the user to specify mesh refinement that follows the front of the contaminant. The mesh then relaxes dynamically as need be to optimize solution time.