A mathematical model of coupled heat and moisture flow in soils have been developed. The model includes algorithms for phase change of soil moisture and frost heave and permits several types of boundary and initial conditions. The finite element method of weighted residuals (Galerkin procedure) was chosen to simulate the spatial regime and the Crank-Nicholson method was used for the time domain portion of the model. To facilitate evaluation of the model, the heat and moisture fluxes were essentially decoupled; moisture flux was then simulated accurately, as were heat flux and frost heave in a laboratory test. Comparision of the simulated and experimental data illustrates the importance of unsaturated hydraulic conductivity. It is one parameter which is difficult to measure and for which only a few laboratory test results are available. Therefore, unsaturated hyraulic conductivities calculated in the computer model may be a significant source of error in calculations of frost heave. The algorithm incorporating effects of surcharge and overburden was inconclusively evaluated. Time-dependent frost penetration and frost heave in laboratory specimens were closely simulated with the model. After 10 days of simulation, the computed frost heave was about 2.3 cm va 2.0 and 2.8 cm in two tests. Frost penetration was computed as 15 cm and was measured at 12.0 cm and 12.2 cm in the two laboratory samples after 10 days.

Media Info

  • Features: Appendices; Figures; Tables;
  • Pagination: 49 p.

Subject/Index Terms

Filing Info

  • Accession Number: 00312458
  • Record Type: Publication
  • Report/Paper Numbers: CRREL 80-10
  • Files: TRIS, USDOT
  • Created Date: Aug 5 1980 12:00AM