16 Octobre – Soutenance Thèse Lucas Lestandi
09 h30 Bat. A ESNCBP, 16 Rue Pey Berlan 33600 PESSAC, Amphi 1
Subject : Low rank approximation techniques and Reduced Order Modeling applied to some fluid dynamics problems
M. Mejdi AZAIEZ
Université de Bordeaux
Directeur de these
M. Gianluigi ROZZA
M. Elias CUETO
Université de Saragosse
M. Tomas CHACON REBOLLO
Université de Séville
CoDirecteur de these
M. Luc MIEUSSENS
Mme Marianne BERINGHIER
M. Samuele RUBINO
Université de Séville
In the last decades, numerical simulation has experienced tremendous improvements driven by massive growth of computing power. Exascale computing has been achieved this year and will allow solving ever more complex problems. But such large systems produce colossal amounts of data which leads to its own difficulties. Moreover, many engineering problems such as multiphysics or optimisation and control, require far more power that any computer architecture could achieve within the current scientific computing paradigm. In this thesis, we propose to shift the paradigm in order to break the curse of dimensionality by introducing decomposition and building reduced order models (ROM) for complex fluid flows.
First, we propose an extended review of data reduction techniques and intends to bridge between applied mathematics community and the computational mechanics one. Thus, founding bivariate separation is studied, including discussions on the equivalence of proper orthogonal decomposition (POD, continuous framework) and singular value decomposition (SVD, discrete matrices). Then a wide review of tensor formats and their approximation is proposed. Such work has already been provided in the literature but either on separate papers or into a purely applied mathematics framework. Here, we offer to the data enthusiast scientist a comparison of Canonical, Tucker, Hierarchical and Tensor train formats including their approximation algorithms. Their relative benefits are studied both theoretically and numerically thanks to the python library pydecomp that was developed during this thesis. A careful analysis of the link between continuous and discrete methods is performed. Finally, we conclude that for most applications ST-HOSVD is best when the number of dimensions d lower than four and TT-SVD (or their POD equivalent) when d grows larger.
Then we focus on a complex fluid dynamics flow, in particular the singular lid driven cavity at high Reynolds number. This flow exhibits a series of Hopf bifurcation which are known to be hard to capture accurately which is why a detailed analysis was performed both with classical tools and POD. Once this flow has been characterized, time-scaling, a new “physics based” interpolation ROM is presented on internal and external flows. This methods gives encouraging results while excluding recent advanced developments in the area such as EIM or Grassmann manifold interpolation.
ROM, Model reduction, Tensor decomposition , Interpolation