Pub Date : 2026-01-09DOI: 10.1016/j.enganabound.2025.106625
Chenliang Li, Donglin Guo, Huihua Zhang
Magneto-electro-elastic (MEE) materials are pivotal in a multitude of fields. In this work, the numerical manifold method (NMM) is innovatively extended to establish 2-D numerical models for perforated MEE solids. The unique dual-cover system, namely, the mathematical cover and the physical cover, enables the NMM to discretize the physical domain with non-conforming mathematical covers straightforwardly. The governing equations and the boundary conditions for hole problems of MEE materials are firstly introduced. Then, by taking into account the governing equations, boundary conditions, and the NMM field approximations, the NMM global discrete equations are derived using the weighted residual method. Through three benchmark examples, the precision of the proposed method is verified, and it is then applied to study two more complex cases, where the effects of hole configurations, loading conditions, and polarization directions on the field quantities of perforated MEE materials are further examined.
{"title":"Modeling of magneto-electro-elastic solids with complex cutouts by the numerical manifold method","authors":"Chenliang Li, Donglin Guo, Huihua Zhang","doi":"10.1016/j.enganabound.2025.106625","DOIUrl":"10.1016/j.enganabound.2025.106625","url":null,"abstract":"<div><div>Magneto-electro-elastic (MEE) materials are pivotal in a multitude of fields. In this work, the numerical manifold method (NMM) is innovatively extended to establish 2-D numerical models for perforated MEE solids. The unique dual-cover system, namely, the mathematical cover and the physical cover, enables the NMM to discretize the physical domain with non-conforming mathematical covers straightforwardly. The governing equations and the boundary conditions for hole problems of MEE materials are firstly introduced. Then, by taking into account the governing equations, boundary conditions, and the NMM field approximations, the NMM global discrete equations are derived using the weighted residual method. Through three benchmark examples, the precision of the proposed method is verified, and it is then applied to study two more complex cases, where the effects of hole configurations, loading conditions, and polarization directions on the field quantities of perforated MEE materials are further examined.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"184 ","pages":"Article 106625"},"PeriodicalIF":4.1,"publicationDate":"2026-01-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145928752","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-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-01-09","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-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-01-09","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-01-07DOI: 10.1016/j.enganabound.2025.106628
Yonghui Li , Mengyao Ma , Dinghao Zhang , Jiawei He
Soil seepage erosion caused by leaks in underground drainage pipelines is one of the key factors leading to ground subsidence. This study establishes a numerical model coupling computational fluid dynamics with the discrete element method (CFD-DEM) to conduct a microscopic investigation of seepage erosion phenomena in unsaturated silt. This numerical model employs the Richards equation to describe seepage processes in unsaturated silt and combines it with the VG model to determine parameters such as the soil's water capacity, saturation, and hydraulic conductivity. Furthermore, based on existing formulas for calculating the shear strength of unsaturated silt, a dynamic updating algorithm for particle cohesion strength was developed to account for saturation effects, thereby reproducing the moisture-induced degradation characteristics of silt. By comparing with laboratory tests, the numerical model was validated for its accuracy in simulating pore water pressure evolution, wetting front propagation rates, bonding strength updates, and soil particle loss processes. Finally, based on numerical simulation results, this study analyzed changes in formation water pressure, contact force chains, and pipeline stress conditions during the erosion-loss process. It revealed the intrinsic mechanism by which leakage-induced erosion of underground drainage pipelines triggers ground collapse at the microscopic level.
{"title":"A CFD-DEM method for simulating the erosion evolution of unsaturated silt induced by leakage through underground pipeline defect","authors":"Yonghui Li , Mengyao Ma , Dinghao Zhang , Jiawei He","doi":"10.1016/j.enganabound.2025.106628","DOIUrl":"10.1016/j.enganabound.2025.106628","url":null,"abstract":"<div><div>Soil seepage erosion caused by leaks in underground drainage pipelines is one of the key factors leading to ground subsidence. This study establishes a numerical model coupling computational fluid dynamics with the discrete element method (CFD-DEM) to conduct a microscopic investigation of seepage erosion phenomena in unsaturated silt. This numerical model employs the Richards equation to describe seepage processes in unsaturated silt and combines it with the VG model to determine parameters such as the soil's water capacity, saturation, and hydraulic conductivity. Furthermore, based on existing formulas for calculating the shear strength of unsaturated silt, a dynamic updating algorithm for particle cohesion strength was developed to account for saturation effects, thereby reproducing the moisture-induced degradation characteristics of silt. By comparing with laboratory tests, the numerical model was validated for its accuracy in simulating pore water pressure evolution, wetting front propagation rates, bonding strength updates, and soil particle loss processes. Finally, based on numerical simulation results, this study analyzed changes in formation water pressure, contact force chains, and pipeline stress conditions during the erosion-loss process. It revealed the intrinsic mechanism by which leakage-induced erosion of underground drainage pipelines triggers ground collapse at the microscopic level.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"184 ","pages":"Article 106628"},"PeriodicalIF":4.1,"publicationDate":"2026-01-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145928767","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-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-01-07","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-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-01-06","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-01-03DOI: 10.1016/j.enganabound.2025.106626
Yifei Zhang
In this work, the linear-gradient smoothing meshfree Galerkin method (LGSM) is developed to analyze 2D static fracture problems of multiple materials, including isotropic materials, orthotropic materials, and functionally graded materials (FGMs). The discontinuous displacement fields and singular stress fields can be characterized by combining the diffraction method with intrinsic enrichment in the meshfree Galerkin method. However, the meshfree approximation based on intrinsic enrichment complicates numerical integration due to the non-polynomial nature of the shape function. Typically, high-order Gaussian integration is required for domain integration in the meshfree Galerkin method, but this is time-consuming. To address this, we introduce the linear-gradient smoothed integral scheme (LGSI) into the Galerkin weak form to construct the system matrix. LGSI uses only two Gaussian integration points at each boundary of the smoothing domain for the numerical integration. This process can significantly reduce the number of integral points and improves computational efficiency. Additionally, LGSI generates linear smoothed strains using the smoothed gradient technique, thereby providing accurate numerical results. The static stress intensity factors (SIFs) are determined using interaction integration. Through several numerical crack problems, the accuracy and efficiency of the LGSM are validated by comparing the results with those obtained from different methods in existing literature.
{"title":"A linear-gradient smoothing meshfree Galerkin method for 2D static fracture analysis of multiple materials","authors":"Yifei Zhang","doi":"10.1016/j.enganabound.2025.106626","DOIUrl":"10.1016/j.enganabound.2025.106626","url":null,"abstract":"<div><div>In this work, the linear-gradient smoothing meshfree Galerkin method (LGSM) is developed to analyze 2D static fracture problems of multiple materials, including isotropic materials, orthotropic materials, and functionally graded materials (FGMs). The discontinuous displacement fields and singular stress fields can be characterized by combining the diffraction method with intrinsic enrichment in the meshfree Galerkin method. However, the meshfree approximation based on intrinsic enrichment complicates numerical integration due to the non-polynomial nature of the shape function. Typically, high-order Gaussian integration is required for domain integration in the meshfree Galerkin method, but this is time-consuming. To address this, we introduce the linear-gradient smoothed integral scheme (LGSI) into the Galerkin weak form to construct the system matrix. LGSI uses only two Gaussian integration points at each boundary of the smoothing domain for the numerical integration. This process can significantly reduce the number of integral points and improves computational efficiency. Additionally, LGSI generates linear smoothed strains using the smoothed gradient technique, thereby providing accurate numerical results. The static stress intensity factors (SIFs) are determined using interaction integration. Through several numerical crack problems, the accuracy and efficiency of the LGSM are validated by comparing the results with those obtained from different methods in existing literature.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"183 ","pages":"Article 106626"},"PeriodicalIF":4.1,"publicationDate":"2026-01-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145884140","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-01-02DOI: 10.1016/j.enganabound.2025.106613
Wenjie Luo , Chaorong Li , Xudong Ling , Yilan Xue , L.O. Mubashiru
The CNN-based conditional encoder extracts prior knowledge for medical image segmentation to support the diffusion model. However, convolution operations are limited in capturing cross-channel and spatial dependencies, leading to the loss of crucial conditional information. Moreover, during the denoising UNet’s decoder, insufficient fusion of cross-layer features during up-sampling yields blurry feature maps and weak adaptability of subsequent convolutions. To tackle these challenges, we propose the Spatial and Channel Mixing Attention (SCMA) for the encoder and the Multi-Scale Feature Modulation Residual Module (MFMRM) for decoder feature fusion. MFMRM comprises the Multi-Scale Context Modulation Module (MCMM) and a Residual Dual Convolution Module (RDCM), adaptively integrating multi-resolution features to enhance representation and improve robustness to input variations. Furthermore, by introducing the Kolmogorov–Arnold Network (KAN) to optimize SCMA, we obtain KSCMA, which mitigates the curse of dimensionality and strengthens the representation of critical features. Experiments on ultrasound thyroid nodule, MRI brain tumor, and a self-constructed invasive breast cancer dataset demonstrate that our approach outperforms existing methods in segmentation accuracy. Our project is open source and available on GitHub at: https://github.com/lwj018/MSMedDiff-1.
{"title":"Enhancing cancer segmentation using conditional diffusion networks with KAN optimized attention and adaptive feature fusion","authors":"Wenjie Luo , Chaorong Li , Xudong Ling , Yilan Xue , L.O. Mubashiru","doi":"10.1016/j.enganabound.2025.106613","DOIUrl":"10.1016/j.enganabound.2025.106613","url":null,"abstract":"<div><div>The CNN-based conditional encoder extracts prior knowledge for medical image segmentation to support the diffusion model. However, convolution operations are limited in capturing cross-channel and spatial dependencies, leading to the loss of crucial conditional information. Moreover, during the denoising UNet’s decoder, insufficient fusion of cross-layer features during up-sampling yields blurry feature maps and weak adaptability of subsequent convolutions. To tackle these challenges, we propose the Spatial and Channel Mixing Attention (SCMA) for the encoder and the Multi-Scale Feature Modulation Residual Module (MFMRM) for decoder feature fusion. MFMRM comprises the Multi-Scale Context Modulation Module (MCMM) and a Residual Dual Convolution Module (RDCM), adaptively integrating multi-resolution features to enhance representation and improve robustness to input variations. Furthermore, by introducing the Kolmogorov–Arnold Network (KAN) to optimize SCMA, we obtain KSCMA, which mitigates the curse of dimensionality and strengthens the representation of critical features. Experiments on ultrasound thyroid nodule, MRI brain tumor, and a self-constructed invasive breast cancer dataset demonstrate that our approach outperforms existing methods in segmentation accuracy. Our project is open source and available on GitHub at: <span><span>https://github.com/lwj018/MSMedDiff-1</span><svg><path></path></svg></span>.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"183 ","pages":"Article 106613"},"PeriodicalIF":4.1,"publicationDate":"2026-01-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145884142","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 : 2025-12-31DOI: 10.1016/j.enganabound.2025.106618
Zhuoheng Wang , Wenxuan Xie , Junseok Kim , Yibao Li
Reliability-based topology optimization is the approach that incorporates uncertainty quantification into topology optimization to achieve robustness and reliability. In this paper, we propose a reliability-based topology optimization algorithm that integrates a phase-field model with uncertainty quantification techniques. Material uncertainty in Young’s modulus is modeled as a random field and reduced via the Karhunen–Loève expansion. A sequential optimization and reliability assessment strategy is adopted to decouple the optimization and reliability evaluation for reducing computational complexity. The topology optimization problem is solved by using the coupled finite element and finite difference approach, in which time integration is performed using a second-order Crank–Nicolson scheme. Reliability analysis is conducted by using the inverse reliability analysis method, with the limit-state function approximated by a stochastic response surface method. Monte Carlo simulations are performed to validate the accuracy of the computed failure probabilities. Numerical results confirm the efficiency and robustness of the proposed method in generating reliable structures under material uncertainty.
{"title":"Efficient phase field structural design algorithm for reliability-based topology optimization with material uncertainties","authors":"Zhuoheng Wang , Wenxuan Xie , Junseok Kim , Yibao Li","doi":"10.1016/j.enganabound.2025.106618","DOIUrl":"10.1016/j.enganabound.2025.106618","url":null,"abstract":"<div><div>Reliability-based topology optimization is the approach that incorporates uncertainty quantification into topology optimization to achieve robustness and reliability. In this paper, we propose a reliability-based topology optimization algorithm that integrates a phase-field model with uncertainty quantification techniques. Material uncertainty in Young’s modulus is modeled as a random field and reduced via the Karhunen–Loève expansion. A sequential optimization and reliability assessment strategy is adopted to decouple the optimization and reliability evaluation for reducing computational complexity. The topology optimization problem is solved by using the coupled finite element and finite difference approach, in which time integration is performed using a second-order Crank–Nicolson scheme. Reliability analysis is conducted by using the inverse reliability analysis method, with the limit-state function approximated by a stochastic response surface method. Monte Carlo simulations are performed to validate the accuracy of the computed failure probabilities. Numerical results confirm the efficiency and robustness of the proposed method in generating reliable structures under material uncertainty.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"183 ","pages":"Article 106618"},"PeriodicalIF":4.1,"publicationDate":"2025-12-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145884143","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 : 2025-12-31DOI: 10.1016/j.enganabound.2025.106615
Zhongting Xu , Zhengjie Sun , Shengliang Zhang
This paper proposes a novel meshless Galerkin energy-preserving method for conservative evolutionary partial differential equations (PDEs). By combining the meshless Galerkin approach for spatial discretization with the average vector field method, we establish a general framework for constructing meshless energy-preserving schemes applicable to a wide range of PDEs. Moreover, we present explicit schemes for the Korteweg–de Vries and nonlinear Schrödinger equations, analyzing their energy conservation properties and convergence results. A salient feature of our method is the use of problem-driven kernels, with periodic kernels for one-dimensional periodic problems and localized kernels for two-dimensional ones. Numerical results demonstrate the high accuracy of our method on scattered nodes and highlight the benefits of the meshless approach. It also achieves superior energy and mass preservation compared to mesh-dependent methods.
{"title":"Meshless energy-conserving schemes for conservative partial differential equations using kernel-based Galerkin methods","authors":"Zhongting Xu , Zhengjie Sun , Shengliang Zhang","doi":"10.1016/j.enganabound.2025.106615","DOIUrl":"10.1016/j.enganabound.2025.106615","url":null,"abstract":"<div><div>This paper proposes a novel meshless Galerkin energy-preserving method for conservative evolutionary partial differential equations (PDEs). By combining the meshless Galerkin approach for spatial discretization with the average vector field method, we establish a general framework for constructing meshless energy-preserving schemes applicable to a wide range of PDEs. Moreover, we present explicit schemes for the Korteweg–de Vries and nonlinear Schrödinger equations, analyzing their energy conservation properties and convergence results. A salient feature of our method is the use of problem-driven kernels, with periodic kernels for one-dimensional periodic problems and localized kernels for two-dimensional ones. Numerical results demonstrate the high accuracy of our method on scattered nodes and highlight the benefits of the meshless approach. It also achieves superior energy and mass preservation compared to mesh-dependent methods.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"183 ","pages":"Article 106615"},"PeriodicalIF":4.1,"publicationDate":"2025-12-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145884139","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号