A hydrological modelling effort using MODFLOW was undertaken in order to determine the relative importance of some of the factors influencing hydraulic gradient reversals in peatlands. Model domains were of two types, large raised bog type (LRBT) and kettle bog type (KBT), and were made to undergo various levels of meteorological forcing (water deficit). Substrate, too, was varied in order to determine its importance on reversals. Domain-wide reversals were successfully simulated in LRBT systems, but not in KBT systems. Although simulated flow patterns matched field-observed patterns, both pre- and post-drought, simulated reversals occurred more quickly than in the field. This may be due to insufficiently distributed parameters, such as hydraulic conductivity. Reversals were easily terminated by simulating non-drought conditions. In the LRBT system, reversal duration decreased, and time-to-reversal increased, with a decrease in drought severity. Increasing drought severity in KBT systems had the opposite effect on the duration of semi-reversed flow patterns, suggesting a possibly different/additional mechanism for flow reversals in KBT systems. Hydraulic conductivity had an appreciable effect on flow reversal evolution, though neither changing porosity, nor differences in catotelm layering had a great effect.

No description available.

This research proposes an Infrastructure to model complex systems for hydrological modeling. Currently, the three main hydrological packages are: i) SEAWAT (modeling groundwater flow); ii) HECRAS (modeling surface water flow); iii) HEC-HMS (modeling atmospheric water flow). Each of these models is self-contained and has a different timescale and simulation speed. Consequently, any integrated model will only run as fast as the slowest of the models. This makes it difficult to provide reliable and dynamic information on water levels and water availability for a given geographical region in a timely manner. The goal of this research is to facilitate the integration of multiple hydrological models from different hydrological packages by applying Electronic Design Automation (EDA) methodologies, including System Level Design (SLD) methodology, SystemC-AMS language, Python language and libraries (numpy, Statsmodels, and ctypes). The EDA methodology brings in the additional advantage of significantly improved simulation speed. The Infrastructure to Model Complex Systems applications is demonstrated using the following SEAWAT benchmark problems: i) Case 1; ii) Henry; iii) Elder problem. Simulation results from the aforementioned benchmarks are analyzed and discussed. Lastly, future research work is presented. / Includes bibliography. / Dissertation (Ph.D.)--Florida Atlantic University, 2014. / FAU Electronic Theses and Dissertations Collection

Contamination of ground water by non-aqueous phase liquids (NAPLs) has received increasing attention. The most common approach to numerical modeling of NAPL movement through the unsaturated zone is the use of the finite difference or finite element methods to solve a set of partial differential equations derived from Darcy's law and the continuity equations (Abriola and Pinder, 1985; Kaluarachchi and Parker, 1989). These methods work well in many settings, but have given little insights as to why certain non-ideal flow phenomena will occur. The network modeling method, which considers flow at the pore-scale, was used in this study to better understand macroscopic flow phenomena in porous media. Pore-scale network models approximate porous medium as a connected network of pores and channels. Two and three-dimensional pore-scale network models were constructed in this study. A uniform statistical distribution was assumed to represent the random arrangement of pore and tube sizes. Both hysteristic scanning curves and intermediate fluid distribution are studied. The simulation results suggested that network models may be used to predict the characteristic curves of three-phase systems. The results also suggested that three-dimensional models are necessary to study the three-phase problems. Two-dimensional models do not provide realistic results as evidenced by their inability to provide scale-invariant representation of flow processes. The network sizes used in this study ranged from 10 x 5 (50) to 156 x 78 (12168) pores for two-dimensional and from 10 x 5 x 5 (250) to 100 x 50 x 5 (25000) pores for three-dimensional domains. The domain size of 100 x 50 x 5 pores was large enough to provide descriptions independent of the domain scale. The one important limitation of network models is the considerable computational requirements. The use of very high speed computers is essential. Except for this limitation, the network model provides an invaluable technique to study fluid transport mechanisms in the vadose zone. / Graduation date: 1996

"January 14, 2004" Includes bibliographical references (leaves 201-213)

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. / Thesis (Ph.D.)--School of Mathematical Sciences (Applied Mathematics), 2004.