Hydrologic transport and geochemical reaction model

Introduction to HYDROGEOCHEM
HYDROGEOCHEM Special Features

HYDROGEOCHEM is a coupled model of hydrologic transport and geochemical reaction in saturated-unsaturated media. It is designed to simulate transient and/or steady-state transport of Na, aqueous components and transient and/or steady-state mass balance of Ns adsorbent components and ion-exchange sites. Along the transport path, HYDROGEOCHEM computes the species distribution of N component species, Mx complexed species, My adsorbed species. Mz ion-exchanged species, and Mp potentially precipitated species. The physical, hydrological and chemical settings are as follows.

HYDROGEOCHEM is the only commercially-available model for the simulation of reactive multispecies-multicomponent chemical transport through saturated-unsaturated media. It is not a path model; it is a true transport model coupled with homogeneous and heterogeneous geochemical reactions.


The purpose of HYDROGEOCHEM is to simulate transient and/or steady-state density-dependent flow fields and temperature distribution and to compute transient and/or steady-state distribution of reactive multispecies chemical concentrations in subsurface media.

HYDROGEOCHEM computes and predicts the distribution of pressure head, moisture content, flow velocity, and total head over a three-dimensional plane in either completely saturated, completely unsaturated, partially unsaturated, or partially saturated subsurface media. It also computes and predicts the spatial-temporal distribution of multi-chemical components. The media may consist of as many types of soils and geologic units as desired with different material properties. Each soil type may be isotropic or anisotropic. The processes governing the distribution of chemical distribution include (1) geochemical equilibrium of aqueous complexation, reduction-oxidation, sorption, and precipitation and dissolution, and (2) hydrological transport by flow advection, dispersion, and effect of unsaturation.

The generalized Richards' equation and Darcy's law governing pressure distribution and water flow in saturated-unsaturated media are simulated with the Galerkin finite-element method subject to appropriate initial and four types of boundary conditions. The hydrological transport equations (a set of PDEs) are derived based on the principle of conservation of mass, and the geochemical equilibrium equations (a set of AEs) are derived based on the mass balance and mass action. The coupled set of PDEs and AEs are simulated with either the conventional finite-element methods or the hybrid Langrangian-Eulerian finite-element method with peak capturing scheme subject to appropriate initial and four types of boundary conditions. Hexahedral elements, triangular prism, and tetrahedral elements are used to facilitate the discretization of the region of interest.


The special features of HYDROGEOCHEM are its flexibility and versatility in modeling as wide a range of problems as possible. The model is designed to (1) treat heterogeneous and anisotropic media, (2) consider spatially and temporally-distributed as well as point sources/sinks, (3) accept the prescribed initial conditions or obtain initial conditions by simulating the steady-state version of the system under consideration, (4) deal with prescribed transient concentrations distributed over a Dirichlet boundary, (5) handle time-dependent fluxes over variable boundaries, (6) deal with time-dependent total fluxes over Cauchy boundaries, (7) include the off-diagonal dispersion coefficient tensor components in the governing equation for dealing with cases when the coordinate system does not coincide with the principal directions of the dispersion coefficient tensor, (8) provide two options for treating the mass matrix - consistent and lumping, (9) give three options (exact relaxation, under- and over-relaxation) for estimating the nonlinear matrix, (10) include two options (direct solution with Gaussian elimination method and successive point iterations) for solving the linearized matrix equations, (11) include both quadrilateral and triangular elements to facilitate the discretization of the region, (12) automatically reset time step size when boundary conditions or sources/sinks change abruptly, and (13) include simultaneous chemical processes of aqueous complexation, precipitation/dissolution, adsorption, ion exchange, redox, and acid-base reactions.


(1) Geometry in terms of nodes and elements, and boundaries in terms of nodes and segments; (2) soil properties including (a) saturated hydraulic conductivities or permeabilities; (b) compressibility of water and the media, respectively; (c) bulk density; (d) three soil characteristic curves for each type of soil or geologic unit which are the retention curve, relative conductivity vs head curve, and water capacity curve; (e) effect porosity; and (f) dispersivities, and effective molecular diffusion coefficient for each soil type or geologic unit; (3) initial distribution of pressure head over the region of interest; (4) net precipitation, allowed ponding depth, potential evaporation, and allowed minimum pressure head in the soil; (5) prescribed head on Dirichlet boundaries; (6) prescribed fluxes on Cauchy and/or Neumann boundaries; (7) artificial withdrawals or injections of water; (8) number of chemical components as well as chemical species and their thermodynamic data base; (9) artificial sources/sinks of water and all chemical components; (10) prescribed total concentrations of all chemical components on Dirichlet boundaries; (11) prescribed fluxes of all chemical components on variable boundaries; and (12) initial distribution of all chemical component concentrations. All inputs in items 4 through 11 can be time-dependent or constant with time.


(1) pressure head, total head, moisture content, and flow velocity over the two-dimensional grid at any desired time; (2) water fluxes through all types of boundaries and amount of water accumulated in the media at any desired time; (3) distribution of total analytical concentrations, total dissolved concentrations, total sorbed concentrations, total precipitated concentrations, and free ion concentrations of all chemical components over a three-dimensional grid at any desired time; (4) amount of waste fluxes through the variable boundary; and (5) equivalent kds as a function of time and space in the region of interest.

HYDROGEOCHEM Requirements: 486/Pentium (Pentium recommended) with 16 MB RAM and FORTRAN Compiler, any Workstation, e.g., IBM RS6000, DEC Alpha, Silicon Graphics, Sun SparcStation, and HP 9000 Series.

Incredible NetworksΤελευταία Ενημέρωση 27 Ιουλίου 2004 - Last Revised on July 27th 2004
Προηγούμενη σελίδα
Copyright 1998-2005 MP & Associates - Φορμίωνος 119-121 - ΑΘΗΝΑ 16121 - Τηλ: (210) 7600955 - Fax: (210) 7600956