Pub Date : 2025-02-12DOI: 10.1016/j.cma.2025.117809
Riccardo Alberini , Michele Terzano , Gerhard A. Holzapfel , Andrea Spagnoli
Advanced simulations of the mechanical behavior of soft tissues frequently rely on structure-based constitutive models, including smeared descriptions of collagen fibers. Among them, the so-called Discrete Fiber Dispersion (DFD) modeling approach is based on a discrete integration of the fiber-strain energy over all the fiber directions. In this paper, we review the theoretical framework of the DFD model, including a derivation of the stress and stiffness tensors required for the finite element implementation. Specifically, their expressions for incompressible plane stress problems are obtained. The use of a Lebedev quadrature, built exploiting the octahedral symmetry, is then proposed, illustrating the particular choice adopted for the orientation of the integration points. Next, the convergence of this quadrature scheme is assessed by means of three numerical benchmark tests, highlighting the advantages with respect to other angular integration methods available in the literature. Finally, using the implemented model, we analyze the mechanical properties of the Z-plasty, a technique commonly used in reconstructive skin surgery, considering multiple geometrical configurations, orientations of the fibers, and levels of skin prestress. The results are presented in the form of mechanical quantities relevant to surgical practice.
{"title":"A discrete fiber dispersion model with octahedral symmetry quadrature for mechanical analyses of skin corrective surgeries","authors":"Riccardo Alberini , Michele Terzano , Gerhard A. Holzapfel , Andrea Spagnoli","doi":"10.1016/j.cma.2025.117809","DOIUrl":"10.1016/j.cma.2025.117809","url":null,"abstract":"<div><div>Advanced simulations of the mechanical behavior of soft tissues frequently rely on structure-based constitutive models, including smeared descriptions of collagen fibers. Among them, the so-called Discrete Fiber Dispersion (DFD) modeling approach is based on a discrete integration of the fiber-strain energy over all the fiber directions. In this paper, we review the theoretical framework of the DFD model, including a derivation of the stress and stiffness tensors required for the finite element implementation. Specifically, their expressions for incompressible plane stress problems are obtained. The use of a Lebedev quadrature, built exploiting the octahedral symmetry, is then proposed, illustrating the particular choice adopted for the orientation of the integration points. Next, the convergence of this quadrature scheme is assessed by means of three numerical benchmark tests, highlighting the advantages with respect to other angular integration methods available in the literature. Finally, using the implemented model, we analyze the mechanical properties of the Z-plasty, a technique commonly used in reconstructive skin surgery, considering multiple geometrical configurations, orientations of the fibers, and levels of skin prestress. The results are presented in the form of mechanical quantities relevant to surgical practice.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117809"},"PeriodicalIF":6.9,"publicationDate":"2025-02-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143386409","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-02-12DOI: 10.1016/j.cma.2025.117812
Jingjing Yang , Pei Zhou , Lin Lan , Chun-Gang Zhu
We propose a novel isogeometric collocation method (IGA-C) for trimmed domains, using weighted extended B-splines (WEB-splines). Our approach employs the implicit representations of the trimming boundaries to construct the weighted basis, which allows for the subsequent calculations based on a parametrization over a single tensor-product patch, despite the nontrivial shape of the domain. The stabilization of the weighted basis is accomplished by means of extension. We present the classification criterion and the calculation procedure for the extension coefficients of the inner B-splines. The utilization of WEB-splines in analysis enables a natural application of the Dirichlet boundary conditions. Additionally, we adopt a hybrid collocation–Galerkin approach to impose the Neumann boundary conditions on the trimming boundaries. Our proposed method combines the advantages of WEB-splines and IGA-C in terms of straightforward implementation, well-conditioned system matrices and high computational efficiency, which we illustrate by numerical tests on 2D and 3D trimmed geometries. The numerical results further demonstrate that our methodology guarantees the same convergence rates as IGA-C.
{"title":"A hybrid isogeometric collocation method on implicitly trimmed domains","authors":"Jingjing Yang , Pei Zhou , Lin Lan , Chun-Gang Zhu","doi":"10.1016/j.cma.2025.117812","DOIUrl":"10.1016/j.cma.2025.117812","url":null,"abstract":"<div><div>We propose a novel isogeometric collocation method (IGA-C) for trimmed domains, using weighted extended B-splines (WEB-splines). Our approach employs the implicit representations of the trimming boundaries to construct the weighted basis, which allows for the subsequent calculations based on a parametrization over a single tensor-product patch, despite the nontrivial shape of the domain. The stabilization of the weighted basis is accomplished by means of extension. We present the classification criterion and the calculation procedure for the extension coefficients of the inner B-splines. The utilization of WEB-splines in analysis enables a natural application of the Dirichlet boundary conditions. Additionally, we adopt a hybrid collocation–Galerkin approach to impose the Neumann boundary conditions on the trimming boundaries. Our proposed method combines the advantages of WEB-splines and IGA-C in terms of straightforward implementation, well-conditioned system matrices and high computational efficiency, which we illustrate by numerical tests on 2D and 3D trimmed geometries. The numerical results further demonstrate that our methodology guarantees the same convergence rates as IGA-C.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117812"},"PeriodicalIF":6.9,"publicationDate":"2025-02-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143386412","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-02-12DOI: 10.1016/j.cma.2025.117808
Jonathan Stollberg , Tarun Gangwar , Oliver Weeger , Dominik Schillinger
We present a new framework for the simultaneous optimization of both the topology as well as the relative density grading of cellular structures and materials, also known as lattices. Due to manufacturing constraints, the optimization problem falls into the class of mixed-integer nonlinear programming problems. Since no algorithm is capable of solving these problems in polynomial time, we obtain a relaxed problem from a multiplicative split of the relative density and a penalization approach. The sensitivities of the objective function are derived such that any gradient-based solver might be applied for the iterative update of the design variables. In a next step, we introduce a material model that is parametric in the design variables of interest and suitable to describe the isotropic deformation behavior of quasi-stochastic lattices. For that, we derive and implement further physical constraints and enhance a physics-augmented neural network from the literature that was formulated initially for rhombic materials. Finally, to illustrate the applicability of the method, we incorporate the material model into our computational framework and exemplary optimize two-and three-dimensional benchmark structures as well as a complex aircraft component.
{"title":"Multiscale topology optimization of functionally graded lattice structures based on physics-augmented neural network material models","authors":"Jonathan Stollberg , Tarun Gangwar , Oliver Weeger , Dominik Schillinger","doi":"10.1016/j.cma.2025.117808","DOIUrl":"10.1016/j.cma.2025.117808","url":null,"abstract":"<div><div>We present a new framework for the simultaneous optimization of both the topology as well as the relative density grading of cellular structures and materials, also known as lattices. Due to manufacturing constraints, the optimization problem falls into the class of mixed-integer nonlinear programming problems. Since no algorithm is capable of solving these problems in polynomial time, we obtain a relaxed problem from a multiplicative split of the relative density and a penalization approach. The sensitivities of the objective function are derived such that any gradient-based solver might be applied for the iterative update of the design variables. In a next step, we introduce a material model that is parametric in the design variables of interest and suitable to describe the isotropic deformation behavior of quasi-stochastic lattices. For that, we derive and implement further physical constraints and enhance a physics-augmented neural network from the literature that was formulated initially for rhombic materials. Finally, to illustrate the applicability of the method, we incorporate the material model into our computational framework and exemplary optimize two-and three-dimensional benchmark structures as well as a complex aircraft component.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117808"},"PeriodicalIF":6.9,"publicationDate":"2025-02-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143386408","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
This paper proposes and studies two extensions of applying hp-variational physics-informed neural networks, more precisely the FastVPINNs framework, to convection-dominated convection–diffusion–reaction problems. First, a term in the spirit of a SUPG stabilization is included in the loss functional and a network architecture is proposed that predicts spatially varying stabilization parameters. Having observed that the selection of the indicator function in hard-constrained Dirichlet boundary conditions has a big impact on the accuracy of the computed solutions, the second novelty is the proposal of a network architecture that learns good parameters for a class of indicator functions. Numerical studies show that both proposals lead to noticeably more accurate results than approaches that can be found in the literature.
{"title":"Improving hp-variational physics-informed neural networks for steady-state convection-dominated problems","authors":"Thivin Anandh , Divij Ghose , Himanshu Jain , Pratham Sunkad , Sashikumaar Ganesan , Volker John","doi":"10.1016/j.cma.2025.117797","DOIUrl":"10.1016/j.cma.2025.117797","url":null,"abstract":"<div><div>This paper proposes and studies two extensions of applying hp-variational physics-informed neural networks, more precisely the FastVPINNs framework, to convection-dominated convection–diffusion–reaction problems. First, a term in the spirit of a SUPG stabilization is included in the loss functional and a network architecture is proposed that predicts spatially varying stabilization parameters. Having observed that the selection of the indicator function in hard-constrained Dirichlet boundary conditions has a big impact on the accuracy of the computed solutions, the second novelty is the proposal of a network architecture that learns good parameters for a class of indicator functions. Numerical studies show that both proposals lead to noticeably more accurate results than approaches that can be found in the literature.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117797"},"PeriodicalIF":6.9,"publicationDate":"2025-02-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143386407","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
A finite-element method dependent adjoint-based procedure to determine the temperature field of structures based on measured displacements/strains and a set of standard loads is developed and tested. Given a series of force and deformation measurements, the temperature field is obtained by minimizing the adequately weighted differences between the measured and computed values. Three numerical examples — a Plate With a Hole, a Bridge, and a Hoover Dam example — each with multiple sensors distributed in different configurations, demonstrate the procedure’s capabilities. A target temperature distribution is prescribed in all cases, and the displacement sensor data is recorded. The optimization algorithm (here, steepest descent with Barzilai–Borwein step) uses this data to optimize the temperatures such that the same deformation is obtained at the sensor locations. Vertex Morphing is used as a filter to mitigate the ill-conditioning. Results show that the proposed approach can accurately reconstruct the target thermal distribution, especially when more sensors are used. Additionally, it is observed that the sensors do not need to be positioned in the region of interest; the method remains effective as long as the sensors can detect changes related to that area. A comparison with standard spatial interpolation techniques, namely, k-nearest neighbors and ordinary and universal kriging, is performed using temperature sensors in the same configurations. The proposed approach performs remarkably better than the interpolation techniques with a reduction in the norm of up to 41.3%, 93.9%, and 41.3%, for the Plate With a Hole, the Bridge, and the Dam examples, respectively.
{"title":"Adjoint-based recovery of thermal fields from displacement or strain measurements","authors":"Talhah Shamshad Ali Ansari , Rainald Löhner , Roland Wüchner , Harbir Antil , Suneth Warnakulasuriya , Ihar Antonau , Facundo Airaudo","doi":"10.1016/j.cma.2025.117818","DOIUrl":"10.1016/j.cma.2025.117818","url":null,"abstract":"<div><div>A finite-element method dependent adjoint-based procedure to determine the temperature field of structures based on measured displacements/strains and a set of standard loads is developed and tested. Given a series of force and deformation measurements, the temperature field is obtained by minimizing the adequately weighted differences between the measured and computed values. Three numerical examples — a Plate With a Hole, a Bridge, and a Hoover Dam example — each with multiple sensors distributed in different configurations, demonstrate the procedure’s capabilities. A target temperature distribution is prescribed in all cases, and the displacement sensor data is recorded. The optimization algorithm (here, steepest descent with Barzilai–Borwein step) uses this data to optimize the temperatures such that the same deformation is obtained at the sensor locations. Vertex Morphing is used as a filter to mitigate the ill-conditioning. Results show that the proposed approach can accurately reconstruct the target thermal distribution, especially when more sensors are used. Additionally, it is observed that the sensors do not need to be positioned in the region of interest; the method remains effective as long as the sensors can detect changes related to that area. A comparison with standard spatial interpolation techniques, namely, k-nearest neighbors and ordinary and universal kriging, is performed using temperature sensors in the same configurations. The proposed approach performs remarkably better than the interpolation techniques with a reduction in the <span><math><mrow><mi>L</mi><mn>2</mn></mrow></math></span> norm of up to 41.3%, 93.9%, and 41.3%, for the Plate With a Hole, the Bridge, and the Dam examples, respectively.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117818"},"PeriodicalIF":6.9,"publicationDate":"2025-02-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143386410","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-02-12DOI: 10.1016/j.cma.2025.117807
R. Fu , D. Xiao , A.G. Buchan , X. Lin , Y. Feng , G. Dong
Reduced order modelling is popular and widely used in engineering as it has the potential to gain several orders of magnitude CPU speedup for simulations with different parameters such as different initial or boundary conditions. This work presents a new parametric non-linear non-intrusive reduced-order model (P-NLNIROM) for the fluid problems, which extends the capabilities for non-linear NIROM Fuet al. (2023) on parametric problems. The Deep Auto-Encoder (DAE), Deep Residual Learning Neural network (ResNet), and deep transfer learning are used to construct the P-NLNIROM. The DAE is developed to project or map the original high-dimensional dynamical systems into a much lower dimensional nonlinear reduced latent space and the deep ResNet is used to construct a set of functions that represents the relationships between model input parameters (such as initial or boundary conditions) and extracted representations in the latent or reduced space (reduced representation of fluid dynamics). Transfer learning is used to extend the predictive capability of the model for different parameters. The novelty of this work lies in that ResNet and transfer learning are used to predict different parametric conditions for the AutoEncoder-based, non-linear, non-intrusive reduced-order model (NLNIROM). Transfer learning expands the predictive range of parametric space and makes the transferred P-NLNIROM perform well with much less data. The capability of this new P-NLNIROM is illustrated numerically by two test cases: a lock exchange, and a flow past a cylinder. The results obtained show that the P-NLNIROM performs well and the transferred model shows more promising results under new parameter conditions.
{"title":"A parametric non-linear non-intrusive reduce-order model using deep transfer learning","authors":"R. Fu , D. Xiao , A.G. Buchan , X. Lin , Y. Feng , G. Dong","doi":"10.1016/j.cma.2025.117807","DOIUrl":"10.1016/j.cma.2025.117807","url":null,"abstract":"<div><div>Reduced order modelling is popular and widely used in engineering as it has the potential to gain several orders of magnitude CPU speedup for simulations with different parameters such as different initial or boundary conditions. This work presents a new parametric non-linear non-intrusive reduced-order model (P-NLNIROM) for the fluid problems, which extends the capabilities for non-linear NIROM Fuet al. (2023) on parametric problems. The Deep Auto-Encoder (DAE), Deep Residual Learning Neural network (ResNet), and deep transfer learning are used to construct the P-NLNIROM. The DAE is developed to project or map the original high-dimensional dynamical systems into a much lower dimensional nonlinear reduced latent space and the deep ResNet is used to construct a set of functions that represents the relationships between model input parameters (such as initial or boundary conditions) and extracted representations in the latent or reduced space (reduced representation of fluid dynamics). Transfer learning is used to extend the predictive capability of the model for different parameters. The novelty of this work lies in that ResNet and transfer learning are used to predict different parametric conditions for the AutoEncoder-based, non-linear, non-intrusive reduced-order model (NLNIROM). Transfer learning expands the predictive range of parametric space and makes the transferred P-NLNIROM perform well with much less data. The capability of this new P-NLNIROM is illustrated numerically by two test cases: a lock exchange, and a flow past a cylinder. The results obtained show that the P-NLNIROM performs well and the transferred model shows more promising results under new parameter conditions.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117807"},"PeriodicalIF":6.9,"publicationDate":"2025-02-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143386411","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-02-11DOI: 10.1016/j.cma.2025.117816
Francesco Duronio , Andrea Di Mascio
Many fluid flows of engineering interest involve high-density ratios, the formation of shock waves, and the presence of multiple chemical species. In these situations, obtaining accurate results depends on the computational fluid dynamics (CFD) algorithm’s ability to ensure stable and robust time and space discretization of the governing equations, effectively minimizing numerical diffusion.
In this work, we develop and evaluate a low-dissipation, high-order solver for the Navier–Stokes equations applicable to compressible, multi-species flows within the OpenFOAM library. We introduce a new shock sensor specifically designed for high-order discretization methods, which guarantees accuracy up to the fourth order in both space and time when using unstructured grids. Additionally, this approach is compatible with real-gas equations of state.
We evaluated the solver performances on three test cases, each chosen for its relevance to real-world engineering problems: a 1D multi-specie shock tube, a 2D shock tube where Richtmyer–Meshkov instability develops as a consequence of bubble-shock interaction, and, finally, the flow through the complex 3D geometry of a high-pressure fuel injector used in propulsion applications, to investigate its performances on unstructured grids. The results confirm the effectiveness of the proposed sensor in scenarios characterized by field discontinuities, shock waves, high-density ratio, and distorted grids. This means that our solver can accurately simulate complex fluid flows in engineering applications where these conditions are often encountered.
The developed high-order scheme and shock sensor were implemented in an OpenFoam solver called rhoCubic4kFoam.
{"title":"Development of an efficient shock sensor for high-order multi-species compressible flow solvers on unstructured grids","authors":"Francesco Duronio , Andrea Di Mascio","doi":"10.1016/j.cma.2025.117816","DOIUrl":"10.1016/j.cma.2025.117816","url":null,"abstract":"<div><div>Many fluid flows of engineering interest involve high-density ratios, the formation of shock waves, and the presence of multiple chemical species. In these situations, obtaining accurate results depends on the computational fluid dynamics (CFD) algorithm’s ability to ensure stable and robust time and space discretization of the governing equations, effectively minimizing numerical diffusion.</div><div>In this work, we develop and evaluate a low-dissipation, high-order solver for the Navier–Stokes equations applicable to compressible, multi-species flows within the OpenFOAM library. We introduce a new shock sensor specifically designed for high-order discretization methods, which guarantees accuracy up to the fourth order in both space and time when using unstructured grids. Additionally, this approach is compatible with real-gas equations of state.</div><div>We evaluated the solver performances on three test cases, each chosen for its relevance to real-world engineering problems: a 1D multi-specie shock tube, a 2D shock tube where Richtmyer–Meshkov instability develops as a consequence of bubble-shock interaction, and, finally, the flow through the complex 3D geometry of a high-pressure fuel injector used in propulsion applications, to investigate its performances on unstructured grids. The results confirm the effectiveness of the proposed sensor in scenarios characterized by field discontinuities, shock waves, high-density ratio, and distorted grids. This means that our solver can accurately simulate complex fluid flows in engineering applications where these conditions are often encountered.</div><div>The developed high-order scheme and shock sensor were implemented in an OpenFoam solver called <em>rhoCubic4kFoam</em>.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117816"},"PeriodicalIF":6.9,"publicationDate":"2025-02-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143386404","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-02-11DOI: 10.1016/j.cma.2025.117811
Maria Strazzullo , Francesco Ballarin , Traian Iliescu , Tomás Chacón Rebollo
The evolve-filter (EF) model is a filter-based numerical stabilization for under-resolved convection-dominated flows. EF is a simple, modular, and effective strategy for both full-order models (FOMs) and reduced-order models (ROMs). It is well-known, however, that when the filter radius is too large, EF can be overdiffusive and yield inaccurate results. To alleviate this, EF is usually supplemented with a relaxation step. The relaxation parameter, however, is very sensitive with respect to the model parameters. In this paper, we propose a novel strategy to alleviate the EF overdiffusivity. Specifically, we leverage the variational multiscale (VMS) framework to separate the large resolved scales from the small resolved scales in the evolved velocity, and we use the filtered small scales to correct the large scales. Furthermore, in the new VMS-EF strategy, we use two different approaches to decompose the evolved velocity: the VMS Evolve-Filter-Filter-Correct (VMS-EFFC) and the VMS Evolve-Postprocess-Filter-Correct (VMS-EPFC) algorithms. The new VMS-based algorithms yield significantly more accurate results than the standard EF in both the FOM and the ROM simulations of a flow past a cylinder at Reynolds number Re = 1000.
{"title":"Variational multiscale evolve and filter strategies for convection-dominated flows","authors":"Maria Strazzullo , Francesco Ballarin , Traian Iliescu , Tomás Chacón Rebollo","doi":"10.1016/j.cma.2025.117811","DOIUrl":"10.1016/j.cma.2025.117811","url":null,"abstract":"<div><div>The evolve-filter (EF) model is a filter-based numerical stabilization for under-resolved convection-dominated flows. EF is a simple, modular, and effective strategy for both full-order models (FOMs) and reduced-order models (ROMs). It is well-known, however, that when the filter radius is too large, EF can be overdiffusive and yield inaccurate results. To alleviate this, EF is usually supplemented with a relaxation step. The relaxation parameter, however, is very sensitive with respect to the model parameters. In this paper, we propose a novel strategy to alleviate the EF overdiffusivity. Specifically, we leverage the variational multiscale (VMS) framework to separate the large resolved scales from the small resolved scales in the evolved velocity, and we use the filtered small scales to correct the large scales. Furthermore, in the new VMS-EF strategy, we use two different approaches to decompose the evolved velocity: the VMS Evolve-Filter-Filter-Correct (VMS-EFFC) and the VMS Evolve-Postprocess-Filter-Correct (VMS-EPFC) algorithms. The new VMS-based algorithms yield significantly more accurate results than the standard EF in both the FOM and the ROM simulations of a flow past a cylinder at Reynolds number Re = 1000.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117811"},"PeriodicalIF":6.9,"publicationDate":"2025-02-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143386403","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-02-11DOI: 10.1016/j.cma.2025.117814
Chen Yang , Qianqian Yu
In force reconstruction, multi-source incomplete information makes it difficult for traditional methods to model and solve the problem accurately, especially with noise or measurement errors. Inspired by multi-objective optimization, this paper proposes a novel set theory-based regularization approach (STR) to enhance adaptability to uncertainties and improve reconstruction accuracy and robustness. The nominal force inversion is constituted and then extended into the interval uncertainty framework, and an effective orthogonal sampling-based interval prediction method is proposed to analyze the coupled effect of force inversion and uncertainty propagation. Once the uncertainty level of the complex structure is known, this prediction method can accurately and quickly estimate the fluctuation bound of the identified force. Enlightened by the completely same logic between the regularization method used in force inversion and multi-objective optimization problem, namely, simultaneously satisfying the norm minimization of the solution and the residual parameter, this study develops a novel multi-objective optimization-inspired set theory-based regularization parameter selection method. This method incorporates the interval dominance relationship to select the most competitive regularization parameter under interval uncertainties. Therefore, an accurate reconstuction framework with bounded uncertainties is finally proposed and verified by two numerical examples.
{"title":"Multi-objective optimization-inspired set theory-based regularization approach for force reconstruction with bounded uncertainties","authors":"Chen Yang , Qianqian Yu","doi":"10.1016/j.cma.2025.117814","DOIUrl":"10.1016/j.cma.2025.117814","url":null,"abstract":"<div><div>In force reconstruction, multi-source incomplete information makes it difficult for traditional methods to model and solve the problem accurately, especially with noise or measurement errors. Inspired by multi-objective optimization, this paper proposes a novel set theory-based regularization approach (STR) to enhance adaptability to uncertainties and improve reconstruction accuracy and robustness. The nominal force inversion is constituted and then extended into the interval uncertainty framework, and an effective orthogonal sampling-based interval prediction method is proposed to analyze the coupled effect of force inversion and uncertainty propagation. Once the uncertainty level of the complex structure is known, this prediction method can accurately and quickly estimate the fluctuation bound of the identified force. Enlightened by the completely same logic between the regularization method used in force inversion and multi-objective optimization problem, namely, simultaneously satisfying the norm minimization of the solution and the residual parameter, this study develops a novel multi-objective optimization-inspired set theory-based regularization parameter selection method. This method incorporates the interval dominance relationship to select the most competitive regularization parameter under interval uncertainties. Therefore, an accurate reconstuction framework with bounded uncertainties is finally proposed and verified by two numerical examples.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"438 ","pages":"Article 117814"},"PeriodicalIF":6.9,"publicationDate":"2025-02-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143386401","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-02-09DOI: 10.1016/j.cma.2025.117825
Mohamed Abdel-Basset , Reda Mohamed , Mohamed Abouhawwash
This study presents a new nature-inspired metaheuristic algorithm, known as the fungal growth optimizer (FGO), which is inspired by fungal growth behavior in nature. Fungal growth behavior includes hyphal growth, branching, and spore germination. Hyphal growth behavior replicates hyphal extension and chemotropism to precisely explore the search space, and thus, reach and exploit nutrient-rich regions. This behavior provides a variety of search patterns in FGO, promoting its performance against stagnation into local optima and slow convergence speed. The branching behavior replicates how new hyphal branches from the side of an existing hypha explore the surrounding regions in search of more nutrients, promoting the exploratory operator throughout the optimization process. The final behavior is spore germination, which represents how existing hyphae explore new environments to reach safer and nutrient-rich areas. When spores land in an environment that is rich in moisture and nutrition, they germinate and grow. FGO assumes that spores will land in a random position at the beginning of the optimization process to promote the exploratory operator. As the optimization process is exceeded, this random position is transformed into a position between the best-so-far solution and a random position, promoting the exploitation operator while preventing premature convergence. FGO is evaluated against four well-known Congress on Evolutionary Computation (CEC) benchmarks (CEC2020, CEC2017, CEC2014, and CEC2022) and eleven engineering design problems. In addition, it is compared with fifteen recently proposed algorithms and eleven highly-performing algorithms, such as L-SHADE, LSHADE-cnEpSin, AL-SHADE, mantis search algorithm (MSA), IMODE, AGSK, SOMA_T3A, HyDE-DF, modified LSHADE-SPACMA, SHADE, and LSHADE-SPACMA, to demonstrate its superiority. According to the experimental results, FGO outperforms or is competitive with all of the compared algorithms for the majority of the test functions, implying that it is a high-performing optimizer and a powerful alternative technique for dealing with complex optimization problems. The FGO source code is available on this link
{"title":"Fungal growth optimizer: A novel nature-inspired metaheuristic algorithm for stochastic optimization","authors":"Mohamed Abdel-Basset , Reda Mohamed , Mohamed Abouhawwash","doi":"10.1016/j.cma.2025.117825","DOIUrl":"10.1016/j.cma.2025.117825","url":null,"abstract":"<div><div>This study presents a new nature-inspired metaheuristic algorithm, known as the fungal growth optimizer (FGO), which is inspired by fungal growth behavior in nature. Fungal growth behavior includes hyphal growth, branching, and spore germination. Hyphal growth behavior replicates hyphal extension and chemotropism to precisely explore the search space, and thus, reach and exploit nutrient-rich regions. This behavior provides a variety of search patterns in FGO, promoting its performance against stagnation into local optima and slow convergence speed. The branching behavior replicates how new hyphal branches from the side of an existing hypha explore the surrounding regions in search of more nutrients, promoting the exploratory operator throughout the optimization process. The final behavior is spore germination, which represents how existing hyphae explore new environments to reach safer and nutrient-rich areas. When spores land in an environment that is rich in moisture and nutrition, they germinate and grow. FGO assumes that spores will land in a random position at the beginning of the optimization process to promote the exploratory operator. As the optimization process is exceeded, this random position is transformed into a position between the best-so-far solution and a random position, promoting the exploitation operator while preventing premature convergence. FGO is evaluated against four well-known Congress on Evolutionary Computation (CEC) benchmarks (CEC2020, CEC2017, CEC2014, and CEC2022) and eleven engineering design problems. In addition, it is compared with fifteen recently proposed algorithms and eleven highly-performing algorithms, such as L-SHADE, LSHADE-cnEpSin, AL-SHADE, mantis search algorithm (MSA), IMODE, AGSK, SOMA_T3A, HyDE-DF, modified LSHADE-SPACMA, SHADE, and LSHADE-SPACMA, to demonstrate its superiority. According to the experimental results, FGO outperforms or is competitive with all of the compared algorithms for the majority of the test functions, implying that it is a high-performing optimizer and a powerful alternative technique for dealing with complex optimization problems. The FGO source code is available on this link</div><div><u>https://drive.mathworks.com/sharing/7b881d79-c7cb-4b64-bdfa-99ab7f57d984</u></div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"437 ","pages":"Article 117825"},"PeriodicalIF":6.9,"publicationDate":"2025-02-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143372897","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}