We focus on reducing the computational costs associated with the hydrodynamic stability of solutions of the incompressible Navier–Stokes equations for a Newtonian and viscous fluid in contraction–expansion channels. In particular, we are interested in studying steady bifurcations, occurring when non-unique stable solutions appear as physical and/or geometric control parameters are varied. The formulation of the stability problem requires solving an eigenvalue problem for a partial differential operator. An alternative to this approach is the direct simulation of the flow to characterize the asymptotic behavior of the solution. Both approaches can be extremely expensive in terms of computational time. We propose to apply Reduced Order Modeling (ROM) techniques to reduce the demanding computational costs associated with the detection of a type of steady bifurcations in fluid dynamics. 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 regurgitant mitral valve orifice shape.

In this work we review a reduced basis method for the solution of uncertainty quantification problems. Based on the basic setting of an elliptic partial differential equation with random input, we introduce the key ingredients of the reduced basis method, including proper orthogonal decomposition and greedy algorithms for the construction of the reduced basis functions, a priori and a posteriori error estimates for the reduced basis approximations, as well as its computational advantages and weaknesses in comparison with a stochastic collocation method [I. Babuška, F. Nobile, and R. Tempone, *SIAM Rev.*, 52 (2010), pp. 317--355]. We demonstrate its computational efficiency and accuracy for a benchmark problem with parameters ranging from a few to a few hundred dimensions. Generalizations to more complex models and applications to uncertainty quantification problems in risk prediction, evaluation of statistical moments, Bayesian inversion, and optimal control under uncertainty are also presented to illustrate how to use the reduced basis method in practice. Further challenges, advancements, and research opportunities are outlined.

Read More: http://epubs.siam.org/doi/abs/10.1137/151004550

This monograph addresses the state of the art of reduced order methods for modeling and computational reduction of complex parametrized systems, governed by ordinary and/or partial differential equations, with a special emphasis on real time computing techniques and applications in computational mechanics, bioengineering and computer graphics.

Several topics are covered, including: design, optimization, and control theory in real-time with applications in engineering; data assimilation, geometry registration, and parameter estimation with special attention to real-time computing in biomedical engineering and computational physics; real-time visualization of physics-based simulations in computer science; the treatment of high-dimensional problems in state space, physical space, or parameter space; the interactions between different model reduction and dimensionality reduction approaches; the development of general error estimation frameworks which take into account both model and discretization effects.

This book is primarily addressed to computational scientists interested in computational reduction techniques for large scale differential problems.

10areduced order methods, MOR, ROM, POD, RB, greedy, CFD, Numerical Analysis1 aQuarteroni, Alfio1 aRozza, Gianluigi uhttps://www.math.sissa.it/publication/reduced-order-methods-modeling-and-computational-reduction01557nas a2200133 4500008004100000245009400041210006900135260001700204520109300221100001501314700002201329700002101351856005101372 2014 en d00aA weighted empirical interpolation method: A priori convergence analysis and applications0 aweighted empirical interpolation method A priori convergence ana bEDP Sciences3 aWe extend the classical empirical interpolation method [M. Barrault, Y. Maday, N.C. Nguyen and A.T. Patera, An empirical interpolation method: application to efficient reduced-basis discretization of partial differential equations. Compt. Rend. Math. Anal. Num. 339 (2004) 667-672] to a weighted empirical interpolation method in order to approximate nonlinear parametric functions with weighted parameters, e.g. random variables obeying various probability distributions. A priori convergence analysis is provided for the proposed method and the error bound by Kolmogorov N-width is improved from the recent work [Y. Maday, N.C. Nguyen, A.T. Patera and G.S.H. Pau, A general, multipurpose interpolation procedure: the magic points. Commun. Pure Appl. Anal. 8 (2009) 383-404]. We apply our method to geometric Brownian motion, exponential Karhunen-Loève expansion and reduced basis approximation of non-affine stochastic elliptic equations. We demonstrate its improved accuracy and efficiency over the empirical interpolation method, as well as sparse grid stochastic collocation method.1 aChen, Peng1 aQuarteroni, Alfio1 aRozza, Gianluigi uhttp://urania.sissa.it/xmlui/handle/1963/3502101318nas a2200121 4500008004100000245007900041210006900120520083000189100002001019700002201039700002101061856011401082 2013 eng d00aFree Form Deformation Techniques Applied to 3D Shape Optimization Problems0 aFree Form Deformation Techniques Applied to 3D Shape Optimizatio3 aThe purpose of this work is to analyse and study an efficient parametrization technique for a 3D shape optimization problem. After a brief review of the techniques and approaches already available in literature, we recall the Free Form Deformation parametrization, a technique which proved to be efficient and at the same time versatile, allowing to manage complex shapes even with few parameters. We tested and studied the FFD technique by establishing a path, from the geometry definition, to the method implementation, and finally to the simulation and to the optimization of the shape. In particular, we have studied a bulb and a rudder of a race sailing boat as model applications, where we have tested a complete procedure from Computer-Aided-Design to build the geometrical model to discretization and mesh generation.1 aKoshakji, Anwar1 aQuarteroni, Alfio1 aRozza, Gianluigi uhttps://www.math.sissa.it/publication/free-form-deformation-techniques-applied-3d-shape-optimization-problems01701nas a2200157 4500008004100000245007600041210006900117300001800186490000700204520113900211100002001350700002101370700002001391700002201411856011001433 2013 eng d00aReduced basis method for parametrized elliptic optimal control problems0 aReduced basis method for parametrized elliptic optimal control p aA2316–A23400 v353 aWe propose a suitable model reduction paradigm-the certified reduced basis method (RB)-for the rapid and reliable solution of parametrized optimal control problems governed by partial differential equations. In particular, we develop the methodology for parametrized quadratic optimization problems with elliptic equations as a constraint and infinite-dimensional control variable. First, we recast the optimal control problem in the framework of saddle-point problems in order to take advantage of the already developed RB theory for Stokes-type problems. Then, the usual ingredients of the RB methodology are called into play: a Galerkin projection onto a low-dimensional space of basis functions properly selected by an adaptive procedure; an affine parametric dependence enabling one to perform competitive offline-online splitting in the computational procedure; and an efficient and rigorous a posteriori error estimate on the state, control, and adjoint variables as well as on the cost functional. Finally, we address some numerical tests that confirm our theoretical results and show the efficiency of the proposed technique.1 aNegri, Federico1 aRozza, Gianluigi1 aManzoni, Andrea1 aQuarteroni, Alfio uhttps://www.math.sissa.it/publication/reduced-basis-method-parametrized-elliptic-optimal-control-problems00548nas a2200133 4500008004100000245009200041210006900133260001000202100001800212700002000230700002200250700002100272856012100293 2013 en d00aA Reduced Computational and Geometrical Framework for Inverse Problems in Haemodynamics0 aReduced Computational and Geometrical Framework for Inverse Prob bSISSA1 aLassila, Toni1 aManzoni, Andrea1 aQuarteroni, Alfio1 aRozza, Gianluigi uhttps://www.math.sissa.it/publication/reduced-computational-and-geometrical-framework-inverse-problems-haemodynamics00568nas a2200133 4500008004100000245010500041210006900146260001000215100001800225700002000243700002200263700002100285856012800306 2013 en d00aA reduced-order strategy for solving inverse Bayesian identification problems in physiological flows0 areducedorder strategy for solving inverse Bayesian identificatio bSISSA1 aLassila, Toni1 aManzoni, Andrea1 aQuarteroni, Alfio1 aRozza, Gianluigi uhttps://www.math.sissa.it/publication/reduced-order-strategy-solving-inverse-bayesian-identification-problems-physiological01509nas 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-elliptic01375nas a2200145 4500008004100000245010300041210006900144300001600213490000700229520080500236100001501041700002201056700002101078856013001099 2013 eng d00aA weighted reduced basis method for elliptic partial differential equations with random input data0 aweighted reduced basis method for elliptic partial differential a3163–31850 v513 aIn this work we propose and analyze a weighted reduced basis method to solve elliptic partial differential equations (PDEs) with random input data. The PDEs are first transformed into a weighted parametric elliptic problem depending on a finite number of parameters. Distinctive importance of the solution at different values of the parameters is taken into account by assigning different weights to the samples in the greedy sampling procedure. A priori convergence analysis is carried out by constructive approximation of the exact solution with respect to the weighted parameters. Numerical examples are provided for the assessment of the advantages of the proposed method over the reduced basis method and the stochastic collocation method in both univariate and multivariate stochastic problems.1 aChen, Peng1 aQuarteroni, Alfio1 aRozza, Gianluigi uhttps://www.math.sissa.it/publication/weighted-reduced-basis-method-elliptic-partial-differential-equations-random-input-data01501nas a2200157 4500008004100000245010600041210006900147260003100216520095600247653002301203100001801226700002001244700002201264700002101286856003601307 2012 en d00aBoundary control and shape optimization for the robust design of bypass anastomoses under uncertainty0 aBoundary control and shape optimization for the robust design of bCambridge University Press3 aWe review the optimal design of an arterial bypass graft following either a (i) boundary optimal control approach, or a (ii) shape optimization formulation. The main focus is quantifying and treating the uncertainty in the residual flow when the hosting artery is not completely occluded,\\r\\nfor which the worst-case in terms of recirculation e ffects is inferred to correspond to a strong ori fice flow through near-complete occlusion. A worst-case optimal control approach is applied to the steady\\r\\nNavier-Stokes equations in 2D to identify an anastomosis angle and a cu ed shape that are robust with respect to a possible range of residual \\r\\nflows. We also consider a reduced order modelling framework\\r\\nbased on reduced basis methods in order to make the robust design problem computationally feasible. The results obtained in 2D are compared with simulations in a 3D geometry but without model\\r\\nreduction or the robust framework.10ashape optimization1 aLassila, Toni1 aManzoni, Andrea1 aQuarteroni, Alfio1 aRozza, Gianluigi uhttp://hdl.handle.net/1963/633701643nas a2200157 4500008004100000245012600041210006900167260001300236520109700249653002201346100001801368700002001386700002201406700002101428856003601449 2012 en d00aGeneralized reduced basis methods and n-width estimates for the approximation of the solution manifold of parametric PDEs0 aGeneralized reduced basis methods and nwidth estimates for the a bSpringer3 aThe set of solutions of a parameter-dependent linear partial di fferential equation with smooth coe fficients typically forms a compact manifold in a Hilbert space. In this paper we review the generalized reduced basis method as a fast computational tool for the uniform approximation of the solution manifold. We focus on operators showing an affi ne parametric dependence, expressed as a linear combination of parameter-independent operators through some smooth, parameter-dependent scalar functions. In the case that the parameter-dependent operator has a dominant term in its affi ne expansion, one can prove the existence of exponentially convergent uniform approximation spaces for the entire solution manifold. These spaces can be constructed without any assumptions on the parametric regularity of the manifold \\r\\nonly spatial regularity of the solutions is required. The exponential convergence rate is then inherited by the generalized reduced basis method. We provide a numerical example related to parametrized elliptic\\r\\nequations con rming the predicted convergence rates.10asolution manifold1 aLassila, Toni1 aManzoni, Andrea1 aQuarteroni, Alfio1 aRozza, Gianluigi uhttp://hdl.handle.net/1963/634001580nas 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-hemodynamics00524nas a2200145 4500008004100000245006800041210006800109300001400177490000700191100001800198700001700216700002300233700001900256856010300275 2004 eng d00aCalculation of impulsively started incompressible viscous flows0 aCalculation of impulsively started incompressible viscous flows a877–9020 v461 aMarra, Andrea1 aMola, Andrea1 aQuartapelle, Luigi1 aRiviello, Luca uhttps://www.math.sissa.it/publication/calculation-impulsively-started-incompressible-viscous-flows