{"title":"Adaptive mixed virtual element method for the fourth-order singularly perturbed problem","authors":"Jian Meng , Xu Qian , Jiali Qiu , Jingmin Xia","doi":"10.1016/j.jcp.2025.113738","DOIUrl":null,"url":null,"abstract":"<div><div>The singularly perturbed theory mainly arises in the system of differential equations with the small enough perturbed parameters acting on the highest-order derivatives. In this paper, we introduce the adaptive mixed virtual element method for the fourth-order singularly perturbed problem and the associated eigenvalue problem. It allows to apply the <span><math><msup><mrow><mi>H</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span>-conforming virtual elements to discrete the continuous spaces and reduces the total number of required degrees of freedom of the <span><math><msup><mrow><mi>H</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span>-conforming virtual element method. Basically, the great flexibility of virtual element method becomes appealing in mesh refinement because the locally mesh post-processing to remove hanging nodes is never needed. This naturally motivates us to develop an <em>a posteriori</em> error estimate for the model problem. Based on the numerical solutions, the interior and edge residual terms, and the error terms related to the inconsistency of the virtual element scheme, the error estimators applied to adaptively refine meshes are constructed and then proved to be equivalent to numerical errors under the balanced energy norms. Moreover, we also consider the approximation method for the fourth-order singularly perturbed eigenvalue problem in two-dimensional space. Analogous with the source problem, we not only discuss the boundedness of the eigenfunctions, but also present the upper bound for the error of the approximated eigenvalues by these error estimators. Necessitated by supporting the theoretical analysis, representative numerical examples are reported. We show that the current numerical method converges at the optimal rate uniformly with respect to the singularly perturbed parameters when using the adaptive polygonal meshes.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"524 ","pages":"Article 113738"},"PeriodicalIF":3.8000,"publicationDate":"2025-01-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Computational Physics","FirstCategoryId":"101","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S002199912500021X","RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
引用次数: 0
Abstract
The singularly perturbed theory mainly arises in the system of differential equations with the small enough perturbed parameters acting on the highest-order derivatives. In this paper, we introduce the adaptive mixed virtual element method for the fourth-order singularly perturbed problem and the associated eigenvalue problem. It allows to apply the -conforming virtual elements to discrete the continuous spaces and reduces the total number of required degrees of freedom of the -conforming virtual element method. Basically, the great flexibility of virtual element method becomes appealing in mesh refinement because the locally mesh post-processing to remove hanging nodes is never needed. This naturally motivates us to develop an a posteriori error estimate for the model problem. Based on the numerical solutions, the interior and edge residual terms, and the error terms related to the inconsistency of the virtual element scheme, the error estimators applied to adaptively refine meshes are constructed and then proved to be equivalent to numerical errors under the balanced energy norms. Moreover, we also consider the approximation method for the fourth-order singularly perturbed eigenvalue problem in two-dimensional space. Analogous with the source problem, we not only discuss the boundedness of the eigenfunctions, but also present the upper bound for the error of the approximated eigenvalues by these error estimators. Necessitated by supporting the theoretical analysis, representative numerical examples are reported. We show that the current numerical method converges at the optimal rate uniformly with respect to the singularly perturbed parameters when using the adaptive polygonal meshes.
期刊介绍:
Journal of Computational Physics thoroughly treats the computational aspects of physical problems, presenting techniques for the numerical solution of mathematical equations arising in all areas of physics. The journal seeks to emphasize methods that cross disciplinary boundaries.
The Journal of Computational Physics also publishes short notes of 4 pages or less (including figures, tables, and references but excluding title pages). Letters to the Editor commenting on articles already published in this Journal will also be considered. Neither notes nor letters should have an abstract.