HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS		QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

This Article Figures Only Full Text Free Full Text (PDF) Free Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in ISI Web of Science Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via ISI Web of Science (10) Citing Articles via Google Scholar Google Scholar Articles by Hansson, K. Articles by van Genuchten, M. Th. Search for Related Content PubMed Articles by Hansson, K. Articles by van Genuchten, M. Th. GeoRef GeoRef Citation Agricola Articles by Hansson, K. Articles by van Genuchten, M. Th. Related Collections Heat Transport Numerical Solutions Coupled Flow/Transport Models

ORIGINAL RESEARCH

Water Flow and Heat Transport in Frozen Soil

Numerical Solution and Freeze–Thaw Applications

Klas Hansson^*,a, Jirka imnek^b, Masaru Mizoguchi^c, Lars-Christer Lundin^a and Martinus Th. van Genuchten^d

^a Dep. of Earth Sciences, Uppsala Univ., Uppsala, Sweden
^b Dep. of Environmental Science, Univ. of California, Riverside, CA
^c Dep. of Biological and Environmental Engineering, Univ. of Tokyo, Japan
^d George E. Brown, Jr. Salinity Laboratory, USDA-ARS, Riverside, CA
^* Corresponding author (klas.hansson{at}hyd.uu.se).

Received 19 August 2003.

A new method is presented to account for phase changes in a fully implicit numerical model for coupled heat transport and variably saturated water flow involving conditions both above and below zero temperature. The method is based on a mixed formulation for both water flow and heat transport similar to the approach commonly used for the Richards equation. The approach enabled numerically stable, energy- and mass-conservative solutions. The model was evaluated by comparing predictions with data from laboratory column freezing experiments. These experiments involved 20-cm long soil columns with an internal diameter of 8 cm that were exposed at the top to a circulating fluid with a temperature of –6°C. Water and soil in the columns froze from the top down during the experiment, with the freezing process inducing significant water redistribution within the soil. A new function is proposed to better describe the dependency of the thermal conductivity on the ice and water contents of frozen soils. Predicted values of the total water content compared well with measured values. The model proved to be numerically stable also for a hypothetical road problem involving simultaneous heat transport and water flow. The problem was simulated using measured values of the surface temperature for the duration of almost 1 yr. Since the road was snow-plowed during winter, surface temperatures varied more rapidly, and reached much lower values, than would have been the case under a natural snow cover. The numerical experiments demonstrate the ability of the code to cope with rapidly changing boundary conditions and very nonlinear water content and pressure head distributions in the soil profile.

Abbreviations: SVAT, soil–vegetation–atmosphere transfer

This article has been cited by other articles:

				J. Simunek, M. Th. van Genuchten, and M. Sejna Development and Applications of the HYDRUS and STANMOD Software Packages and Related Codes Vadose Zone J., May 27, 2008; 7(2): 587 - 600. [Abstract] [Full Text] [PDF]

				R. P. Daanen, D. Misra, and H. Epstein Active-Layer Hydrology in Nonsorted Circle Ecosystems of the Arctic Tundra Vadose Zone J., October 8, 2007; 6(4): 694 - 704. [Abstract] [Full Text] [PDF]