Applied Numerical Mathematics最新文献

{"title":"Energy stability and error estimate of the RKMK2e scheme for the extended Fisher–Kolmogorov equation","authors":"Haifeng Wang,&nbsp;Yan Wang,&nbsp;Hong Zhang,&nbsp;Songhe Song","doi":"10.1016/j.apnum.2025.01.014","DOIUrl":"10.1016/j.apnum.2025.01.014","url":null,"abstract":"<div><div>In this paper, we develop a second-order accurate scheme for the extended Fisher–Kolmogorov (EFK) equation and investigate its global-in-time energy stability and convergence. The proposed scheme uses the Fourier spectral collocation method in space and the stabilization Runge–Kutta–Munthe–Kaas-2e (RKMK2e) method for temporal approximation. To demonstrate the global-in-time energy stability of the proposed scheme, we first verify that, under the assumption that all numerical solutions are uniformly bounded, the scheme is energy stable when using a sufficiently large stabilization parameter. Then, to establish the uniform-in-time boundedness of the numerical solutions, we fully utilize the nonlinear operator estimates and discrete Sobolev embedding in each stage of the scheme. Moreover, we conduct an optimal rate convergence analysis with a sufficient regularity assumption for the exact solution. Several numerical examples are presented to validate the accuracy, computational efficiency, and energy stability of the proposed scheme.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"212 ","pages":"Pages 60-76"},"PeriodicalIF":2.2,"publicationDate":"2025-01-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143155010","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":"Refined variational iteration method through optimal linear operator and Lagrange multiplier: Local/global convergence","authors":"Tapas Roy,&nbsp;Dilip Kumar Maiti","doi":"10.1016/j.apnum.2025.01.013","DOIUrl":"10.1016/j.apnum.2025.01.013","url":null,"abstract":"<div><div>Our focus is on refining the Variational Iteration Method (VIM), which will henceforth be called R-VIM. We begin by introducing a general linear operator and consequently the Lagrange's multiplier; then, minimizing the residual error, we subsequently find the optimal linear operator as well as the best fitting Lagrange's multiplier for a given differential equation. By means of Banach Fixed Point Theorem, the necessary conditions for the convergence of the solutions of VIM, Optimal VIM, and R-VIM are established involving the Lagrange's multiplier; it has been demonstrated (and subsequently numerically certified) that solutions of VIM and Optimal VIM locally converge. While R-VIM solutions, with the aid of necessary and sufficient criteria, have been shown to converge globally for appropriate values of the unknown parameters involved in the general linear operator. We make sure that the convergence of the VIM solution depends critically on the optimal selection of the linear operator and, consequently, the best suited Lagrange's multiplier.</div><div>Our bird's eye view would be to differentiate between the contributions of the homotopy equation in homotopy-based approaches and Lagrange's multiplier-based correction functional of the VIM for the fastest convergence. Furthermore, in order to achieve fast convergence, we additionally examine if the optimal VIM's regulating parameter offers any further benefit above proposed R-VIM. We aim to apply R-VIM to problems with multiple solutions to ensure that all solutions are accurately identified. Our proposed method has been proven superior and more widely applicable through both theoretical and numerical illustrations. However, it also has its limits. Hence, we suggest the implementation of the multistep R-VIM, referred to as MR-VIM, as a means to solve chaotic dynamical systems.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"212 ","pages":"Pages 1-28"},"PeriodicalIF":2.2,"publicationDate":"2025-01-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143154831","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 unconditionally stable second-order scheme for Maxwell's equations in the Cole–Cole dispersive medium","authors":"Jingjing Xiao ,&nbsp;Desong Kong","doi":"10.1016/j.apnum.2025.01.009","DOIUrl":"10.1016/j.apnum.2025.01.009","url":null,"abstract":"<div><div>Coupling the averaged L1 scheme and the Crank–Nicolson scheme for temporal derivatives, we study Maxwell's equations in the Cole–Cole dispersive medium. A rigorous analysis is carried out to show that the proposed scheme is unconditionally stable and has a second-order convergence in time for sufficiently smooth solutions. We avoid using the mathematical induction method, which simplifies the analysis than the existing schemes. A fully discrete scheme with a finite difference method at Yee's grid is proposed. Numerical experiments are carefully designed to illustrate our theoretical analysis.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"211 ","pages":"Pages 211-227"},"PeriodicalIF":2.2,"publicationDate":"2025-01-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143141931","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":"Optimal L2 error estimates of two structure-preserving finite element methods for Schrödinger-Boussinesq equations","authors":"Houchao Zhang,&nbsp;Junjun Wang,&nbsp;Xueran Gong","doi":"10.1016/j.apnum.2025.01.012","DOIUrl":"10.1016/j.apnum.2025.01.012","url":null,"abstract":"<div><div>This paper is concerned with the construction and analysis of two structure-preserving finite element approximation schemes for solving one and two dimensional coupled Schrödinger-Boussinesq (SBq) equations. Firstly, two finite element approximation schemes are developed and the total mass and energy preserving properties of the proposed schemes are demonstrated in the discrete sense. Secondly, by use of the innovative cut-off function method, an auxiliary fully-discrete system is established, which leads to the unique solvability with the Brouwer's fixed-pointed theorem and the optimal <span><math><msup><mrow><mi>L</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span> error estimates for the proposed nonlinear scheme. Finally, linearized iterative algorithms are showed rigorously and extensive numerical results are presented to check the theoretical analysis of the structure-preserving numerical methods.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"211 ","pages":"Pages 193-210"},"PeriodicalIF":2.2,"publicationDate":"2025-01-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143141887","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":"Gradient descent-based parameter-free methods for solving coupled matrix equations and studying an application in dynamical systems","authors":"Akbar Shirilord,&nbsp;Mehdi Dehghan","doi":"10.1016/j.apnum.2025.01.011","DOIUrl":"10.1016/j.apnum.2025.01.011","url":null,"abstract":"<div><div>This paper explores advanced gradient descent-based parameter-free methods for solving coupled matrix equations and examines their applications in dynamical systems. We focus on the coupled matrix equations<span><span><span><math><mrow><mo>{</mo><mtable><mtr><mtd><mi>A</mi><mi>X</mi><mo>+</mo><mi>Y</mi><mi>B</mi><mo>=</mo><mi>C</mi><mo>,</mo></mtd></mtr><mtr><mtd><mi>D</mi><mi>X</mi><mo>+</mo><mi>Y</mi><mi>E</mi><mo>=</mo><mi>F</mi><mo>,</mo></mtd></mtr></mtable></mrow></math></span></span></span> where <span><math><mi>A</mi><mo>,</mo><mi>D</mi><mo>∈</mo><msup><mrow><mi>R</mi></mrow><mrow><mi>m</mi><mo>×</mo><mi>m</mi></mrow></msup><mo>,</mo><mi>B</mi><mo>,</mo><mi>E</mi><mo>∈</mo><msup><mrow><mi>R</mi></mrow><mrow><mi>n</mi><mo>×</mo><mi>n</mi></mrow></msup><mo>,</mo><mi>C</mi><mo>,</mo><mi>F</mi><mo>∈</mo><msup><mrow><mi>R</mi></mrow><mrow><mi>m</mi><mo>×</mo><mi>n</mi></mrow></msup></math></span> are given matrices, and <span><math><mi>X</mi><mo>,</mo><mi>Y</mi><mo>∈</mo><msup><mrow><mi>R</mi></mrow><mrow><mi>m</mi><mo>×</mo><mi>n</mi></mrow></msup></math></span> are the unknown matrices to be determined. We propose a novel gradient descent-based approach with parameter-free, enhancing convergence through an accelerated technique related to momentum methods. A comprehensive analysis of the convergence and characteristics of these methods is provided. Our convergence analysis demonstrates that if the spectrum of a block matrix constructed from the matrices <em>A</em>, <em>B</em>, <em>D</em>, and <em>E</em> is confined within a horizontal ellipse in the complex plane, centered at <span><math><mo>(</mo><mn>0</mn><mo>,</mo><mn>0</mn><mo>)</mo></math></span> with a major axis length of 3 and a minor axis length of 1, then the accelerated momentum method will converge to the exact solution of the discussed model. The numerical results indicate that proposed methods significantly improve efficiency, showing faster convergence and reduced computational time compared to traditional approaches. Additionally, we apply these methods to linear dynamic systems, demonstrating their effectiveness in real-world scenarios.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"212 ","pages":"Pages 29-59"},"PeriodicalIF":2.2,"publicationDate":"2025-01-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143155008","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 posteriori error analysis of a mixed finite element method for the stationary convective Brinkman–Forchheimer problem","authors":"Sergio Caucao ,&nbsp;Gabriel N. Gatica ,&nbsp;Luis F. Gatica","doi":"10.1016/j.apnum.2025.01.007","DOIUrl":"10.1016/j.apnum.2025.01.007","url":null,"abstract":"<div><div>We consider a Banach spaces-based mixed variational formulation that has been recently proposed for the nonlinear problem given by the stationary convective Brinkman–Forchheimer equations, and develop a reliable and efficient residual-based <em>a posteriori</em> error estimator for the 2D and 3D versions of the associated mixed finite element scheme. For the reliability analysis, we utilize the global inf-sup condition of the problem, combined with appropriate small data assumptions, a stable Helmholtz decomposition in nonstandard Banach spaces, and the local approximation properties of the Raviart–Thomas and Clément interpolants. In turn, inverse inequalities, the localization technique based on bubble functions in local <span><math><msup><mrow><mi>L</mi></mrow><mrow><mi>p</mi></mrow></msup></math></span>-spaces, and known results from previous works, are the main tools yielding the efficiency estimate. Finally, several numerical results confirming the theoretical properties of the estimator and illustrating the performance of the associated adaptive algorithm are reported. In particular, the case of flow through a 2D porous medium with fracture networks is considered.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"211 ","pages":"Pages 158-178"},"PeriodicalIF":2.2,"publicationDate":"2025-01-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143141885","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 robust mapping spectral method for elastic equations in curved fan-shaped domains","authors":"Yuling Guo ,&nbsp;Zhongqing Wang ,&nbsp;Chao Zhang","doi":"10.1016/j.apnum.2025.01.008","DOIUrl":"10.1016/j.apnum.2025.01.008","url":null,"abstract":"<div><div>In this paper, we introduce a robust mapping Legendre spectral-Galerkin method for addressing elastic problems in simply-connected, fan-shaped domains with curved boundaries. By employing a polar coordinate transformation, we map the fan-shaped domain onto a rectangle, which transforms the original elastic equation into a variable coefficient equation. We then develop a Legendre spectral-Galerkin scheme for this variable coefficient equation. Additionally, we demonstrate the optimal convergence of the numerical solution in the <span><math><msup><mrow><mi>H</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span>-norm as the Lamé coefficient <em>λ</em> remains bounded. Numerical examples illustrate the high accuracy and robustness of our method, even as <em>λ</em> approaches infinity.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"211 ","pages":"Pages 144-157"},"PeriodicalIF":2.2,"publicationDate":"2025-01-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143141888","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":"Variable-time-step weighted IMEX FEMs for nonlinear evolution equations","authors":"Meng Li ,&nbsp;Dan Wang ,&nbsp;Junjun Wang ,&nbsp;Xiaolong Zhao","doi":"10.1016/j.apnum.2025.01.005","DOIUrl":"10.1016/j.apnum.2025.01.005","url":null,"abstract":"<div><div>In this paper, the variable-time-step weighted implicit-explicit (IMEX) finite element methods (FEMs) are developed for some types of nonlinear real- or complex-valued evolution equations. Extensive research is conducted on the discrete orthogonal convolution (DOC) kernels and the discrete complementary convolution (DCC) kernels of the variable-time-step weighted IMEX scheme, elucidating their crucial properties in both real- and complex-valued scenarios. We prove that the scheme exhibits optimal convergence without any restrictions on the time-space step ratio. At last, several numerical examples are provided to demonstrate our theoretical results. With the weighted parameter <span><math><mi>θ</mi><mo>=</mo><mn>1</mn></math></span>, the scheme in this work can degenerate into a special case: variable-time-step two-step backward differentiation formula (BDF2) scheme, and the convergence analysis in this special case was introduced in Liao et al. (2020) <span><span>[15]</span></span> and Liao et al. (2021) <span><span>[29]</span></span>.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"211 ","pages":"Pages 123-143"},"PeriodicalIF":2.2,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143141889","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":"Analysis of average bound preserving time-implicit discretizations for convection-diffusion-reaction equation","authors":"Fengna Yan ,&nbsp;Yinhua Xia","doi":"10.1016/j.apnum.2025.01.006","DOIUrl":"10.1016/j.apnum.2025.01.006","url":null,"abstract":"<div><div>We propose a high-order average bound-preserving limiter for implicit backward differentiation formula (BDF) and local discontinuous Galerkin (LDG) discretizations applied to convection-diffusion-reaction equations. Our approach first imposes cell average bounds of the numerical solution using the Karush-Kuhn-Tucker (KKT) limiter and then enforces pointwise bounds with an explicit bound-preserving limiter. This method reduces the number of constraints compared to using only the KKT system to directly ensure pointwise bounds, resulting in a relatively small system of nonlinear equations to solve at each time step. We prove the unique solvability of the proposed average bound-preserving BDF-LDG discretizations. Furthermore, we establish the stability and optimal error estimates for the second-order average bound-preserving BDF2-LDG discretization. The unique solvability and stability are derived by transforming the KKT-limited cell average bounds-preserving LDG discretizations into a variational inequality. The error estimates are derived using the cell average bounds-preserving inequality constraints. Numerical results are presented to validate the accuracy and effectiveness of the proposed method in preserving the bounds.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"211 ","pages":"Pages 103-122"},"PeriodicalIF":2.2,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143141891","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":"Singular fractional integro-differential equations with non-constant coefficients","authors":"Kaido Lätt,&nbsp;Arvet Pedas,&nbsp;Hanna Britt Soots","doi":"10.1016/j.apnum.2025.01.004","DOIUrl":"10.1016/j.apnum.2025.01.004","url":null,"abstract":"<div><div>We investigate a class of singular fractional integro-differential equations with non-constant coefficients. After reformulating the original problem as a cordial Volterra integral equation we study the unique solvability of the underlying problem. We construct a numerical method based on collocation techniques to find an approximation to the solution of the original problem and analyse the convergence and the convergence order of the proposed method. Additionally, we present the results of some numerical experiments.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"211 ","pages":"Pages 179-192"},"PeriodicalIF":2.2,"publicationDate":"2025-01-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143142752","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}