{"title":"The RBF-FD method for solving the time-fractional viscoelastic wave propagation in irregular domains","authors":"Feng Wang, Yan Zhu, Sihao Jia, Xu Guo","doi":"10.1016/j.enganabound.2025.106121","DOIUrl":"10.1016/j.enganabound.2025.106121","url":null,"abstract":"<div><div>The time-fractional viscoelastic wave equation plays a crucial role in geophysical exploration by accurately modeling wave attenuation and velocity dispersion in Earth’s media. However, solving this equation is challenging due to the stress–strain relationship governed by the Caputo fractional derivative of small orders and the complexity of irregular surface topographies. The requirement for significant memory and computational resources when dealing with small fractional orders limits the efficiency of traditional methods. Conventional approaches, which rely on horizontal reference planes, fixed-step grids, and stair-step approximations for irregular surfaces, often lead to staircase scattering and reduced accuracy. To address these challenges, this study proposes a numerical algorithm based on the Radial Basis Function-Finite Difference (RBF-FD) method for simulating time-fractional viscoelastic waves in irregular domains. The meshless nature of the RBF-FD method allows for flexible node distribution, making it well-suited for complex interfaces. Additionally, a short-memory algorithm is implemented to efficiently solve the stress–strain relationship governed by the fractional derivative. Several numerical experiments are presented to validate the accuracy and efficiency of the proposed scheme.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"173 ","pages":"Article 106121"},"PeriodicalIF":4.2,"publicationDate":"2025-02-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143372126","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":"Meshfree analysis of brain tumor growth under various treatment plans considering the mechanical stress effect","authors":"Amir Khosravifard, Ghazaleh Ansari","doi":"10.1016/j.enganabound.2025.106151","DOIUrl":"10.1016/j.enganabound.2025.106151","url":null,"abstract":"<div><div>Brain tumors are among the deadliest types of cancer and the most challenging to treat. Predicting the growth behavior of tumors can help physicians choose the best treatment program. Herein, a numerical technique based on the meshless radial point interpolation methods is presented for simulating the growth of brain tumors under the effects of radiotherapy and chemotherapy. In this work, the stress field induced in the brain due to tumor growth and its effect on tumor diffusion is analyzed. In the meshfree method, to obtain the system of discrete equations, the weak-form of the coupled nonlinear system of governing equations is developed. Part of the simulations is based on clinical data obtained from other research; additionally, data found in the literature is used to validate the results. The effects of various conventional treatment programs, with and without considering the effect of the stress field in brain tissue, are analyzed and compared. Furthermore, the order of radiotherapy and chemotherapy treatments is investigated. It is shown that when the stress effect is considered, the points with maximum stress are where tumor growth is highest. After the cell density reaches its maximum value at these points, growth is transferred to the surrounding areas.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"173 ","pages":"Article 106151"},"PeriodicalIF":4.2,"publicationDate":"2025-02-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143348206","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}
Songlin Hu , Jianwen Liang , Zhenning Ba , Zhongxian Liu
{"title":"Scattering of spherical P-waves by three-dimensional cavity in an elastic half-space","authors":"Songlin Hu , Jianwen Liang , Zhenning Ba , Zhongxian Liu","doi":"10.1016/j.enganabound.2025.106137","DOIUrl":"10.1016/j.enganabound.2025.106137","url":null,"abstract":"<div><div>This study adopts the indirect boundary integral equation method (IBIEM) to solve the scattering of spherical P-waves by a three-dimensional (3D) cavity in an elastic half-space. Specifically, the free field of the spherical wave is obtained by the method of full space superposition. Based on the single-layer potential theory, the scattered field is constructed using concentrated force sources applied on the fictitious wave source surface. Our method’s numerical accuracy and stability are verified by comparing it against existing results. Additionally, considering a spherical cavity in a half-space as an example, this study investigates the influence of the wave source orientation, incident wave frequency, distance between the wave source and the cavity, and cavity depth on the surface displacement and dynamic stress concentration factor (DSCF) on the cavity surface. The results indicate significant differences in the spatial distribution characteristics of surface displacement and DSCF on the cavity surface for different wave source orientations and cavity depths. As the incident frequency increases, the spatial oscillation of surface displacement near the cavity intensifies, and the DSCF gradually decreases. As the wave source approaches the cavity, the amplification effect of surface displacement near the cavity becomes more apparent. At the same time, the maximum DSCF shows significant non-monotonic variation, and its position also changes accordingly. When the distance between the wave source and the cavity is large, the spherical wave can be approximated as a plane wave.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"173 ","pages":"Article 106137"},"PeriodicalIF":4.2,"publicationDate":"2025-02-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143308580","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}
Redouane El Kadmiri , Youssef Belaasilia , Abdelaziz Timesli
{"title":"A hybrid FEM-Meshless approach for dynamic analysis of homogeneous and inhomogeneous structures","authors":"Redouane El Kadmiri , Youssef Belaasilia , Abdelaziz Timesli","doi":"10.1016/j.enganabound.2025.106135","DOIUrl":"10.1016/j.enganabound.2025.106135","url":null,"abstract":"<div><div>In this article, a hybrid FEM-MESHLESS method for dynamic analysis of homogeneous and inhomogeneous structures is developed. This hybrid method has already been shown by El Kadmiri et al (2021, 2022, 2024) for static analysis of homogeneous and inhomogeneous structures. To verify the proposed method for studying the dynamic response of homogeneous and inhomogeneous structures, a study based on small deformations theory is presented by comparison with the results of finite element method and those of analytical solutions. In addition, numerical tests are carried out to demonstrate the reliability and performance of the proposed algorithm by establishing a convergence study using regular and irregular discretization schemes.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"173 ","pages":"Article 106135"},"PeriodicalIF":4.2,"publicationDate":"2025-02-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143308581","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 method based on the method of fundamental solution for time harmonic electromagnetic field with a three-dimensional elastic body","authors":"Yao Sun, Jiaxin Chen","doi":"10.1016/j.enganabound.2024.106056","DOIUrl":"10.1016/j.enganabound.2024.106056","url":null,"abstract":"<div><div>In this paper, we propose a numerical formula to calculate time-harmonic electromagnetic field interacting with three-dimensional elastic body. The formula is based on the method of fundamental solutions. Firstly, we perform Helmholtz decomposition on the displacement field. The problem will transform into a coupled bounded problem including a scaler Helmholtz equation, a vector Helmholtz equation and a Maxwell equation. Then, we use the method of fundamental solutions to solve the new problem. Finally, we provide some examples to demonstrate the effectiveness of the proposed method. We construct the exact solutions for the boundary value problem to verify the accuracy and present a comparative study with the Galerkin scheme.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"171 ","pages":"Article 106056"},"PeriodicalIF":4.2,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142793484","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}
Peng Yin , Xu-Chang Liu , Jin-Shui Yang , Yao-Yao Xu , Shuang Li , Xiao-Bin Lu , Lin-Zhi Wu
{"title":"Numerical study of flow-induced noise response of lightweight composite sandwich sail based on the boundary element method","authors":"Peng Yin , Xu-Chang Liu , Jin-Shui Yang , Yao-Yao Xu , Shuang Li , Xiao-Bin Lu , Lin-Zhi Wu","doi":"10.1016/j.enganabound.2024.106059","DOIUrl":"10.1016/j.enganabound.2024.106059","url":null,"abstract":"<div><div>The submarine's sail, as the largest appendage structure, is more susceptible to turbulence induced vibrations during medium to high-speed navigation, making it a critical area for the generation of flow-induced noise, significantly impacting the stealth and safety of submarine. Considering the excellent mechanical properties and high damping characteristics of lightweight composite sandwich structures, by combining large eddy simulation with vibro-acoustic coupling methods based on boundary element method, under the premise of verifying the accuracy of the numerical methods, a series of three-dimensional dynamic numerical models are established to investigate the flow-induced noise response of the novel composite sandwich sail. The results indicate that the overall sound power level of composite sandwich sail is reduced by approximately 8.9 dB compared to steel structure. The maximum sound power level of composite sandwich sail is lower than the steel with equal areal density. The sound pressure of the sail with buoyant material is lower than that of foam and PVC with the same damping. This work can provide theoretical support for the design methods of new lightweight, multifunctional sail structures.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"171 ","pages":"Article 106059"},"PeriodicalIF":4.2,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142793487","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 IGN-BEM coupled model for the interaction between fully nonlinear waves and 2D floating bodies over variable topography","authors":"Gao-fei Su , Ying Gou , Bin Teng , Ming Zhao","doi":"10.1016/j.enganabound.2024.106077","DOIUrl":"10.1016/j.enganabound.2024.106077","url":null,"abstract":"<div><div>A two-dimensional time domain coupled model is developed to analyze the interaction between fully nonlinear waves and floating bodies over variable topography. The whole calculation domain is divided into an inner domain close to the structure and two outer domains far away from the structure. The fully nonlinear free surface boundary conditions are used in each sub-domain. Irrotational Green-Naghdi (IGN) equations are applied to compute the wave motion in the outer domains, which are solved by the finite element method (FEM). The Laplace equation is solved by the boundary element method (BEM) in the inner domain. The outer and inner domains are coupled through the overlapping regions. The experimental data of waves propagating over a submerged breakwater and the interaction between shallow-water waves and a box fixed on the still-water surface are used to verify the rationality and accuracy of the coupled model. The coupled model is applied to compute the wave exciting force and the motion response of a barge over a sloping terrain. The influence of the terrain height on wave forces and barge motions is studied.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"171 ","pages":"Article 106077"},"PeriodicalIF":4.2,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142825004","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 total Lagrangian‒Riemann SPH method with MUSCL reconstruction for large elastic‒plastic deformation and fracture simulation","authors":"Longkui Chen , Zhanming Wang , Shenghong Huang","doi":"10.1016/j.enganabound.2024.106083","DOIUrl":"10.1016/j.enganabound.2024.106083","url":null,"abstract":"<div><div>The smooth particle hydrodynamics (SPH) method possesses inherent advantages in simulating large deformations, fractures and crack propagations in solids. However, challenging issues, including tensile instability and numerical oscillations, persist. Total Lagrangian smooth particle hydrodynamics (TLSPH) was proposed to eliminate tensile instability by applying the kernel approximation consistently in the reference configuration; however, the artificial viscosity model has to be added to reduce the numerical oscillation induced by shock and other contact discontinuity simulations, which severely decreases its accuracy and robustness. Motivated by the advantages of both TLSPH and Riemann-SPH of the ULSPH frame, a second-order solid Riemann scheme is constructed on the basis of the Monotone Upwind-Centered Scheme for Conservation Laws (MUSCL) reconstruction and incorporated into the total Lagrangian SPH (TLSPH) framework. The resulting MUSCL-TLSPH method is designed for solving dynamic elastic‒plastic structural impact problems, including large deformations and fractures. This method effectively overcomes the challenges faced by traditional SPH approaches, eliminating the need to introduce artificial stresses related to tunable parameters to maintain computational stability. Finally, the accuracy and robustness of the MUSCL-TLSPH method presented in this paper are verified through a series of numerical validations.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"171 ","pages":"Article 106083"},"PeriodicalIF":4.2,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142901668","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}
Limei Zhang , Yueping Yin , Hong Zheng , Sainan Zhu , Nan Zhang
{"title":"Singularity treatments in transient confined seepage using numerical manifold method","authors":"Limei Zhang , Yueping Yin , Hong Zheng , Sainan Zhu , Nan Zhang","doi":"10.1016/j.enganabound.2024.106100","DOIUrl":"10.1016/j.enganabound.2024.106100","url":null,"abstract":"<div><div>The numerical manifold method (NMM) is proposed for analysis of the two-dimensional transient confined seepage flow problems with singular corner points. To deal with the singularity of corner points, the asymptotic expansion of the solution in the vicinity of corner points is incorporated into the local approximations of the relevant physical patches of the NMM, while the constant local approximation is assigned to the other patches far from the singularity points. Then, the NMM discrete formulation for the initial – boundary value problem for transient seepage flow is deduced based on the Galerkin approximation. For time integration, the backward time integration scheme is adopted. The accuracy and effectiveness of the proposed method are demonstrated in typical examples involving homogeneous, heterogeneous, and anisotropic material. Comparing with constant local approximations to all the patches, the proposed method can better reflect the strong singularity of corner points.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"171 ","pages":"Article 106100"},"PeriodicalIF":4.2,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142929263","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 multiscale and multiphysical numerical approach for sandwich multiphase hybrid fiber plates with smart composite facesheets","authors":"Duy-Khuong Ly , Huy-Cuong Vu-Do , Chanachai Thongchom , T. Nguyen-Thoi","doi":"10.1016/j.enganabound.2025.106134","DOIUrl":"10.1016/j.enganabound.2025.106134","url":null,"abstract":"<div><div>This study introduces a comprehensive multiscale and multiphysical numerical approach for analyzing sandwich three-phase nanocomposite plate with multiferroic facesheets in its upper and lower surfaces. The proposed research investigates the zigzag effect and quasi-3D sinusoidal shear deformation, capturing the complex interactions between the core and multiferroic facesheets across multiple physical fields. A distinct feature of the three-phase polymer/CNT/fiber material is the embedding of Carbon Nanotube (CNT) nanofiller within the matrix phase, enhancing the overall properties of the carbon fiber composite. Micromechanical models for three-phase systems are employed to determine the effective elastic properties of the composite core. A unified numerical approach is developed to address the global and local behavior of the structure, capturing the mechanical, electrical, and magnetic coupling effects inherent in multiferroic materials. This model utilizes isogeometric analysis for high-fidelity representation, ensuring precise geometric accuracy and smooth continuity, and incorporates Eringen’s nonlocal strain gradient multiferroic theory to account for size effects. The zigzag effect is characterized by a multiscale kinematic description, where the displacement field is represented by the superposition of coarse and fine contributions. Numerical simulations validate the model, demonstrating its effectiveness in predicting the mechanical, electrical, and magnetic responses of the smart composite plates. This work offers a robust tool for the design and optimization of advanced composite structures in engineering applications.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"173 ","pages":"Article 106134"},"PeriodicalIF":4.2,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143077701","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}