Engineering Analysis with Boundary Elements最新文献

{"title":"A linear fracture constitutive model for the two-dimensional finite-discrete element method (FDEM) and its parameters calibration procedure","authors":"Chengzeng Yan ,&nbsp;Du Han ,&nbsp;Hong Zheng ,&nbsp;Tie Wang ,&nbsp;Sajid Ali","doi":"10.1016/j.enganabound.2025.106228","DOIUrl":"10.1016/j.enganabound.2025.106228","url":null,"abstract":"<div><div>The parameters calibration of the original fracture constitutive model for the joint element in the finite-discrete element method (FDEM) is very complex. To simplify the parameter calibration procedure, we firstly propose a linear fracture constitutive model for the two-dimensional finite-discrete element method (FDEM). Based on the simplified constitutive model for the joint element, the proposed calibration procedure only involves adjusting two parameters (<em>r_o</em> and <em>r_s</em>), which represent the post-peak strain to pre-peak strain ratio and have a clear physical meaning. Other rock mechanical parameters can directly take from the experimental values. Then, the correctness of the proposed fracture constitutive model is verified by a series of rock numerical tests. Finally, the calibration procedure of the input parameters for the Jinping Baishan Group marble is provided by using the linear fracture constitutive model. The proposed constitutive model can significantly simplify the parameter calibration process and improve the computational efficiency, which enables FDEM better applied to practical engineering.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"176 ","pages":"Article 106228"},"PeriodicalIF":4.2,"publicationDate":"2025-04-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143767946","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":"Boundary element analysis for MHD Brinkman flow around circular cylinders inside a microchannel exhibiting wall roughness","authors":"Vishal Chhabra ,&nbsp;Chandra Shekhar Nishad ,&nbsp;Manoj Sahni ,&nbsp;Vineet Kumar Chaurasiya","doi":"10.1016/j.enganabound.2025.106235","DOIUrl":"10.1016/j.enganabound.2025.106235","url":null,"abstract":"<div><div>Inspired by the dynamics of blood flow around clots, emboli, and drug capsules in blood vessels, this study introduces a hydrodynamic model describing the steady, pressure-driven flow of a viscous, incompressible fluid past multiple equally sized circular cylinders within a rectangular microchannel. To create a permeable environment, the microchannel is embedded with a uniform, isotropic porous medium. To simulate roughness on the permeable walls, alternating Navier slip and no-slip boundary conditions are imposed, maintaining the same phase configuration. The system operates at a low Reynolds number and is influenced by an external magnetic field. The Brinkman equations, which govern the flow through the porous domain, are solved using the boundary element method (BEM). The hydrodynamics of the proposed model, with two instances of slip-periodicity, are studied extensively. Increasing Darcy number induces a more permeable porous medium, causing a significant reduction in flow resistance, particularly in regions farther from the channel boundaries. The Lorentz force, which is most effective at generating drag when perpendicular to the flow, becomes less impactful as the inclination angle of the magnetic field increases. The shear stress is minimized at the points at which Navier's slip and no-slip boundary conditions coincide. The proposed model enhances microfluidic systems for precise drug delivery, optimizes lab-on-chip devices for diagnostics, and improves fluid dynamics in porous systems like heat exchangers and filtration processes. It also supports the development of medical devices that better simulate natural blood flow, advancing the efficiency of artificial organs and implants.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"176 ","pages":"Article 106235"},"PeriodicalIF":4.2,"publicationDate":"2025-04-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143760311","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":"Implicit-splitting physics-driven particle relaxation for enhancement of Eulerian SPH","authors":"Bo Zhang ,&nbsp;Zhentong Wang ,&nbsp;Decheng Wan ,&nbsp;Xiangyu Hu","doi":"10.1016/j.enganabound.2025.106239","DOIUrl":"10.1016/j.enganabound.2025.106239","url":null,"abstract":"<div><div>Physics-driven particle relaxation, driven by either constant background pressure or the kernel gradient correction (KGC) matrix, has been proposed to generate isotropic and body-fitted particle distributions for complex geometries while ensuring zero-order consistency in smoothed particle hydrodynamics (SPH). However, this relaxation process often encounters challenges, such as a low decay rate of residuals and difficulties in achieving convergence to minor zero-order consistency errors, particularly for three-dimensional complex geometries and scenarios with a small ratio of smoothing length <span><math><mi>h</mi></math></span> to particle spacing <span><math><mrow><mi>Δ</mi><mi>x</mi></mrow></math></span>. This limitation hinders the development of Eulerian SPH (ESPH), as its accuracy depends on the initial particle configurations, while the smoothing length determines the computational cost. In this work, we introduce an implicit-splitting approach to improve the relaxation process, aiming to obtain isotropic particle distributions with negligible zero-order consistency errors at reduced <span><math><mrow><mi>h</mi><mo>/</mo><mi>Δ</mi><mi>x</mi></mrow></math></span> values. This approach provides the initial particle distribution for ESPH, enhancing its computational efficiency by reducing the number of neighboring particles. Extensive relaxation examples demonstrate that the proposed method significantly reduces the relaxation residual and achieves substantially smaller zero-order consistency errors. Subsequently, different incompressible numerical examples based on relaxed particles have been validated using Eulerian SPH. A value of <span><math><mrow><mi>h</mi><mo>/</mo><mi>Δ</mi><mi>x</mi></mrow></math></span> smaller than 1.0 was selected to reduce neighboring particles, and the reverse kernel gradient correction (RKGC) was adopted to ensure first-order consistency. The numerical results consistently show good accuracy and smoother contours than those obtained from the finite volume method (FVM) implemented in an SPH framework based on unstructured meshes.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"176 ","pages":"Article 106239"},"PeriodicalIF":4.2,"publicationDate":"2025-04-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143747891","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":"A novel heterogeneous CPU/GPGPU-accelerated 3D CDEM and its application to modeling deep roadway excavation","authors":"Junguang Huang ,&nbsp;Yiming Zhang ,&nbsp;Chun Feng ,&nbsp;Huanning Hu ,&nbsp;Minjie Wen","doi":"10.1016/j.enganabound.2025.106224","DOIUrl":"10.1016/j.enganabound.2025.106224","url":null,"abstract":"<div><div>To improve the accuracy and computational efficiency of the CDEM for deep coal mine roadway excavation modeling, this study proposes a heterogeneous CPU/GPGPU-accelerated solver that integrates a mixed continuous–discontinuous media algorithm. The solver employs an explicit time integration method combined with a modular approach for 3D tetrahedral solid finite elements and fracturable penalty springs, which model rock fracture behavior and the transition from continuum to discontinuum in rock masses. To maximize computational efficiency, the solver uses a hybrid CPU/GPGPU framework with SIMD parallel techniques, achieving up to 600-fold speedup on a single GPGPU. The solver’s accuracy is validated for both quasi-static and dynamic problems, and its scalability across different hardware accelerators is demonstrated. Modeling results from the 22nd mining area of the Quandian coal mine show significant shear deformation and crack evolution in the soft rock, particularly at the intersection of the roof slab and sidewall, where stress concentration and large deformation are most pronounced. These findings validate the efficiency and reliability of the proposed method for simulating and analyzing underground excavation processes.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"177 ","pages":"Article 106224"},"PeriodicalIF":4.2,"publicationDate":"2025-04-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143759161","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 simple and efficient three-dimensional spring element model for pore seepage problems","authors":"Jing Li ,&nbsp;Xinguang Zhu ,&nbsp;Chun Feng ,&nbsp;Minjie Wen ,&nbsp;Yiming Zhang","doi":"10.1016/j.enganabound.2025.106225","DOIUrl":"10.1016/j.enganabound.2025.106225","url":null,"abstract":"<div><div>This study introduces a novel spring element model for efficient simulation of nonlinear seepage in porous media. The model discretizes the simulation domain into tetrahedral elements and constructs orthogonal Three-dimensional permeability networks within each element, establishing a quantitative relationship between pipe flow and nodal pressure differences. By developing a mathematical model linking network flow to nodal pressure differences, the method enables precise allocation of pipe flow in the local coordinate system and accurate transformation to the global coordinate system, thereby determining nodal flow and velocity. The Three-Dimensional Seepage Spring Element Method (3D-SSEM) simplifies the element flow matrix in finite element analysis to three essential pipe permeability stiffness values, thereby reducing computational complexity. Coupled with parallel computing strategies, the algorithm achieves significant improvements in computational efficiency and memory usage. The method is validated through four numerical examples, demonstrating high efficiency and accuracy in solving saturated-unsaturated seepage problems. Compared with analytical solutions and other numerical methods, it exhibits superior convergence and reduced solution time while maintaining precision. Additionally, the method effectively simulates complex coupled processes in large-scale real-world environments, offering robust support for practical engineering design optimization.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"176 ","pages":"Article 106225"},"PeriodicalIF":4.2,"publicationDate":"2025-04-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143747888","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":"Weak voltage signal enhancement for accurate image reconstruction in the craniocerebral EIT","authors":"Yanyan Shi ,&nbsp;Hanxiao Dou ,&nbsp;Meng Wang ,&nbsp;Hao Su ,&nbsp;Feng Fu","doi":"10.1016/j.enganabound.2025.106237","DOIUrl":"10.1016/j.enganabound.2025.106237","url":null,"abstract":"<div><div>As a promising imaging technique, electrical impedance tomography (EIT) is used to reflect the conductivity distribution variation of human tissues. Different from the lung EIT, the application of the craniocerebral EIT is challenging. This is attributed to the fact that the skull with high resistivity greatly restricts the injected current from flowing into the brain tissue. Consequently, the measured voltage signal is very weak causing poor reconstructed images. To solve this problem, a new strategy based on a multi-layer convolutional neural network (CNN) is proposed for weak voltage signal enhancement. Voltage measurements from the three-layer head model and the single-layer head model perform as the input and the output of the network respectively. The trained network is supposed to enhance the voltage data of the three-layer head model. To test the performance of the proposed method, voltage data processed by the multi-layer CNN is compared with the single-layer voltage data. Besides, comparisons are also made in the case of noise interruption and when the skull thickness varies. The results demonstrate that the processed voltage data is almost consistent with the single-layer voltage data. Compared with the image reconstruction with the three-layer voltage data, there is a large improvement when using the proposed method.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"176 ","pages":"Article 106237"},"PeriodicalIF":4.2,"publicationDate":"2025-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143738460","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":"PINNs-MPF: A Physics-Informed Neural Network framework for Multi-Phase-Field simulation of interface dynamics","authors":"Seifallah Elfetni ,&nbsp;Reza Darvishi Kamachali","doi":"10.1016/j.enganabound.2025.106200","DOIUrl":"10.1016/j.enganabound.2025.106200","url":null,"abstract":"<div><div>We present PINNs-MPF framework, an application of Physics-Informed Neural Networks (PINNs) to handle Multi-Phase-Field (MPF) simulations of microstructure evolution. A combination of optimization techniques within PINNs and in direct relation to MPF method are extended and adapted. The numerical resolution is realized through a multi-variable time-series problem by using fully discrete resolution. Within each interval, space, time, and phases/grains are treated separately, constituting discrete subdomains. PINNs-MPF is equipped with an extended multi-networking (parallelization) concept to subdivide the simulation domain into multiple batches, with each batch associated with an independent NN trained to predict the solution. To ensure continuity across the spatio-temporal-phasic subdomains, a Master NN efficiently is to handle interactions among the multiple networks and facilitates the transfer of learning. A pyramidal training approach is proposed to the PINN community as a dual-impact method: to facilitate the initialization of training when dealing with multiple networks, and to unify the solution through an extended transfer of learning. Furthermore, a comprehensive approach is adopted to specifically focus the attention on the interfacial regions through a dynamic meshing process, significantly simplifying the tuning of hyper-parameters, serving as a key concept for addressing MPF problems using machine learning. We perform a set of systematic simulations that benchmark foundational aspects of MPF simulations, i.e., the curvature-driven dynamics of a diffuse interface, in the presence and absence of an external driving force, and the evolution and equilibrium of a triple junction. The proposed PINNs-MPF framework successfully reproduces benchmark tests with high fidelity and Mean Squared Error (MSE) loss values ranging from 10⁻⁶ to 10⁻⁴ compared to ground truth solutions.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"176 ","pages":"Article 106200"},"PeriodicalIF":4.2,"publicationDate":"2025-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143738461","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":"Eigensolutions to functionally graded material plates covered with piezoelectric surface layers","authors":"Pengchong Zhang ,&nbsp;Wentao Wang ,&nbsp;Haohao Xu ,&nbsp;Jia Peng ,&nbsp;Zixuan Zhou","doi":"10.1016/j.enganabound.2025.106240","DOIUrl":"10.1016/j.enganabound.2025.106240","url":null,"abstract":"<div><div>The transverse free vibration analysis of composite intelligent plates constituted by the functionally graded substrate and full size surface-attached piezoelectric laminae is conducted by means of the scaled boundary finite element method (SBFEM) in association with the precise integration algorithm (PIA). It is needful to point out that material coefficients of the functionally graded host layer are altered across the single direction or both in-plane coordinate axes and changed as any form of polynomial or non-polynomial formulae. In the suggested methodology, only the outer surface parallel with the <em>x</em>-O-<em>y</em> coordinate plane is required to be meshed with two dimensional spectral elements. Moreover, three and four degrees of freedom for non-homogeneous and piezoelectric laminae respectively are adopted, which is convenient to decrease the calculation expense and increase the computational efficiency. Supported by the scaled boundary coordinate system <em>z</em>-<em>η</em>-<em>ζ</em>, the governing equation of the hybrid plate is expressed as an ordinary differential one. Aided by the highly accurate PIA, the stiffness matrix of each layer from the analytical matrix exponential function can be acquired. Eigensolutions to laminated coupling plates are derived from the eigenvalue equation composed of stiffness and mass matrices. From tabular and graphical comparisons with available data provided by open works, the brilliant correctness and broad utilization of the exploited procedure are revealed. Finally, numerical exercises are implemented to examine impacts of gradient functions, constraint conditions, gradation parameters and aspect ratios on free vibration responses of composite smart plates.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"176 ","pages":"Article 106240"},"PeriodicalIF":4.2,"publicationDate":"2025-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143747892","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":"Efficient BEM for thin-walled inhomogeneous potential problems: Theory and MATLAB code","authors":"Xiaotong Gao ,&nbsp;Yan Gu ,&nbsp;Bo Yu ,&nbsp;Wenzhen Qu ,&nbsp;Haodong Ma","doi":"10.1016/j.enganabound.2025.106241","DOIUrl":"10.1016/j.enganabound.2025.106241","url":null,"abstract":"<div><div>The traditional boundary element method (BEM) often faces challenges in efficiently solving inhomogeneous problems, particularly in thin-walled geometries, due to the need for domain discretization and the handling of nearly singular integrals. In this study, we propose an efficient hybrid algorithm that combines the BEM with physics-informed neural networks (PINNs) to solve inhomogeneous potential problems in thin-walled structures. The approach transforms inhomogeneous equations into equivalent homogeneous ones by subtracting a closed-form particular solution, which is derived using the learning capabilities of PINNs. This methodology not only simplifies the problem formulation but also enhances computational efficiency by eliminating the need for domain discretization, making it particularly well-suited for thin-walled geometries. Additionally, the scaled coordinate transformation BEM, a recently developed technique for solving domain integrals, is also employed for comparative analysis. Finally, a nonlinear coordinate transformation is employed to effectively regularize nearly singular integrals, which are critical in BEM for thin structures. The proposed method achieves accurate and reliable results with a small number of boundary elements, even for structures with extremely small thickness-to-length ratios, as low as 10⁻⁹. This makes the method highly suitable for modeling thin films and thin-walled structures, particularly in the context of advanced smart materials. The unique contribution of this work lies in the integration of PINNs with BEM to tackle challenges specific to thin-walled inhomogeneous problems, offering a more efficient and accurate solution compared to traditional BEM-based method.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"176 ","pages":"Article 106241"},"PeriodicalIF":4.2,"publicationDate":"2025-03-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143738462","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 alternative dual reciprocity BEM for P-SV wave propagation problems: A comparative study","authors":"Pouya Kavandi ,&nbsp;Mehdi Panji ,&nbsp;Navid Ganjian ,&nbsp;Jafar Asgari Marnani","doi":"10.1016/j.enganabound.2025.106238","DOIUrl":"10.1016/j.enganabound.2025.106238","url":null,"abstract":"<div><div>This research introduces a dual reciprocity boundary element method (BEM) designed to analyze the transient scattering of vertically travelling incident <em>P</em>-<em>SV</em> waves. By using static fundamental solutions and appropriate predictor operations, the domain inertia integrals from the equilibrium equation were transformed into boundary integral equations. The computable format of the integral equations was achieved by incorporating the effects of free-field displacements into the equations. After coding the formulation, the proposed technique's accuracy and efficiency were assessed through various wave scattering problems, such as an encased round hole, a half-circle depression, and a double-peaked mound, all exposed to vertically incoming <em>P</em>-<em>SV</em> waves. An extensive comparison was conducted with the full-plane time-domain boundary element technique available in the literature, focusing on accuracy and computation time. The findings indicated that despite the complexity and necessity for interior points in dual reciprocity models, it is more advantageous than the traditional full-plane time-domain approach, as it substantially reduces analysis time.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"176 ","pages":"Article 106238"},"PeriodicalIF":4.2,"publicationDate":"2025-03-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143735323","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}