In 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.

1 aDemo, Nicola1 aTezzele, Marco1 aRozza, Gianluigi uhttps://arxiv.org/abs/2006.0728200473nas a2200121 4500008004100000245004900041210004900090300001000139490000700149100002000156700002100176856015400197 2020 eng d00aSpecial Issue on Reduced Order Models in CFD0 aSpecial Issue on Reduced Order Models in CFD a91-920 v341 aPerotto, Simona1 aRozza, Gianluigi uhttps://www.scopus.com/inward/record.uri?eid=2-s2.0-85084258805&doi=10.1080%2f10618562.2020.1756497&partnerID=40&md5=d9316aad9ba95f244e07379318ebbcba01342nas a2200145 4500008004100000245010000041210007100141300001200212490000800224520076800232100002101000700002101021700002101042856013301063 2020 eng d00aA spectral element reduced basis method for navier–stokes equations with geometric variations0 aspectral element reduced basis method for navier–stokes equation a561-5710 v1343 aWe consider the Navier-Stokes equations in a channel with a narrowing of varying height. The model is discretized with high-order spectral element ansatz functions, resulting in 6372 degrees of freedom. The steady-state snapshot solutions define a reduced order space through a standard POD procedure. The reduced order space allows to accurately and efficiently evaluate the steady-state solutions for different geometries. In particular, we detail different aspects of implementing the reduced order model in combination with a spectral element discretization. It is shown that an expansion in element-wise local degrees of freedom can be combined with a reduced order modelling approach to enhance computational times in parametric many-query scenarios.

1 aHess, Martin, W.1 aQuaini, Annalisa1 aRozza, Gianluigi uhttps://www.math.sissa.it/publication/spectral-element-reduced-basis-method-navier%E2%80%93stokes-equations-geometric-variations01963nas a2200145 4500008004100000245009800041210006900139300001400208490000700222520138100229100001701610700001701627700002101644856015201665 2020 eng d00aStabilized reduced basis methods for parametrized steady Stokes and Navier–Stokes equations0 aStabilized reduced basis methods for parametrized steady Stokes a2399-24160 v803 aIt is well known in the Reduced Basis approximation of saddle point problems that the Galerkin projection on the reduced space does not guarantee the inf–sup approximation stability even if a stable high fidelity method was used to generate snapshots. For problems in computational fluid dynamics, the lack of inf–sup stability is reflected by the inability to accurately approximate the pressure field. In this context, inf–sup stability is usually recovered through the enrichment of the velocity space with suitable supremizer functions. The main goal of this work is to propose an alternative approach, which relies on the residual based stabilization techniques customarily employed in the Finite Element literature, such as Brezzi–Pitkaranta, Franca–Hughes, streamline upwind Petrov–Galerkin, Galerkin Least Square. In the spirit of offline–online reduced basis computational splitting, two such options are proposed, namely offline-only stabilization and offline–online stabilization. These approaches are then compared to (and combined with) the state of the art supremizer enrichment approach. Numerical results are discussed, highlighting that the proposed methodology allows to obtain smaller reduced basis spaces (i.e., neglecting supremizer enrichment) for which a modified inf–sup stability is still preserved at the reduced order level.

1 aAli, Shafqat1 aBallarin, F.1 aRozza, Gianluigi uhttps://www.scopus.com/inward/record.uri?eid=2-s2.0-85083340115&doi=10.1016%2fj.camwa.2020.03.019&partnerID=40&md5=7ace96eee080701acb04d8155008dd7d02446nas a2200121 4500008004100000245014200041210006900183520189700252100001902149700001702168700002102185856011802206 2019 eng d00aShape optimization through proper orthogonal decomposition with interpolation and dynamic mode decomposition enhanced by active subspaces0 aShape optimization through proper orthogonal decomposition with 3 aWe 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.

1 aTezzele, Marco1 aDemo, Nicola1 aRozza, Gianluigi uhttps://www.scopus.com/inward/record.uri?eid=2-s2.0-85075390244&partnerID=40&md5=3e1f2e9a2539d34594caff13766c94b801397nas a2200193 4500008004100000245006200041210006000103260003800163490000800201520080600209100002101015700002101036700002501057700001801082700001601100700002801116700002201144856003701166 2019 eng d00aA Spectral Element Reduced Basis Method in Parametric CFD0 aSpectral Element Reduced Basis Method in Parametric CFD bSpringer International Publishing0 v1263 aWe consider the Navier-Stokes equations in a channel with varying Reynolds numbers. The model is discretized with high-order spectral element ansatz functions, resulting in 14 259 degrees of freedom. The steady-state snapshot solu- tions define a reduced order space, which allows to accurately evaluate the steady- state solutions for varying Reynolds number with a reduced order model within a fixed-point iteration. In particular, we compare different aspects of implementing the reduced order model with respect to the use of a spectral element discretization. It is shown, how a multilevel static condensation in the pressure and velocity boundary degrees of freedom can be combined with a reduced order modelling approach to enhance computational times in parametric many-query scenarios.

1 aHess, Martin, W.1 aRozza, Gianluigi1 aRadu, Florin, Adrian1 aKumar, Kundan1 aBerre, Inga1 aNordbotten, Jan, Martin1 aPop, Iuliu, Sorin uhttps://arxiv.org/abs/1712.0643201451nas a2200133 4500008004100000245006200041210006000103300001200163490000800175520093900183100002101122700002101143856015301164 2019 eng d00aA spectral element reduced basis method in parametric CFD0 aspectral element reduced basis method in parametric CFD a693-7010 v1263 aWe consider the Navier-Stokes equations in a channel with varying Reynolds numbers. The model is discretized with high-order spectral element ansatz functions, resulting in 14,259 degrees of freedom. The steady-state snapshot solutions define a reduced order space, which allows to accurately evaluate the steady-state solutions for varying Reynolds number with a reduced order model within a fixed-point iteration. In particular, we compare different aspects of implementing the reduced order model with respect to the use of a spectral element discretization. It is shown, how a multilevel static condensation (Karniadakis and Sherwin, Spectral/hp element methods for computational fluid dynamics, 2nd edn. Oxford University Press, Oxford, 2005) in the pressure and velocity boundary degrees of freedom can be combined with a reduced order modelling approach to enhance computational times in parametric many-query scenarios.

1 aHess, Martin, W.1 aRozza, Gianluigi uhttps://www.scopus.com/inward/record.uri?eid=2-s2.0-85060005503&doi=10.1007%2f978-3-319-96415-7_64&partnerID=40&md5=d1a900db8ddb92cd818d797ec212a4c601912nas a2200157 4500008004100000245009800041210006900139260003000208520136800238100001701606700001901623700002101642700002201663700002101685856004801706 2018 eng d00aShape Optimization by means of Proper Orthogonal Decomposition and Dynamic Mode Decomposition0 aShape Optimization by means of Proper Orthogonal Decomposition a aTrieste, ItalybIOS Press3 aShape 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.1 aDemo, Nicola1 aTezzele, Marco1 aGustin, Gianluca1 aLavini, Gianpiero1 aRozza, Gianluigi uhttp://ebooks.iospress.nl/publication/4922901774nas a2200169 4500008004100000245006400041210006100105520122900166100001801395700001801413700001401431700001701445700001701462700001901479700002101498856008501519 2018 eng d00aSRTP 2.0 - The evolution of the safe return to port concept0 aSRTP 20 The evolution of the safe return to port concept3 aIn 2010 IMO (International Maritime Organisation) introduced new rules in SOLAS with the aim of intrinsically increase the safety of passenger ships. This requirement is achieved by providing safe areas for passengers and essential services for allowing ship to Safely Return to Port (SRtP). The entry into force of these rules has changed the way to design passenger ships. In this respect big effort in the research has been done by industry to address design issues related to the impact on failure analysis of the complex interactions among systems. Today the research activity is working to bring operational matters in the design stage. This change of research focus was necessary because human factor and the way to operate the ship itself after a casualty on board may have a big impact in the design of the ship/systems. Also the management of the passengers after a casualty is becoming a major topic for safety. This paper presents the state of the art of Italian knowledge in the field of system engineering applied to passenger ship address to safety improvement and design reliability. An overview of present tools and methodologies will be offered together with future focuses in the research activity.

1 aCangelosi, D.1 aBonvicini, A.1 aNardo, M.1 aMola, Andrea1 aMarchese, A.1 aTezzele, Marco1 aRozza, Gianluigi uhttps://www.math.sissa.it/publication/srtp-20-evolution-safe-return-port-concept01436nas a2200145 4500008004100000245011100041210006900152300001400221490000600235520085400241100001401095700001701109700002101126856014301147 2018 eng d00aStabilized weighted reduced basis methods for parametrized advection dominated problems with random inputs0 aStabilized weighted reduced basis methods for parametrized advec a1475-15020 v63 aIn this work, we propose viable and eficient strategies for stabilized parametrized advection dominated problems, with random inputs. In particular, we investigate the combination of the wRB (weighted reduced basis) method for stochastic parametrized problems with the stabilized RB (reduced basis) method, which is the integration of classical stabilization methods (streamline/upwind Petrov-Galerkin (SUPG) in our case) in the ofine-online structure of the RB method. Moreover, we introduce a reduction method that selectively enables online stabilization; this leads to a sensible reduction of computational costs, while keeping a very good accuracy with respect to high-fdelity solutions. We present numerical test cases to assess the performance of the proposed methods in steady and unsteady problems related to heat transfer phenomena.

1 aTorlo, D.1 aBallarin, F.1 aRozza, Gianluigi uhttps://www.scopus.com/inward/record.uri?eid=2-s2.0-85058246502&doi=10.1137%2f17M1163517&partnerID=40&md5=6c54e2f0eb727cb85060e988486b8ac800682nas a2200145 4500008004100000245009900041210006900140260001000209520018600219100001700405700002000422700002200442700002100464856005100485 2015 en d00aSupremizer stabilization of POD-Galerkin approximation of parametrized Navier-Stokes equations0 aSupremizer stabilization of PODGalerkin approximation of paramet bWiley3 aIn this work, we present a stable proper orthogonal decomposition–Galerkin approximation for parametrized steady incompressible Navier–Stokes equations with low Reynolds number.1 aBallarin, F.1 aManzoni, Andrea1 aQuarteroni, Alfio1 aRozza, Gianluigi uhttp://urania.sissa.it/xmlui/handle/1963/3470101619nas a2200145 4500008004100000245010700041210006900148260001300217520111600230100001701346700002001363700002101383700001801404856005101422 2014 en d00aShape Optimization by Free-Form Deformation: Existence Results and Numerical Solution for Stokes Flows0 aShape Optimization by FreeForm Deformation Existence Results and bSpringer3 aShape optimization problems governed by PDEs result from many applications in computational fluid dynamics. These problems usually entail very large computational costs and require also a suitable approach for representing and deforming efficiently the shape of the underlying geometry, as well as for computing the shape gradient of the cost functional to be minimized. Several approaches based on the displacement of a set of control points have been developed in the last decades, such as the so-called free-form deformations. In this paper we present a new theoretical result which allows to recast free-form deformations into the general class of perturbation of identity maps, and to guarantee the compactness of the set of admissible shapes. Moreover, we address both a general optimization framework based on the continuous shape gradient and a numerical procedure for solving efficiently three-dimensional optimal design problems. This framework is applied to the optimal design of immersed bodies in Stokes flows, for which we consider the numerical solution of a benchmark case study from literature.1 aBallarin, F.1 aManzoni, Andrea1 aRozza, Gianluigi1 aSalsa, Sandro uhttp://urania.sissa.it/xmlui/handle/1963/3469801201nas a2200133 4500008004100000245007800041210006900119300001100188490000800199520070800207100001900915700002100934856011200955 2014 eng d00aStabilized reduced basis method for parametrized advection-diffusion PDEs0 aStabilized reduced basis method for parametrized advectiondiffus a1–180 v2743 aIn this work, we propose viable and efficient strategies for the stabilization of the reduced basis approximation of an advection dominated problem. In particular, we investigate the combination of a classic stabilization method (SUPG) with the Offline-Online structure of the RB method. We explain why the stabilization is needed in both stages and we identify, analytically and numerically, which are the drawbacks of a stabilization performed only during the construction of the reduced basis (i.e. only in the Offline stage). We carry out numerical tests to assess the performances of the ``double'' stabilization both in steady and unsteady problems, also related to heat transfer phenomena.

1 aPacciarini, P.1 aRozza, Gianluigi uhttps://www.math.sissa.it/publication/stabilized-reduced-basis-method-parametrized-advection-diffusion-pdes01104nas a2200121 4500008004100000245016100041210006900202300001600271520058500287100001900872700002100891856007000912 2014 eng d00aStabilized reduced basis method for parametrized scalar advection-diffusion problems at higher Péclet number: Roles of the boundary layers and inner fronts0 aStabilized reduced basis method for parametrized scalar advectio a5614–56243 aAdvection-dominated problems, which arise in many engineering situations, often require a fast and reliable approximation of the solution given some parameters as inputs. In this work we want to investigate the coupling of the reduced basis method - which guarantees rapidity and reliability - with some classical stabilization techiques to deal with the advection-dominated condition. We provide a numerical extension of the results presented in [1], focusing in particular on problems with curved boundary layers and inner fronts whose direction depends on the parameter.

1 aPacciarini, P.1 aRozza, Gianluigi uhttps://infoscience.epfl.ch/record/203327/files/ECCOMAS_PP_GR.pdf01509nas a2200145 4500008004100000245009700041210006900138300001600207490000700223520095300230100001501183700002201198700002101220856012201241 2013 eng d00aStochastic optimal robin boundary control problems of advection-dominated elliptic equations0 aStochastic optimal robin boundary control problems of advectiond a2700–27220 v513 aIn this work we deal with a stochastic optimal Robin boundary control problem constrained by an advection-diffusion-reaction elliptic equation with advection-dominated term. We assume that the uncertainty comes from the advection field and consider a stochastic Robin boundary condition as control function. A stochastic saddle point system is formulated and proved to be equivalent to the first order optimality system for the optimal control problem, based on which we provide the existence and uniqueness of the optimal solution as well as some results on stochastic regularity with respect to the random variables. Stabilized finite element approximations in physical space and collocation approximations in stochastic space are applied to discretize the optimality system. A global error estimate in the product of physical space and stochastic space for the numerical approximation is derived. Illustrative numerical experiments are provided.1 aChen, Peng1 aQuarteroni, Alfio1 aRozza, Gianluigi uhttps://www.math.sissa.it/publication/stochastic-optimal-robin-boundary-control-problems-advection-dominated-elliptic01580nas a2200145 4500008004100000245008700041210006900128260001000197520093400207653011301141100001501254700002201269700002101291856012201312 2012 en d00aSimulation-based uncertainty quantification of human arterial network hemodynamics0 aSimulationbased uncertainty quantification of human arterial net bWiley3 aThis work aims at identifying and quantifying uncertainties from various sources in human cardiovascular\r\nsystem based on stochastic simulation of a one dimensional arterial network. A general analysis of\r\ndifferent uncertainties and probability characterization with log-normal distribution of these uncertainties\r\nis introduced. Deriving from a deterministic one dimensional fluid structure interaction model, we establish\r\nthe stochastic model as a coupled hyperbolic system incorporated with parametric uncertainties to describe\r\nthe blood flow and pressure wave propagation in the arterial network. By applying a stochastic collocation\r\nmethod with sparse grid technique, we study systemically the statistics and sensitivity of the solution with\r\nrespect to many different uncertainties in a relatively complete arterial network with potential physiological\r\nand pathological implications for the first time.10auncertainty quantification, mathematical modelling of the cardiovascular system, fluid-structure interaction1 aChen, Peng1 aQuarteroni, Alfio1 aRozza, Gianluigi uhttps://www.math.sissa.it/publication/simulation-based-uncertainty-quantification-human-arterial-network-hemodynamics