John W. Molson | Geology Reports

USGS2018

Groundwater flow and heat transport for systems undergoing freeze-thaw: Intercomparison of numerical simulators for 2D test cases

In high-elevation, boreal and arctic regions, hydrological processes and associated water bodies can be strongly influenced by the distribution of permafrost. Recent field and modeling studies indicate that a fully-coupled multidimensional thermo-hydraulic approach is required to accurately model the evolution of these permafrost-impacted landscapes and groundwater systems. However, the relatively new and complex numerical codes being developed for coupled non-linear freeze-thaw systems require verification. This issue is addressed by means of an intercomparison of thirteen numerical codes for two-dimensional test cases with several performance metrics (PMs). These codes comprise a wide range of numerical approaches, spatial and temporal discretization strategies, and computational efficiencies. Results suggest that the codes provide robust results for the test cases considered and that minor discrepancies are explained by computational precision. However, larger discrepancies are observed for some PMs resulting from differences in the governing equations, discretization issues, or in the freezing curve used by some codes.

Advances in Water Resources↗

USGSDate not supplied

Characterizing the two-dimensional thermal conductivity distribution in a sand and gravel aquifer

Both hydrologic and thermal transport properties play a significant role in the movement of heat through permeable sedimentary material; however, the thermal conductivity is rarely characterized in detailed spatial resolution. As part of a study of the movement of thermal plumes through a sand and gravel aquifer, we have constructed a two-dimensional profile of thermal conductivity. This work consisted of: (i) measuring the thermal conductivity of the soil solids, λ s , for the main stratigraphic units using the steady-state divided-bar apparatus and estimating conductivity from mineral composition; (ii) measuring the volumetric water content and porosity using crosshole ground-penetrating radar; (iii) evaluating four models used to predict the apparent thermal conductivity, λ, of variably saturated soils and selecting the best model using the information-theoretic approach, (iv) calculating the λ field on a 0.25-m square cell grid using measured data and the selected model, and (v) simulating thermal transport within the two-dimensional domain using a finite element numerical model. The apparent thermal conductivity in the saturated aquifer ranges from 2.14 to 2.69 W m −1 K −1 with a mean of 2.42 W m −1 K −1 Numerical simulations show that the heterogeneous thermal conductivity field results in increased thermal dispersion that is most pronounced at the plume front. Our values for λ and λ s may be used for glacial soils with similar mineralogy and texture. Our methods may also be used at other sites to construct the thermal conductivity distribution.

Ontario↗