{"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":null,"pages":null},"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":null,"pages":null},"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":null,"pages":null},"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":null,"pages":null},"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":null,"pages":null},"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":null,"pages":null},"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":null,"pages":null},"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":null,"pages":null},"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}
{"title":"Tethered flexible polymer under oscillatory linear flow","authors":"A. Lamura","doi":"10.1016/j.apnum.2024.07.009","DOIUrl":"10.1016/j.apnum.2024.07.009","url":null,"abstract":"<div><p>The non-equilibrium structural and dynamical properties of a flexible polymer tethered to a reflecting wall and subject to oscillatory linear flow are studied by numerical simulations. Polymer is confined in two dimensions and is modeled as a bead-spring chain immersed in a fluid described by the Brownian multiparticle collision dynamics. At high strain, the polymer is stretched along the flow direction following the applied flow, then recoils at flow inversion before flipping and elongate again. When strain is reduced, it may happen that the chain recoils without flipping when the applied shear changes sign. Conformations are analyzed and compared to stiff polymers revealing more compact patterns at low strains and less stretched configurations at high strain. The dynamics is investigated by looking at the center-of-mass motion which shows a frequency doubling along the direction normal to the external flow. The center-of-mass correlation function is characterized by smaller amplitudes when reducing bending rigidity.</p></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":null,"pages":null},"PeriodicalIF":2.2,"publicationDate":"2024-07-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141817929","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":"Numerical solution of a hydrodynamic model with cavitation using finite difference method at arbitrary meshes","authors":"A. García , M. Negreanu , F. Ureña , A.M. Vargas","doi":"10.1016/j.apnum.2024.07.007","DOIUrl":"10.1016/j.apnum.2024.07.007","url":null,"abstract":"<div><p>In this paper, we investigate the implementation of the finite difference method on arbitrary meshes in conjunction with a fixed-point algorithm for the lubrication problem of a journal bearing with cavitation, considering the Elrod-Adams model. We establish numerical properties of the generalized finite difference scheme and provide several illustrative examples.</p></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":null,"pages":null},"PeriodicalIF":2.2,"publicationDate":"2024-07-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0168927424001855/pdfft?md5=b0297d53fcb365d790dc28df2d29a8cb&pid=1-s2.0-S0168927424001855-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141839198","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}