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

1 aHess, Martin1 aQuaini, Annalisa1 aRozza, Gianluigi uhttps://arxiv.org/abs/1901.0370801466nas a2200121 4500008004100000245010700041210006900148520102800217100002001245700002101265700002101286856003701307 2020 eng d00aA Reduced Order technique to study bifurcating phenomena: application to the Gross-Pitaevskii equation0 aReduced Order technique to study bifurcating phenomena applicati3 aWe 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. journal = SIAM Journal on Scientific Computing

1 aPichi, Federico1 aQuaini, Annalisa1 aRozza, Gianluigi uhttps://arxiv.org/abs/1907.0708201763nas a2200145 4500008004100000245010400041210006900145300001000214490000800224520128300232100002301515700002101538700002101559856003701580 2019 eng d00aA Finite Volume approximation of the Navier-Stokes equations with nonlinear filtering stabilization0 aFinite Volume approximation of the NavierStokes equations with n a27-450 v1873 aWe 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.

1 aGirfoglio, Michele1 aQuaini, Annalisa1 aRozza, Gianluigi uhttps://arxiv.org/abs/1901.0525102124nas a2200169 4500008004100000245008400041210006900125300001200194490000800206520160300214100001701817700002101834700002101855700002101876700002001897856003701917 2019 eng d00aA Localized Reduced-Order Modeling Approach for PDEs with Bifurcating Solutions0 aLocalized ReducedOrder Modeling Approach for PDEs with Bifurcati a379-4030 v3513 aReduced-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.

1 aHess, Martin1 aAlla, Alessandro1 aQuaini, Annalisa1 aRozza, Gianluigi1 aGunzburger, Max uhttps://arxiv.org/abs/1807.08851