{"title":"High-order spline finite element method for solving time-dependent electromagnetic waves","authors":"Imad El-Barkani , Imane El-Hadouti , Mohamed Addam , Mohammed Seaid","doi":"10.1016/j.apnum.2024.08.002","DOIUrl":"10.1016/j.apnum.2024.08.002","url":null,"abstract":"<div><p>In this paper we propose a high-order spline finite element method for solving a class of time-dependent electromagnetic waves and its associated frequency-domain approach. A Fourier transform and its inverse are used for the time integration of the wave problem. The spatial discretization is performed using a partitioned mesh with tensorial spline functions to form bases of the discrete solution in the variational finite element space. Quadrature methods such as the Gauss-Hermite quadrature are implemented in the inverse Fourier transform to compute numerical solutions of the time-dependent electromagnetic waves. In the present study we carry out a rigorous convergence analysis and establish error estimates for the wave solution in the relevant norms. We also provide a full algorithmic description of the method and assess its performance by solving several test examples of time-dependent electromagnetic waves with known analytical solutions. The method is shown to verify the theoretical estimates and to provide highly accurate and efficient simulations. We also evaluate the computational performance of the proposed method for solving a problem of wave transmission through non-homogeneous materials. The obtained computational results confirm the excellent convergence, high accuracy and applicability of the proposed spline finite element method.</p></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"206 ","pages":"Pages 48-74"},"PeriodicalIF":2.2,"publicationDate":"2024-08-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141979538","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"The virtual element method for a contact problem with wear and unilateral constraint","authors":"Bangmin Wu , Fei Wang , Weimin Han","doi":"10.1016/j.apnum.2024.08.004","DOIUrl":"10.1016/j.apnum.2024.08.004","url":null,"abstract":"<div><p>This paper is dedicated to the numerical solution of a mathematical model that describes frictional quasistatic contact between an elastic body and a moving foundation, with the wear effect on the contact interface of the moving foundation due to friction. The mathematical problem is a system consisting of a time-dependent quasi-variational inequality and an integral equation. The numerical method is based on the use of the virtual element method (VEM) for the spatial discretization of the variational inequality and a variable step-size left rectangle integration formula for the integral equation. The existence and uniqueness of a numerical solution are shown, and optimal order error estimates are derived for both the displacement and the wear function for the lowest order VEM. Numerical results are presented to demonstrate the efficiency of the method and to illustrate the numerical convergence orders.</p></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"206 ","pages":"Pages 29-47"},"PeriodicalIF":2.2,"publicationDate":"2024-08-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141979540","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Convergence analysis of higher-order approximation of singularly perturbed 2D semilinear parabolic PDEs with non-homogeneous boundary conditions","authors":"Narendra Singh Yadav , Kaushik Mukherjee","doi":"10.1016/j.apnum.2024.08.001","DOIUrl":"10.1016/j.apnum.2024.08.001","url":null,"abstract":"<div><p>This article focuses on developing and analyzing an efficient higher-order numerical approximation of singularly perturbed two-dimensional semilinear parabolic convection-diffusion problems with time-dependent boundary conditions. We approximate the governing nonlinear problem by an implicit fitted mesh method (FMM), which combines an alternating direction implicit scheme in the temporal direction together with a higher-order finite difference scheme in the spatial directions. Since the solution possesses exponential boundary layers, a Cartesian product of piecewise-uniform Shishkin meshes is used to discretize in space. To begin our analysis, we establish the stability corresponding to the continuous nonlinear problem, and obtain a-priori bounds for the solution derivatives. Thereafter, we pursue the stability analysis of the discrete problem, and prove <em>ε</em>-uniform convergence in the maximum-norm. Next, for enhancement of the temporal accuracy, we use the Richardson extrapolation technique solely in the temporal direction. In addition, we investigate the order reduction phenomenon naturally occurring due to the time-dependent boundary data and propose a suitable approximation to tackle this effect. Finally, we present the computational results to validate the theoretical estimates.</p></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"206 ","pages":"Pages 210-246"},"PeriodicalIF":2.2,"publicationDate":"2024-08-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142040287","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A modified PRP conjugate gradient method for unconstrained optimization and nonlinear equations","authors":"Haijuan Cui","doi":"10.1016/j.apnum.2024.07.014","DOIUrl":"10.1016/j.apnum.2024.07.014","url":null,"abstract":"<div><p>A modified Polak Ribiere Polyak(PRP) conjugate gradient(CG) method is proposed for solving unconstrained optimization problems. The search direction generated by this method satisfies sufficient descent condition at each iteration and this method inherits one remarkable property of the standard PRP method. Under the standard Armijo line search, the global convergence and the linearly convergent rate of the presented method is established. Some numerical results are given to show the effectiveness of the proposed method by comparing with some existing CG methods.</p></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"205 ","pages":"Pages 296-307"},"PeriodicalIF":2.2,"publicationDate":"2024-08-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141935229","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A second-order structure-preserving discretization for the Cahn-Hilliard/Allen-Cahn system with cross-kinetic coupling","authors":"Aaron Brunk , Herbert Egger , Oliver Habrich","doi":"10.1016/j.apnum.2024.07.016","DOIUrl":"10.1016/j.apnum.2024.07.016","url":null,"abstract":"<div><p>We study the numerical solution of a Cahn-Hilliard/Allen-Cahn system with strong coupling through state and gradient dependent non-diagonal mobility matrices. A fully discrete approximation scheme in space and time is proposed which preserves the underlying gradient flow structure and leads to dissipation of the free-energy on the discrete level. Existence and uniqueness of the discrete solution is established and relative energy estimates are used to prove optimal convergence rates in space and time under minimal smoothness assumptions. Numerical tests are presented for illustration of the theoretical results and to demonstrate the viability of the proposed methods.</p></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"206 ","pages":"Pages 12-28"},"PeriodicalIF":2.2,"publicationDate":"2024-08-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S016892742400196X/pdfft?md5=c541e871e3b0cf5488415210f7fead7d&pid=1-s2.0-S016892742400196X-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141935227","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Numerical simulation for an initial-boundary value problem of time-fractional Klein-Gordon equations","authors":"Zaid Odibat","doi":"10.1016/j.apnum.2024.07.015","DOIUrl":"10.1016/j.apnum.2024.07.015","url":null,"abstract":"<div><p>This paper mainly presents numerical solutions to an initial-boundary value problem of the time-fractional Klein-Gordon equations. We developed a numerical scheme with the help of the finite difference methods and the predictor-corrector methods to find numerical solutions of the considered problems. The proposed scheme is based on discretizing the considered problems with respect to spatial and temporal domains. Numerical results are derived for some illustrative problems, and the outputs are compared with the exact solution in the integer order case. The solution behavior and 3D graphics of the discussed problems are demonstrated using the proposed scheme. Finally, the proposed scheme, which does not require solving large systems of linear equations, can be extended and modified to handle other classes of time-fractional PDEs.</p></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"206 ","pages":"Pages 1-11"},"PeriodicalIF":2.2,"publicationDate":"2024-08-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141935228","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"An efficient iterative method for solving the graph regularization Q-weighted nonnegative matrix factorization problem in multi-view clustering","authors":"Chunmei Li, Dan Tian, Xuefeng Duan, Naya Yang","doi":"10.1016/j.apnum.2024.07.010","DOIUrl":"10.1016/j.apnum.2024.07.010","url":null,"abstract":"<div><p>In this paper, we consider the graph regularization Q-weighted nonnegative matrix factorization problem in multi-view clustering. Based on the Q-weighted norm property, this problem is transformed into the minimization problem of the trace function. The necessary condition for the existence of a solution is given. The proximal alternating nonnegative least squares method and its acceleration method are designed to solve it. The convergence theorem is also given. The feasibility and effectiveness of the proposed methods are verified by numerical experiments.</p></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"205 ","pages":"Pages 255-266"},"PeriodicalIF":2.2,"publicationDate":"2024-08-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141935042","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Multi-step Hermite-Birkhoff predictor-corrector schemes","authors":"Arjun Thenery Manikantan, Jochen Schütz","doi":"10.1016/j.apnum.2024.07.011","DOIUrl":"10.1016/j.apnum.2024.07.011","url":null,"abstract":"<div><p>In this study, we introduce a multi-step multi-derivative predictor-corrector time integration scheme analogous to the schemes in Schütz et al. (2022) <span><span>[13]</span></span>, incorporating a multi-step quadrature rule. We conduct stability analysis up to order eight and optimize the schemes to achieve <span><math><mi>A</mi><mo>(</mo><mi>α</mi><mo>)</mo></math></span>-stability for large <em>α</em>. Numerical experiments are performed on ordinary differential equations exhibiting diverse stiffness conditions, as well as on partial differential equations showcasing non-linearity and higher-order terms. Results demonstrate the convergence and flexibility of the proposed schemes across diverse situations.</p></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"205 ","pages":"Pages 281-295"},"PeriodicalIF":2.2,"publicationDate":"2024-07-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141935044","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A rotational velocity-correction projection method for the Micropolar Navier-Stokes equations","authors":"Zhiyong Si , Ziyi Li , Leilei Wei","doi":"10.1016/j.apnum.2024.07.013","DOIUrl":"10.1016/j.apnum.2024.07.013","url":null,"abstract":"<div><p>In this paper, we introduce a velocity correction projection method for the Micropolar Navier-Stokes Equations. The velocity correction method are adopted to approximate the time derivative term, stability analysis and error estimation of the first-order semi-discrete scheme are proved. At the same time, the optimal error estimate using the technique of dual norm are obtained. In this way, the divergence free of the velocity <strong>u</strong> can be conserved. Finally, the numerical results show the method has an optimal convergence order. The numerical results are consistent with our theoretical analysis, and our method is effective.</p></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"205 ","pages":"Pages 267-280"},"PeriodicalIF":2.2,"publicationDate":"2024-07-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141935043","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A reduced-order two-grid method based on POD technique for the semilinear parabolic equation","authors":"Junpeng Song, Hongxing Rui","doi":"10.1016/j.apnum.2024.07.012","DOIUrl":"10.1016/j.apnum.2024.07.012","url":null,"abstract":"<div><p>In the conventional two-grid (TG) method, the nonlinear system on the fine grid is transformed into a nonlinear subsystem on the coarse grid and a linear subsystem on the fine grid to reduce computational costs. It has been successfully applied in various fields. Nonetheless, its computational efficiency remains relatively low. For this, we develop a novel reduced-order two-grid (ROTG) method with less degrees of freedom for solving the semilinear parabolic equation. For the two subsystems mentioned, the proper orthogonal decomposition (POD) technique is utilized to substantially reduce degrees of freedom. An a priori error estimate for the ROTG scheme is derived. Finally, we conduct several numerical tests to observe the ROTG method's behavior and verify the theoretical analysis.</p></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"205 ","pages":"Pages 240-254"},"PeriodicalIF":2.2,"publicationDate":"2024-07-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141935045","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}