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
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
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
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 |
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 |
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
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) |
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
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)
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
Instabilité de précession �(2006-2009, 2012-2015)
11
Principaux collaborateurs |
C. Eloy, P. Meunier |
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
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,
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é
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é
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
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 »
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
Instabilité de flambage / plissement�(2012-2015)
19
Principaux collaborateurs |
J. Dunkel, P. Reis, N. Stoop, D. Terwagne, F. Lopez Jimenez |
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
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
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
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
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
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 (???)
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
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
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
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
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
Instabilité fluide-élastique en faisceau
31
Phase 1
Phase 2
STABLE
INSTABLE
STABLE
INSTABLE
Diphasique avec
deux phases identiques
Instabilité fluide-élastique en faisceau
32
Diphasique avec
deux phases identiques
Diphasique
?
STABLE
STABLE
INSTABLE
INSTABLE
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
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
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
Perspectives de recherche
36
Thèse: quantification des incertitudes et analyse de sensibilité pour les vibrations de structures minces sous écoulement axial
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
Merci de votre attention
38
39
40
ANNEXES
41
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)
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
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)
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
Activités de recherche
46
47
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
Aeronautic application
49
Asymmetry 🡪 Precession
Liquid propellant 🡪 Fluid motion
Torque 🡪 angle amplification 🡪 change of trajectory
Flying object in precession
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
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
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
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)
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
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
57
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)
2D substrate problem
59
60
2. Modeling
Solution ordre 0
No substrate
Film much stiffer than substrate
Same limits for
Instabilité de flambage / plissement
61
Curvature has no effect on the wavenumber
Discretization of the wavelength due to the periodic closing conditions of the ring
62
Curvature has no effect on the critical stress
63
3. Results
Comparison to plane case
Curvature has no effect on the wavelenth
64
3. Results
Phase diagram
Transition has the same scaling law for small inner radius
65
The aspect ratio determines the surface pattern
Rapport d’aspect,
Instabilité de flambage / plissement
Diagramme de bifurcation
Diagramme de phase
Instabilité de surface
avec P. Reis, J. Dunkel et al. (MIT)
Trainée
Depressur.
Mode hexagonal
Mode mixte
Mode labyrinthe
Rapport d’aspect
Stress
augmente
Rapport de rigidité
Théorie vs. Expériences
Min.
SH
Mode hexagonal
Mode mixte
Mode labyrinthe
augmente
68
69
70
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
Instabilité fluide-élastique
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
Instabilité fluide-élastique
74
Instabilité fluide-élastique
Lines are regressions of the form
with A and B two fitting parameters
corresponding to the thresholds in single-phase flow
Instabilité fluide-élastique
76
Fitted law for the instability threshold
Instabilité fluide-élastique
77
9 independent parameters (13 for two-phase)
Tube
Fluid
Flow
Geometry
Single-phase cross flow
Lift direction
Instabilité fluide-élastique
78
For a given bundle and aspect ratio
The instability threshold is an iso-value of
Total damping
Instabilité fluide-élastique
79
Shaker
Instability threshold
Direct method (variation of the Stokes number)
New law for the instability threshold in single-phase flow
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
Instabilité fluide-élastique
81
Bubbly
Annular
Plug/slug
Churn
Map of the flow regimes
Instability threshold
Instabilité fluide-élastique
82
CONCLUSIONS
Total damping
Theory
Experiments
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
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
Self-added mass
Fluid forces
Fluid forces
Cross-added mass
increases with confinement
86
Fluid force terms
|Added mass coefficients|
increases with confinement
Isolated cylinder limits
Cylinders repelled from each other
D’Alembert’s paradox
87
Sinusoidal vibrations
and
colinear
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
Schéma de la cuve d’un REP. Thèse M. Kocher
Espace annulaire
Admission
Evacuation
Panier interne
Plenum inférieur
Enveloppe de coeur