{"title":"General enrichments of stable GFEM for interface problems: Theory and extreme learning machine construction","authors":"Dongmei Wang , Hengguang Li , Qinghui Zhang","doi":"10.1016/j.apnum.2025.03.009","DOIUrl":null,"url":null,"abstract":"<div><div>Generalized finite element methods (GFEMs), when applied to interface problems (IPs), need to be enriched with special functions to enhance approximation accuracy. These functions include distance functions, one-side distance functions, level set functions, and exponential forms of level set function. For the IP with geometrically complex interface curves, computation of the distance function or level set function could be challenging, and algorithms of computational geometry are usually involved. Moreover, theoretical analysis on optimal convergence of the GFEM enriched by these functions has not been fully investigated. In this study we propose a general enrichment scheme, based on which all the aforementioned enrichments can be viewed as special examples. We prove that a stable GFEM (SGFEM) with such a new enrichment scheme reaches the optimal convergence rate. Most importantly, the new scheme provides an instruction to construct machine learning (ML) based enrichments, which advances the ability of GFEM to handle geometrically complex interfaces. Two ML methods, deep neural network (DNN) and extreme learning machine (ELM), are studied. Among them, the ELM is highly suggested because it exhibits high accuracy for the interface curve with complex geometries. The learning dimension for the ML is one dimension less than that of the domain so that the proposed ML algorithm can be implemented efficiently. The numerical experiments demonstrate that the SGFEM with the ELM enrichment achieves the optimal convergence rates for the IP, as predicted theoretically.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"214 ","pages":"Pages 143-159"},"PeriodicalIF":2.2000,"publicationDate":"2025-04-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Applied Numerical Mathematics","FirstCategoryId":"100","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0168927425000716","RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"MATHEMATICS, APPLIED","Score":null,"Total":0}
引用次数: 0
Abstract
Generalized finite element methods (GFEMs), when applied to interface problems (IPs), need to be enriched with special functions to enhance approximation accuracy. These functions include distance functions, one-side distance functions, level set functions, and exponential forms of level set function. For the IP with geometrically complex interface curves, computation of the distance function or level set function could be challenging, and algorithms of computational geometry are usually involved. Moreover, theoretical analysis on optimal convergence of the GFEM enriched by these functions has not been fully investigated. In this study we propose a general enrichment scheme, based on which all the aforementioned enrichments can be viewed as special examples. We prove that a stable GFEM (SGFEM) with such a new enrichment scheme reaches the optimal convergence rate. Most importantly, the new scheme provides an instruction to construct machine learning (ML) based enrichments, which advances the ability of GFEM to handle geometrically complex interfaces. Two ML methods, deep neural network (DNN) and extreme learning machine (ELM), are studied. Among them, the ELM is highly suggested because it exhibits high accuracy for the interface curve with complex geometries. The learning dimension for the ML is one dimension less than that of the domain so that the proposed ML algorithm can be implemented efficiently. The numerical experiments demonstrate that the SGFEM with the ELM enrichment achieves the optimal convergence rates for the IP, as predicted theoretically.
期刊介绍:
The purpose of the journal is to provide a forum for the publication of high quality research and tutorial papers in computational mathematics. In addition to the traditional issues and problems in numerical analysis, the journal also publishes papers describing relevant applications in such fields as physics, fluid dynamics, engineering and other branches of applied science with a computational mathematics component. The journal strives to be flexible in the type of papers it publishes and their format. Equally desirable are:
(i) Full papers, which should be complete and relatively self-contained original contributions with an introduction that can be understood by the broad computational mathematics community. Both rigorous and heuristic styles are acceptable. Of particular interest are papers about new areas of research, in which other than strictly mathematical arguments may be important in establishing a basis for further developments.
(ii) Tutorial review papers, covering some of the important issues in Numerical Mathematics, Scientific Computing and their Applications. The journal will occasionally publish contributions which are larger than the usual format for regular papers.
(iii) Short notes, which present specific new results and techniques in a brief communication.