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

• Record Type: Journal Article
• Authors: Grenier, Christophe; Anbergen, Hauke; Bense, Victor; Chanzy, Quentin; Coon, Ethan; Collier, Nathaniel; Costard, François; Ferry, Michel; Frampton, Andrew; Frederick, Jennifer; Gonçalvès, Julio; Holmén, Johann; Jost, Anne; Kokh, Samuel; Kurylyk, Barret; McKenzie, Jeffrey; Molson, John; Mouche, Emmanuel; Orgogozo, Laurent; Pannetier, Romain; Rivière, Agnès; Roux, Nicolas; Rühaak, Wolfram; Scheidegger, Johanna; Selroos, Jan-Olof; Therrien, René; Vidstrand, Patrik; Voss, Clifford
• Is Part Of: Advances in water resources. Volume 114(2018)
• Page Start: 196
• Page End: 218
• Publication Date: 2018-04
• Subjects: Numerical simulation -- Code benchmarking -- Thermo-hydrological coupling -- Permafrost -- Sharp interface problems; Hydrology -- Periodicals; Hydrodynamics -- Periodicals; Hydraulic engineering -- Periodicals; 551.48
• DOI: 10.1016/j.advwatres.2018.02.001
• Languages: English
• ISSNs: 0309-1708

Abstract

Highlights: Coupled Thermo-Hydrological codes including latent heat for Cryohydrogeological codes require some level of evaluation. Two 2D benchmark cases were developed for intercomparison of TH codes. Intercomparison results suggest that the codes provide robust results for the test cases considered. Minor discrepancies found in the intercomparison results were studied raising implementation issues. Abstract: 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.