{"title":"针对 Barenblatt-Biot 孔弹性模型的新型混合有限元法分析","authors":"","doi":"10.1016/j.camwa.2024.07.011","DOIUrl":null,"url":null,"abstract":"<div><p>In this work, we study the locking-free numerical method for a Barenblatt-Biot poroelastic model. When solving by the continuous Galerkin mixed finite element method, the model exists two kind of locking phenomena for special physical parameters. To overcome these locking phenomena, we introduce new variables to reformulate the original problem into a new problem, which exists a built-in mechanism to keep the continuous Galerkin mixed finite element method stable. It can be regarded as a generalized Stokes problem for two given <em>β</em> and <em>η</em> and two diffusion problems for a given <em>δ</em>. The generalized Stokes problem can be adopted by some stable solver and the diffusion problems can be solved by continuous Galerkin finite element. Moreover, the existence and uniqueness of weak solution is proved by using the standard Galerkin method and combining with priori estimates and some invariant quantities. After that, we design fully discrete time-stepping schemes to use mixed finite element method with <span><math><msub><mrow><mi>P</mi></mrow><mrow><mn>2</mn></mrow></msub><mo>−</mo><msub><mrow><mi>P</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>−</mo><msub><mrow><mi>P</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>−</mo><msub><mrow><mi>P</mi></mrow><mrow><mn>1</mn></mrow></msub></math></span> element pairs for space variables and backward Euler method for time variable, and analyze the coupled and decoupled time stepping methods based on the proposed scheme. The optimal convergence order is obtained in both space and time. Finally, some numerical examples are presented to show the optimal convergence rates about variables and the robustness of proposed method with respect to <em>ν</em>, and to verify that there is no locking phenomenon.</p></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":null,"pages":null},"PeriodicalIF":2.9000,"publicationDate":"2024-07-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"Analysis of new mixed finite element method for a Barenblatt-Biot poroelastic model\",\"authors\":\"\",\"doi\":\"10.1016/j.camwa.2024.07.011\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"<div><p>In this work, we study the locking-free numerical method for a Barenblatt-Biot poroelastic model. When solving by the continuous Galerkin mixed finite element method, the model exists two kind of locking phenomena for special physical parameters. To overcome these locking phenomena, we introduce new variables to reformulate the original problem into a new problem, which exists a built-in mechanism to keep the continuous Galerkin mixed finite element method stable. It can be regarded as a generalized Stokes problem for two given <em>β</em> and <em>η</em> and two diffusion problems for a given <em>δ</em>. The generalized Stokes problem can be adopted by some stable solver and the diffusion problems can be solved by continuous Galerkin finite element. Moreover, the existence and uniqueness of weak solution is proved by using the standard Galerkin method and combining with priori estimates and some invariant quantities. After that, we design fully discrete time-stepping schemes to use mixed finite element method with <span><math><msub><mrow><mi>P</mi></mrow><mrow><mn>2</mn></mrow></msub><mo>−</mo><msub><mrow><mi>P</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>−</mo><msub><mrow><mi>P</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>−</mo><msub><mrow><mi>P</mi></mrow><mrow><mn>1</mn></mrow></msub></math></span> element pairs for space variables and backward Euler method for time variable, and analyze the coupled and decoupled time stepping methods based on the proposed scheme. The optimal convergence order is obtained in both space and time. Finally, some numerical examples are presented to show the optimal convergence rates about variables and the robustness of proposed method with respect to <em>ν</em>, and to verify that there is no locking phenomenon.</p></div>\",\"PeriodicalId\":55218,\"journal\":{\"name\":\"Computers & Mathematics with Applications\",\"volume\":null,\"pages\":null},\"PeriodicalIF\":2.9000,\"publicationDate\":\"2024-07-22\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"Computers & Mathematics with Applications\",\"FirstCategoryId\":\"100\",\"ListUrlMain\":\"https://www.sciencedirect.com/science/article/pii/S0898122124003109\",\"RegionNum\":2,\"RegionCategory\":\"数学\",\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q1\",\"JCRName\":\"MATHEMATICS, APPLIED\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"Computers & Mathematics with Applications","FirstCategoryId":"100","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0898122124003109","RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"MATHEMATICS, APPLIED","Score":null,"Total":0}
Analysis of new mixed finite element method for a Barenblatt-Biot poroelastic model
In this work, we study the locking-free numerical method for a Barenblatt-Biot poroelastic model. When solving by the continuous Galerkin mixed finite element method, the model exists two kind of locking phenomena for special physical parameters. To overcome these locking phenomena, we introduce new variables to reformulate the original problem into a new problem, which exists a built-in mechanism to keep the continuous Galerkin mixed finite element method stable. It can be regarded as a generalized Stokes problem for two given β and η and two diffusion problems for a given δ. The generalized Stokes problem can be adopted by some stable solver and the diffusion problems can be solved by continuous Galerkin finite element. Moreover, the existence and uniqueness of weak solution is proved by using the standard Galerkin method and combining with priori estimates and some invariant quantities. After that, we design fully discrete time-stepping schemes to use mixed finite element method with element pairs for space variables and backward Euler method for time variable, and analyze the coupled and decoupled time stepping methods based on the proposed scheme. The optimal convergence order is obtained in both space and time. Finally, some numerical examples are presented to show the optimal convergence rates about variables and the robustness of proposed method with respect to ν, and to verify that there is no locking phenomenon.
期刊介绍:
Computers & Mathematics with Applications provides a medium of exchange for those engaged in fields contributing to building successful simulations for science and engineering using Partial Differential Equations (PDEs).