Engineering Analysis with Boundary Elements最新文献

{"title":"Isogeometric methods for thermal analysis with spatially varying thermal conductivity under general boundary and other constraints","authors":"Zulfiqar Ali ,&nbsp;Weiyin Ma","doi":"10.1016/j.enganabound.2025.106130","DOIUrl":"10.1016/j.enganabound.2025.106130","url":null,"abstract":"<div><div>This paper presents some results on steady-state thermal analysis with variable thermal conductivity under general boundary conditions and other internal constraints using isogeometric methods. Non-Uniform Rational B-splines (NURBS) serve as basis functions for representing both the geometry of the physical domains and the solution. While both isogeometric collocation method and Galarkin formulation are discussed for facilitating comparisons, the main emphasis of the presented work is on isogeometric collocation method (IGA-C) for thermal analysis. To obtain the final solution, the respective partial differential equation (PDE) is discretized in its strong form at a number of collocation sites in IGA-C, as opposed to Galerkin formulations that involve a costly process of numerical integration in building up the system equations. The proposed method on IGA-C for thermal analysis can be easily implemented due to the simplicity of IGA-C in setting up the system equations. In addition to general boundary conditions of the respective PDE, other arbitrary constraints can also be easily incorporated into the final system of equations for producing desired solutions. Numerical examples with different kinds of spatially varying thermal conductivity along with other additional constraints and heat sources are provided to demonstrate the effectiveness of the proposed methods. The results show that the proposed methods are capable of conveniently handling arbitrary boundary and other additional constraints when solving thermal PDEs and can produce stable and accurate solutions with expected convergence.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"172 ","pages":"Article 106130"},"PeriodicalIF":4.2,"publicationDate":"2025-01-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143077673","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":"Electromagnetic scattering sensitivity analysis for perfectly conducting objects in TM polarization with isogeometric BEM","authors":"Leilei Chen ,&nbsp;Chengmiao Liu ,&nbsp;Haojie Lian ,&nbsp;Wenxiang Gu","doi":"10.1016/j.enganabound.2025.106126","DOIUrl":"10.1016/j.enganabound.2025.106126","url":null,"abstract":"<div><div>This study proposes a sensitivity analysis framework for Transverse Magnetic polarized electromagnetic scattering problems, with a focus on Perfectly Electric Conductors (PEC). To enable seamless integration of Computer-Aided Design and Computer-Aided Engineering, the isogeometric boundary element method based on the Galerkin scheme is employed. This method utilizes Non-Uniform Rational B-splines to represent geometries and unknown physical fields. A direct differentiation formulation is derived to evaluate the sensitivity of the Radar Cross Section (RCS) with respect to shape design variables and incident wave angles. The accuracy and effectiveness of the algorithm are validated through numerical examples.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"172 ","pages":"Article 106126"},"PeriodicalIF":4.2,"publicationDate":"2025-01-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143049713","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 cell-based smoothed radial point interpolation method applied to lower bound limit analysis of thin plates","authors":"Shenshen Chen ,&nbsp;Hao Dong ,&nbsp;Xing Wei ,&nbsp;Fengtao Liu","doi":"10.1016/j.enganabound.2025.106123","DOIUrl":"10.1016/j.enganabound.2025.106123","url":null,"abstract":"<div><div>This paper proposes a novel numerical method based on the cell-based smoothed radial point interpolation method (CS-RPIM) combined with second-order cone programming to perform lower bound limit analysis of elastic-perfectly-plastic thin plates, using only deflection as nodal variable. The problem domain is initially discretized using a simple triangular background mesh, where each triangular cell is subsequently subdivided into multiple smoothing domains. Shape functions are formulated using the radial point interpolation method, allowing direct imposition of essential boundary conditions for deflection. Rotational constraints are conveniently handled through the construction of smoothed curvatures. By utilizing a generalized gradient smoothing technique, complex domain integrals are simplified into boundary integrals over the smoothing domains, thus eliminating the need to compute second-order derivatives of the shape functions. The virtual work principle is employed to enforce the equilibrium conditions for the self-equilibrated residual moment field in a weak sense. The von Mises yield conditions are expressed as conic constraints and the resulting optimization problems are solved using highly efficient primal-dual interior point solvers. Numerical examples demonstrate that it is feasible and effective to conduct lower bound limit analysis of thin plates using the proposed CS-RPIM and second-order cone programming.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"172 ","pages":"Article 106123"},"PeriodicalIF":4.2,"publicationDate":"2025-01-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143049714","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":"Bi-material V-notch fracture analysis in functionally graded materials","authors":"C.Y. Fu ,&nbsp;Y. Yang ,&nbsp;P.H. Wen ,&nbsp;J. Sladek ,&nbsp;V. Sladek","doi":"10.1016/j.enganabound.2025.106127","DOIUrl":"10.1016/j.enganabound.2025.106127","url":null,"abstract":"<div><div>The finite block method (FBM) in the Cartesian coordinate system is developed to deal with the problems of the bi-materials V-notches in functionally graded materials (FGM) under static and dynamic loads. The first partial differential matrix is established via Lagrange series. Higher-order derivatives can be deduced from the first order partial differential matrix directly. In order to obtain the high accurate at the V-notch tip, the asymptotic expansions of the stress and displacement around the notch tip are introduced with a singular polygonal core technique. For the dynamic problems, the Laplace transform method with Durbin inverse algorithm is utilized. The degrees of accuracy and convergence of the FBM are demonstrated through four case studies. Comparisons are implemented with the finite element method (FEM) results.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"172 ","pages":"Article 106127"},"PeriodicalIF":4.2,"publicationDate":"2025-01-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143049715","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":"Green’s function representation and numerical approximation of the two-dimensional stochastic Stokes equation","authors":"Jie Zhu ,&nbsp;Yujun Zhu ,&nbsp;Ju Ming ,&nbsp;Xiaoming He","doi":"10.1016/j.enganabound.2025.106117","DOIUrl":"10.1016/j.enganabound.2025.106117","url":null,"abstract":"<div><div>This paper investigates the two-dimensional unsteady Stokes equation with general additive noise. The primary contribution is the derivation of the relevant estimate of Green’s tensor, which provides a fundamental representation for the solution of this stochastic equation. We demonstrate the crucial role of Green’s function in understanding the stability and perturbation characteristics of the stochastic Stokes system. Furthermore, we analyze the convergence properties of the Euler–Maruyama (EM) scheme for temporal discretization and derive error estimates for a Galerkin finite element discretization using the Taylor–Hood method for spatial approximation. This work provides a strong convergence of order <span><math><mrow><mi>O</mi><mfenced><mrow><mi>h</mi><msup><mrow><mrow><mo>(</mo><mi>Δ</mi><mi>t</mi><mo>)</mo></mrow></mrow><mrow><mo>−</mo><mfrac><mrow><mn>1</mn></mrow><mrow><mn>2</mn></mrow></mfrac></mrow></msup><mo>+</mo><mi>Δ</mi><mi>t</mi></mrow></mfenced></mrow></math></span> of the velocity in the <span><math><mrow><msup><mrow><mi>L</mi></mrow><mrow><mn>2</mn></mrow></msup><mrow><mo>(</mo><mn>0</mn><mo>,</mo><mi>T</mi><mo>;</mo><msup><mrow><mi>L</mi></mrow><mrow><mn>2</mn></mrow></msup><mrow><mo>(</mo><mi>D</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow></math></span> norm and <span><math><mrow><mi>O</mi><mfenced><mrow><mi>h</mi><msup><mrow><mrow><mo>(</mo><mi>Δ</mi><mi>t</mi><mo>)</mo></mrow></mrow><mrow><mo>−</mo><mfrac><mrow><mn>3</mn></mrow><mrow><mn>2</mn></mrow></mfrac></mrow></msup><mo>+</mo><mi>Δ</mi><mi>t</mi></mrow></mfenced></mrow></math></span> of the pressure in the <span><math><mrow><msup><mrow><mi>L</mi></mrow><mrow><mn>2</mn></mrow></msup><mrow><mo>(</mo><mn>0</mn><mo>,</mo><mi>T</mi><mo>;</mo><msup><mrow><mi>H</mi></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mi>D</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow></math></span> norm based on the Green tensor approach. These results contribute to a deeper understanding of the stochastic behavior of fluid dynamics systems, paving the way for improved theoretical modeling and more accurate numerical simulations in diverse fields such as meteorology, oceanography, and engineering applications.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"172 ","pages":"Article 106117"},"PeriodicalIF":4.2,"publicationDate":"2025-01-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143049716","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 fully cell-based immersed smoothed finite element method with the mean value coordinate projection using quadrilateral elements for fluid-structure interaction","authors":"Shuhao Huo ,&nbsp;Hengzhi Wang ,&nbsp;Zhipeng Li ,&nbsp;Zhiqiang Li ,&nbsp;Chen Jiang ,&nbsp;Guirong Liu","doi":"10.1016/j.enganabound.2025.106122","DOIUrl":"10.1016/j.enganabound.2025.106122","url":null,"abstract":"<div><div>In this work, an effective and stable immersed cell-based smoothed finite element method (ICS-FEM) together with mean value coordinate (MVC) projection using quadrilateral elements is presented for 2D fluid-structure interaction (FSI) problems. In an immersed-based algorithm, the entire system can be divided into three components: large-deformed nonlinear structure, incompressible viscous fluid, and fictitious fluid for FSI force. A characteristic-based split (CBS) CS-FEM solver is developed to solve the Navier-Stokes (N-S) equation. An explicit total Lagrange (TL) CS-FEM solver is utilized to depict the large deformation of the elastic structure. The problem domains are discretized into a set of quadrilateral elements and four cell-based smoothing domains are constructed for each element. The cell-based smoothing operation instead of the traditional mapping operation of isoparametric elements is executed to all gradient-related terms for both structures and fluids. A method based on the mean value coordinates is used for the interpolation process and the search criterion. Four numerical examples are utilized to illustrate the advance of the proposed method, including computational precision, convergence rate, robustness, etc.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"172 ","pages":"Article 106122"},"PeriodicalIF":4.2,"publicationDate":"2025-01-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143050030","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 new boundary integral method for investigating the roughness scaling law of heterogeneous interfacial fracture","authors":"Wei Du ,&nbsp;Xiaohua Zhao ,&nbsp;Wei Jiang ,&nbsp;Yongcheng Guo ,&nbsp;Jinping Fu ,&nbsp;Zhen Wang","doi":"10.1016/j.enganabound.2025.106131","DOIUrl":"10.1016/j.enganabound.2025.106131","url":null,"abstract":"<div><div>A new semi-analytical and semi-numerical approach is proposed to investigate the scaling law of in-plane roughness due to the fracture of a heterogeneous interface involving spatial correlation of disorders. The model is considered as a composite structure composed of two cantilever rectangular plates bonded with an interfacial layer. Based on the theory of solid mechanics, the dynamic process of interfacial fracture is derived analytically and reduced to two coupled integral equations, which further become a system of linear algebraic equations after discretizing the interface to a set of prismatic elements. Numerical simulations present that the morphology of interfacial fracture fronts in all cases show self-affine scaling properties with the roughness exponent in the range (0.36,0.70), depending on stiffness ratio of laminate structure and heterogeneous properties of interface. Remarkably, the present results cover most of the exponent values observed in previous experiments and provides strong evidence that it is the microstructure and heterogeneous properties that mainly control the roughness.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"172 ","pages":"Article 106131"},"PeriodicalIF":4.2,"publicationDate":"2025-01-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143049717","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":"Seismic response analysis of marine undulating sites based on indirect boundary element method","authors":"Zhong-Xian Liu,&nbsp;Xiang Liu,&nbsp;Tian-Chun Ai,&nbsp;Jia-Wei Zhao,&nbsp;Lei Huang","doi":"10.1016/j.enganabound.2025.106125","DOIUrl":"10.1016/j.enganabound.2025.106125","url":null,"abstract":"<div><div>The seafloor is mostly made up of soft silt, and seismic waves collide with particles before scattering during their propagation. Moreover, the ocean terrain includes basins, seamounts, islands and reefs, contributing to the intricate propagation of seismic waves in seawater. This study proposes a two-dimensional wave simulation algorithm for the marine seismic site effect based on the indirect boundary element method (IBEM). The \"seafloor bedrock - complex site - seawater\" system is used to simulate the scattering characteristics of seismic waves. The method proposed in this study can significantly improve the computing performance while drastically reducing the computation and storage requirements. Porosity has a significant impact on the displacement amplitude factor (DAF) and peak ground acceleration (PGA), with increases in porosity leading to over 30 % growth in both metrics. When the terrain becomes more complex, the resulting seismic amplification effects are significantly enhanced, causing DAF and PGA to increase by 72.0 % and 32.9 %, respectively. A comprehensive analysis reveals that porosity not only affects the propagation characteristics of seismic waves but also exacerbates site amplification effects through complex topographic conditions. Therefore, the influence of topographic factors on the site effect can be considered during the design and construction of relevant ocean engineering to guarantee the safety and stability of the structure.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"172 ","pages":"Article 106125"},"PeriodicalIF":4.2,"publicationDate":"2025-01-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143049718","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":"Modeling of rarefied gas flows in streamwise periodic channels: Application of coupled constitutive relations and the method of fundamental solutions","authors":"Ankit Farkya,&nbsp;Anirudh Singh Rana","doi":"10.1016/j.enganabound.2024.106108","DOIUrl":"10.1016/j.enganabound.2024.106108","url":null,"abstract":"<div><div>Periodic structures are ubiquitous in nature and engineering, offering unique properties that inspire a range of applications. This paper explores the mathematical modeling of periodic structures in rarefied gas flows using the coupled constitutive relations (CCR) model. The method of fundamental solutions (MFS), known for its meshfree nature and computational efficiency, is utilized as a numerical tool. The streamwise periodic boundary conditions are incorporated into the MFS for modeling two-dimensional flow in periodically patterned channels. We validate the developed CCR-MFS framework with analytical solutions for force-driven Poiseuille and Couette flow. The error analysis is also performed to determine the optimal singularity location. Furthermore, we simulate the flow in channels with periodic patterns by varying the accommodation coefficient. This allows us to analyze the effects of patterning and accommodation coefficients in the Maxwell model of boundary conditions. Effects of patterning on mass flux, energy flux, and average friction coefficients are also presented for the force-driven flow in patterned channels. Our simulations demonstrate the potential of the mathematical and computational techniques to enhance the performance and functionality of a range of technological applications.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"172 ","pages":"Article 106108"},"PeriodicalIF":4.2,"publicationDate":"2025-01-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142989686","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":"Extension of three-dimensional discontinuous deformation analysis for solid block motions in predefined fluid field","authors":"Xinyan Peng ,&nbsp;Xuanmei Fan ,&nbsp;Pengcheng Yu ,&nbsp;Guangqi Chen ,&nbsp;Mingyao Xia ,&nbsp;Yingbin Zhang ,&nbsp;Xiao Cheng ,&nbsp;Chao Liang","doi":"10.1016/j.enganabound.2025.106124","DOIUrl":"10.1016/j.enganabound.2025.106124","url":null,"abstract":"<div><div>Solid–fluid numerical simulations involving open channels are usually complicated, especially for large solid displacements. An extended three–dimensional discontinuous deformation analysis (3D DDA) method incorporating depth-integrated two-dimensional fluid dynamics was proposed to evaluate solid movement considering fluid actions. In this method, two types of fluid forces on solid blocks, buoyancy and drag forces, are calculated under complex fluid surface conditions and behave as external forces on solid blocks. The accuracy of the extended 3D DDA method in calculating the block buoyancy and fluid drag forces was validated using a block acceleration model. A hypothetical open-channel model involving non-uniform and asymmetric terrain was studied to show the extended 3D DDA method is applicable to study sliding mass deposition in general cases. A hypothetical bridge unseating case threatened by flood and floating objects was then studied, which showed that both the uprising flood and floating objects could cause serious damage to a bridge in a flood event. Finally, the deposition of sliding mass in the Yangjiagou landslide was simulated, demonstrating the influence of fluid dynamics on the deposition characteristics. The results showed that the extended 3D DDA method is adaptable and efficient in analyzing discrete block systems involving open-channel fluid action.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"172 ","pages":"Article 106124"},"PeriodicalIF":4.2,"publicationDate":"2025-01-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142989684","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}