Link to the University of Pittsburgh Homepage
Link to the University Library System Homepage Link to the Contact Us Form

Multiblock modeling of flow in porous media and applications

Pencheva, Gergina (2007) Multiblock modeling of flow in porous media and applications. Doctoral Dissertation, University of Pittsburgh. (Unpublished)

Primary Text

Download (6MB) | Preview


We investigate modeling flow in porous media in multiblock domain. Mixed finite element methods are used for subdomain discretizations. Physically meaningful boundary conditions are imposed on the non-matching interfaces via mortar finite element spaces. We investigate the pollution effect of nonmatching grids error on the numerical solution away from interfaces. We prove that most of the error in the velocity occurs along the interfaces, and that high accuracy is preserved in the interior of the subdomains. In case of discontinuous coefficients, the pollution from the singularity affects the accuracy in the whole domain. </br></br> We investigate the upscaling error resulting when fine resolution data is approximated on a very coarse scale. Extending work of Wheeler and Yotov, we incorporate this upscaling error in an a posteriori error estimator for the pressure, velocity and mortar pressure. We employ a non-overlapping domain decomposition method reducing the global system to one that is solved iteratively via a preconditioned conjugate gradient method. This approach is suitable for parallel implementation. The balancing domain decomposition method for mixed finite elements following Cowsar, Mandel, and Wheeler is extended to the case of mortar mixed finite elements on non-matching multiblock grids. The algorithm involves solution of a mortar interface problem with one local Dirichlet solve and one local Neumann solve on each iteration. A coarse solve is used to guarantee consistency and to provide global exchange of information. Quasi-optimal condition number bounds independent of the jump in coefficients are derived. </br></br> We finally consider multiscale mortar mixed finite element discretizations for single and two phase flows. We show optimal convergence and some superconvergence in the fine scale for the solution and its flux. We also derive efficient and reliable a posteriori error estimators suitable for adaptive mesh refinement. </br></br> We have incorporated the above methods into a parallel multiblock simulator on unstructured prismatic meshes employing a non-overlapping domain decomposition algorithm and mortar spaces. Numerical experiments are presented confirming all theoretical results.


Social Networking:
Share |


Item Type: University of Pittsburgh ETD
Status: Unpublished
CreatorsEmailPitt UsernameORCID
ETD Committee:
TitleMemberEmail AddressPitt UsernameORCID
Committee ChairYotov, Ivanyotov@math.pitt.eduYOTOV
Committee MemberSussman, Myronsussmanm@math.pitt.eduMMS125
Committee MemberRiviere,
Committee MemberLayton, William Jwjl@pitt.eduWJL
Date: 22 June 2007
Date Type: Completion
Defense Date: 4 January 2007
Approval Date: 22 June 2007
Submission Date: 18 December 2006
Access Restriction: No restriction; Release the ETD for access worldwide immediately.
Institution: University of Pittsburgh
Schools and Programs: Dietrich School of Arts and Sciences > Mathematics
Degree: PhD - Doctor of Philosophy
Thesis Type: Doctoral Dissertation
Refereed: Yes
Uncontrolled Keywords: ; balancing domain decomposition
Other ID:, etd-12182006-141339
Date Deposited: 10 Nov 2011 20:11
Last Modified: 15 Nov 2016 13:54


Monthly Views for the past 3 years

Plum Analytics

Actions (login required)

View Item View Item