Techniques for improving numerical modeling of water flow in variably saturated, heterogenous media.

Persistent Link:
http://hdl.handle.net/10150/187427
Title:
Techniques for improving numerical modeling of water flow in variably saturated, heterogenous media.
Author:
Pan, Lehua.
Issue Date:
1995
Publisher:
The University of Arizona.
Rights:
Copyright © is held by the author. Digital access to this material is made possible by the University Libraries, University of Arizona. Further transmission, reproduction or presentation (such as public display or performance) of protected items is prohibited except with permission of the author.
Abstract:
We introduce a simple nonlinear transformation. The new approach was compared to Kirkland's φ-based method and Ross' p-based method as well as Celia's h-based method. All algorithms use the modified Picard iteration with a linear solver of the Cholesky preconditioned conjugate gradient method (2-D) and with Thomas' method (1-D). An adaptively optimum relaxation procedure was also proposed to improve the efficiency and the robustness for 2-D cases. To test the effects of the transformation and the optimal relaxation. a total of 14 different 1-D and 2-D infiltration cases were considered. The results show that the newly introduced P(t)-based method is the most effective of the four methods and the simplest of the three transformation methods. It is several orders of magnitude faster than the h-based method and also reduces the truncation error using the same grid. The proposed adaptively optimum relaxation procedure, enhances the efficiency and robustness of all four methods to different degrees. Combination of the P(t)-based method with the adaptively optimum relaxation procedure results in a very efficient and robust algorithm. Modeling water flow in variably saturated, porous media is important in many branches of science and engineering. Highly nonlinear relationships between water content and hydraulic conductivity and soil-water pressure result in very steep wetting fronts causing numerical problems. These include poor efficiency when modeling water infiltration into very dry porous media, and numerical oscillation near a steep wetting front when a traditional mass-distributed finite element method is used.
Type:
text; Dissertation-Reproduction (electronic)
Degree Name:
Ph.D.
Degree Level:
doctoral
Degree Program:
Soil, Water and Environmental Science; Graduate College
Degree Grantor:
University of Arizona
Committee Chair:
Wierenga, Peter J.

Full metadata record

DC FieldValue Language
dc.language.isoenen_US
dc.titleTechniques for improving numerical modeling of water flow in variably saturated, heterogenous media.en_US
dc.creatorPan, Lehua.en_US
dc.contributor.authorPan, Lehua.en_US
dc.date.issued1995en_US
dc.publisherThe University of Arizona.en_US
dc.rightsCopyright © is held by the author. Digital access to this material is made possible by the University Libraries, University of Arizona. Further transmission, reproduction or presentation (such as public display or performance) of protected items is prohibited except with permission of the author.en_US
dc.description.abstractWe introduce a simple nonlinear transformation. The new approach was compared to Kirkland's φ-based method and Ross' p-based method as well as Celia's h-based method. All algorithms use the modified Picard iteration with a linear solver of the Cholesky preconditioned conjugate gradient method (2-D) and with Thomas' method (1-D). An adaptively optimum relaxation procedure was also proposed to improve the efficiency and the robustness for 2-D cases. To test the effects of the transformation and the optimal relaxation. a total of 14 different 1-D and 2-D infiltration cases were considered. The results show that the newly introduced P(t)-based method is the most effective of the four methods and the simplest of the three transformation methods. It is several orders of magnitude faster than the h-based method and also reduces the truncation error using the same grid. The proposed adaptively optimum relaxation procedure, enhances the efficiency and robustness of all four methods to different degrees. Combination of the P(t)-based method with the adaptively optimum relaxation procedure results in a very efficient and robust algorithm. Modeling water flow in variably saturated, porous media is important in many branches of science and engineering. Highly nonlinear relationships between water content and hydraulic conductivity and soil-water pressure result in very steep wetting fronts causing numerical problems. These include poor efficiency when modeling water infiltration into very dry porous media, and numerical oscillation near a steep wetting front when a traditional mass-distributed finite element method is used.en_US
dc.typetexten_US
dc.typeDissertation-Reproduction (electronic)en_US
thesis.degree.namePh.D.en_US
thesis.degree.leveldoctoralen_US
thesis.degree.disciplineSoil, Water and Environmental Scienceen_US
thesis.degree.disciplineGraduate Collegeen_US
thesis.degree.grantorUniversity of Arizonaen_US
dc.contributor.chairWierenga, Peter J.en_US
dc.contributor.committeememberWarrick, A. W.en_US
dc.contributor.committeememberBrusseau, Marken_US
dc.contributor.committeememberYeh, T.-C. Jimen_US
dc.contributor.committeememberNeuman, Shlomo P.en_US
dc.identifier.proquest9623319en_US
All Items in UA Campus Repository are protected by copyright, with all rights reserved, unless otherwise indicated.