Adaptation de la POD aux

domaines fluides mobiles :

application à la réduction de
modèle par POD en IFS
LIBERGE Erwan
[email protected]

Directeur de thèse : A. Hamdouni

Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 1/18
Introduction
– Intérêt de la réduction de modèles par la
méthode POD
– Système dynamique de faible dimension
– Etude de stabilité
– Contrôle actif
– Temps de calcul réduit et diminution du stockage de données
– Développement croissant en turbulence depuis
une décénnie
– Initié en turbulence par Lumley(1967)
– En mécanique des structure la POD est un
domaine récent d’investigation
– Trindade (2005), Sarkar et Paidousis (2004), Epureanu (2004)
– Cela débute en interaction fluide structure
Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 2/18
Plan
– Introduction
– Décomposition orthogonale aux valeurs propres
– problème aux valeurs propres associé à la POD
– Propriétés des fonctions issues de la POD
– La snaphot POD
– Application aux domaines mobiles
– Problématique
– Introduction d’un domaine de référence
– Test de la base POD
– Application à la réduction de modèle
Test sur un cas monodimensionel

– Conclusion
Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 3/18
Décomposition orthogonale aux
valeurs propres (POD)

Analyse en composantes principales

–
(ACP) en analyse de données

– Introduite en mécanique des fluides par Lumley (1967)

→identification des structures cohérentes
– Construction d’une base "physique" optimale au sens
énergétique
– Essor fin des années 80
nombreuses applications : canal, cavité entraînée, disque
contra-rotatif, couche de mélange, jet ...
(Sirovich et al., Adaptation
Cazemier ...)domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 4/18
de la POD aux
Problème aux valeurs propres
associé à la POD
– Soit u (x, t) le champ de vitesses (supposé aléatoire) à
valeurs réelles, on cherche une fonction φ (déterministe),
dans H = L2 (Ω), telle que
 D E
 max (u, Φ)2
Φ∈H (1)
 (Φ, Φ) = 1

(•, •) correspond au produit scalaire de H et h•iest un

opérateur de moyenne statistique.
– Trouver λ ∈ R et φ ∈ L2 (Ω) tels que
Z
R (x, y) φ (y) = λφ (x) (2)
Ω

où R (x, y) = hu (x, t) ⊗ u (y, t)i est le tenseur de corrélations

spatial Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 5/18
Propriétés des fonctions issues
de la POD
– λn ∈ R+ λ1 > λ2 > . . . > λn−1 > λn

– (φn , φm ) = δnm , φn vérifient les CL et div (φn ) = 0

– Toute réalisation du champ u (x, t) peut s’écrire :

+∞
X
u (x, t) = an (t) φn (x) au sens de L2 (Ω)
n=1

où an (t) = (u (•, t) , φn )
N
X N
X
– Energie : EN = han (t) an (t)i = λn
n=1 n=1

– Meilleure approximation au sens énergétique

Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 6/18
La Snapshot POD (Sirovitch,
1987)
M : nombre de réalisation de l’écoulement
Nm le nombre de noeud
nc Le nombre de composantes
M << Nm nc On recherche les coefficient temporels a sous la
forme :
M
X
Φ (x) = a (tk ) u (x, tk )
k=1

Ce qui revient à résoudre le problème aux valeurs propres

suivant :
M
X 1
(u (ti ) , u (tk )) a (tk ) = λa (ti ) pour i = 1 . . . M
k=1
M
Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 7/18
Application aux domaines mo-
biles

– uΓI = VΓI
– FF = σF .n
avec n normale extérieure

Problématique
– Base spatiale alors que les domaines sont mobiles.

– Difficulté à connaître le domaine relatif à chaque modes

Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 8/18
Application aux domaines mo-
biles
Calcul de la matrice Snapshot :

C = [cij ]
1
où cij = M
(u (ti , •) , u (tj , •))

Problème :
Definition du produit scalaire (•, •) pour le calcul de C ?
les domaine fluides à l’instant i et j ne sont pas identiques.

Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 8/18
Application aux domaines mo-
biles
Solutions envisagées :

– Avoir une évolution des modes propres au cours du temps

(base adaptative ?)

– Utiliser un domaine de référence englobant toutes les

positions du domaines mobiles.

Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 8/18
Utilisation d’un domaine de ré-
férence
– Domaine de référence fixe Ω comportant une discrétisation
fixe au cours du temps.

Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 9/18
Utilisation d’un domaine de ré-
férence
– Domaine de référence fixe Ω comportant une discrétisation
fixe au cours du temps.
– Utilisation d’une fonction caractéristique du domaine
mobile 
 1 si ~x ∈ Ω (t)
F
χF (~x, t) = (3)
 0 si ~x ∈ ΩS (t)

Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 9/18
Utilisation d’un domaine de ré-
férence
– Domaine de référence fixe Ω comportant une discrétisation
fixe au cours du temps.
– Utilisation d’une fonction caractéristique du domaine
mobile 
 1 si ~x ∈ Ω (t)
F
χF (~x, t) = (4)
 0 si ~x ∈ ΩS (t)

– Interpolation des champs du domaines mobile (vitesse,

pression,. . .) sur le maillage fixe.

 v (x, t) if x ∈ Ω (t)
f F
v (x, t) = (5)
 vs (x, t) if x ∈ ΩS (t)
Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 9/18
Interpolation sur le maillage de
référence
On considère le problème 2D d’un solide rigide immergé
dans un fluide incompressible

ΩS domaine solide
ΩF domaine fluide
ΓI Interface fluide solide

Fluide incompressible de Reynolds 200. La simulation à été

effectuée sous Castem pour une durée de 5.55s

Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 10/18
Interpolation sur le maillage de
référence

discrétisation initial
Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 10/18
Interpolation sur le maillage de
référence

Noeuds composants le
maillage initial

Discrétisation du domaine
de référence

Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 10/18
Interpolation sur le maillage de
référence

Superposition des deux discrétisations

Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 10/18
Interpolation sur le maillage de
référence
Le noeud N ∈ au maillage
de reférence

ri :distance entre le noeud

i du maillage mobile et le
noeud N

Interpolation inverse (Shepard method, 1968)

U1 U2 U3 U4
r1
+ r2
+ r3
+ r4
UN = 1
r1
+ r12 + 1
r3
+ 1
r4

Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 10/18
Fonction caratéristique
On définit la fonction caratéristique « par maille » de la
manière suivante :

 1 si les 4 noeuds ∈ ΩF (t)


χF (i, t) = Taux d’occupation de la maille par le domaine fluide


0 si les 4 noeuds ∈ ΩS (t)

(6)

Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 11/18
Fonction caratéristique

Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 11/18
Test de la base POD
Efficacité de la base POD à retranscrire le phénomène testé
sur le champ de vitesse total interpolé (sur
Ω = ΩF (t) ∪ ΩS (t))
Avec 5 modes on capture 99.999% de l’énergie totale
mode i eigenvalue λi % energy
1 254 99.837
2 0.343 99.9722
3 5.553e − 2 99.994
4 1.100e − 2 99.998
5 2.674e − 3 99.999

Contribution énergétique des 5 premiers modes

Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 12/18
Test de la base POD
La reconstruction du champ de vitesse sur la base POD
tronquée à Nmodes modes est évaluée en norme L2

an (t) = (u (•, t) , Φn )
NX
modes
U (x, t) = an (t) Φn (x)
n=1

kU − UN modes kL∞ (0,T ;L2 (Ω))

Res (N modes) = (7)
kU kL∞ (0,T ;L2 (Ω))

Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 12/18
Test de la base POD

Erreur de reconstruction inférieure à 1% à partir de 5 modes

Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 12/18
Application à la réduction du
système dynamique
Formulation pour une structure rigide en ALE


∂u
ρ + ρ(u − w).∇u = −∇p + µ4u


 ∂t


∇·u = 0 dans ΩF (t)
dV



 H
 + [V, HV ] = Fs + FF pour ΩS (t)
dt
avec w la vitesse du maillage
Formulation variationnelle
V ∗ et v ∗ compatibles
   
∂vf ∗ ∗ dV
ρ + ρvf .∇vf , v − (µ4vf , v ) + H ,V ∗
∂t dt
+ h[V, HV ] , V ∗ i = hFs , V ∗ i
Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 13/18
Application à la réduction du
système dynamique
Formulation variationnelle
V ∗ et v ∗ compatibles
   
∂vf ∗ ∗ dV
ρ + ρvf .∇vf , v − (µ4vf , v ) + H ,V ∗
∂t dt
+ h[V, HV ] , V ∗ i = hFs , V ∗ i

V ∗ et v ∗ compatibles + action réaction

Z
hFf , V ∗ i + σ · n · v ∗ dγ = 0
ΓI

Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 13/18
Système réduit
En choissisant v ∗ = Φm
Pour n = 1..N :


N X N N Z
 dan X X ∂vs
+ Bnmp am ap + Anm am =


 .Φn dx
dt Ω ∂t
m=1 p=1 m=1
 H dV + [V, HV ] = Fs + FF in ΩS (t)



dt
où
Z Z
Anm = µ ∇Φm .∇Φn XΩF dx + 2µ (∇Φm .Φn ) dΓ
Z Ω ΓI

Bnmp = Φm .∇Φp .Φn XΩF dx

Ω

Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 14/18
Application à un cas monodi-
mensionel

F IG . 1 – schématisation du domaine d’étude

Le fluide est modélisé par l’équation de Burgers résolue en
ALE couplée à une équation de ressort

Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 15/18
Contribution énergétique des
modes
mode i eigenvalue λi % energy
1 18.255 91.648
2 1.515 99.256
3 8.393e − 2 99.677
4 6.079e − 2 99.982
5 2.347e − 3 99.994
6 7.804e − 4 99.999

Valeurs propres et contribution énergétique des modes POD

calculés pour le champ de vitesse sur le maillage de référence
Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 16/18
Contribution énergétique des
modes

Modes POD 1 et 2 obtenus sur le maillage de référence

Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 16/18
Performance du système dyna-
mique
N
X
On injecte vN (x, t) = an (t) Φn (x) dans le système
n=1
dynamique

Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 17/18
Performance du système dyna-
mique

Résidu issue de la POD directe et du système réduit

Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 17/18
Performance du système dyna-
mique

Solution numérique et solution du système réduit à différents

instants
Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 17/18
Performance du système dyna-
mique

F IG . 2 – temporal coefficients a1 ,a2 ,a3 ,a4

Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 17/18
Conclusions
– La méthode de réduction de modèle par la POD
semble permettre une bonne reproduction pour
peu de modes.

– Perspectives :
– Finir le système réduit 2D
– Recherche d’une base a priori (Chinesta,
Ryckelynck,Verdon)

Adaptation de la POD aux domaines fluides mobiles : application à la réduction de modèle par POD en IFS – p. 18/18