Vertically averaged multi-constituent flow simulations of geological CO2 sequestration — Stabilized finite element methods and quadratic elements

IF 4.4 2区数学 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Mathematics and Computers in Simulation Pub Date : 2025-03-31 DOI:10.1016/j.matcom.2025.03.023

Chris Ladubec , Robert Gracie

{"title":"Vertically averaged multi-constituent flow simulations of geological CO2 sequestration — Stabilized finite element methods and quadratic elements","authors":"Chris Ladubec , Robert Gracie","doi":"10.1016/j.matcom.2025.03.023","DOIUrl":null,"url":null,"abstract":"<div><div>This paper compares several stabilized Finite Element Methods (FEMs) at solving the vertically averaged multi-constituent flow for CO<sub>2</sub> sequestration. Non-physical oscillations can occur when using the Galerkin FEM. Linear elements have been previously studied and are not able to provide a complete representation of the stabilization term in the consistent stabilized FEMs. Quadratic elements solve this problem and also allow additional formulations to be compared.</div><div>Stability is not guaranteed when using implicit time integration on non-linear advection dominated problems. Therefore, several stabilization methods are considered: the Streamline-Upwind (SU and SU-<span><math><mrow><mn>8</mn><mi>τ</mi></mrow></math></span>) method, the Streamline-Upwind Petrov–Galerkin (SUPG) method, the Galerkin Least Squares (GLS) method, the Subgrid Scales (SGS) method, and the Hughes Variational Method (HVM).</div><div>Quadratic elements and consistent formulations improve stability when discretizing in time first. Inconsistent methods (SU and SU-<span><math><mrow><mn>8</mn><mi>τ</mi></mrow></math></span>) can perform well with linear or quadratic elements, however they are especially susceptible to over-diffusion. Suitable stabilization methods can be obtained for various injection scenarios, by selecting an appropriate combination of stabilization method, stabilization parameter, Courant number and element size.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"235 ","pages":"Pages 114-131"},"PeriodicalIF":4.4000,"publicationDate":"2025-03-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Mathematics and Computers in Simulation","FirstCategoryId":"100","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0378475425000977","RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}

引用次数: 0

Abstract

This paper compares several stabilized Finite Element Methods (FEMs) at solving the vertically averaged multi-constituent flow for CO₂ sequestration. Non-physical oscillations can occur when using the Galerkin FEM. Linear elements have been previously studied and are not able to provide a complete representation of the stabilization term in the consistent stabilized FEMs. Quadratic elements solve this problem and also allow additional formulations to be compared.

Stability is not guaranteed when using implicit time integration on non-linear advection dominated problems. Therefore, several stabilization methods are considered: the Streamline-Upwind (SU and SU-

8 τ

) method, the Streamline-Upwind Petrov–Galerkin (SUPG) method, the Galerkin Least Squares (GLS) method, the Subgrid Scales (SGS) method, and the Hughes Variational Method (HVM).

Quadratic elements and consistent formulations improve stability when discretizing in time first. Inconsistent methods (SU and SU-

8 τ

) can perform well with linear or quadratic elements, however they are especially susceptible to over-diffusion. Suitable stabilization methods can be obtained for various injection scenarios, by selecting an appropriate combination of stabilization method, stabilization parameter, Courant number and element size.

查看原文本刊更多论文

求助全文

约1分钟内获得全文求助全文

来源期刊

Mathematics and Computers in Simulation 数学-计算机：跨学科应用

CiteScore

8.90

自引率

4.30%

发文量

335

审稿时长

54 days

期刊介绍： The aim of the journal is to provide an international forum for the dissemination of up-to-date information in the fields of the mathematics and computers, in particular (but not exclusively) as they apply to the dynamics of systems, their simulation and scientific computation in general. Published material ranges from short, concise research papers to more general tutorial articles. Mathematics and Computers in Simulation, published monthly, is the official organ of IMACS, the International Association for Mathematics and Computers in Simulation (Formerly AICA). This Association, founded in 1955 and legally incorporated in 1956 is a member of FIACC (the Five International Associations Coordinating Committee), together with IFIP, IFAV, IFORS and IMEKO. Topics covered by the journal include mathematical tools in: •The foundations of systems modelling •Numerical analysis and the development of algorithms for simulation They also include considerations about computer hardware for simulation and about special software and compilers. The journal also publishes articles concerned with specific applications of modelling and simulation in science and engineering, with relevant applied mathematics, the general philosophy of systems simulation, and their impact on disciplinary and interdisciplinary research. The journal includes a Book Review section -- and a "News on IMACS" section that contains a Calendar of future Conferences/Events and other information about the Association.