Reduced-order models at work: Industry and Medicine

Europe/Paris
Ada Lovelace (Inria Bordeaux Sud Ouest)

Ada Lovelace

Inria Bordeaux Sud Ouest

200 avenue de la ville tour, 33405 TALENCE
Angelo Iollo (Université de Bordeaux), Mejdi Azaiez (Bordeaux INP), Michel Bergmann (INRIA Bordeaux), Tommaso Taddei (Inria Bordeaux)
Description

 

The workshop aims at tightening the link between fundamental research on Model Reduction and clinical and industrial applications. Topics of the presentations will include reduced basis methods for parameterized systems, Proper Orthogonal Decomposition, nonlinear compression methods, data assimilation and geometry inference methods.

This workshop is supported by: Inria - EDF - Bordeaux INP - Université de Bordeaux - CNRS- IdEx - HORIZON H2020

Participants
  • Agustina Giuliodori Picco
  • Alberto Remigi
  • Alejandro Bandera
  • amaury deldicque
  • Andrea Ferrero
  • Andrey Popov
  • Angelo Iollo
  • Angélique Ponçot
  • Anthony Patera
  • Bernd Noack
  • Bernd Noack
  • Catharina Czech
  • Duc Khai Nguyen Pham
  • Eki AGOUZAL
  • Ettorche Rayen
  • Florian Bernard
  • georges Oppenheim
  • Gianluigi Rozza
  • Gianluigi Rozza
  • Giulia Sambataro
  • Guilhem Ferté
  • Haysam Telib
  • Jonathan Pham
  • Jose Raul Bravo Martinez
  • Juan Valverde
  • Karl Maroun
  • Laurent Cordier
  • Lei Zhang
  • luca abergo
  • Luca Berti
  • Ludovica Saccaro
  • Maksim Elizarev
  • Manuel Bernardino del Pino
  • Marie BILLAUD FRIESS
  • Matthew Zahr
  • Mejdi Azaiez
  • Michel Bergmann
  • Michele Giuliano Carlino
  • Mirco Ciallella
  • NADJIMA ISSA
  • Pedro Morel Rosa
  • Pierre Mollo
  • RABAH Mohammed
  • Rémi Abgrall
  • Saigre Thomas
  • Samuele Rubino
  • Sean Reiter
  • Simone Camarri
  • Spenser Anderson
  • Thomas Milcent
  • Tommaso Taddei
  • Tomás Chacón Rebollo
  • Traian Iliescu
  • Victor Zucatti da Silva
  • William Wolf
  • Yassine Chahid
    • 9:00 AM
      Welcome and coffee Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE
    • General session: Maday + Vignon-Clémentel Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE

      Maday

    • 1
      Some new developments on the non invasive reduced basis method Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE

      The efficient implementation of reduced basis methods relying on a high fidelity discretization method to compute the elements of the reduced basis, requires to enter within the code, offline, so that the online solution can be produced very rapidly. Since this is sometimes impossible, in particular for codes used in industrial framework, we have proposed a Non Invasive alternative where the code is used at two stages : a) offline to build the reduced basis, b) online to, first, compute a coarse approximation using the code with few degrees of freedom (thus more rapidly than with the high fidelity requirement) then by processing this coarse solution to improve the accuracy and fulfil the high fidelity requirement at much mower cost. This method that first appeared with its numerical analysis in [1] for finite element approximations has been generalised in various directions (finite volume, truncation of the domain, different new rectifications) and implemented in a library in the frame of a joined project with large and medium size companies [2].
      In this paper I will present some of the new concepts in this drection

      [1] Rachida Chakir, Yvon Maday :A two-grid finite-element/reduced basis scheme for the approximation of the solution of parameter dependent PDE. In 9e Colloque national en calcul des structures. ISO 690 Y. (2009, May).
      [2] Elise Grosjean : Variations and further developments on the Non-Intrusive Reduced Basis two-grid method, PhD Thesis at Sorbonne Université - 2022

      Elise Grosjean and Yvon Maday

      Speaker: Prof. Yvon Maday (UPMC Paris 6)
    • 2
      Hemodynamics examples of reduced order models for clinical applications Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE
      Speaker: Prof. Irène Vignon-Clémentel (Inria Saclay)
    • 12:10 PM
      Lunch Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE
    • General session: Iliescu + Mathelin + Enchery Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE

      Maday

      • 1:40 PM
        Welcome and coffee
    • 3
      ROM Closures and Stabilizations for Under-Resolved Turbulent Flows Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE

      In this talk, I will survey reduced order model (ROM) closures and
      stabilizations for under-resolved turbulent flows. Over the past
      decade, several closure and stabilization strategies have been
      developed to tackle the ROM inaccuracy in the convection-dominated,
      under-resolved regime, i.e., when the number of degrees of freedom is
      too small to capture the complex underlying dynamics. I will present
      regularized ROMs, which are stabilizations that employ spatial
      filtering to alleviate the spurious numerical oscillations generally
      produced by standard ROMs in the convection-dominated, under-resolved
      regime. I will also survey three classes of ROM closures, i.e.,
      correction terms that increase the ROM accuracy: (i) functional
      closures, which are based on physical insight; (ii) structural
      closures, which are developed by using mathematical arguments; and
      (iii) data-driven closures, which leverage available data. Throughout
      my talk, I will highlight the impact made by data on classical
      numerical methods over the past decade. I will also emphasize the role
      played by physical constraints in data-driven modeling of ROM closures
      and stabilizations.

      Speaker: Prof. Traian Iliescu (Virginia Tech)
    • 4
      Field reconstruction using manifold learning and structure-preserving metrics Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE

      The problem of estimating the state of a physical system is ubiquitous in science. However observations are always limited so that the high-dimensional state cannot be observed and the associated mathematical problem is ill-posed. Popular workarounds include dimension reduction and regularization by imposing some structure to the class of elements in which the estimation is sought.
      We here rely on a purely data-driven approach and learn the map between extended measurements and the nonlinear manifold the state vector lies on. Specifically, we use embedding to address the non-Markovianity of the raw measurements. Combined with multi-kernel learning, it results in high-dimensional measurement features. The state vector nonlinear manifold is approximated and the map from measurement features to the estimated state is the solution to a Sylvester equation.
      The methodology is illustrated with the estimation of a fluid flow field from a few wall-mounted sensors.

      Speaker: Dr Lionel Mathelin (LIMSI-CNRS)
    • 5
      Use of reduced basis techniques for two-phase flows in porous media Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE
      Speaker: Dr Guillaume Enchery (IFPEN)
    • 4:10 PM
      Coffee break Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE
    • Poster session Salle de détente (Institut de mathématiques de Bordeaux)

      Salle de détente

      Institut de mathématiques de Bordeaux

      351 cours de la libération, 33405 TALENCE Cédex
    • General session: Farhat + Zahr Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE

      Maday

    • 6
      Quadratic Approximation Manifold for Mitigating the Kolmogorov Barrier in Nonlinear Projection-Based Model Order Reduction Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE

      A quadratic approximation manifold is presented for performing nonlinear, projection-based, model order reduction (PMOR). It constitutes a departure from the traditional affine subspace approximation aimed at mitigating the Kolmogorov barrier for nonlinear PMOR, particularly for convection-dominated transport problems. It builds on the data-driven approach underlying the traditional construction of projection-based reduced-order models (PROMs); is application-independent; is linearization-free; and therefore is robust for highly nonlinear problems. Most importantly, this approximation leads to quadratic PROMs that deliver the same accuracy as their traditional counterparts using however the square root of their dimension. The computational advantages of the proposed high-order approach to nonlinear PMOR over the traditional approach are highlighted for the detached-eddy simulation-based prediction of the Ahmed body turbulent wake flow, which is a popular CFD benchmark problem in the automotive industry. For a fixed accuracy level, these advantages include: a reduction of the total offline computational cost by a factor greater than five; a reduction of its online wall clock time by a factor greater than 32; and a reduction of the wall clock time of the underlying high-dimensional model by a factor greater than two orders of magnitude.

      co-author: Joshua Barnett

      Speaker: Prof. Charbel Farhat (Stanford)
    • 7
      Model reduction of convection-dominated partial differential equations via optimization-based implicit feature tracking Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE

      Partial differential equations (PDEs) that model convection-dominated phenomena often arise in engineering practice and scientific applications, ranging from the study of high-speed, turbulent flow over vehicles to wave propagation through solid media. The solutions of these equations are characterized by local features or disturbances that propagate throughout the domain as time evolves or a system parameter varies. Numerical methods to approximate these solutions require stabilization and fine, usually adaptive, grids to adequately resolve the local features, which lead to expensive discretizations with a large number of degrees of freedom. Projection-based model reduction methods tend to be ineffective in reducing the computational cost of such problems due to a slowly decaying Kolmogorov n-width of the solution manifold.

      To avoid the fundamental linear reducibility limitation associated with convection-dominated problems, we construct a nonlinear approximation by composing a low-dimensional linear space with a parametrized domain mapping. The linear space is constructed using the method of snapshots and POD; prior to compression, each snapshot is composed with a mapping that causes its local features to align (same spatial location) with the corresponding features in all other snapshots. The parametrized domain mapping is chosen such that the local features present in the linear space deform to the corresponding features in the solution being approximated, effectively removing the convection-dominated nature of the problem. The domain mapping is determined implicitly through the solution of a residual minimization problem, rather than relying on explicit sensing/detection. We provide numerous numerical experiments to demonstrate the effectivity of the proposed method on benchmark problems from computational fluid dynamics.

      Speaker: Prof. Mattew Zahr (Notre Dame)
    • 10:40 AM
      Coffee break Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE
    • 8
      From handcrafted Galerkin models to automated cluster models Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE

      Reduced-order models (ROM) are of paramount importance for physical understanding, data
      compression, estimation, control and optimization. Over a century ago, simple dynamical models
      of coherent structures have been facilitated by stability theory (Orr-Sommerfeld equation
      1907) and by vortex models (von K arm an 1911). Data-driven reduced-order modeling of coherent
      structures has been pioneered by Aubry et al. (1988) with a celebrated POD model of
      the turbulent boundary layer. Since then, machine learning [1] has significantly simplified and
      enriched the spectrum of possibilities for data-driven ROM.
      In this talk, we exemplify different ROM approaches for the fluidic pinball [2], the wake flow
      behind a cluster of three parallel cylinders on an equilateral triangle pointing upstream. The
      flow may be actuated by rotating cylinders. First, the transition scenario of the unforced fluidic
      pinball is modeled with a five-mode mean-field Galerkin model. This model comprises
      successive Hopf and pitchfork bifucations, which are typical for a number of wake flows. Second,
      a cluster-based network model (CNM) [3, 4] is presented describing the fluidic pinball
      wake with actuation as free input, employing thousand differently actuated pinball simulations.
      CNM yields a robust dynamics from a fully automatable procedure. Finally, perspectives of
      ROM for common tasks of data management are given.
      [1] BRUNTON, S. L., NOACK, B. R. & KOUMOUTASKOS, P. 2020 Ann. Rev. Fluid Mech. 52:477–
      508.
      [2] DENG, N., NOACK, B. R., MORZY N SKI, M., & PASTUR, L. R. 2020 Low-order model for
      successive bifurcations of the fluidic pinball. J. Fluid Mech. 884, A37:1–41.
      [3] LI, H., FERNEX, D., SEMAAN, R., TAN, J., MORZY N SKI, M. & NOACK, B. R. 2021 Clusterbased
      network model. J. Fluid Mech. 906, A21:1–41.
      [4] FERNEX, D., NOACK, B. R. & SEMAAN, R. 2021 Cluster-based network model—From snapshots
      to complex dynamical systems. Science Advances (online).

      Bernd R. Noack 1, Nan Deng1, Daniel Fernex2
      Luc Pastur3, Richard Semaan4 Marek Morzy  nski5
      1 Harbin Institute of Technology, Shenzhen, China,
      2 EPFL, Switzerland,
      3 ENSTA, Paris,
      4 TU Braunschweig, Germany,
      5 Poznan University of Technology, Poland

      Speaker: Prof. Bernd Noack (Harbin Institute ot Technology, Shenzhen, China)
    • General session: Noack + Chacón Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE

      Maday

    • 9
      Reduced basis Smagorinsky turbulence models Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE
      Speaker: Prof. Tomás Chacón (Instituto de Matemáticas de la Universidad de Sevilla)
    • 12:50 PM
      Lunch Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE
    • 10
      Fear and loathing in Aerodynamics: Are ROMs chimeras or game-changers? Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE

      Over the last years promising advances in reduced order modelling and machine learning have been reported by the scientific community.
      Within this talk we will try to project (some of) these findings on different industrial settings to discuss their potential.
      We will try to make this evaluation using all metrices that are relevant in industry to identify applications where ROMs can have a real impact.

      Speaker: Dr Haysam Telib (Optimad)
    • Industrial Presentations: perspectives for next ten years: Telib + Bernard + Ferté Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE
    • 11
      TBD Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE
      Speaker: Dr Florian Bernard (Nurea)
    • 12
      TBD Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE
      Speaker: Dr Guilhem Ferté (EDF)
    • 4:00 PM
      Coffee break Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE
    • Round table Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE
    • 13
      Beyond PCA by explicitely taking into account system symmetries Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE

      Linear principal component analysis (PCA) experiences an increase in the dimensionality of the latent space when it is applied to configurations that exhibit symmetries. In this study, we introduce a novel machine learning embedding, which uses spatial transformer networks and siamese networks to account for continuous and discrete symmetries, respectively. This embedding, which we term symmetry-aware PCA, will be applied to three configurations: Burger's equation exhibiting a continuous translation symmetry, flow in sudden expansion, a discrete reflexional symmetry, and Kolmogorov Flow which combines discrete shift-reflect and continuous translation symmetries. We will show a drastic increase in the number of modes required to represent given trajectories.

      Simon Kneer, Taraneh Sayadi, Denis Sipp, Peter J. Schmid, Georgios Rigas

      Speaker: Prof. Denis Sipp (ONERA)
    • General session: Sipp + Lombardi Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE

      Maday

    • 14
      Tensor methods for high-dimensional problems and model reduction Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE

      We present several contributions related to Tensor methods for high-dimensional problems and discuss how they are inherently related to model reduction.
      In particular, we will show several ways to introduce a principle of adaptivity, making tensor representations more suitable to parsimoniously represent certain solutions sets. In the last part of the talk, a contribution on a possible way to exploit a tensor representation in state estimation is presented. In this, we show that variational and sequential state estimation methods can be derived after casting state estimation as an optimal recovery problem, using tensors to have a space-time representation of the solutions set.

      Speaker: Dr Damiano Lombardi (Inria Paris)
    • 11:10 AM
      Coffee break Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE
    • General session: Rozza + Perotto Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE

      Maday

    • 15
      Model order reduction for CFD: state of the art, advances in applications Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE
      Speaker: Prof. Gianlugi Rozza (SISSA)
    • 16
      POD- and RB-Hierarchical Model Reduction Techniques in a Parametrized Setting Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE

      Different methods have been proposed in the scientific panorama to offer a compromise between modeling accuracy and computational efficiency. Reduced order models represent a widespread solution in such a direction. In this presentation, we focus on the Hierarchical Model (HiMod) reduction technique, which proved to be an effective approach to discretize CFD configurations where a principal horizontal dynamics overwhelms the transverse ones (e.g., when modeling hemodynamics or signal propagation in waveguides). In particular, we address the generalization of the HiMod procedure to a parametrized setting. We propose two alternative approaches, which combine HiMod with Proper Orthogonal Decomposition (POD) and the Reduced Basis (RB) method, respectively. The two strategies will be analyzed and cross-compared in order to identify the associated pros and cons.

      Speaker: Prof. Simona Perotto (Politecnico Milano)
    • 1:20 PM
      Lunch Ada Lovelace

      Ada Lovelace

      Inria Bordeaux Sud Ouest

      200 avenue de la ville tour, 33405 TALENCE