Premium
Locally conservative groundwater flow in the continuous G alerkin method using 3‐ D prismatic patches
Author(s) -
Wu Qiang,
Zhao Yingwang,
Lin YuFeng F.,
Xu Hua
Publication year - 2016
Publication title -
water resources research
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.863
H-Index - 217
eISSN - 1944-7973
pISSN - 0043-1397
DOI - 10.1002/2016wr018967
Subject(s) - streamlines, streaklines, and pathlines , conservation of mass , sink (geography) , advection , groundwater flow , geology , flow (mathematics) , finite element method , geometry , groundwater , algorithm , mathematics , computer science , mechanics , engineering , geotechnical engineering , aquifer , physics , structural engineering , geography , cartography , thermodynamics
A new procedure has been developed to improve the velocity field computed by the continuous Galerkin finite element method (CG). It enables extending the postprocessing algorithm proposed by Cordes and Kinzelbach (1992) to three‐dimensional (3‐D) models by using prismatic patches for saturated groundwater flow. This approach leverages a dual mesh to preserve local mass conservation and provides interpolated velocities based on consistent fluxes. To develop this 3‐D approach, a triangular conservative patch is introduced by computing not only advection fluxes, but also vertical infiltrations, storage changes, and other sink or source terms. This triangular patch is then used to develop a prismatic patch, which consists of subprisms in two layers. By dividing a single two‐layer patch into two separate one‐layer patches, two dimensional (2‐D) algorithms can be applied to compute velocities. As a consequence, each subelement is able to preserve local mass conservation. A hypothetical 3‐D model is used to evaluate the precision of streamlines and flow rates generated by this approach and the FEFLOW simulation program.