Reduced formulation of a steady fluid-structure interaction problem with parametric coupling

We propose a two-fold approach to model reduction of fluid-structure interaction. The state equations for the fluid are solved with reduced basis methods. These are model reduction methods for parametric partial differential equations using well-chosen snapshot solutions in order to build a set of global basis functions. The other reduction is in terms of the geometric complexity of the moving fluid-structure interface. We use free-form deformations to parameterize the perturbation of the flow channel at rest configuration. As a computational example we consider a steady fluid-structure interaction problem: an incmpressible Stokes flow in a channel that has a flexible wall.Comment: 10 pages, 3 figure

Nonlinear model reduction via discrete empirical interpolation

This thesis proposes a model reduction technique for nonlinear dynamical systems based upon combining Proper Orthogonal Decomposition (POD) and a new method, called the Discrete Empirical Interpolation Method (DEIM). The popular method of Galerkin projection with POD basis reduces dimension in the sense that far fewer variables are present, but the complexity of evaluating the nonlinear term generally remains that of the original problem. DEIM, a discrete variant of the approach from [11], is introduced and shown to effectively overcome this complexity issue. State space error estimates for POD-DEIM reduced systems are also derived. These [Special characters omitted.] error estimates reflect the POD approximation property through the decay of certain singular values and explain how the DEIM approximation error involving the nonlinear term comes into play. An application to the simulation of nonlinear miscible flow in a 2-D porous medium shows that the dynamics of a complex full-order system of dimension 15000 can be captured accurately by the POD-DEIM reduced system of dimension 40 with a factor of [Special characters omitted.] (1000) reduction in computational time

Reduced Models for Optimal Control, Shape Optimization and Inverse Problems in Haemodynamics

The objective of this thesis is to develop reduced models for the numerical solution of optimal control, shape optimization and inverse problems. In all these cases suitable functionals of state variables have to be minimized. State variables are solutions of a partial differential equation (PDE), representing a constraint for the minimization problem. The solution of these problems induce large computational costs due to the numerical discretization of PDEs and to iterative procedures usually required by numerical optimization (many-query context). In order to reduce the computational complexity, we take advantage of the reduced basis (RB) approximation for parametrized PDEs, once the state problem has been reformulated in parametrized form. This method enables a rapid and reliable approximation of parametrized PDEs by constructing low-dimensional, problem-specific approximation spaces. In case of PDEs defined over domains of variable shapes (e.g. in shape optimization problems) we need to introduce suitable, low-dimensional shape parametrization techniques in order to tackle the geometrical complexity. Free-Form Deformations and Radial-Basis Functions techniques have been analyzed and successfully applied with this aim. We analyze the reduced framework built by coupling these tools and apply it to the solution of optimal control and shape optimization problems. Robust optimization problems under uncertain conditions are also taken into consideration. Moreover, both deterministic and Bayesian frameworks are set in order to tackle inverse identification problems. As state equations, we consider steady viscous flow problems described by Stokes or Navier-Stokes equations, for which we provide a detailed analysis and construction of RB approximation and a posteriori error estimation. Several numerical test cases are also illustrated to show efficacy and reliability of RB approximations. We exploit this general reduced framework to solve some optimization and inverse problems arising in haemodynamics. More specifically, we focus on the optimal design of cardiovascular prostheses, such as bypass grafts, and on inverse identification of pathological conditions or flow/shape features in realistic parametrized geometries, such as carotid artery bifurcations