{"title":"Unconditionally Energy Stable and First-Order Accurate Numerical Schemes for the Heat Equation with Uncertain Temperature-Dependent Conductivity","authors":"Fiordilino, J. A., Winger, M.","doi":"10.4208/ijnam2023-1035","DOIUrl":"https://doi.org/10.4208/ijnam2023-1035","url":null,"abstract":"In this paper, we present first-order accurate numerical methods for solution of the heat equation with uncertain temperature-dependent thermal conductivity. Each algorithm yields a shared coefficient matrix for the ensemble set improving computational efficiency. Both mixed and Robin-type boundary conditions are treated. In contrast with alternative, related methodologies, stability and convergence are unconditional. In particular, we prove unconditional, energy stability and optimal-order error estimates. A battery of numerical tests are presented to illustrate both the theory and application of these algorithms.","PeriodicalId":50301,"journal":{"name":"International Journal of Numerical Analysis and Modeling","volume":"44 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135219537","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A Novel Deep Neural Network Algorithm for the Helmholtz Scattering Problem In the Unbounded Domain","authors":"Andy L Yang","doi":"10.4208/ijnam2023-1032","DOIUrl":"https://doi.org/10.4208/ijnam2023-1032","url":null,"abstract":". In this paper, we develop a novel meshless, ray-based deep neural network algorithm for solving the high-frequency Helmholtz scattering problem in the unbounded domain. While our recent work [44] designed a deep neural network method for solving the Helmholtz equation over (cid:12)nite bounded domains, this paper deals with the more general and di(cid:14)cult case of unbounded regions. By using the perfectly matched layer method, the original mathematical model in the unbounded domain is transformed into a new format of second-order system in a (cid:12)nite bounded domain with simple homogeneous Dirichlet boundary conditions. Compared with the Helmholtz equation in the bounded domain, the new system is equipped with variable coe(cid:14)cients. Then, a deep neural network algorithm is designed for the new system, where the rays in various random directions are used as the basis of the numerical solution. Various numerical examples have been carried out to demonstrate the accuracy and e(cid:14)ciency of the proposed numerical method. The proposed method has the advantage of easy implementation and meshless while maintaining high accuracy. To the best of the author’s knowledge, this is the (cid:12)rst deep neural network method to solve the Helmholtz equation in the unbounded domain.","PeriodicalId":50301,"journal":{"name":"International Journal of Numerical Analysis and Modeling","volume":"7 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135143752","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A Conforming Dg Method for the Biharmonic Equation on Polytopal Meshes","authors":"Xiu Ye null, Shangyou Zhang","doi":"10.4208/ijnam2023-1037","DOIUrl":"https://doi.org/10.4208/ijnam2023-1037","url":null,"abstract":"A conforming discontinuous Galerkin finite element method is introduced for solving the biharmonic equation. This method, by its name, uses discontinuous approximations and keeps simple formulation of the conforming finite element method at the same time. The ultra simple formulation of the method will reduce programming complexity in practice. Optimal order error estimates in a discrete $H^2$ norm is established for the corresponding finite element solutions. Error estimates in the $L^2$ norm are also derived with a sub-optimal order of convergence for the lowest order element and an optimal order of convergence for all high order of elements. Numerical results are presented to confirm the theory of convergence.","PeriodicalId":50301,"journal":{"name":"International Journal of Numerical Analysis and Modeling","volume":"12 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135195304","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Champike Attanayake, So-Hsiang Chou null, Quanling Deng
{"title":"High-Order Enriched Finite Element Methods for Elliptic Interface Problems with Discontinuous Solutions","authors":"Champike Attanayake, So-Hsiang Chou null, Quanling Deng","doi":"10.4208/ijnam2023-1038","DOIUrl":"https://doi.org/10.4208/ijnam2023-1038","url":null,"abstract":"Elliptic interface problems whose solutions are $C^0$ continuous have been well studied over the past two decades. The well-known numerical methods include the strongly stable generalized finite element method (SGFEM) and immersed FEM (IFEM). In this paper, we study numerically a larger class of elliptic interface problems where their solutions are discontinuous. A direct application of these existing methods fails immediately as the approximate solution is in a larger space that covers discontinuous functions. We propose a class of high-order enriched unfitted FEMs to solve these problems with implicit or Robin-type interface jump conditions. We design new enrichment functions that capture the imposed discontinuity of the solution while keeping the condition number from fast growth. A linear enriched method in 1D was recently developed using one enrichment function and we generalized it to an arbitrary degree using two simple discontinuous one-sided enrichment functions. The natural tensor product extension to the 2D case is demonstrated. Optimal order convergence in the $L^2$ and broken $H^1$-norms are established. We also establish superconvergence at all discretization nodes (including exact nodal values in special cases). Numerical examples are provided to confirm the theory. Finally, to prove the efficiency of the method for practical problems, the enriched linear, quadratic, and cubic elements are applied to a multi-layer wall model for drug-eluting stents in which zero-flux jump conditions and implicit concentration interface conditions are both present.","PeriodicalId":50301,"journal":{"name":"International Journal of Numerical Analysis and Modeling","volume":"12 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135195306","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
P.-A. Gourdain, M. Evans, H. R. Hasson, J. R. Young, I. West-Abdallah null, M. B. Adams
{"title":"A Fully Implicit Method Using Nodal Radial Basis Functions to Solve the Linear Advection Equation","authors":"P.-A. Gourdain, M. Evans, H. R. Hasson, J. R. Young, I. West-Abdallah null, M. B. Adams","doi":"10.4208/ijnam2023-1018","DOIUrl":"https://doi.org/10.4208/ijnam2023-1018","url":null,"abstract":"","PeriodicalId":50301,"journal":{"name":"International Journal of Numerical Analysis and Modeling","volume":"10 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135421665","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Sean Ingimarson, Monika Neda, Leo G. Rebholz, Jorge Reyes null, An Vu
{"title":"Improved Long Time Accuracy for Projection Methods for Navier-Stokes Equations Using Emac Formulation","authors":"Sean Ingimarson, Monika Neda, Leo G. Rebholz, Jorge Reyes null, An Vu","doi":"10.4208/ijnam2023-1008","DOIUrl":"https://doi.org/10.4208/ijnam2023-1008","url":null,"abstract":"","PeriodicalId":50301,"journal":{"name":"International Journal of Numerical Analysis and Modeling","volume":"64 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"136370904","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A Sharp $alpha$-Robust $L1$ Scheme on Graded Meshes for Two-Dimensional Time Tempered Fractional Fokker-Planck Equation","authors":"Can Wang, Weihua Deng null, Xiangong Tang","doi":"10.4208/ijnam2023-1033","DOIUrl":"https://doi.org/10.4208/ijnam2023-1033","url":null,"abstract":"","PeriodicalId":50301,"journal":{"name":"International Journal of Numerical Analysis and Modeling","volume":"283 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135195302","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A Posteriori Error Estimates for a Local Discontinuous Galerkin Approximation of Semilinear Second-Order Elliptic Problems on Cartesian Grids","authors":"Mahboub Baccouch","doi":"10.4208/ijnam2023-1034","DOIUrl":"https://doi.org/10.4208/ijnam2023-1034","url":null,"abstract":". In this paper, we design and analyze new residual-type a posteriori error estimators for the local discontinuous Galerkin (LDG) method applied to semilinear second-order elliptic problems in two dimensions of the type (cid:0) ∆ u = f ( x ;u ). We use our recent superconvergence results derived in Commun. Appl. Math. Comput. (2021) to prove that the LDG solution is superconvergent with an order p +2 towards the p -degree right Radau interpolating polynomial of the exact solution, when tensor product polynomials of degree at most p are considered as basis for the LDG method. Moreover, we show that the global discretization error can be decomposed into the sum of two errors. The first error can be expressed as a linear combination of two ( p +1)- degree Radau polynomials in the x - and y (cid:0) directions. The second error converges to zero with order p + 2 in the L 2 -norm. This new result allows us to construct a posteriori error estimators of residual type. We prove that the proposed a posteriori error estimators converge to the true errors in the L 2 -norm under mesh refinement at the optimal rate. The order of convergence is proved to be p + 2. We further prove that our a posteriori error estimates yield upper and lower bounds for the actual error. Finally, a series of numerical examples are presented to validate the theoretical results and numerically demonstrate the convergence of the proposed a posteriori error estimators.","PeriodicalId":50301,"journal":{"name":"International Journal of Numerical Analysis and Modeling","volume":"12 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135195303","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A Finite Volume Element Solution Based on Postprocessing Technique Over Arbitrary Convex Polygonal Meshes","authors":"Yanlong Zhang null, Yanhui Zhou","doi":"10.4208/ijnam2023-1026","DOIUrl":"https://doi.org/10.4208/ijnam2023-1026","url":null,"abstract":". A special (cid:12)nite volume element method based on postprocessing technique is proposed to solve the anisotropic di(cid:11)usion problem on arbitrary convex polygonal meshes. The shape function of polygonal (cid:12)nite element method is constructed by Wachspress generalized barycentric coordinate, and by adding some element-wise bubble functions to the (cid:12)nite element solution, we get a new (cid:12)nite volume element solution that satis(cid:12)es the local conservation law on a certain dual mesh. The postprocessing algorithm only needs to solve a local linear algebraic system on each primary cell, so that it is easy to implement. More interesting is that, a general construction of the bubble functions is introduced on each polygonal cell, which enables us to prove the existence and uniqueness of the post-processed solution on arbitrary convex polygonal meshes with full anisotropic di(cid:11)usion tensor. The optimal H 1 and L 2 error estimates of the post-processed solution are also obtained. Finally, the local conservation property and convergence of the new polygonal (cid:12)nite volume element solution are veri(cid:12)ed by numerical experiments.","PeriodicalId":50301,"journal":{"name":"International Journal of Numerical Analysis and Modeling","volume":"1 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135143757","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}