This contribution describes the implementation of a data-driven shape optimization pipeline in a naval architecture application. We adopt reduced order models in order to improve the efficiency of the overall optimization, keeping a modular and equation-free nature to target the industrial demand. We applied the above mentioned pipeline to a realistic cruise ship in order to reduce the total drag. We begin by defining the design space, generated by deforming an initial shape in a parametric way using free form deformation. The evaluation of the performance of each new hull is determined by simulating the flux via finite volume discretization of a two-phase (water and air) fluid. Since the fluid dynamics model can result very expensive—especially dealing with complex industrial geometries—we propose also a dynamic mode decomposition enhancement to reduce the computational cost of a single numerical simulation. The real-time computation is finally achieved by means of proper orthogonal decomposition with Gaussian process regression technique. Thanks to the quick approximation, a genetic optimization algorithm becomes feasible to converge towards the optimal shape.

%B Bolletino dell Unione Matematica Italiana %V 14 %P 211-230 %G eng %R 10.1007/s40574-020-00263-4 %0 Unpublished Work %D 2021 %T AN EXTENDED PHYSICS INFORMED NEURAL NETWORK FOR PRELIMINARY ANALYSIS OF PARAMETRIC OPTIMAL CONTROL PROBLEMS %A Nicola Demo %A Maria Strazzullo %A Gianluigi Rozza %G eng %0 Journal Article %J Journal of Marine Science and Engineering %D 2021 %T Hull Shape Design Optimization with Parameter Space and Model Reductions, and Self-Learning Mesh Morphing %A Nicola Demo %A Marco Tezzele %A Andrea Mola %A Gianluigi Rozza %XIn the field of parametric partial differential equations, shape optimization represents a challenging problem due to the required computational resources. In this contribution, a data-driven framework involving multiple reduction techniques is proposed to reduce such computational burden. Proper orthogonal decomposition (POD) and active subspace genetic algorithm (ASGA) are applied for a dimensional reduction of the original (high fidelity) model and for an efficient genetic optimization based on active subspace property. The parameterization of the shape is applied directly to the computational mesh, propagating the generic deformation map applied to the surface (of the object to optimize) to the mesh nodes using a radial basis function (RBF) interpolation. Thus, topology and quality of the original mesh are preserved, enabling application of POD-based reduced order modeling techniques, and avoiding the necessity of additional meshing steps. Model order reduction is performed coupling POD and Gaussian process regression (GPR) in a data-driven fashion. The framework is validated on a benchmark ship.

%B Journal of Marine Science and Engineering %V 9 %P 185 %G eng %U https://www.mdpi.com/2077-1312/9/2/185 %R 10.3390/jmse9020185 %0 Generic %D 2021 %T The Neural Network shifted-Proper Orthogonal Decomposition: a Machine Learning Approach for Non-linear Reduction of Hyperbolic Equations %A Davide Papapicco %A Nicola Demo %A Michele Girfoglio %A Giovanni Stabile %A Gianluigi Rozza %G eng %0 Journal Article %J Software Impacts %D 2021 %T PyGeM: Python Geometrical Morphing %A Marco Tezzele %A Nicola Demo %A Andrea Mola %A Gianluigi Rozza %K Free form deformation %K Geometrical morphing %K Inverse distance weighting %K Python %K Radial basis functions interpolation %X PyGeM is an open source Python package which allows to easily parametrize and deform 3D object described by CAD files or 3D meshes. It implements several morphing techniques such as free form deformation, radial basis function interpolation, and inverse distance weighting. Due to its versatility in dealing with different file formats it is particularly suited for researchers and practitioners both in academia and in industry interested in computational engineering simulations and optimization studies. %B Software Impacts %V 7 %P 100047 %G eng %R 10.1016/j.simpa.2020.100047 %0 Journal Article %J SIAM Journal on Scientific Computing %D 2021 %T A supervised learning approach involving active subspaces for an efficient genetic algorithm in high-dimensional optimization problems %A Nicola Demo %A Marco Tezzele %A Gianluigi Rozza %XIn this work, we present an extension of the genetic algorithm (GA) which exploits the active subspaces (AS) property to evolve the individuals on a lower dimensional space. In many cases, GA requires in fact more function evaluations than others optimization method to converge to the optimum. Thus, complex and high-dimensional functions may result intractable with the standard algorithm. To address this issue, we propose to linearly map the input parameter space of the original function onto its AS before the evolution, performing the mutation and mate processes in a lower dimensional space. In this contribution, we describe the novel method called ASGA, presenting differences and similarities with the standard GA method. We test the proposed method over n-dimensional benchmark functions – Rosenbrock, Ackley, Bohachevsky, Rastrigin, Schaffer N. 7, and Zakharov – and finally we apply it to an aeronautical shape optimization problem.

%B SIAM Journal on Scientific Computing %V 43 %G eng %U https://arxiv.org/abs/2006.07282 %N 3 %& B831 %R https://doi.org/10.1137/20M1345219 %0 Conference Paper %B Proceedings of the 6th European Conference on Computational Mechanics: Solids, Structures and Coupled Problems, ECCM 2018 and 7th European Conference on Computational Fluid Dynamics, ECFD 2018 %D 2020 %T Advances in reduced order methods for parametric industrial problems in computational fluid dynamics %A Gianluigi Rozza %A M.H. Malik %A Nicola Demo %A Marco Tezzele %A Michele Girfoglio %A Giovanni Stabile %A Andrea Mola %XReduced order modeling has gained considerable attention in recent decades owing to the advantages offered in reduced computational times and multiple solutions for parametric problems. The focus of this manuscript is the application of model order reduction techniques in various engineering and scientific applications including but not limited to mechanical, naval and aeronautical engineering. The focus here is kept limited to computational fluid mechanics and related applications. The advances in the reduced order modeling with proper orthogonal decomposition and reduced basis method are presented as well as a brief discussion of dynamic mode decomposition and also some present advances in the parameter space reduction. Here, an overview of the challenges faced and possible solutions are presented with examples from various problems.

%B Proceedings of the 6th European Conference on Computational Mechanics: Solids, Structures and Coupled Problems, ECCM 2018 and 7th European Conference on Computational Fluid Dynamics, ECFD 2018 %G eng %U https://www.scopus.com/inward/record.uri?eid=2-s2.0-85075395686&partnerID=40&md5=fb0b1a3cfdfd35a104db9921bc9be675 %0 Journal Article %J Advanced Modeling and Simulation in Engineering Sciences %D 2020 %T Enhancing CFD predictions in shape design problems by model and parameter space reduction %A Marco Tezzele %A Nicola Demo %A Giovanni Stabile %A Andrea Mola %A Gianluigi Rozza %XIn this work we present an advanced computational pipeline for the approximation and prediction of the lift coefficient of a parametrized airfoil profile. The non-intrusive reduced order method is based on dynamic mode decomposition (DMD) and it is coupled with dynamic active subspaces (DyAS) to enhance the future state prediction of the target function and reduce the parameter space dimensionality. The pipeline is based on high-fidelity simulations carried out by the application of finite volume method for turbulent flows, and automatic mesh morphing through radial basis functions interpolation technique. The proposed pipeline is able to save 1/3 of the overall computational resources thanks to the application of DMD. Moreover exploiting DyAS and performing the regression on a lower dimensional space results in the reduction of the relative error in the approximation of the time-varying lift coefficient by a factor 2 with respect to using only the DMD.

%B Advanced Modeling and Simulation in Engineering Sciences %V 7 %G eng %U https://arxiv.org/abs/2001.05237 %N 40 %R https://doi.org/10.1186/s40323-020-00177-y %0 Journal Article %J Lecture Notes in Computational Science and Engineering %D 2020 %T Reduced order isogeometric analysis approach for pdes in parametrized domains %A Fabrizio Garotta %A Nicola Demo %A Marco Tezzele %A Massimo Carraturo %A Alessandro Reali %A Gianluigi Rozza %XIn this contribution, we coupled the isogeometric analysis to a reduced order modelling technique in order to provide a computationally efficient solution in parametric domains. In details, we adopt the free-form deformation method to obtain the parametric formulation of the domain and proper orthogonal decomposition with interpolation for the computational reduction of the model. This technique provides a real-time solution for any parameter by combining several solutions, in this case computed using isogeometric analysis on different geometrical configurations of the domain, properly mapped into a reference configuration. We underline that this reduced order model requires only the full-order solutions, making this approach non-intrusive. We present in this work the results of the application of this methodology to a heat conduction problem inside a deformable collector pipe.

%B Lecture Notes in Computational Science and Engineering %V 137 %P 153-170 %G eng %U https://www.scopus.com/inward/record.uri?eid=2-s2.0-85089615035&doi=10.1007%2f978-3-030-48721-8_7&partnerID=40&md5=7b15836ae65fa28dcfe8733788d7730c %R 10.1007/978-3-030-48721-8_7 %0 Conference Paper %B 8th International Conference on Computational Methods in Marine Engineering, MARINE 2019 %D 2019 %T A complete data-driven framework for the efficient solution of parametric shape design and optimisation in naval engineering problems %A Nicola Demo %A Marco Tezzele %A Andrea Mola %A Gianluigi Rozza %XIn the reduced order modeling (ROM) framework, the solution of a parametric partial differential equation is approximated by combining the high-fidelity solutions of the problem at hand for several properly chosen configurations. Examples of the ROM application, in the naval field, can be found in [31, 24]. Mandatory ingredient for the ROM methods is the relation between the high-fidelity solutions and the parameters. Dealing with geometrical parameters, especially in the industrial context, this relation may be unknown and not trivial (simulations over hand morphed geometries) or very complex (high number of parameters or many nested morphing techniques). To overcome these scenarios, we propose in this contribution an efficient and complete data-driven framework involving ROM techniques for shape design and optimization, extending the pipeline presented in [7]. By applying the singular value decomposition (SVD) to the points coordinates defining the hull geometry — assuming the topology is inaltered by the deformation —, we are able to compute the optimal space which the deformed geometries belong to, hence using the modal coefficients as the new parameters we can reconstruct the parametric formulation of the domain. Finally the output of interest is approximated using the proper orthogonal decomposition with interpolation technique. To conclude, we apply this framework to a naval shape design problem where the bulbous bow is morphed to reduce the total resistance of the ship advancing in calm water.

%B 8th International Conference on Computational Methods in Marine Engineering, MARINE 2019 %G eng %U https://www.scopus.com/inward/record.uri?eid=2-s2.0-85075342565&partnerID=40&md5=d76b8a1290053e7a84fb8801c0e6bb3d %0 Conference Paper %B VIII International Conference on Computational Methods in Marine Engineering %D 2019 %T A complete data-driven framework for the efficient solution of parametric shape design and optimisation in naval engineering problems %A Nicola Demo %A Marco Tezzele %A Andrea Mola %A Gianluigi Rozza %XIn the reduced order modeling (ROM) framework, the solution of a parametric partial differential equation is approximated by combining the high-fidelity solutions of the problem at hand for several properly chosen configurations. Examples of the ROM application, in the naval field, can be found in [31, 24]. Mandatory ingredient for the ROM methods is the relation between the high-fidelity solutions and the parameters. Dealing with geometrical parameters, especially in the industrial context, this relation may be unknown and not trivial (simulations over hand morphed geometries) or very complex (high number of parameters or many nested morphing techniques). To overcome these scenarios, we propose in this contribution an efficient and complete data-driven framework involving ROM techniques for shape design and optimization, extending the pipeline presented in [7]. By applying the singular value decomposition (SVD) to the points coordinates defining the hull geometry –- assuming the topology is inaltered by the deformation –-, we are able to compute the optimal space which the deformed geometries belong to, hence using the modal coefficients as the new parameters we can reconstruct the parametric formulation of the domain. Finally the output of interest is approximated using the proper orthogonal decomposition with interpolation technique. To conclude, we apply this framework to a naval shape design problem where the bulbous bow is morphed to reduce the total resistance of the ship advancing in calm water.

%B VIII International Conference on Computational Methods in Marine Engineering %G eng %U https://arxiv.org/abs/1905.05982 %0 Journal Article %J Comptes Rendus - Mecanique %D 2019 %T A non-intrusive approach for the reconstruction of POD modal coefficients through active subspaces %A Nicola Demo %A Marco Tezzele %A Gianluigi Rozza %XReduced order modeling (ROM) provides an efficient framework to compute solutions of parametric problems. Basically, it exploits a set of precomputed high-fidelity solutions—computed for properly chosen parameters, using a full-order model—in order to find the low dimensional space that contains the solution manifold. Using this space, an approximation of the numerical solution for new parameters can be computed in real-time response scenario, thanks to the reduced dimensionality of the problem. In a ROM framework, the most expensive part from the computational viewpoint is the calculation of the numerical solutions using the full-order model. Of course, the number of collected solutions is strictly related to the accuracy of the reduced order model. In this work, we aim at increasing the precision of the model also for few input solutions by coupling the proper orthogonal decomposition with interpolation (PODI)—a data-driven reduced order method—with the active subspace (AS) property, an emerging tool for reduction in parameter space. The enhanced ROM results in a reduced number of input solutions to reach the desired accuracy. In this contribution, we present the numerical results obtained by applying this method to a structural problem and in a fluid dynamics one.

%B Comptes Rendus - Mecanique %V 347 %P 873-881 %G eng %U https://www.scopus.com/inward/record.uri?eid=2-s2.0-85075379471&doi=10.1016%2fj.crme.2019.11.012&partnerID=40&md5=dcb27af39dc14dc8c3a4a5f681f7d84b %R https://doi.org/10.1016/j.crme.2019.11.012 %0 Conference Paper %B 8th International Conference on Computational Methods in Marine Engineering, MARINE 2019 %D 2019 %T Shape optimization through proper orthogonal decomposition with interpolation and dynamic mode decomposition enhanced by active subspaces %A Marco Tezzele %A Nicola Demo %A Gianluigi Rozza %XWe propose a numerical pipeline for shape optimization in naval engineering involving two different non-intrusive reduced order method (ROM) techniques. Such methods are proper orthogonal decomposition with interpolation (PODI) and dynamic mode decomposition (DMD). The ROM proposed will be enhanced by active subspaces (AS) as a pre-processing tool that reduce the parameter space dimension and suggest better sampling of the input space. We will focus on geometrical parameters describing the perturbation of a reference bulbous bow through the free form deformation (FFD) technique. The ROM are based on a finite volume method (FV) to simulate the multi-phase incompressible flow around the deformed hulls. In previous works we studied the reduction of the parameter space in naval engineering through AS [38, 10] focusing on different parts of the hull. PODI and DMD have been employed for the study of fast and reliable shape optimization cycles on a bulbous bow in [9]. The novelty of this work is the simultaneous reduction of both the input parameter space and the output fields of interest. In particular AS will be trained computing the total drag resistance of a hull advancing in calm water and its gradients with respect to the input parameters. DMD will improve the performance of each simulation of the campaign using only few snapshots of the solution fields in order to predict the regime state of the system. Finally PODI will interpolate the coefficients of the POD decomposition of the output fields for a fast approximation of all the fields at new untried parameters given by the optimization algorithm. This will result in a non-intrusive data-driven numerical optimization pipeline completely independent with respect to the full order solver used and it can be easily incorporated into existing numerical pipelines, from the reference CAD to the optimal shape.

%B 8th International Conference on Computational Methods in Marine Engineering, MARINE 2019 %G eng %U https://www.scopus.com/inward/record.uri?eid=2-s2.0-85075390244&partnerID=40&md5=3e1f2e9a2539d34594caff13766c94b8 %0 Conference Proceedings %B The 28th International Ocean and Polar Engineering Conference %D 2018 %T An efficient shape parametrisation by free-form deformation enhanced by active subspace for hull hydrodynamic ship design problems in open source environment %A Nicola Demo %A Marco Tezzele %A Andrea Mola %A Gianluigi Rozza %K Active subspaces %K Boundary element method %K Dynamic mode decomposition %K Fluid structure interaction %K Free form deformation %K Fully nonlinear potential %K Numerical towing tank %X In this contribution, we present the results of the application of a parameter space reduction methodology based on active subspaces to the hull hydrodynamic design problem. Several parametric deformations of an initial hull shape are considered to assess the influence of the shape parameters considered on the hull total drag. The hull resistance is typically computed by means of numerical simulations of the hydrodynamic flow past the ship. Given the high number of parameters involved - which might result in a high number of time consuming hydrodynamic simulations - assessing whether the parameters space can be reduced would lead to considerable computational cost reduction. Thus, the main idea of this work is to employ the active subspaces to identify possible lower dimensional structures in the parameter space, or to verify the parameter distribution in the position of the control points. To this end, a fully automated procedure has been implemented to produce several small shape perturbations of an original hull CAD geometry which are then used to carry out high-fidelity flow simulations and collect data for the active subspaces analysis. To achieve full automation of the open source pipeline described, both the free form deformation methodology employed for the hull perturbations and the solver based on unsteady potential flow theory, with fully nonlinear free surface treatment, are directly interfaced with CAD data structures and operate using IGES vendor-neutral file formats as input files. The computational cost of the fluid dynamic simulations is further reduced through the application of dynamic mode decomposition to reconstruct the steady state total drag value given only few initial snapshots of the simulation. The active subspaces analysis is here applied to the geometry of the DTMB-5415 naval combatant hull, which is which is a common benchmark in ship hydrodynamics simulations. %B The 28th International Ocean and Polar Engineering Conference %I International Society of Offshore and Polar Engineers %C Sapporo, Japan %G eng %U https://www.onepetro.org/conference-paper/ISOPE-I-18-481 %0 Journal Article %J The Journal of Open Source Software %D 2018 %T EZyRB: Easy Reduced Basis method %A Nicola Demo %A Marco Tezzele %A Gianluigi Rozza %B The Journal of Open Source Software %V 3 %P 661 %G eng %U https://joss.theoj.org/papers/10.21105/joss.00661 %R 10.21105/joss.00661 %0 Conference Paper %B Technology and Science for the Ships of the Future: Proceedings of NAV 2018: 19th International Conference on Ship & Maritime Research %D 2018 %T Model Order Reduction by means of Active Subspaces and Dynamic Mode Decomposition for Parametric Hull Shape Design Hydrodynamics %A Marco Tezzele %A Nicola Demo %A Mahmoud Gadalla %A Andrea Mola %A Gianluigi Rozza %X We present the results of the application of a parameter space reduction methodology based on active subspaces (AS) to the hull hydrodynamic design problem. Several parametric deformations of an initial hull shape are considered to assess the influence of the shape parameters on the hull wave resistance. Such problem is relevant at the preliminary stages of the ship design, when several flow simulations are carried out by the engineers to establish a certain sensibility with respect to the parameters, which might result in a high number of time consuming hydrodynamic simulations. The main idea of this work is to employ the AS to identify possible lower dimensional structures in the parameter space. The complete pipeline involves the use of free form deformation to parametrize and deform the hull shape, the full order solver based on unsteady potential flow theory with fully nonlinear free surface treatment directly interfaced with CAD, the use of dynamic mode decomposition to reconstruct the final steady state given only few snapshots of the simulation, and the reduction of the parameter space by AS, and shared subspace. Response surface method is used to minimize the total drag. %B Technology and Science for the Ships of the Future: Proceedings of NAV 2018: 19th International Conference on Ship & Maritime Research %I IOS Press %C Trieste, Italy %G eng %U http://ebooks.iospress.nl/publication/49270 %R 10.3233/978-1-61499-870-9-569 %0 Journal Article %J The Journal of Open Source Software %D 2018 %T PyDMD: Python Dynamic Mode Decomposition %A Nicola Demo %A Marco Tezzele %A Gianluigi Rozza %B The Journal of Open Source Software %V 3 %P 530 %G eng %U https://joss.theoj.org/papers/734e4326edd5062c6e8ee98d03df9e1d %R 10.21105/joss.00530 %0 Conference Paper %B Technology and Science for the Ships of the Future: Proceedings of NAV 2018: 19th International Conference on Ship & Maritime Research %D 2018 %T Shape Optimization by means of Proper Orthogonal Decomposition and Dynamic Mode Decomposition %A Nicola Demo %A Marco Tezzele %A Gianluca Gustin %A Gianpiero Lavini %A Gianluigi Rozza %X Shape optimization is a challenging task in many engineering fields, since the numerical solutions of parametric system may be computationally expensive. This work presents a novel optimization procedure based on reduced order modeling, applied to a naval hull design problem. The advantage introduced by this method is that the solution for a specific parameter can be expressed as the combination of few numerical solutions computed at properly chosen parametric points. The reduced model is built using the proper orthogonal decomposition with interpolation (PODI) method. We use the free form deformation (FFD) for an automated perturbation of the shape, and the finite volume method to simulate the multiphase incompressible flow around the deformed hulls. Further computational reduction is done by the dynamic mode decomposition (DMD) technique: from few high dimensional snapshots, the system evolution is reconstructed and the final state of the simulation is faithfully approximated. Finally the global optimization algorithm iterates over the reduced space: the approximated drag and lift coefficients are projected to the hull surface, hence the resistance is evaluated for the new hulls until the convergence to the optimal shape is achieved. We will present the results obtained applying the described procedure to a typical Fincantieri cruise ship. %B Technology and Science for the Ships of the Future: Proceedings of NAV 2018: 19th International Conference on Ship & Maritime Research %I IOS Press %C Trieste, Italy %G eng %U http://ebooks.iospress.nl/publication/49229 %& 212 %R 10.3233/978-1-61499-870-9-212 %0 Conference Paper %B International Conference on Parallel Processing and Applied Mathematics %D 2015 %T Experience on vectorizing lattice Boltzmann kernels for multi-and many-core architectures %A Calore, Enrico %A Nicola Demo %A Schifano, Sebastiano Fabio %A Tripiccione, Raffaele %B International Conference on Parallel Processing and Applied Mathematics %I Springer %P 53–62 %G eng