The majority of the most common physical phenomena can be described using partial differential equations (PDEs). However, they are very often characterized by strong nonlinearities. Such features lead to the coexistence of multiple solutions studied by the bifurcation theory. Unfortunately, in practical scenarios, one has to exploit numerical methods to compute the solutions of systems of PDEs, even if the classical techniques are usually able to compute only a single solution for any value of a parameter when more branches exist. In this work we implemented an elaborated deflated continuation method, that relies on the spectral element method (SEM) and on the reduced basis (RB) one, to efficiently compute bifurcation diagrams with more parameters and more bifurcation points. The deflated continuation method can be obtained combining the classical continuation method and the deflation one: the former is used to entirely track each known branch of the diagram, while the latter is exploited to discover the new ones. Finally, when more than one parameter is considered, the efficiency of the computation is ensured by the fact that the diagrams can be computed during the online phase while, during the offline one, one only has to compute one-dimensional diagrams. In this work, after a more detailed description of the method, we will show the results that can be obtained using it to compute a bifurcation diagram associated with a problem governed by the Navier-Stokes equations.

%B Advances in Computational Mathematics %G eng %U https://arxiv.org/abs/1912.06089 %0 Conference Paper %B QUIET Selected Contributions %D 2020 %T Non-Intrusive Polynomial Chaos Method Applied to Problems in Computational Fluid Dynamics with a Comparison to Proper Orthogonal Decomposition %A Saddam Hijazi %A Giovanni Stabile %A Andrea Mola %A Gianluigi Rozza %E van Brummelen, Harald %E Corsini, Alessandro %E Perotto, Simona %E Rozza, Gianluigi %XIn this work, Uncertainty Quantification (UQ) based on non-intrusive Polynomial Chaos Expansion (PCE) is applied to the CFD problem of the flow past an airfoil with parameterized angle of attack and inflow velocity. To limit the computational cost associated with each of the simulations required by the non-intrusive UQ algorithm used, we resort to a Reduced Order Model (ROM) based on Proper Orthogonal Decomposition (POD)-Galerkin approach. A first set of results is presented to characterize the accuracy of the POD-Galerkin ROM developed approach with respect to the Full Order Model (FOM) solver (OpenFOAM). A further analysis is then presented to assess how the UQ results are affected by substituting the FOM predictions with the surrogate ROM ones.

%B QUIET Selected Contributions %I Springer International Publishing %G eng %U https://arxiv.org/abs/1901.02285 %& Non-Intrusive Polynomial Chaos Method Applied to Problems in Computational Fluid Dynamics [...] %0 Journal Article %J International Journal of Computational Fluid Dynamics %D 2020 %T Reduced Basis Model Order Reduction for Navier-Stokes equations in domains with walls of varying curvature %A Hess, Martin %A Quaini, Annalisa %A Rozza, Gianluigi %XWe consider the Navier-Stokes equations in a channel with a narrowing and walls of varying curvature. By applying the empirical interpolation method to generate an affine parameter dependency, the offline-online procedure can be used to compute reduced order solutions for parameter variations. The reduced order space is computed from the steady-state snapshot solutions by a standard POD procedure. The model is discretised with high-order spectral element ansatz functions, resulting in 4752 degrees of freedom. The proposed reduced order model produces accurate approximations of steady-state solutions for a wide range of geometries and kinematic viscosity values. The application that motivated the present study is the onset of asymmetries (i.e., symmetry breaking bifurcation) in blood flow through a regurgitant mitral valve, depending on the Reynolds number and the valve shape. Through our computational study, we found that the critical Reynolds number for the symmetry breaking increases as the wall curvature increases.

%B International Journal of Computational Fluid Dynamics %V 34 %P 119-126 %G eng %U https://arxiv.org/abs/1901.03708 %R 10.1080/10618562.2019.1645328 %0 Journal Article %J SIAM Journal on Scientific Computing %D 2020 %T A Reduced Order technique to study bifurcating phenomena: application to the Gross-Pitaevskii equation %A Pichi, Federico %A Quaini, Annalisa %A Rozza, Gianluigi %XWe propose a computationally efficient framework to treat nonlinear partial differential equations having bifurcating solutions as one or more physical control parameters are varied. Our focus is on steady bifurcations. Plotting a bifurcation diagram entails computing multiple solutions of a parametrized, nonlinear problem, which can be extremely expensive in terms of computational time. In order to reduce these demanding computational costs, our approach combines a continuation technique and Newton's method with a Reduced Order Modeling (ROM) technique, suitably supplemented with a hyper-reduction method. To demonstrate the effectiveness of our ROM approach, we trace the steady solution branches of a nonlinear Schrödinger equation, called Gross-Pitaevskii equation, as one or two physical parameters are varied. In the two parameter study, we show that our approach is 60 times faster in constructing a bifurcation diagram than a standard Full Order Method.

%B SIAM Journal on Scientific Computing %G eng %U https://arxiv.org/abs/1907.07082 %R https://doi.org/10.1137/20M1313106 %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 Demo, Nicola %A Tezzele, Marco %A Mola, Andrea %A Rozza, Gianluigi %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 Conference Paper %B VIII International Conference on Computational Methods in Marine Engineering %D 2019 %T Efficient Reduction in Shape Parameter Space Dimension for Ship Propeller Blade Design %A Mola, Andrea %A Tezzele, Marco %A Gadalla, Mahmoud %A Valdenazzi, Federica %A Grassi, Davide %A Padovan, Roberta %A Rozza, Gianluigi %XIn this work, we present the results of a ship propeller design optimization campaign carried out in the framework of the research project PRELICA, funded by the Friuli Venezia Giulia regional government. The main idea of this work is to operate on a multidisciplinary level to identify propeller shapes that lead to reduced tip vortex-induced pressure and increased efficiency without altering the thrust. First, a specific tool for the bottom-up construction of parameterized propeller blade geometries has been developed. The algorithm proposed operates with a user defined number of arbitrary shaped or NACA airfoil sections, and employs arbitrary degree NURBS to represent the chord, pitch, skew and rake distribution as a function of the blade radial coordinate. The control points of such curves have been modified to generate, in a fully automated way, a family of blade geometries depending on as many as 20 shape parameters. Such geometries have then been used to carry out potential flow simulations with the Boundary Element Method based software PROCAL. Given the high number of parameters considered, such a preliminary stage allowed for a fast evaluation of the performance of several hundreds of shapes. In addition, the data obtained from the potential flow simulation allowed for the application of a parameter space reduction methodology based on active subspaces (AS) property, which suggested that the main propeller performance indices are, at a first but rather accurate approximation, only depending on a single parameter which is a linear combination of all the original geometric ones. AS analysis has also been used to carry out a constrained optimization exploiting response surface method in the reduced parameter space, and a sensitivity analysis based on such surrogate model. The few selected shapes were finally used to set up high fidelity RANS simulations and select an optimal shape.

%B VIII International Conference on Computational Methods in Marine Engineering %G eng %U https://arxiv.org/abs/1905.09815 %0 Journal Article %J Computers & Fluids %D 2019 %T A Finite Volume approximation of the Navier-Stokes equations with nonlinear filtering stabilization %A Girfoglio, Michele %A Quaini, Annalisa %A Rozza, Gianluigi %XWe consider a Leray model with a nonlinear differential low-pass filter for the simulation of incompressible fluid flow at moderately large Reynolds number (in the range of a few thousands) with under-refined meshes. For the implementation of the model, we adopt the three-step algorithm Evolve-Filter-Relax (EFR). The Leray model has been extensively applied within a Finite Element (FE) framework. Here, we propose to combine the EFR algorithm with a computationally efficient Finite Volume (FV) method. Our approach is validated against numerical data available in the literature for the 2D flow past a cylinder and against experimental measurements for the 3D fluid flow in an idealized medical device, as recommended by the U.S. Food and Drug Administration. We will show that for similar levels of mesh refinement FV and FE methods provide significantly different results. Through our numerical experiments, we are able to provide practical directions to tune the parameters involved in the model. Furthermore, we are able to investigate the impact of mesh features (element type, non-orthogonality, local refinement, and element aspect ratio) and the discretization method for the convective term on the agreement between numerical solutions and experimental data.

%B Computers & Fluids %V 187 %P 27-45 %G eng %U https://arxiv.org/abs/1901.05251 %R 10.1016/j.compfluid.2019.05.001 %0 Journal Article %J Computer Methods in Applied Mechanics and Engineering %D 2019 %T A Localized Reduced-Order Modeling Approach for PDEs with Bifurcating Solutions %A Hess, Martin %A Alla, Alessandro %A Quaini, Annalisa %A Rozza, Gianluigi %A Gunzburger, Max %XReduced-order modeling (ROM) commonly refers to the construction, based on a few solutions (referred to as snapshots) of an expensive discretized partial differential equation (PDE), and the subsequent application of low-dimensional discretizations of partial differential equations (PDEs) that can be used to more efficiently treat problems in control and optimization, uncertainty quantification, and other settings that require multiple approximate PDE solutions. In this work, a ROM is developed and tested for the treatment of nonlinear PDEs whose solutions bifurcate as input parameter values change. In such cases, the parameter domain can be subdivided into subregions, each of which corresponds to a different branch of solutions. Popular ROM approaches such as proper orthogonal decomposition (POD), results in a global low-dimensional basis that does no respect not take advantage of the often large differences in the PDE solutions corresponding to different subregions. Instead, in the new method, the k-means algorithm is used to cluster snapshots so that within cluster snapshots are similar to each other and are dissimilar to those in other clusters. This is followed by the construction of local POD bases, one for each cluster. The method also can detect which cluster a new parameter point belongs to, after which the local basis corresponding to that cluster is used to determine a ROM approximation. Numerical experiments show the effectiveness of the method both for problems for which bifurcation cause continuous and discontinuous changes in the solution of the PDE.

%B Computer Methods in Applied Mechanics and Engineering %V 351 %P 379-403 %G eng %U https://arxiv.org/abs/1807.08851 %R 10.1016/j.cma.2019.03.050 %0 Journal Article %J Computer Methods in Applied Mechanics and Engineering %D 2019 %T A reduced basis approach for PDEs on parametrized geometries based on the shifted boundary finite element method and application to a Stokes flow %A Karatzas, Efthymios N %A Stabile, Giovanni %A Nouveau, Leo %A Scovazzi, Guglielmo %A Rozza, Gianluigi %XWe propose a model order reduction technique integrating the Shifted Boundary Method (SBM) with a POD-Galerkin strategy. This approach allows to treat more complex parametrized domains in an efficient and straightforward way. The impact of the proposed approach is threefold. First, problems involving parametrizations of complex geometrical shapes and/or large domain deformations can be efficiently solved at full-order by means of the SBM, an unfitted boundary method that avoids remeshing and the tedious handling of cut cells by introducing an approximate surrogate boundary. Second, the computational effort is further reduced by the development of a reduced order model (ROM) technique based on a POD-Galerkin approach. Third, the SBM provides a smooth mapping from the true to the surrogate domain, and for this reason, the stability and performance of the reduced order basis are enhanced. This feature is the net result of the combination of the proposed ROM approach and the SBM. Similarly, the combination of the SBM with a projection-based ROM gives the great advantage of an easy and fast to implement algorithm considering geometrical parametrization with large deformations. The transformation of each geometry to a reference geometry (morphing) is in fact not required. These combined advantages will allow the solution of PDE problems more efficiently. We illustrate the performance of this approach on a number of two-dimensional Stokes flow problems.

%B Computer Methods in Applied Mechanics and Engineering %V 347 %P 568–587 %G eng %U https://arxiv.org/abs/1807.07790 %R 10.1016/j.cma.2018.12.040 %0 Journal Article %D 2019 %T Reduced basis approaches for parametrized bifurcation problems held by non-linear Von Kármán equations %A Pichi, Federico %A Rozza, Gianluigi %XThis work focuses on the computationally efficient detection of the buckling phenomena and bifurcation analysis of the parametric Von Kármán plate equations based on reduced order methods and spectral analysis. The computational complexity - due to the fourth order derivative terms, the non-linearity and the parameter dependence - provides an interesting benchmark to test the importance of the reduction strategies, during the construction of the bifurcation diagram by varying the parameter(s). To this end, together the state equations, we carry out also an analysis of the linearized eigenvalue problem, that allows us to better understand the physical behaviour near the bifurcation points, where we lose the uniqueness of solution. We test this automatic methodology also in the two parameter case, understanding the evolution of the first buckling mode. journal = Journal of Scientific Computing

%V 81 %P 112–135 %G eng %U https://arxiv.org/abs/1804.02014 %R 10.1007/s10915-019-01003-3 %0 Conference Paper %B VIII International Conference on Computational Methods in Marine Engineering %D 2019 %T Shape optimization through proper orthogonal decomposition with interpolation and dynamic mode decomposition enhanced by active subspaces %A Tezzele, Marco %A Demo, Nicola %A Rozza, Gianluigi %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 VIII International Conference on Computational Methods in Marine Engineering %G eng %U https://arxiv.org/abs/1905.05483 %0 Book Section %B Numerical Mathematics and Advanced Applications - ENUMATH 2017 %D 2019 %T A Spectral Element Reduced Basis Method in Parametric CFD %A Hess, Martin W. %A Rozza, Gianluigi %E Radu, Florin Adrian %E Kumar, Kundan %E Berre, Inga %E Nordbotten, Jan Martin %E Pop, Iuliu Sorin %XWe 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.

%B Numerical Mathematics and Advanced Applications - ENUMATH 2017 %I Springer International Publishing %V 126 %G eng %U https://arxiv.org/abs/1712.06432 %& A Spectral Element Reduced Basis Method in Parametric CFD %R 10.1007/978-3-319-96415-7_64 pages = 693–701 %0 Book Section %B Numerical Methods for PDEs %D 2018 %T Reduced Basis Approximation and A Posteriori Error Estimation: Applications to Elasticity Problems in Several Parametric Settings %A Huynh, D. B. P. %A Pichi, Federico %A Rozza, Gianluigi %B Numerical Methods for PDEs %V 15 %G eng %U https://link.springer.com/chapter/10.1007/978-3-319-94676-4_8 %R https://doi.org/10.1007/978-3-319-94676-4_8