Numerical Methods for Two-phase Flow in Rigid and Deformable Porous Media

dc.contributor.advisorRiviere, Beatrice
dc.creatorShen, Boqian
dc.date.accessioned2022-09-26T19:13:17Z
dc.date.available2022-11-01T05:01:15Z
dc.date.created2022-05
dc.date.issued2022-04-22
dc.date.submittedMay 2022
dc.date.updated2022-09-26T19:13:17Z
dc.description.abstractThe thesis focuses on developing numerical schemes for two-phase flow in rigid and deformable porous media problems. We present a stable and efficient sequential Discontinuous Galerkin (DG) method for solving the linear poroelasticity equations, which characterize two-phase flow within a deformable porous media. More precisely, we approximate the pressure of the wetting phase, the pressure of the non-wetting phase, and the displacement of the solid skeleton in three dimensions by a high-order interior penalty discontinuous Galerkin (IPDG) spatial discretization combined with a backward Euler discretization in time. The proposed work is based on previous developments in single fluid flow in deformable porous media. The numerical scheme solves the coupled equations sequentially while keeping each equation implicitly with respect to its unknown. The equations are fully decoupled in this sequential approach, which significantly reduces the computational cost compared to the implicit and iterative approaches. Numerical experiments show the convergence of the scheme is optimal. Finally, we apply the sequential DG scheme to a variety of physical problems with realistic data including common benchmarks, heterogeneous porous media with discontinuous permeability, porosity and capillary pressure, and porous media subjected to load. The second part of this thesis proposes an adaptive hybrid numerical scheme for solving two-phase flow in rigid porous media problems. The spatial discretization for transport phenomena problem in heterogeneous porous media requires locally mass conservative methods, such as finite volume methods and discontinuous Galerkin methods. The numerical scheme uses discontinuous Galerkin methods in regions of interest where high accuracy is needed and uses finite volume methods in the rest of the domain. The proposed schemes take advantage of the high accuracy of the discontinuous Galerkin method due to its local mesh adaptivity and local choice of polynomial degree. Finite volume methods are only first-order accurate but computationally efficient and robust for general geometries with structured mesh. We develop adaptive indicators to dynamically identify the regions of each method. By using such an adaptive indicator we are able to find the optimal balance between accuracy and computational cost.
dc.embargo.terms2022-11-01
dc.format.mimetypeapplication/pdf
dc.identifier.citationShen, Boqian. "Numerical Methods for Two-phase Flow in Rigid and Deformable Porous Media." (2022) Diss., Rice University. <a href="https://hdl.handle.net/1911/113384">https://hdl.handle.net/1911/113384</a>.
dc.identifier.urihttps://hdl.handle.net/1911/113384
dc.language.isoeng
dc.rightsCopyright is held by the author, unless otherwise indicated. Permission to reuse, publish, or reproduce the work beyond the bounds of fair use or other exemptions to copyright law must be obtained from the copyright holder.
dc.subjecttwo-phase poroelasticity
dc.subjectsequential implicit
dc.subjectdiscontinuous Galerkin
dc.subjectheterogeneities
dc.subjectmultinumerics
dc.subjectFinite Volume
dc.subjectcoupled flow
dc.subjectgeomechanics
dc.titleNumerical Methods for Two-phase Flow in Rigid and Deformable Porous Media
dc.typeThesis
dc.type.materialText
thesis.degree.departmentComputational and Applied Mathematics
thesis.degree.disciplineEngineering
thesis.degree.grantorRice University
thesis.degree.levelDoctoral
thesis.degree.nameDoctor of Philosophy
Files
Original bundle
Now showing 1 - 2 of 2
Loading...
Thumbnail Image
Name:
SHEN-DOCUMENT-2022.pdf
Size:
48.6 MB
Format:
Adobe Portable Document Format
Loading...
Thumbnail Image
Name:
CertificateOfCompletion.pdf
Size:
117.62 KB
Format:
Adobe Portable Document Format
License bundle
Now showing 1 - 2 of 2
No Thumbnail Available
Name:
PROQUEST_LICENSE.txt
Size:
5.84 KB
Format:
Plain Text
Description:
No Thumbnail Available
Name:
LICENSE.txt
Size:
2.6 KB
Format:
Plain Text
Description: