Ce papier décrit l'application de la méthode de réduction a priori (APR) dans plusieurs cas académiques. Cet algorithme, qui consiste a déterminer de façon a priori une base d'un écoulement, utilise la décomposition de Karhunen-Loève et des sous-espaces de Krylov, et a été appliquée dans le cas de l'équation de convection-diffusion 2D, ainsi que dans le cas de la cavité entraînée. Couplée à un système dynamique, cette méthode présente l'avantage d'être extrémement rapide et permet un réajustement dynamique de la base lors de modifications de paramètres contrôlant l'écoulement

ALLERY, Cyrille

HAMDOUNI, Aziz

VERDON, Nicolas

I-Revues

18 ème Congrès Français de Mécanique Grenoble, 27-31 août 2007
Résolution numérique des équations de transfert par une méthode de
réduction a priori
Nicolas Verdon, Cyrille Allery & Aziz Hamdouni
Université de La Rochelle
Laboratoire d’Etude des Phénomènes de Transfert Appliqués au Bâtiment
Avenue Michel Crépeau F-17042 La Rochelle cedex 1
nverdon@univ-lr.fr
Résumé :
Ce papier décrit l’application de la méthode de réduction a priori (APR) à plusieurs cas académiques. Cet algo-
rithme, qui consiste a déterminer sans échantillonage préalable une base d’un écoulement, utilise la décomposi-
tion de Karhunen-Loève et des sous-espaces de Krylov, et a été appliquée dans le cas de l’équation de convection-
diffusion 2D, ainsi que dans le cas de la cavité entraînée. Couplée à un système dynamique, cette méthode présente
l’avantage d’être extrémement rapide et permet un réajustement dynamique de la base lors de modifications de
paramètres contrôlant l’écoulement.
Abstract :
This paper describes the application of the A Priori Reduction method (APR) in several academic test cases. This
algorithm, which consists in determining a basis of a flow without prerequisite knowledge, uses the Karhunen-
Loève Decomposition and Krylov subspaces, and has been been applied to the convection-diffusion equation and
also for the driven cavity configuration. Coupled with a dynamical system, this method has the important ad-
vantage to be computationally very fast and allows a dynamical update of the basis when control parameters are
changing.
Mots-clefs :
Modèles d’Ordre Réduit ; Navier-Sokes ; Système Dynamique
1 Introduction
La résolution numérique des équations de transfert a une importance grandissante du point
de vue fondamental mais aussi à un niveau industriel. Cependant, certaines méthodes requièrent
à l’heure actuel un effort de calcul trop important pour prétendre à des applications indus-
trielles. Par exemple, le contrôle actif nécessite des temps de calcul très courts pour permettre
des adaptations du système mis en jeu en temps réel. Dans ce contexte on trouve bon nombre
de méthodes de réduction de modèles. Dans le cas de la mécanique des fluides, qui présente des
applications variées, la méthode intensivement utilisée est la décomposition orthogonale aux
valeurs propres (POD) qui possède l’avantage de construire une base optimale du point de vue
énergétique. Plus de propriétés et de détails théoriques sur cette méthode peuvent être trouvés
dans Holmes et al. (1996). Associée à un système dynamique obtenu par une méthode Galer-
kin, cette approche permet de représenter l’évolution d’un écoulement sur des temps bien plus
longs que ceux nécessaires à la construction de la dite base. Cependant, il est connu que le sys-
tème dynamique n’est valable que pour un certain temps, voir Sirisup et Karniadakis (2004),
avant que la solution ne diverge. De même, un changement dans les paramètres de contrôle de
l’écoulement, à savoir le nombre de Reynolds ou la modification d’une condition aux limites,
n’est pas bien représenté par le système dynamique. La base calculée n’est effectivement va-
1
18 ème Congrès Français de Mécanique Grenoble, 27-31 août 2007
lable que pour une gamme étroite de ces paramètres.
Pour s’affranchir des problèmes liés à la POD, nous proposons une méthode itérative de construc-
tion de la base. Cette approche présente l’avantage de réajuster la base lorsqu’elle n’est plus
valable (temps trop longs, changement de paramètres de contrôle, etc...).
Le papier s’organise de la façon suivante : dans un premier temps la méthode est décrite pour
une équation de transfert quelconque et dans un deuxième temps les résultats obtenus pour la
convection-diffusion ainsi que pour les équations de Navier-Stokes dans la configuration de la
cavité entraînée seront présentés.
2 Présentation de la méthode APR
La méthode APR étant itérative, il convient d’initialiser le processus de résolution, c’est-à-
dire ici de s’imposer une base initiale. Comme il a déjà été mis en évidence par Ryckelynck
(2006), le choix de cette base n’est pas un élément très important pour la qualité de la solution,
mais détermine la rapidité et le nombre d’itérations que l’algorithme devra effectuer. En effet,
choisir une base plus physique requiert moins d’itérations de calcul APR pour faire converger
la solution. On appelle cette base initiale ψ0. La méthode étant expliquée en détails dans Verdon
et al. (2006), seules ses grandes lignes seront expliquées dans les paragraphes suivants.
2.1 Construction du système réduit
On considère donc le problème physique 1 que l’on doit résoudre, que l’on note :
L(u, p) = b (1)
La discrétisation de l’équation (1) conduit, après linéarisation si nécessaire, au système d’équa-
tions linéaires suivant :
Aζ = B (2)
où A est la matrice issue de la discrétisation de l’équation (1), de taille N , ζ est le vecteur
inconnues et enfin B le second membre, incluant les conditions aux limites et les données de la
solution calculée au pas de temps antérieur.
Notons maintenant n le nombre de vecteurs de base de ψ (n ≪ N). Le système réduit associé
au modèle complet (2) est obtenu en écrivant l’inconnue ζ sur la base ψ
ζ(x, t) =
n∑
i=1
ψi(x)ai(t) (3)
et en projetant le système obtenu sur cette même base. Le système réduit dont les nouvelles
inconnues sont les coefficients temporels a(t) s’écrit de la façon suivante :
F(ai(t)) = 0 ∀t ∈ T ∀i = 1, · · · , n (4)
où T est l’intervalle de temps sur lequel on désire construire la base de la solution. Il faut donc
résoudre ce système réduit pour tous les pas de temps de l’intervalle T . La solution complète
est alors reconstruite en utilisant l’expression de sa décomposition (3) sur la base ψ. La qualité
de la solution est alors évaluée sur tout l’intervalle en considérant le résidu du modèle complet
r = Aζ − B. Si la solution n’est pas satisfaisante, la base est alors adaptée comme expliqué
dans le paragraphe suivant.
1équation de convection-diffusion, équations de Navier-Stokes, etc...
2
18 ème Congrès Français de Mécanique Grenoble, 27-31 août 2007
2.2 Adaptation et enrichissement de la base
Lors de l’itération courante, on a résolu l’équation réduite sur tout l’intervalle de temps,
c’est-à-dire que l’on a déterminé les coefficients temporels a(t), ∀t ∈ T . Pour la phase d’amé-
lioration de la base, on réalise une décomposition de Karhunen-Loève sur cet ensemble. On est
donc ammené à résoudre le problème aux valeurs propres suivant :
CΦ = λΦ (5)
où C est la matrice de corrélation définie de la façon suivante :
Cij =
Nt∑
l=1
ai(tl)aj(tl) (6)
avec Nt le nombre de pas de temps. On obtient ainsi n vecteurs propres correspondant aux
n valeurs propres solutions de cette équation. On considère alors les η plus grandes valeurs
propres, définies telles que λη > µsel · λ1 2 afin de ne garder que les modes propres dont la
contributions dans l’écoulement est significative de la dynamique observée sur l’intervalle de
temps. La matrice d’amélioration V est alors formée de la façon suivante : ses vecteurs colonnes
sont les η vecteurs propres conservés, c’est-à-dire V = [Φ1 Φ2 · · ·Φη]. La nouvelle base est
alors améliorée à l’aide de cette matriceV :
ψ˜ = ψ · V (7)
Ensuite, il faut également tenir compte de l’erreur que l’on a comise à l’itération k, c’est pour-
quoi on conserve le premier résidu du modèle complet r = Aζ−B qui ne satisfait pas le critère
de convergence et on construit l’espace de Krylov3 d’ordre m associé à ce résidu :
Km = {r,Ar,A
2r, · · · ,Am−1r} (8)
A la base améliorée ψ˜, on rajoute les vecteurs du sous-espace de Krylov que l’on vient de cal-
culer, c’est ce qu’on appelle la phase d’enrichissement de la base, c’est-à-dire le moment où
l’on prend en considération les erreurs de l’itération précédente. Dans la suite, on prend m = 1.
La base ainsi mise à jour, c’est-à-dire améliorée par la décomposition de Karhunen-Loève et
enrichie par les vecteurs de Krylov, est utilisée pour l’itération suivante. Si elle permet d’obte-
nir une solution suffisamment précise, elle sera conservée par la suite pour la construction du
système dynamique. Dans le cas contraire, elle sera mise à jour jusqu’à ce que la qualité de la
solution calculée soit jugée satisfaisante.
3 Système dynamique
Une fois que l’on a obtenu une base qui représente correctement la solution du problème
considéré sur l’intervalle de temps étudié, on espère simuler la solution sur un intervalle de
temps beaucoup plus long que celui nécessaire au calcul de la base. Pour ce faire, on construit
un système dynamique en effectuant une projection de Galerkin de l’équation sur la base de
l’écoulement.
En utilisant l’expression (3) de la décomposition de la vitesse sur la base ψ, et par projection,
on obtient le système dynamique suivant :
dai
dt
=
n∑
j=1
Bijaj +
n∑
j=1
n∑
k=1
Cijkajak +Di pour i = 1, · · · , n (9)
2où les valeurs propres sont rangées dans l’ordre décroissant et µsel est le coefficient de sélection que l’on se fixe
3Pour plus de détails sur les méthodes de Krylov, le lecteur pourra se référer à Saad (1996)
3
18 ème Congrès Français de Mécanique Grenoble, 27-31 août 2007
où l’on a pour l’équation de convection-diffusion :
Bij = (ψn,−αψm + ν∆ψm) Cijk = 0 et Di = 0 (10)
et pour les équations de Navier-Stokes :
Bij = (ψn, ν∆ψm) Cijk = −(ψn,∇ψm) et Di =
∫
∂Ω
P · ψnn dS (11)
La résolution du système d’équations différentielles (9) par un algorithme de Runge-Kutta est
très rapide en comparaison avec la résolution des équations pendant la précédente phase. La
solution complète est alors calculée à l’aide de (3) sur des temps bien plus long que le temps
nécessaire pour la méthode APR.
4 Résultats
Cette partie présente les deux cas tests que nous avons étudiés, à savoir l’équation de convec-
tion diffusion et les équations de Navier-Stokes pour le cas d’une cavité entraînée.
4.1 Cas 1 : équation de convection-diffusion
L’équation de convection-diffusion a été étudiée dans le domaine carré 0 ≤ x, y ≤ 2 avec la
connaissance de la solution analytique donnée dans Tian et Ge (2007)4 :
φ(x, y, t) =
1
4t+ 1
exp
[
−
(x− αt− 0.5)2
Γ(4t+ 1)
−
(y − αt− 0.5)2
Γ(4t+ 1)
]
(12)
Les conditions aux limites de type Dirichlet et les conditions initiales ont été obtenues à partir
de cette solution analytique. Par ailleurs, la base initiale de notre algorithme correspond aux
conditions initiales du problème.
Plusieurs expérimentations numériques ont été effectuées sur cet exemple.
Dans un premier temps, nous avons testé la faculté de réadaptation de la base APR lorsque le
système dynamique diverge, et dans une deuxième approche, on a étudié l’adaptation rapide de
la base pour un changement de paramètre de contrôle, à savoir lorsque Γ varie.
Dans tout ce qui suit, l’intervalle de calcul est noté T et correspond à une durée T = 0.0125s.
Dans un premier temps les calculs on été menés sur un intervalle T à partir de t = 0s afin de
déterminer une première base APR. A partir de cette base, comportant 5 vecteurs, on construit
un système dynamique que l’on résout sur 10T . On constate alors qu’à partir de t = 6T , la base
APR ne représente plus correctement la dynamique du système et l’erreur calculée par rapport
à la solution analytique dépasse les 5%5 pour atteindre 22% à t = 10T . On effectue alors un
nouveau calcul APR sur un intervalle T à partir de t = 0.075s et avec la nouvelle base obtenue,
on construit le nouveau système dynamique correspondant qui est résolu jusqu’à t = 10T c’est-
à-dire l’intervalle que l’on voulait étudier. L’erreur à la fin du calcul est alors aux alentours de
5% (voir fig. 1), ce qui est fortement amélioré par rapport aux performances de la première base
et le temps de calcul nécessaire a été divisé par deux par rapport au modèle complet.
Dans une deuxième approche, on a testé la sensibilité de la base à un changement de para-
mètre, ici Γ. On repart de la première base calculée sur l’intervalle T et on construit le système
dynamique correspondant, mais cette fois-ci avec Γ = 0.02 alors que l’itération d’APR a été ef-
fectuée pour Γ = 0.01. On se rend alors compte que même sur l’intervalle T où a été calculée la
4avec α = 0.8 et Γ = 0.01
5on juge arbitrairement que la base n’est plus valable lorsque l’erreur dépasse 5%
4
18 ème Congrès Français de Mécanique Grenoble, 27-31 août 2007
(a) Solution analytique (b) Modèle complet (c) Modèle réduit + SD
FIG. 1 – Résultats du cas 1 : isolignes de φ pour t = 0.0125s et t = 0.125s
base, le système dynamique n’est pas capable de reproduire fidèlement la solution puisque l’on
trouve une erreur égale à 22.16% pour t = T . L’idée est alors de relancer la méthode APR, en
considérant comme base initiale la première base calculée. En seulement 3 itérations APR, on
arrive à une base, comprenant 4 vecteurs, qui représente correctement la solution, c’est-à-dire
que le système dynamique formé à partir de cette base ne donne plus qu’une erreur d’environ
1.9% à t = 0.0125s.
Par ailleurs, la modification de la base n’a nécessité que 6s de temps de calcul alors que la
construction de la première base en a requis environ 20s, ce qui présente un intérêt en vue
d’applications en contrôle actif.
4.2 Cas 2 : équations de Navier-Stokes
Le cas étudié ici est la cavité entraînée unitaire (0 ≤ x, y ≤ 1) représentée sur la figure 2
pour un nombre de Reynolds Re = 1 avec, initialement, le fluide au repos. Ce cas a été traîté
u = 0
v = 0
u = 1
v = 0
u = 0v = 0
u = 0
v = 0
FIG. 2 – Configuration de la cavité entraînée 2D
sans le système dynamique mais nous a tout de même permis de valider la méthode APR pour
les équations de Navier-Stokes. On constate que la méthode APR permet de retrouver la solu-
tion stationnaire (fig. 3) pour un nombre de Reynolds valantRe = 1 en environ deux fois moins
de temps que le modèle complet, c’est-à-dire 413.3s de calcul contre 781, 8s pour le modèle
complet, et une erreur relativement faible, comme présenté sur la figure 3(c).
5
18 ème Congrès Français de Mécanique Grenoble, 27-31 août 2007
(a) Modèle complet (b) Modèle réduit
0
0.2
0.4
0.6
0.8
1
0
0.2
0.4
0.6
0.8
1
−0.01
−0.005
0
0.005
0.01
0.015
0.02
0.025
0.03
 
xy
 
Er
re
ur
(c) Erreur entre modèle réduit et modèle complet
FIG. 3 – Résultats du cas 2 : isolignes de la vitesse pourRe = 1
Cette étude est une étude préliminaire à un cas instationnaire qui sera présenté lors de la com-
munication orale.
5 Conclusions
Dans ce travail, la méthode APR a été appliquée avec succès à deux configurations : l’équa-
tion de convection-diffusion instationnaire et les équations de Navier-Stokes dans le cas de la
cavité entraînée à Re = 1. Les performances qui sont ici mises en avant sont d’une part, le
gain de temps comparé à un calcul avec un modèle complet, puis une facilité et une rapidité
d’adaptation de la base APR, d’une part lorsque le système dynamique diverge, et d’autre part
lorsqu’un paramètre de contrôle change. La prochaine étape du travail consiste donc à appliquer
cette méthode à des écoulements turbulents.
Références
Holmes P., Lumley J.L., and Berkooz G., Turbulence, Coherent Structures, Dynamical Systems
and Symmetry, Cambridge University Press, Cambridge, ISBN 0 521 55142 0, 1996.
Saad Y., Iterative Methods for Sparse Linear Systems, 2nd Edition, 2000.
Ryckelynck D., A priori hyperreduction method : an adaptive approach, Journal of Computa-
tional Physics, Volume 202, Issue 1, pp. 346-366, 1er janvier 2005.
Verdon N., Allery C., Ryckelynck D., Hamdouni A., An adaptive ROM approach for solving
transfer equations, Revue européenne de mécanique numérique, Volume 15, N. 5, 2006
Sirisup S., Karniadakis G. E., A spectral viscosity method for correcting the long term behavior
of POD models, Journal of Computational Physics, Volume 194, p. 92-116, 2004
Tian Z. F., Ge Y. B. 2007 A fourth-order compact ADI method for solving two-dimensional
unsteady convection-diffusion problems Journal of Computational and Applied Mathematics
Volume 198, pp. 268-286
6


Résolution numérique des équations de transfert par une méthode de réduction a priori

http://documents.irevues.inist.fr/bitstream/2042/16505/1/CFM2007-0797.pdf

Résolution numérique des équations de transfert par une méthode de réduction a priori

Abstract

Similar works

Full text

Available Versions

I-Revues