{"title":"Convergence properties of the radial basis function-finite difference method on specific stencils with applications in solving partial differential equations","authors":"Fazlollah Soleymani , Shengfeng Zhu","doi":"10.1016/j.enganabound.2024.106026","DOIUrl":"10.1016/j.enganabound.2024.106026","url":null,"abstract":"<div><div>We consider the problem of approximating a linear differential operator on several specific stencils using the radial basis function method in the finite difference scheme. We prove a linear convergence order on a non-equispaced five-point stencil. Then, we discuss how the convergence rate can be boosted up to the second-order on an equispaced stencil. Moreover, we show that including additional nearby nodes (six to twelve) in the stencil does not improve the convergence rate, thus increasing the computational load without enhancing convergence. To overcome this limitation, we propose a stencil that accelerates the convergence up to four using a nine-point stencil, unlike existing approaches which are based on thirteen-point equispaced stencils to achieve such an order of convergence. To support our findings, we conduct numerical experiments by solving Poisson equations and a parabolic problem.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"169 ","pages":"Article 106026"},"PeriodicalIF":4.2,"publicationDate":"2024-11-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142650900","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 radial basis function-compact finite difference method for Sobolev equation arising from fluid dynamics","authors":"Mohammad Ilati","doi":"10.1016/j.enganabound.2024.106020","DOIUrl":"10.1016/j.enganabound.2024.106020","url":null,"abstract":"<div><div>In this article, a new high-order, local meshless technique is presented for numerically solving multi-dimensional Sobolev equation arising from fluid dynamics. In the proposed method, Hermite radial basis function (RBF) interpolation technique is applied to approximate the operators of the model over local stencils. This leads to compact RBF generated finite difference (RBF-FD) formula, which provides a significant improvement in the accuracy and computational efficiency. In the first stage of the proposed method, the time discretization is performed by Crank–Nicolson finite difference scheme along with temporal Richardson extrapolation technique. In the second stage, the space dimension is discretized by applying the local radial basis function-compact finite difference (RBF-CFD) method. By performing some numerical simulations and comparing the results with existing methods, the high accuracy and computational efficiency of the proposed method are clearly demonstrated. The numerical results show that the presented method has fourth-order accuracy in both space and time dimensions. Finally, it can be concluded that the proposed method is a suitable alternative to the existing numerical techniques for the Sobolev model.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"169 ","pages":"Article 106020"},"PeriodicalIF":4.2,"publicationDate":"2024-11-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142650901","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 merging approach for hole identification with the NMM and WOA-BP cooperative neural network in heat conduction problem","authors":"X.L. Ji, H.H. Zhang, S.Y. Han","doi":"10.1016/j.enganabound.2024.106042","DOIUrl":"10.1016/j.enganabound.2024.106042","url":null,"abstract":"<div><div>Defect identification is an important issue in structural health monitoring. Herein, originated from inverse techniques, a merging approach is established by the numerical manifold method (NMM) and whale optimization algorithm-back propagation (WOA-BP) cooperative neural network to identify hole defects in heat conduction problems. On the one hand, the NMM can simulate varying hole configurations on a fixed mathematical cover, which eases the generation of “big data” for the training of neural network to a large extent. On the other hand, the WOA, a global optimization algorithm, is adopted to optimize the initial weights and thresholds of the BP neural network to alleviate its frequently encountered local optimum phenomenon. The boundary temperatures of sampling points by the NMM and the associated hole geometries are used for the learning of WOA-BP neural network, which is then applied to predict the hole defects. Numerical examples concerning the detection of circular/ elliptical holes demonstrate that the proposed method possesses higher accuracy and satisfying robustness in holes prediction compared with standard BP network under the same condition. The present work provides a convenient pathway and great potential in application of structural health monitoring.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"169 ","pages":"Article 106042"},"PeriodicalIF":4.2,"publicationDate":"2024-11-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142650898","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":"Fluid flow simulation with an H2-accelerated Boundary-Domain Integral Method","authors":"J. Tibaut , J. Ravnik , M. Schanz","doi":"10.1016/j.enganabound.2024.106015","DOIUrl":"10.1016/j.enganabound.2024.106015","url":null,"abstract":"<div><div>The development of new numerical methods for fluid flow simulations is challenging but such tools may help to understand flow problems better. Here, the Boundary-Domain Integral Method is applied to simulate laminar fluid flow governed by a dimensionless velocity–vorticity formulation of the Navier–Stokes equation. The Reynolds number is chosen in all examples small enough to ensure laminar flow conditions. The false transient approach is utilized to improve stability.</div><div>As all boundary element methods, the Boundary-Domain Integral Method has a quadratic complexity. Here, the <span><math><msup><mrow><mi>H</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span>-methodology is applied to obtain an almost linear complexity. This acceleration technique is not only applied to the boundary only part but more important to the domain related part of the formulation. The application of the <span><math><msup><mrow><mi>H</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span>-methodology does not allow to use the rigid body method to determine the singular integrals and the integral free term as done until now. It is shown how to apply the technique of Guigiani and Gigante to handle the strongly singular integrals in this application. Further, a parametric study shows the influence of the introduced approximation parameters. For this purpose the example of a lid driven cavity is utilized. The second example demonstrates the performance of the proposed method by simulating the Hagen–Poiseuille flow in a pipe. The third example considers the flow around a rigid cylinder to show the behavior of the method for an unstructured grid. All examples show that the proposed method results in an almost linear complexity as the mathematical analysis promisses.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"169 ","pages":"Article 106015"},"PeriodicalIF":4.2,"publicationDate":"2024-11-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142650902","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Underwater acoustic scattering of multiple elastic obstacles using T-matrix method","authors":"Yuzheng Yang , Qiang Gui , Yingbin Chai , Wei Li","doi":"10.1016/j.enganabound.2024.106028","DOIUrl":"10.1016/j.enganabound.2024.106028","url":null,"abstract":"<div><div>In this paper, the T-matrix method is applied to investigate the monostatic and bistatic far-field acoustic scattering patterns of underwater elastic multi-obstacles, which is a semi-analytical method and its results can be used to verify the accuracy of various numerical methods. The T-matrix formula for underwater multi-obstacle acoustic scattering is obtained by utilizing the addition theorem of the spherical harmonic function. Furthermore, an iterative algorithm is introduced to quickly solve the separation matrix in the addition theorem. The investigation into the far-field acoustic scattering characteristics of a pair of solid elastic spheres covers a full range of scattering angles, revealing that the resonance structure presented in the scattering spectra is attributed to the Rayleigh wave, specular reflections and Franz waves. Finally, experimental validation demonstrates good agreement between the numerical simulation results and experimental results.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"169 ","pages":"Article 106028"},"PeriodicalIF":4.2,"publicationDate":"2024-11-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142650998","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":"Geometrically nonlinear topology optimization of porous structures","authors":"Yongfeng Zheng, Rongna Cai, Jiawei He, Zihao Chen","doi":"10.1016/j.enganabound.2024.106014","DOIUrl":"10.1016/j.enganabound.2024.106014","url":null,"abstract":"<div><div>Porous structures are extensively used in engineering, and current designs of porous structures are constructed based on linear assumptions. In engineering, deformation cannot be ignored, so it is necessary to consider the effect of geometric nonlinearity in structural design. For the first time, this paper performs the geometrically nonlinear topology optimization of porous structures. This paper presents the theory of geometric nonlinear analysis, the bi-directional evolutionary method is used to search for the topological configurations of porous structures, the number of structural holes is determined by the number of periodicities. Furthermore, the optimization equation, sensitivity analysis, and optimization process are provided in detail. Lastly, four numerical examples are investigated to discuss the influence of geometric nonlinearity on the design of porous structures, such as comparisons between geometric nonlinear and linear design, the ability of geometric nonlinear design to resist cracks, changes in load amplitude and position and 3D porous designs. The conclusions drawn can provide strong reference for the design of high-performance porous structures.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"169 ","pages":"Article 106014"},"PeriodicalIF":4.2,"publicationDate":"2024-11-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142592833","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}
Aman Garg , Li Li , Weiguang Zheng , Mohamed-Ouejdi Belarbi , Roshan Raman
{"title":"Free vibration behaviour of bio-inspired helicoidal laminated composite panels of revolution under thermal conditions: Multi-output machine learning approach","authors":"Aman Garg , Li Li , Weiguang Zheng , Mohamed-Ouejdi Belarbi , Roshan Raman","doi":"10.1016/j.enganabound.2024.106024","DOIUrl":"10.1016/j.enganabound.2024.106024","url":null,"abstract":"<div><div>The present work aims to study the free vibration behaviour of bio-inspired helicoidal laminated composite spherical, toroid, and conical shell panels using a single-output Support Vector Machine (SVM) algorithm trained in the chassis of parabolic shear deformation theory under thermal conditions. Different helicoidal lamination schemes are adopted, such as Fibonacci, semi-circular, exponential, recursive, and linear helicoidal schemes. Temperature-dependent material properties are adopted. The effect of the geometry of the shell, temperature, and lamination scheme on the free vibration behaviour of spherical, toroid, and conical shell panels is studied. Also, the mode shapes are obtained using different multi-output SVM surrogate in which the displacements are obtained at different locations and are predicted to obtain the fundamental mode shape. The trained surrogate model can predict the values of fundamental frequency and mode shapes much faster than the parabolic shear deformation theory.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"169 ","pages":"Article 106024"},"PeriodicalIF":4.2,"publicationDate":"2024-11-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142592834","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}
Bin Li , Huayu Liu , Jian Liu , Miao Cui , Xiaowei Gao , Jun Lv
{"title":"A novel weak-form meshless method based on Lagrange interpolation for mechanical analysis of complex thin plates","authors":"Bin Li , Huayu Liu , Jian Liu , Miao Cui , Xiaowei Gao , Jun Lv","doi":"10.1016/j.enganabound.2024.106021","DOIUrl":"10.1016/j.enganabound.2024.106021","url":null,"abstract":"<div><div>In this paper, a novel weak-form meshless method, Galerkin Free Element Collocation Method (GFECM), is proposed for the mechanical analysis of thin plates. This method assimilates the benefits of establishing spatial partial derivatives by isoparametric elements and forming coefficient matrices node by node, which makes the calculation more convenient and stable. The pivotal aspect of GFECM is that the surrounding nodes can be freely chosen as collocation elements, which can adapt to irregular node distribution and suitable for complex models. Meanwhile, each collocation element is used as a Lagrange isoparametric element individually, which can easily construct high-order elements and improve the calculation accuracy, especially for high-order partial differential equations such as the Kirchhoff plate bending problem. In order to obtain the weak-form of the governing equation, the Galerkin form of the governing equation is constructed based on the virtual work principle and variational method. In addition, due to the Lagrange polynomials possessing the Kronecker delta property as shape functions, it can accurately impose boundary conditions compared with traditional meshless methods that use rational functions. Several numerical examples are proposed to verify the correctness and effectiveness of the proposed method in thin plate bending problems.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"169 ","pages":"Article 106021"},"PeriodicalIF":4.2,"publicationDate":"2024-11-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142592835","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":"Fourth-order phase field modelling of brittle fracture with strong form meshless method","authors":"Izaz Ali , Gašper Vuga , Boštjan Mavrič , Umut Hanoglu , Božidar Šarler","doi":"10.1016/j.enganabound.2024.106025","DOIUrl":"10.1016/j.enganabound.2024.106025","url":null,"abstract":"<div><div>This study aims to find a solution for crack propagation in 2D brittle elastic material using the local radial basis function collocation method. The staggered solution of the fourth-order phase field and mechanical model is structured with polyharmonic spline shape functions augmented with polynomials. Two benchmark tests are carried out to assess the performance of the method. First, a non-cracked square plate problem is solved under tensile loading to validate the implementation by comparing the numerical and analytical solutions. The analysis shows that the iterative process converges even with a large loading step, whereas the non-iterative process requires smaller steps for convergence to the analytical solution. In the second case, a single-edge cracked square plate subjected to tensile loading is solved, and the results show a good agreement with the reference solution. The effects of the incremental loading, length scale parameter, and mesh convergence for regular and scattered nodes are demonstrated. This study presents a pioneering attempt to solve the phase field crack propagation using a strong-form meshless method. The results underline the essential role of the represented method for an accurate and efficient solution to crack propagation. It also provides valuable insights for future research towards more sophisticated material models.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"169 ","pages":"Article 106025"},"PeriodicalIF":4.2,"publicationDate":"2024-11-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142592655","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Accelerated boundary integral analysis of energy eigenvalues for confined electron states in quantum semiconductor heterostructures","authors":"J.D. Phan , A.-V. Phan","doi":"10.1016/j.enganabound.2024.106012","DOIUrl":"10.1016/j.enganabound.2024.106012","url":null,"abstract":"<div><div>This paper presents a novel and efficient approach for the computation of energy eigenvalues in quantum semiconductor heterostructures. Accurate determination of the electronic states in these heterostructures is crucial for understanding their optical and electronic properties, making it a key challenge in semiconductor physics. The proposed method is based on utilizing series expansions of zero-order Bessel functions to numerically solve the Schrödinger equation using boundary integral method for bound electron states in a computationally efficient manner. To validate the proposed technique, the approach was applied to address issues previously explored by other research groups. The results clearly demonstrate the computational efficiency and high precision of the approach. Notably, the proposed technique significantly reduces the computational time compared to the conventional method for searching the energy eigenvalues in quantum structures.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"169 ","pages":"Article 106012"},"PeriodicalIF":4.2,"publicationDate":"2024-11-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142578666","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}