{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2025,10,12]],"date-time":"2025-10-12T03:22:49Z","timestamp":1760239369024,"version":"build-2065373602"},"reference-count":50,"publisher":"MDPI AG","issue":"4","license":[{"start":{"date-parts":[[2020,11,13]],"date-time":"2020-11-13T00:00:00Z","timestamp":1605225600000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":["Computation"],"abstract":"<jats:p>Carbon capture, utilization, and storage (CCUS) describes a set of technically viable processes to separate carbon dioxide (CO2) from industrial byproduct streams and inject it into deep geologic formations for long-term storage. Legacy wells located within the spatial domain of new injection and production activities represent potential pathways for fluids (i.e., CO2 and aqueous phase) to leak through compromised components (e.g., through fractures or micro-annulus pathways). The finite element (FE) method is a well-established numerical approach to simulate the coupling between multi-phase fluid flow and solid phase deformation interactions that occur in a compromised well system. We assumed the spatial domain consists of a three-phases system: a solid, liquid, and gas phase. For flow in the two fluids phases, we considered two sets of primary variables: the first considering capillary pressure and gas pressure (PP) scheme, and the second considering liquid pressure and gas saturation (PS) scheme. Fluid phases were coupled with the solid phase using the full coupling (i.e., monolithic coupling) and iterative coupling (i.e., sequential coupling) approaches. The challenge of achieving numerical stability in the coupled formulation in heterogeneous media was addressed using the mass lumping and the upwinding techniques. Numerical results were compared with three benchmark problems to assess the performance of coupled FE solutions: 1D Terzaghi\u2019s consolidation, Liakopoulos experiments, and the Kueper and Frind experiments. We found good agreement between our results and the three benchmark problems. For the Kueper and Frind test, the PP scheme successfully captured the observed experimental response of the non-aqueous phase infiltration, in contrast to the PS scheme. These exercises demonstrate the importance of fluid phase primary variable selection for heterogeneous porous media. We then applied the developed model to the hypothetical case of leakage along a compromised well representing a heterogeneous media. Considering the mass lumping and the upwinding techniques, both the monotonic and the sequential coupling provided identical results, but mass lumping was needed to avoid numerical instabilities in the sequential coupling. Additionally, in the monolithic coupling, the magnitude of primary variables in the coupled solution without mass lumping and the upwinding is higher, which is essential for the risk-based analyses.<\/jats:p>","DOI":"10.3390\/computation8040098","type":"journal-article","created":{"date-parts":[[2020,11,13]],"date-time":"2020-11-13T08:44:02Z","timestamp":1605257042000},"page":"98","update-policy":"https:\/\/doi.org\/10.3390\/mdpi_crossmark_policy","source":"Crossref","is-referenced-by-count":4,"title":["Effect of Computational Schemes on Coupled Flow and Geo-Mechanical Modeling of CO2 Leakage through a Compromised Well"],"prefix":"10.3390","volume":"8","author":[{"ORCID":"https:\/\/orcid.org\/0000-0002-1798-1162","authenticated-orcid":false,"given":"Mohammad","family":"Islam","sequence":"first","affiliation":[{"name":"National Energy Technology Laboratory (NETL), US Department of Energy, Pittsburgh, PA 15236, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Nicolas","family":"Huerta","sequence":"additional","affiliation":[{"name":"Pacific Northwest National Laboratory, Richland, WA 99352, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Robert","family":"Dilmore","sequence":"additional","affiliation":[{"name":"National Energy Technology Laboratory (NETL), US Department of Energy, Pittsburgh, PA 15236, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"1968","published-online":{"date-parts":[[2020,11,13]]},"reference":[{"key":"ref_1","doi-asserted-by":"crossref","unstructured":"Volk, T. (2008). CO2 Rising: The World\u2019s Greatest Environmental Challenge, MIT Press.","DOI":"10.7551\/mitpress\/7945.001.0001"},{"key":"ref_2","doi-asserted-by":"crossref","unstructured":"Nordbotten, J.M., and Celia, M.A. (2012). Geological Storage of CO2 Modeling Approaches for Large-Scale Simulation, John Wiley and Sons.","DOI":"10.1002\/9781118137086"},{"key":"ref_3","doi-asserted-by":"crossref","first-page":"953","DOI":"10.1016\/S0196-8904(99)00149-1","article-title":"Sequestration of CO2 in geological media: criteria and approach for site selection in response to climate change","volume":"41","author":"Bachu","year":"2000","journal-title":"Energy Convers. Manag."},{"key":"ref_4","doi-asserted-by":"crossref","first-page":"S70","DOI":"10.1016\/j.ijggc.2012.12.028","article-title":"Wellbore integrity and corrosion of carbon steel in CO2 geologic storage environments: A literature review","volume":"16","author":"Choi","year":"2013","journal-title":"Int. J. Greenh. Gas Control"},{"key":"ref_5","doi-asserted-by":"crossref","first-page":"292","DOI":"10.1016\/j.ijggc.2015.06.014","article-title":"Recent advances in risk assessment and risk management of geologic CO2 storage","volume":"40","author":"Pawar","year":"2015","journal-title":"Int. J. Greenh. Gas Control"},{"key":"ref_6","doi-asserted-by":"crossref","first-page":"168","DOI":"10.1016\/j.compgeo.2016.10.001","article-title":"Investigation of wellbore microannulus permeability under stress via experimental wellbore mock-up and finite element modeling","volume":"83","author":"Gomez","year":"2017","journal-title":"Comput. Geotech."},{"key":"ref_7","doi-asserted-by":"crossref","first-page":"13","DOI":"10.2118\/124719-PA","article-title":"How To Prevent Loss of Zonal Isolation Through a Comprehensive Analysis of Microannulus Formation","volume":"26","author":"Bois","year":"2011","journal-title":"SPE Drill. Completion"},{"key":"ref_8","doi-asserted-by":"crossref","first-page":"707","DOI":"10.1007\/s00254-004-1073-5","article-title":"Spatial characterization of the location of potentially leaky wells penetrating a deep saline aquifer in a mature sedimentary basin","volume":"46","author":"Gasda","year":"2004","journal-title":"Environ. Geol."},{"key":"ref_9","unstructured":"Rubin, E.S., Keith, D.W., Gilboy, C.F., Wilson, M., Morris, T., Gale, J., and Thambimuthu, K. (2005). Quantitative estimation of CO2 leakage from geological storage: Analytical models, numerical models, and data needs. Greenhouse Gas Control Technologies 7, Elsevier Science Ltd."},{"key":"ref_10","doi-asserted-by":"crossref","first-page":"38","DOI":"10.1016\/j.petrol.2014.05.024","article-title":"Three-dimensional analysis of cement sheath integrity around Wellbores","volume":"121","author":"Wang","year":"2014","journal-title":"J. Pet. Sci. Eng."},{"key":"ref_11","unstructured":"United States Environmental Protection Agency (2008). Determination of the Mechanical Integrity of Injection Wells: United States Environmental Protection Agency, Region 5: Underground Injection Control (UIC) Branch, Regional Guidance #5."},{"key":"ref_12","doi-asserted-by":"crossref","unstructured":"Cavanagh, P.H., Johnson, C.R., le Roy-Delage, S., DeBruijn, G.G., Cooper, I., Guillot, D.J., Bulte, H., and Dargaud, B. (2007, January 20\u201322). Self-Healing Cement\u2014Novel Technology to Achieve Leak-Free Wells. Proceedings of the SPE\/IADC Drilling Conference, Amsterdam, The Netherlands.","DOI":"10.2523\/105781-MS"},{"key":"ref_13","doi-asserted-by":"crossref","first-page":"107","DOI":"10.1306\/eg.04231414004","article-title":"Seepage pathway assessment for natural gas to shallow groundwater during well stimulation, in production, and after abandonment","volume":"21","author":"Dusseault","year":"2014","journal-title":"Environ. Geosci."},{"key":"ref_14","doi-asserted-by":"crossref","first-page":"525","DOI":"10.1007\/s10706-011-9491-0","article-title":"The Geomechanics of CO2 Storage in Deep Sedimentary Formations","volume":"30","author":"Rutqvist","year":"2012","journal-title":"Geotech. Geol. Eng."},{"key":"ref_15","doi-asserted-by":"crossref","first-page":"200","DOI":"10.1016\/j.ijggc.2014.01.003","article-title":"Numerical modeling of the pumping tests at the Ketzin pilot site for CO2 injection: Model calibration and heterogeneity effects","volume":"22","author":"Chen","year":"2014","journal-title":"Int. J. Greenh. Gas Control"},{"key":"ref_16","first-page":"58","article-title":"Retention and distribution of moisture in soils after infiltration has ceased","volume":"10","author":"Liakopoulos","year":"1965","journal-title":"Hydrol. Sci. J."},{"key":"ref_17","unstructured":"Lewis, R.W., and Schrefler, B.A. (1998). The Finite Element Method in the Static and Dynamic Deformation and Consolidation of Porous Media, John Wiley & Sons. [2nd ed.]."},{"key":"ref_18","doi-asserted-by":"crossref","unstructured":"Helmig, R. (1997). Multiphase Flow and Transport Processes in the Subsurface: A Contribution to the Modeling of Hydrosystems, Springer.","DOI":"10.1007\/978-3-642-60763-9"},{"key":"ref_19","doi-asserted-by":"crossref","first-page":"1049","DOI":"10.1029\/91WR00266","article-title":"Two-phase flow in heterogeneous porous media: 1. Model development","volume":"27","author":"Kueper","year":"1991","journal-title":"Water Resour. Res."},{"key":"ref_20","doi-asserted-by":"crossref","first-page":"461","DOI":"10.1016\/S0309-1708(98)00022-0","article-title":"Practical implementation of the fractional flow approach to multi-phase flow simulation","volume":"22","author":"Binning","year":"1999","journal-title":"Adv. Water Resour."},{"key":"ref_21","doi-asserted-by":"crossref","unstructured":"Kolditz, O., Gorke, U.-J., Shao, H., and Wang, W. (2012). Thermo-Hydro-Mechanical-Chemical Processes in Fractured Porous Media: Benchmarks and Examples, Springer.","DOI":"10.1007\/978-3-642-27177-9"},{"key":"ref_22","unstructured":"Zienkiewicz, O.C., Taylor, R.L., and Zhu, J.Z. (2005). The Finite Element Method: Its Basis and Fundamentals, Elsevier. [6th ed.]."},{"key":"ref_23","doi-asserted-by":"crossref","unstructured":"Chen, Z., Huan, G., and Ma, Y. (2006). Computational Methods for Multiphase Flows in Porous Media, Society for Industrial and Applied Mathematics.","DOI":"10.1137\/1.9780898718942"},{"key":"ref_24","unstructured":"Bacon, D., Zheng, L., Yonkofski, C., Demirkanli, I., Rose, K., Zhou, Q., and Yang, Y.-M. (2018). National Risk Assessment Partnership\u2014Application of Risk Assessment Tools and Methodologies to Synthetic and Field Data (PNNL-SA-137185)."},{"key":"ref_25","unstructured":"Brooks, R.H., and Corey, A.T. (1964). Hydraulic Properties of Porous Media, Hydrology Papers 3, Colorado State University."},{"key":"ref_26","doi-asserted-by":"crossref","unstructured":"Coussy, O. (2004). Poromechanics, John Wiley and Sons.","DOI":"10.1002\/0470092718"},{"key":"ref_27","doi-asserted-by":"crossref","first-page":"589","DOI":"10.1007\/s12665-012-1546-x","article-title":"OpenGeoSys: an open-source initiative for numerical simulation of thermo-hydro-mechanical\/chemical (THM\/C) processes in porous media","volume":"67","author":"Kolditz","year":"2012","journal-title":"Environ. Earth Sci."},{"key":"ref_28","doi-asserted-by":"crossref","first-page":"1309","DOI":"10.1002\/nme.2579","article-title":"Gmsh: A 3-D finite element mesh generator with built-in pre- and post-processing facilities","volume":"79","author":"Geuzaine","year":"2009","journal-title":"Int. J. Numer. Methods Eng."},{"key":"ref_29","doi-asserted-by":"crossref","unstructured":"Bear, J., and Bachmat, Y. (1990). Introduction to Modeling of Transport Phenomena in Porous Media, Kluwer Academic Publishers.","DOI":"10.1007\/978-94-009-1926-6"},{"key":"ref_30","doi-asserted-by":"crossref","unstructured":"De Boer, R. (2005). Trends in Continuum Mechanics in Porous Media, Springer.","DOI":"10.1007\/1-4020-3144-0"},{"key":"ref_31","doi-asserted-by":"crossref","unstructured":"Truesdell, C. (1984). Rational Thermodynamics, Springer. [2nd ed.].","DOI":"10.1007\/978-1-4612-5206-1"},{"key":"ref_32","doi-asserted-by":"crossref","first-page":"369","DOI":"10.1007\/BF00165780","article-title":"Reservoir transport equations by compositional approach","volume":"4","author":"Panday","year":"1989","journal-title":"Transp. Porous Media"},{"key":"ref_33","unstructured":"Aziz, K., and Settari, A. (1979). Petroleum Reservoir Simulation, Applied Science Publishers Ltd."},{"key":"ref_34","unstructured":"Yip, S. (2005). Finite difference, finite element and finite volume methods for partial differential equations. Chapter 8.2. Handbook of Materials Modeling, Springer."},{"key":"ref_35","doi-asserted-by":"crossref","unstructured":"Daehlen, M., and Tveito, A. (1997). Finite Element Methods for Two-Phase Flow in Heterogeneous Porous Media. Chapter 10. Numerical Methods and Software Tools in Industrial Mathematics, Springer.","DOI":"10.1007\/978-1-4612-1984-2"},{"key":"ref_36","doi-asserted-by":"crossref","first-page":"3357","DOI":"10.1029\/94WR02046","article-title":"A comparison of Picard and Newton iteration in the numerical solution of multidimensional variably saturated flow problems","volume":"30","author":"Paniconi","year":"1994","journal-title":"Water Resour. Res."},{"key":"ref_37","doi-asserted-by":"crossref","unstructured":"Saad, Y. (2003). Iterative Methods for Sparse Linear Systems, The Society for Industrial and Applied Mathematics.","DOI":"10.1137\/1.9780898718003"},{"key":"ref_38","unstructured":"Wood, W.L. (1990). Practical Time-Stepping Schemes, Oxford University Press."},{"key":"ref_39","doi-asserted-by":"crossref","unstructured":"Bear, J., and Corapcioglu, M.Y. (1984). The Finite element method in porous media flow. Fundamentals of Transport Phenomena in Porous Media, Martinus Nijhoff Publishers.","DOI":"10.1007\/978-94-009-6175-3"},{"key":"ref_40","doi-asserted-by":"crossref","unstructured":"Terzaghi, K. (1943). Theoretical Soil Mechanics, John Wiley and Sons.","DOI":"10.1002\/9780470172766"},{"key":"ref_41","doi-asserted-by":"crossref","first-page":"359","DOI":"10.1016\/0045-7825(92)90193-N","article-title":"Improved accuracy in finite element analysis of Biot\u2019s consolidation problem","volume":"95","author":"Murad","year":"1992","journal-title":"Comput. Methods Appl. Mech. Eng."},{"key":"ref_42","unstructured":"Liakopoulos, A.C. (1965). Transient Flow through Unsaturated Porous Media, University of California."},{"key":"ref_43","first-page":"318","article-title":"Capillary conduction of liquids through porous mediums","volume":"1","author":"Richards","year":"1931","journal-title":"J. Appl. Phys."},{"key":"ref_44","doi-asserted-by":"crossref","first-page":"3223","DOI":"10.1016\/S0045-7825(00)00390-X","article-title":"A fully coupled dynamic model for two-phase fluid flow in deformable porous media","volume":"190","author":"Schrefler","year":"2001","journal-title":"Comput. Methods Appl. Mech. Eng."},{"key":"ref_45","unstructured":"Ravi, K., Bosma, M., and Gastebled, O. (May, January 30). Safe and Economic Gas Wells through Cement Design for Life of the Well. Proceedings of the SPE Gas Technology Symposium, Calgary, AB, Canada."},{"key":"ref_46","unstructured":"NIST-WebBook (2019). NIST Chemistry WebBook, NIST Standard Reference Database."},{"key":"ref_47","doi-asserted-by":"crossref","unstructured":"Islam, M.N., Bunger, A.P., Huerta, N., and Dilmore, R. (2019). Bentonite Extrusion into Near-Borehole Fracture. Geosciences, 9.","DOI":"10.3390\/geosciences9120495"},{"key":"ref_48","unstructured":"Lubliner, J. (1990). Plasticity Theory, Dover Publications."},{"key":"ref_49","doi-asserted-by":"crossref","first-page":"744","DOI":"10.1061\/(ASCE)0733-9437(1989)115:4(744)","article-title":"On the Correspondence between Brooks-Corey and van Genuchten Models","volume":"115","author":"Lenhard","year":"1989","journal-title":"J. Irrig. Drain. Eng."},{"key":"ref_50","doi-asserted-by":"crossref","first-page":"892","DOI":"10.2136\/sssaj1980.03615995004400050002x","article-title":"A closed form equation for predicting the hydrauclic conductivity of unsaturated soils","volume":"44","year":"1980","journal-title":"Soil Sci. Soc. Am. J."}],"container-title":["Computation"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/www.mdpi.com\/2079-3197\/8\/4\/98\/pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2025,10,11]],"date-time":"2025-10-11T10:32:59Z","timestamp":1760178779000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.mdpi.com\/2079-3197\/8\/4\/98"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2020,11,13]]},"references-count":50,"journal-issue":{"issue":"4","published-online":{"date-parts":[[2020,12]]}},"alternative-id":["computation8040098"],"URL":"https:\/\/doi.org\/10.3390\/computation8040098","relation":{},"ISSN":["2079-3197"],"issn-type":[{"type":"electronic","value":"2079-3197"}],"subject":[],"published":{"date-parts":[[2020,11,13]]}}}