{"id":{"repo_id":"uiuc","oai_identifier":"oai:www.ideals.illinois.edu:2142/97245"},"canonical_url":"https://search.dev.ndltd.org/etd/uiuc/oai:www.ideals.illinois.edu:2142/97245","repository":{"repo_id":"uiuc","name":"University of Illinois - Urbana-Champaign","base_url":"https://www.ideals.illinois.edu/oai-pmh"},"display":{"title":"Development of a high-order accurate reservoir simulator using spectral element method","abstract":"Reservoir simulation serves as an important tool for reservoir management to predict and optimize the future performance of a reservoir. Modeling multiphase fluid flow in porous media can be computationally challenging due to heterogeneous geologic features and sharp, moving fluid interfaces. Low-order convergence of conventional reservoir simulators, often based on finite volume or finite element methods, can be computationally prohibitive to fully resolve the physics of flows with sharp fronts. A high-order numerical scheme with strong stability properties and robustness to accurately resolve the physics of reservoir flows is thus highly desirable. The current research work focuses on the case of two-phase immiscible, incompressible fluid flow in oil reservoirs and seeks to develop a high-order accurate numerical scheme. Governing equations for two-phase incompressible flow in porous media, derived in terms of pressure p and saturation s, are a coupled system of partial differential equations (PDEs) with elliptic-parabolic nature in general. Under certain conditions, the governing equations can become elliptic-hyperbolic in nature. Computational challenges in numerically solving PDEs with discontinuous solutions or discontinuous data are identified and various numerical schemes reviewed for meeting desired goals of accuracy and stability. In particular, numerical comparisons are given for convective hyperbolic PDEs using (a) a first-order upwind finite volume method, (b) a hybrid 1st/2nd order Jameson-Schmidt-Turkel scheme, (c) explicit Runge-Kutta Discontinuous Galerkin method of Cockburn and Shu, and (d) an implicit Discontinuous Galerkin method. Based on the conclusions drawn from these schemes, a fully coupled implicit Discontinuous Galerkin spectral element method (DGSEM) is proposed to solve the two non-linear governing equations in pressure and saturation. The proposed DGSEM scheme is capable of using high-order spectral elements and controlled amount of artificial diffusion to achieve stable, robust numerical solutions. Spectral/hp-convergence of DGSEM scheme is demonstrated for quarter five-spot pattern flow in homogeneous reservoirs. Superior performance capabilities of the DGSEM scheme in resolving heterogeneous geologic features are demonstrated. A novel approach for using volumetric source terms for injection/production wells on spectral element grids not required to be conforming with geometric or geologic features is presented. Numerical results are presented for simulations with different types of geologic heterogeneities and injection/production patterns for 2D reservoirs.","abstract_html":"Reservoir simulation serves as an important tool for reservoir management to predict and optimize the future performance of a reservoir. Modeling multiphase fluid flow in porous media can be computationally challenging due to heterogeneous geologic features and sharp, moving fluid interfaces. Low-order convergence of conventional reservoir simulators, often based on finite volume or finite element methods, can be computationally prohibitive to fully resolve the physics of flows with sharp fronts. A high-order numerical scheme with strong stability properties and robustness to accurately resolve the physics of reservoir flows is thus highly desirable. The current research work focuses on the case of two-phase immiscible, incompressible fluid flow in oil reservoirs and seeks to develop a high-order accurate numerical scheme. Governing equations for two-phase incompressible flow in porous media, derived in terms of pressure p and saturation s, are a coupled system of partial differential equations (PDEs) with elliptic-parabolic nature in general. Under certain conditions, the governing equations can become elliptic-hyperbolic in nature. Computational challenges in numerically solving PDEs with discontinuous solutions or discontinuous data are identified and various numerical schemes reviewed for meeting desired goals of accuracy and stability. In particular, numerical comparisons are given for convective hyperbolic PDEs using (a) a first-order upwind finite volume method, (b) a hybrid 1st/2nd order Jameson-Schmidt-Turkel scheme, (c) explicit Runge-Kutta Discontinuous Galerkin method of Cockburn and Shu, and (d) an implicit Discontinuous Galerkin method. Based on the conclusions drawn from these schemes, a fully coupled implicit Discontinuous Galerkin spectral element method (DGSEM) is proposed to solve the two non-linear governing equations in pressure and saturation. The proposed DGSEM scheme is capable of using high-order spectral elements and controlled amount of artificial diffusion to achieve stable, robust numerical solutions. Spectral/hp-convergence of DGSEM scheme is demonstrated for quarter five-spot pattern flow in homogeneous reservoirs. Superior performance capabilities of the DGSEM scheme in resolving heterogeneous geologic features are demonstrated. A novel approach for using volumetric source terms for injection/production wells on spectral element grids not required to be conforming with geometric or geologic features is presented. Numerical results are presented for simulations with different types of geologic heterogeneities and injection/production patterns for 2D reservoirs.","abstract_has_math":false,"creators":["Taneja, Ankur"],"institution":"University of Illinois at Urbana-Champaign","degree_name":"Ph.D.","degree_level":"Dissertation","degree_discipline":"Chemical Engineering","degree_department":null,"school":null,"contributors":["Higdon, Jonathan","Sing, Charles","Fischer, Paul","Tortorelli, Daniel"],"advisors":[],"committee_chairs":[],"committee_members":[],"year":2017,"date_issued":"2017-08-10T19:14:21Z","date_published":"2017-08-10T19:14:21Z","updated_at":"2026-07-22T22:24:32Z","subjects":["Discontinuous galerkin method","Spectral element method","Reservoir simulation","Reservoir management"],"languages":["en"],"rights":["Copyright 2017 Ankur Taneja"],"rights_urls":[],"identifier_entries":[]},"links":{"outbound_url":"http://hdl.handle.net/2142/97245","outbound_label":"Handle","outbound_source":"dc:identifier"},"metadata_groups":[{"id":"people","label":"People","entries":[{"key":"dc:contributor","label":"Contributor","values":["Higdon, Jonathan","Sing, Charles","Fischer, Paul","Tortorelli, Daniel"]},{"key":"dc:creator","label":"Author","values":["Taneja, Ankur"]}]},{"id":"academic_context","label":"Academic Context","entries":[{"key":"dc:date","label":"Dc Date","values":["2017-08-10T19:14:21Z","2017-01-09","2017-05"]},{"key":"dc:type","label":"Dc Type","values":["text"]},{"key":"thesis:degree_discipline","label":"Discipline","values":["Chemical Engineering"]},{"key":"thesis:degree_level","label":"Degree Level","values":["Dissertation"]},{"key":"thesis:degree_name","label":"Degree Name","values":["Ph.D."]},{"key":"thesis:institution_name","label":"Thesis Institution Name","values":["University of Illinois at Urbana-Champaign"]}]},{"id":"subjects_keywords","label":"Subjects and Keywords","entries":[{"key":"dc:subject","label":"Dc Subject","values":["Discontinuous galerkin method","Spectral element method","Reservoir simulation","Reservoir management"]}]},{"id":"language_rights","label":"Language and Rights","entries":[{"key":"dc:language","label":"Dc Language","values":["en"]},{"key":"dc:rights","label":"Dc Rights","values":["Copyright 2017 Ankur Taneja"]}]},{"id":"identifiers","label":"Identifiers","entries":[{"key":"dc:identifier","label":"Identifier","values":["http://hdl.handle.net/2142/97245"]}]},{"id":"additional","label":"Additional Metadata","entries":[{"key":"dc:description","label":"Description","values":["Reservoir simulation serves as an important tool for reservoir management to predict and optimize the future performance of a reservoir. Modeling multiphase fluid flow in porous media can be computationally challenging due to heterogeneous geologic features and sharp, moving fluid interfaces. Low-order convergence of conventional reservoir simulators, often based on finite volume or finite element methods, can be computationally prohibitive to fully resolve the physics of flows with sharp fronts. A high-order numerical scheme with strong stability properties and robustness to accurately resolve the physics of reservoir flows is thus highly desirable. The current research work focuses on the case of two-phase immiscible, incompressible fluid flow in oil reservoirs and seeks to develop a high-order accurate numerical scheme. Governing equations for two-phase incompressible flow in porous media, derived in terms of pressure p and saturation s, are a coupled system of partial differential equations (PDEs) with elliptic-parabolic nature in general. Under certain conditions, the governing equations can become elliptic-hyperbolic in nature. Computational challenges in numerically solving PDEs with discontinuous solutions or discontinuous data are identified and various numerical schemes reviewed for meeting desired goals of accuracy and stability. In particular, numerical comparisons are given for convective hyperbolic PDEs using (a) a first-order upwind finite volume method, (b) a hybrid 1st/2nd order Jameson-Schmidt-Turkel scheme, (c) explicit Runge-Kutta Discontinuous Galerkin method of Cockburn and Shu, and (d) an implicit Discontinuous Galerkin method. Based on the conclusions drawn from these schemes, a fully coupled implicit Discontinuous Galerkin spectral element method (DGSEM) is proposed to solve the two non-linear governing equations in pressure and saturation. The proposed DGSEM scheme is capable of using high-order spectral elements and controlled amount of artificial diffusion to achieve stable, robust numerical solutions. Spectral/hp-convergence of DGSEM scheme is demonstrated for quarter five-spot pattern flow in homogeneous reservoirs. Superior performance capabilities of the DGSEM scheme in resolving heterogeneous geologic features are demonstrated. A novel approach for using volumetric source terms for injection/production wells on spectral element grids not required to be conforming with geometric or geologic features is presented. Numerical results are presented for simulations with different types of geologic heterogeneities and injection/production patterns for 2D reservoirs.","Submission original under an indefinite embargo labeled 'Open Access'. The submission was exported from vireo on 2017-08-10 without embargo terms","The student, Ankur Taneja, accepted the attached license on 2017-01-05 at 15:45.","The student, Ankur Taneja, submitted this Dissertation for approval on 2017-01-05 at 15:57.","This Dissertation was approved for publication on 2017-01-09 at 12:45.","DSpace SAF Submission Ingestion Package generated from Vireo submission #10530 on 2017-08-10 at 13:37:09","Made available in DSpace on 2017-08-10T19:14:21Z (GMT). No. of bitstreams: 16 TANEJA-DISSERTATION-2017.pdf: 5082263 bytes, checksum: 2d9a460dd00da5e3d2311fc22b1c2bd2 (MD5) fig_5-10a.mp4: 911371 bytes, checksum: 2b360cd42aa7c910884e1ddfd70df8d2 (MD5) fig_5-10b.mp4: 932161 bytes, checksum: f3abc8900add6207f1f2392f63ecec6d (MD5) fig_5-10c.mp4: 907018 bytes, checksum: 941653eb1e37814a08ed580ec854b202 (MD5) fig_5-15.mp4: 962290 bytes, checksum: dfe00f97ca6d56cc67a7498067baf367 (MD5) fig_5-16.mp4: 1023313 bytes, checksum: 1841c6f34eeca9b63cc0682321f8e054 (MD5) fig_5-17.mp4: 1035174 bytes, checksum: 49a36fb4fa89d52c1566d681ebfa0ac6 (MD5) fig_5-18.mp4: 1050292 bytes, checksum: f3852e4aee55c66f318ff2d36a410931 (MD5) fig_5-19.mp4: 986214 bytes, checksum: b0dc9e059633b04a6a6f8e0fe009101d (MD5) fig_5-20.mp4: 995058 bytes, checksum: 3744ee5ca4808d361c2eb692418aec11 (MD5) fig_5-21.mp4: 1017429 bytes, checksum: af7166a79bb4bee365181400c61decef (MD5) fig_5-23a.mp4: 820013 bytes, checksum: 7a16da165ee6a943f86a7fdee687c026 (MD5) fig_5-23b.mp4: 842424 bytes, checksum: f50ef44f3f6810dbf1eda876f8c1cab5 (MD5) thesis_source_files_Jan6_2017.zip: 5578397 bytes, checksum: 7aae054a36e8100afa0cfc7e818bd13e (MD5) LICENSE.txt: 4209 bytes, checksum: d75528ad94658afd92673b992059ebca (MD5) PROQUEST_LICENSE.txt: 4555 bytes, checksum: c51c4d4c7c40669c4b9e15ad8a212a28 (MD5) Previous issue date: 2017-01-09"]},{"key":"dc:format","label":"Dc Format","values":["application/pdf"]},{"key":"dc:title","label":"Title","values":["Development of a high-order accurate reservoir simulator using spectral element method"]}]}],"canonical_facts":{"dc:contributor":["Higdon, Jonathan","Sing, Charles","Fischer, Paul","Tortorelli, Daniel"],"dc:creator":["Taneja, Ankur"],"dc:date":["2017-08-10T19:14:21Z","2017-01-09","2017-05"],"dc:description":["Reservoir simulation serves as an important tool for reservoir management to predict and optimize the future performance of a reservoir. Modeling multiphase fluid flow in porous media can be computationally challenging due to heterogeneous geologic features and sharp, moving fluid interfaces. Low-order convergence of conventional reservoir simulators, often based on finite volume or finite element methods, can be computationally prohibitive to fully resolve the physics of flows with sharp fronts. A high-order numerical scheme with strong stability properties and robustness to accurately resolve the physics of reservoir flows is thus highly desirable. The current research work focuses on the case of two-phase immiscible, incompressible fluid flow in oil reservoirs and seeks to develop a high-order accurate numerical scheme. Governing equations for two-phase incompressible flow in porous media, derived in terms of pressure p and saturation s, are a coupled system of partial differential equations (PDEs) with elliptic-parabolic nature in general. Under certain conditions, the governing equations can become elliptic-hyperbolic in nature. Computational challenges in numerically solving PDEs with discontinuous solutions or discontinuous data are identified and various numerical schemes reviewed for meeting desired goals of accuracy and stability. In particular, numerical comparisons are given for convective hyperbolic PDEs using (a) a first-order upwind finite volume method, (b) a hybrid 1st/2nd order Jameson-Schmidt-Turkel scheme, (c) explicit Runge-Kutta Discontinuous Galerkin method of Cockburn and Shu, and (d) an implicit Discontinuous Galerkin method. Based on the conclusions drawn from these schemes, a fully coupled implicit Discontinuous Galerkin spectral element method (DGSEM) is proposed to solve the two non-linear governing equations in pressure and saturation. The proposed DGSEM scheme is capable of using high-order spectral elements and controlled amount of artificial diffusion to achieve stable, robust numerical solutions. Spectral/hp-convergence of DGSEM scheme is demonstrated for quarter five-spot pattern flow in homogeneous reservoirs. Superior performance capabilities of the DGSEM scheme in resolving heterogeneous geologic features are demonstrated. A novel approach for using volumetric source terms for injection/production wells on spectral element grids not required to be conforming with geometric or geologic features is presented. Numerical results are presented for simulations with different types of geologic heterogeneities and injection/production patterns for 2D reservoirs.","Submission original under an indefinite embargo labeled 'Open Access'. The submission was exported from vireo on 2017-08-10 without embargo terms","The student, Ankur Taneja, accepted the attached license on 2017-01-05 at 15:45.","The student, Ankur Taneja, submitted this Dissertation for approval on 2017-01-05 at 15:57.","This Dissertation was approved for publication on 2017-01-09 at 12:45.","DSpace SAF Submission Ingestion Package generated from Vireo submission #10530 on 2017-08-10 at 13:37:09","Made available in DSpace on 2017-08-10T19:14:21Z (GMT). No. of bitstreams: 16 TANEJA-DISSERTATION-2017.pdf: 5082263 bytes, checksum: 2d9a460dd00da5e3d2311fc22b1c2bd2 (MD5) fig_5-10a.mp4: 911371 bytes, checksum: 2b360cd42aa7c910884e1ddfd70df8d2 (MD5) fig_5-10b.mp4: 932161 bytes, checksum: f3abc8900add6207f1f2392f63ecec6d (MD5) fig_5-10c.mp4: 907018 bytes, checksum: 941653eb1e37814a08ed580ec854b202 (MD5) fig_5-15.mp4: 962290 bytes, checksum: dfe00f97ca6d56cc67a7498067baf367 (MD5) fig_5-16.mp4: 1023313 bytes, checksum: 1841c6f34eeca9b63cc0682321f8e054 (MD5) fig_5-17.mp4: 1035174 bytes, checksum: 49a36fb4fa89d52c1566d681ebfa0ac6 (MD5) fig_5-18.mp4: 1050292 bytes, checksum: f3852e4aee55c66f318ff2d36a410931 (MD5) fig_5-19.mp4: 986214 bytes, checksum: b0dc9e059633b04a6a6f8e0fe009101d (MD5) fig_5-20.mp4: 995058 bytes, checksum: 3744ee5ca4808d361c2eb692418aec11 (MD5) fig_5-21.mp4: 1017429 bytes, checksum: af7166a79bb4bee365181400c61decef (MD5) fig_5-23a.mp4: 820013 bytes, checksum: 7a16da165ee6a943f86a7fdee687c026 (MD5) fig_5-23b.mp4: 842424 bytes, checksum: f50ef44f3f6810dbf1eda876f8c1cab5 (MD5) thesis_source_files_Jan6_2017.zip: 5578397 bytes, checksum: 7aae054a36e8100afa0cfc7e818bd13e (MD5) LICENSE.txt: 4209 bytes, checksum: d75528ad94658afd92673b992059ebca (MD5) PROQUEST_LICENSE.txt: 4555 bytes, checksum: c51c4d4c7c40669c4b9e15ad8a212a28 (MD5) Previous issue date: 2017-01-09"],"dc:format":["application/pdf"],"dc:identifier":["http://hdl.handle.net/2142/97245"],"dc:language":["en"],"dc:rights":["Copyright 2017 Ankur Taneja"],"dc:subject":["Discontinuous galerkin method","Spectral element method","Reservoir simulation","Reservoir management"],"dc:title":["Development of a high-order accurate reservoir simulator using spectral element method"],"dc:type":["text"],"thesis:degree_discipline":["Chemical Engineering"],"thesis:degree_level":["Dissertation"],"thesis:degree_name":["Ph.D."],"thesis:institution_name":["University of Illinois at Urbana-Champaign"]},"updated_at":"2026-07-22T22:24:32Z"}