Pub Date : 2026-02-01Epub Date: 2025-12-31DOI: 10.1016/j.enganabound.2025.106624
Yu Tian , Fu-Ren Ming , Hao Chen , Xiang-Li Fang , Ping-Ping Wang
The metal jet formed by a conventional shaped charge has a high penetration depth but a narrow damage range. To expand the damage range of the metal jet, an annular shaped charge can be employed. In this paper, the underwater explosions of annular shaped charges are simulated by a graphic processing unit accelerated axisymmetric Riemann-SPH method, and the accuracy is verified by experiments. The hole-cutting effect using an annular shaped charge in an underwater explosion is analyzed, and the characteristics of the annular shaped charge are compared with a spherical segment shaped charge. Furthermore, the effects of the liner thickness on the annular jet are explored. It is revealed that the damage mode inflicted on the plate by an explosive formed projectile (EFP) is impact penetration, whereas that by an annular jet is cutting. The plate’s breach size caused by the annular jet reaches 1.79 times the charge radius and 3.15 times that of EFP. Our findings also reveal that there is an optimal dimensionless maximum liner thickness (/) to maximize the breach size , which is 0.12 for a charge mass of 60kg. This paper can provide support for the optimization design of the shaped charge.
{"title":"Numerical study of hole-cutting effect using annular shaped charges in underwater explosions","authors":"Yu Tian , Fu-Ren Ming , Hao Chen , Xiang-Li Fang , Ping-Ping Wang","doi":"10.1016/j.enganabound.2025.106624","DOIUrl":"10.1016/j.enganabound.2025.106624","url":null,"abstract":"<div><div>The metal jet formed by a conventional shaped charge has a high penetration depth but a narrow damage range. To expand the damage range of the metal jet, an annular shaped charge can be employed. In this paper, the underwater explosions of annular shaped charges are simulated by a graphic processing unit accelerated axisymmetric Riemann-SPH method, and the accuracy is verified by experiments. The hole-cutting effect using an annular shaped charge in an underwater explosion is analyzed, and the characteristics of the annular shaped charge are compared with a spherical segment shaped charge. Furthermore, the effects of the liner thickness on the annular jet are explored. It is revealed that the damage mode inflicted on the plate by an explosive formed projectile (EFP) is impact penetration, whereas that by an annular jet is cutting. The plate’s breach size caused by the annular jet reaches 1.79 times the charge radius and 3.15 times that of EFP. Our findings also reveal that there is an optimal dimensionless maximum liner thickness <span><math><mi>λ</mi></math></span> (<span><math><mrow><mi>λ</mi><mo>=</mo><mi>T</mi></mrow></math></span>/<span><math><msub><mrow><mi>R</mi></mrow><mrow><mi>z</mi></mrow></msub></math></span>) to maximize the breach size <span><math><msub><mrow><mi>d</mi></mrow><mrow><mi>p</mi></mrow></msub></math></span>, which is 0.12 for a charge mass of 60kg. This paper can provide support for the optimization design of the shaped charge.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"183 ","pages":"Article 106624"},"PeriodicalIF":4.1,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145884238","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-02-01Epub Date: 2025-12-19DOI: 10.1016/j.enganabound.2025.106603
Pelin Senel , Daniel Lesnic , Andreas Karageorghis
We consider a geometric inverse problem which requires detecting an unknown obstacle, e.g. a submarine or an aquatic mine, submerged in a Stokes slow viscous stationary flow of an incompressible fluid. In two-dimensions, the problem is formulated in terms of the biharmonic streamfunction in an unbounded domain which is approximated using the Trefftz method and the method of fundamental solutions (MFS). This is, apparently, the first time the Trefftz method and the MFS are applied for the solution of the biharmonic equation in an unbounded domain. We first examine direct problems and then consider inverse problems. The unknown obstacle is determined by employing a nonlinear Tikhonov regularization procedure. Numerical results are presented and discussed.
{"title":"Meshless methods for the detection of an obstacle submerged in a two-dimensional Stokes flow","authors":"Pelin Senel , Daniel Lesnic , Andreas Karageorghis","doi":"10.1016/j.enganabound.2025.106603","DOIUrl":"10.1016/j.enganabound.2025.106603","url":null,"abstract":"<div><div>We consider a geometric inverse problem which requires detecting an unknown obstacle, e.g. a submarine or an aquatic mine, submerged in a Stokes slow viscous stationary flow of an incompressible fluid. In two-dimensions, the problem is formulated in terms of the biharmonic streamfunction in an unbounded domain which is approximated using the Trefftz method and the method of fundamental solutions (MFS). This is, apparently, the first time the Trefftz method and the MFS are applied for the solution of the biharmonic equation in an unbounded domain. We first examine direct problems and then consider inverse problems. The unknown obstacle is determined by employing a nonlinear Tikhonov regularization procedure. Numerical results are presented and discussed.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"183 ","pages":"Article 106603"},"PeriodicalIF":4.1,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145784768","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-02-01Epub Date: 2025-12-16DOI: 10.1016/j.enganabound.2025.106601
Chunquan Li, Runxin Cao, Yuanhao Zheng, Hongyan Huang, Ming Zhang
The microchannel heat sink with gradient array of ribs and pin fins (MCHS-GDRPF) integrates two types of thermal microstructures, rib and pin fin, which effectively enhances the thermal performance of the microchannel and improves the temperature uniformity. To enhance the thermal performance of the MCHS-GDRPF, machine learning methods are employed to predict and optimize performance based on structural design parameters. Four machine learning methods were evaluated to construct the thermal performance prediction model. The Support Vector Regression (SVR) model, demonstrating the highest accuracy, was then integrated with the fast Non-dominated Sorting Genetic Algorithm-II (NSGA-II) to perform multi-objective optimization of thermal resistance, temperature inhomogeneity, and friction coefficient. Subsequently, the Pareto front solution set was derived, and the optimal structural parameters were determined via the TOPSIS method. Finally, a thermal performance analysis was conducted for MCHS-GDRPF with four distinct structural parameter configurations. Comparative evaluation at Reynolds number 622 indicates three key advantages of the MCHS-GDRPF: (1) 28 K lower maximum surface temperature, (2) 72% improvement in temperature uniformity, and (3) 1.397 PEC value, surpassing conventional microchannel performance.
{"title":"A machine learning-driven multi-objective optimization framework for advanced microchannel heat sinks with gradient rib-pin fin arrays","authors":"Chunquan Li, Runxin Cao, Yuanhao Zheng, Hongyan Huang, Ming Zhang","doi":"10.1016/j.enganabound.2025.106601","DOIUrl":"10.1016/j.enganabound.2025.106601","url":null,"abstract":"<div><div>The microchannel heat sink with gradient array of ribs and pin fins (MCHS-GDRPF) integrates two types of thermal microstructures, rib and pin fin, which effectively enhances the thermal performance of the microchannel and improves the temperature uniformity. To enhance the thermal performance of the MCHS-GDRPF, machine learning methods are employed to predict and optimize performance based on structural design parameters. Four machine learning methods were evaluated to construct the thermal performance prediction model. The Support Vector Regression (SVR) model, demonstrating the highest accuracy, was then integrated with the fast Non-dominated Sorting Genetic Algorithm-II (NSGA-II) to perform multi-objective optimization of thermal resistance, temperature inhomogeneity, and friction coefficient. Subsequently, the Pareto front solution set was derived, and the optimal structural parameters were determined via the TOPSIS method. Finally, a thermal performance analysis was conducted for MCHS-GDRPF with four distinct structural parameter configurations. Comparative evaluation at Reynolds number 622 indicates three key advantages of the MCHS-GDRPF: (1) 28 K lower maximum surface temperature, (2) 72% improvement in temperature uniformity, and (3) 1.397 PEC value, surpassing conventional microchannel performance.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"183 ","pages":"Article 106601"},"PeriodicalIF":4.1,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145784769","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-02-01Epub Date: 2025-12-13DOI: 10.1016/j.enganabound.2025.106600
Haoran Yan , Yunpeng Lu , Hong Song , Yan Wang , Guiyong Zhang
This study develops a multiple-relaxation-time multiphase lattice Boltzmann flux solver (MRT-MLBFS) for the simulation of incompressible multiphase flows involving Newtonian and non-Newtonian fluids with large density ratio and rheological contrasts. The method integrates a phase-field model for interface capturing and a finite-volume-based LBFS framework to solve the macroscopic Navier–Stokes equations. By combining a MRT formulation with a flux reconstruction strategy, the method effectively suppresses numerical instabilities while reducing artificial dissipation, thereby ensuring robust simulations under extreme density and viscosity contrasts. Additionally, the algorithm is implemented in a fully explicit scheme and accelerated on GPU, which leads to a significant gain in computational efficiency, achieving speedups exceeding two orders of magnitude compared to CPU implementations. Through a series of benchmark cases, including droplet on plate, Rayleigh–Taylor instability, and droplet spreading on thin film, MRT-MLBFS is shown to capture complex interfacial evolution and nonlinear rheological effects with high accuracy and stability. This work establishes MRT-MLBFS as a reliable solution strategy for non-Newtonian multiphase flows with broad applicability.
{"title":"Multiphase lattice Boltzmann flux solver for non-Newtonian power-law fluid flows with high efficiency and stability","authors":"Haoran Yan , Yunpeng Lu , Hong Song , Yan Wang , Guiyong Zhang","doi":"10.1016/j.enganabound.2025.106600","DOIUrl":"10.1016/j.enganabound.2025.106600","url":null,"abstract":"<div><div>This study develops a multiple-relaxation-time multiphase lattice Boltzmann flux solver (MRT-MLBFS) for the simulation of incompressible multiphase flows involving Newtonian and non-Newtonian fluids with large density ratio and rheological contrasts. The method integrates a phase-field model for interface capturing and a finite-volume-based LBFS framework to solve the macroscopic Navier–Stokes equations. By combining a MRT formulation with a flux reconstruction strategy, the method effectively suppresses numerical instabilities while reducing artificial dissipation, thereby ensuring robust simulations under extreme density and viscosity contrasts. Additionally, the algorithm is implemented in a fully explicit scheme and accelerated on GPU, which leads to a significant gain in computational efficiency, achieving speedups exceeding two orders of magnitude compared to CPU implementations. Through a series of benchmark cases, including droplet on plate, Rayleigh–Taylor instability, and droplet spreading on thin film, MRT-MLBFS is shown to capture complex interfacial evolution and nonlinear rheological effects with high accuracy and stability. This work establishes MRT-MLBFS as a reliable solution strategy for non-Newtonian multiphase flows with broad applicability.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"183 ","pages":"Article 106600"},"PeriodicalIF":4.1,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145730708","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-02-01Epub Date: 2025-12-10DOI: 10.1016/j.enganabound.2025.106595
Y.C. Shiah
A common concern in the use of thin laminated composites is the potential debonding caused by interlaminar stresses of adhesives. However, their simulations by conventional schemes present modelling difficulties stemming from the thinness of each laminate and adhesive layer. For that, conventional simulations usually neglect the presence of thin adhesives on interfaces. Due to such negligence, any potential debonding from the failure of adhesives cannot be properly assessed. This paper presents efficiently modelling of the interlaminar stresses in 3D thin laminated composites when subjected to the inertial loads of self-weight or rotations using the boundary element method (BEM). By applying the FG-Squircular mapping (introduced by Manuel Fernandez-Guasti in 1992), the transformed boundary singular integrals due to the inertial loads are completely regularized. By the BEM modelling, thin laminated composites bonded with thin sheets of adhesives on interfaces can thus be modelled using very coarse meshes.
{"title":"BEM simulation of the interlaminar stresses of 3D thin composites subjected to inertial loads","authors":"Y.C. Shiah","doi":"10.1016/j.enganabound.2025.106595","DOIUrl":"10.1016/j.enganabound.2025.106595","url":null,"abstract":"<div><div>A common concern in the use of thin laminated composites is the potential debonding caused by interlaminar stresses of adhesives. However, their simulations by conventional schemes present modelling difficulties stemming from the thinness of each laminate and adhesive layer. For that, conventional simulations usually neglect the presence of thin adhesives on interfaces. Due to such negligence, any potential debonding from the failure of adhesives cannot be properly assessed. This paper presents efficiently modelling of the interlaminar stresses in 3D thin laminated composites when subjected to the inertial loads of self-weight or rotations using the boundary element method (BEM). By applying the FG-Squircular mapping (introduced by Manuel Fernandez-Guasti in 1992), the transformed boundary singular integrals due to the inertial loads are completely regularized. By the BEM modelling, thin laminated composites bonded with thin sheets of adhesives on interfaces can thus be modelled using very coarse meshes.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"183 ","pages":"Article 106595"},"PeriodicalIF":4.1,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145730700","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-02-01Epub Date: 2025-12-03DOI: 10.1016/j.enganabound.2025.106553
Ruiping Niu , Shijie Zhao , Xinglong Lu , Qiuxia Fan , Lei han Wang
This paper proposes, for the first time, an adaptive cell-based smoothed finite element method (A-CS-FEM) based on arbitrary polygonal elements for thermo-mechanical coupling problems. By utilizing mean value coordinates, the proposed model accommodates non-convex polygons without self-intersection, enabling robust handling of diverse polygonal elements. The approach integrates constrained Delaunay triangulation with adaptive techniques to partition triangulation-cell-based smoothing domains, naturally ensuring the positivity condition for a normed space without additional stabilization, while maintaining the high-quality meshes. The gradient smoothing technique in CS-FEM eliminates the coordinate mapping inherent in traditional polygonal finite element methods, because it requires only the shape function values along the segments of cell smoothing domains instead of the shape function derivatives. Numerical results demonstrate that A-CS-FEM significantly improves the quality of smoothing domains for complex geometries with arbitrary convex and concave polygonal discretization, thereby achieving high-precision solutions for both displacement and temperature.
{"title":"An adaptive cell-based smoothed finite element method with arbitrary polygonal elements for coupled thermo-mechanical analysis","authors":"Ruiping Niu , Shijie Zhao , Xinglong Lu , Qiuxia Fan , Lei han Wang","doi":"10.1016/j.enganabound.2025.106553","DOIUrl":"10.1016/j.enganabound.2025.106553","url":null,"abstract":"<div><div>This paper proposes, for the first time, an adaptive cell-based smoothed finite element method (A-CS-FEM) based on arbitrary polygonal elements for thermo-mechanical coupling problems. By utilizing mean value coordinates, the proposed model accommodates non-convex polygons without self-intersection, enabling robust handling of diverse polygonal elements. The approach integrates constrained Delaunay triangulation with adaptive techniques to partition triangulation-cell-based smoothing domains, naturally ensuring the positivity condition for a normed <span><math><msubsup><mi>G</mi><mi>h</mi><mn>1</mn></msubsup></math></span> space without additional stabilization, while maintaining the high-quality meshes. The gradient smoothing technique in CS-FEM eliminates the coordinate mapping inherent in traditional polygonal finite element methods, because it requires only the shape function values along the segments of cell smoothing domains instead of the shape function derivatives. Numerical results demonstrate that A-CS-FEM significantly improves the quality of smoothing domains for complex geometries with arbitrary convex and concave polygonal discretization, thereby achieving high-precision solutions for both displacement and temperature.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"183 ","pages":"Article 106553"},"PeriodicalIF":4.1,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145685098","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-02-01Epub Date: 2025-12-02DOI: 10.1016/j.enganabound.2025.106569
Zefeng Liu , Jinshuai Bai , Yuantong Gu , Ping Xiang
This paper introduces a novel computational framework that integrates physics-informed neural network (PINN) with generalized layerwise theory (LW) for the bending analysis of laminated composite plates. The framework leverages the approximation capability of deep neural networks while incorporating the physical constraints from LW theory to accurately capture the displacement fields of laminated composite plates as well as the shear stresses variations along the thickness direction. The framework is validated using various laminated plate configurations and loading conditions, with results showing excellent agreement with the meshless radial point interpolation method (RPIM), as well as other published solutions. These results highlight the potential of the PINN framework to enhance the predictive bending analysis of laminated composite plates, positioning it as a promising alternative for laminated composite structures.
{"title":"Physics-informed neural network based on layerwise theory for bending analysis of laminated plates","authors":"Zefeng Liu , Jinshuai Bai , Yuantong Gu , Ping Xiang","doi":"10.1016/j.enganabound.2025.106569","DOIUrl":"10.1016/j.enganabound.2025.106569","url":null,"abstract":"<div><div>This paper introduces a novel computational framework that integrates physics-informed neural network (PINN) with generalized layerwise theory (LW) for the bending analysis of laminated composite plates. The framework leverages the approximation capability of deep neural networks while incorporating the physical constraints from LW theory to accurately capture the displacement fields of laminated composite plates as well as the shear stresses variations along the thickness direction. The framework is validated using various laminated plate configurations and loading conditions, with results showing excellent agreement with the meshless radial point interpolation method (RPIM), as well as other published solutions. These results highlight the potential of the PINN framework to enhance the predictive bending analysis of laminated composite plates, positioning it as a promising alternative for laminated composite structures.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"183 ","pages":"Article 106569"},"PeriodicalIF":4.1,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145645759","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-02-01Epub Date: 2025-12-02DOI: 10.1016/j.enganabound.2025.106561
Dianjian Ruan , Zhanheng Chen , Daming Yuan
Bilateral obstacle problems are fundamental in the study of partial differential equations (PDEs) and variational inequalities, with significant applications in optimal control, elasticity, and material deformation under constraints. However, numerically solving these problems is challenging due to the inherent nonlinearities and the presence of free boundaries that evolve with complex contact dynamics. Conventional discretization methods, including finite element and finite difference approaches, often struggle to balance accuracy with computational efficiency, especially when dealing with irregular geometries or the need for adaptive resolution. In the present work we introduce a meshless method that overcomes these challenges by combining the modified finite-particle method (MFPM) for discretization with the Picard iteration technique for solving the result piecewise linear system. The proposed technique employs adaptive stencil selection to guarantee a result linear system with a moderate condition number. An adaptive meshless refinement method enhances the free boundary resolution, particularly in capturing the unknown free boundary a priori. Numerical experiments confirm the method’s flexibility and robustness across a range of node layouts – including Cartesian grids, PNP nodes, and Halton points – demonstrating its potential as an effective tool for solving bilateral obstacle problems and broadening the applicability of PDE and variational inequality models.
{"title":"A modified finite particle method with adaptive strategy for solving bilateral obstacle problems","authors":"Dianjian Ruan , Zhanheng Chen , Daming Yuan","doi":"10.1016/j.enganabound.2025.106561","DOIUrl":"10.1016/j.enganabound.2025.106561","url":null,"abstract":"<div><div>Bilateral obstacle problems are fundamental in the study of partial differential equations (PDEs) and variational inequalities, with significant applications in optimal control, elasticity, and material deformation under constraints. However, numerically solving these problems is challenging due to the inherent nonlinearities and the presence of free boundaries that evolve with complex contact dynamics. Conventional discretization methods, including finite element and finite difference approaches, often struggle to balance accuracy with computational efficiency, especially when dealing with irregular geometries or the need for adaptive resolution. In the present work we introduce a meshless method that overcomes these challenges by combining the modified finite-particle method (MFPM) for discretization with the Picard iteration technique for solving the result piecewise linear system. The proposed technique employs adaptive stencil selection to guarantee a result linear system with a moderate condition number. An adaptive meshless refinement method enhances the free boundary resolution, particularly in capturing the unknown free boundary a priori. Numerical experiments confirm the method’s flexibility and robustness across a range of node layouts – including Cartesian grids, PNP nodes, and Halton points – demonstrating its potential as an effective tool for solving bilateral obstacle problems and broadening the applicability of PDE and variational inequality models.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"183 ","pages":"Article 106561"},"PeriodicalIF":4.1,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145657536","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-02-01Epub Date: 2025-12-08DOI: 10.1016/j.enganabound.2025.106573
Jinlian Ren , Wei Zhang , Weigang Lu , Tao Jiang
In this paper, a novel stabilization term (ST) is introduced into the nonlinear constitutive model, and a corrected smoothed particle hydrodynamics approximation scheme incorporating the stabilization term, which is termed as C-SPH+ST, is developed to address the High Weissenberg Number Problem (HWNP) in viscoelastic fluids based on the Oldroyd-B constitutive model. To establish a high-accuracy, enhanced-stability algorithm for solving the HWNP of viscoelastic fluids, we perform the improvements based on a corrected SPH method in our previous work which enforces the corrected kernel gradient, particle shifting and density reinitialization techniques in the original SPH framework. Meanwhile, a novel stabilization term, distinct from existing artificial stress terms, is proposed and integrated into the C-SPH+noST method. In numerical experiments, the validity of the proposed stabilized corrected SPH scheme is first verified by simulating classical viscoelastic flows with HWNP, with results compared against reference solutions. The proposed stable particle scheme is then applied to predict the complex nonlinear behavior of viscoelastic Wannier flow, and focusing primarily on the influence of high Weissenberg Number on unstable viscoelastic flow. Additionally, the effects of the parameters in the given stabilization term and other techniques are also analyzed and discussed. All numerical results indicate that the proposed C-SPH+ST algorithm is a robust tool for addressing the HWNP in viscoelastic fluids, and can efficiently prevent non-physical particle clustering and accurately capture the complex phenomena in viscoelastic flows with high Weissenberg numbers.
{"title":"C-SPH+ST: A novel stabilized corrected smoothed particle hydrodynamics scheme for high Weissenberg number problem of viscoelastic fluids","authors":"Jinlian Ren , Wei Zhang , Weigang Lu , Tao Jiang","doi":"10.1016/j.enganabound.2025.106573","DOIUrl":"10.1016/j.enganabound.2025.106573","url":null,"abstract":"<div><div>In this paper, a novel stabilization term (ST) is introduced into the nonlinear constitutive model, and a corrected smoothed particle hydrodynamics approximation scheme incorporating the stabilization term, which is termed as C-SPH+ST, is developed to address the High Weissenberg Number Problem (HWNP) in viscoelastic fluids based on the Oldroyd-B constitutive model. To establish a high-accuracy, enhanced-stability algorithm for solving the HWNP of viscoelastic fluids, we perform the improvements based on a corrected SPH method in our previous work which enforces the corrected kernel gradient, particle shifting and density reinitialization techniques in the original SPH framework. Meanwhile, a novel stabilization term, distinct from existing artificial stress terms, is proposed and integrated into the C-SPH+noST method. In numerical experiments, the validity of the proposed stabilized corrected SPH scheme is first verified by simulating classical viscoelastic flows with HWNP, with results compared against reference solutions. The proposed stable particle scheme is then applied to predict the complex nonlinear behavior of viscoelastic Wannier flow, and focusing primarily on the influence of high Weissenberg Number on unstable viscoelastic flow. Additionally, the effects of the parameters in the given stabilization term and other techniques are also analyzed and discussed. All numerical results indicate that the proposed C-SPH+ST algorithm is a robust tool for addressing the HWNP in viscoelastic fluids, and can efficiently prevent non-physical particle clustering and accurately capture the complex phenomena in viscoelastic flows with high Weissenberg numbers.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"183 ","pages":"Article 106573"},"PeriodicalIF":4.1,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145697382","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-02-01Epub Date: 2025-12-24DOI: 10.1016/j.enganabound.2025.106602
P. Phung-Van , P.T. Hung , M. Abdel Wahab , Chien H. Thai
This study investigates the nanoscale free vibration behavior of a novel smart sandwich nanoplate integrating a functionally graded triply periodic minimal surface (FG-TPMS) core combined with magneto-electro-elastic face sheets. Unlike conventional sandwich structures, this design leverages the mechanical efficiency of FG-TPMS architectures with coupled piezoelectric–piezomagnetic responses of magneto-electro-elastic materials to enhance vibrational performance compared to traditional sandwich composites. The governing equations are formulated using nonlocal strain gradient theory to accurately capture small-scale effects, while isogeometric analysis is employed to ensure high precision and continuity in numerical simulations. Additionally, the displacement approximation is constructed using a newly developed Chebyshev shear deformation theory, which provides improved representation of shear effects in nanoscale plates. The findings demonstrate that synergistic interaction between FG-TPMS topologies, magneto-electro-elastic face sheets and small-scale effects significantly influences natural frequencies. Moreover, this study shows that a larger length scale parameter increases stiffness and raises frequencies, while a higher nonlocal parameter lowers stiffness and reduces frequencies. And the magnetic field strengthens the nanoplate, whereas the electric voltage weakens it. These results offer valuable insights into the dynamic analysis of smart nanostructures with potential applications in aerospace, biomedical engineering and vibration control systems.
{"title":"Piezoelectric and piezomagnetic effects on functionally graded triply periodic minimal surface smart sandwich nanoscale plates using Chebyshev shear deformation theory","authors":"P. Phung-Van , P.T. Hung , M. Abdel Wahab , Chien H. Thai","doi":"10.1016/j.enganabound.2025.106602","DOIUrl":"10.1016/j.enganabound.2025.106602","url":null,"abstract":"<div><div>This study investigates the nanoscale free vibration behavior of a novel smart sandwich nanoplate integrating a functionally graded triply periodic minimal surface (FG-TPMS) core combined with magneto-electro-elastic face sheets. Unlike conventional sandwich structures, this design leverages the mechanical efficiency of FG-TPMS architectures with coupled piezoelectric–piezomagnetic responses of magneto-electro-elastic materials to enhance vibrational performance compared to traditional sandwich composites. The governing equations are formulated using nonlocal strain gradient theory to accurately capture small-scale effects, while isogeometric analysis is employed to ensure high precision and continuity in numerical simulations. Additionally, the displacement approximation is constructed using a newly developed Chebyshev shear deformation theory, which provides improved representation of shear effects in nanoscale plates. The findings demonstrate that synergistic interaction between FG-TPMS topologies, magneto-electro-elastic face sheets and small-scale effects significantly influences natural frequencies. Moreover, this study shows that a larger length scale parameter increases stiffness and raises frequencies, while a higher nonlocal parameter lowers stiffness and reduces frequencies. And the magnetic field strengthens the nanoplate, whereas the electric voltage weakens it. These results offer valuable insights into the dynamic analysis of smart nanostructures with potential applications in aerospace, biomedical engineering and vibration control systems.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"183 ","pages":"Article 106602"},"PeriodicalIF":4.1,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145823739","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号