1 of 89

Institut Polytechnique de Paris��Habilitation à Diriger des Recherches

Analyse de quelques problèmes d’instabilité en mécanique

par Romain LAGRANGE

le 12/12/2023

Commission d’examen

M. Jean Marc CHOMAZ

Mme. Lydia BOUROUIBA

M. Bruno COCHELIN

M. Stéphane ETIENNE

M. Xavier AMANDOLESE

Mme. Marianna BRAZA

M. Arnaud LAZARUS

Président

Rapporteur

Rapporteur

Rapporteur

Examinateur

Examinatrice

Examinateur

2 of 89

Plan

2

Dossier personnel

Positionnement de mon activité de recherche

Instabilité d’un fluide dans un cylindre en précession

Instabilité fluide-élastique d’un faisceau de tubes

Instabilité de flambage / plissement

Perspectives de recherche

Hors programme

Instabilité de flottement

Instabilité de flambage sur substrat non-linéaire

3 of 89

Plan

3

Dossier personnel

Positionnement de mon activité de recherche

Instabilité d’un fluide dans un cylindre en précession

Instabilité fluide-élastique d’un faisceau de tubes

Instabilité de flambage / plissement

Perspectives de recherche

Hors programme

Instabilité de flottement

Instabilité de flambage sur substrat non-linéaire

Instabilité de plissement 3D

4 of 89

Dossier personnel

4

Education

Année

Ecole

Diplôme

Spécialité

2009

Aix-Marseille Université

Doctorat

Mécanique des fluides

(Instabilité de précession)

2006

Aix-Marseille Université

Master 2

Mécanique des fluides

(Instabilité de flottement)

2006

Centrale Marseille

Ingénieur

Mécanique des solides et des structures

Parcours professionnel

Année

Entreprise / Université

Poste

Recherche

2015-…

CEA

Ingénieur-Chercheur

Instabilité fluide-élastique

2012-2015

MIT

Instructeur

Instabilité de plissement

Instabilité de précession

2010-2011

IFPEN

Post-doctorat

Instabilité de flambage sur substrat non-linéaire

5 of 89

Dossier personnel

5

Enseignements (depuis 2012)

En cours

Domaine

Fonction

Ecole / Université

2022-…

MMC

Responsable de module

Cours de 15h

ENPC

2020-…

Mathématiques

Colleur

EPF

2017-…

Mathématiques

Examinateur

Concours Bachelor, X

Passés

Domaine

Cours

TD

Niveau

Ecole / Université

2019-2021

Mathématiques

Concours Centrale

2017-2018

Mathématiques

X

L1

CNAM Versailles

2018

Mathématiques

X

L1

IUT Orsay

2017

MMC

X

M1

UPMC

2016-2018

Mathématiques

X

L1

Polytech-UPMC

2015-2018

Dynamique non-linéaire et instabilités

X

Cours de 24h

M2

UPMC

2012-2015

Mathématiques

X

Cours de 50h

X

Under.

MIT

6 of 89

Dossier personnel

6

Journal

Nombre d’articles

Journal of Fluids and Structures

6

Journal of Fluid Mechanics

4

Journal of Applied Mechanics

2

Comptes Rendus de Mécanique

2

Theoretical and Applied Mechanics Letters

2

Nature Materials

1

Physics of Fluids

1

Physical Review Letters

1

Journal of the Mechanics and Physics of Solids

1

Proceedings of the National Academy of Sciences

1

International Journal of Mechanical Sciences

1

European Journal of Mechanics - B/Fluids

1

Mechanics and Industry

1

24

(19 depuis thèse en 2009)

(10 ces 5 dernières années)

Publications

Cinquantaine de communications dans congrès nationaux et int.

Monographie CEA: vibrations sous écoulement

Membre du comité d’organisation du FIV2022

7 of 89

Dossier personnel

7

Activité d’encadrement

Etudiant

Ecole

Année

Master

Thèse

L. Plado Costante

ENPC

2023

X

M. Kocher

IPP

2021-2024

X

C. Bazin

IPP

2020-2023

X

L. Lorand

ESTACA

2022

X

E. Diaferia

ENSAM

2021

X

E. Beaugendre

POLYTECH

2021

X

A. Epalle

INSA

2021

X

D. Panunzio

ENSAM

2020

X

M. Monteverde

ENSTA

2020

X

D. Zegarra

ESILV

2020

X

S. Ambriola

ENSAM

2019

X

S. Tsitkov

MIT

2013 (1 an)

X

A. Nasto

MIT

2013

X

11

2

(+ 1 en 2024)

8 of 89

Plan

8

Dossier personnel

Positionnement de mon activité de recherche

Instabilité d’un fluide dans un cylindre en précession

Instabilité fluide-élastique d’un faisceau de tubes

Instabilité de flambage / plissement

Perspectives de recherche

Hors programme

Instabilité de flottement

Instabilité de flambage sur substrat non-linéaire

Instabilité de plissement 3D

9 of 89

Positionnement de mon activité de recherche

9

Modélisation (fluides / solides)

Adim.

Maths appli. (recherche , analyse stab.)

Système dynamique 1D

Diagramme de bifurcation

Diagrammes de phase

Nouveaux critères

Etudes paramétriques Comparaisons expérience / numérique

Approche générale

Problème fondamental / industriel

Bifurcation en

(Type selle)

10 of 89

Plan

10

Dossier personnel

Positionnement de mon activité de recherche

Instabilité d’un fluide dans un cylindre en précession

Instabilité fluide-élastique d’un faisceau de tubes

Instabilité de flambage / plissement

Perspectives de recherche

Hors programme

Instabilité de flottement

Instabilité de flambage sur substrat non-linéaire

Instabilité de plissement 3D

11 of 89

Instabilité de précession �(2006-2009, 2012-2015)

11

Principaux collaborateurs

C. Eloy, P. Meunier

12 of 89

Instabilité de précession

12

Cylindre en précession

θ

Ω2

Ω1

H

2R

Six paramètres

Longueur caractéristique

Temps caractéristique

Quatre nombres sans dimension

Rapport d’aspect

Nombre de Rossby (convectif / Coriolis)

Rapport de fréquences

Nombre de Reynolds

13 of 89

Instabilité de précession

13

Navier-Stokes

Solution ordre

Modes de Kelvin forcés

+

+

+

Cisaillement

Vorticité Mode 1

Vorticité Mode 2

Mode 1

Mode 2

Relation de dispersion

Nombre d’onde axial,

Rapport de fréquence,

Résonance du mode 1

Amplitude du mode1

Rapport de fréquence,

14 of 89

Instabilité de précession

14

1st res. Mode 1

PIV vorticité axiale, écoulement instable

Mécanisme de résonance triadique

Mode de base

Mode 1 de l’instabilité

Mode 2 de l’instabilité

15 of 89

Instabilité de précession

15

Ecoulement de base

Modes de Kelvin libres

(Modes de l’instabilité)

Navier-Stokes (perturbation 2 temps)

Seuil d’instabilité (analytique)

Ordre 0 de NS : modes de Kelvin libres

Ordre 1 de NS : équations d’amplitude

+

+

Diagramme de stabilité

16 of 89

Instabilité de précession

16

Introduction d’un mode géostrophique (amplitude )

Attracteur étrange dans l’espace des phases

et

Diagramme de bifurcation

Nombre de Reynolds

Bif. 1

Bif. 2

17 of 89

Instabilité de précession

17

Visualisation du mélange d’un colorant dans un cylindre en précession.

Crédit : P. Meunier, brevet soft-mixer

Avant instabilité: pas de mélange

Après instabilité: mélange « optimal »

18 of 89

Plan

18

Dossier personnel

Positionnement de mon activité de recherche

Instabilité d’un fluide dans un cylindre en précession

Instabilité fluide-élastique d’un faisceau de tubes

Instabilité de flambage / plissement

Perspectives de recherche

Hors programme

Instabilité de flottement

Instabilité de flambage sur substrat non-linéaire

Instabilité de plissement 3D

19 of 89

Instabilité de flambage / plissement�(2012-2015)

19

Principaux collaborateurs

J. Dunkel, P. Reis, N. Stoop, D. Terwagne, F. Lopez Jimenez

20 of 89

Instabilité de flambage / plissement

20

Configuration plane, sous compression uni-axiale

Influence d’une courbure du substrat?

Flambage

Plissement

Diagramme de phase

Wang (2008)

21 of 89

21

1. Problem formulation

Rapport d’aspect

Taille de la cavité

Rapport de rigidités

Schéma d’un anneau sur un substrat 2D

Sept paramètres

Déplacements

Instabilité de flambage / plissement

Trois nombres sans dimension

22 of 89

22

2. Modeling

Développement de Koiter

Equilibres statiques

Equations d’Euler-Lagrange

Substrat (élasticité linéaire 2D, Airy, continuité déplacement)

Energie potentielle du film

Instabilité de flambage / plissement

Cinématique Euler-Bernoulli, déplacement ligne moyenne

Energie de déformation élastique

Travail des efforts mécaniques extérieurs (substrat)

et

23 of 89

23

2. Modeling

Ordre 1

Instabilité de flambage / plissement

Ordre 0

Mode global (m=2) ou mode de plissement (m>2)

Solution axisymétrique

Transition

Diagramme de stabilité

Rapport d’aspect,

Pression critique adim.,

Diagramme de phase

Rapport d’aspect,

24 of 89

24

Chargement isotrope du film

augmente

Mode hexagonal

Mode mixte

Mode labyrinthique

Instabilité de flambage / plissement

Equation d’Euler-Lagrange

Cinématique de Koiter, déplacement feuille moyenne

Diagramme de phase

Modes de l’instabilité

Théorie

Expériences

25 of 89

Instabilité de flambage / plissement

25

Motivation

Contrôle actif du coefficient de trainée via un changement de la rugosité de la surface

Trainée

Dépressurisation

Ecoulement

Dispositif expérimental Crédit: D. Terwagne

Stable

Instable

Dépressurisation,

Coefficient de trainée,

Coefficient de trainée

Voiture aérodynamique (???)

26 of 89

Plan

26

Dossier personnel

Positionnement de mon activité de recherche

Instabilité d’un fluide dans un cylindre en précession

Instabilité fluide-élastique d’un faisceau de tubes

Instabilité de flambage / plissement

Perspectives de recherche

Hors programme

Instabilité de flottement

Instabilité de flambage sur substrat non-linéaire

Instabilité de plissement 3D

27 of 89

Instabilité fluide-élastique en faisceau�(2015-…)

27

Principaux collaborateurs

P. Piteau, X. Delaune, D. Panunzio, J. Antunes,

Framatome

Echangeur de chaleur

Vibrations des tubes en U sous l’effet d’un écoulement diphasique

Vidéo Framatome

28 of 89

Instabilité fluide-élastique en faisceau

28

Tube monté sur deux lames flexibles

Ecoulement transverse

Monophasique ou diphasique (eau-air)

Vibrations en portance ou trainée

Système de contrôle actif

DIVA (5x3 tubes, pas carré)

29 of 89

29

1. Problem formulation

Treize paramètres

Ecoulement diphasique transverse

Vibrations en portance, pas carré

Instabilité fluide-élastique en faisceau

Neuf nombres sans dimension

Nombre de Stokes

Nombre de Reynolds

Nombre de Scruton

Pas réduit

Rapport d’aspect

Nombre de Webber

30 of 89

Instabilité fluide-élastique en faisceau

30

Seuil d’instabilité

STABLE

INSTABLE

Coefficient d’amortissement

Faisceau, tube, fluide fixés

Force fluide-élastique

Equation du mouvement

31 of 89

Instabilité fluide-élastique en faisceau

31

Phase 1

Phase 2

STABLE

INSTABLE

STABLE

INSTABLE

Diphasique avec

deux phases identiques

32 of 89

Instabilité fluide-élastique en faisceau

32

Diphasique avec

deux phases identiques

Diphasique

?

STABLE

STABLE

INSTABLE

INSTABLE

33 of 89

Instabilité fluide-élastique en faisceau

33

FINAL RESULT

Diagramme de stabilité monophasique

(méthode directe)

Stable

Instable

Stable

Instable

Test expérimental de la loi

Diagramme de stabilité diphasique

(méthode indirecte)

3 jeux de lames

34 of 89

Plan

34

Dossier personnel

Positionnement de mon activité de recherche

Instabilité d’un fluide dans un cylindre en précession

Instabilité fluide-élastique d’un faisceau de tubes

Instabilité de flambage / plissement

Perspectives de recherche

Hors programme

Instabilité de flottement

Instabilité de flambage sur substrat non-linéaire

Instabilité de plissement 3D

35 of 89

Perspectives de recherche

35

Instabilité fluide-élastique,

faisceau à pas triangulaire

Modélisation des forces fluides diphasiques en faisceau à pas triangulaire

Modélisation des forces turbulentes diphasiques (portance, trainée, large gamme de taux de vide), faisceau à pas carré

Recherche orientée vers l’industrie (Framatome, EDF, projets Européens: GO-VIKING, NARSIS+)

Spectre des forces turbulentes diphasiques, pas carré, portance

Expérience vs. Neptune-CFD. Thèse C. Bazin

De Langre

Axisa

Lagrange

(20 paramètres)

Spectre adimensionné des forces turbulentes diphasiques, faisceau à pas carré, portance

Expériences

Fréquence adimensionnée

Spectre adimensionné

Recherche amont appliquée

Développement et validation de l’IFS monophasique dans TrioCFD (code interne CEA) sur base de modèles simplifiés et de R&D

36 of 89

Perspectives de recherche

36

Thèse: quantification des incertitudes et analyse de sensibilité pour les vibrations de structures minces sous écoulement axial

37 of 89

Perspectives de recherche

37

IFS cylindres concentriques, vibrations imposées, fluide au repos

IFS, cylindres excentriques, vibrations imposées, fluide au repos

Pression adimensionnée

Masse ajoutée propre

Masse ajoutée croisée

38 of 89

Merci de votre attention

38

39 of 89

39

40 of 89

40

41 of 89

ANNEXES

41

42 of 89

42

IRPHE

Marseille

(Master, Thèse 2006-2009)

IFPEN

Lyon

(Post-doc 2010-2011)

MIT

Cambridge, USA

(Instructeur 2012-2015)

Mécanique des fluides

Mécanique des structures

Interactions fluides-structures

Instabilité de battement

(Master)

(C. Eloy, L. Schouveiler)

Instabilité de précession

(Thèse)

(P. Meunier, C. Eloy, F. Nadal)

Transition turbulence précession

(1 stage)

(1 stage)

IFS class project

Centrale

Marseille

(2003-2006)

Instabilité de plissement

(P. Reis, J.Dunkel)

Flambage sur substrat NL

(D. Averbuch)

CEA

Saclay

(Ingénieur-chercheur 2015-…)

Instabilité fluide-élastique en faisceau

(9 stages, 2 thèses en cours)

Instabilité précession hors RES

(P. Meunier, C. Eloy)

43 of 89

Dynamique non-linéaire

et chaos

Niveau

Nombre d’heures

Cours

UPMC

Master 2

24 (3)

MMC

UPMC

Licence / CPGE

19 (2)

Mathématiques

MIT

Undergraduate

42 (2)

CNAM

Tensor calculus

MIT

3

42

Travaux dirigés

Mathématiques

UPMC

ORSAY

Licence / CPGE

IUT

UPMC

Master 1

19 (2)

22

14

Expériences de l’enseignement

Differential equations

Multivariable calculus

MIT

MIT

Undergraduate

Undergraduate

52 (2)

52 (2)

Licence

Etablissement

Total HETD

83.5

2017-2018

182.5

2016-2017

44 of 89

Dossier personnel

44

Automne

Janvier

Printemps

2012-2013

Multivariable calculus

X

Multivariable calculus

2013-2014

Nonlinear dynamics and chaos

X

Differential equations

2014-2015

Nonlinear dynamics and chaos

Tensor calculus

(PDF disponibles)

Differential equations

Recitations : 4h / semaine + Office hours + grading

Lectures : 3h / semaine + Office hours + PSets hebdo. + Examens + grading

Enseigné en Master 2 à l’UPMC (2016-2019)

Enseignements au MIT (2012-2015)

45 of 89

Dossier personnel

45

Enseignement à l’ENPC (2022-…)

Séance

Chapitre

Test en ligne

1

Algèbre et analyse tensorielle

X

2

Cinématique et déformations (HPP)

X

3

Lois fondamentales et contraintes

X

4

Loi de comportement d’un matériau élastique linéaire homogène isotrope

X

5

Exercice 1 (réservoir sphérique sous pression, approche en déplacement)

6

Exercice 2 (poutre en flexion, approche en contrainte)

Mise à niveau en mécanique des milieux continus

(40-45 étudiants, 6 séances de 2h30)

Rédaction d’un manuscrit de cours (disponible ici) et quizz en ligne

46 of 89

Activités de recherche

46

47 of 89

47

48 of 89

Instabilité de précession

48

Précession de la Terre

Structure du noyau de la Terre

Noyau externe : fer liquide (80 %)

Conducteur électrique 🡪 Champ magnétique

49 of 89

Aeronautic application

49

Asymmetry 🡪 Precession

Liquid propellant 🡪 Fluid motion

Torque 🡪 angle amplification 🡪 change of trajectory

Flying object in precession

50 of 89

50

Linear amplitude

Amplitude of the first Kelvin mode

Good prediction of the structure and amplitude of the modes

At a resonance: viscous and weakly nonlinear theory

51 of 89

Axial vorticity of modes (PIV)

51

Different Kelvin modes observed depending on the frequency ratio

A Kelvin mode is made of several counter rotativive vortices

Mode 2

Mode 3

Mode 5

Mode 1

52 of 89

52

1st res.

2nd res.

Mode 3

Mode 1

Mode 2

Mode 5

Mode 1

Mode 2

Mode 3

Mode 5

Good prediction of the resonance frequencies

If , the mode fits the cylinder height: resonance,

Instabilité de précession

53 of 89

Instabilité de précession

53

Navier-Stokes

Modes de Kelvin forcés

+

+

+

1st res.

2nd res.

Mode 1

Mode 2

Solution Euler avec conditions aux limites, ordre

Cisaillement

McEwan (1970)

54 of 89

Instabilité de précession

54

Base flow

(forced mode)

Free Kelvin mode

(instability modes)

Geostrophic mode

Ecoulement total

Amplitude equations

Base flow linearly unstable for

55 of 89

55

Stable

Unstable

Fixed point

Unstable

Limit cycle

1st bif. : the flow becomes unstable. Transcritical bif.

2nd bif. : the fixed point becomes unstable. Supercritical Hopf bif.

Good prediction of the Hopf bif threshold and period of the limit cycle

Instabilité de précession

56 of 89

56

57 of 89

57

58 of 89

Instabilité de flambage / plissement

58

Uniaxial compression

Equi-biaxial compression

Allen (1969)

Slight initial curvature imperfection may play a role in mode selection

Audoly,

Boudaoud (2008)

59 of 89

2D substrate problem

59

60 of 89

60

2. Modeling

Solution ordre 0

No substrate

Film much stiffer than substrate

Same limits for

Instabilité de flambage / plissement

61 of 89

61

Curvature has no effect on the wavenumber

Discretization of the wavelength due to the periodic closing conditions of the ring

62 of 89

62

Curvature has no effect on the critical stress

63 of 89

63

3. Results

Comparison to plane case

Curvature has no effect on the wavelenth

64 of 89

64

3. Results

Phase diagram

Transition has the same scaling law for small inner radius

65 of 89

65

The aspect ratio determines the surface pattern

Rapport d’aspect,

Instabilité de flambage / plissement

Diagramme de bifurcation

Diagramme de phase

66 of 89

Instabilité de surface

avec P. Reis, J. Dunkel et al. (MIT)

Trainée

Depressur.

Mode hexagonal

Mode mixte

Mode labyrinthe

Rapport d’aspect

Stress

    • 3 nombres sans dimension:

augmente

Rapport de rigidité

67 of 89

Théorie vs. Expériences

    • Energie potentielle:

Min.

    • Swift-Hohenberg:
    • Stabilité des modes:
    • Equations d’amplitude:

SH

Mode hexagonal

Mode mixte

Mode labyrinthe

augmente

68 of 89

68

69 of 89

69

70 of 89

70

71 of 89

Dimensionless numbers (Single phase)

71

Five dimensionless numbers

Aspect ratio

Pitch ratio

Stokes number

Reynolds number

Mass ratio

Flexible tube in a rigid tube bundle, single phase flow

Nine parameters

Characteristic length

Characteristic time

Characteristic mass

72 of 89

Instabilité fluide-élastique

73 of 89

Instabilité fluide-élastique

73

3 sets of blades : L2C (24 Hz in water)

L4C (18 Hz in water)

L1C (11 Hz in water)

Blade L4C (18 Hz in water)

0,0051

0,0174

77,29

5,36

1,356

2,77

0,0278

0,1008

78,37

3,64

1,16

0,0295

0,1070

78,39

2,93

0,44

0,0313

0,1134

78,37

2,18

-0,30

0,0347

0,1261

78,42

0,79

-1,51

5 to 6 points per VF

VF in the range

74 of 89

Instabilité fluide-élastique

74

75 of 89

Instabilité fluide-élastique

Lines are regressions of the form

with A and B two fitting parameters

corresponding to the thresholds in single-phase flow

76 of 89

Instabilité fluide-élastique

76

Fitted law for the instability threshold

77 of 89

Instabilité fluide-élastique

77

9 independent parameters (13 for two-phase)

Tube

Fluid

Flow

Geometry

Single-phase cross flow

Lift direction

78 of 89

Instabilité fluide-élastique

78

For a given bundle and aspect ratio

The instability threshold is an iso-value of

Total damping

79 of 89

Instabilité fluide-élastique

79

Shaker

Instability threshold

Direct method (variation of the Stokes number)

New law for the instability threshold in single-phase flow

80 of 89

Instabilité fluide-élastique

80

CONCLUSIONS

New law for the instability threshold (single flexible tube in a square bundle, lift direction, air-water)

New law for the instability threshold

81 of 89

Instabilité fluide-élastique

81

Bubbly

Annular

Plug/slug

Churn

Map of the flow regimes

Instability threshold

82 of 89

Instabilité fluide-élastique

82

CONCLUSIONS

Total damping

Theory

Experiments

83 of 89

83

Conclusions

Precessing instability

Wrinkling instability

Curvature plays a role in the pattern selection

Fluid-elastic instability in GV

Potential theory not effective

Instability threshold (in single phase flow) depends on

Triadic resonance

84 of 89

84

2D potential model (no incident flow)

Set position of

to align cylinders with iso-contours

Instantaneous position of cylinders

Bipolar coordinates

Coordinates of a fluid particle

Update the coordinate system with the position of cylinders

Aspect ratio

Separation distance

85 of 89

85

Self-added mass

Fluid forces

Fluid forces

Cross-added mass

increases with confinement

86 of 89

86

Fluid force terms

|Added mass coefficients|

increases with confinement

Isolated cylinder limits 

Cylinders repelled from each other

D’Alembert’s paradox

87 of 89

87

Sinusoidal vibrations

and

colinear

88 of 89

Objectif 2028 : réaliser le couplage MANTA / TrioCFD pour disposer d’un code fiable, prédictif et qualifié pour l’IFS

Mieux comprendre la physique de l’IFS et améliorer sa modélisation (forces élastiques, stabilité…)

Disposer d’un code interne pour les études R&D et les calculs industriels (GV, RJH, assemblages combustibles…)

Axe 1

Développement et validation de l’IFS dans TrioCFD,

sur base de modèles simplifiés et de R&D

Axe 2

Développement et validation du couplage MANTA / TrioCFD

Modèles de turbulence et lois de paroi

Nouvelle cinématique maillage ALE vs. forts déplacements

Etudes paramétriques, nouvelles lois, validation CFD-théorie

Couplage interne fluide-structure

Réalisation du couplage MANTA / TrioCFD

Calcul cible (26) du plan de développement MANTA

(2026-2028)

Calculs de R&D en écoulement axial (assemblage, RJH,…) et transverses (faisceaux de tubes, …) avec

modèles mécaniques simples

Perspectives de recherche

89 of 89

89

Schéma de la cuve d’un REP. Thèse M. Kocher

Espace annulaire

Admission

Evacuation

Panier interne

Plenum inférieur

Enveloppe de coeur