{"title":"Deep learning-complex variable meshless method for inverse Poisson problems","authors":"Wenna He, Xiaofeng Liu, Heng Cheng","doi":"10.1016/j.enganabound.2025.106495","DOIUrl":"10.1016/j.enganabound.2025.106495","url":null,"abstract":"<div><div>This study proposes a deep learning-complex variable meshless method (DL-CVMM) framework that integrates deep neural networks (DNNs) with an improved complex variable element-free Galerkin (ICVEFG) method for solving inverse Poisson problem. The framework takes 2D coordinates as input and predicts source terms via forward propagation in DNNs. These predicted source terms are then incorporated into the ICVEFG discretization scheme to reconstruct the physical field. The inverse problem is formulated as an optimization problem by minimizing the empirical risk function in the problem domain between the reconstructed and observed values. This framework leverages DNNs for source term prediction, harnessing their generalization and learning capabilities, while employing the ICVEFG method for efficient physical field reconstruction. A key advantage of ICVEFG compared to the element-free Galerkin (EFG) method is its reduction in the number of unknown coefficients, reducing the matrix order in the shape function computations. Rigorous validations on three 2D Poisson inverse problem benchmarks demonstrate that the DL-CVMM framework achieves good convergence and computational accuracy.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"180 ","pages":"Article 106495"},"PeriodicalIF":4.1,"publicationDate":"2025-10-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145266832","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":"Domain decomposition peridynamic differential operator for solving singularly perturbed reaction–diffusion problems","authors":"Jun Wang , Chunlei Ruan , Feifei Zhou , Yun Chen","doi":"10.1016/j.enganabound.2025.106500","DOIUrl":"10.1016/j.enganabound.2025.106500","url":null,"abstract":"<div><div>A domain decomposition peridynamic differential operator (DD-PDDO) method for singularly perturbed reaction–diffusion problems is constructed. The method first uses the background mesh, such as Shishkin mesh here, to divide the computational domain into nonoverlapping subdomains; then uses PDDO to solve the problem in each subdomain independently; finally uses the continuity of the first-order derivative at the interface as a constraint mechanism to combine the solution of subdomains together. Present method overcomes the disadvantages of the weight function in the original PDDO for determining the interaction domain size at the interface transition points which have a large grid ratio. Numerical examples verify that the DD-PDDO method has stronger stability and higher accuracy than the PDDO method under conditions of large grid ratio.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"180 ","pages":"Article 106500"},"PeriodicalIF":4.1,"publicationDate":"2025-10-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145266831","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":"New general solutions and MFS methodology of Stokes equations","authors":"Chein-Shan Liu , Tai-Wen Hsu , Chia-Cheng Tsai","doi":"10.1016/j.enganabound.2025.106497","DOIUrl":"10.1016/j.enganabound.2025.106497","url":null,"abstract":"<div><div>Three-dimensional (3D) Stokes equations are reformulated to be the third-order partial differential equations, with four specific solutions being derived in Theorems 1–4. Then a third-order method of fundamental solutions (MFS) to solve the Stokes flow problems is developed. The Papkovich–Neuber solution is proven with an easier manner, which needs four 3D harmonic functions. The new solution with one 3D harmonic function and three 2D in-plane harmonic functions is more saving. Three important methods in Theorems 5–7 are proven to seek new solutions of the Stokes equations by means of a biharmonic potential function or a biharmonic vector; they are complete through a lengthy proof. An effort is made in Theorem 8 for providing a complete general solution, and the Slobodyanskii general solution is re-derived via a simple way; both of them are presented in terms of three harmonic functions in the Cartesian coordinates. The new solutions under a point force are employed to generate the Stokeslet as an application. Five fresh numerical methods are developed which automatically satisfy the incompressibility condition. A reduced MFS together with the Papkovich–Neuber solution are merged into an unsymmetric Stokeslet method. Two MFS with dipole source are also derived. To explore the efficiency and accuracy of the proposed numerical methods, some examples including a benchmark problem are tested.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"180 ","pages":"Article 106497"},"PeriodicalIF":4.1,"publicationDate":"2025-10-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145266834","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 meshless generalized finite difference method for the closed-loop geothermal system","authors":"Xiaotong Han , Lina Song , Cong Xie , Xiaoming He","doi":"10.1016/j.enganabound.2025.106482","DOIUrl":"10.1016/j.enganabound.2025.106482","url":null,"abstract":"<div><div>This study develops a meshless method to solve the closed-loop geothermal system. Unlike mesh-based methods that require additional interface terms, this method handles complex interface conditions in a simple and direct way by employing the Taylor series expansion in combination with the moving least squares approach. Moreover, to overcome the pressure oscillation problem in the model, an effective mixed boundary condition is introduced, without making any changes to the governing equations. Through the numerical cases, the method demonstrates the high accuracy and effectively simulates the target model.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"180 ","pages":"Article 106482"},"PeriodicalIF":4.1,"publicationDate":"2025-10-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145267627","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":"Employ a multigrid algorithm to solve the shape-transformation phase-field model","authors":"Jilong He","doi":"10.1016/j.enganabound.2025.106501","DOIUrl":"10.1016/j.enganabound.2025.106501","url":null,"abstract":"<div><div>This paper presents an original three-dimensional shape-transformation phase-field model solved by a highly efficient, optimally complex multigrid algorithm. The model extends the classical Allen–Cahn equation through the novel introduction of a static coupling term, <span><math><mrow><mi>α</mi><msup><mrow><mrow><mo>|</mo><msup><mrow><mi>ϕ</mi></mrow><mrow><mn>2</mn></mrow></msup><mo>−</mo><mn>1</mn><mo>|</mo></mrow></mrow><mrow><mn>2</mn></mrow></msup><mrow><mo>(</mo><mi>ψ</mi><mo>−</mo><mi>ϕ</mi><mo>)</mo></mrow></mrow></math></span>, which forces alignment to a target phase while preserving interface sharpness via a localization factor that vanishes at equilibrium (<span><math><mrow><mi>ϕ</mi><mo>=</mo><mo>±</mo><mn>1</mn></mrow></math></span>). We develop a customized V-cycle multigrid solver with four key components: (1) Gauss–Seidel relaxation for high-frequency error damping, (2) residual computation, (3) restriction operators, and (4) prolongation operators. For stability and efficiency, an operator-splitting technique decouples nonlinear diffusion and coupling terms, solved sequentially. Rigorous theoretical analysis proves <span><math><msup><mrow><mi>L</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span>-stability under a time-step constraint and establishes Q-linear convergence of Newton-multigrid iterations. Numerical experiments demonstrate <span><math><mrow><mo>></mo><mn>99</mn><mtext>%</mtext></mrow></math></span> similarity in topology-changing transformations, confirming the method’s robustness and computational superiority.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"180 ","pages":"Article 106501"},"PeriodicalIF":4.1,"publicationDate":"2025-10-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145267629","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}
Carlos Friedrich Loeffler, Leone Bernardo Florindo, Luciano de Oliveira Castro Lara, Lucas Silveira Campos
{"title":"Mitigating the requirement of poles in incompressible diffusive-advective heat transfer problems using the decomposition of derivatives and the multiple reciprocity techniques","authors":"Carlos Friedrich Loeffler, Leone Bernardo Florindo, Luciano de Oliveira Castro Lara, Lucas Silveira Campos","doi":"10.1016/j.enganabound.2025.106498","DOIUrl":"10.1016/j.enganabound.2025.106498","url":null,"abstract":"<div><div>The Derivatives Decomposition Technique is presented, applied in the context of the Direct Interpolation Boundary Element Method, to solve the incompressible diffusive-advective heat transfer equation. Despite achieving better accuracy, the primary objective is to develop an accurate model that uses the decomposition technique, aiming to mitigate the need for internal poles to solve the domain integral related to advective effects. Thus, the model is derived in conjunction with the well-known Multiple Reciprocity Method, which is first applied to the diffusive-advective equation. Within the contributions of this research, the significant role of the derivatives technique must be highlighted, as the velocity field is usually expressed in terms of Cartesian coordinates, which commonly requires radial approximations that compromise its precision. However, such components can be expressed in terms of standard and tangential components. Given that the derivatives of radial functions along the boundary are more precise than the normal derivatives, and that these normal derivatives can be included and solved directly in the matrix equation, the use of the Multiple Reciprocity Method becomes viable and advantageous, also in diffusive-advective models. This combined methodology yields a model that mitigates the requirement for internal poles for diffusive-advective problems, particularly at low Péclet numbers.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"180 ","pages":"Article 106498"},"PeriodicalIF":4.1,"publicationDate":"2025-10-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145266833","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-informed neural network for barrier option pricing in coupled financial quantitative system with varying interest rate and volatility","authors":"Yu Chen , Xing Lü , Hao Tian , Rui-Heng Li","doi":"10.1016/j.enganabound.2025.106457","DOIUrl":"10.1016/j.enganabound.2025.106457","url":null,"abstract":"<div><div>Accurate pricing of barrier options is essential for facilitating informed investment decisions, optimizing resource allocation, and promoting market stability. The dynamics of interest rate and volatility significantly influence the barrier option pricing. Traditional equations associated with constant parameters may fail to capture these complexities. In this paper, the underlying asset is assumed to follow an extended geometric Brownian motion incorporating varying interest rate and volatility, and then a coupled pricing system for the up-and-out call option is derived based on the Kolmogorov forward equation and backward equation, enabling the analysis of volatility fluctuations. Higher volatility indicates a greater level of risk, which may correspond to higher potential returns. The fusion framework, the physics-informed neural network (PINN), is introduced to solve this coupled system, consisting of two subnetworks: one dedicated to estimating the expected values of barrier option prices, and another for capturing the volatility surface associated with the option prices. Experimental results based on the closing price data of CSI 300ETF options show that PINN offers an effective and efficient framework for evaluating the prices of financial derivatives, achieving high precision and interpretability, even in cases where closed-form analytical solutions are unavailable.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"180 ","pages":"Article 106457"},"PeriodicalIF":4.1,"publicationDate":"2025-10-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145266830","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}
Zaizheng Ge , Haifeng Gao , Changjun Zheng , Haojie Lian , Qun Yan , Hongwei Zhou , Xiao Han , Zhe Xu , Chuanxing Bi
{"title":"A level set-based topology optimization method for the design of finite phononic crystals","authors":"Zaizheng Ge , Haifeng Gao , Changjun Zheng , Haojie Lian , Qun Yan , Hongwei Zhou , Xiao Han , Zhe Xu , Chuanxing Bi","doi":"10.1016/j.enganabound.2025.106492","DOIUrl":"10.1016/j.enganabound.2025.106492","url":null,"abstract":"<div><div>This study presents a level set-based topology optimization method for finite phononic crystals. To minimize the wave transmission on the boundary of the output domain, the adjoint method is employed to derive the topological derivatives. The boundary element method is then utilized to efficiently solve both the original and adjoint elastic dynamic problems, allowing for effective handling of boundary conditions. The boundaries of the material region are described by the zero-contour line of the level set function. Within a finite periodic unit, the optimization is performed, and numerical examples are provided for different frequency ranges. The results show that the proposed method can significantly reduce the wave transmission at a specific frequency and exhibit good vibration isolation performance. The vibration isolation effect can be increased by increasing the number of optimized layers, but the attenuation band gap is not necessarily enlarged. A parameter is also investigated for its influence on structural complexity, demonstrating that appropriate adjustment can balance geometric complexity and engineering feasibility. This method provides an effective means for the optimal design of phononic crystals and has the potential to be extended to a wider range of elastic metamaterial design problems.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"180 ","pages":"Article 106492"},"PeriodicalIF":4.1,"publicationDate":"2025-10-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145266829","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 development and analysis of a mesh-free technique for 1D and 2D time-fractional Kuramoto-Sivashinsky equation","authors":"Farzaneh Safari , Mojtaba Fardi","doi":"10.1016/j.enganabound.2025.106494","DOIUrl":"10.1016/j.enganabound.2025.106494","url":null,"abstract":"<div><div>The time-fractional Kuramoto–Sivashinsky equation (TFKSE) involving Caputo derivatives provides a model that describes pattern formation on flame fronts. We develop a numerical scheme to construct a time discretization for the TFKSE based on the <span><math><mrow><mi>L</mi><mn>1</mn></mrow></math></span> method for the Caputo fractional derivative. Furthermore, the proposed backward substitution method (BSM), which employs trigonometric basis functions, is used for numerical approximation. The first part of the BSM involves the approximation of the boundary data and the second part approximates the correction functions using basis functions. The numerical analysis of the scheme is presented, and the problem is discussed in one- and two-dimensional domains.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"180 ","pages":"Article 106494"},"PeriodicalIF":4.1,"publicationDate":"2025-10-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145267628","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}
A. Cheraghback , M. Botshekanan Dehkordi , Y. Kiani
{"title":"Dynamic snap-through phenomena in pseudoelastic shape memory alloy arch-beams","authors":"A. Cheraghback , M. Botshekanan Dehkordi , Y. Kiani","doi":"10.1016/j.enganabound.2025.106487","DOIUrl":"10.1016/j.enganabound.2025.106487","url":null,"abstract":"<div><div>This work deals with the dynamic snap-through phenomena in shape memory alloy (SMA) arch-beams taking into account the pseudoelastic effect of SMAs, a type of material nonlinearity, for the first time. For this aim, the Lagoudas formulation is used to model the pseudoelastic effect of SMAs. The SMA arch-beams are modeled employing the Timoshenko beam theory. By assuming the von Karman nonlinear strains, the governing equations of motion are derived while are coupled with the nonlinear phase transformation equations of SMAs. In this regard, the differential quadrature method (DQM) is employed to solve the nonlinear governing equations and the Newmark method is implemented to integrate these equations in the time domain. Meanwhile, the return mapping algorithm along with the Newton–Raphson technique is used to overcome the nonlinearities of the problem. In this work, the Budiansky criterion is employed to detect the dynamic buckling load of the SMA arch-beam. Under the dynamic buckling load, the arch-beams made of elastic materials exhibit multiple subsequent snap-throughs. A very interesting finding in this study is that, the material nonlinearity or pseudoelastic effect of SMAs prevents the subsequent dynamic snap-throughs in the SMA arch-beams, such that after the first snap phenomenon, the response of the system reaches a stable condition. The reason for this fact is the dissipation of energy due to the phase transformation in SMAs. In this regard, some novel results for pinned-pinned and fixed-fixed boundaries for different geometrical parameters have been presented.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"180 ","pages":"Article 106487"},"PeriodicalIF":4.1,"publicationDate":"2025-10-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145266828","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}