Club CastemParis 20/11/2006 Calcul d´instabilités en interaction fluide-structure

Paris 20/11/2006Club Castem

Calcul d´instabilités

en

interaction fluide-structure


Plan de la présentation

1. Introduction: instabilités fluide-élastiques

2. Essais AMOVI

3. Modélisation

4. Résultats

5. Conclusion

2 5 4 3 1

Paris 20/11/2006Club Castem 2 5 4 3 1

Instabilités fluide-élastiques



2 5 4 3 1





1 5 4 3 2


•Effets inertiels :

masse ajoutée

•Effets visqueux :

amortissement ajouté

•Amortissement structurel


2 5 4 3 1


•Masse, amortissement et rigidité ajoutés

•Amortissement et rigidité de signe quelconque


2 5 4 3 1


Essais AMOVI

1 5 4 3 2


Essais AMOVI

•fair = 14.3 Hz

•D = 12.15 mm

•Un seul tube flexible

•Enregistrement du déplacement

•Analyse fréquentielle

Amortissement et fréquence en fonction de la vitesse

1 5 4 3 2


Modélisation : formalisme ALE

Déplacement du maillage.

n1 1

2 2

2n n

ts

Règle du trapèze

1 5 4 2 3


Modélisation : maillage

Maillage: environ 150 000 éléments quadratiques

Δx = 2. 10-4 m (pas moins d’un élt dans CL) x/D~10-2

Calcul 2D laminaire incompressible

1 5 4 2 3


Modélisation : équations

0

.

div u

uu grad u grad p u

t

.. .

.. .

( ) ( ) ( ) ( )

( ) ( ) ( ) 0

i i i

i i i i i i

s i s i s i i

s a i s a i s a i

q t q t q t t

q t q t q t

M C K F

M M C C K K

. .s

su

t

p n n n

Fluide:

Structure:

Interface:

1 5 4 2 3


Résultats

Fluide sans écoulement


Résultats: fluide sans écoulement

1 5 3 2 4


Résultats : fluide sans écoulement

1 5 3 2 4


Résultats

Fluide sous écoulement

Calcul de la maquette AMOVI


Résultats : calcul AMOVI sous écoulement

(a) Re = 1789 (b) Re = 2786

(c) Re = 5964

1 5 3 2 4


0,00

0,50

1,00

1,50

2,00

2,50

3,00

0,025 0,03 0,035 0,04 0,045 0,05 0,055 0,06 0,065

V (m/s)

Am

ort

issm

en

t ré

du

it

Amovi Calcul numérique

-1,50

0,50

2,50

4,50

6,50

8,50

10,50

12,50

0 0,02 0,04 0,06 0,08 0,1 0,12 0,14 0,16

V (m/s)F

réq

uen

ce f

Amovi Calcul numérique

Résultats : comparaison avec AMOVI

1 5 3 2 4


( 1)( ) ( )t tA A h Ch O h

( ) (a)

( / 2) (b)2

( / 4) (c)4

t

t

t

A h A Ch

hA h A C

hA h A C

Résultats : Étude de convergence temporelle (cas a)

Extrapolation de Richardson (1) :

1 5 3 2 4



A(h/4)-A(h/2)ln : ln(2)

A(h/2)-A(h)t

( / 2) 2 ( )

1 2

t

t

A h A hA

Extrapolation de Richardson (2) :

•Ordre :

•Valeur Convergée:

1 5 3 2 4



1 5 3 2 4

Ordre de ~ 1 ; Valeur Convergée : - 2 % d’erreur !!


Résultats : Analyse des champs de pressions (cas b)

1 5 3 2 4


Résultats : Analyse des champs de pression-vitesse (cas b)

1 5 3 2 4


Résultats

Fluide sous écoulement

Calcul à Re plus faible


Résultats : 0 < Re < 500 0 < Vr < 5

R

VV

fD

1 5 3 2 4


2,00

2,05

2,10

2,15

2,20

2,25

2,30

2,35

0,0025 0,005 0,0075 0,0113 0,015 0,0187 0,02 0,0225 0,025

V (m/s)

Fré

qu

ence

f-3

-2

-1

0

1

2

3

4

0,0025 0,005 0,0075 0,0113 0,015 0,0187 0,02 0,0225 0,025

V (m/s)A

mo

rtis

sm

en

t ξ

Résultats : 0 < Re < 500 0 < Vr < 5

1 5 3 2 4


Résultats : 0 < Re < 500 0 < Vr < 5

1 5 3 2 4


Conclusion et perspectives

• Calculs AMOVI :

Bons résultats à bas Reynolds (Richardson)

Problèmes à plus haut Re (régime écoulement ? Limites de l’approche 2D ?)

• Cas théorique :

Bons résultats qualitatifs, pas de comparaison possible

• Perspectives :

Compréhension de la physique (animations, étude paramétrique…). Comparaison PIV

1 4 3 2 5

Documents

Club CastemParis 20/11/2006 Calcul d´instabilités en interaction fluide-structure