ICE-1 is a research model for the analysis of coupled flow of water, heat and solute in unsaturated, partially frozen soils, including heave effects. The model is based on solving simultaneously equations representing the suction-water content relationship, the hydraulic conductivity function, equation of state for energy of the liquid phase, approximate time integral of Darcy's law, approximate time integral of the Fournier heat diffusion equation, and the mass balance of solutes. The soil column is divided into elements, each one centimeter thick. It is assumed that the soil hydraulic and thermal properties are known, as well as the initial condition. The soil surface temperature should be specified at all times. The code employs a temperature distribution that decreases as a linear function of time with an upper limit on the (negative) temperature. Two conditions may be specified at the lower boundary: either a zero water flow boundary or a flux boundary. A specified temperature boundary condition is assume at the lower boundary. The program computes and displays for given depth, temperature, liquid water content, total water content, salt content, heaving pressure and overburden pressure. It also computes latent heat and total heave.

The public domain, DOS-based program supports interactive data entry and correction, graphic display of results, and writing results to a file for printing and inclusion in reports. ICE-1 is distributed complete with QuickBASIC source code and executable image. The program documentation consists of installation instructions, theory, and an example problem.

SYSTEM REQUIREMENTS

Intel 80i86 based computer, 640 Kb conventional RAM, DOS 3.0 or higher, CGA or EGA graphics; math coprocessor recommended.

Developers: A.I. El-Kadi (IGWMC) and J.W. Cary (Battelle Pacific Northwest Laboratory)