{"title":"The two-grid decoupled method for hybrid-dimensional fracture models based on the discontinuous Galerkin method","authors":"Shuangshuang Chen, Yuna Xu, Longchao Jin","doi":"10.1016/j.camwa.2025.07.029","DOIUrl":"10.1016/j.camwa.2025.07.029","url":null,"abstract":"<div><div>In this paper, we consider a kind of hybrid-dimensional fracture models, in which fractures are treated as <span><math><mo>(</mo><mi>n</mi><mo>−</mo><mn>1</mn><mo>)</mo></math></span>-dimensional objects in a <em>n</em>-dimensional porous medium, and the interaction between fractures and surrounding domain is taken into account. For the model with one single fracture, the discontinuous Galerkin (DG) method is firstly proposed, and the optimal order error estimate in the discrete <span><math><msup><mrow><mi>H</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span>-norm for the pressure is derived. Moreover, a two-grid decoupled method is considered, which uses a coarse grid rough approximations of the interface variables to decouple the mixed domain problem on a fine grid, and the optimal error estimate is also analyzed. The proposed DG method and two-grid decoupled method are then extended to a model with complex intersecting fractures, as well as the corresponding theoretical analysis. Numerical examples with both one single fracture and intersecting fractures are all presented to verity the accuracy of the considered methods.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"195 ","pages":"Pages 296-321"},"PeriodicalIF":2.9,"publicationDate":"2025-07-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144713795","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}
Patrick Ciarlet , Minh Hieu Do , Mario Gervais , François Madiot
{"title":"A posteriori error estimates for the DD+L2 jumps method on the neutron diffusion equations","authors":"Patrick Ciarlet , Minh Hieu Do , Mario Gervais , François Madiot","doi":"10.1016/j.camwa.2025.07.026","DOIUrl":"10.1016/j.camwa.2025.07.026","url":null,"abstract":"<div><div>We analyze <em>a posteriori</em> error estimates for the discretization of the neutron diffusion equations with a Domain Decomposition Method, the so-called DD+<span><math><msup><mrow><mi>L</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span> jumps method. We provide guaranteed and locally efficient estimators on a base block equation, the one-group neutron diffusion equation. Classically, one introduces a Lagrange multiplier to account for the jumps on the interface. This Lagrange multiplier is used for the reconstruction of the physical variables. Remarkably, no reconstruction of the Lagrange multiplier is needed to achieve the optimal <em>a posteriori</em> estimates.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"195 ","pages":"Pages 349-365"},"PeriodicalIF":2.9,"publicationDate":"2025-07-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144713796","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":"Equilibrated flux a posteriori error estimates and adaptivity for nonlinear and doubly degenerate elliptic problems","authors":"Elyes Ahmed , Saber Amdouni","doi":"10.1016/j.camwa.2025.07.019","DOIUrl":"10.1016/j.camwa.2025.07.019","url":null,"abstract":"<div><div>In this work, we establish a posteriori error analysis for a class of nonlinear and doubly degenerate elliptic equations, including the Stefan problem and both fast and slow diffusion in porous media. Our analysis leverages equilibrated flux reconstructions, yielding guaranteed and fully computable upper bounds on an energy-type error norm with local efficiency. These bounds are fully computable and robust, remaining unaffected by the degree of nonlinearity and degeneracy. The proposed estimators steer an adaptive solver that dynamically switches between nonlinear solvers to optimize iterations. The adaptive algorithm addresses discretization, regularization, and linearization errors. When Newton's method encounters convergence failures, the algorithm adaptively switches to the L-scheme solver. This solver precomputes the stabilization parameter <span><math><mi>L</mi><mo>></mo><mn>0</mn></math></span> in an offline phase to approximate the Jacobian. The efficiency of the adaptive algorithm is demonstrated through four prototypical examples, showcasing its effective error control and significant computational savings.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"195 ","pages":"Pages 239-264"},"PeriodicalIF":2.9,"publicationDate":"2025-07-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144703834","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}
Mohsen Jalalian , Manochehr Kazemi , Mohammad Esmael Samei
{"title":"Solving the second kind Volterra-Fredholm type of two-dimensional integral equations on non-rectangular domains via radial basis functions","authors":"Mohsen Jalalian , Manochehr Kazemi , Mohammad Esmael Samei","doi":"10.1016/j.camwa.2025.07.027","DOIUrl":"10.1016/j.camwa.2025.07.027","url":null,"abstract":"<div><div>This research, introduces a method to solve two-dimensional nonlinear Volterra-Fredholm integral equations with non-rectangular domains, numerically, based on radial basis functions. The method doesn't need a background mesh or cell structure in the domain. In the approach, all of the integrals are estimated using Gauss-Legendre quadrature formula. Error analysis and the rate of convergence of this method are also investigated. Numerical examples are included to demonstrate the validity and efficiency of this method.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"195 ","pages":"Pages 265-279"},"PeriodicalIF":2.9,"publicationDate":"2025-07-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144703833","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":"Performance of a vertex-centred block-coupled finite volume methodology for small-strain static elastoplasticity","authors":"Federico Mazzanti , Philip Cardiff","doi":"10.1016/j.camwa.2025.07.020","DOIUrl":"10.1016/j.camwa.2025.07.020","url":null,"abstract":"<div><div>This article assesses the performance of a vertex-centred block-coupled finite volume methodology for static small-strain elastoplasticity based on a Newton-Raphson algorithm. The proposed coupled solution algorithm is compared with vertex-centred and cell-centred segregated solution procedures in terms of accuracy, efficiency and robustness for quasi-static problems. This coupled methodology is then verified on four test cases: a patch test to check its ability to reproduce first-order polynomials; the method of manufactured solutions to check the order of convergence of displacement and stress; a 2-D perforated plate and a 3-D narrow T-member under elastoplastic conditions to examine the accuracy of the solutions on non-trivial test cases. The proposed methodology demonstrates high accuracy and a computational efficiency of up to 28× speedups over the standard cell-centred segregated approach.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"195 ","pages":"Pages 212-238"},"PeriodicalIF":2.9,"publicationDate":"2025-07-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144694977","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 PDE model for pricing European options with transaction costs under the 3/2 non-affine stochastic volatility model","authors":"Jianguo Tan , Jiling Cao","doi":"10.1016/j.camwa.2025.07.014","DOIUrl":"10.1016/j.camwa.2025.07.014","url":null,"abstract":"<div><div>In this paper, we study the problem of pricing European options with transaction costs under the 3/2 non-affine stochastic volatility model. First, we derive a nonlinear partial differential equations (PDE) model for pricing European options by using the expectation of transaction costs in a small time interval. It is worth to mention that the nonlinear PDE degenerates into the corresponding pricing PDE under the 3/2 stochastic volatility model when the transaction cost rate is set to zero. Then, we solve the nonlinear PDE numerically by using the finite difference method. Finally, we present numerical simulations and sensitivity analysis to illustrate both the consistency and the impact of transaction costs on option pricing.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"196 ","pages":"Pages 246-262"},"PeriodicalIF":2.9,"publicationDate":"2025-07-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144687010","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":"Denoising algorithm of the modified Gray-Scott model on non-uniform grids","authors":"Shanshan Ge , Jian Wang","doi":"10.1016/j.camwa.2025.07.016","DOIUrl":"10.1016/j.camwa.2025.07.016","url":null,"abstract":"<div><div>This paper introduces an denoising algorithm based on the Gray-Scott (GS) model, which adopts non-uniform grids and highly adaptive numerical strategy to achieve effective denoising for two-dimensional (2D) and three-dimensional (3D) models. This method adjusts the grid density based on regional characteristics: dense grids are used in the area of interest to improve local accuracy, and sparse grids are used in non-critical areas to enhance the overall computational efficiency. In order to enhance the robustness and shape preserving ability, we modify the original GS model by replacing <span><math><mo>(</mo><mi>F</mi><mo>+</mo><mi>k</mi><mo>)</mo><mi>v</mi></math></span> with <span><math><mo>(</mo><mi>F</mi><mo>+</mo><mi>k</mi><mo>)</mo><mo>(</mo><mi>v</mi><mo>−</mo><msub><mrow><mi>v</mi></mrow><mrow><mn>0</mn></mrow></msub><mo>)</mo></math></span>, which effectively relives the volume shrinkage and shape distortion problems, and achieves the retention of structural details and effective suppression of noise. Numerical experiments show that our algorithm can maintain key geometric features while improving the smoothness, and has excellent denoising performance and extensive applicability.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"195 ","pages":"Pages 191-211"},"PeriodicalIF":2.9,"publicationDate":"2025-07-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144686130","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":"Unconditionally energy gradient stable numerical scheme for Cahn–Hilliard equation with high-order degenerate mobility","authors":"Gyeonggyu Lee , Seunggyu Lee","doi":"10.1016/j.camwa.2025.07.012","DOIUrl":"10.1016/j.camwa.2025.07.012","url":null,"abstract":"<div><div>In this study, we consider the Cahn–Hilliard equation with high-order degenerate mobility, where the mobility function depends on the concentration field. As a Wasserstein gradient flow of the Ginzburg–Landau free energy, the system satisfies both mass preservation and energy dissipation requirements. A novel numerical scheme based on a linear stabilized splitting method proposed satisfying mass conservation, unique solvability, and unconditional energy gradient stability while drawing connections with the porous medium equation. Extensive numerical experiments validate the theoretical findings, including mass conservation, energy stability, and accuracy, in both space and time.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"196 ","pages":"Pages 263-287"},"PeriodicalIF":2.9,"publicationDate":"2025-07-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144687009","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}
N. Rasin , H. Hirol , A.R. Zainal Abidin , M.H. Mokhtaram , M.A. Mohd Noor , A.Y. Mohd Yassin
{"title":"NURBS-Divergence-Meshless (NDM) method","authors":"N. Rasin , H. Hirol , A.R. Zainal Abidin , M.H. Mokhtaram , M.A. Mohd Noor , A.Y. Mohd Yassin","doi":"10.1016/j.camwa.2025.07.011","DOIUrl":"10.1016/j.camwa.2025.07.011","url":null,"abstract":"<div><div>In this study, we present a new formulation that applies the divergence theorem within the NURBS parameter space to discretize governing equations using meshless techniques. We refer to this method as the NURBS-Divergence-Meshless (NDM) method. In this method, NURBS not only represents the domain exactly but also simplifies both the expression and the process of line integration, benefiting from the rectangular structure of the parameter space.</div><div>The primary advantage of NDM lies in its use of the divergence theorem, which simplifies the discretization of governing equations by ensuring that the line integral yields the required number of equations. It also enables the direct imposition of boundary conditions in interpolation schemes lacking the Kronecker delta property, such as moving least squares (MLS).</div><div>To illustrate the foundational concept, the discretization and solution of a simple scalar problem are presented. The equation is discretized using either the MLS or radial basis function (RBF) interpolation, with the divergence theorem providing the resulting algebraic equations. Convergence studies show that the non-overlapping NDM-MLS performs best, exhibiting all the anticipated favorable characteristics.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"196 ","pages":"Pages 151-171"},"PeriodicalIF":2.9,"publicationDate":"2025-07-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144679050","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 adaptive conforming mixed method with high-order elements for clustered eigenvalues of the Helmholtz transmission problem with index of refraction n(x)=1","authors":"Jinhua Feng, Shixi Wang, Hai Bi, Yidu Yang","doi":"10.1016/j.camwa.2025.07.017","DOIUrl":"10.1016/j.camwa.2025.07.017","url":null,"abstract":"<div><div>In this paper, we investigate the conforming mixed method for the Helmholtz transmission eigenvalue problem in cases where the scatterers have the same permeability in 2D or the same sound speed in 3D as the surrounding medium. Based on the work of Cakoni et al. (2009) <span><span>[22]</span></span> and Liu et al. (2023) <span><span>[23]</span></span>, we study the approximation of high-order elements for the clustered eigenvalues of the problem. We present an a priori error estimate, derive an a posteriori error estimator, and prove the reliability of the estimator for the approximate eigenvalues and the approximate eigenspaces. Numerical experiments show that our method is efficient and can compute the approximate eigenvalues with high precision.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"196 ","pages":"Pages 201-217"},"PeriodicalIF":2.9,"publicationDate":"2025-07-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144679088","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}