Flow of groundwater with variable density and viscosity, Atikokan Research Area, Canada

Research output: Contribution to journalArticle

10 Citations (Scopus)

Abstract

Flow of groundwater with variable density and viscosity was simulated at the Atikokan Research Area (ARA) in northwestern Ontario, Canada. An empirical viscosity-concentration equation was modified to include total-dissolved-solids (TDS) data from the ARA. The resulting equation was used successfully to estimate reasonably accurate viscosity values over the expected range of temperature and concentration, in comparison with experimental values derived for sodium chloride solutions. A three-dimensional finite-element code, MOTIF, developed by Atomic Energy of Canada Limited, was used in the simulations. The inclusion of the effects of depth-increasing temperature and TDS-dependent fluid-density distribution, while maintaining only a temperature-dependent viscosity relationship in a simulation, resulted in a more penetrative flow against expected buoyancy effects (i.e., the physics of the system was not honored). Accounting for concentration in the viscosity equation caused water to be less penetrative and more in accordance with the expected physics of the system. A conclusion is that fluid concentration should be considered simultaneously in calculating the density and viscosity of a fluid during modeling of variable-density flow in areas underlain by fluids with high TDS. Results of simulations suggest that both flow directions and magnitudes should be employed simultaneously during the calibration of a model. Large-scale groundwater movement in the ARA may be analyzed with carefully selected vertical no-flow boundaries. By incorporating the geothermal temperature gradient, groundwater recharge increases by 12%; thus, this gradient plays a significant role in groundwater flow at the ARA. Variability in the fluid concentration at the ARA neither decreases nor increases recharge into the groundwater system. The hypothesis that an isolated continuous regional flow system may exist at depth in the ARA is not supported by these simulations.

Original languageEnglish
Pages (from-to)193-203
Number of pages11
JournalHydrogeology Journal
Volume6
Issue number2
DOIs
StatePublished - 1 Jan 1998

Fingerprint

viscosity
groundwater
fluid
simulation
recharge
physics
temperature
sodium chloride
density current
temperature gradient
buoyancy
groundwater flow
calibration
modeling
energy
water
effect

Cite this

@article{8619198d87b54d2b88bd494c41c69253,
title = "Flow of groundwater with variable density and viscosity, Atikokan Research Area, Canada",
abstract = "Flow of groundwater with variable density and viscosity was simulated at the Atikokan Research Area (ARA) in northwestern Ontario, Canada. An empirical viscosity-concentration equation was modified to include total-dissolved-solids (TDS) data from the ARA. The resulting equation was used successfully to estimate reasonably accurate viscosity values over the expected range of temperature and concentration, in comparison with experimental values derived for sodium chloride solutions. A three-dimensional finite-element code, MOTIF, developed by Atomic Energy of Canada Limited, was used in the simulations. The inclusion of the effects of depth-increasing temperature and TDS-dependent fluid-density distribution, while maintaining only a temperature-dependent viscosity relationship in a simulation, resulted in a more penetrative flow against expected buoyancy effects (i.e., the physics of the system was not honored). Accounting for concentration in the viscosity equation caused water to be less penetrative and more in accordance with the expected physics of the system. A conclusion is that fluid concentration should be considered simultaneously in calculating the density and viscosity of a fluid during modeling of variable-density flow in areas underlain by fluids with high TDS. Results of simulations suggest that both flow directions and magnitudes should be employed simultaneously during the calibration of a model. Large-scale groundwater movement in the ARA may be analyzed with carefully selected vertical no-flow boundaries. By incorporating the geothermal temperature gradient, groundwater recharge increases by 12{\%}; thus, this gradient plays a significant role in groundwater flow at the ARA. Variability in the fluid concentration at the ARA neither decreases nor increases recharge into the groundwater system. The hypothesis that an isolated continuous regional flow system may exist at depth in the ARA is not supported by these simulations.",
author = "Duke Ophori",
year = "1998",
month = "1",
day = "1",
doi = "10.1007/s100400050144",
language = "English",
volume = "6",
pages = "193--203",
journal = "Hydrogeology Journal",
issn = "1431-2174",
publisher = "Springer Heidelberg",
number = "2",

}

Flow of groundwater with variable density and viscosity, Atikokan Research Area, Canada. / Ophori, Duke.

In: Hydrogeology Journal, Vol. 6, No. 2, 01.01.1998, p. 193-203.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Flow of groundwater with variable density and viscosity, Atikokan Research Area, Canada

AU - Ophori, Duke

PY - 1998/1/1

Y1 - 1998/1/1

N2 - Flow of groundwater with variable density and viscosity was simulated at the Atikokan Research Area (ARA) in northwestern Ontario, Canada. An empirical viscosity-concentration equation was modified to include total-dissolved-solids (TDS) data from the ARA. The resulting equation was used successfully to estimate reasonably accurate viscosity values over the expected range of temperature and concentration, in comparison with experimental values derived for sodium chloride solutions. A three-dimensional finite-element code, MOTIF, developed by Atomic Energy of Canada Limited, was used in the simulations. The inclusion of the effects of depth-increasing temperature and TDS-dependent fluid-density distribution, while maintaining only a temperature-dependent viscosity relationship in a simulation, resulted in a more penetrative flow against expected buoyancy effects (i.e., the physics of the system was not honored). Accounting for concentration in the viscosity equation caused water to be less penetrative and more in accordance with the expected physics of the system. A conclusion is that fluid concentration should be considered simultaneously in calculating the density and viscosity of a fluid during modeling of variable-density flow in areas underlain by fluids with high TDS. Results of simulations suggest that both flow directions and magnitudes should be employed simultaneously during the calibration of a model. Large-scale groundwater movement in the ARA may be analyzed with carefully selected vertical no-flow boundaries. By incorporating the geothermal temperature gradient, groundwater recharge increases by 12%; thus, this gradient plays a significant role in groundwater flow at the ARA. Variability in the fluid concentration at the ARA neither decreases nor increases recharge into the groundwater system. The hypothesis that an isolated continuous regional flow system may exist at depth in the ARA is not supported by these simulations.

AB - Flow of groundwater with variable density and viscosity was simulated at the Atikokan Research Area (ARA) in northwestern Ontario, Canada. An empirical viscosity-concentration equation was modified to include total-dissolved-solids (TDS) data from the ARA. The resulting equation was used successfully to estimate reasonably accurate viscosity values over the expected range of temperature and concentration, in comparison with experimental values derived for sodium chloride solutions. A three-dimensional finite-element code, MOTIF, developed by Atomic Energy of Canada Limited, was used in the simulations. The inclusion of the effects of depth-increasing temperature and TDS-dependent fluid-density distribution, while maintaining only a temperature-dependent viscosity relationship in a simulation, resulted in a more penetrative flow against expected buoyancy effects (i.e., the physics of the system was not honored). Accounting for concentration in the viscosity equation caused water to be less penetrative and more in accordance with the expected physics of the system. A conclusion is that fluid concentration should be considered simultaneously in calculating the density and viscosity of a fluid during modeling of variable-density flow in areas underlain by fluids with high TDS. Results of simulations suggest that both flow directions and magnitudes should be employed simultaneously during the calibration of a model. Large-scale groundwater movement in the ARA may be analyzed with carefully selected vertical no-flow boundaries. By incorporating the geothermal temperature gradient, groundwater recharge increases by 12%; thus, this gradient plays a significant role in groundwater flow at the ARA. Variability in the fluid concentration at the ARA neither decreases nor increases recharge into the groundwater system. The hypothesis that an isolated continuous regional flow system may exist at depth in the ARA is not supported by these simulations.

UR - http://www.scopus.com/inward/record.url?scp=0012064525&partnerID=8YFLogxK

U2 - 10.1007/s100400050144

DO - 10.1007/s100400050144

M3 - Article

AN - SCOPUS:0012064525

VL - 6

SP - 193

EP - 203

JO - Hydrogeology Journal

JF - Hydrogeology Journal

SN - 1431-2174

IS - 2

ER -