Pub Date : 2026-03-01Epub Date: 2026-01-19DOI: 10.1016/j.enganabound.2026.106648
Ruiping Niu, Yi Cai, Chengtao Wu
Computer-assisted thermal ablation therapy techniques depend on precise and effective soft tissue temperature prediction. This paper proposes an accurate explicit dynamics face-based smoothed finite element algorithm (FS-FEM) for thermoelastic analysis of soft tissue, grounded in large deformation thermoselasticity and total Lagrangian formulation. A tightly coupled model is proposed to capture the interactive behavior: (i) bioheat transfer with tissue deformation, (ii) tissue deformation due to thermal expansion. Then the weakened weak form of the presented coupled model is derived, and the explicit dynamic face-based smoothed finite element approach is formulated. Finally, the effectiveness and compatibility of the established methodology are verified using a human liver modelling with publicly available CT scan datasets to illustrate a clinically pertinent scenario of thermal ablation for hepatocellular carcinoma. The results demonstrate that our proposed numerical algorithm efficiently solves the simulation of liver thermal ablation, with temperature and tissue deformation predictions being more accurate than those obtained by the finite element method (FEM).
{"title":"An explicit dynamic face-based smoothed finite element approach to thermoelastic modeling in thermal ablation therapy","authors":"Ruiping Niu, Yi Cai, Chengtao Wu","doi":"10.1016/j.enganabound.2026.106648","DOIUrl":"10.1016/j.enganabound.2026.106648","url":null,"abstract":"<div><div>Computer-assisted thermal ablation therapy techniques depend on precise and effective soft tissue temperature prediction. This paper proposes an accurate explicit dynamics face-based smoothed finite element algorithm (FS-FEM) for thermoelastic analysis of soft tissue, grounded in large deformation thermoselasticity and total Lagrangian formulation. A tightly coupled model is proposed to capture the interactive behavior: (i) bioheat transfer with tissue deformation, (ii) tissue deformation due to thermal expansion. Then the weakened weak form of the presented coupled model is derived, and the explicit dynamic face-based smoothed finite element approach is formulated. Finally, the effectiveness and compatibility of the established methodology are verified using a human liver modelling with publicly available CT scan datasets to illustrate a clinically pertinent scenario of thermal ablation for hepatocellular carcinoma. The results demonstrate that our proposed numerical algorithm efficiently solves the simulation of liver thermal ablation, with temperature and tissue deformation predictions being more accurate than those obtained by the finite element method (FEM).</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"184 ","pages":"Article 106648"},"PeriodicalIF":4.1,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146000869","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}
Pub Date : 2026-03-01Epub Date: 2026-01-06DOI: 10.1016/j.enganabound.2025.106612
Jiaxin Li , Shuihuai Yang , Xing Wei , Linlin Sun , Yue Yu
This paper presents the first attempt to apply the generalized finite difference method (GFDM) for the plane strain analysis of magneto-electro-elastic (MEE) materials. The computational domain is discretized into a set of nodes, where each node is associated with a local influence domain comprising its neighboring nodes. The variables at the nodes within each local influence domain are approximated using Taylor series expansion. A moving least squares method is then employed to establish a relationship between the partial derivatives of the variables and their values at the nodes in the local domain. Finally, the governing equations and boundary conditions are transformed into a system of linear equations. To improve the property of the coefficient matrix, each linear equation is divided by the maximum material parameter present in that equation. The feasibility and accuracy of the GFDM are validated through three test cases, with results compared against those from the meshless local Petrov–Galerkin, radial point interpolation, local radial basis function, boundary element, and finite element methods.
{"title":"Analysis of plane problems in magneto-electro-elastic media using the generalized finite difference method","authors":"Jiaxin Li , Shuihuai Yang , Xing Wei , Linlin Sun , Yue Yu","doi":"10.1016/j.enganabound.2025.106612","DOIUrl":"10.1016/j.enganabound.2025.106612","url":null,"abstract":"<div><div>This paper presents the first attempt to apply the generalized finite difference method (GFDM) for the plane strain analysis of magneto-electro-elastic (MEE) materials. The computational domain is discretized into a set of nodes, where each node is associated with a local influence domain comprising its neighboring nodes. The variables at the nodes within each local influence domain are approximated using Taylor series expansion. A moving least squares method is then employed to establish a relationship between the partial derivatives of the variables and their values at the nodes in the local domain. Finally, the governing equations and boundary conditions are transformed into a system of linear equations. To improve the property of the coefficient matrix, each linear equation is divided by the maximum material parameter present in that equation. The feasibility and accuracy of the GFDM are validated through three test cases, with results compared against those from the meshless local Petrov–Galerkin, radial point interpolation, local radial basis function, boundary element, and finite element methods.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"184 ","pages":"Article 106612"},"PeriodicalIF":4.1,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145904246","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}
Pub Date : 2026-03-01Epub Date: 2026-01-17DOI: 10.1016/j.enganabound.2026.106642
Weichen Sun , Qiang Xie , Kai Wu , Zhilin Cao , Hexing Zhang , Xiang Fu , Yuxin Ban
This paper proposes a novel Anisotropic Lattice Spring Model (ALSM) to address the limitations of traditional LSMs in simulating anisotropic deformation of layered rocks such as schist, shale, and slate. The model introduces two key innovations: (1) normal-tangential coupled spring bonds that form a full stiffness matrix with non-zero off-diagonal terms, overcoming the inherent Poisson's ratio constraints of LSMs; and (2) a dual-constraint stiffness matching method that incorporates stress symmetry into the energy equivalence framework, ensuring physical consistency in macroscopic constitutive relations and enabling prescise mapping between microscopic parameters and macroscopic elastic constants. Verification results show that the ALSM significantly broadens the achievable range of Poisson's ratios and accurately simulates strongly anisotropic rocks. In uniaxial compression and Brazilian disc tests, the ALSM yields apparent elastic moduli consistent with theory/experiment, captures non-monotonic anisotropy, and markedly reduces displacement/stress errors compared to traditional models. The model allows direct input of macroscopic anisotropic parameters, avoiding complex microscopic calibration. The established elastic homogenization framework supports further study of anisotropic fracture and offers a new approach for bottom-up design of anisotropic materials.
{"title":"A novel anisotropic lattice spring model for elastically-homogeneous modelling of layered rocks","authors":"Weichen Sun , Qiang Xie , Kai Wu , Zhilin Cao , Hexing Zhang , Xiang Fu , Yuxin Ban","doi":"10.1016/j.enganabound.2026.106642","DOIUrl":"10.1016/j.enganabound.2026.106642","url":null,"abstract":"<div><div>This paper proposes a novel Anisotropic Lattice Spring Model (ALSM) to address the limitations of traditional LSMs in simulating anisotropic deformation of layered rocks such as schist, shale, and slate. The model introduces two key innovations: (1) normal-tangential coupled spring bonds that form a full stiffness matrix with non-zero off-diagonal terms, overcoming the inherent Poisson's ratio constraints of LSMs; and (2) a dual-constraint stiffness matching method that incorporates stress symmetry into the energy equivalence framework, ensuring physical consistency in macroscopic constitutive relations and enabling prescise mapping between microscopic parameters and macroscopic elastic constants. Verification results show that the ALSM significantly broadens the achievable range of Poisson's ratios and accurately simulates strongly anisotropic rocks. In uniaxial compression and Brazilian disc tests, the ALSM yields apparent elastic moduli consistent with theory/experiment, captures non-monotonic anisotropy, and markedly reduces displacement/stress errors compared to traditional models. The model allows direct input of macroscopic anisotropic parameters, avoiding complex microscopic calibration. The established elastic homogenization framework supports further study of anisotropic fracture and offers a new approach for bottom-up design of anisotropic materials.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"184 ","pages":"Article 106642"},"PeriodicalIF":4.1,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145995149","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}
Pub Date : 2026-03-01Epub Date: 2026-01-14DOI: 10.1016/j.enganabound.2026.106643
Anjali Singh, Rajen Kumar Sinha
The conservative Allen–Cahn (CAC) equation is a second-order nonlinear partial differential equation that models phase separation in binary mixtures while preserving the total volume. Physics-informed neural networks (PINNs) have demonstrated considerable success in approximating solutions to various classes of partial differential equations; however, their application to CAC models remains challenging. These difficulties stem from the presence of a small interfacial parameter in the nonlinear term and highly nonlinear mass-correction terms , , which significantly degrade the approximation accuracy and mass conservation properties of standard PINNs. In this work, we propose a novel hybrid mass-constrained physics-informed neural network (Mc-PINN) framework for efficiently and accurately solving three types of CAC models in convex polygonal domains. The proposed method integrates deep learning with operator-splitting strategies to decompose the original CAC equations into simpler subproblems. One subproblem admits an analytical solution, while the other is solved using the Mc-PINN scheme. To further enhance efficiency, a Metropolis–Hastings based adaptive sampling strategy is introduced. In addition, we derive error estimates for the proposed method applied to all three CAC models. Numerical experiments demonstrate the robustness, accuracy, and mass-conservation capability of the proposed framework.
{"title":"Mc-PINN for solving Conservative Allen–Cahn equations using Metropolis–Hasting based sampling","authors":"Anjali Singh, Rajen Kumar Sinha","doi":"10.1016/j.enganabound.2026.106643","DOIUrl":"10.1016/j.enganabound.2026.106643","url":null,"abstract":"<div><div>The conservative Allen–Cahn (CAC) equation is a second-order nonlinear partial differential equation that models phase separation in binary mixtures while preserving the total volume. Physics-informed neural networks (PINNs) have demonstrated considerable success in approximating solutions to various classes of partial differential equations; however, their application to CAC models remains challenging. These difficulties stem from the presence of a small interfacial parameter <span><math><mi>ϵ</mi></math></span> in the nonlinear term <span><math><mrow><msup><mrow><mi>F</mi></mrow><mrow><mo>′</mo></mrow></msup><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow></mrow></math></span> and highly nonlinear mass-correction terms <span><math><msub><mrow><mi>G</mi></mrow><mrow><mi>i</mi></mrow></msub></math></span>, <span><math><mrow><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>3</mn></mrow></math></span>, which significantly degrade the approximation accuracy and mass conservation properties of standard PINNs. In this work, we propose a novel hybrid mass-constrained physics-informed neural network (Mc-PINN) framework for efficiently and accurately solving three types of CAC models in convex polygonal domains. The proposed method integrates deep learning with operator-splitting strategies to decompose the original CAC equations into simpler subproblems. One subproblem admits an analytical solution, while the other is solved using the Mc-PINN scheme. To further enhance efficiency, a Metropolis–Hastings based adaptive sampling strategy is introduced. In addition, we derive error estimates for the proposed method applied to all three CAC models. Numerical experiments demonstrate the robustness, accuracy, and mass-conservation capability of the proposed framework.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"184 ","pages":"Article 106643"},"PeriodicalIF":4.1,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145980172","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}
Pub Date : 2026-03-01Epub Date: 2026-01-07DOI: 10.1016/j.enganabound.2025.106629
Ahmed Mostafa Shaaban , Elena Atroshchenko , Steffen Marburg
A discontinuous isogeometric boundary element method is proposed for the acoustic model, which is based on the Helmholtz time-harmonic wave propagation equation, as part of acoustic–structural interaction problems. The discontinuous formulation is combined with an offset-based collocation scheme, which enhances accuracy over the continuous approach while simplifying integration, the computation of free terms, and the representation of highly distorted elements in pole-based models. The dynamic structural model, representing the second component of the interaction problem, is formulated using isogeometric Reissner–Mindlin shell theory, which is particularly effective for modeling thin curved structures. The acoustic and structural models are directly coupled through conforming numerical meshes on the interaction surface. Isogeometric analysis is applied for both the acoustic and structural formulations, as well as to the coupling scheme, due to its ability to represent exact geometries and its superior accuracy compared to conventional numerical approaches.
Numerical examples are presented to assess the performance of the proposed solution, with results compared against available analytical solutions and previously reported outcomes based on alternative coupling strategies.
{"title":"Discontinuous isogeometric boundary elements for direct acoustic–structural coupling with Reissner–Mindlin shells","authors":"Ahmed Mostafa Shaaban , Elena Atroshchenko , Steffen Marburg","doi":"10.1016/j.enganabound.2025.106629","DOIUrl":"10.1016/j.enganabound.2025.106629","url":null,"abstract":"<div><div>A discontinuous isogeometric boundary element method is proposed for the acoustic model, which is based on the Helmholtz time-harmonic wave propagation equation, as part of acoustic–structural interaction problems. The discontinuous formulation is combined with an offset-based collocation scheme, which enhances accuracy over the continuous approach while simplifying integration, the computation of free terms, and the representation of highly distorted elements in pole-based models. The dynamic structural model, representing the second component of the interaction problem, is formulated using isogeometric Reissner–Mindlin shell theory, which is particularly effective for modeling thin curved structures. The acoustic and structural models are directly coupled through conforming numerical meshes on the interaction surface. Isogeometric analysis is applied for both the acoustic and structural formulations, as well as to the coupling scheme, due to its ability to represent exact geometries and its superior accuracy compared to conventional numerical approaches.</div><div>Numerical examples are presented to assess the performance of the proposed solution, with results compared against available analytical solutions and previously reported outcomes based on alternative coupling strategies.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"184 ","pages":"Article 106629"},"PeriodicalIF":4.1,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145904247","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}
Pub Date : 2026-03-01Epub Date: 2026-01-09DOI: 10.1016/j.enganabound.2025.106627
Youjie Xu , Rui Yong , Lianjin Zhang , Fei Zhang , Zhenglin Mao , Xixiang Liu
The multiple multi-stage fractured horizontal wells (MFHWs) interference has been identified in tight gas reservoirs, but current method and model cannot achieve quantitative inter-well interference evaluation of non-uniform fractured region. Therefore, MFHWs mathematical model is established with consideration of multiple sub-region permeability difference. Coupling method of boundary element and source function is employed to solved the mathematical model. Inter-well interference factor is defined and used to evaluate inter-well interference degree. Total interference factor minimum value is used to determine the optimal rate ratio. The result shows that larger production time and well distance leads to small inter-well interference factor, but they have no influence optimal rate ratio. The optimal rate ratio will increase with the increasing of number of adjacent well fractures. If permeability of center region is larger than that of other region, inter-well interference factor and total interference factor will decrease and optimal rate ratio will increase. The model and method can evaluate inter-well interference degree and optimal rate ratio quantitatively, which provides guidance for horizontal well fracturing and parameter optimization design.
{"title":"An inter-well interference quantitative evaluation approach of multiple multi-stage fractured horizontal well with non-uniform simulated reservoirs volume in tight gas reservoirs","authors":"Youjie Xu , Rui Yong , Lianjin Zhang , Fei Zhang , Zhenglin Mao , Xixiang Liu","doi":"10.1016/j.enganabound.2025.106627","DOIUrl":"10.1016/j.enganabound.2025.106627","url":null,"abstract":"<div><div>The multiple multi-stage fractured horizontal wells (MFHWs) interference has been identified in tight gas reservoirs, but current method and model cannot achieve quantitative inter-well interference evaluation of non-uniform fractured region. Therefore, MFHWs mathematical model is established with consideration of multiple sub-region permeability difference. Coupling method of boundary element and source function is employed to solved the mathematical model. Inter-well interference factor is defined and used to evaluate inter-well interference degree. Total interference factor minimum value is used to determine the optimal rate ratio. The result shows that larger production time and well distance leads to small inter-well interference factor, but they have no influence optimal rate ratio. The optimal rate ratio will increase with the increasing of number of adjacent well fractures. If permeability of center region is larger than that of other region, inter-well interference factor and total interference factor will decrease and optimal rate ratio will increase. The model and method can evaluate inter-well interference degree and optimal rate ratio quantitatively, which provides guidance for horizontal well fracturing and parameter optimization design.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"184 ","pages":"Article 106627"},"PeriodicalIF":4.1,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145928754","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}
Pub Date : 2026-03-01Epub Date: 2026-01-09DOI: 10.1016/j.enganabound.2025.106621
Ziheng Li, Hong Zheng, Shouyang Huang
The original Discontinuous Deformation Analysis (DDA) faces limitations due to its linear displacement assumption and contact solution approach, making it challenging to handle nonlinear contact problems involving localized small deformations, particularly in cases like low-velocity impact on composite materials. This study presents an multi-point constraints enhanced DDA (MPC-DDA) specifically designed for modeling low-velocity impact scenarios in composites. The proposed modification incorporates two key features: (1) removal of the penalty spring mechanism preserves residual block penetration states in OC iterative results, (2) Iterative residuals are reconstructed as surface indentations through multi-point constraints, enabling full historical contact trace visualization. A novel contact algorithm utilizing virtual entrance points is developed, enabling MPC-DDA to dynamically capture the master and slave points on the contact boundary. This approach effectively represents surface indentations on composite materials through virtual entrance point displacements, significantly reducing the computational complexity associated with resolving localized small deformations characteristic of original DDA implementations. The MPC-DDA method has been implemented in a Matlab program, and its accuracy has been validated by comparing the computational results with those published in existing literature. Comparative analyses demonstrate that MPC-DDA outperforms original DDA in addressing nonlinear contact problems, particularly in terms of solution accuracy and computational efficiency.
{"title":"Multi-point constraints enhanced discontinuous deformation analysis for low-velocity impact on low composite materials","authors":"Ziheng Li, Hong Zheng, Shouyang Huang","doi":"10.1016/j.enganabound.2025.106621","DOIUrl":"10.1016/j.enganabound.2025.106621","url":null,"abstract":"<div><div>The original Discontinuous Deformation Analysis (DDA) faces limitations due to its linear displacement assumption and contact solution approach, making it challenging to handle nonlinear contact problems involving localized small deformations, particularly in cases like low-velocity impact on composite materials. This study presents an multi-point constraints enhanced DDA (MPC-DDA) specifically designed for modeling low-velocity impact scenarios in composites. The proposed modification incorporates two key features: (1) removal of the penalty spring mechanism preserves residual block penetration states in OC iterative results, (2) Iterative residuals are reconstructed as surface indentations through multi-point constraints, enabling full historical contact trace visualization. A novel contact algorithm utilizing virtual entrance points is developed, enabling MPC-DDA to dynamically capture the master and slave points on the contact boundary. This approach effectively represents surface indentations on composite materials through virtual entrance point displacements, significantly reducing the computational complexity associated with resolving localized small deformations characteristic of original DDA implementations. The MPC-DDA method has been implemented in a Matlab program, and its accuracy has been validated by comparing the computational results with those published in existing literature. Comparative analyses demonstrate that MPC-DDA outperforms original DDA in addressing nonlinear contact problems, particularly in terms of solution accuracy and computational efficiency.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"184 ","pages":"Article 106621"},"PeriodicalIF":4.1,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145928753","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}
Pub Date : 2026-03-01Epub Date: 2026-01-13DOI: 10.1016/j.enganabound.2026.106641
Yue Zhuo , Kai Chen , Degao Zou , Shanlin Tian , Shiyong Wu , Shi Zhang
Fracture analysis of anti-seepage systems in earth-rock dams involves pronounced material nonlinearity, frictional contact behavior, and convergence difficulties, which remain long-standing challenges in geotechnical engineering. In this study, A nonlinear 3D SBFEM–PFM formulation was created by embedding the phase-field model (PFM) into the scaled boundary finite element method (SBFEM) and implementing an intra-element phase-field interpolation scheme. Base on independently developed GEODYNA, a unified and parallelized framework was established through various coupling schemes, including FEM-SBFEM, phase-field, and nonlinear contact algorithm. The method is validated against a classical benchmark and subsequently applied to the world’s highest asphalt concrete core rockfill dam (ACCRD) on the overburden to simulate full-process cracking of gallery. Cracks were identified along the inner surface of the gallery and on the outer surfaces of both banks, with reservoir impoundment exhibiting opposing effects on crack width. Additionally, the severity of structural damage was significantly influenced by the interface characteristics between the gallery and surrounding geomaterial. This study signifies the inaugural implementation of PFM to achieve a full-process simulation of stress accumulation, crack initiation, and progressive propagation within dam anti-seepage systems. High risk zones were precisely identified, and the practical optimization measure was suggested, providing an innovative and effective approach for structural analysis in related geotechnical engineering applications.
{"title":"Crack analysis of the foundation gallery within an asphalt concrete core dam based on 3D SBFEM-PFM","authors":"Yue Zhuo , Kai Chen , Degao Zou , Shanlin Tian , Shiyong Wu , Shi Zhang","doi":"10.1016/j.enganabound.2026.106641","DOIUrl":"10.1016/j.enganabound.2026.106641","url":null,"abstract":"<div><div>Fracture analysis of anti-seepage systems in earth-rock dams involves pronounced material nonlinearity, frictional contact behavior, and convergence difficulties, which remain long-standing challenges in geotechnical engineering. In this study, A nonlinear 3D SBFEM–PFM formulation was created by embedding the phase-field model (PFM) into the scaled boundary finite element method (SBFEM) and implementing an intra-element phase-field interpolation scheme. Base on independently developed GEODYNA, a unified and parallelized framework was established through various coupling schemes, including FEM-SBFEM, phase-field, and nonlinear contact algorithm. The method is validated against a classical benchmark and subsequently applied to the world’s highest asphalt concrete core rockfill dam (ACCRD) on the overburden to simulate full-process cracking of gallery. Cracks were identified along the inner surface of the gallery and on the outer surfaces of both banks, with reservoir impoundment exhibiting opposing effects on crack width. Additionally, the severity of structural damage was significantly influenced by the interface characteristics between the gallery and surrounding geomaterial. This study signifies the inaugural implementation of PFM to achieve a full-process simulation of stress accumulation, crack initiation, and progressive propagation within dam anti-seepage systems. High risk zones were precisely identified, and the practical optimization measure was suggested, providing an innovative and effective approach for structural analysis in related geotechnical engineering applications.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"184 ","pages":"Article 106641"},"PeriodicalIF":4.1,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145956606","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}
Pub Date : 2026-03-01Epub Date: 2026-01-15DOI: 10.1016/j.enganabound.2026.106646
Gagan Sahoo, Harekrushna Behera, Tai-Wen Hsu
Capillary–gravity waves, influenced by both surface tension and gravity, interact strongly with marine structures, especially in the presence of uniform currents. Despite extensive studies on wave scattering by porous structures, the combined effects of surface tension, current, and porous barriers over a porous bottom remain insufficiently explored. This study examines the scattering of such waves by two thin surface-piercing porous barriers in the presence of a uniform current over a porous sea bed. A linear wave–structure interaction model is solved numerically through a hybrid Boundary Element-Finite Difference Method (BEM–FDM) and analytically through an eigenfunction expansion combined with a least-squares approach. The hybrid BEM–FDM efficiently handles higher-order boundary conditions that cannot be directly addressed by conventional BEM, while the analytical method eliminates the need for eigenfunction orthogonality and explicit mode coupling. The effects of surface tension, current velocity and direction, porous effect parameters of barriers as well as bottom, barrier length and spacing between them on reflection, transmission, and energy dissipation are analyzed. Results show that surface tension enhances reflection and dissipation while reducing transmission. Current direction strongly affects scattering: following currents enhance transmission, whereas opposing currents increase reflection and dissipation. Longer barriers and larger porous-effect parameters of both porous barriers and porous bottom enhance energy dissipation, while spacing between porous barriers induce interference driven oscillations.
{"title":"Scattering of capillary-gravity waves by surface-piercing porous barriers in the presence of uniform current over a porous sea bed","authors":"Gagan Sahoo, Harekrushna Behera, Tai-Wen Hsu","doi":"10.1016/j.enganabound.2026.106646","DOIUrl":"10.1016/j.enganabound.2026.106646","url":null,"abstract":"<div><div>Capillary–gravity waves, influenced by both surface tension and gravity, interact strongly with marine structures, especially in the presence of uniform currents. Despite extensive studies on wave scattering by porous structures, the combined effects of surface tension, current, and porous barriers over a porous bottom remain insufficiently explored. This study examines the scattering of such waves by two thin surface-piercing porous barriers in the presence of a uniform current over a porous sea bed. A linear wave–structure interaction model is solved numerically through a hybrid Boundary Element-Finite Difference Method (BEM–FDM) and analytically through an eigenfunction expansion combined with a least-squares approach. The hybrid BEM–FDM efficiently handles higher-order boundary conditions that cannot be directly addressed by conventional BEM, while the analytical method eliminates the need for eigenfunction orthogonality and explicit mode coupling. The effects of surface tension, current velocity and direction, porous effect parameters of barriers as well as bottom, barrier length and spacing between them on reflection, transmission, and energy dissipation are analyzed. Results show that surface tension enhances reflection and dissipation while reducing transmission. Current direction strongly affects scattering: following currents enhance transmission, whereas opposing currents increase reflection and dissipation. Longer barriers and larger porous-effect parameters of both porous barriers and porous bottom enhance energy dissipation, while spacing between porous barriers induce interference driven oscillations.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"184 ","pages":"Article 106646"},"PeriodicalIF":4.1,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145980173","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}
This paper presents a spectral method for the efficient evaluation of Green’s functions in three-dimensional anisotropic thermoelastic and thermomagnetoelectroelastic problems. The method expands the Green’s function kernel in spherical harmonics, reducing its integral representation to a finite sum containing only odd- or even-degree harmonic coefficients. This eliminates the need for mesh-based evaluation and interpolation, significantly improving computational speed and numerical robustness. The approach requires precomputation of only a small set of spectral coefficients, after which Green’s functions and their spatial derivatives can be evaluated rapidly at arbitrary target points. Moreover, the precomputed coefficients depend only on the material properties, can be calculated to the desired accuracy, and can be reused to analyze different solid geometries composed of the same material. The resulting formulation is well suited for boundary element methods (BEM) and other integral equation schemes. Numerical experiments demonstrate fast spectral convergence of the spherical harmonic series, while performance benchmarks show that the proposed approach reduces the overall BEM computation time for the considered problems by approximately a factor of two. The method provides an efficient and scalable tool for simulating complex multiphysics phenomena in anisotropic solids.
{"title":"A spectral approach for fast evaluation of 3D thermomagnetoelectroelastic Green’s functions and their derivatives in the boundary element method","authors":"Viktoriya Pasternak , Heorhiy Sulym , Andrii Korniichuk , Iaroslav Pasternak","doi":"10.1016/j.enganabound.2026.106647","DOIUrl":"10.1016/j.enganabound.2026.106647","url":null,"abstract":"<div><div>This paper presents a spectral method for the efficient evaluation of Green’s functions in three-dimensional anisotropic thermoelastic and thermomagnetoelectroelastic problems. The method expands the Green’s function kernel in spherical harmonics, reducing its integral representation to a finite sum containing only odd- or even-degree harmonic coefficients. This eliminates the need for mesh-based evaluation and interpolation, significantly improving computational speed and numerical robustness. The approach requires precomputation of only a small set of spectral coefficients, after which Green’s functions and their spatial derivatives can be evaluated rapidly at arbitrary target points. Moreover, the precomputed coefficients depend only on the material properties, can be calculated to the desired accuracy, and can be reused to analyze different solid geometries composed of the same material. The resulting formulation is well suited for boundary element methods (BEM) and other integral equation schemes. Numerical experiments demonstrate fast spectral convergence of the spherical harmonic series, while performance benchmarks show that the proposed approach reduces the overall BEM computation time for the considered problems by approximately a factor of two. The method provides an efficient and scalable tool for simulating complex multiphysics phenomena in anisotropic solids.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"184 ","pages":"Article 106647"},"PeriodicalIF":4.1,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145962064","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}
扫码关注我们
求助内容：
应助结果提醒方式：
京公网安备 11010802042870号