HydroGeoSphere

HydroGeoSphere

Developer(s)	Aquanty
Stable release	2013
Operating system	Windows
Type	Hydrogeology software
Website	http://www.aquanty.com/

HydroGeoSphere (HGS) is a 3D control-volume finite element groundwater model, and is based on a rigorous conceptualization of the hydrologic system comprising of surface and subsurface flow regimes. The model is designed to take into account all key components of the hydrologic cycle. For each time step, the model solves surface and subsurface flow, solute and energy transport equations simultaneously, and provides a complete water and solute balance.

Governing Equations

In order to accomplish the integrated analysis, HydroGeoSphere utilizes a rigorous, mass conservative modeling approach that fully couples the surface flow and transport equations with the 3-D, variably saturated subsurface flow and transport equations. This approach is significantly more robust than previous conjunctive approaches that rely on linkage of separate surface and subsurface modeling codes.

Groundwater Flow

HydroGeoSphere assumes that the subsurface flow equation in a porous medium is always solved during a simulation, either for fully-saturated or variably-saturated flow conditions. The subsurface flow equation can be expanded to incorporate discrete fractures, a second interacting porous continuum, wells, tile drains and surface flow. The following assumptions are made for subsurface flow:

• The fluid is essentially incompressible.
• The porous medium and fractures (or macropores), if present, are non-deformable.
• The system is under isothermal conditions.
• The air phase is infinitely mobile.

The Richards’ equation is used to describe three-dimensional transient subsurface flow in a variably-saturated porous medium:

$-\nabla \cdot (w_{m}{\textbf {q}})+\sum \Gamma _{eq}\pm Q=w_{m}{\frac {\partial }{\partial t}}(\theta _{s}S_{w})$

The fluid flux, ${\textbf {q}}$, is represented by the Darcy's law shown as:

${\textbf {q}}=-{\textbf {K}}\cdot k_{r}\nabla (\psi +z)$

where

$w_{m}$ is the volumetric fraction of the total porosity occupied by the porous medium

$\Gamma _{ex}$ is the internal fluid exchange rate (e.g. surface water, wells, and tile drains)

$Q$ is the external fluid outside of the model domain

$\theta _{s}$ is the saturated water content

$S_{w}$ is the degree of saturation

${\textbf {K}}$ is the hydraulic conductivity tensor

$k_{r}$ is the relative permeability of the medium calculated as a function of saturation

$\psi$ is the pressure head

$z$ is the elevation head

Surface Water Flow

Areal surface water flow is represented in HydroGeoSphere by a two-dimensional depth-averaged flow equation, which is the diffusion-wave approximation of the Saint Venant equation for surface water flow. HydroGeoSphere's surface water flow component is implemented with the following assumptions:

• Depth-averaged flow velocities
• Hydrostatic pressure distribution vertically
• Mild slope
• Dominant bottom shear stresses.

The surface flow components are solved by the following three equations, which are given by the following mass balance equation:

${\frac {\partial \phi _{o}h_{o}}{\partial t}}+{\frac {\partial {\bar {v}}_{xo}d_{o}}{\partial x}}+{\frac {\partial {\bar {v}}_{yo}d_{o}}{\partial y}}+d_{o}\Gamma _{o}\pm Q_{o}=0$

coupled with the momentum equations, neglecting inertia terms, for the x-direction:

${\bar {v}}_{ox}=-K_{ox}{\frac {\partial h_{o}}{\partial x}}$

and for the y-direction:

${\bar {v}}_{oy}=-K_{oy}{\frac {\partial h_{o}}{\partial y}}$

where

$\phi _{o}$ is the surface flow domain porosity

$h_{o}$ is the water surface elevation

${\bar {v}}_{xo}$ and ${\bar {v}}_{yo}$ are the vertically averaged flow velocities in the x and y directions

$d_{o}$ is the depth of surface water flow

$\Gamma _{o}$ is the internal fluid exchange

$Q_{o}$ is the external fluid exchange.

The surface conductances, $K_{ox}$ and $K_{oy}$ are approximated by either the Manning or Chezy equation.

Solute Transport

Three-dimensional transport of solutes is described by the modified reactive transport advective-dispersion equation:

$-\nabla \cdot w_{m}({\textbf {q}}C-\theta _{s}S_{w}{\textbf {D}}\nabla C)+[w_{m}\theta _{s}S_{w}R\lambda C]_{par}+\sum {\Omega _{ex}+Q_{c}}=w_{m}\left[{\frac {\partial (\theta _{s}S_{w}RC)}{\partial t}}+\theta _{s}S_{w}R\lambda C\right]$

where

$C$ is the solute concentration

$\lambda$ is the first-order decay constant

$Q_{c}$ is the external source or sink term

$\Omega$ is the internal solute transfer between domains

$R$ is the retardation factor

${\textbf {D}}$ is the diffusion coefficient.

$par$ designates parent species for the case of a decay chain

Heat Transport

Graf [2005] incorporated heat transport within the saturated-zone flow regime into HydroGeoSphere together with temperature-dependent fluid properties, such as viscosity and density. The model’s capability was successfully demonstrated for the case of thermohaline flow and transport in porous and fractured porous media [Graf and Therrien, 2007]. This work extends the model’s capability to include thermal energy transport in the unsaturated zone and in the surface water, which is considered a key step in the linkage between the atmospheric and hydrologic systems. Surface heat fluxes from atmospheric inputs are an important source/sink of thermal energy, especially to the surface water system. As such, surface heat fluxes across the land surface were also incorporated into HydroGeoSphere. A complete description of the physical processes and governing flow and solute transport equations that form the basis of HydroGeoSphere can be found in Therrien et al. [2007] and therefore will not be presented here.

The general equation for variably saturated subsurface thermal energy transport following Molson et al. [1992] is given by:

$-\nabla \cdot {\Big (}{\textbf {q}}\rho _{w}c_{w}T-(k_{b}+c_{b}\rho _{b}{\textbf {D}})\nabla T{\Big )}+\Omega _{o}\pm Q_{T}={\frac {\partial \rho _{b}c_{b}T}{\partial t}}$

where

$\rho$ is the density

$c$ is the heat capacity

$T$ is the temperature of the bulk subsurface

$k$ is the thermal conductivity

$D$ is the thermal dispersion term

$Q_{T}$ is the thermal source/sink

$\Omega _{o}$ is the thermal interactions between the surface and subsurface

$Q_{T}$ is the external thermal interactions

Surface-Subsurface Coupling

Integrated surface/subsurface flow with evapotranspiration and precipitation processes.

The 2-D areal surface flow modules of HydroGeoSphere follow the same conventions for spatial and temporal discretizations as those used by the subsurface modules. The surface flow equation is solved on a 2-D finite-element mesh stacked upon a subsurface grid when solving for both domains (i.e. the x- and y-locations of nodes are the same for each layer of nodes). For superposition, the grid generated for the subsurface domain is mirrored areally for the surface flow nodes, with surface flow node elevations corresponding to the top elevation of the topmost active layer of the subsurface grid. Note that surface flow node elevations may vary substantially to conform with topography. However, the assumptions of small slope inherent in the diffusion-wave equation will not allow for modeling of inertial effects.

The discretized surface equation is coupled with the 3-D subsurface flow equation via superposition (common node approach) or via leakage through a surficial skin layer (dual node approach). For both approaches, fully implicit coupling of the surface and subsurface flow regimes provides an integral view of the movement of water, as opposed to the traditional division of surface and subsurface regimes. Flux across the land surface is, therefore, a natural internal process allowing water to move between the surface and subsurface flow systems as governed by local flow hydrodynamics, instead of using physically artificial boundary conditions at the interface. When the subsurface connection is provided via superposition, HydroGeoSphere adds the surface flow equation terms for the 2-D surface mesh to those of the top layer of subsurface nodes. In that case, the fluid exchange flux, which contains leakance term does not need to be explicitly defined.

Features

The HGS model, developed jointly at the University of Waterloo and the University of Laval, is a three-dimensional control-volume finite element simulator which is designed to simulate the entire terrestrial portion of the hydrologic cycle. It uses a globally implicit approach to simultaneously solve the 2D diffusive-wave equation and the 3D form of Richards’ equation. HGS also dynamically integrates key components of the hydrologic cycle such as evaporation from bare soil and water bodies, vegetation-dependent transpiration with root uptake, snowmelt and soil freeze/thaw. Features such as macro pores, fractures, and tile drains can either be incorporated discretely or using a dual-porosity, dual permeability formulation.

HydroGeoSphere is a powerful numerical simulator specifically developed for supporting water resource and engineering projects pertaining to hydrologic systems with surface and subsurface flow and transport components. In terms of simulation capability and computational aspects, the HydroGeoSphere code has the following attributes:

Fluid flow • Complete hydrologic cycle modeling using detailed physics of surface and subsurface flow in one integrated code. The surface regime can be represented as 2-D areal flow for the entire surface or as 2-D runoff into 1-D channels. The subsurface regime consists of 3-D unsaturated/saturated flow. Both regimes naturally interact with each other through considerations of the physics of flow between them. • Physically-based accounting of all components of the hydrologic cycle water budget. • Accurate delineation and tracking of the water table position, taking into account flow in the unsaturated zone, delayed yield and vertical flow components. • Handling of non-ponding or prescribed ponding recharge conditions. • Automatic handling of seepage face conditions at the land-atmosphere interface. • Automatic and correct apportioning of the total flow rate of a multi-layer well to the well nodes, including the simulation of water flow and solute/temperature mixing within the water column in the well. • Accommodation of wellbore storage. CHAPTER 1. INTRODUCTION 6 • Arbitrarycombinationsofporous,discretely-fractured,dual-porosityanddual-permeability media for the subsurface. Mass and thermal energy transport • Capability of modeling non-reactive and reactive chemical species transport in the associated surface and subsurface flow fields, including solute interactions between the surface and subsurface flow regimes. • Calculation of temperatures in the surface and subsurface flow regimes as driven by air temperature and incoming solar radiation, accounting for land surface-atmospheric thermal interactions. • Accurate handling of fluid and mass/thermal energy exchanges between fractures and matrix including matrix diffusion effects and solute/thermal energy advection in the matrix. • Straight or branching decay chains representing degradation reactions. Numerical methods • The fully-implicit coupling approach used by the code provides for a robust mass conserved solution scheme, which is essential for systems with strong interactions between regimes. • Parallelization employing OpenMP. • Advanced computational algorithms and a flexible, user-programmable interface allow the code to perform unprecedented, fully-integrated, 3-D simulation/animation on a personal computer. • Fluid, solute and energy balance tracking. • Unstructured finite-element grids. • Axi-symmetric grid option. • 7-point finite-difference option. • 8-node block or 6-node prism elements, 4-node tetrahedral elements, 3- and 4-node plate elements for fractures and surface water, and 2-node line elements for wells, storm and sanitary sewers, water supply mains and tile drains or other types of linear infrastructure features. • Adaptive time-stepping schemes with automatic generation and control of time steps, sub-timing option at select nodes where flow and transport is rapid. • Straightforward organization and control of simulation output. CHAPTER 1. INTRODUCTION 7 • Robust and efficient ILU-preconditioned iterative sparse-matrix solver. • Robust and efficient Newton-Raphson linearization option. • Flexible pre- and post-processing capabilities. For field applications and research investigations, HydroGeoSphere can be used to perform event-based and continuous simulations on widely varying spatial scales ranging from single soil column profiles to large-scale basins, which may include several catchments. Examples of field applications of HydroGeoSphere include: • Integrated water resource assessment. • Watershed hydrologic analysis, including impacts of land-use or climate-change impacts on both surface and subsurface water. • Floodplain hydrologic analysis. • Fluvial hydraulic analysis. • Contaminant migration and fate in both surface and subsurface water.

Some features include:

• 2D overland flow
• 2D stream network flow
• 3D variably saturated flow
• 3D flow in fractures and conduits
• Advection/dispersion modeling
• Reactions
• Density dependent flow
• Land Surface Processes

Other Groundwater Models

• MODFLOW
• FEFLOW

External links

• HGS Homepage
• HGS Review Article