Parallel finite element algorithm for transient flow in bounded randomly heterogeneous domains

Persistent Link:
http://hdl.handle.net/10150/280226
Title:
Parallel finite element algorithm for transient flow in bounded randomly heterogeneous domains
Author:
Ye, Ming
Issue Date:
2002
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 consider the effect of randomness of hydraulic conductivities K(x) on numerical predictions, without resorting to either Monte Carlo simulation, of transient flow in bounded domains driven by random source, initial and boundary terms. Our aim is to allow optimum unbiased prediction of hydraulic heads h(x, t) and fluxes q(x,t) by means of their respective ensemble moments, <h( x,t)>c and < q(x,t)>c, conditioned on measurements of K(x). These predictors have been shown by Tartakovsky and Neuman (1998) to satisfy exactly a space-time nonlocal (integro-differential) conditional mean flow equation in which < q(x,t)>c is generally non-Darcian. Exact nonlocal equations have been obtained for second conditional moments of head and flux that serve as measures of predictive uncertainty. The authors developed recursive closure approximations for the first and second conditional moment equations through expansion in powers of a small parameter σᵧ , which represents the standard estimation error of ln K(x). The authors explored the possibility of localizing the exact moment equations in real, Laplace- and/or infinite Fourier-transformed domains. In this paper we show how to solve recursive closure approximations of nonlocal first and second conditional moment equations numerically, to first order in σ²ᵧ, in a bounded two-dimensional domain. Our solution is based on Laplace transformation of the moment equations, parallel finite element solution in the complex Laplace domain, and numerical inversion of the solution from the Laplace to the real time domain. We present a detailed comparison between numerical solutions of nonlocal and localized moment equations, and Monte Carlo simulations, under superimposed mean-uniform and convergent flow regimes in two dimensions. The results are shown to compare very well for variances σ²ᵧ as large as 4. The degree to which parallelization enhances computational efficiency is explored.
Type:
text; Dissertation-Reproduction (electronic)
Keywords:
Hydrology.; Engineering, Civil.
Degree Name:
Ph.D.
Degree Level:
doctoral
Degree Program:
Graduate College; Hydrology and Water Resources
Degree Grantor:
University of Arizona
Advisor:
Neuman, Shlomo P.

Full metadata record

DC FieldValue Language
dc.language.isoen_USen_US
dc.titleParallel finite element algorithm for transient flow in bounded randomly heterogeneous domainsen_US
dc.creatorYe, Mingen_US
dc.contributor.authorYe, Mingen_US
dc.date.issued2002en_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 consider the effect of randomness of hydraulic conductivities K(x) on numerical predictions, without resorting to either Monte Carlo simulation, of transient flow in bounded domains driven by random source, initial and boundary terms. Our aim is to allow optimum unbiased prediction of hydraulic heads h(x, t) and fluxes q(x,t) by means of their respective ensemble moments, <h( x,t)>c and < q(x,t)>c, conditioned on measurements of K(x). These predictors have been shown by Tartakovsky and Neuman (1998) to satisfy exactly a space-time nonlocal (integro-differential) conditional mean flow equation in which < q(x,t)>c is generally non-Darcian. Exact nonlocal equations have been obtained for second conditional moments of head and flux that serve as measures of predictive uncertainty. The authors developed recursive closure approximations for the first and second conditional moment equations through expansion in powers of a small parameter σᵧ , which represents the standard estimation error of ln K(x). The authors explored the possibility of localizing the exact moment equations in real, Laplace- and/or infinite Fourier-transformed domains. In this paper we show how to solve recursive closure approximations of nonlocal first and second conditional moment equations numerically, to first order in σ²ᵧ, in a bounded two-dimensional domain. Our solution is based on Laplace transformation of the moment equations, parallel finite element solution in the complex Laplace domain, and numerical inversion of the solution from the Laplace to the real time domain. We present a detailed comparison between numerical solutions of nonlocal and localized moment equations, and Monte Carlo simulations, under superimposed mean-uniform and convergent flow regimes in two dimensions. The results are shown to compare very well for variances σ²ᵧ as large as 4. The degree to which parallelization enhances computational efficiency is explored.en_US
dc.typetexten_US
dc.typeDissertation-Reproduction (electronic)en_US
dc.subjectHydrology.en_US
dc.subjectEngineering, Civil.en_US
thesis.degree.namePh.D.en_US
thesis.degree.leveldoctoralen_US
thesis.degree.disciplineGraduate Collegeen_US
thesis.degree.disciplineHydrology and Water Resourcesen_US
thesis.degree.grantorUniversity of Arizonaen_US
dc.contributor.advisorNeuman, Shlomo P.en_US
dc.identifier.proquest3073277en_US
dc.identifier.bibrecord.b43479224en_US
All Items in UA Campus Repository are protected by copyright, with all rights reserved, unless otherwise indicated.