{"title":"Unconditionally maximum principle preserving scheme for the Allen–Cahn equation with a logarithmic free energy","authors":"Junxiang Yang , Sangkwon Kim , Junseok Kim","doi":"10.1016/j.camwa.2025.07.032","DOIUrl":"10.1016/j.camwa.2025.07.032","url":null,"abstract":"<div><div>We present a novel computational scheme that unconditionally satisfies the maximum principle to solve the Allen–Cahn (AC) equation with a logarithmic Flory–Huggins potential. The proposed scheme uses an operator splitting approach integrated with a frozen coefficient technique without using any stabilization term. Due to the use of the frozen coefficient method, we can apply a closed-form solution to the highly nonlinear term. This combination allows the development of a scheme that rigorously maintains the maximum principle throughout the computation. We provide a detailed analytical proof of the discrete maximum principle for the proposed scheme to ensure its theoretical robustness. To validate the high performance, we conduct computational experiments, which confirm the unconditional stability and demonstrate the method's effectiveness in preserving the maximum principle. These numerical results highlight the proposed scheme's potential as a reliable tool for accurately solving the AC equation.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"195 ","pages":"Pages 366-375"},"PeriodicalIF":2.5,"publicationDate":"2025-07-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144721592","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}
Adam Sky , Michael Neunteufel , Jack S. Hale , Andreas Zilian
{"title":"Formulae and transformations for simplicial tensorial finite elements via polytopal templates","authors":"Adam Sky , Michael Neunteufel , Jack S. Hale , Andreas Zilian","doi":"10.1016/j.camwa.2025.07.028","DOIUrl":"10.1016/j.camwa.2025.07.028","url":null,"abstract":"<div><div>We introduce a unified method for constructing the basis functions of a wide variety of partially continuous tensor-valued finite elements on simplices using polytopal templates. These finite element spaces are essential for achieving well-posed discretisations of mixed formulations of partial differential equations that involve tensor-valued functions, such as the Hellinger–Reissner formulation of linear elasticity. In our proposed polytopal template method, the basis functions are constructed from template tensors associated with the geometric polytopes (vertices, edges, faces etc.) of the reference simplex and any scalar-valued <span><math><msup><mrow><mi>H</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span>-conforming finite element space. From this starting point we can construct the Regge, Hellan–Herrmann–Johnson, Pechstein–Schöberl, Hu–Zhang, Hu–Ma–Sun and Gopalakrishnan–Lederer–Schöberl elements. Because the Hu–Zhang element and the Hu–Ma–Sun element cannot be mapped from the reference simplex to a physical simplex via standard double Piola mappings, we also demonstrate that the polytopal template tensors can be used to define a consistent mapping from a reference simplex even to a non-affine simplex in the physical mesh. Finally, we discuss the implications of element regularity with two numerical examples for the Reissner–Mindlin plate problem.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"195 ","pages":"Pages 322-348"},"PeriodicalIF":2.9,"publicationDate":"2025-07-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144713999","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":"Numerical analysis of a discontinuous Galerkin method for time-fractional Navier-Stokes-Fokker-Planck equations with weakly singular solutions","authors":"Dong Liu , Weihua Deng","doi":"10.1016/j.camwa.2025.07.031","DOIUrl":"10.1016/j.camwa.2025.07.031","url":null,"abstract":"<div><div>In this paper, a class of time-fractional Navier-Stokes-Fokker-Planck equations (TF-NSFPEs) describing position of particles with anomalous diffusion in viscous incompressible fluids are proposed. We develop the symmetric interior penalty discontinuous Galerkin (IPDG) method for TF-NSFPEs. The <em>L</em>1 method in the time on graded mesh is used for the reason that the solution of the time fractional Fokker-Planck equation (TFFPE) usually has a weak singularity near the initial time. The stability and the optimal error estimates of the IPDG semi-discrete scheme are proved by using the discrete fractional Grönwall inequality. Further, based on these results, the fully discrete optimal error estimates are obtained. Finally, some numerical experiments are performed to justify the effectiveness of the theoretical results.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"195 ","pages":"Pages 280-295"},"PeriodicalIF":2.9,"publicationDate":"2025-07-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144713797","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 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}