Anatoly Alikhanov , Huiran Bai , Jian Huang , Aleksei Tyrylgin , Yin Yang
{"title":"Multiscale model reduction for the time fractional thermoporoelasticity problem in fractured and heterogeneous media","authors":"Anatoly Alikhanov , Huiran Bai , Jian Huang , Aleksei Tyrylgin , Yin Yang","doi":"10.1016/j.cam.2024.116157","DOIUrl":null,"url":null,"abstract":"<div><p>In this paper, we consider the time fractional thermoporoelasticity problem in fractured and heterogeneous media. The mathematical model with a time memory formalism is described by a coupled system of equations for pressure, temperature and displacements. We use an implicit finite difference approximation for temporal discretization. We present a fine grid approximation based on the finite element method and Discrete Fracture Model (DFM) for two-dimensional model problems. Further, we use the Generalized Multiscale Finite Element Method (GMsFEM) for coarse grid approximation. The primary concept behind the proposed method is to streamline the complexity inherent in the thermoporoelasticity problem. Given that our model equation incorporates multiple fractional powers, leading to multiple unknowns with memory effects, we aim to address this intricacy by optimizing the problem’s dimensionality. As a result, the solution is sought on a coarse grid, a strategic choice that not only simplifies the computational cost but also contributes to significant time savings. We present numerical results for the two-dimensional model problems in heterogeneous fractured porous media. We derive relative errors between the reference fine grid solution and the multiscale solution for different numbers of multiscale basis functions. The results confirm that the proposed method is able to achieve good accuracy with a few degrees of freedoms on the coarse grid.</p></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"455 ","pages":"Article 116157"},"PeriodicalIF":2.1000,"publicationDate":"2024-07-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Computational and Applied Mathematics","FirstCategoryId":"100","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0377042724004060","RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"MATHEMATICS, APPLIED","Score":null,"Total":0}
引用次数: 0
Abstract
In this paper, we consider the time fractional thermoporoelasticity problem in fractured and heterogeneous media. The mathematical model with a time memory formalism is described by a coupled system of equations for pressure, temperature and displacements. We use an implicit finite difference approximation for temporal discretization. We present a fine grid approximation based on the finite element method and Discrete Fracture Model (DFM) for two-dimensional model problems. Further, we use the Generalized Multiscale Finite Element Method (GMsFEM) for coarse grid approximation. The primary concept behind the proposed method is to streamline the complexity inherent in the thermoporoelasticity problem. Given that our model equation incorporates multiple fractional powers, leading to multiple unknowns with memory effects, we aim to address this intricacy by optimizing the problem’s dimensionality. As a result, the solution is sought on a coarse grid, a strategic choice that not only simplifies the computational cost but also contributes to significant time savings. We present numerical results for the two-dimensional model problems in heterogeneous fractured porous media. We derive relative errors between the reference fine grid solution and the multiscale solution for different numbers of multiscale basis functions. The results confirm that the proposed method is able to achieve good accuracy with a few degrees of freedoms on the coarse grid.
期刊介绍:
The Journal of Computational and Applied Mathematics publishes original papers of high scientific value in all areas of computational and applied mathematics. The main interest of the Journal is in papers that describe and analyze new computational techniques for solving scientific or engineering problems. Also the improved analysis, including the effectiveness and applicability, of existing methods and algorithms is of importance. The computational efficiency (e.g. the convergence, stability, accuracy, ...) should be proved and illustrated by nontrivial numerical examples. Papers describing only variants of existing methods, without adding significant new computational properties are not of interest.
The audience consists of: applied mathematicians, numerical analysts, computational scientists and engineers.