Francesco Ballarin

Shape Optimization by Free-Form Deformation: Existence Results and Numerical Solution for Stokes Flows

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

Fast simulations of patient-specific haemodynamics of coronary artery bypass grafts based on a POD-Galerkin method and a vascular shape parametrization

In this work a reduced-order computational framework for the study of haemodynamics in three-dimensional patient-specific configurations of coronary artery bypass grafts dealing with a wide range of scenarios is proposed. We combine several efficient algorithms to face at the same time both the geometrical complexity involved in the description of the vascular network and the huge computational cost entailed by time dependent patient-specific flow simulations. Medical imaging procedures allow to reconstruct patient-specific configurations from clinical data. A centerlines-based parametrization is proposed to efficiently handle geometrical variations. POD-Galerkin reduced-order models are employed to cut down large computational costs. This computational framework allows to characterize blood flows for different physical and geometrical variations relevant in the clinical practice, such as stenosis factors and anastomosis variations, in a rapid and reliable way. Several numerical results are discussed, highlighting the computational performance of the proposed framework, as well as its capability to carry out sensitivity analysis studies, so far out of reach. In particular, a reduced-order simulation takes only a few minutes to run, resulting in computational savings of 99 % of CPU time with respect to the full-order discretization. Moreover, the error between full-order and reduced-order solutions is also studied, and it is numerically found to be less than 1 % for reduced-order solutions obtained with just O(100) online degrees of freedom.

Isogeometric analysis-based reduced order modelling for incompressible linear viscous flows in parametrized shapes

In this work we provide a combination of isogeometric analysis with reduced order modelling techniques, based on proper orthogonal decomposition, to guarantee computational reduction for the numerical model, and with free-form deformation, for versatile geometrical parametrization. We apply it to computational fluid dynamics problems considering a Stokes flow model. The proposed reduced order model combines efficient shape deformation and accurate and stable velocity and pressure approximation for incompressible viscous flows, computed with a reduced order method. Efficient offline–online computational decomposition is guaranteed in view of repetitive calculations for parametric design and optimization problems. Numerical test cases show the efficiency and accuracy of the proposed reduced order model.

Reduced-Order Semi-Implicit Schemes for Fluid-Structure Interaction Problems

POD–Galerkin reduced-order models (ROMs) for fluid-structure interaction problems (incompressible fluid and thin structure) are proposed in this paper. Both the high-fidelity and reduced-order methods are based on a Chorin-Temam operator-splitting approach. Two different reduced-order methods are proposed, which differ on velocity continuity condition, imposed weakly or strongly, respectively. The resulting ROMs are tested and compared on a representative haemodynamics test case characterized by wave propagation, in order to assess the capabilities of the proposed strategies.

On a certified smagorinsky reduced basis turbulence model

In this work we present a reduced basis Smagorinsky turbulence model for steady flows. We approximate the nonlinear eddy diffusion term using the empirical interpolation method (cf. [M. A. Grepl et al., ESAIM Math. Model. Numer. Anal., 41 (2007), pp. 575--605; Barrault et al., C. R. Acad. Sci. Paris Ser. I Math., 339 (2004), pp. 667--672]) and the velocity-pressure unknowns by an independent reduced-basis procedure. This model is based upon an a posteriori error estimation for a Smagorinsky turbulence model. The theoretical development of the a posteriori error estimation is based on [S. Deparis, SIAM J. Sci. Comput., 46 (2008), pp. 2039--2067] and [A. Manzoni, ESAIM Math. Model. Numer. Anal., 48 (2014), pp. 1199--1226], according to the Brezzi--Rappaz--Raviart stability theory, and adapted for the nonlinear eddy diffusion term. We present some numerical tests, programmed in FreeFem++ (cf. [F. Hecht, J. Numer. Math., 20 (2012), pp. 251--265]), in which we show a speedup on the computation by factor larger...

Advances in geometrical parametrization and reduced order models and methods for computational fluid dynamics problems in applied sciences and engineering: overview and perspectives

Several problems in applied sciences and engineering require reduction techniques in order to allow computational tools to be employed in the daily practice, especially in iterative procedures such as optimization or sensitivity analysis. Reduced order methods need to face increasingly complex problems in computational mechanics, especially into a multiphysics setting. Several issues should be faced: stability of the approximation, efficient treatment of nonlinearities, uniqueness or possible bifurcations of the state solutions, proper coupling between fields, as well as offline-online computing, computational savings and certification of errors as measure of accuracy. Moreover, efficient geometrical parametrization techniques should be devised to efficiently face shape optimization problems, as well as shape reconstruction and shape assimilation problems. A related aspect deals with the management of parametrized interfaces in multiphysics problems, such as fluid-structure interaction problems, and also a domain decomposition based approach for complex parametrized networks. We present some illustrative industrial and biomedical problems as examples of recent advances on methodological developments.

Model Reduction for Parametrized Optimal Control Problems in Environmental Marine Sciences and Engineering

We propose reduced order methods as a suitable approach to face parametrized optimal control problems governed by partial differential equations, with applications in en- vironmental marine sciences and engineering. Environmental parametrized optimal control problems are usually studied for different configurations described by several physical and/or geometrical parameters representing different phenomena and structures. The solution of parametrized problems requires a demanding computational effort. In order to save com- putational time, we rely on reduced basis techniques as a reliable and rapid tool to solve parametrized problems. We introduce general parametrized linear quadratic optimal control problems, and the saddle-point structure of their optimality system. Then, we propose a POD-Galerkin reduction of the optimality system. Finally, we test the resulting method on two environmental applications: a pollutant control in the Gulf of Trieste, Italy and a solution tracking governed by quasi-geostrophic equations, in its linear and nonlinear version, describing North Atlantic Ocean dynamic. The two experiments underline how reduced order methods are a reliable and convenient tool to manage several environmental optimal control problems, for different mathematical models, geographical scale as well as physical meaning.

A Weighted POD Method for Elliptic PDEs with Random Inputs

In this work we propose and analyze a weighted proper orthogonal decomposition method to solve elliptic partial differential equations depending on random input data, for stochastic problems that can be transformed into parametric systems. The algorithm is introduced alongside the weighted greedy method. Our proposed method aims to minimize the error in a

Supremizer stabilization of POD–Galerkin approximation of parametrized steady incompressible Navier–Stokes equations

L^2