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 th...

Full description

Bibliographic Details
Published in:Advances in Water Resources
Main 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
Language:unknown
Published: 2023
Subjects:
Online Access:http://www.osti.gov/servlets/purl/1464021
https://www.osti.gov/biblio/1464021
https://doi.org/10.1016/j.advwatres.2018.02.001
Description
Summary: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.