THE VELOCITY ORIENTED APPROACH REVISITED
Journal Title: Biuletyn Państwowego Instytutu Geologicznego - Year 2010, Vol 441, Issue 441
Abstract
A great deal of hydrogeological situations requires an extremely accurate calculation of the 3-dimensional groundwater discharge rates in the subsoil. Examples are: hydrology of wetlands, water balances of aquatic ecosystems depending on groundwater recharge, river-groundwater interaction, advective transport of pollution underneath waste disposal sites, particle trajectories in aquifer-aquitard systems with contrasting heterogeneities and many others. Numerical determination of the vertical groundwater velocity is a notoriously difficult problem. In nature this component may be two or three orders of magnitude smaller than the horizontal velocity components. In such cases application of Darcy’s law to the numerically calculated hydraulic heads obtained from a finite difference or finite element model may lead to relatively inaccurate vertical velocities. More specifically, when estimating vertical velocity components in cases where the Dupuit approximation – negligible vertical head gradient – holds, numerical differentiation of hydraulic heads yields zero vertical velocity. In the 1980s of the last century Zijl and Nawalany proposed to invert the order of calculating the velocity field by eliminating the head from Darcy’s law and to consider the Darcy velocity as the primary variable. For 2-dimensional flow this was already common practice and the challenge was a 3-dimensional extension, which was called the Velocity Oriented Approach (VOA). In two dimensions such methods were conventionally based on a stream function as primary variable. However, at that time application of a 3D stream function was not feasible and, therefore, the Darcy velocity itself was considered as the primary variable. This approach has been proven to yield a high accuracy for all three components of the specific discharge, including the relatively small vertical component, especially in cases where the subsoil is smoothly heterogeneous in the horizontal directions. In the 1990s the mixed-hybrid finite element method was developed. The physical interpretation of this method shows the way how to liberate the VOA from its smoothness requirement by introduction of a practical applicable 3D stream function. In conclusion, the velocity oriented approach indicates a change in paradigm regarding the accurate calculation of specific discharge in groundwater flow.
Authors and Affiliations
Marek NAWALANY, Wouter ZIJL
Keywords
Related Articles
- EP ID EP64691
- DOI -
- Views 31
- Downloads 0
RECORD OF CHANGES IN THE OLIGOCENE–MIOCENE SEDIMENTS OF THE MENILITE–KROSNO SERIES OF THE SKOLE UNIT BASED ON CALCAREOUS NANNOPLANKTON STUDIES – BIOSTRATIGRAPHY AND PALAEOGEOGRAPHICAL IMPLICATIONS (POLISH OUTER CARPATHIANS)
The biostratigraphic investigation included the Menilite–Krosno Series of the Skole Unit overlying the Globigerina Marls. Seven calcareous nannoplankton zones (sensu Martini, 1971) were distinguished in these sediments:...
OCENA PRACY I ZASIĘGU ODDZIAŁYWANIA „DUŻEGO” UJĘCIA WÓD PODZIEMNYCH PO 40 LATACH UŻYTKOWANIA
W artykule przedstawiono wyniki kompleksowych badań hydrogeologicznych w celu oceny sprawności 40-letnich studni, poprawności wyznaczonych w latach 70. XX w. wielkości zasobów eksploatacyjnych „dużego” ujęcia i zasięgu j...
CHANGES OF RECHARGE CONDITIONS OF INFLOW W VI-32 AT WIELICZKA SALT MINE IN THE LIGHT OF LONG-TERM ISOTOPIC INVESTIGATIONS
In a paper the results of long-term isotopic investigations are presented. Both stable isotopes ([sup]18[/sup]O and [sup]2[/sup]H) and tritium concentrations were measured in water appearing in the Z-32 chamber at the V...
GEOSTATISTICS AS A TOOL SUPPORTING RESEARCH ON MINERAL DEPOSITS OF POLAND – MAIN DIRECTIONS OF APPLICATIONS
The history of development of geostatistical methods in the world and their applications in Poland are shortly presented. Main directions of the use of geostatistics in the evaluation of mineral deposits of Poland are di...
MODEL RESEARCH ON THE EFFECTIVENESS OF TECHNICAL METHODS OF LIMITING THE EXPANSION OF THE DEPRESSION CONE CAUSED BY DEWATERING OF THE PLANNED GUBIN OPEN-PIT BROWN COAL MINE
Due to the planned exploitation of the Gubin lignite deposit by opencast method, a numerical hydrogeological model was built. The main goal of the model was to forecast the development of the cone of depression caused by...