Yang Yu, Chongshuai Wang, Ean Tat Ooi, Sundararajan Natarajan, Haitian Yang
{"title":"A TPAA-SBFE Based Partitioning Algorithm for Structural Viscoelastic Dynamic Analysis With Cyclic Symmetry","authors":"Yang Yu, Chongshuai Wang, Ean Tat Ooi, Sundararajan Natarajan, Haitian Yang","doi":"10.1002/nme.70053","DOIUrl":"https://doi.org/10.1002/nme.70053","url":null,"abstract":"<div>\u0000 \u0000 <p>An innovative and efficient reduced order algorithm is proposed for the dynamic analysis of 2-D and 3-D viscoelastic structures featuring cyclic symmetry, which is applicable to complete cyclically symmetric structures with or without common nodes, as well as incomplete cyclically symmetric structures. The proposed algorithm is developed on a recursive platform based on TPAA-SBFEM and inherits all the advantages of TPAA-SBFEM. For complete cyclically symmetric structures with or without common nodes, the global stiffness and mass matrices are proved to be block-circulant and can be generated via a cyclically symmetric part instead of the whole domain. And then a recursive partitioning algorithm is rendered, such that the solution scale and computational cost can be efficiently reduced via solving a series of independent problems with small scales. For incomplete cyclically symmetric structures, the complete cyclic symmetry and Woodbury formula are utilized to reduce computational expense. Four numerical examples are provided to elucidate the effectiveness and efficiency of the proposed approach with satisfactory results.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 11","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-05-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144171785","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Topology Optimization of Arc Unit Cells by Extending Energy-Based Homogenization and Parametric Level Set Methods","authors":"Gang Chen, Xiangyu Chen, Qiwen Zeng, Hang Yang","doi":"10.1002/nme.70044","DOIUrl":"https://doi.org/10.1002/nme.70044","url":null,"abstract":"<div>\u0000 \u0000 <p>This paper extends the energy-based homogenization method (EBHM) to evaluate the effective elastic properties of two-dimensional arc unit cells, where the periodic boundary conditions are constructed based on the polar strain–displacement relationships, and the displacement field is solved by the polar finite element method. Combining the extended EBHM with the parametric level set method (PLSM) and the augmented Lagrangian multiplier method (ALM), it is able to design the topological configurations of arc unit cells for seeking superior effective elastic properties. In addition, the level set function is first filtered by compactly supported radial basis functions (CSRBFs) in order to eliminate small holes in the configuration. Numerical examples are performed to demonstrate the advantages of the proposed method. The periodic boundary condition performs well when the circumferential angle of the arc unit cell is in the range of [0.0002, 0.2]. Results indicate that the effective elastic properties of optimized arc unit cells override those of geometrically transformed optimal square unit cells. This superiority mainly benefits from the fact that unisymmetrical configurations are allowed in arc unit cells but not in square unit cells. Results also imply that the capability of nucleation of new holes is retained even though small holes are eliminated by filtering. It contributes to the trade-off between performance and manufacturability.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 11","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-05-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144171784","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Zhuo Deng, Shi-Xuan Liu, Song Cen, Ming Sun, Yan Shang
{"title":"Corotational Unsymmetric Membrane Element Formulation for Geometric Nonlinear Analysis of Flexoelectric Solids Within the Consistent Couple Stress Theory","authors":"Zhuo Deng, Shi-Xuan Liu, Song Cen, Ming Sun, Yan Shang","doi":"10.1002/nme.70045","DOIUrl":"https://doi.org/10.1002/nme.70045","url":null,"abstract":"<div>\u0000 \u0000 <p>In this work, the corotational (CR) penalty membrane element formulation for geometric nonlinear analysis of flexoelectric solids with large displacement but small deformation is proposed. To achieve this, the two-dimensional couple stress-based flexoelectric model is first extended to the CR form where the effect of rigid body rotation is excluded and, accordingly, the curvature that governs the electromechanical behavior of flexoelectricity is defined as the derivatives of the elastic rotation in the CR configuration. Next, a quadrilateral 8-node element is constructed based on the CR method. In the development, with the use of the penalty function method, the independent nodal rotation degrees of freedom are utilized to approximate the elastic rotation field, ensuring the higher-order continuity requirement for the displacement test function in a weak sense. In addition, to fully improve the element performance, the trial function of CR force-stress is formulated based on the polynomial force-stress functions, which are continually updated with reference to the latest deformed configuration. As demonstrated by the numerical benchmark examples, the new element can efficiently and accurately predict the geometric nonlinear electromechanical response of slender flexoelectric structures and capture the size effect very well.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 11","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-05-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144171591","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Qifan Zhou, Yingqing Guo, Guicai Li, Kejie Xu, Kun Wang
{"title":"Multimodal Hybrid Aero-Engine Mechanical Wear Fault Diagnosis Algorithm Based on Two-Channel Data Input Types","authors":"Qifan Zhou, Yingqing Guo, Guicai Li, Kejie Xu, Kun Wang","doi":"10.1002/nme.70046","DOIUrl":"https://doi.org/10.1002/nme.70046","url":null,"abstract":"<div>\u0000 \u0000 <p>Aero-engines are complex and sophisticated systems combining mechanical, thermal, and fluidic domains. Abnormal wear of mechanical components is becoming more prevalent due to severe changes in flight conditions and the external environment, which may lead to drastic performance degradation and accidents. Therefore, the diagnosis of such wear and tear faults is urgent, and based on this need, more researchers and scholars are focusing their attention on it. To address the current shortcomings of fault diagnosis algorithms that only rely on one-dimensional datasets or two-dimensional image analysis and the low accuracy of final fault identification, an innovative hybrid algorithm is proposed in this study. The algorithm integrates one-dimensional time series data and two-dimensional image data, converts the one-dimensional dataset into a two-dimensional image dataset through the Gramian Angle Field technique, and subsequently uses a dual-channel GRU-CNN (Convolutional Neural Network-Gated Recurrent Unit) algorithm model designed for fault diagnosis, which can simultaneously analyze and map the features and fault modes of the one-dimensional dataset and the two-dimensional image. In order to extract features with richer semantic information and stronger discriminative ability, a multimodal fusion technique is employed, which successfully addresses the limitations of the wear-and-tear feature distributions of the two datasets using the cross-extraction fusion method and combines the advantages of both in terms of trend distributions of the time series and edge feature distributions of the image sequences, respectively. The best fault diagnosis results were achieved by using the strong mapping relationship between the saliency feature expression and the fault modes. The final analysis shows that the recognition rate of typical mechanical wear of aero-engines exceeds 97%, thus achieving the desired goal.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 10","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-05-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144126027","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Shape Optimization of Supercapacitor Electrode to Maximize Charge Storage","authors":"Jiajie Li, Shenggao Zhou, Shengfeng Zhu","doi":"10.1002/nme.70052","DOIUrl":"https://doi.org/10.1002/nme.70052","url":null,"abstract":"<div>\u0000 \u0000 <p>This work proposes a shape optimization approach for electrode morphology to maximize charge storage in supercapacitors. The ionic distributions and electric potential are modeled by the steady-state Poisson–Nernst–Planck system. Shape sensitivity analysis is performed to derive the Eulerian derivative in both volumetric and boundary expressions. An optimal electrode morphology is obtained through gradient flow algorithms. The steady-state Poisson–Nernst–Planck system is efficiently solved by the Gummel fixed-point scheme with finite-element discretization, in which exponential coefficients with large variation are tackled with inverse averaging techniques. Extensive numerical experiments are performed to demonstrate the effectiveness of the proposed optimization model and corresponding numerical methods in enhancing charge storage in supercapacitors. It is expected that the proposed shape optimization approach provides a promising tool in the design of electrode morphology from a perspective of charge storage enhancement.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 10","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-05-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144117960","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Matteo Giacomini, Davide Cortellessa, Luan M. Vieira, Ruben Sevilla, Antonio Huerta
{"title":"A Hybrid Pressure Formulation of the Face-Centred Finite Volume Method for Viscous Laminar Incompressible Flows","authors":"Matteo Giacomini, Davide Cortellessa, Luan M. Vieira, Ruben Sevilla, Antonio Huerta","doi":"10.1002/nme.70037","DOIUrl":"https://doi.org/10.1002/nme.70037","url":null,"abstract":"<p>This work presents a hybrid pressure face-centred finite volume (FCFV) solver to simulate steady-state incompressible Navier-Stokes flows. The method leverages the robustness, in the incompressible limit, of the hybridisable discontinuous Galerkin paradigm for compressible and weakly compressible flows to derive the formulation of a novel, low-order face-based discretization. The incompressibility constraint is enforced in a weak sense by introducing an inter-cell mass flux, defined in terms of a new, hybrid variable that represents the pressure at the cell faces. This results in a new hybridization strategy where cell variables (velocity, pressure, and deviatoric strain rate tensor) are expressed as a function of velocity and pressure at the barycentre of the cell faces. The hybrid pressure formulation provides first-order convergence of all variables, including the stress, without the need for gradient reconstruction, thus being less sensitive to cell type, stretching, distortion, and skewness than traditional low-order finite volume solvers. Numerical benchmarks of Navier-Stokes flows at low and moderate Reynolds numbers, in two and three dimensions, are presented to evaluate the accuracy and robustness of the method. In particular, the hybrid pressure formulation outperforms the FCFV method when convective effects are relevant, achieving accurate predictions on significantly coarser meshes.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 10","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-05-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.70037","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144100881","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Thermodynamically Consistent Modeling of Compressible Two-Phase Flow in Porous Media","authors":"Junkai Wang, Xiaolin Zhou, Qiaolin He","doi":"10.1002/nme.70050","DOIUrl":"https://doi.org/10.1002/nme.70050","url":null,"abstract":"<div>\u0000 \u0000 <p>In this paper, we develop a thermodynamically consistent model for compressible two-phase flow in porous media according to the Reynolds transport theorem and the second law of thermodynamics. The Helmholtz free energy and saturation free energy are introduced to characterize fluid compressibility and capillary pressure effects, respectively. Especially, we use molar density instead of pressure as one of the primary variables and derive a Darcy-type momentum equation with chemical potential gradient as the main driving force. Furthermore, the discrete total free energy dissipation and prior error estimates of the proposed numerical scheme are derived. Numerical results are presented to validate the accuracy, stability, and effectiveness of our proposed method.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 10","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-05-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144091844","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A Unified Algorithm Framework for Phase-Field Based Multimaterial Topology Optimization With Various Graded Interfaces","authors":"Qian Yu, Qing Xia, Yibao Li, Chao Yang","doi":"10.1002/nme.70048","DOIUrl":"https://doi.org/10.1002/nme.70048","url":null,"abstract":"<div>\u0000 \u0000 <p>In this work, we present a unified algorithm framework for multimaterial topology optimization with various graded interfaces based on the phase field method. In the framework, the topology optimization problem is transformed into solving a set of linear partial differential equations (PDEs). We then introduce a scalar variable to grade the material property, resulting in a complex nonlinear interpolation operator in the effective elastic tensor. Following that, we define an auxiliary variable to substitute all the nonlinear transformations for a linear elliptic equation system. A second-order accurate Crank–Nicolson scheme is applied on the reformulated system, in which all high-order nonlinear terms are treated in a semi-explicit fashion. We prove that the proposed algorithm framework is unconditionally energy stable and demonstrate its stability as well as accuracy by numerical examples. A series of benchmark problems with different interfacial behaviors in topological design are investigated to verify the effectiveness of our method. The sensitivity of different parameters in the model is analyzed to evaluate their effects on the resulting structure.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 10","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-05-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144091845","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"The Force Method Version of Goodman's Joint Element with Convergence Proof","authors":"Hong Zheng, Jia-han Jiang, Ming-kai Sun","doi":"10.1002/nme.70047","DOIUrl":"https://doi.org/10.1002/nme.70047","url":null,"abstract":"<div>\u0000 \u0000 <p>Interfaces widely present in nature and civil engineering are the most fundamental elements causing discontinuous deformation and failure. The Goodman joint element is incorporated into the vast majority of commercial and proprietary finite element programs due to its simplicity. However, the numerical properties of the Goodman element are not ideal. A lot of effort has gone into enhancing its numerical properties, only to make some empirical suggestions for adjusting mesh and computation parameters, yet there are still many examples that fail to reach convergence or incur drastic numerical oscillations of contact stress. Considering that the Goodman element is implemented within the framework of the displacement method, and poorlyposed problems in the displacement method are usually well-posed problems in the force method, this study proposes a force method version of the Goodman element, abbreviated as FMVGE. The primal unknowns in FMVGE are the contact stresses on the interface rather than nodal deformations. The core problem in FMVGE represents the contact conditions in the form of a quasi-variational inequality (QVI). By employing process iteration rather than state iteration used in the displacement method, at the same time, FMVGE precisely satisfies the contact conditions, with convergence being theoretically guaranteed. Analysis of classic examples and engineering cases indicates that the numerical properties of FMVGE are significantly superior to the widely adopted master–slave block method. Appendices A and B provide the proof of the convergence of the FMVGE solution and the core Matlab code for solving the QVI, respectively.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 10","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-05-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143950296","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Improved Conservative Phase Field Method for Contact Line Motion With Solidification","authors":"Mingguang Shen, Ben Q. Li","doi":"10.1002/nme.70036","DOIUrl":"https://doi.org/10.1002/nme.70036","url":null,"abstract":"<p>Phase field method is appealing due to its implicit capture of the interface. It mainly contains two types, one being the Allen–Cahn and the other being the Cahn-Hilliard. The latter conserves mass, but is fourth order, not as amenable to numerical methods as the former. Recently, a conservative Allen–Cahn phase field method has been put forward, improving the boundedness of the phase field. The new method conserves mass well and is only second order. This paper couples the conservative Allen–Cahn phase field model and the Navier–Stokes equation and the heat balance equation. Interfacial tension is expressed in a potential form. The model is discretized using a finite difference method on a half-staggered grid, and the velocity–pressure coupling is decoupled using an explicit projection method. The heat balance equation is solved using an iterative method. The model was tested against the classic static bubble case and was applied to impact in various conditions. Moreover, the model is parallelized using shared memory parallelism, OpenMP.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 9","pages":""},"PeriodicalIF":2.7,"publicationDate":"2025-05-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.70036","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143944930","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}