{"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":"https://doi.org/10.1016/j.enganabound.2024.106083","url":null,"abstract":"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.","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"5 1","pages":""},"PeriodicalIF":3.3,"publicationDate":"2024-12-24","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}
{"title":"Various near singularity regularization methods derived from distance transformations in 3D boundary element method","authors":"Yukai Jin, Yidan Zhang","doi":"10.1016/j.enganabound.2024.106094","DOIUrl":"https://doi.org/10.1016/j.enganabound.2024.106094","url":null,"abstract":"This paper focuses on applying non-linear transformations for near singularity regularization combined distance transformations. In the previous methods, near singularities are usually considered in only the polar direction, ignoring those in the circular direction, which leads to low accuracy when calculating nearly singular integrals of narrow element or when the projection point is located near the element end. In this paper, the near singularities are traced firstly based on the distance function, by which the distance function can be constructed in two ways. The general form of nearly singular integrals in the two directions is extracted. Then, several non-linear transformations are introduced about removal of the near singularities in one direction. In our method, the only one directional methods are combined to solve the nearly singular integrals. Finally, comparisons of the results by the combined distance transformations show that by employing the non-linear transformations in both directions, more stable and accurate results can be obtained especially for nearly singular integrals of the narrow elements.","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"6 1","pages":""},"PeriodicalIF":3.3,"publicationDate":"2024-12-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142901667","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}
Ali Habibirad, Yadollah Ordokhani, Omid Baghani, Hadis Azin
{"title":"A local meshless numerical scheme based on the radial point interpolation for the generalized time-fractional Allen–Cahn equation","authors":"Ali Habibirad, Yadollah Ordokhani, Omid Baghani, Hadis Azin","doi":"10.1016/j.enganabound.2024.106058","DOIUrl":"https://doi.org/10.1016/j.enganabound.2024.106058","url":null,"abstract":"This research has been conducted to investigate a numerical solution for the Allen–Cahn equation featuring the generalized fractional time derivative. The finite difference method is employed to discretize the equation in the time variable. Subsequently, an error estimate is derived for the proposed method in <mml:math altimg=\"si1.svg\" display=\"inline\"><mml:msub><mml:mrow><mml:mi mathvariant=\"script\">L</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi><mml:mo>,</mml:mo><mml:mi>μ</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant=\"bold\">q</mml:mi></mml:mrow></mml:msub></mml:math> space. Furthermore, a meshless technique based on radial point interpolation is used to discretize the problem in spatial variables. Through these procedures, the equation is transformed into a system of linear equations at each time step. The method’s effectiveness for solving this equation is demonstrated by three examples on both regular and irregular domains. These examples illustrate that the current method has a high level of accuracy and efficiency for solving the given problem.","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"32 1","pages":""},"PeriodicalIF":3.3,"publicationDate":"2024-12-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142874805","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":"Numerical analysis of basin response using Indirect Boundary Element Method (IBEM) for dip-slip sources","authors":"Zhonghan Liu, Zhenning Ba, Jingxuan Zhao, Jiaqi Niu","doi":"10.1016/j.enganabound.2024.106081","DOIUrl":"https://doi.org/10.1016/j.enganabound.2024.106081","url":null,"abstract":"Accurate modeling of basin structures and quantitative analysis of basin amplification effects are critical for seismologists and engineers. The Indirect Boundary Element Method (IBEM), developed from the Boundary Element Method (BEM), is particularly well-suited for these tasks due to its capability to manage layers with lateral inhomogeneities. However, current IBEM studies mostly focus on wavefields generated by plane wave incidence or point sources emitting pure spherical P or S waves, which often fail to capture the complex behavior of real seismic events. To address this limitation, we propose an updated IBEM scheme that simulates wavefields from dip-slip sources using the displacement-stress discontinuity vector. We validate this updated IBEM through degenerated models and sources, confirming its robustness in handling dip-slip sources. Further numerical tests reveal several features of basin amplification effects from dip-slip sources: (1) Sources with varying dips cause substantial differences in-basin peak numbers and amplitudes compared to simplified sources. For the surface responses, the ratio of maximum amplitudes (in-basin to out-of-basin) can exceed 10 with asymmetric distributed peaks. (2) The basin amplification effect for dip-slip sources shows similar frequency and velocity contrast dependencies compared to other simplified sources. Our study shows the IBEM scheme is effective and advantageous in basin amplification effect analysis, especially for different source configurations under the same model, with great potential for further applications for seismology, earthquake engineering, and hazard control.","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"16 1","pages":""},"PeriodicalIF":3.3,"publicationDate":"2024-12-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142825005","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":"https://doi.org/10.1016/j.enganabound.2024.106077","url":null,"abstract":"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.","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"27 1","pages":""},"PeriodicalIF":3.3,"publicationDate":"2024-12-12","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}
V. Gutiérrez-Posada, I. Ubero-Martínez, J. Cifuentes-Rodrǵuez, J. Vallepuga-Espinosa
{"title":"BEM formulation for simulating heat dissipation in microelectronic packaging with point heat sources","authors":"V. Gutiérrez-Posada, I. Ubero-Martínez, J. Cifuentes-Rodrǵuez, J. Vallepuga-Espinosa","doi":"10.1016/j.enganabound.2024.106076","DOIUrl":"https://doi.org/10.1016/j.enganabound.2024.106076","url":null,"abstract":"This work presents a new and robust formulation for studying the effect of point heat sources on three-dimensional thermomechanical contact problems. The aim of this work is to accurately analyze heat dissipation in microchips with known heat sources. To achieve this, the Boundary Element Method (BEM) has been used to calculate the thermomechanical influence coefficients. The traditional BEM has been extended in the proposed formulation to account for the complex interactions between the microchip and heatsink, which result from the inclusion of point heat sources. This includes the effects of microcontacts, microgaps, and thermal contact resistance. To solve the problem’s inherent non-linearity, a dual iterative solution scheme has been proposed. The formulation is validated by a numerical example and is shown to be capable of accurately capturing traction variations in the contact zone and the influence of heat sources on the contact zone temperature. The formulation is then applied to study the thermomechanical behavior of microelectronic packages, particularly those with low heat dissipation, under the influence of simulated point heat sources.","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"16 1","pages":""},"PeriodicalIF":3.3,"publicationDate":"2024-12-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142825049","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 volume compensation model for multi-resolution moving particle method simulating free-surface flow","authors":"Xiaoxing Liu, Kai Wang, Shunhua Chen, Shuai Zhang","doi":"10.1016/j.enganabound.2024.106080","DOIUrl":"https://doi.org/10.1016/j.enganabound.2024.106080","url":null,"abstract":"This study presents a novel volume compensation model for multi-resolution moving particle method simulating free surface flows. The volume-compensation model is developed to conserve volume when simulating free surface flow using multi-resolution particles, a topic that has been rarely discussed for multi-resolution simulations in previous literature. The free surface is reconstructed by a linear polynomial, enabling the volume of fluid occupied to be easily calculated. Then particles on the free surface are shifted according to the volume change to conserve volume. In the multi-resolution model, the particle size is capable of adapting in an adaptive manner in accordance with its position. In the refined area, the large particles split into four smaller child particles. Child particles near to the large particles are further merged into middle-size particle to avoid the instability caused by the contact of child particles and parent particles. In order to ensure the accuracy of the calculations, the LSMPS framework is employed. A series of numerical tests were conducted to validate the capability of the developed method using multi-resolution particles.","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"116 1","pages":""},"PeriodicalIF":3.3,"publicationDate":"2024-12-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142825129","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 rational kernel function selection for Galerkin meshfree methods through quantifying relative interpolation errors","authors":"Like Deng, Dongdong Wang","doi":"10.1016/j.enganabound.2024.106062","DOIUrl":"https://doi.org/10.1016/j.enganabound.2024.106062","url":null,"abstract":"Although kernel functions play a pivotal role in meshfree approximation, the selection of kernel functions is often experience-based and lacks a theoretical basis. As an attempt to resolve this issue, a rational matching between kernel functions and nodal supports is proposed in this work for Galerkin meshfree methods, where the quadratic through quintic B-spline kernel functions are particularly investigated. The foundation of this rational matching is the design of an efficient quantification of relative interpolation errors. The proposed relative interpolation error measures are not problem-dependent and can be easily and efficiently evaluated. More importantly, these relative interpolation error measures effectively reflect the variation of the real interpolation errors for meshfree approximation, which essentially control the solution accuracy of the Galerkin meshfree formulation with consistent numerical integration. Consequently, an optimal selection of kernel functions that match the nodal supports of meshfree approximation can be readily realized via minimizing the relative interpolation errors of meshfree approximation. The efficacy of the proposed rational matching between kernel functions and nodal supports is well demonstrated by meshfree numerical solutions.","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"210 1","pages":""},"PeriodicalIF":3.3,"publicationDate":"2024-12-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142793486","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":"Numerical simulation of fracture and breakage behaviors in rock disks containing pre-defects with an improved non-local model","authors":"Shijun Zhao, Liang Kong, Qing Zhang, Xinbo Zhao, Wei Xu","doi":"10.1016/j.enganabound.2024.106061","DOIUrl":"https://doi.org/10.1016/j.enganabound.2024.106061","url":null,"abstract":"The characterization and understanding of cracking propagation behaviors in non-uniform geological structures are crucial for predicting the mechanical response of rock-like materials under varying loading conditions. In this study, an improved Peridynamics (PD) model with degree of heterogeneity characterized by random pre-breaking \"bond\" ratio is introduced to capture the intricacies of crack initiation, and propagation in Brazilian disk (BD) specimens with different pre-cracks was explored. The crack propagation paths and the whole process of breakage of specimens with single pre-crack or pre-crack system were obtained. Next, the splitting mechanism and bearing capacity of the BD specimens with different inclinations of pre-crack were investigated through the improved PD model. Finally, the results of the PD simulations based on non-local actions were compared with the experimental ones in terms of crack propagation and breakage patterns. The simulation results are in good agreement with the experimental results, which can help explain the observed behaviors of cracking propagation and growth behaviors. The simulation results show that the improved PD model can accurately describe the whole process of micro-cracks initiation, aggregation, macro-crack generation, propagation and breakage process. This study provides valuable insights for improving simulations and contributing to the field of geomechanics.","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"14 1","pages":""},"PeriodicalIF":3.3,"publicationDate":"2024-12-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142825136","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":"https://doi.org/10.1016/j.enganabound.2024.106056","url":null,"abstract":"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.","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"20 1","pages":""},"PeriodicalIF":3.3,"publicationDate":"2024-12-07","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}