Directly to content
  1. Publishing |
  2. Search |
  3. Browse |
  4. Recent items rss |
  5. Open Access |
  6. Jur. Issues |
  7. DeutschClear Cookie - decide language by browser settings

Discontinuous Galerkin based Geostatistical Inversion of Stationary Flow and Transport Processes in Groundwater

Ngo, Adrian Q. T.

[thumbnail of Diss_A.Ngo_2015.pdf]
Preview
PDF, English
Download (7MB) | Terms of use

Citation of documents: Please do not cite the URL that is displayed in your browser location input, instead use the DOI, URN or the persistent URL below, as we can guarantee their long-time accessibility.

Abstract

The hydraulic conductivity of a confined aquifer is estimated using the quasi-linear geostatistical approach (QLGA), based on measurements of dependent quantities such as the hydraulic head or the arrival time of a tracer. This requires the solution of the steady-state groundwater flow and solute transport equations, which are coupled by Darcy's law. The standard Galerkin finite element method (FEM) for the flow equation combined with the streamline diffusion method (SDFEM) for the transport equation is widely used in the hydrogeologists' community. This work suggests to replace the first by the two-point flux cell-centered finite volume scheme (CCFV) and the latter by the Discontinuous Galerkin (DG) method. The convection-dominant case of solute (contaminant) transport in groundwater has always posed a special challenge to numerical schemes due to non-physical oscillations at steep fronts. The performance of the DG method is experimentally compared to the SDFEM approach with respect to numerical stability, accuracy and efficient solvability of the occurring linear systems. A novel method for the reduction of numerical under- and overshoots is presented as a very efficient alternative to local mesh refinement. The applicability and software-technical integration of the CCFV/DG combination into the large-scale inversion scheme mentioned above is realized. The high-resolution estimation of a synthetic hydraulic conductivity field in a 3-D real-world setting is simulated as a showcase on Linux high performance computing clusters. The setup in this showcase provides examples of realistic flow fields for which the solution of the convection-dominant transport problem by the streamline diffusion method fails.

Translation of abstract (German)

Die hydraulische Leitfähigkeit in einem gespannten Grundwasserleiter wird mit dem QLGA Inversionsverfahren (Quasi-linearer geostatistischer Ansatz) geschätzt, indem man Messungen von abhängigen Größen wie etwa des hydraulischen Drucks und der Ankunftszeit eines eingeleiteten Tracers (z.B. gefärbter Indikator) auswertet. Dabei müssen die stationären Gleichungen für Strömung und Transport, die über Darcys Gesetz gekoppelt sind, numerisch gelöst werden. Bei Hydrogeologen ist die Standard Finite Elemente Methode (FEM) zur Lösung der Strömungsgleichung und die Methode der Stromlinien-Diffusion (SDFEM) zur Lösung der Transportgleichung weit verbreitet. In dieser praxisorientierten Arbeit wird dieser Ansatz durch das zellzentrierte Finite-Volumen-Verfahren (CCFV) kombiniert mit dem unstetigen Galerkin-Verfahren (DG) ersetzt. Der konvektionsdominante Stofftransport im Grundwasser stellt seit jeher eine Herausforderung an numerische Methoden dar. Das DG-Verfahren wird mit dem SDFEM-Verfahren verglichen, und zwar hinsichtlich Stabilität, Genauigkeit und der effizienten Lösbarkeit der auftretenden linearen Gleichungssysteme. Eine neue Methode zur Dämpfung von numerischen Unter- und Überschwinger wird als eine sehr effiziente Alternative zur lokalen Gitterverfeinerung präsentiert. Die Anwendbarkeit der Kombination CCFV/DG auf das erwähnte großskalige Inversionsverfahren wird realisiert. In einer dreidimensionalen Versuchsumgebung wird ein synthetisch generiertes Parameterfeld mit hoher Auflösung geschätzt. Wir demonstrieren damit eine parallele Inversion, die auf Linux HPC-Clusters durchgeführt werden kann. Die in diesem realistischen Beispiel auftretenden Strömungsfelder induzieren Transportprobleme, die im konvektionsdominanten Fall mit der Methode der Stromlinien-Diffusion nicht mehr lösbar sind.

Document type: Dissertation
Supervisor: Bastian, Prof. Dr. Peter
Date of thesis defense: 13 February 2015
Date Deposited: 10 Mar 2015 09:45
Date: 2015
Faculties / Institutes: The Faculty of Mathematics and Computer Science > Dean's Office of The Faculty of Mathematics and Computer Science
DDC-classification: 004 Data processing Computer science
510 Mathematics
550 Earth sciences
Controlled Keywords: convectiondominated, downstream, adaptive, Gauss-Newton, adjoint, FFT, Bayes, DG, SUPG, FEM
About | FAQ | Contact | Imprint |
OA-LogoDINI certificate 2013Logo der Open-Archives-Initiative