PA05 SYSIF– METHODES
NUMERIQUES
Application numérique des équations de Barré de
Saint-Venant et Euler
ARAÚJO DA SILVA João Victor
PEREIRA DO NASCIMENTO, Lucas
08/11/2019 ARAUJO DA SILVA, João Victor
PEREIRA DO NASCIMENTO, Lucas
Table des matières
1. Introduction................................................................................................. 2
2. Barré de Saint-Venant .................................................................................3
2.1. La validation du code ...............................................................................3
2.2. Les modifications introduites ..................................................................3
2.3.1. Exploitation des résultats et discussion ............................................5
3. Euler............................................................................................................... 6
3.1. Modèles physiques .................................................................................. 6
3.1.1. Le Tube de Sod ..................................................................................... 6
3.1.2. L’interaction entre 4 ondes de choc planes ..................................... 6
3.2. Les modifications introduites : Tude de Sod ..........................................7
3.2.1. La validation du code ........................................................................ 8
3.3. Les modifications introduites : problème 2D ........................................ 11
3.3.1. Exploitation des résultats et discussion ............................................. 12
4. Conclusion................................................................................................... 13
PA 05 – Application numérique des équations de Barré de Saint-Venant et Euler 1
08/11/2019 ARAUJO DA SILVA, João Victor
PEREIRA DO NASCIMENTO, Lucas
1. Introduction
Dans le cadre de ce projet, il nous a été proposé d'étudier des différents cas liés à
l'écoulement des fluides. Dans un premier temps, nous nous appuierons sur l'écoulement
d'un film de liquide de faible épaisseur, dont la modélisation sera faite par les équations de
Barré de Saint-Venant.
L'objectif est de simuler numériquement la rupture d'un barrage circulaire avec un fond
plat en utilisant le schéma de Rusanov pour valider les calculs avec un article donnée.
Ensuite, il faut créer une nouvelle géométrie du fond et commenter son influence.
Dans un deuxième temps, nous modifierons l'approche pour pouvoir passer à un problème
concernant les équations d'Euler pour les gaz. Tout d'abord, nous nous intéressons à l'étude
d'un problème 1D, le tube de choc, plus spécifiquement celui de Sod, avec des conditions
initiales données.
Il faudra alors réadapter le codage pour obtenir des courbes d'intérêt (pression, masse
volumique, ...) et les analyser. A la fin, nous passerons à un problème 2D, toujours en
utilisant les équations d'Euler, pour étudier la structure finale résultant de l'interaction
entre 4 ondes de choc.
Pour ce faire, nous utiliserons le logiciel Python pour pouvoir modifier un code déjà existant
et l'adapter aux différents cas d'étude, aussi bien que le VisIt pour pouvoir obtenir des
visualisations plus riches et améliorer la compréhension du problème.
PA 05 – Application numérique des équations de Barré de Saint-Venant et Euler 2
08/11/2019 ARAUJO DA SILVA, João Victor
PEREIRA DO NASCIMENTO, Lucas
2. Barré de Saint-Venant
2.1. La validation du code
La validation de notre code sera faite par comparaison des résultats obtenus
numériquement selon notre code et ceux présentés par l'article de base.
Les géométries et les conditions initiales et de contour étant identiques, on analysera les
projections des 2 résultats sur le domaine à un temps spécifique (t = 3s) afin d'observer la
formation d'une vague en amont de l´écoulement du barrage. Nous avons aussi choisi
d’utiliser un CFL de 0.4.
Notre Article (Ricchiuto et al 2009)
Maillage 600x600 Maillage 1000x1000 Maillage 10113x10113
Tableau 1 - Comparaison entre les resultats obtenus
Les similarités des résultats peuvent être observées par les courbes obtenues. Les
différences des allures peuvent être attribuées à plusieurs facteurs comme, par exemple, la
différente méthode utilisée (éléments finis) et donc aussi au maillage du domaine.
[Link] modifications introduites
Les modifications ont été surtout au niveau du fond et du terme source. Par rapport au
premier, il a fallu définir une paramétrisation pour une géométrie désirée pour pouvoir
obtenir des valeurs du vecteur z, qui était une matrice carré nulle auparavant. Nous
obtenons, alors :
PA 05 – Application numérique des équations de Barré de Saint-Venant et Euler 3
08/11/2019 ARAUJO DA SILVA, João Victor
PEREIRA DO NASCIMENTO, Lucas
Figure 1 - Modifications introduites pour le fond
En plus, a été aussi nécessaire de définir ses dérivées partielles par rapport à x et y, pour
pouvoir obtenir le terme source. Dans ce cas, les modifications ont été les suivantes :
Figure 2 - Modifications au niveau du terme source
Nous pouvons aussi commenter la modification faite au niveau de la fonction d’obtention
des fichiers Tecplot, où nous avons sauvegardé la somme du fond + niveau d’eau dans la
variable_0, suivi du fond dans la variable_1, pour pouvoir visualiser les deux surfaces sur
VisIt. On a alors :
Figure 3 - Modification du codage pour affichage sur VisIt
Pour finir, on a imposé une condition de micro-pluie pour créer une couche d’eau permanente et
artificielle qui éviterait d’avoir des divisions par zéro au niveau des flux :
PA 05 – Application numérique des équations de Barré de Saint-Venant et Euler 4
08/11/2019 ARAUJO DA SILVA, João Victor
PEREIRA DO NASCIMENTO, Lucas
2.3.1. Exploitation des résultats et discussion
En ce qui concerne les tests réalisés, nous avons choisi d’utiliser un maillage de 350x350
avec 3600 iterations, pour pouvoir étudier les phénomènes liés à l’écoulement de l’eau
jusqu’à un moment proche de la stabilisation. A propos du fluide, la hauter initialle a été
𝑑𝑜𝑚(𝑥)
de 40m par rapport à la rampe dans une petite zone comprise entre 0 ≤ 𝑥 ≤ , pour
20
pouvoir bien voir l’effet de la pesanteur sur le système.
Par rapport au fond, il a été modélisé comme une rampe avec un obstacle sous forme de
paraboloïde elliptique inversé centré au milieu, avec une base dont le rayon vaut 25m. Ci-
dessous, nous montrons les effets de l’écoulement pour quelques iterations.
Iteration : 0 Iteration : 1600 Iteration : 3600
Tableau 2 - Evolution au cours du temps de l'écoulement de l'eau
Il est alors possible de vérifier que l’obstacle empêche bien l’écoulement de l’eau, en la
faisant s’étaler vers les côtés, vers l’eau (en suivant le format de l’obstacle) et vers l’arrière.
En plus, on voit bien que les conditions limites ont été bien définies, car l’eau né dépasse
pas les limites établies par le parallélépipède, au contraire, elle est réflechie.
Il est aussi possible de déduire que le fluide commence à se stabiliser, une fois que les
couches qui se forment au-dessus commencent à reduire son épaisseur, et le niveau d’eau
réfletie est aussi plus faible.
PA 05 – Application numérique des équations de Barré de Saint-Venant et Euler 5
08/11/2019 ARAUJO DA SILVA, João Victor
PEREIRA DO NASCIMENTO, Lucas
3. Euler
3.1. Modèles physiques
3.1.1. Le Tube de Sod
Tout d’abord, on s’intéressera à un cas spécifique du Tube de Choc, un problème 1D régit
par les équations d’Euler, avec des conditions initiales suivantes :
𝑢𝐿 = 0 𝑢𝑅 = 0
𝑃𝑜𝑢𝑟 𝑙𝑎 𝑚𝑜𝑖𝑡𝑖é 𝑔𝑎𝑢𝑐ℎ𝑒: {𝑝𝐿 = 1 𝑃𝑜𝑢𝑟 𝑙𝑎 𝑚𝑜𝑖𝑡𝑖é 𝑑𝑟𝑜𝑖𝑡𝑒: { 𝑝𝑅 = 0.1
𝜌𝐿 = 1 𝜌𝑅 = 0.125
De manière succincte, l’essai consiste à remplir un tube en métal avec un gaz dans deux
états différents (comprimé en amont et à basse pression en aval), séparés par une
membrane, qui sera prudemment déchirée, en générant une onde qui se propagera vers
l’aval.
Cette propagation créera des zones spécifiques qui feront varier les propriétés du gaz lors
de son homogénéisation, soit sous forme continue (dans la zone de détente), soit sous
forme de discontinuité (dans les zones de contact et de choc). C’est donc ce comportement
qui sera le cas de notre étude. Il est possible de mieux visualiser la physique du problème
dans la figure ci-dessous :
Figure 4 - Description du tube de Sod
3.1.2. L’interaction entre 4 ondes de choc planes
Dans un deuxième temps, nous nous intéresserons donc à un problème 2D pour étudier
le résultat de l’interaction entre 4 ondes de choc planes, toujours en ayant une zone de
détente qui se propagera de manière déphasée. Les conditions initiales sont alors :
PA 05 – Application numérique des équations de Barré de Saint-Venant et Euler 6
08/11/2019 ARAUJO DA SILVA, João Victor
PEREIRA DO NASCIMENTO, Lucas
𝑣1 = 0 𝑣2 = 0
𝑢 =0 𝑢2 = 1.206
Ω1 : 𝑦 > 0.8 𝑒𝑡 𝑥 > 0.8 → { 1 Ω2 : 𝑦 > 0.8 𝑒𝑡 𝑥 < 0.8 → {
𝑝1 = 1.5 𝑝2 = 0.3
𝜌1 = 1.5 𝜌2 = 0.5323
𝑣3 = 1.206 𝑣4 = 1.206
𝑢3 = 1.206 𝑢4 = 0
Ω3 : 𝑦 < 0.8 𝑒𝑡 𝑥 < 0.8 → { Ω4 : 𝑦 < 0.8 𝑒𝑡 𝑥 > 0.8 → {
𝑝3 = 0.029 𝑝4 = 0.3
𝜌3 = 0.138 𝜌4 = 0.5323
3.2. Les modifications introduites : Tude de Sod
Au niveau des modifications introduites, on a enlevé tout ce qui concernait le fond (donc
le terme source, les variables d’altitude etc). Ensuite, nous nous sommes appuyés sur les
équations d’Euler pour déclarer les vecteurs contenant les variables d’intérêt (la pression,
la masse volumique, les vitesses, l’Energie, le nombre de Mach), en s’utilisant du fait que le
gaz étudié peut être modélisé comme un gaz parfait, et donc les équations suivantes sont
valides :
𝛾×𝑃
𝑃 = (𝛾 − 1) × 𝜌 × 𝑒 , 𝑎2 =
𝜌
En plus, on a maintenant 3 valeurs propres, soit 𝑢 + 𝑎, 𝑢, 𝑢 − 𝑎. Dû à ces faits, nous avons
changé les expressions des vecteurs w, f, g et ceux des flux physiques comme suit :
Figure 5 - Calcul du vecteur w dans le temps initial
Figure 6 - Calcul du flux physique en x
PA 05 – Application numérique des équations de Barré de Saint-Venant et Euler 7
08/11/2019 ARAUJO DA SILVA, João Victor
PEREIRA DO NASCIMENTO, Lucas
Figure 7 - Calcul du flux numérique en x
A la fin, on a aussi actualisé le vecteur pression après chaque pas de temps comme suit :
Figure 8 - Actualisation de la pression
3.2.1. La validation du code
Pour pouvoir valider les résultats obtenus, on a décidé de comparer les courbes numériques
avec la résolution analytique. Pour ce faire, nous avons stipulé des temps différents pour
pouvoir caractériser son influence lors de la précision du schéma. Nous avons utilisé un
maillage de 498x3, en obtenant :
PA 05 – Application numérique des équations de Barré de Saint-Venant et Euler 8
08/11/2019 ARAUJO DA SILVA, João Victor
PEREIRA DO NASCIMENTO, Lucas
Résultat numérique
t = 0.1 s t = 0.2 s
Résultat analytique
t = 0.1s t = 0.2s
Tableau 3 - Comparaison entre les résultats numériques et analytiques
Avec le but d’étudier les résultats, on a calculé les erreurs, pour 3 variables d’intérêt
(pression, masse volumique et la vitesse u). La formule utilisé a été la suivante :
𝐸𝑟𝑟𝑒𝑢𝑟_𝑣𝑎𝑟 = (𝑣𝑎𝑙𝑒𝑢𝑟_𝑎𝑛𝑎𝑙𝑦𝑡𝑖𝑞𝑢𝑒_𝑣𝑎𝑟 − 𝑣𝑎𝑙𝑒𝑢𝑟_𝑛𝑢𝑚é𝑟𝑖𝑞𝑢𝑒_𝑣𝑎𝑟) × 100
PA 05 – Application numérique des équations de Barré de Saint-Venant et Euler 9
08/11/2019 ARAUJO DA SILVA, João Victor
PEREIRA DO NASCIMENTO, Lucas
Erreurs
t = 0.1 s
𝑒𝑟𝑟𝑒𝑢𝑟_ max_𝐩 = 20.7% 𝑒𝑟𝑟𝑒𝑢𝑟_ max_𝐫𝐡𝐨 = 15.1% 𝑒𝑟𝑟𝑒𝑢𝑟_ max_𝐮 = 91.8%
t = 0.1 s
𝑒𝑟𝑟𝑒𝑢𝑟_ max_𝐩 = 31% 𝑒𝑟𝑟𝑒𝑢𝑟_ max_𝐫𝐡𝐨 = 23.3% 𝑒𝑟𝑟𝑒𝑢𝑟_ max_𝐮 = 92.7%
Tableau 4 - Calcul et répresentation des erreurs
On voit que les erreurs absolues, dans la majorité des cas, ont été inférieures à 31% (et donc
l’erreur relative doit être de l’ordre de 5% ou inférieur), ce qui montre la précision de
l’approche numérique effectuée.
Cependant, en ce qui concerne la vitesse, on voit une erreur très faible sauf dans une petite
zone pour les deux cas, et cela est dû au fait qu’il y a un décalage des courbes numériques
et analytiques dû à la discontinuité, ce qui fait que la courbe analytique devienne nulle
avant que la courbe numérique, et donc l’erreur prend la valeur de numérique entièrement
dans ces zones.
En plus, on peut aussi voir, bien définies, les 5 zones qui se forment après le choc au niveau
de la pression et de la masse volumique :
PA 05 – Application numérique des équations de Barré de Saint-Venant et Euler 10
08/11/2019 ARAUJO DA SILVA, João Victor
PEREIRA DO NASCIMENTO, Lucas
Propriétés après choc
Pression Masse volumique Vitesse u
3.3. Les modifications introduites : Le problème 2D
En passant du tube de choc au problème 2D d’interaction d’ondes planes, il nous a fallu
changer la déclaration des conditions initiales pour bien délimiter les 4 zones formées ainsi
que décommenter l’appel des fonctions de flux physique et numérique et la variable sigmay
dans le code principal. Ces fonctions étaient déjà prêtes puisque le code a été conçu pour
des applications 2D depuis le début. Les modifications des conditions initiales sont les
suivantes :
Figure 9 - Déclaration des nouvelles conditions limites
PA 05 – Application numérique des équations de Barré de Saint-Venant et Euler 11
08/11/2019 ARAUJO DA SILVA, João Victor
PEREIRA DO NASCIMENTO, Lucas
3.3.1. Exploitation des résultats et discussion
Pour étudier ce phénomène, on a utilisé un maillage de 1200x1200 avec 5000 itérations. Le
résultat obtenu a été le suivant :
Figure 11 - Résultat numérique obtenu par notre code Figure 10 - Résultat de l'article Liang Pan et al., 2016
On constate alors que l’allure commence à se former de manière satisfaisante, cependant le
nombre d’itérations n’a pas été suffisant pour que le phénomène ait lieu de manière complète.
En plus, l’article de Pan a utilisé un schéma d’ordre 2 (avec un maillage de 1000x1000).
Dû à ce fait, on a essayé d’augmenter le maillage pour pouvoir approximer le résultat (comme
l’erreur est élevée au carré d’une ordre à l’outre, on a fait le même avec le maillage, mais il
fallait utiliser 106 × 106 ). Malheureusement, par des raisons d’infrastructure, il n’a pas été
possible d’effectuer des calculs plus précis, puisque nos ordinateurs ne subissaient pas, et le
temps de calcul serait trop élevée.
PA 05 – Application numérique des équations de Barré de Saint-Venant et Euler 12
08/11/2019 ARAUJO DA SILVA, João Victor
PEREIRA DO NASCIMENTO, Lucas
4. Conclusion
Dans le cadre du projet réalisé jusqu’ici, nous avons pu nous familiariser avec les différents
aspects d’implémentation des deux ensembles d’équations (Euler et Saint-Venant). La
compréhension et l’utilisation de l’équation d’advection linéaires, les différents types de
schémas et flux numériques et leurs impacts sur la précision, la stabilité et les coûts de
calcul de nos codes.
Pour la première partie, avec la modification proposée, l’on a pu évaluer les influences de
la géométrie choisie pour le fond sur le terme source. Nous avons aussi observé le besoin
de l’insertion d’une couche d’eau artificielle constante sur tout le fond afin d’éviter une
erreur de division par zéro, vue que la variable zéro “wc[h,i,j]” est le dénominateur des
plusieurs d’autres expressions.
Deuxièmement, après avoir réalisé les modifications nécessaires pour évaluer le tube de Sod
dans une approche d’Euler et testé le programme pour le deux cas (1D et 2D), l’utilisation
d’une discrétisation de premier ordre s’est avérée assez couteuse en termes de calcul à cause
du nombre important de points requis afin d’obtenir des erreurs plus faibles par rapport
aux résultats de Pan. La mise en œuvre d’un schéma d’ordre supérieur pourrait donc réduire
considérablement les couts de calculs et conséquemment le temps total demandé.
PA 05 – Application numérique des équations de Barré de Saint-Venant et Euler 13