{"title":"Spatiotemporal dynamics deduced by nonlocal delay competition in a diffusive Lotka-Volterra population model","authors":"Xiaosong Tang , Jiaxin Shen , Xinchang Wang , Zhaoyun Zeng , Jingwen Zhu","doi":"10.1016/j.camwa.2025.07.008","DOIUrl":"10.1016/j.camwa.2025.07.008","url":null,"abstract":"<div><div>In this article, under the influence of nonlocal delay competition, we are devoted to investigating spatiotemporal dynamics of a diffusive Lotka-Volterra population model. According to the different values of some parameters, the original model may be changed into Lotka-Volterra predator-prey model, cooperative population model, or competition population model. Then, through the characteristic equation analysis, we find that nonlocal delay competition can deduce the presence of Turing-Hopf bifurcation when the original model is Lotka-Volterra cooperative population model or competition population model. However, when the original model is Lotka-Volterra predator-prey model, nonlocal delay competition cannot deduce the presence of Turing-Hopf bifurcation, but delay can deduce stability switches and the presence of Hopf bifurcation. Moreover, in the known literatures, to our knowledge, reaction-diffusion model with nonlocal delay competition has been investigated more rarely, which implies that our results in this paper are new. Finally, by presenting some numerical calculations and simulations, we obtain the rich results of stable spatially homogeneous periodic solutions, spatially steady state solutions, spatially inhomogeneous periodic solutions and stability switches deduced by delay.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"196 ","pages":"Pages 172-182"},"PeriodicalIF":2.9,"publicationDate":"2025-07-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144679086","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":"Physics-constrained deep kernel learning for inverse problems with noisy data","authors":"Zhenjie Tang, Li He","doi":"10.1016/j.camwa.2025.07.022","DOIUrl":"10.1016/j.camwa.2025.07.022","url":null,"abstract":"<div><div>We propose a novel physics-constrained deep kernel learning (PCDKL) to estimate physical parameters and learn forward solutions for problems described by partial differential equations (PDEs) and noisy data. In this framework, a Gaussian Process (GP) with a deep kernel is constructed to model the forward solution. The posterior function samples from the GP serve as surrogates for the PDE solution. These GP posterior samples are constrained by two likelihoods: one to fit the noisy observations and the other to enforce conformity with the governing equation. To efficiently and effectively infer the deep kernel and physical parameters, we develop a stochastic estimation algorithm based on the evidence lower bound (ELBO), which serves as a posterior regularization objective function. The effectiveness of the proposed PCDKL is demonstrated through a systematic comparison with a Bayesian physics-informed neural network (B-PINN), a state-of-the-art method for solving inverse problems in PDEs with noisy observations. Our experiments show that PCDKL not only achieves forward solutions with informative uncertainty estimates comparable to B-PINN, but also yields accurate estimates of physical parameters. These results suggest that PCDKL has significant potential for uncertainty quantification in forward solutions and accurate physical parameter estimation, making it valuable for practical applications.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"196 ","pages":"Pages 135-150"},"PeriodicalIF":2.9,"publicationDate":"2025-07-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144670241","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}
Cai Mingchao , Li Jingzhi , Li Ziliang , Liu Qiang
{"title":"An efficient iterative decoupling method for thermo-poroelasticity based on a four-field formulation","authors":"Cai Mingchao , Li Jingzhi , Li Ziliang , Liu Qiang","doi":"10.1016/j.camwa.2025.07.021","DOIUrl":"10.1016/j.camwa.2025.07.021","url":null,"abstract":"<div><div>This paper studies the thermo-poroelasticity model. By introducing an intermediate variable, we transform the original three-field model into a four-field model. Building upon this four-field model, we present both a coupled finite element method and a decoupled iterative finite element method. We prove the stability and optimal convergence of the coupled finite element method. Furthermore, we establish the convergence of the decoupled iterative method. This paper focuses primarily on analyzing the iterative decoupled algorithm. It demonstrates that the algorithm's convergence does not require any additional assumptions about physical parameters or stabilization parameters. Numerical results are provided to demonstrate the effectiveness and theoretical validity of these new methods.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"195 ","pages":"Pages 139-160"},"PeriodicalIF":2.9,"publicationDate":"2025-07-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144671079","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 family of even-order edge-oriented nonconforming finite elements in 2D with efficient locally conservative flux reconstruction","authors":"Gwanghyun Jo , Hyeokjoo Park","doi":"10.1016/j.camwa.2025.07.024","DOIUrl":"10.1016/j.camwa.2025.07.024","url":null,"abstract":"<div><div>In this paper, we propose a new family of edge-oriented even-order nonconforming finite elements in 2D. The proposed element has fewer degrees of freedom compared to the existing edge-oriented nonconforming elements, and preserves some attractive properties of the Crouzeix-Raviart element, such as the optimal approximation capability and the inf-sup stability. We also present an efficient <span><math><mi>H</mi><mo>(</mo><mi>div</mi><mo>)</mo></math></span>-conforming flux reconstruction for the finite element discretization by the proposed element, which is possible due to its edge-oriented degrees of freedom.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"196 ","pages":"Pages 127-134"},"PeriodicalIF":2.9,"publicationDate":"2025-07-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144672318","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":"Vertex-centered control-volume mimetic finite difference methods","authors":"Rainer Helmig , Martin Schneider , Ivan Yotov","doi":"10.1016/j.camwa.2025.07.018","DOIUrl":"10.1016/j.camwa.2025.07.018","url":null,"abstract":"<div><div>We develop a new class of vertex-centered control-volume mimetic finite difference methods on polytopal meshes for second order elliptic equations. The schemes are based on a mixed velocity-pressure formulation. The pressure is constant on dual mesh control-volumes constructed around the primary mesh vertices. The normal velocity is constant on the faces of the control-volumes, resulting in local mass conservation over the control-volumes. We consider both symmetric velocity integration rules constructed over the control-volumes, as well as non-symmetric quadrature rules constructed over sub-volumes obtained by the intersection of primary and dual elements. The latter choice allows for explicit gradient construction and local multipoint flux elimination within the primary elements, resulting in a positive definite vertex-centered pressure system. On simplicial, quadrilateral or hexahedral meshes, these local flux methods are closely related, and in some cases equivalent, to the classical vertex-centered control-volume finite element methods based on piecewise polynomial finite element basis functions for the pressure. The mimetic finite difference framework is utilized to analyze the well posedness and accuracy of the proposed methods. We establish first order convergence for the pressure and the velocity in the discrete mimetic norms, as well as second order pressure superconvergence in the case of symmetric quadrature rules. A series of numerical experiments illustrates the convergence properties of the methods on problems with varying degree of anisotropy, heterogeneity, and grid complexity in two and three dimensions.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"196 ","pages":"Pages 104-126"},"PeriodicalIF":2.9,"publicationDate":"2025-07-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144672322","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 weak Galerkin finite element method for the two-dimensional Burgers' equation","authors":"Xiaoxiao Chen, Yanli Chen","doi":"10.1016/j.camwa.2025.07.010","DOIUrl":"10.1016/j.camwa.2025.07.010","url":null,"abstract":"<div><div>We propose a weak Galerkin finite element method for solving the two-dimensional Burgers' equation. Based on a special weak variational formulation, both semi-discrete and fully-discrete numerical schemes are established and analyzed. Further, we prove the existence and uniqueness of the fully discrete solution and derive the optimal error estimates in the discrete <span><math><msup><mrow><mi>H</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span> norm and the <span><math><msup><mrow><mi>L</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span> norm, respectively. Finally, some numerical examples are presented to support our theoretical analysis.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"196 ","pages":"Pages 84-103"},"PeriodicalIF":2.9,"publicationDate":"2025-07-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144654978","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}
Anouk Zandbergen , Tycho van Noorden , Alexander Heinlein
{"title":"Improving pseudo-time stepping convergence for CFD simulations with neural networks","authors":"Anouk Zandbergen , Tycho van Noorden , Alexander Heinlein","doi":"10.1016/j.camwa.2025.07.006","DOIUrl":"10.1016/j.camwa.2025.07.006","url":null,"abstract":"<div><div>Computational fluid dynamics (CFD) simulations of viscous fluids described by the stationary Navier–Stokes equations are considered. Depending on the Reynolds number of the flow, the Navier–Stokes equations may exhibit a highly nonlinear behavior. The system of nonlinear equations resulting from the discretization of the Navier–Stokes equations can be solved using nonlinear iteration methods, such as Newton's method. However, fast quadratic convergence is typically only obtained in a local neighborhood of the solution, and for many configurations, the classical Newton iteration does not converge at all. In such cases, so-called globalization techniques may help to improve convergence.</div><div>In this paper, pseudo-time stepping (also known as pseudo-transient continuation) is employed in order to improve nonlinear convergence. The classical algorithm is enhanced by a neural network model that is trained to predict a local pseudo-time step. Generalization of the novel approach is facilitated by predicting the local pseudo-time step separately on each element using only local information on a patch of adjacent elements as input. Numerical results for standard benchmark problems, including flow over a backward facing step geometry and Couette flow, show the performance of the machine learning-enhanced globalization approach; as the software for the simulations, the CFD Module of COMSOL Multiphysics® is employed.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"196 ","pages":"Pages 64-83"},"PeriodicalIF":2.9,"publicationDate":"2025-07-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144632942","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}
Jingjie Cheng , Qing Xia , Binhu Xia , Junseok Kim , Yibao Li
{"title":"Decoupled, efficient and structure-preserving numerical scheme for a non-isothermal phase field sintering model","authors":"Jingjie Cheng , Qing Xia , Binhu Xia , Junseok Kim , Yibao Li","doi":"10.1016/j.camwa.2025.07.007","DOIUrl":"10.1016/j.camwa.2025.07.007","url":null,"abstract":"<div><div>In order to investigate the impact of temperature in the selective laser sintering industrial process, we focus on developing the non-isothermal phase-field sintering model which combines the phase field equations with the thermodynamic framework. We employ a Lagrange multiplier method to handle nonlinear terms and a second-order accurate scheme is carried out using the backward differentiation formula framework. The system is transformed into linear quadratic harmonic equations. We only need to solve one nonlinear equation by Newton iteration in order to update the Lagrange multiplier. The discrete scheme is rigorously proved to follow the second law of thermodynamics. Numerical simulations are conducted to demonstrate the deformation of grains in the system under the influence of thermal driving, as well validate the stability, accuracy, and efficiency of our method.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"196 ","pages":"Pages 49-63"},"PeriodicalIF":2.9,"publicationDate":"2025-07-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144632941","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}
Yuan Xu , Siu-Long Lei , Hai-Wei Sun , Zhi-Wei Fang
{"title":"A preconditioned fast finite volume method for two-dimensional conservative space-fractional diffusion equations on non-uniform meshes","authors":"Yuan Xu , Siu-Long Lei , Hai-Wei Sun , Zhi-Wei Fang","doi":"10.1016/j.camwa.2025.07.005","DOIUrl":"10.1016/j.camwa.2025.07.005","url":null,"abstract":"<div><div>We develop a finite volume (FV) method to solve two-dimensional (2D) conservative space-fractional diffusion equations (SFDEs) on non-uniform meshes via a sum-of-exponentials (SOE) technique. The SOE technique is used to approximate the spatial kernel, resulting in a numerically stable scheme. To further improve efficiency, a preconditioned fast Krylov subspace iterative method is exploited to obtain the numerical solution. The matrix-vector multiplication can be performed in <span><math><mi>O</mi><mo>(</mo><mi>M</mi><msup><mrow><mi>log</mi></mrow><mrow><mn>2</mn></mrow></msup><mo></mo><mi>M</mi><mo>)</mo></math></span> operations for a matrix of size <em>M</em>. Numerical experiments confirm the effectiveness of the proposed algorithm in terms of the computational time and the number of iterations.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"196 ","pages":"Pages 30-48"},"PeriodicalIF":2.9,"publicationDate":"2025-07-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144614815","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}
Magdalena Pabisz , Askold Vilkha , Dominika Ciupek , Maciej Woźniak , Maciej Paszyński
{"title":"Concurrency of the exponential integrators for the glioblastoma brain tumor simulations","authors":"Magdalena Pabisz , Askold Vilkha , Dominika Ciupek , Maciej Woźniak , Maciej Paszyński","doi":"10.1016/j.camwa.2025.07.003","DOIUrl":"10.1016/j.camwa.2025.07.003","url":null,"abstract":"<div><div>We employ the trace theory as a tool to investigate the concurrency of the novel implementation of the exponential integrators method using the multinomial theorem. The exponential integrators code is applied for the glioblastoma brain tumor simulations. The paper introduces a solid theoretical background for trace theory-based concurrency analysis and uses an important example of brain tumor simulations to validate its applicability. We perform our simulations on available public MRI scan databases. The spatial discretization uses the finite difference method. The temporal discretization uses the exponential integrators method. We employ the trace theory to identify basic undividable tasks and their dependency relation. Based on the coloring of the task dependency graph, we localize sets of tasks that can be executed fully in parallel. From the results of the trace theory-based concurrency analysis, we designed an efficient MATLAB code parallelized for the GPGPU. We test the implementation by running 1000 time steps of the glioblastoma brain tumor simulation over the three-dimensional brain data with 194 x 248 x 256 pixels. We measure the scalability of the code, reaching 32 times speedup, using 128 x 128 x 128 computational mesh that can be computed on RTX 2080 Super GPU card equipped with 3072 processing cores and 8 GB of memory.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"196 ","pages":"Pages 13-29"},"PeriodicalIF":2.9,"publicationDate":"2025-07-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144614814","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}