{"title":"Decoupling technology for systems of evolutionary equations","authors":"P.N. Vabishchevich","doi":"10.1016/j.camwa.2025.04.022","DOIUrl":"10.1016/j.camwa.2025.04.022","url":null,"abstract":"<div><div>Numerical methods of approximate solution of the Cauchy problem for coupled systems of evolution equations are considered. Separating simpler subproblems for individual components of the solution achieves simplification of the problem at a new level in time. The decoupling method, a significant approach to simplifying the problem, is based on the decomposition of the problem's operator matrix. The approximate solution is constructed based on the linear composition of solutions to auxiliary problems. The paper investigates decoupling variants based on extracting the diagonal part of the operator matrix and the lower and upper triangular submatrices. The study introduces a new decomposition approach, which involves splitting the operator matrix into rows and columns. The composition stage utilizes various variants of splitting schemes, showcasing the versatility of the approach. In additive operator-difference schemes, we can distinguish explicit-implicit schemes, factorized schemes for two-component splitting, and regularized schemes for general multi-component splitting. The study of stability of two- and three-level decoupling composition schemes is carried out using the theory of stability (correctness) of operator-difference schemes for finite-dimensional Hilbert spaces. The theoretical results of the decoupling technique for systems of evolution equations are illustrated on a test two-dimensional problem for a coupled system of two diffusion equations with inhomogeneous self- and cross-diffusion coefficients.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"191 ","pages":"Pages 105-128"},"PeriodicalIF":2.9,"publicationDate":"2025-04-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143882097","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 positive meshless finite difference scheme for scalar conservation laws with adaptive artificial viscosity driven by fault detection","authors":"Cesare Bracco , Oleg Davydov , Carlotta Giannelli , Alessandra Sestini","doi":"10.1016/j.camwa.2025.04.006","DOIUrl":"10.1016/j.camwa.2025.04.006","url":null,"abstract":"<div><div>We present a meshless finite difference method for multivariate scalar conservation laws that generates positive schemes satisfying a local maximum principle on irregular nodes and relies on artificial viscosity for shock capturing. Coupling two different numerical differentiation formulas and the adaptive selection of the sets of influence allows to meet a local CFL condition without any <em>a priori</em> time step restriction. The artificial viscosity term is chosen in an adaptive way by applying it only in the vicinity of the sharp features of the solution identified by an algorithm for fault detection on scattered data. Numerical tests demonstrate a robust performance of the method on irregular nodes and advantages of adaptive artificial viscosity. The accuracy of the obtained solutions is comparable to that for standard monotone methods available only on Cartesian grids.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"190 ","pages":"Pages 103-121"},"PeriodicalIF":2.9,"publicationDate":"2025-04-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143882997","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}
Abu Naser Sarker , Ronald D. Haynes , Michael Robertson
{"title":"A moving mesh method for pitting corrosion of heterogeneous materials","authors":"Abu Naser Sarker , Ronald D. Haynes , Michael Robertson","doi":"10.1016/j.camwa.2025.04.005","DOIUrl":"10.1016/j.camwa.2025.04.005","url":null,"abstract":"<div><div>An adaptive moving mesh or <em>r</em>-refinement method for the numerical approximation of pitting corrosion in heterogeneous materials is designed and applied to the problem of pitting corrosion in metals. The pitting corrosion is described by Laplace's equation with a moving boundary where the moving boundary problem is coupled with the partial differential equations describing the mesh movement. We show that the numerical approach is able to track evolving pit geometry for complicated materials with varying crystallography, corrosion-resistant inclusions, and material voids.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"191 ","pages":"Pages 48-59"},"PeriodicalIF":2.9,"publicationDate":"2025-04-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143874705","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 improved boundary knot method with fictitious points for solving high-order Helmholtz-type PDEs","authors":"L. Liu, L.L. Zhang, M. Lei, R.P. Niu","doi":"10.1016/j.camwa.2025.04.017","DOIUrl":"10.1016/j.camwa.2025.04.017","url":null,"abstract":"<div><div>An improved boundary knot method (IBKM) is proposed to enhance the performance of BKM in solving homogeneous high-order Helmholtz-type partial differential equations. Compared with the classical BKM where the sources are always placed on the physical boundary as collocation points, the new sources named fictitious points are now placed on multi-layer extended pseudo boundaries. This modification leads to higher accuracy without the loss of efficiency. Numerical examples are presented to demonstrate the superiority of the IBKM.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"191 ","pages":"Pages 36-47"},"PeriodicalIF":2.9,"publicationDate":"2025-04-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143869978","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":"Quasi-uniform unconditional superconvergent error estimates of FEMs for the time-dependent singularly perturbed Bi-wave problem modeling d-wave superconductors","authors":"Yanmi Wu , Dongyang Shi","doi":"10.1016/j.camwa.2025.04.018","DOIUrl":"10.1016/j.camwa.2025.04.018","url":null,"abstract":"<div><div>For the fourth order time-dependent singularly perturbed Bi-wave equation modeling <em>d</em>-wave superconductors, the implicit Backward Euler (BE) and Crank-Nicolson (CN) schemes of Galerkin finite element method (FEM) are studied by Bonner-Fox-Shmite element. Then the quasi-uniform and unconditional superconvergent error estimates of orders <span><math><mi>O</mi><mo>(</mo><msup><mrow><mi>h</mi></mrow><mrow><mn>3</mn></mrow></msup><mo>+</mo><mi>τ</mi><mo>)</mo></math></span> and <span><math><mi>O</mi><mo>(</mo><msup><mrow><mi>h</mi></mrow><mrow><mn>3</mn></mrow></msup><mo>+</mo><msup><mrow><mi>τ</mi></mrow><mrow><mn>2</mn></mrow></msup><mo>)</mo></math></span> (<em>h</em>, the spatial parameter, and <em>τ</em>, the time step) in the energy norm are derived respectively for the above schemes, which are independent of the negative powers of the perturbation parameter appearing in the model. Finally, some numerical results are provided to verify the theoretical analysis.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"191 ","pages":"Pages 24-35"},"PeriodicalIF":2.9,"publicationDate":"2025-04-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143869977","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 virtual element method for a convective Brinkman-Forchheimer problem coupled with a heat equation","authors":"Danilo Amigo , Felipe Lepe , Enrique Otárola , Gonzalo Rivera","doi":"10.1016/j.camwa.2025.04.015","DOIUrl":"10.1016/j.camwa.2025.04.015","url":null,"abstract":"<div><div>We develop a virtual element method to solve a convective Brinkman-Forchheimer problem coupled with a heat equation. This coupled model may allow for thermal diffusion and viscosity as a function of temperature. Under standard discretization assumptions and appropriate assumptions on the data, we prove the well posedness of the proposed numerical scheme. We also derive optimal error estimates under suitable regularity assumptions for the solution and appropriate assumptions on the data. We conclude with a series of numerical tests performed on different families of meshes that complement the theoretical findings.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"191 ","pages":"Pages 1-23"},"PeriodicalIF":2.9,"publicationDate":"2025-04-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143860614","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 spectral Galerkin method with higher order time discretization for the nonlinear stochastic Fisher's type equation driven by multiplicative noise","authors":"Huanrong Li , Rushuang Yang","doi":"10.1016/j.camwa.2025.04.016","DOIUrl":"10.1016/j.camwa.2025.04.016","url":null,"abstract":"<div><div>This paper primarily focuses on developing a high-order-in-time spectral Galerkin approximation method for nonlinear stochastic Fisher's type equations driven by multiplicative noise. For this reason, we first design an improved discretization scheme in time based on the Milstein method, and then propose a spectral Galerkin approximation method in space. We analyze the <span><math><msup><mrow><mi>H</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span> stability and <span><math><msup><mrow><mi>L</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span> stability estimations for numerical solutions of the proposed fully discrete spectral Galerkin approximation formulation under reasonable assumptions about the multiplicative noise function <span><math><mi>g</mi><mo>(</mo><mi>u</mi><mo>)</mo></math></span> and the nonlinear multiplicative function <span><math><mi>f</mi><mo>(</mo><mi>u</mi><mo>(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo>)</mo><mo>)</mo></math></span>. Additionally, we achieve nearly optimal error convergence orders in both space and time. Especially, the time convergence order almost reaches 1 under the <span><math><msup><mrow><mi>L</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span> norm. Finally, several numerical experiments are carried out for the stochastic Fisher's type models to validate all the theoretical results, and it can also be seen that the numerical results are consistent with the physical properties of the Fisher's equation.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"190 ","pages":"Pages 90-102"},"PeriodicalIF":2.9,"publicationDate":"2025-04-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143858303","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":"High order difference schemes for nonlinear Riesz space variable-order fractional diffusion equations","authors":"Qiu-Ya Wang","doi":"10.1016/j.camwa.2025.04.010","DOIUrl":"10.1016/j.camwa.2025.04.010","url":null,"abstract":"<div><div>This article aims at studying new finite difference methods for one-dimensional and two-dimensional nonlinear Riesz space variable-order (VO) fractional diffusion equations. In the presented model, fractional derivatives are defined in the Riemann-Liouville type. Based on 4-point weighted-shifted-Grünwald-difference (4WSGD) operators for Riemann-Liouville constant-order fractional derivatives, which have a free parameter and have at least third order accuracy, we derive variable-order 4WSGD operators for space Riesz variable-order fractional derivatives. In order that the fully discrete schemes exhibit robust stability and can handle the nonlinear term efficiently, we employ the implicit Euler (IE) method to discretize the time derivative, which leads to IE-4WSGD schemes. The stability and convergence properties of the IE-4WSGD schemes are analyzed theoretically. Additionally, a parameter selection strategy is derived for 4WSGD schemes and banded preconditioners are put forward to accelerate the GMRES methods for solving the discretization linear systems. Numerical results demonstrate the effectiveness of the proposed IE-4WSGD schemes and preconditioners.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"188 ","pages":"Pages 221-243"},"PeriodicalIF":2.9,"publicationDate":"2025-04-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143854843","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 transient free surface seepage flow using numerical manifold method","authors":"Zhen Jia, Hong Zheng","doi":"10.1016/j.camwa.2025.04.011","DOIUrl":"10.1016/j.camwa.2025.04.011","url":null,"abstract":"<div><div>In the analysis of transient seepage flow with free surfaces, not only the free surfaces but also the boundary conditions vary with time, introducing significant challenges to those traditional mesh-based numerical methods. Although the numerical manifold method (NMM) has shown great advantages in tracking time-independent free surface seepage flow due to its dual cover systems – the mathematical cover and the physical cover, in the analysis of transient free surface seepage flow it will encounter the inheritance issue of degrees of freedom between two consecutive time steps, which is still an open issue for all the partition of unity based methods such as the extended finite element method (XFEM) and the generalized finite element method (GFEM). It is shown in this study that the issue can be easily overcome if a different discretization order from the classical discretization order is adopted, <em>i.e.</em>, time discretization is carried out before to spatial discretization. By analyzing typical transient seepage examples, the positions of the transient free surfaces predicted by the proposed method are excellently consistent with analytical solutions or experimental results. At the same time, it also points out the errors and possible consequences of some literature concerning the handling of sudden drops in upstream water level. The results demonstrate that the proposed procedure not only effectively predicts the evolution of free surfaces but also accurately addresses transient seepage problems, including those with complex drainage systems.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"189 ","pages":"Pages 129-143"},"PeriodicalIF":2.9,"publicationDate":"2025-04-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143839089","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}
Yiyi Liu , Xueqing Teng , Xiaoqiang Yan , Hong Zhang
{"title":"A second-order, unconditionally invariant-set-preserving scheme for the FitzHugh-Nagumo equation","authors":"Yiyi Liu , Xueqing Teng , Xiaoqiang Yan , Hong Zhang","doi":"10.1016/j.camwa.2025.04.013","DOIUrl":"10.1016/j.camwa.2025.04.013","url":null,"abstract":"<div><div>In this paper, we present and analyze a second-order exponential time differencing Runge–Kutta (ETDRK2) scheme for the FitzHugh-Nagumo equation. Utilizing a second-order finite-difference space discretization, we derive the fully discrete numerical scheme by incorporating both the stabilization technique and the ETDRK2 scheme for temporal approximation. The smallest invariant set of the FitzHugh-Nagumo equation is presented. We demonstrate that the proposed scheme unconditionally preserves the invariant set without any time-step constraint. The convergence in both time and space is verified to achieve second-order accuracy. Numerical experiments are carried out to illustrate the efficiency, stability, and structure-preserving property of the proposed scheme.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"189 ","pages":"Pages 161-175"},"PeriodicalIF":2.9,"publicationDate":"2025-04-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143844179","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}