www.the-cryosphere-discuss.net/4/1243/2010/ doi:10.5194/tcd-4-1243-2010 © Author(s) 2010. This work is distributed under the Creative Commons Attribution 3.0 License. An energy-conserving model of freezing variably-saturated soil 1Department of Civil and Environmental Engineering, University of Trento, Trento, Italy 2Environment Canada – National Hydrology Research Centre, Saskatoon, Saskatchewan, Canada 3Department of Geography, University of Zurich, Winterthurerstrasse 190 Zurich, Switzerland *now at: Mountain-eering srl, Via Siemens 19 Bolzano, Italy Abstract. In this paper we provide a method for solving the energy equation in freezing soil. The solver is linked with the solution of Richards equation, and therefore able to approximate water movement near the liquid-solid phase transition. The equations show non-linear characteristics causing oscillatory behavior in the solution close to the phase transition, when normal methods of iterative integration, as Newton or Picard, are used. Thus, a globally convergent Newton method has been implemented to achieve convergence. The method is tested by comparison with an analytical solution to the Stefan problem and by comparison with experimental data derived from literature. Discussion Paper (PDF, 1871 KB) Interactive Discussion (Closed, 6 Comments) Final Revised Paper (TC) Special Issue Citation: Dall'Amico, M., Endrizzi, S., Gruber, S., and Rigon, R.: An energy-conserving model of freezing variably-saturated soil, The Cryosphere Discuss., 4, 1243-1276, doi:10.5194/tcd-4-1243-2010, 2010. Bibtex EndNote Reference Manager XML |