Aquifer Modeling by Numerical Methods Applied to an Arizona Groundwater Basin
KeywordsAquifers -- Arizona -- Mathematical models.
Groundwater -- Arizona -- Mathematical models.
Aquifers -- Mathematical models.
Groundwater -- Mathematical models.
MetadataShow full item record
RightsCopyright © Arizona Board of Regents
Collection InformationThis title from the Hydrology & Water Resources Technical Reports collection is made available by the Department of Hydrology & Atmospheric Sciences and the University Libraries, University of Arizona. If you have questions about titles in this collection, please contact firstname.lastname@example.org.
AbstractFLUMP, a recently developed mixed explicit -implicit finite -element program, was calibrated against a data base obtained from a portion of the Tucson Basin aquifer, Arizona, and represents its first application to a real -world problem. Two previous models for the same region were constructed (an electric analog and a finite -difference model) in which calibration was based on prescribed flux boundary conditions along stream courses and mountain fronts. These fluxes are not directly measured and estimates are subject to large uncertainties. In contrast, boundary conditions used in the calibration of FLUMP were prescribed hydraulic heads obtained from direct measurement. At prescribed head boundaries FLUMP computed time - varying fluxes representing subsurface lateral flow and recharge along streams. FLUMP correctly calculated fluctuations in recharge along the Santa Cruz River due to fluctuations in storm runoff and sewage effluent release rates. FLUMP also provided valuable insight into distributions of recharge, discharge, and subsurface flow in the study area.Properties of FLUMP were compared with those of two other programs in current use: ISOQUAD, a finite -element program developed by Pinder and Frind (1972), and a finite- difference program developed by the U.S. Geological Survey (Trescott, et al., 1976). It appears that FLUMP can handle a larger class of problems than the other two programs, including those in which the boundary conditions and aquifer parameters vary arbitrarily with time and /or head. FLUMP also has the ability to solve explicitly when accuracy requires small time steps, the ability to solve explicitely in certain parts of the flow region while solving implicitly in other parts, flexibility in mesh design and numbering of nodes, computation of internal as well as external fluxes, and global as well as local mass balance checks at each time step.
Series/Report no.Technical Reports on Natural Resource Systems, No. 32
Showing items related by title, author, creator and subject.
Modeling System Reliability For Digital Preservation: Model Modification and Four-Copy Model StudyHan, Yan; Chan, Chi Pak (British Library, 2008)Research has been studied to evaluate the reliability of storage media and the reliability of a computer backup system. In this paper, we use the Continuous Time Markov Chain to model and analyze the reliability of a computer backup system. We propose a modified model from that of the Constantopoulos, Doerr and Petraki . We analyze the difference, show computational results, and propose new input parameters (e.g. time to repair) for the model from our experience. Further we developed a four-copy data model to test if it fulfills the sample reliability rate set by the RLG-NARA. The modeling process can be applied to construct models for computer preservation systems using different storage media. The reliability of constructed models can be calculated so that preservation institutions can have quantitative data to decide their preservation strategies.
The effects of surface roughness and stress on lattice gas models using kinetic Monte Carlo modelingJackson, Kenneth A.; Bentz, Daniel N. (The University of Arizona., 2003)This thesis reports kinetic Monte Carlo computer simulations using a lattice gas model conducted on a variety of systems. These studies may be divided into two main categories: rod eutectics and related surface roughening, and surface morphology changes due to local stresses. The first grouping is a study of irregular rod eutectic systems. Simulations of directional solidification of rod eutectic systems were conducted using a model similar to the spin one Ising model. Growth of the rods was initiated from columns of pure A atoms embedded in a matrix of B atoms. The growth characteristics of the eutectic depend on the location of the surface roughening transition for the two phases. The surface roughening transition was determined using fluctuation dissipation theory, where the kinetic behavior of the interface is related to a characteristic time of fluctuations about an equilibrium position. These times were determined by time correlations. Results show a sharp transition in the kinetic behavior of the interface as a function of Jackson's alpha factor. This is the first time this method has been used to locate the surface roughening transition. An applied temperature gradient supplied the restoring force for the interface. The roughening behavior of binary alloys was also examined and compared to pure component systems. The second set of simulations reported here examine the effects of local stress on surface morphology. The weakening of bonds due to the dislocation stress field was studied as the origin of the formation of etch-pits at dislocations. Atoms from a diamond cubic lattice were irreversibly removed with a probability which depends on an local surface configuration as well as on the local stress developed from its physical location with respect to a dislocation in the lattice. In accordance with experimental observations, both faceted and non-faceted dislocation etch-pits have been observed. Simulations of crystal surfaces near equilibrium have reproduced direct experimental results using atomic force microscopy (AFM). The probe tip interacts with the shape and the motion of step edges, and the motion of a step is retarded in the vicinity of the tip.
Discrete stochastic arbitrage models for pricing options on short and long term yields: Tests of the Ho and Lee and the Black, Derman and Toy models.Raj, Mahendra.; Sternberg, Joel S.; Dyl, Edward A.; Kwan, Simon H. (The University of Arizona., 1992)The Chicago Board of Options Exchange introduced the short-term and the long-term options on interest rates in June, 1989. This paper develops versions of the Ho and Lee and the Black, Derman and Toy models, two of the most popular arbitrage based discrete stochastic term structure models, and conducts joint tests of the models and the option markets. Both parametric and non-parametric tests are employed. The data consist of the prices of these options and the term structure of interest rates from June, 1989 to June, 1992. The prices of the short and long-term interest options have been obtained from the Wall Street Journals for the period. The term structure of interest rates has been obtained as the prices of treasury bills and the seven, ten and thirty year treasury strips. It is found that the version of Black, Derman and Toy model used in the study is superior to the Ho and Lee model. The Black, Derman and Toy explains the option prices to some extent as seen from the significant R-square values but it does tend to underprice both the short and long-term options slightly. However, since the tests are joint tests of the model and the short and long-term interest rate option markets, it may be that though the model is pricing the options correctly, the markets are not efficient. This may be especially true in view of the fact that the trading in these options has been, at times, quite thin. The Ho and Lee model fails to significantly explain both the short and long-term option prices. The call prices are found to be reasonably close to the actual prices whereas the put prices are grossly over priced. As indicated by the theoretical concepts the unstable and sensitive nature of δ and the tendency for the interest rates at the extreme nodes to assume unrealistic values causes this problem.