Please use this identifier to cite or link to this item: https://hdl.handle.net/2440/37924
Type: Thesis
Title: Numerical Accuracy of Variable-Density Groundwater Flow and Solute Transport Simulations
Author: Woods, Juliette
Issue Date: 2004
School/Discipline: School of Mathematical Sciences : Applied Mathematics
Abstract: The movement of a fluid and solute through a porous medium is of great practical interest because this describes the spread of contaminants through an aquifer. Many contaminants occur at concentrations sufficient to alter the density of the fluid, in which case the physics is typically modelled mathematically by a pair of coupled, nonlinear partial differential equations. There is disagreement as to the exact form of these governing equations. Codes aiming to solve some version of the governing equations are typically tested against the Henry and Elder benchmark problems. Neither benchmark has an analytic solution, so in practice they are treated as exercises in inter code comparison. Different code developers define the boundary conditions of the Henry problem differently, and the Elder problems results are poorly understood. The Henry, Elder and some other problems are simulated on several different codes, which produce widely-varying results. The existing benchmarks are unable to distinguish which code, if any, simulates the problems correctly, illustrating the benchmarks' limitations. To determine whether these discrepancies might be due to numerical error, one popular code, SUTRA, is considered in detail. A numerical analysis of a special case reveals that SUTRA is numerically dispersive. This is confirmed using the Gauss pulse test, a benchmark that does have an analytic solution. To further explain inter code discrepancies, a testcode is developed which allows a choice of numerical methods. Some of the methods are based on SUTRA's while others are finite difference methods of varying levels of accuracy. Simulations of the Elder problem reveal that the benchmark is extremely sensitive to the choice of solution method: qualitative differences are seen in the flow patterns. Finally, the impact of numerical error on a real-world application, the simulation of saline disposals, is considered. Saline disposal basins are used to store saline water away from rivers and agricultural land in parts of Australia. Existing models of disposal basins are assessed in terms of their resemblance to real fieldsite conditions, and in terms of numerical error. This leads to the development of a new model which aims to combine verisimilitude with numerical accuracy.
Advisor: Teubner, Michael David
Dissertation Note: Thesis (Ph.D.)--School of Mathematical Sciences (Applied Mathematics), 2004.
Keywords: SUTRA (Computer file), Groundwater flow Computer simulation, Water salinization Measurement, Groundwater Mathematical models, Finite element method Computer programs, Computer interfaces
Appears in Collections:Research Theses

Files in This Item:
File Description SizeFormat 
01front.pdf89.29 kBAdobe PDFView/Open
02chapter1-chapter5.pdf2.48 MBAdobe PDFView/Open
03chapter6-bibliography.pdf2.32 MBAdobe PDFView/Open


Items in DSpace are protected by copyright, with all rights reserved, unless otherwise indicated.