{"title":"Temporal splitting with Explicit-Implicit-Null method for nonlinear generalized multiscale subdiffusion equations","authors":"Zhengya Yang , Yating Wang , Wing Tat Leung","doi":"10.1016/j.cam.2025.116885","DOIUrl":"10.1016/j.cam.2025.116885","url":null,"abstract":"<div><div>In this paper, we propose an efficient temporal splitting scheme coupled with explicit-implicit-null (EIN) method for nonlinear time fractional multiscale diffusion equations with generalized kernel functions. To avoid the iterative method for nonlinear problems at each time step, the idea of EIN is to add and subtract a linear term to the original equation and apply implicit and explicit time marching scheme to the linear damping term and other terms, respectively. To handle the multiscale feature in the linear part, we further introduce a partially explicit temporal splitting scheme and construct two multiscale subspaces by means of the nonlocal multi-continua method to speed up the computation. The splitting schemes treat implicitly the subspace containing the high contrast and explicitly the other. We perform a rigorous stability analysis of the proposed algorithm and give the stability condition. Several numerical experiments are presented to confirm our theoretical results and demonstrate that the proposed algorithm achieves high accuracy while reducing computational cost for high contrast problems.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"473 ","pages":"Article 116885"},"PeriodicalIF":2.1,"publicationDate":"2025-07-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144570487","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 scalar/vector stabilized hybrid mixed method for Darcy–Stokes flows","authors":"Iury Igreja","doi":"10.1016/j.cam.2025.116879","DOIUrl":"10.1016/j.cam.2025.116879","url":null,"abstract":"<div><div>This work proposes a novel stabilized mixed hybrid finite element formulation for simulating coupled Stokes–Darcy flows. This approach achieves a significant reduction in global degrees of freedom compared to hybrid formulations that employ vector-valued Lagrange multipliers for the velocity field in both the free fluid and porous media regions. This advantage is realized by employing a vector multiplier for the Stokes velocity field and a scalar multiplier for the Darcy normal velocity component. This specific choice not only allows the natural simulation of heterogeneous porous media, but also naturally satisfies the interface conditions between the Stokes and Darcy domains. The linear system is constructed by condensing the velocity and pressure degrees of freedom at the element level, resulting in a problem solely defined by the scalar/vector multipliers. Furthermore, the global inf–sup stability of the method is supported by theoretical results from the literature. Numerical experiments in two and three dimensions evaluate the method’s accuracy, convergence rates, and robustness across homogeneous and heterogeneous porous media, as well as scenarios involving a spatially varying Beavers–Joseph–Saffman slip coefficient. The results confirm that, compared to hybrid formulations using vector multipliers in both the Stokes and Darcy subdomains, the proposed method achieves comparable or improved accuracy in homogeneous porous media and substantially reduces the number of degrees of freedom. In addition, the results demonstrate that the proposed method effectively handles heterogeneous porous media and maintains stability under spatially varying interface conditions, highlighting its suitability for complex coupled flow simulations.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"473 ","pages":"Article 116879"},"PeriodicalIF":2.1,"publicationDate":"2025-07-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144564111","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 dynamic mode decomposition-based Kalman filter for Bayesian inverse problem of nonlinear dynamical systems","authors":"Yuming Ba , Qiuqi Li , Zhiwen Zhang","doi":"10.1016/j.cam.2025.116880","DOIUrl":"10.1016/j.cam.2025.116880","url":null,"abstract":"<div><div>Ensemble Kalman filter (EnKF) method has been widely used in parameter estimation of the dynamic models, which needs to compute the forward model repeatedly. For nonlinear parameterized PDEs, constructing an accurate reduced order model is extremely challenging. To accelerate the posterior exploration efficiently, building surrogates of the forward models is necessary. In this paper, the dynamic mode decomposition (DMD) coupled with the weighted and interpolated nearest-neighbors (wiNN) algorithm is introduced to construct the surrogates for nonlinear dynamical systems. This extends the applicability of DMD to parameterized problems. Moreover, a low rank approximation of Kalman gain is used to EnKF, which can avoid the ensemble degenerate from the singularity of the covariance matrix. Finally, we apply the proposed method to nonlinear parameterized PDEs for the two-dimensional fluid flow and investigate their Bayesian inverse problems. The results are presented to show the applicability and efficiency of the proposed EnKF with DMD-wiNN method by taking account of parameters in nonlinear diffusion functions, nonlinear reaction functions and source functions.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"473 ","pages":"Article 116880"},"PeriodicalIF":2.1,"publicationDate":"2025-07-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144587803","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}
Dongmei Yu , Tiange Ma , Gehao Zhang , Cairong Chen , Deren Han
{"title":"Predefined-time robust gradient neural network for solving absolute value equations","authors":"Dongmei Yu , Tiange Ma , Gehao Zhang , Cairong Chen , Deren Han","doi":"10.1016/j.cam.2025.116887","DOIUrl":"10.1016/j.cam.2025.116887","url":null,"abstract":"<div><div>A predefined-time robust gradient neural network (PRGNN) for solving absolute value equations is investigated. PRGNN achieves predefined-time convergence and exhibits complete resilience against bounded vanishing or bounded non-vanishing noise. Compared with several existing continuous-time models, PRGNN can achieve better noise-tolerance performance under large constant noise. Numerical simulations demonstrate our claims.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"473 ","pages":"Article 116887"},"PeriodicalIF":2.1,"publicationDate":"2025-07-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144597528","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 novel fractional-order model and analysis of cancer-immune system interaction in an avascular environment with an efficient control mechanism","authors":"Hardik Joshi , Mehmet Yavuz","doi":"10.1016/j.cam.2025.116888","DOIUrl":"10.1016/j.cam.2025.116888","url":null,"abstract":"<div><div>The interaction between cancer and the immune system plays a vital role in tumor development, progression, and treatment, as well as in the advancement of effective cancer therapies. This paper explores the dynamic interaction between immune cells and cancerous cells, highlighting the critical role of the immune system in cancer prevention and treatment. A fractional-order model is developed to examine the influence of T-helper cells, cytotoxic T cells, B cells, and antibodies on cancerous cells. The model is thoroughly analyzed for feasibility and solution positivity. Additionally, the existence and uniqueness criteria, along with possible equilibrium points, are established. Conditions for both local and global asymptotic stability of equilibrium points are derived. Finally, numerical simulations validate the theoretical findings and emphasize key parameters affecting cancerous cell dynamics. The findings demonstrate that cytotoxic T cells and antibodies are vital in targeting and eliminating cancer cells, thereby strengthening the immune response. Additionally, the memory effect inherent in the fractional-order derivative profoundly shapes the system’s dynamics.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"473 ","pages":"Article 116888"},"PeriodicalIF":2.1,"publicationDate":"2025-07-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144579889","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":"Constructing orthogonal maximin distance uniform projection designs for computer experiments","authors":"A.M. Elsawah","doi":"10.1016/j.cam.2025.116902","DOIUrl":"10.1016/j.cam.2025.116902","url":null,"abstract":"<div><div>There is a significant need for computer experiments to study and model complex physical systems. In both computer and physical experiments, constructing experimental designs with good space-filling and column-orthogonality properties is crucial. While maximin distance designs and uniform designs ensure space-filling in full-dimensional spaces, they lack guarantees for low-dimensional projections. Uniform projection designs address this gap by ensuring space-filling properties in low-dimensional subspaces. Orthogonal designs enable efficient factor screening in Gaussian processes and ensure uncorrelated estimates of main effects in linear models. However, constructing such optimal designs remains challenging. A design that combines these advantages would outperform individual approaches. This paper fills this gap by proposing seven novel theoretical techniques for constructing orthogonal maximin distance uniform projection designs. The proposed designs demonstrate superior performance as the number of factors increases, making them particularly well-suited for surrogate modeling and linear trend estimation in high-dimensional Gaussian processes. Comparative studies show that the proposed techniques outperform existing methods.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"473 ","pages":"Article 116902"},"PeriodicalIF":2.1,"publicationDate":"2025-07-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144579891","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":"FEM-implicit/explicit surface mesh discrete schemes for 2D-3C time-dependent shallow shell problem","authors":"Rongfang Wu , Xiaoqin Shen , Ying Liu , Yumin Cheng","doi":"10.1016/j.cam.2025.116878","DOIUrl":"10.1016/j.cam.2025.116878","url":null,"abstract":"<div><div>In this work, we propose the finite element method (FEM) coupled implicit/explicit time difference schemes for the two-dimensional three-component (2D-3C) shallow shell surface equation, which features the second-order time variable and fourth-order space variable. Concretely, the higher-order spatial variable is discretized by employing a <span><math><msup><mrow><mi>C</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span>-conforming element; meanwhile, the existence and uniqueness of semi-discrete solutions for the shallow shell model are analyzed using eigenvalues, and the optimal error estimate is obtained. Then, the time variable is discretized using the Newmark format. We prove that the space–time fully discretization scheme achieves first-order convergence when the parameter <span><math><mrow><mi>γ</mi><mo>≠</mo><mn>0</mn><mo>.</mo><mn>5</mn></mrow></math></span>. In addition, we adopt the leapfrog format to improve the convergence, which achieves second-order accuracy within small time steps. Finally, numerical experiments are conducted to validate the stability and efficiency of numerical algorithms.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"473 ","pages":"Article 116878"},"PeriodicalIF":2.1,"publicationDate":"2025-07-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144587802","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 a priori error analysis of a thermoelastic problem with history dependence on the mechanical and thermal components","authors":"N. Bazarra , J.R. Fernández , R. Quintanilla","doi":"10.1016/j.cam.2025.116883","DOIUrl":"10.1016/j.cam.2025.116883","url":null,"abstract":"<div><div>Here, we provide an a priori error analysis of a thermoelastic problem, with the Moore–Gibson–Thompson (MGT) equation, where the history dependence is assumed on both the mechanical and thermal parts. An existence and uniqueness result, and the exponential stability, are recalled. Then, a fully discrete approximation is introduced by using the finite element method and the implicit Euler scheme, to approximate the spatial variable and the time derivatives, respectively. A discrete stability property is proved and a main a priori error estimates result is obtained. The linear convergence of the approximations is deduced under suitable regularity conditions. Finally, we perform some one- and two-dimensional simulations to show the accuracy of the algorithm, the exponential decay of the discrete energy and the behavior of the solution.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"473 ","pages":"Article 116883"},"PeriodicalIF":2.1,"publicationDate":"2025-07-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144570902","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":"Nonlinear quartic quasi-interpolating splines for piecewise smooth function approximation","authors":"Francesc Aràndiga , Paola Lamberti , Sara Remogna","doi":"10.1016/j.cam.2025.116890","DOIUrl":"10.1016/j.cam.2025.116890","url":null,"abstract":"<div><div>Quasi-interpolation based on spline approximation methods is used in numerous applications. A quartic quasi-interpolating spline is a piecewise polynomial of degree four satisfying <span><math><msup><mrow><mi>C</mi></mrow><mrow><mn>3</mn></mrow></msup></math></span> continuity and fifth-order approximation, if the function to be approximated is sufficiently smooth. However, if the function has jump discontinuities, we observe that the Gibbs phenomenon appears when approximating near discontinuities. In this paper, we present nonlinear modifications of such a spline, based on weighted essentially non-oscillatory (WENO) techniques to avoid this phenomenon near discontinuities and, at the same time, maintain the fifth-order accuracy in smooth regions. We also provide some numerical and graphical tests confirming the theoretical results.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"474 ","pages":"Article 116890"},"PeriodicalIF":2.1,"publicationDate":"2025-07-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144704809","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}
José M. Ramón , Juan Ruiz-Álvarez , Dionisio F. Yáñez
{"title":"Non-linear Partition of Unity method","authors":"José M. Ramón , Juan Ruiz-Álvarez , Dionisio F. Yáñez","doi":"10.1016/j.cam.2025.116891","DOIUrl":"10.1016/j.cam.2025.116891","url":null,"abstract":"<div><div>This paper introduces the Non-linear Partition of Unity Method (NL-PUM), a novel technique integrating Radial Basis Function interpolation and Weighted Essentially Non-Oscillatory algorithms. As far as we know, this is the first time that a non-linear PUM method is introduced in the literature. The main advance of this proposal is providing an algorithm that keeps the properties of the PUM at smooth zones, while introducing a non-linear modification close to the discontinuities to avoid oscillations. This is done automatically computing estimations of the smoothness of the data and replacing the PUM by Sheppard method when all the data is affected by a discontinuity. Thus, the computation of smoothness indicators and the use of compactly supported base functions ensure precision in regions affected by the presence of discontinuities. Error bounds are calculated and validate the effectiveness of the new method, showing improved interpolation capabilities at discontinuity regions as well as at smooth zones. A battery of experiments is presented to check the theoretical results provided.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"473 ","pages":"Article 116891"},"PeriodicalIF":2.1,"publicationDate":"2025-07-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144579892","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}