首页 > 最新文献

Engineering Analysis with Boundary Elements最新文献

英文 中文
Point Jacobi-type preconditioning and parameter tuning for Calderon-preconditioned Burton–Miller method in transmission problems 传输问题中Calderon-preconditioned Burton-Miller方法的点jacobi型预处理和参数整定
IF 4.1 2区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2026-02-01 Epub Date: 2025-12-03 DOI: 10.1016/j.enganabound.2025.106572
Keigo Tomoyasu, Hiroshi Isakari
It was recently demonstrated that the boundary element method based on the Burton–Miller formulation (BM-BEM), widely used for solving exterior problems, can be adapted to solve transmission problems efficiently. This adaptation utilises Calderon’s identities to improve the spectral properties of the underlying integral operator. Consequently, most eigenvalues of the squared BEM coefficient matrix, i.e., the collocation-discretised version of the operator, cluster at a few points in the complex plane. When these clustering points are closely packed, the resulting linear system is well-conditioned and can be solved efficiently using the generalised minimal residual (GMRES) method with only a few iterations. However, when multiple materials with significantly different material constants are involved, some eigenvalues become separated, deteriorating the conditioning. To address this, we propose an enhanced Calderon-preconditioned BM-BEM with two strategies. First, we apply a preconditioning scheme inspired by the point Jacobi method. Second, we tune the BM parameters to improve the conditioning of the coefficient matrix. Both strategies leverage a newly derived analytical expression for the eigenvalue clustering points of the relevant operator. Numerical experiments demonstrate that the proposed method, combining both strategies, is particularly efficient for solving scattering problems involving composite penetrable materials with high contrast in material properties.
近年来的研究表明,基于Burton-Miller公式的边界元法(BM-BEM)可以有效地适用于求解外部问题。这种自适应利用卡尔德隆恒等式来改善底层积分算子的谱性质。因此,平方的BEM系数矩阵的大多数特征值,即算子的配位离散版本,聚集在复平面上的几个点上。当这些聚类点紧密聚集时,得到的线性系统是条件良好的,并且可以使用广义最小残差(GMRES)方法,只需少量迭代即可有效地求解。然而,当涉及多个材料且材料常数显著不同时,一些特征值会分离,使条件恶化。为了解决这个问题,我们提出了一个具有两种策略的增强型卡尔德龙预置BM-BEM。首先,我们采用了一种受点Jacobi法启发的预处理方案。其次,我们调整了BM参数,以改善系数矩阵的条件。这两种策略都利用了相关算子的特征值聚类点的新导出的解析表达式。数值实验表明,本文提出的方法结合了这两种策略,特别有效地解决了材料性能具有高对比度的复合可穿透材料的散射问题。
{"title":"Point Jacobi-type preconditioning and parameter tuning for Calderon-preconditioned Burton–Miller method in transmission problems","authors":"Keigo Tomoyasu,&nbsp;Hiroshi Isakari","doi":"10.1016/j.enganabound.2025.106572","DOIUrl":"10.1016/j.enganabound.2025.106572","url":null,"abstract":"<div><div>It was recently demonstrated that the boundary element method based on the Burton–Miller formulation (BM-BEM), widely used for solving exterior problems, can be adapted to solve transmission problems efficiently. This adaptation utilises Calderon’s identities to improve the spectral properties of the underlying integral operator. Consequently, most eigenvalues of the squared BEM coefficient matrix, i.e., the collocation-discretised version of the operator, cluster at a few points in the complex plane. When these clustering points are closely packed, the resulting linear system is well-conditioned and can be solved efficiently using the generalised minimal residual (GMRES) method with only a few iterations. However, when multiple materials with significantly different material constants are involved, some eigenvalues become separated, deteriorating the conditioning. To address this, we propose an enhanced Calderon-preconditioned BM-BEM with two strategies. First, we apply a preconditioning scheme inspired by the point Jacobi method. Second, we tune the BM parameters to improve the conditioning of the coefficient matrix. Both strategies leverage a newly derived analytical expression for the eigenvalue clustering points of the relevant operator. Numerical experiments demonstrate that the proposed method, combining both strategies, is particularly efficient for solving scattering problems involving composite penetrable materials with high contrast in material properties.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"183 ","pages":"Article 106572"},"PeriodicalIF":4.1,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145685097","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}
引用次数: 0
An RBF–FD method for pricing under the Bates model: Handling stochastic volatility and jump processes 贝茨模型下的RBF-FD定价方法:处理随机波动和跳跃过程
IF 4.1 2区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2026-02-01 Epub Date: 2025-12-30 DOI: 10.1016/j.enganabound.2025.106622
Yicheng Liu , Yilin Li , Tao Liu
This work introduces a meshless computational methodology for the valuation of European-style options under the Bates stochastic volatility jump-diffusion framework. The governing dynamics are described by a two-dimensional partial integro-differential equation (PIDE) that incorporates both mixed derivative operators and a nonlocal integral component. The numerical strategy is constructed within the Radial Basis Function-Finite Difference (RBF–FD) paradigm, employing a modified multiquadric kernel that facilitates the analytical determination of localized differentiation weights. Spatial discretization is achieved through the method of lines, where the integral operator is approximated by means of quadrature rules combined with accurate interpolation techniques. Numerical experiments reveal that the proposed solver delivers superior accuracy while lowering computational cost in comparison with classical schemes.
本文介绍了在贝茨随机波动跳-扩散框架下欧式期权估值的无网格计算方法。控制动力学由包含混合导数算子和非局部积分分量的二维偏积分微分方程(PIDE)描述。数值策略是在径向基函数-有限差分(RBF-FD)范式内构建的,采用改进的多重二次核,便于解析确定局部微分权值。空间离散化是通过线的方法实现的,其中积分算子是通过正交规则结合精确的插值技术来逼近的。数值实验表明,与经典格式相比,该算法具有更高的精度和更低的计算成本。
{"title":"An RBF–FD method for pricing under the Bates model: Handling stochastic volatility and jump processes","authors":"Yicheng Liu ,&nbsp;Yilin Li ,&nbsp;Tao Liu","doi":"10.1016/j.enganabound.2025.106622","DOIUrl":"10.1016/j.enganabound.2025.106622","url":null,"abstract":"<div><div>This work introduces a meshless computational methodology for the valuation of European-style options under the Bates stochastic volatility jump-diffusion framework. The governing dynamics are described by a two-dimensional partial integro-differential equation (PIDE) that incorporates both mixed derivative operators and a nonlocal integral component. The numerical strategy is constructed within the Radial Basis Function-Finite Difference (RBF–FD) paradigm, employing a modified multiquadric kernel that facilitates the analytical determination of localized differentiation weights. Spatial discretization is achieved through the method of lines, where the integral operator is approximated by means of quadrature rules combined with accurate interpolation techniques. Numerical experiments reveal that the proposed solver delivers superior accuracy while lowering computational cost in comparison with classical schemes.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"183 ","pages":"Article 106622"},"PeriodicalIF":4.1,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145884138","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}
引用次数: 0
A time-domain method of auxiliary sources for efficient analysis of transient electromagnetic scattering by moderately conductive cylinders 中等导电性圆柱瞬变电磁散射有效分析的辅助源时域方法
IF 4.1 2区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2026-02-01 Epub Date: 2025-12-29 DOI: 10.1016/j.enganabound.2025.106605
Minas Kouroublakis , Nikolaos L. Tsitsas , Yehuda Leviatan
This paper presents a time-domain implementation of the Method of Auxiliary Sources (MAS) combined with the Standard Impedance Boundary Condition (SIBC) for electromagnetic scattering problems involving cylindrical scatterers with finite but moderate conductivity. The proposed approach focuses on solving the two-dimensional problem using a first-order SIBC, which is valid when the conductivity is sufficiently higher than the maximum spectral frequency times the dielectric permittivity of the scatterer. This regime includes moderately conductive materials — such as carbon-based composites, conductive polymers, and doped dielectrics — that are increasingly used in real-world radio-frequency applications, including wearable electronics, electromagnetic interference shielding, and biomedical sensors. Under the above validity conditions, the interaction between the incident wave and the scatterer is dominated by surface effects, allowing for an efficient and accurate modeling strategy without the need to compute internal fields. The theoretical formulation of the time-domain MAS-SIBC method is developed, followed by extensive numerical testing on various geometries whose cross section is a closed curve. Such geometries include circular, elliptical, super-circular, rounded-triangular, and inverted-elliptical scatterers. A planar geometry is also tested. All results are validated against analytical solutions and commercial frequency-domain solvers, demonstrating the accuracy and practical potential of the proposed method. The findings suggest that time-domain MAS-SIBC offers a promising and computationally efficient approach for modeling scattering from materials even with moderate conductivity.
本文提出了一种结合标准阻抗边界条件(SIBC)的辅助源法(MAS)的时域实现方法,用于求解电导率有限但中等的圆柱形散射体的电磁散射问题。该方法着重于利用一阶SIBC求解二维问题,当电导率足够高于散射体的最大频谱频率乘以介电常数时,该方法有效。该体系包括中等导电性材料,如碳基复合材料、导电聚合物和掺杂电介质,这些材料越来越多地用于现实世界的射频应用,包括可穿戴电子产品、电磁干扰屏蔽和生物医学传感器。在上述有效条件下,入射波与散射体之间的相互作用主要由表面效应主导,因此无需计算内部场即可实现高效准确的建模策略。建立了时域MAS-SIBC方法的理论公式,并对各种截面为封闭曲线的几何形状进行了大量的数值试验。这些几何形状包括圆形、椭圆形、超圆形、圆三角形和倒椭圆形散射体。一个平面几何也进行了测试。所有结果都与解析解和商用频域解进行了验证,证明了所提出方法的准确性和实用潜力。研究结果表明,时域MAS-SIBC为中等导电性材料的散射建模提供了一种有前途的、计算效率高的方法。
{"title":"A time-domain method of auxiliary sources for efficient analysis of transient electromagnetic scattering by moderately conductive cylinders","authors":"Minas Kouroublakis ,&nbsp;Nikolaos L. Tsitsas ,&nbsp;Yehuda Leviatan","doi":"10.1016/j.enganabound.2025.106605","DOIUrl":"10.1016/j.enganabound.2025.106605","url":null,"abstract":"<div><div>This paper presents a time-domain implementation of the Method of Auxiliary Sources (MAS) combined with the Standard Impedance Boundary Condition (SIBC) for electromagnetic scattering problems involving cylindrical scatterers with finite but moderate conductivity. The proposed approach focuses on solving the two-dimensional problem using a first-order SIBC, which is valid when the conductivity is sufficiently higher than the maximum spectral frequency times the dielectric permittivity of the scatterer. This regime includes moderately conductive materials — such as carbon-based composites, conductive polymers, and doped dielectrics — that are increasingly used in real-world radio-frequency applications, including wearable electronics, electromagnetic interference shielding, and biomedical sensors. Under the above validity conditions, the interaction between the incident wave and the scatterer is dominated by surface effects, allowing for an efficient and accurate modeling strategy without the need to compute internal fields. The theoretical formulation of the time-domain MAS-SIBC method is developed, followed by extensive numerical testing on various geometries whose cross section is a closed curve. Such geometries include circular, elliptical, super-circular, rounded-triangular, and inverted-elliptical scatterers. A planar geometry is also tested. All results are validated against analytical solutions and commercial frequency-domain solvers, demonstrating the accuracy and practical potential of the proposed method. The findings suggest that time-domain MAS-SIBC offers a promising and computationally efficient approach for modeling scattering from materials even with moderate conductivity.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"183 ","pages":"Article 106605"},"PeriodicalIF":4.1,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145884137","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}
引用次数: 0
A high-order meshless method for the Allen-Cahn phase-field model Allen-Cahn相场模型的一种高阶无网格方法
IF 4.1 2区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2026-01-01 Epub Date: 2025-12-01 DOI: 10.1016/j.enganabound.2025.106559
Yuqian Xu, Wentao Ma
Accurately capturing interface dynamics in phase-field models remains a significant challenge, especially in the presence of narrow interfacial layers and geometrically complex domains. In this work, we propose a high-order meshless method based on the Generalized Finite Difference Method (GFDM) to efficiently solve the Allen-Cahn (AC) equation. The spatial discretization is constructed via Taylor expansions combined with moving least-squares approximations, achieving arbitrary-order accuracy without requiring mesh generation. Moreover, the method supports localized node refinement near interfaces and direct boundary treatment in irregular geometries, significantly enhancing computational efficiency and geometric adaptability. A Crank-Nicolson (CN) scheme is employed for time discretization to preserve the energy dissipation property of the model. Extensive numerical experiments, including mean curvature flow and phase separation in both regular and irregular domains, demonstrate the proposed method’s accuracy, energy stability, and geometric adaptability.
在相场模型中准确捕获界面动力学仍然是一个重大挑战,特别是在存在窄界面层和几何复杂区域的情况下。在本文中,我们提出了一种基于广义有限差分法(GFDM)的高阶无网格方法来有效地求解Allen-Cahn (AC)方程。空间离散化是通过泰勒展开结合移动最小二乘逼近,实现任意阶精度而不需要网格生成。此外,该方法支持界面附近的局部节点细化和不规则几何的直接边界处理,显著提高了计算效率和几何适应性。采用Crank-Nicolson (CN)格式进行时间离散,以保持模型的能量耗散特性。大量的数值实验,包括规则和不规则域的平均曲率流动和相分离,证明了该方法的准确性、能量稳定性和几何适应性。
{"title":"A high-order meshless method for the Allen-Cahn phase-field model","authors":"Yuqian Xu,&nbsp;Wentao Ma","doi":"10.1016/j.enganabound.2025.106559","DOIUrl":"10.1016/j.enganabound.2025.106559","url":null,"abstract":"<div><div>Accurately capturing interface dynamics in phase-field models remains a significant challenge, especially in the presence of narrow interfacial layers and geometrically complex domains. In this work, we propose a high-order meshless method based on the Generalized Finite Difference Method (GFDM) to efficiently solve the Allen-Cahn (AC) equation. The spatial discretization is constructed via Taylor expansions combined with moving least-squares approximations, achieving arbitrary-order accuracy without requiring mesh generation. Moreover, the method supports localized node refinement near interfaces and direct boundary treatment in irregular geometries, significantly enhancing computational efficiency and geometric adaptability. A Crank-Nicolson (CN) scheme is employed for time discretization to preserve the energy dissipation property of the model. Extensive numerical experiments, including mean curvature flow and phase separation in both regular and irregular domains, demonstrate the proposed method’s accuracy, energy stability, and geometric adaptability.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"182 ","pages":"Article 106559"},"PeriodicalIF":4.1,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145657539","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}
引用次数: 0
Semi-analytical Discrete Green’s Function Method for cathodic protection problems 阴极保护问题的半解析离散格林函数法
IF 4.1 2区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2026-01-01 Epub Date: 2025-11-26 DOI: 10.1016/j.enganabound.2025.106544
Kevin D. Cole , Barış Çetin , Barbaros Çetin
Cathodic protection (CP) is a vital technique to prevent corrosion in submerged or intermittently exposed metallic structures. While numerical methods such as finite element and boundary element methods are commonly used to model CP systems, they often involve high computational costs, especially under nonlinear boundary conditions. This paper introduces a semi-analytical approach based on the Discrete Green’s Function Method (DGFM) for solving cathodic protection problems with both linear and nonlinear boundary conditions. The DGFM combines the accuracy of analytical methods with the flexibility of numerical techniques, enabling efficient and precise solutions. A graph theory-based framework is employed to generate the discrete Green’s functions from the Laplacian matrix, allowing the electric potential distribution to be calculated through matrix multiplications. The method is verified against analytical solutions for a linear test problem and verified through comparison with COMSOL Multiphysics simulations for a nonlinear CP model involving realistic polarization curves for zinc and steel. Results demonstrate that DGFM significantly reduces computational time while maintaining high accuracy, making it a promising tool for CP design and optimization.
阴极保护(CP)是防止水下或间歇暴露金属结构腐蚀的重要技术。虽然有限元和边界元等数值方法通常用于模拟CP系统,但它们往往涉及较高的计算成本,特别是在非线性边界条件下。本文介绍了一种基于离散格林函数法的半解析方法,用于求解具有线性和非线性边界条件的阴极保护问题。DGFM结合了分析方法的准确性和数值技术的灵活性,实现了高效和精确的解决方案。采用基于图论的框架,由拉普拉斯矩阵生成离散格林函数,通过矩阵乘法计算电势分布。通过一个线性测试问题的解析解验证了该方法的有效性,并通过COMSOL Multiphysics仿真对锌和钢的实际极化曲线非线性CP模型进行了验证。结果表明,DGFM显著减少了计算时间,同时保持了较高的精度,使其成为CP设计和优化的一个有前途的工具。
{"title":"Semi-analytical Discrete Green’s Function Method for cathodic protection problems","authors":"Kevin D. Cole ,&nbsp;Barış Çetin ,&nbsp;Barbaros Çetin","doi":"10.1016/j.enganabound.2025.106544","DOIUrl":"10.1016/j.enganabound.2025.106544","url":null,"abstract":"<div><div>Cathodic protection (CP) is a vital technique to prevent corrosion in submerged or intermittently exposed metallic structures. While numerical methods such as finite element and boundary element methods are commonly used to model CP systems, they often involve high computational costs, especially under nonlinear boundary conditions. This paper introduces a semi-analytical approach based on the Discrete Green’s Function Method (DGFM) for solving cathodic protection problems with both linear and nonlinear boundary conditions. The DGFM combines the accuracy of analytical methods with the flexibility of numerical techniques, enabling efficient and precise solutions. A graph theory-based framework is employed to generate the discrete Green’s functions from the Laplacian matrix, allowing the electric potential distribution to be calculated through matrix multiplications. The method is verified against analytical solutions for a linear test problem and verified through comparison with COMSOL Multiphysics simulations for a nonlinear CP model involving realistic polarization curves for zinc and steel. Results demonstrate that DGFM significantly reduces computational time while maintaining high accuracy, making it a promising tool for CP design and optimization.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"182 ","pages":"Article 106544"},"PeriodicalIF":4.1,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145598942","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}
引用次数: 0
Coupled role of mineral heterogeneity and cross fractures in the macro–meso mechanical behavior of granite: Insights from image-informed numerical modeling 矿物非均质性和交叉裂缝在花岗岩宏观细观力学行为中的耦合作用:来自图像信息数值模拟的见解
IF 4.1 2区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2026-01-01 Epub Date: 2025-11-29 DOI: 10.1016/j.enganabound.2025.106564
Tingting Liu , Shenghao Yang , Luyang Ding , Xiaohan Xie , Hui Shen , Xinping Li
Tectonic evolution induces discontinuities and mineral heterogeneity in rock masses, complicating mechanical behavior. This study examines how mineral composition affects the response of granite with cross fractures, a key issue for underground engineering stability. The discrete element method (DEM) was implemented using the particle flow code in two dimensions (PFC2D). Parameters were calibrated through a digital image processing grain-based model (DIP–GBM), on which a random distribution properties grain-based model (RDP–GBM) was established to simulate mechanical response, fracture evolution, and energy characteristics under varied mineral compositions. Results show that fracture connectivity (ω) is the primary control: uniaxial compressive strength (UCS) decreases with increasing connectivity, while the elastic modulus (E) increases slightly at ω = 0.141 before declining. Mineral composition exerts a secondary influence, with higher biotite content and feldspar-to-quartz (λ) ratios reducing strength. In intact granite, failure is governed by mineral heterogeneity, whereas in cross-fractured granite, fracture geometry dominates. Low-stress region exists in granite with cross fractures and is linked to crack evolution. With increasing connectivity, total and elastic strain energy decrease, while biotite-rich rocks show reduced total energy. These findings highlight the coupled effects of fractures and minerals and provide a basis for stability evaluation in underground engineering.
构造演化导致岩体中的不连续和矿物非均质性,使力学行为复杂化。本文研究了矿物成分如何影响具有交叉裂缝的花岗岩的响应,这是地下工程稳定性的一个关键问题。采用二维粒子流代码(PFC2D)实现离散元法(DEM)。通过数字图像处理颗粒模型(DIP-GBM)标定参数,在此基础上建立随机分布属性颗粒模型(RDP-GBM),模拟不同矿物成分下的力学响应、裂缝演化和能量特征。结果表明,裂缝连通性(ω)是主要控制因素,单轴抗压强度(UCS)随连通性的增加而降低,而弹性模量(E)在ω = 0.141时略有增加,然后下降。矿物组成是次要影响,较高的黑云母含量和长石与石英(λ)比降低了强度。在完整的花岗岩中,破坏是由矿物非均质性决定的,而在交叉断裂的花岗岩中,裂缝几何形状占主导地位。低应力区存在于具有交叉裂缝的花岗岩中,与裂缝演化有关。随着连通性的增加,总应变能和弹性应变能降低,而富含黑云母的岩石总应变能降低。这些发现突出了裂缝与矿物的耦合作用,为地下工程稳定性评价提供了依据。
{"title":"Coupled role of mineral heterogeneity and cross fractures in the macro–meso mechanical behavior of granite: Insights from image-informed numerical modeling","authors":"Tingting Liu ,&nbsp;Shenghao Yang ,&nbsp;Luyang Ding ,&nbsp;Xiaohan Xie ,&nbsp;Hui Shen ,&nbsp;Xinping Li","doi":"10.1016/j.enganabound.2025.106564","DOIUrl":"10.1016/j.enganabound.2025.106564","url":null,"abstract":"<div><div>Tectonic evolution induces discontinuities and mineral heterogeneity in rock masses, complicating mechanical behavior. This study examines how mineral composition affects the response of granite with cross fractures, a key issue for underground engineering stability. The discrete element method (DEM) was implemented using the particle flow code in two dimensions (PFC2D). Parameters were calibrated through a digital image processing grain-based model (DIP–GBM), on which a random distribution properties grain-based model (RDP–GBM) was established to simulate mechanical response, fracture evolution, and energy characteristics under varied mineral compositions. Results show that fracture connectivity (<em>ω</em>) is the primary control: uniaxial compressive strength (UCS) decreases with increasing connectivity, while the elastic modulus (<em>E</em>) increases slightly at <em>ω</em> = 0.141 before declining. Mineral composition exerts a secondary influence, with higher biotite content and feldspar-to-quartz (<em>λ</em>) ratios reducing strength. In intact granite, failure is governed by mineral heterogeneity, whereas in cross-fractured granite, fracture geometry dominates. Low-stress region exists in granite with cross fractures and is linked to crack evolution. With increasing connectivity, total and elastic strain energy decrease, while biotite-rich rocks show reduced total energy. These findings highlight the coupled effects of fractures and minerals and provide a basis for stability evaluation in underground engineering.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"182 ","pages":"Article 106564"},"PeriodicalIF":4.1,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145619675","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}
引用次数: 0
A Feature-Enhanced Physics Informed Neural Network for trajectory simulation of charged particles 用于带电粒子轨迹模拟的特征增强物理信息神经网络
IF 4.1 2区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2026-01-01 Epub Date: 2025-11-29 DOI: 10.1016/j.enganabound.2025.106562
Fan Yang, Xuan Liu, Pengbo Wang, Xinheng Li
Trajectory simulation of charged particles is a complex multiphysics problem with dynamic field involved. With the mesh-free essence, Physics-Informed Neural Networks (PINNs) offers high capability in conformal modeling and dynamic field mapping. By integrating governing physical laws with real-time or experimental data, PINNs can accurately capture dynamic field variation and system response, demonstrating substantial potential for field–particle coupling simulations. However, conventional PINNs still face challenges caused by space charge effect and multi-scale field variations. This paper proposes a Feature-Enhanced PINN (FE-PINN) framework for the simulation of charged particle dynamics. For electromagnetic fields with evolving internal sources and complex geometries, Feature-Enhanced PINN eliminates the need for network retraining under varying space charge distributions, and applies targeted enhancement in collocation strategies and network structure to improve convergence and accuracy in domains of localized high-gradient. Built upon FE-PINN electromagnetic field solutions, particle trajectory simulation is achieved by iteratively solving the Poisson’s equation and particle motion equation under a fixed magnetic field. The proposed method is validated using the magnetron injection gun (MIG) of an 800 GHz gyrotron with results compared with that of CST Studio Suite.
带电粒子的轨迹模拟是一个涉及动力学场的复杂多物理场问题。物理信息神经网络(PINNs)具有无网格的本质,在保形建模和动态场映射方面具有很高的能力。通过将控制物理定律与实时或实验数据相结合,pinn可以准确地捕获动态场变化和系统响应,显示出场-粒子耦合模拟的巨大潜力。然而,传统的pin - ns仍然面临着空间电荷效应和多尺度场变化带来的挑战。本文提出了一种特征增强的pin - n (fe - pin)框架,用于模拟带电粒子动力学。对于具有演化的内源和复杂几何形状的电磁场,Feature-Enhanced PINN消除了在不同空间电荷分布下对网络进行再训练的需要,并对配置策略和网络结构进行了有针对性的增强,以提高局部高梯度域的收敛性和准确性。在FE-PINN电磁场解的基础上,通过迭代求解固定磁场下的泊松方程和粒子运动方程实现粒子轨迹仿真。利用800 GHz回旋管的磁控管注射枪(MIG)对该方法进行了验证,并与CST Studio Suite的结果进行了比较。
{"title":"A Feature-Enhanced Physics Informed Neural Network for trajectory simulation of charged particles","authors":"Fan Yang,&nbsp;Xuan Liu,&nbsp;Pengbo Wang,&nbsp;Xinheng Li","doi":"10.1016/j.enganabound.2025.106562","DOIUrl":"10.1016/j.enganabound.2025.106562","url":null,"abstract":"<div><div>Trajectory simulation of charged particles is a complex multiphysics problem with dynamic field involved. With the mesh-free essence, Physics-Informed Neural Networks (PINNs) offers high capability in conformal modeling and dynamic field mapping. By integrating governing physical laws with real-time or experimental data, PINNs can accurately capture dynamic field variation and system response, demonstrating substantial potential for field–particle coupling simulations. However, conventional PINNs still face challenges caused by space charge effect and multi-scale field variations. This paper proposes a Feature-Enhanced PINN (FE-PINN) framework for the simulation of charged particle dynamics. For electromagnetic fields with evolving internal sources and complex geometries, Feature-Enhanced PINN eliminates the need for network retraining under varying space charge distributions, and applies targeted enhancement in collocation strategies and network structure to improve convergence and accuracy in domains of localized high-gradient. Built upon FE-PINN electromagnetic field solutions, particle trajectory simulation is achieved by iteratively solving the Poisson’s equation and particle motion equation under a fixed magnetic field. The proposed method is validated using the magnetron injection gun (MIG) of an 800 GHz gyrotron with results compared with that of CST Studio Suite.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"182 ","pages":"Article 106562"},"PeriodicalIF":4.1,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145613985","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}
引用次数: 0
A 2.5D BEM-based approach in the Bézier–Bernstein space for railway noise prediction and acoustic barrier assessment 基于bsamier - bernstein空间的铁路噪声预测和声障评价的2.5D bem方法
IF 4.1 2区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2026-01-01 Epub Date: 2025-11-27 DOI: 10.1016/j.enganabound.2025.106568
R. Velázquez-Mata , C. Knuth , A. Romero , G. Squicciarini , A. Tadeu , D.J. Thompson , P. Galvín
Noise pollution from railway traffic, primarily caused by rolling noise resulting from the vibrations of the track and wheels, is a major public health concern. While traditional acoustic barriers are effective, they are often visually intrusive, particularly in urban settings. This has led to growing interest in more integrated solutions, such as low, close barriers, which require accurate noise prediction tools. This paper presents a two-and-a-half-dimensional BEM for predicting and mitigating railway noise. The method uses Bézier–Bernstein space to accurately model complex geometries, enhancing noise prediction across different rail profiles. Several rail configurations are compared to evaluate their impact on noise emissions and to support the design of more effective and adaptable barrier solutions. The method is then applied to evaluate the performance of a specific low-height barrier configuration, considering the presence of the vehicle to assess its impact on noise reduction. Numerical predictions are validated through comparison with experimental data and other numerical approaches. Results highlight the importance of accurate source modelling for barrier design and demonstrate the potential of the proposed method as a flexible tool for developing noise mitigation solutions that utilize the barrier’s geometry to improve acoustic performance and support visual integration in urban environments.
铁路交通的噪音污染是一个重大的公共卫生问题,主要是由轨道和车轮振动产生的滚动噪音造成的。虽然传统的隔音屏障是有效的,但它们往往是视觉上的干扰,特别是在城市环境中。这导致人们对更集成的解决方案越来越感兴趣,例如需要精确噪声预测工具的低、封闭屏障。本文提出了一种用于铁路噪声预测和抑制的二维半边界元法。该方法使用bsamzier - bernstein空间来精确模拟复杂的几何形状,增强了不同轨道轮廓的噪声预测。对几种轨道配置进行比较,以评估其对噪声排放的影响,并支持设计更有效和适应性更强的屏障解决方案。然后将该方法应用于评估特定低高度屏障配置的性能,考虑车辆的存在以评估其对降噪的影响。通过与实验数据和其他数值方法的比较,验证了数值预测的正确性。结果强调了准确的声源建模对屏障设计的重要性,并证明了所提出的方法作为开发噪声缓解解决方案的灵活工具的潜力,这些解决方案利用屏障的几何形状来改善声学性能并支持城市环境中的视觉整合。
{"title":"A 2.5D BEM-based approach in the Bézier–Bernstein space for railway noise prediction and acoustic barrier assessment","authors":"R. Velázquez-Mata ,&nbsp;C. Knuth ,&nbsp;A. Romero ,&nbsp;G. Squicciarini ,&nbsp;A. Tadeu ,&nbsp;D.J. Thompson ,&nbsp;P. Galvín","doi":"10.1016/j.enganabound.2025.106568","DOIUrl":"10.1016/j.enganabound.2025.106568","url":null,"abstract":"<div><div>Noise pollution from railway traffic, primarily caused by rolling noise resulting from the vibrations of the track and wheels, is a major public health concern. While traditional acoustic barriers are effective, they are often visually intrusive, particularly in urban settings. This has led to growing interest in more integrated solutions, such as low, close barriers, which require accurate noise prediction tools. This paper presents a two-and-a-half-dimensional BEM for predicting and mitigating railway noise. The method uses Bézier–Bernstein space to accurately model complex geometries, enhancing noise prediction across different rail profiles. Several rail configurations are compared to evaluate their impact on noise emissions and to support the design of more effective and adaptable barrier solutions. The method is then applied to evaluate the performance of a specific low-height barrier configuration, considering the presence of the vehicle to assess its impact on noise reduction. Numerical predictions are validated through comparison with experimental data and other numerical approaches. Results highlight the importance of accurate source modelling for barrier design and demonstrate the potential of the proposed method as a flexible tool for developing noise mitigation solutions that utilize the barrier’s geometry to improve acoustic performance and support visual integration in urban environments.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"182 ","pages":"Article 106568"},"PeriodicalIF":4.1,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145611796","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}
引用次数: 0
Dual integral formulation for acoustic resonances of a two-dimensional damped cavity with a thinpartition 具有薄隔板的二维阻尼腔的声学共振的对偶积分公式
IF 4.1 2区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2026-01-01 Epub Date: 2025-11-27 DOI: 10.1016/j.enganabound.2025.106560
Kue-Hong Chen , Yi-Kui Liu , Jeng-Hong Kao , Jeng-Tzong Chen
This paper derives a dual integral formulation for the Helmholtz equation used to analyze the problem of the acoustic resonances for a damped cavity with a thin, incomplete partition. The presence of damping in the medium results in the incorporation of complex wave numbers in the Helmholtz equation within mathematical models. To address the singular and hypersingular integrals of the four complex-valued wave number kernel functions, the addition theorem is employed to expand them into series representations involving only real wave number. After then, the singular and hypersingular integrals are successfully transformed into a summation of regular integrals in an infinite series using the proposed regularization technique. This study investigates the acoustic resonances of a cavity with a thin partition, focusing on the effect of damping on these resonances. A case with an exact solution is given as a standard benchmark to evaluate the convergence and accuracy of the proposed dual integral equation. Finally, two numerical examples, involving rectangular and submarine-shaped cavities with a thin, incomplete partition, are presented. The acoustic frequencies and modes are determined using the developed program.
本文导出了用于分析具有薄不完全隔板的阻尼腔的声学共振问题的亥姆霍兹方程的对偶积分公式。介质中阻尼的存在导致在数学模型中亥姆霍兹方程中加入复波数。为了求解四个复值波数核函数的奇异积分和超奇异积分,利用加法定理将它们展开为只涉及实波数的级数表示。然后,利用所提出的正则化技术,将奇异积分和超奇异积分成功地转化为无穷级数上的正则积分和。本文研究了具有薄隔板的腔体的声学共振,重点研究了阻尼对这些共振的影响。给出了一个精确解的例子,作为评价对偶积分方程收敛性和准确性的标准基准。最后,给出了矩形腔和不完全薄隔板的潜艇形腔的两个数值算例。使用所开发的程序确定声频率和模态。
{"title":"Dual integral formulation for acoustic resonances of a two-dimensional damped cavity with a thinpartition","authors":"Kue-Hong Chen ,&nbsp;Yi-Kui Liu ,&nbsp;Jeng-Hong Kao ,&nbsp;Jeng-Tzong Chen","doi":"10.1016/j.enganabound.2025.106560","DOIUrl":"10.1016/j.enganabound.2025.106560","url":null,"abstract":"<div><div>This paper derives a dual integral formulation for the Helmholtz equation used to analyze the problem of the acoustic resonances for a damped cavity with a thin, incomplete partition. The presence of damping in the medium results in the incorporation of complex wave numbers in the Helmholtz equation within mathematical models. To address the singular and hypersingular integrals of the four complex-valued wave number kernel functions, the addition theorem is employed to expand them into series representations involving only real wave number. After then, the singular and hypersingular integrals are successfully transformed into a summation of regular integrals in an infinite series using the proposed regularization technique. This study investigates the acoustic resonances of a cavity with a thin partition, focusing on the effect of damping on these resonances. A case with an exact solution is given as a standard benchmark to evaluate the convergence and accuracy of the proposed dual integral equation. Finally, two numerical examples, involving rectangular and submarine-shaped cavities with a thin, incomplete partition, are presented. The acoustic frequencies and modes are determined using the developed program.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"182 ","pages":"Article 106560"},"PeriodicalIF":4.1,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145611795","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}
引用次数: 0
An improved Tchebychev-radial point interpolation method for large deformation analysis of hyperelastic model 一种用于超弹性模型大变形分析的改进tchebychevv -径向点插值方法
IF 4.1 2区 工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY Pub Date : 2026-01-01 Epub Date: 2025-12-01 DOI: 10.1016/j.enganabound.2025.106566
Nha Thanh Nguyen, Thai Van Vu, Pham Toan Thang
In this study, an improvement for the Tchebychev-radial point interpolation method (TRPIM) is proposed for the finite deformation analysis of the hyperelastic structure. The TRPIM shape function is developed by integrating the radial basis function with the Tchebychev polynomial basis function instead of traditional polynomial basis functions, and this combination produces a more accurate interpolation. Both the first and second kinds of Tchebychev polynomial are tested and evaluated for the proposed method. The modification revolves around integrating the Cartesian transformation method into the TRPIM approach to make it a truly meshfree method that does not use background cells system for numerical integration. The nonlinear behavior of hyperelastic models (Neo-Hookean, Mooney–Rivlin, and Ogden) under a finite deformation state is simulated with the total Lagrange formulation and standard Newton–Raphson algorithm. To assess the large deformation behavior of two-dimensional hyperelastic problems, a series of numerical tests are carried out. The reliability and performance of the proposed method is confirmed by comparing its results with those of available solutions.
本文提出了一种改进的tchebychevv -radial点插值法(TRPIM),用于超弹性结构的有限变形分析。将径向基函数与Tchebychev多项式基函数相结合,取代了传统的多项式基函数,得到了TRPIM形状函数,这种组合可以实现更精确的插值。对所提出的方法进行了第一类和第二类切比切夫多项式的检验和评价。改进的核心是将笛卡尔变换方法与TRPIM方法相结合,使其成为一种真正的无网格方法,不使用背景单元系统进行数值积分。采用全拉格朗日公式和标准Newton-Raphson算法模拟了有限变形状态下的超弹性模型(Neo-Hookean、Mooney-Rivlin和Ogden)的非线性行为。为了评估二维超弹性问题的大变形行为,进行了一系列数值试验。通过与已有解的结果比较,验证了所提方法的可靠性和性能。
{"title":"An improved Tchebychev-radial point interpolation method for large deformation analysis of hyperelastic model","authors":"Nha Thanh Nguyen,&nbsp;Thai Van Vu,&nbsp;Pham Toan Thang","doi":"10.1016/j.enganabound.2025.106566","DOIUrl":"10.1016/j.enganabound.2025.106566","url":null,"abstract":"<div><div>In this study, an improvement for the Tchebychev-radial point interpolation method (TRPIM) is proposed for the finite deformation analysis of the hyperelastic structure. The TRPIM shape function is developed by integrating the radial basis function with the Tchebychev polynomial basis function instead of traditional polynomial basis functions, and this combination produces a more accurate interpolation. Both the first and second kinds of Tchebychev polynomial are tested and evaluated for the proposed method. The modification revolves around integrating the Cartesian transformation method into the TRPIM approach to make it a truly meshfree method that does not use background cells system for numerical integration. The nonlinear behavior of hyperelastic models (Neo-Hookean, Mooney–Rivlin, and Ogden) under a finite deformation state is simulated with the total Lagrange formulation and standard Newton–Raphson algorithm. To assess the large deformation behavior of two-dimensional hyperelastic problems, a series of numerical tests are carried out. The reliability and performance of the proposed method is confirmed by comparing its results with those of available solutions.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"182 ","pages":"Article 106566"},"PeriodicalIF":4.1,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145651032","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}
引用次数: 0
期刊
Engineering Analysis with Boundary Elements
全部 Acc. Chem. Res. ACS Applied Bio Materials ACS Appl. Electron. Mater. ACS Appl. Energy Mater. ACS Appl. Mater. Interfaces ACS Appl. Nano Mater. ACS Appl. Polym. Mater. ACS BIOMATER-SCI ENG ACS Catal. ACS Cent. Sci. ACS Chem. Biol. ACS Chemical Health & Safety ACS Chem. Neurosci. ACS Comb. Sci. ACS Earth Space Chem. ACS Energy Lett. ACS Infect. Dis. ACS Macro Lett. ACS Mater. Lett. ACS Med. Chem. Lett. ACS Nano ACS Omega ACS Photonics ACS Sens. ACS Sustainable Chem. Eng. ACS Synth. Biol. Anal. Chem. BIOCHEMISTRY-US Bioconjugate Chem. BIOMACROMOLECULES Chem. Res. Toxicol. Chem. Rev. Chem. Mater. CRYST GROWTH DES ENERG FUEL Environ. Sci. Technol. Environ. Sci. Technol. Lett. Eur. J. Inorg. Chem. IND ENG CHEM RES Inorg. Chem. J. Agric. Food. Chem. J. Chem. Eng. Data J. Chem. Educ. J. Chem. Inf. Model. J. Chem. Theory Comput. J. Med. Chem. J. Nat. Prod. J PROTEOME RES J. Am. Chem. Soc. LANGMUIR MACROMOLECULES Mol. Pharmaceutics Nano Lett. Org. Lett. ORG PROCESS RES DEV ORGANOMETALLICS J. Org. Chem. J. Phys. Chem. J. Phys. Chem. A J. Phys. Chem. B J. Phys. Chem. C J. Phys. Chem. Lett. Analyst Anal. Methods Biomater. Sci. Catal. Sci. Technol. Chem. Commun. Chem. Soc. Rev. CHEM EDUC RES PRACT CRYSTENGCOMM Dalton Trans. Energy Environ. Sci. ENVIRON SCI-NANO ENVIRON SCI-PROC IMP ENVIRON SCI-WAT RES Faraday Discuss. Food Funct. Green Chem. Inorg. Chem. Front. Integr. Biol. J. Anal. At. Spectrom. J. Mater. Chem. A J. Mater. Chem. B J. Mater. Chem. C Lab Chip Mater. Chem. Front. Mater. Horiz. MEDCHEMCOMM Metallomics Mol. Biosyst. Mol. Syst. Des. Eng. Nanoscale Nanoscale Horiz. Nat. Prod. Rep. New J. Chem. Org. Biomol. Chem. Org. Chem. Front. PHOTOCH PHOTOBIO SCI PCCP Polym. Chem.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:604180095
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1