Eunsuh Kim , Heejae Kwon , Sungha Cho , Kyongmin Yeo , Minseok Choi
{"title":"Stabilize physics-informed neural networks for stiff differential equations: Re-spacing layer","authors":"Eunsuh Kim , Heejae Kwon , Sungha Cho , Kyongmin Yeo , Minseok Choi","doi":"10.1016/j.camwa.2025.09.014","DOIUrl":"10.1016/j.camwa.2025.09.014","url":null,"abstract":"<div><div>Approximating the solution of stiff differential equations, which exhibit abrupt changes in certain regions, using physics-informed neural networks (PINNs) is challenging. Typically, training PINNs involves using a larger number of samples concentrated around regions of rapid changes to resolve the sharp gradients. However, this strategy leads to data imbalance, resulting in slower convergence and reduced solution quality. Here, we propose Re-spacing layer (RS-layer) to mitigate these challenges. RS-layer is a pre-trained encoding layer designed to map the skewed distribution of sampling points onto a uniform distribution, maintaining the desirable statistical properties of the input data for effective PINN training. We demonstrate that RS-layer improves PINN training by regularizing the solution gradient in the transformed space. The efficacy of our method is validated through numerical experiments on one-dimensional singularly perturbed equations, the ROBER problem, and the Akzo Nobel problem. Our results show that RS-layer not only accelerates convergence, but also enhances accuracy.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"200 ","pages":"Pages 167-179"},"PeriodicalIF":2.5,"publicationDate":"2025-09-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145110183","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 local weak form of the generalized finite difference method (GFDM) with control volume in heat conduction problems","authors":"Tao Zhang , Xiaofeng Zhou","doi":"10.1016/j.camwa.2025.09.015","DOIUrl":"10.1016/j.camwa.2025.09.015","url":null,"abstract":"<div><div>An explicit local weak form of the Generalized Finite Difference Method (GFDM) is developed for solving heat conduction problems by incorporating the concept of control volumes from the Finite Volume Method (FVM). The control volume is introduced as the local integral domain in the local weak form, where the Heaviside function is employed as the weight function, similar to traditional FVM. However, instead of the cell-centered control volume used in FVM, a vertex-centered control volume is applied for integration. The trial function is constructed using a Generalized Finite Difference approximation based on a Taylor expansion. By applying the divergence theorem, the integral of the governing equation over the control volume is transformed into a boundary integral, thereby reducing the continuity requirements for both the trial function and thermal conductivity. Several numerical examples demonstrate the proposed method's accuracy, stability, and convergence. Additionally, it offers local integration scheme without requiring any background grid.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"199 ","pages":"Pages 203-224"},"PeriodicalIF":2.5,"publicationDate":"2025-09-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145109779","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}
Mourad Hrizi , Ravi Prakash , Antonio André Novotny
{"title":"Approximation of unknown sources in a time fractional PDE by the optimal ones and their reconstruction","authors":"Mourad Hrizi , Ravi Prakash , Antonio André Novotny","doi":"10.1016/j.camwa.2025.09.020","DOIUrl":"10.1016/j.camwa.2025.09.020","url":null,"abstract":"<div><div>In this paper, our focus is on studying a geometric inverse source problem that is governed by two-dimensional time-fractional subdiffusion. The problem involves determining the shape and location of the unknown source's geometrical support from boundary measurements of its associated potential. Firstly, we prove the uniqueness of the inverse problem. In the second phase, we propose a novel reconstruction method that utilizes the coupled complex boundary method (CCBM) to solve the identification problem. The main idea of this method is to approximate the overdetermined problem to a complex boundary value problem with a complex Robin boundary condition coupling the Dirichlet and Neumann boundary conditions. Next, we utilize the imaginary part of the solution throughout the domain to construct a shape cost function, which we then minimize with respect to ball-shaped sources by using a Newton-type topological derivative method to reconstruct the geometrical support of the unknown source.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"200 ","pages":"Pages 117-144"},"PeriodicalIF":2.5,"publicationDate":"2025-09-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145109835","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 bond-based nonlocal anisotropic diffusion model with variable matrix-valued coefficients and its asymptotically compatible meshfree discretization","authors":"Xiaofang Wang , Hao Tian","doi":"10.1016/j.camwa.2025.09.018","DOIUrl":"10.1016/j.camwa.2025.09.018","url":null,"abstract":"<div><div>We propose a bond-based nonlocal anisotropic diffusion model with variable matrix-valued coefficients. In previous studies, the non-ordinary state-based nonlocal diffusion model <span><span>[37]</span></span> has been used effectively for simulating anisotropic diffusion. However, it encounters challenges such as high computational costs, complexities in implementing boundary conditions, and numerical oscillation of zero-energy mode. In this paper, we propose a novel bond-based, nonlocal anisotropic diffusion model, and the key idea is that incorporates a nonlocal operator via a kernel function, integrating matrix-valued diffusion coefficients. The influence region of our model consists of two parts: an elliptical region determined by the variable diffusion coefficient at a material point and an irregular region shaped by the coefficient at neighboring points. Furthermore, we confirm the well-posedness of the proposed model and deduce various properties, such as weak convergence and mass conservation. For computational implementation, we introduce a meshfree method that is shown to be asymptotically compatible and relies on the quadrature rule, which is compatible with the proposed nonlocal diffusion model and can effectively solve the model. To evaluate the precision and efficiency of the model, we performed comprehensive numerical experiments in both two and three dimensions. We have also confirmed the discrete maximum principle through experimental validation simultaneously.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"199 ","pages":"Pages 148-166"},"PeriodicalIF":2.5,"publicationDate":"2025-09-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145109776","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":"hp-version discontinuous Galerkin time-stepping schemes for diffusive-viscous wave equation","authors":"Min Zhang , Zhaonan Dong , Wenjing Yan","doi":"10.1016/j.camwa.2025.09.021","DOIUrl":"10.1016/j.camwa.2025.09.021","url":null,"abstract":"<div><div>This work introduces a fully discrete scheme for the diffusion-viscous wave equation (DVWe) in the second-order formulation, combining <em>hp</em>-DG time-stepping schemes with conforming finite element methods (FEM). Two major theoretical contributions are presented: (1) <em>hp</em>-version a priori error estimates in both the energy-norm and DG-norm, which are optimal in the spatial mesh size <em>h</em>, temporal step size <em>τ</em>, and temporal polynomial order <em>q</em>, yet suboptimal by one order in the spatial polynomial order <em>p</em>. Furthermore, for solutions exhibiting weak singularities in time, exponential convergence in terms of the total number of temporal degrees of freedom is proven using the <em>hp</em>-refinement strategy. (2) An energy decay estimate that offers explicit bounds involving the model parameters, discretization parameters, and the Poincaré inequality constant. A series of numerical experiments are presented to validate the practical performance of the proposed approach.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"200 ","pages":"Pages 145-166"},"PeriodicalIF":2.5,"publicationDate":"2025-09-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145110184","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":"Hybrid active-passive control of multiphase composite shells through high-order Chebyshev frequency response modeling","authors":"Duy-Khuong Ly , Huy-Cuong Vu-Do , Trung-Hau Dang , T. Nguyen-Thoi","doi":"10.1016/j.camwa.2025.09.012","DOIUrl":"10.1016/j.camwa.2025.09.012","url":null,"abstract":"<div><div>This paper presents a high-efficiency approach for controlled frequency response analysis of laminated multi-phase composite shells, employing high-order Chebyshev finite elements. A key feature of the composite material used in this study is its three-phase composition, consisting of a polymer matrix, carbon fiber, and carbon nanotube nanofillers embedded within the matrix phase. The inclusion of carbon nanotube nanofillers significantly enhances the mechanical and dynamic properties of the overall material. The proposed approach addresses shear/membrane locking and the computational expense of discrete layerwise models, ensuring both efficiency and precision in analyzing complex composite shells. By using high-order shape functions derived from Chebyshev polynomials, this approach achieves rapid convergence without sacrificing accuracy. The developed numerical model captures the frequency response of laminated shells across varying geometric configurations, while incorporating the effects of two-layer control patches that introduce active or passive damping. Additionally, the model accounts for the layerwise effect of the composite, allowing for accurate prediction of structural behavior under excited loads. Numerical validation confirms the robustness and versatility of the high-order Chebyshev finite elements, effectively overcoming common computational challenges such as locking and spurious modes, while providing highly accurate and reliable frequency response predictions for practical engineering applications. This study highlights the potential of using advanced numerical techniques in combination with innovative composite materials for improved dynamic performance in aerospace, automotive, and civil engineering applications.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"199 ","pages":"Pages 184-202"},"PeriodicalIF":2.5,"publicationDate":"2025-09-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145109778","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}
Marina Matthaiou , Volker John , Marwa Zainelabdeen
{"title":"Bound-preserving physics-informed neural networks for steady-state convection-diffusion-reaction problems","authors":"Marina Matthaiou , Volker John , Marwa Zainelabdeen","doi":"10.1016/j.camwa.2025.09.009","DOIUrl":"10.1016/j.camwa.2025.09.009","url":null,"abstract":"<div><div>Numerical approximations of solutions of convection-diffusion-reaction problems should take only physically admissible values. Provided that bounds for the admissible values are known, this paper presents several approaches within physics-informed neural networks (PINNs) and <em>hp</em>-variational PINNs (<em>hp</em>-VPINNs) to preserve these bounds for convection-dominated problems. These approaches comprise the inclusion of the requirement for bound preservation in the cost functional, a simple cut-off strategy for the unphysical values, and two methods that enforce bound preservation via the activation function of the output layer of the neural network. Numerical simulations are performed for convection-dominated problems defined in two-dimensional domains. A variety of choices for several hyperparameters is explored. Enforcing bound preservation with the sine activation function in the output layer turned out to be superior to all other methods with respect to the accuracy of the computed solutions, and in particular, the results are much more accurate than those obtained with the standard PINNs and <em>hp</em>-VPINNs.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"199 ","pages":"Pages 167-183"},"PeriodicalIF":2.5,"publicationDate":"2025-09-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145109777","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":"Shape derivative of the Laplacian eigenvalue problem","authors":"Zhengfang Zhang , Lulu Guo , Xiangjing Gao , Weifeng Chen , Xiaoliang Cheng","doi":"10.1016/j.camwa.2025.09.013","DOIUrl":"10.1016/j.camwa.2025.09.013","url":null,"abstract":"<div><div>The Laplacian eigenvalue problem with two different densities is investigated. By the squeeze theorem, the shape derivative of the least eigenvalue on the interface of two different density subdomains is derived. As an application, the minimization of the least eigenvalue with area constraint is considered. The shape derivative of the objective functional is applied as the velocity for the level set method to involve the interface. The numerical results validate that the proposed method is effective to capture the final optimized distribution of two different densities.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"199 ","pages":"Pages 127-147"},"PeriodicalIF":2.5,"publicationDate":"2025-09-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145093733","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 averaged L1 ADI compact difference scheme for the three-dimensional time-fractional mobile/immobile transport equation with weakly singular solutions","authors":"Kai Liu, Haixiang Zhang, Xuehua Yang","doi":"10.1016/j.camwa.2025.09.019","DOIUrl":"10.1016/j.camwa.2025.09.019","url":null,"abstract":"<div><div>In this paper, a three-dimensional (3D) time-fractional mobile/immobile (MIM) transport equation, which incorporates the Caputo time-fractional derivative of order <span><math><mi>α</mi><mo>∈</mo><mo>(</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>)</mo></math></span>, is taken into consideration. The space derivatives are discretized using the compact finite difference approximation, and the Caputo time-fractional derivative is estimated by employing the averaged L1 formula. Combining with corresponding alternating direction implicit (ADI) algorithms, the overall computational cost is reduced significantly. Using the discrete energy analysis methodology, we demonstrate that the suggested method possesses temporal second-order convergence and spatial fourth-order convergence under the regularity assumption. Numerical experiments demonstrate that ADI techniques is effective in computing 3D problems.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"200 ","pages":"Pages 102-116"},"PeriodicalIF":2.5,"publicationDate":"2025-09-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145093845","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}
Dong Yan , Xin-Jie Huang , Guiyuan Ma , Xin-Jiang He
{"title":"Pricing American options with exogenous and endogenous transaction costs","authors":"Dong Yan , Xin-Jie Huang , Guiyuan Ma , Xin-Jiang He","doi":"10.1016/j.camwa.2025.09.008","DOIUrl":"10.1016/j.camwa.2025.09.008","url":null,"abstract":"<div><div>We study an American option pricing problem with liquidity risks and transaction fees. As endogenous transaction costs, liquidity risks of the underlying asset are modeled by a mean-reverting process. Transaction fees are exogenous transaction costs and are assumed to be proportional to the trading amount, with the long-run liquidity level depending on the proportional transaction costs rate. Two nonlinear partial differential equations are established to characterize the option values for the holder and the writer, respectively. To illustrate the impact of these transaction costs on option prices and optimal exercise prices, we apply the alternating direction implicit method to solve the linear complementarity problem numerically. Finally, we conduct model calibration from market data via maximum likelihood estimation, and find that our model incorporating liquidity risks outperforms the Leland model significantly.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"200 ","pages":"Pages 85-101"},"PeriodicalIF":2.5,"publicationDate":"2025-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145093734","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}