PROJET VOLUME FINI :

Étude du Phénomène
DAM BREAK

-2023/2024-
Page 0 sur 43
 Rapport du Projet CFD :
Étude du Phénomène
de Rupture de Barrage
« DAM BREAK »

▪ Réalisé par :
BENSBIK HAITAM
ELBAROUDI NASSIM
ELJANDOUZ BRAHIM
ERRAMY IBRAHIM
OUFETTOUL FARID
▪ Encadré par :
Pr. AMAL ABOULHASSANE

Page 1 sur 43
 DAM BREAK

Sommaire

INTRODUCTION ---------------------------------------------------------- 4
Description Physique et Mathématique ----------------------- 5
Domaine de Calcul et Type d'Écoulement ------------------------------------------- 5
Mise en Équations du Problème Physique --------------------------------------------- 5
Modélisation Numérique --------------------------------------------------------------------- 5
Modèle de Surface Libre ------------------------------------------------------------------ 5
Paramètres de Simulation ---------------------------------------------------------------- 5
Validation et Calibration du Modèle ------------------------------------------------- 5

Définition des Propriétés Physiques ------------------------------ 7

Propriétés du Fluide ------------------------------------------------------------------------- 7
Autres Propriétés et Paramètres -------------------------------------------------------- 7
Justification des Propriétés --------------------------------------------------------------- 7

DEFINITION DES CONDITIONS AUX LIMITES ET DES

CONDITIONS INITIALES --------------------------------------------------- 9
Conditions aux Limites --------------------------------------------------------------------- 9
Conditions Initiales --------------------------------------------------------------------------- 9
Simulation de la Rupture ------------------------------------------------------------------ 9

Justification du Choix du Solveur ------------------------------- 11

Adaptation aux Écoulements Multiphasiques ---------------------------------- 11
Méthode Volume of Fluid (VoF) ------------------------------------------------------ 11
Méthode des Volumes Finis ------------------------------------------------------------ 11
Flexibilité et Personnalisation----------------------------------------------------------- 11

Simulation numérique du problème -------------------------- 12

Simulation du volume sans obstacle -------------------------------------------------- 12
a- Prétraitement : ----------------------------------------------------------------------- 12
b- Traitement : ---------------------------------------------------------------------------- 15
c- Post-traitement : --------------------------------------------------------------------- 20
Simulation du volume avec l'obstacle de double taille ------------------------ 27

Page 2 sur 43
 DAM BREAK

a- Prétraitement : ----------------------------------------------------------------------- 27
b- Traitement ----------------------------------------------------------------------------- 31
c- Post-traitement ----------------------------------------------------------------------- 31

ANNEXE ------------------------------------------------------------------ 41
La méthode VoF (Volume of Fluid) Vs la méthode des volumes finis -- 41

bibliographie ---------------------------------------------------------- 42

Page 3 sur 43
 DAM BREAK

INTRODUCTION
La rupture de barrage est un phénomène hydraulique critique pouvant
entraîner des conséquences catastrophiques pour les communautés et les
infrastructures situées en aval. Lorsqu'un barrage cède, l'eau retenue est
libérée de manière soudaine et violente, provoquant des inondations rapides
et de haute vélocité. Ces événements peuvent entraîner des pertes
humaines et des dégâts matériels importants, rendant impératif une
compréhension approfondie et une modélisation précise de ces
phénomènes pour prévenir et atténuer les risques associés.
La modélisation numérique par dynamique des fluides (CFD) offre un outil
puissant pour simuler les scénarios de rupture de barrage et analyser les
comportements des écoulements complexes résultants. En utilisant des
solveurs avancés et des techniques de simulation, il est possible de reproduire
les conditions réelles et d'étudier les différents facteurs influençant la
propagation des ondes de rupture.
Ce projet vise à modéliser la rupture de barrage en utilisant les logiciels
OpenFOAM et ParaView. L'objectif est de simuler le processus de rupture,
analyser les résultats obtenus et formuler des recommandations basées sur
ces observations. La méthode de volumes finis sera utilisée pour résoudre les
équations de Navier-Stokes, garantissant ainsi une précision et une stabilité
dans les simulations.
Les enjeux de la rupture de barrage ne se limitent pas seulement à la
gestion immédiate des inondations, mais incluent également la planification
à long terme et la mise en place de stratégies d'atténuation des risques. Une
modélisation précise permet non seulement de comprendre les dynamiques
de l'écoulement, mais aussi de concevoir des infrastructures plus résilientes et
de planifier des interventions d'urgence plus efficaces.
En conclusion, l'intérêt de ce sujet réside dans sa capacité à fournir des
informations cruciales pour la prévention des catastrophes hydrauliques et
l'amélioration de la sécurité des ouvrages. La modélisation CFD de la rupture
de barrage permet de simuler des scénarios complexes, d'anticiper les
conséquences potentielles et de proposer des solutions pratiques pour
minimiser les impacts sur les populations et les infrastructures. Ce projet
contribue ainsi à une meilleure compréhension des phénomènes en jeu et à
une gestion plus efficace des risques associés.

Page 4 sur 43
 DAM BREAK

DESCRIPTION PHYSIQUE ET MATHEMATIQUE

Domaine de Calcul et Type d'Écoulement
Pour simplifier la simulation de la rupture de barrage, nous utiliserons une
modélisation en deux dimensions (2D). Le domaine de calcul inclut le
réservoir, la zone de rupture et le lit de la rivière en aval. La modélisation 2D
permet de réduire la complexité et les ressources de calcul nécessaires tout
en capturant les principales dynamiques de l'écoulement.
L'écoulement est considéré comme isotherme, incompressible et en
régime laminaire. Ce type d'écoulement permet une modélisation plus
précise des interactions entre le fluide et les structures environnantes.

Mise en Équations du Problème Physique

Les équations de Navier-Stokes pour un fluide incompressible en régime
laminaire sont utilisées pour modéliser l’écoulement :
∂u 1
+ (u ⋅ ∇)u = − ∇p + ν∇2 u + g
∂t ρ
∇⋅u=0

Où 𝐮 est le vecteur vitesse, 𝑝 est la pression, 𝜌 est la densité du fluide, 𝜈 est

la viscosité cinématique, et 𝐠 est le vecteur gravité.

Modélisation Numérique
Modèle de Surface Libre
La méthode de Volume of Fluid (VoF) est utilisée pour suivre la surface libre
entre l'eau et l'air dans le modèle. Cette méthode est efficace pour capturer
les interfaces de phase, essentielle dans la simulation des phénomènes de
rupture de barrage.

Paramètres de Simulation
Les paramètres de simulation incluent la définition de la géométrie du
domaine, la création du maillage, et la configuration des conditions initiales
et aux limites. Les conditions aux limites typiques incluent des conditions de
non-glissement sur les parois et des conditions de pression à l'entrée et à la
sortie du domaine. Les conditions initiales définissent un champ de vitesse nul
et une distribution de pression hydrostatique dans le réservoir avant la rupture.

Validation et Calibration du Modèle

La validation du modèle 2D est réalisée en comparant les résultats de
simulation avec des données expérimentales ou des études de cas
documentées. Cette étape est cruciale pour ajuster les paramètres du
modèle et garantir la fiabilité des prédictions. Les études antérieures ont
démontré que les modèles 2D peuvent capturer efficacement les

Page 5 sur 43
 DAM BREAK

dynamiques principales des inondations de rupture de barrage tout en étant

moins exigeants en termes de ressources de calcul.

Page 6 sur 43
 DAM BREAK

DEFINITION DES PROPRIETES PHYSIQUES

Pour simuler la rupture de barrage avec précision, il est essentiel de définir
correctement les propriétés physiques du fluide utilisé dans la modélisation.
Ces propriétés influencent directement le comportement de l'écoulement et
les résultats de la simulation.

Propriétés du Fluide
Dans ce projet, nous considérons que le fluide est de l'eau douce, et les
propriétés physiques suivantes sont définies :
• Densité (ρ): 1000 kg/m³
o La densité de l'eau est une propriété cruciale, car elle affecte la
force de flottabilité et les mouvements du fluide sous l'effet de la
gravité.
• Viscosité dynamique (μ): 0.001 Pa.s (ou 1 x 10⁻³ kg/m·s)
o La viscosité dynamique est la mesure de la résistance du fluide à
la déformation sous une contrainte. Pour l'eau, cette valeur est
standard et affecte la manière dont le fluide se dissipe et se
mélange pendant l'écoulement.
• Viscosité cinématique (ν): 1 x 10⁻⁶ m²/s
o La viscosité cinématique est obtenue en divisant la viscosité
dynamique par la densité du fluide. Elle joue un rôle crucial dans
la diffusion de la vitesse dans le fluide.

Autres Propriétés et Paramètres

• Gravité (g) : 9.81 m/s²
o La gravité est une force motrice importante dans la simulation
des écoulements, surtout lors de la rupture d'un barrage où l'eau
est libérée subitement.
• Tension de Surface : En raison de la méthode VoF utilisée pour capturer
la surface libre, la tension de surface entre l'eau et l'air peut également
être un paramètre pertinent à considérer. Cependant, pour simplifier
cette simulation, nous pouvons négliger cet effet ou utiliser une valeur
standard si nécessaire.

Justification des Propriétés

Les valeurs choisies pour la densité, la viscosité dynamique et la viscosité
cinématique sont standard pour l'eau douce à température ambiante. Ces
propriétés permettent de modéliser de manière réaliste l'écoulement de l'eau

Page 7 sur 43
 DAM BREAK

après la rupture du barrage. La prise en compte de ces propriétés physiques

assure une simulation précise et fiable des dynamiques d'écoulement

Page 8 sur 43
 DAM BREAK

DEFINITION DES CONDITIONS AUX LIMITES ET DES

CONDITIONS INITIALES
Conditions aux Limites
Les conditions aux limites définissent les comportements du fluide aux
frontières du domaine de calcul. Voici les principales conditions aux limites
utilisées dans la simulation de rupture de barrage :

▪ Entrée
o Avant la Rupture : Condition de vitesse nulle (paroi fixe).
o Après la Rupture : Condition de débit basé sur la hauteur d'eau et la
vitesse de l'écoulement post-rupture. Cela peut être modélisé par
une condition de vitesse imposée ou une condition de flux sortant.

▪ Sortie
o Condition de Pression : Une condition de pression nulle (pression
atmosphérique) est appliquée pour permettre au fluide de s'écouler
librement hors du domaine. Cela représente un environnement
ouvert où l'eau peut s'échapper sans restriction.

▪ Parois
o Condition de Non-Glissement : Vitesse nulle sur les parois solides
(mur), ce qui signifie que le fluide adhère aux parois sans
mouvement relatif par rapport à celles-ci. Cela est important pour
modéliser l'interaction entre l'eau et les structures solides
environnantes.

Conditions Initiales
Les conditions initiales définissent l'état du fluide au début de la simulation :

▪ Champ de Vitesse Initial

▪ Vitesse Nulle : Au début de la simulation, le champ de vitesse dans
tout le domaine est initialisé à zéro. Cela signifie que le fluide est au
repos avant le début de la rupture du barrage.

▪ Champ de Pression Initial

Distribution Hydrostatique de la Pression : La pression initiale est définie en
fonction de la profondeur de l'eau dans le réservoir, suivant une distribution
hydrostatique : 𝑝 = 𝜌𝑔ℎ où p est la pression, 𝜌 est la densité de l'eau, g est
l'accélération due à la gravité, et h est la hauteur de la colonne d'eau au-
dessus du point considéré.

Simulation de la Rupture
Lorsque la rupture est simulée, les conditions aux limites et initiales
changent pour refléter le nouveau scénario dynamique :

Page 9 sur 43
 DAM BREAK

• La condition de débit à l'entrée est ajustée pour simuler l'eau s'écoulant

rapidement à travers la brèche.
• Les conditions de vitesse et de pression sont mises à jour en temps réel
en fonction des équations de Navier-Stokes et des conditions limites
définies.

Page 10 sur 43
 DAM BREAK

JUSTIFICATION DU CHOIX DU SOLVEUR

Pour ce projet de simulation de rupture de barrage, nous avons choisi le
solveur interFoam dans OpenFOAM. Voici les principales raisons justifiant ce
choix :

Adaptation aux Écoulements Multiphasiques

interFoam est conçu pour simuler les écoulements multiphasiques, où
plusieurs phases de fluides (comme l'eau et l'air) interagissent. Cette capacité
est cruciale pour modéliser la rupture de barrage, où l'interaction entre l'eau
libérée et l'air environnant joue un rôle significatif dans la dynamique de
l'écoulement.

Méthode Volume of Fluid (VoF)

Le solveur utilise la méthode Volume of Fluid (VoF) pour capturer avec
précision la surface libre entre les fluides. Cette méthode est essentielle pour
suivre l'interface entre l'eau et l'air pendant et après la rupture du barrage,
assurant une simulation précise des phénomènes de surface libre.

Méthode des Volumes Finis

interFoam utilise la méthode des volumes finis pour la discrétisation des
équations de Navier-Stokes. Cette approche est bien adaptée aux
simulations de dynamique des fluides, car elle permet de conserver les
propriétés de masse, de quantité de mouvement et d'énergie à travers les
cellules du maillage. Cela garantit une solution stable et précise même dans
les situations de flux complexes comme les ruptures de barrage.

Flexibilité et Personnalisation
OpenFOAM offre une grande flexibilité et permet de personnaliser les
simulations en fonction des besoins spécifiques du projet. Les utilisateurs
peuvent ajuster les paramètres du solveur, modifier les conditions aux limites
et initiales, et adapter les modèles de turbulence pour obtenir des résultats
optimisés et précis. Cette flexibilité est particulièrement bénéfique pour les
projets de recherche et de simulation avancée.

Page 11 sur 43
 DAM BREAK

SIMULATION NUMERIQUE DU PROBLEME

Le dossier de toutes les simulations, ICI.

Simulation du volume sans obstacle

Le dossier de la simulation, ICI.
Pour le maillage raffiné, ICI.

a- Prétraitement :
▪ Création de la géométrie et du maillage :
Le fichier blockMeshDict est utilisé pour définir la géométrie et le maillage de
base pour la simulation de la rupture de digue. Voici une explication détaillée
de la géométrie créée dans ce fichier :
DIMENSIONS DE LA GEOMETRIE :

La géométrie est définie par un ensemble de points appelés « vertices ».

Dans ce cas, les points définissent un volume en trois dimensions avec une
profondeur de la géométrie très faible, simulant ainsi un problème
bidimensionnel (2D).
- Les coordonnées des sommets (vertices) sont définies dans un système
de coordonnées cartésiennes :
vertices
(
(0 0 0)
(1.168 0 0)
(0 1.168 0)
(1.168 1.168 0)
(0 0 0.1)
(1.168 0 0.1)
(0 1.168 0.1)
(1.168 1.168 0.1)
);

➢ Les points (0 0 0), (1.168 0 0), (0 0 0.1) et (1.168 0 0.1)

définissent la base inf. de la géométrie.
➢ Les points (0 1.168 0), (1.168 1.168 0), (0 1.168 0.1) et (1.168
1.168 0.1) définissent la base sup de la géométrie.

DEFINITION DES BLOCS ET MAILLAGE

• Blocs de maillage :
o Le maillage est défini par un ensemble de blocs hexagonaux,
chacun spécifié par les indices des sommets :

Page 12 sur 43
 DAM BREAK

blocks
(
hex (0 1 3 2 4 5 7 6) (50 50 1) simpleGrading (1 1 1)
);

o hex (0 1 3 2 4 5 7 6) : Les indices des sommets qui définissent le

bloc.
o (50 50 1) : La résolution du maillage dans les directions x, y et z.
Ici, le maillage est divisé en 50 cellules dans les directions x et y,
et en 1 cellule dans la direction z.
o simpleGrading (1 1 1) : Indique un maillage uniforme dans toutes
les directions.

DEFINITION DES FRONTIERES

• Conditions aux limites :

o Les frontières (boundaries) sont définies pour indiquer les types de
conditions aux limites appliquées aux différentes faces du
domaine de calcul :

boundary
(
leftWall
{
type wall;
faces
(
(0 4 6 2)
);
}
rightWall
{
type wall;
faces
(
(3 7 5 1)
);
}
lowerWall
{
type wall;
faces
(
(0 1 5 4)
);
}
atmosphere
{
type patch;
faces

Page 13 sur 43
 DAM BREAK

(
(2 6 7 3)
);
}
);

o leftWall, rightWall, lowerWall : Sont définies comme des murs

(wall) indiquant des conditions aux limites de type mur.
o atmosphere : Est défini comme une condition aux limites de type
patch, indiquant une interface ouverte avec l'atmosphère.

DEFINITION DU VOLUME D'EAU INITIAL

Le fichier setFieldsDict est utilisé pour définir les conditions initiales du

champ de volume fractionnel d'eau dans le domaine de calcul. Voici une
explication détaillée du contenu de ce fichier :

• Contenu du fichier setFieldsDict :

1. Valeurs par défaut du champ :
o Le champ par défaut [Link] est initialisé à 0. Cela signifie
que, par défaut, il n'y a pas d'eau dans le domaine sauf dans les
régions spécifiquement définies.
defaultFieldValues
(
volScalarFieldValue [Link] 0
);

2. Définition des régions contenant de l'eau :

o Une région est définie où la fraction volumique d'eau
[Link] est initialisée à 10. Cela correspond au volume initial
d'eau dans le domaine.
o La région est définie par un volume rectangulaire spécifié par
deux coins opposés de la boîte.
regions
(
boxToCell
{
box (0 0 -1) (0.292 0.584 1);
fieldValues
(
volScalarFieldValue [Link] 10
);
}
);

o Explication de la région d'eau

- Volume de la boîte :
La boîte est définie par deux points (0 0 -1) et (0.292 0.584 1).

Page 14 sur 43
 DAM BREAK

▪ (0 0 -1) : Coin inférieur gauche de la boîte.

▪ (0.292 0.584 1) : Coin supérieur droit de la boîte.
Ces points définissent une boîte rectangulaire en 3D qui contient
l'eau initiale.
- Fraction volumique de l'eau :
Dans cette région définie par la boîte, la fraction volumique de
l'eau [Link] est définie à 10. Cela signifie que cette région
est entièrement remplie d'eau au début de la simulation.

▪ Raffinement du maillage :
Le raffinement du maillage est réalisé en modifiant le fichier
blockMeshDict, en doublant la résolution du maillage dans les zones
souhaitées, puis en relançant la simulation.
blocks
(
hex (0 1 3 2 4 5 7 6) (200 200 1) simpleGrading (1 1 1)
);

Le maillage est finement divisé en 200 cellules dans les directions x et y, et

en une seule cellule dans la direction z.

▪ Visualisation des résultats du prétraitement sur

Paraview :
Ici, vous trouver la vidéo de la simulation pour le premier maillage.
Ici, vous trouver la simulation du maillage raffiné.

b- Traitement :
▪ Définir les paramètres de contrôle de la simulation :
PARAMETRES DE TEMPS :

− startTime : Temps initial de la simulation.

− endTime : Temps final de la simulation.
− deltaT : Intervalle de temps pour chaque pas de temps.
startTime 0;
endTime 1;
deltaT 0.001;

FREQUENCE DE SAUVEGARDE DES RESULTATS

− writeControl : Critère de sauvegarde des résultats (par exemple,

basé sur les pas de temps).
− writeInterval : Fréquence à laquelle les résultats sont sauvegardés.
writeControl adjustable;
writeInterval 0.05;

Page 15 sur 43
 DAM BREAK

▪ Les schémas de discrétisation et l’algorithme de

résolution adoptés
SCHEMAS DE DISCRETISATION (FVSCHEMES) :

o Schémas de dérivée temporelle (ddtSchemes) :

ddtSchemes
{
default Euler;
}

• Euler : Utilisation du schéma d'Euler implicite pour la dérivée temporelle.

Ce schéma est simple et robuste, adapté pour les simulations avec des
pas de temps petits.
o Schémas de gradient (gradSchemes) :
gradSchemes
{
default Gauss linear;
}

• Gauss linear : Utilisation du schéma de Gauss avec une interpolation

linéaire pour calculer les gradients. Cela assure une précision de
second ordre pour les gradients, ce qui est important pour capturer les
variations de vitesse et de pression.
o Schémas de divergence (divSchemes) :
divSchemes
{
div(rhoPhi,U) Gauss linearUpwind grad(U);
div(phi,alpha) Gauss vanLeer;
div(phirb,alpha) Gauss linear;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

• div(phi,U) Gauss upwind : Utilisation du schéma upwind pour la

divergence du flux de vitesse. Ce schéma est stable et diffuse les
discontinuités, ce qui est utile pour capturer les chocs dans les
écoulements.
• div(phi,alpha) Gauss vanLeer : Utilisation du schéma de van Leer pour
la fraction volumique de l'interface. Ce schéma est plus précis pour les
écoulements multiphasiques, conservant mieux les interfaces nettes.
• div(phirb,alpha) Gauss interfaceCompression : Utilisation du schéma de
compression d'interface pour maintenir une interface nette entre les
phases.
• div((muEff*dev(T(grad(U))))) Gauss linear : Utilisation du schéma linéaire
pour les termes de diffusion visqueuse.

Page 16 sur 43
 DAM BREAK

o Schémas de Laplacien (laplacianSchemes) :

laplacianSchemes
{
default Gauss linear corrected;
}

• Gauss linear corrected : Utilisation du schéma linéaire corrigé pour les

termes de Laplacien. Ce schéma corrige les erreurs de non-
orthogonalité dans le maillage, améliorant la précision.
o Schémas d'interpolation (interpolationSchemes) :
interpolationSchemes
{
default linear;
}

• linear : Utilisation de l'interpolation linéaire pour les valeurs aux faces

des cellules. Cela assure une transition douce entre les valeurs des
cellules adjacentes.
o Schémas de dérivée normale de surface (snGradSchemes) :
snGradSchemes
{
default corrected;
}

• corrected : Utilisation du schéma corrigé pour les dérivées normales aux

faces, ce qui améliore la précision pour les maillages non orthogonaux.
ALGORITHMES DE RESOLUTION (FVSOLUTION) :

o Solveurs (solvers) :
solvers
{
"[Link].*"
{
nAlphaCorr 2;
nAlphaSubCycles 1;
cAlpha 1;

MULESCorr yes;
nLimiterIter 5;

solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-8;
relTol 0;

Page 17 sur 43
 DAM BREAK

"pcorr.*"
{
solver PCG;
preconditioner DIC;
tolerance 1e-5;
relTol 0;
}

p_rgh
{
solver PCG;
preconditioner DIC;
tolerance 1e-07;
relTol 0.05;
}

p_rghFinal
{
$p_rgh;
relTol 0;
}

U
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-06;
relTol 0;
}
}

Solveur [Link].* :

Le solveur [Link].* utilise plusieurs paramètres spécifiques pour la

méthode VOF (Volume of Fluid) :
❖ nAlphaCorr : Le nombre de corrections pour le champ alpha est fixé
à 2, ce qui permet d'améliorer la précision du calcul de l'interface.
❖ nAlphaSubCycles : Un sous-cycle est utilisé (nAlphaSubCycles 1),
permettant des itérations supplémentaires au sein d'un même pas
de temps.
❖ cAlpha : Coefficient alpha utilisé dans le modèle, ici fixé à 1.
❖ MULESCorr : La correction MULES (Multidimensional Universal Limiter
for Explicit Solution) est activée pour mieux gérer l'interface.

Page 18 sur 43
 DAM BREAK

❖ nLimiterIter : Fixé à 5, ce qui limite le nombre d'itérations de

correction pour stabiliser l'interface.
❖ Solveur : smoothSolver avec le lisseur symGaussSeidel, choisi pour sa
capacité à fournir une convergence douce et stable.
❖ Tolérance : 1e-8, une tolérance très stricte assurant une grande
précision.
❖ Tolérance relative : 0, exigeant une stricte convergence absolue.

Solveur pcorr.* :

Le solveur PCG (Preconditioned Conjugate Gradient) est utilisé pour les

corrections de pression :
❖ Préconditionneur DIC : (Diagonal Incomplete Cholesky) améliore la
convergence du solveur PCG en préconditionnant le système
d'équations.
❖ Tolérance : 1e-5, équilibrant précision et performance.
❖ Tolérance relative : 0, ce qui signifie que la convergence est basée
uniquement sur la tolérance absolue.

Solveur p_rgh :

Le solveur PCG est également utilisé pour la pression modifiée p_rgh :

❖ Préconditionneur DIC : Améliore la convergence en réduisant le

conditionnement du système.
❖ Tolérance : 1e-07, très stricte pour assurer une haute précision.
❖ Tolérance relative : 0.05, permettant un peu de flexibilité pour
accélérer la convergence tout en maintenant une bonne précision.

Solveur p_rghFinal :

Le solveur pour p_rghFinal réutilise les paramètres de p_rgh :

❖ Tolérance relative : 0, exigeant une stricte convergence absolue

pour la solution finale.

Solveur U :

Le solveur pour la vitesse U utilise smoothSolver :

❖ Lisseur symGaussSeidel : Fournit une convergence plus douce et

stable.
❖ Tolérance : 1e-06, assurant une bonne précision sans être
excessivement stricte.
❖ Tolérance relative : 0, exigeant une stricte convergence absolue.

Page 19 sur 43
 DAM BREAK

c- Post-traitement :
▪ Analyse Visuelle des Frames

Voici les principales observations pour les cinq images extraites de la vidéo
(maillage initial, raffiné) de votre simulation de rupture de barrage :

Frame 1 :
o Le point de rupture initial est visible.
o L'eau commence à se déplacer rapidement en aval, créant une
onde de rupture marquée.

Frame 2 :
o L'onde de rupture s'est propagée davantage.
o Les variations de vitesse et de pression sont observables dans la
zone de déplacement rapide de l'eau.

Frame 3 :
o La propagation continue de l'onde.
o Interactions possibles avec le terrain ou les obstacles visibles.

Frame 4 :
o Le flux d'eau se stabilise légèrement après l'onde de rupture
initiale.
o Variations de direction et de vitesse du flux visibles.

Frame 5 :
o L'eau atteint un nouvel équilibre après la propagation de l'onde.
o Les effets résiduels de la rupture sont encore observables.

▪ Champ de pression (ICI)

ANALYSE

• Début de la simulation :
Au début, la pression est élevée dans la région où se trouvait initialement
l'eau, en amont du barrage. Ceci est dû au poids de l'eau qui exerce une
pression hydrostatique sur le fond et les parois du réservoir.
• Milieu de la simulation :
Lorsque le barrage cède, l'eau commence à se déplacer rapidement vers
l'aval. On observe une diminution de la pression dans la région initiale du
réservoir et une augmentation de la pression en aval, où l'eau s'accumule en
rencontrant des obstacles ou le fond du canal.

Page 20 sur 43
 DAM BREAK

La dynamique des vagues de pression peut être observée à travers les

variations de couleurs, indiquant les gradients de pression causés par
l'écoulement turbulent de l'eau.
• Fin de la simulation :
À la fin, l'eau s'est répandue davantage en aval. La pression se stabilise
progressivement, mais des zones de haute pression peuvent persister dans les
régions où l'eau s'accumule ou rencontre des obstacles.
La pression est plus uniformément répartie comparée aux phases initiales,
mais avec des variations locales dues aux résidus de mouvement et de
turbulence de l'eau.
INTERPRETATION

• Distribution de la pression :
La distribution de la pression dans cette simulation montre clairement
l'effet de la rupture de barrage. La chute initiale de pression en amont et
l'augmentation en aval sont caractéristiques de ce type d'événement.
• Écoulement et turbulence :
Les variations rapides et les gradients de pression indiquent un écoulement
très turbulent, typique des ruptures de barrage où l'eau se déplace
rapidement et de manière désordonnée.
• Impacts environnementaux :
Ces résultats peuvent être utilisés pour évaluer l'impact environnemental
d'une rupture de barrage, notamment sur les structures en aval, l'érosion des
sols et la sécurité des infrastructures humaines.
COMPARAISON ENTRE LE MAILLAGE INITIAL ET LE MAILLAGE RAFFINE [VIDEO]

• Début de la simulation :
o Maillage Initial : La pression est concentrée en amont, avec des
gradients de pression moins distincts.
o Maillage Raffiné : La distribution de la pression est plus détaillée
avec des gradients plus précis. Les variations subtiles de la
pression sont mieux capturées.
• Milieu de la simulation :
o Maillage Initial : Les vagues de pression et les gradients sont
visibles mais moins définis. La turbulence est présente mais moins
détaillée.
o Maillage Raffiné : Les vagues de pression et les gradients sont
beaucoup plus nets. La turbulence est mieux représentée, avec
des détails fins dans les variations de pression.

Page 21 sur 43
 DAM BREAK

• Fin de la simulation :
o Maillage Initial : La pression commence à se stabiliser avec une
distribution relativement uniforme.
o Maillage Raffiné : La stabilisation de la pression est plus détaillée.
Les zones de haute pression résiduelle sont plus clairement
définies, et les variations locales sont mieux capturées.
Interprétation des Résultats
• Précision de la Simulation :
o Maillage Initial : Fournit une vue d'ensemble du phénomène mais
manque de précision dans les détails. Utile pour une analyse
rapide et globale.
o Maillage Raffiné : Offre une représentation beaucoup plus
précise des phénomènes de pression et de turbulence. Les
détails fins des variations de pression sont bien capturés, ce qui
permet une analyse plus approfondie.
• Temps de Calcul :
o Maillage Initial : 3.76 secondes, rapide mais moins précis.

o Maillage Raffiné : 300.82 secondes, beaucoup plus long mais

avec une précision nettement supérieure.
• Implications pour les Applications Pratiques :
o Le maillage raffiné est essentiel pour les applications nécessitant
une grande précision, comme la conception d'infrastructures
critiques ou les études environnementales détaillées.
o Le maillage initial peut suffire pour des analyses préliminaires ou
lorsque le temps de calcul est une contrainte majeure.

▪ Analyse du Champ de Vitesse [Vidéo]

ANALYSE

Début de la simulation :
o La vitesse de l'eau est initialement élevée près de la brèche du
barrage. On observe une accélération rapide de l'eau
s'écoulant vers l'aval.
o La vitesse est maximale au point de rupture, indiquant l'effet de
la gravité et de la différence de niveau d'eau.
Milieu de la simulation :
o L'eau s'écoule rapidement vers l'aval, avec des variations de
vitesse importantes. Les zones de turbulence sont visibles là où
l'eau rencontre des obstacles ou des changements de section.

Page 22 sur 43
 DAM BREAK

o Les gradients de vitesse montrent des zones d'accélération et de

décélération, caractéristiques de l'écoulement turbulent induit
par la rupture de barrage.
Fin de la simulation :
o L'eau s'est répandue en aval et la vitesse commence à diminuer.
La distribution de la vitesse devient plus uniforme, bien que des
zones de haute vitesse persistent dans les zones où l'eau continue
de se déplacer rapidement.
o Les fluctuations de vitesse sont moins prononcées, indiquant une
phase de stabilisation du mouvement de l'eau.
INTERPRETATION DES RESULTATS

• Distribution de la Vitesse :
o La distribution initiale montre une forte accélération de l'eau à
travers la brèche. Cette phase est critique pour comprendre la
force initiale de l'écoulement.
o La phase intermédiaire avec des vitesses élevées et des
turbulences intenses est typique des ruptures de barrage,
marquée par des mouvements chaotiques de l'eau.
o La phase finale montre une stabilisation progressive, mais encore
avec des zones de haute vitesse indiquant des résidus de
mouvement turbulent.
• Effets de la Turbulence :
o La turbulence joue un rôle majeur dans la distribution des vitesses.
Les variations soudaines et les gradients élevés indiquent des
mélanges intenses et des forces hydrodynamiques importantes.
o Les zones de turbulence peuvent être cruciales pour l'analyse de
l'érosion des sols, des forces sur les structures en aval, et des
impacts environnementaux.
• Comparaison avec le Champ de Pression :
o Les zones de haute pression dans la simulation de pression
correspondent aux zones de haute vitesse dans le champ de
vitesse. Cela confirme l'interdépendance des deux phénomènes
dans la dynamique de rupture de barrage.
o Les gradients de pression sont souvent associés à des gradients
de vitesse, renforçant l'importance de l'analyse conjointe pour
une compréhension complète.

Page 23 sur 43
 DAM BREAK

COMPARAISON ENTRE LE MAILLAGE INITIAL ET LE MAILLAGE RAFFINE [VIDEO]

Début de la simulation :
o Maillage Initial : La vitesse de l'eau est élevée près de la brèche
du barrage, mais les détails des gradients de vitesse sont moins
précis.
o Maillage Raffiné : La distribution de la vitesse est plus détaillée,
avec des gradients de vitesse mieux définis. Les variations subtiles
de la vitesse sont mieux capturées, permettant une analyse plus
fine des zones de haute accélération.
Milieu de la simulation :
o Maillage Initial : Les variations de vitesse sont visibles mais moins
définies. Les zones de turbulence sont présentes mais moins
détaillées.
o Maillage Raffiné : Les variations de vitesse sont beaucoup plus
nettes, avec des zones de turbulence clairement identifiées. Les
gradients de vitesse montrent une dynamique complexe et
précise de l'écoulement turbulent.
Fin de la simulation :
o Maillage Initial : La vitesse commence à diminuer et devient plus
uniforme, mais les détails sont moins visibles.
o Maillage Raffiné : La stabilisation de la vitesse est plus détaillée,
avec des zones de haute vitesse résiduelle clairement définies.
Les variations locales sont mieux capturées, offrant une vue plus
précise de la stabilisation du mouvement de l'eau.
Interprétation des Résultats
• Précision de la Simulation :
o Maillage Initial : Fournit une vue d'ensemble du champ de vitesse
mais manque de précision dans les détails des gradients de
vitesse et de la turbulence.
o Maillage Raffiné : Offre une représentation beaucoup plus
précise des phénomènes de vitesse et de turbulence. Les détails
fins des variations de vitesse et des zones de haute accélération
sont bien capturés.
• Effets de la Turbulence :
o Maillage Initial : La turbulence est visible mais moins détaillée.

o Maillage Raffiné : La turbulence est représentée de manière plus

réaliste et détaillée, avec des zones de haute turbulence
clairement définies. Cela permet une meilleure compréhension

Page 24 sur 43
 DAM BREAK

des forces hydrodynamiques et des impacts potentiels sur

l'environnement et les structures en aval.
• Temps de Calcul :
o Comme mentionné précédemment, le maillage raffiné prend
beaucoup plus de temps à calculer (300.82 secondes contre 3.76
secondes pour le maillage initial), mais offre une précision
nettement supérieure.

▪ Les lignes de courants [Vidéo]

ANALYSE ET INTERPRETATION DES RESULTATS :

Début de la Rupture :
o La première image montre le début de la rupture du barrage. On
peut observer la formation initiale de vagues et l'écoulement
rapide de l'eau.
o Les lignes d'influence ici montrent les vecteurs de vitesse et la
direction principale de l'écoulement.
Propagation de l'Eau :
o La deuxième image illustre la propagation de l'eau suite à la
rupture. L'écoulement devient turbulent avec des zones de
haute vitesse près du point de rupture.
o Les lignes d'influence indiquent une distribution complexe des
vecteurs de vitesse, montrant des tourbillons et des zones de
recirculation.
Expansion Latérale :
o À ce stade, l'eau commence à se répandre latéralement. Les
lignes d'influence montrent une dispersion des vecteurs de
vitesse, suggérant une réduction de la vitesse dans les zones
éloignées du point de rupture.
o Des zones de basse vitesse sont visibles aux bords de
l'écoulement.
Stabilisation Initiale :
o L'eau continue de s'étendre et commence à atteindre un état
plus stable. Les lignes d'influence montrent une diminution
progressive de la turbulence.

Page 25 sur 43
 DAM BREAK

o Les vecteurs de vitesse indiquent une direction plus uniforme, et

la vitesse globale de l'écoulement commence à diminuer.
Stabilisation Finale :
o À la fin, l'écoulement de l'eau se stabilise davantage. Les lignes
d'influence sont plus espacées, indiquant une vitesse plus faible
et une distribution plus uniforme de l'écoulement.
o Les zones de haute turbulence disparaissent progressivement,
laissant place à un écoulement plus laminaire.

▪ Vecteurs vitesses [Vidéo]

ANALYSE ET INTERPRETATION DES VECTEURS DE VITESSES :

Début de la Rupture :
o La première image montre les vecteurs de vitesse à l'instant initial
de la rupture du barrage.
o Les vecteurs sont orientés vers le bas, indiquant un écoulement
rapide et initial de l'eau depuis le point de rupture. La densité des
vecteurs de vitesse est élevée, ce qui montre une grande
variation de la vitesse dans cette région.
Propagation Initiale de l'Eau :
o La deuxième image illustre la propagation de l'eau quelques
instants après la rupture.
o Les vecteurs de vitesse montrent une direction principalement
horizontale avec une intensité élevée proche du point de
rupture. Il y a également des vecteurs qui montrent des
mouvements verticaux, indiquant une turbulence initiale.
Expansion Latérale :
o À ce stade, l'eau commence à s'étendre latéralement et les
vecteurs de vitesse montrent une dispersion.
o Les vecteurs indiquent une diminution de la vitesse à mesure que
l'eau s'éloigne du point de rupture. Des vecteurs de direction
opposée peuvent être observés, indiquant des tourbillons et des
mouvements de recirculation.
Stabilisation Initiale :

Page 26 sur 43
 DAM BREAK

o L'eau continue de se répandre et commence à atteindre un état

plus stable. Les vecteurs de vitesse sont plus uniformes et
montrent une diminution de l'intensité.
o Les vecteurs indiquent une direction plus régulière et moins de
turbulence, ce qui suggère une stabilisation progressive de
l'écoulement.
Stabilisation Finale :
o À la fin, les vecteurs de vitesse indiquent une stabilisation quasi
complète de l'écoulement de l'eau.
o Les vecteurs sont espacés uniformément avec une direction
constante, montrant une vitesse réduite et un écoulement
laminaire.

Simulation du volume avec l'obstacle de double taille

Le dossier de la simulation, ICI.
Pour le maillage raffiné, ICI.

a- Prétraitement :
▪ Création de la géométrie et du maillage :
Le fichier blockMeshDict est utilisé pour définir la géométrie et le maillage de
base pour la simulation de la rupture de digue. Voici une explication détaillée
de la géométrie créée dans ce fichier :
DIMENSIONS DE LA GEOMETRIE :

Le fichier blockMeshDict définit la géométrie et le maillage pour la

simulation dans OpenFOAM. Voici les principales sections du fichier :
- Les coordonnées des sommets (vertices) sont définies dans un système
de coordonnées cartésiennes :
vertices
(
(0 0 0)
(0.584 0 0)
(0.632 0 0)
(1.168 0 0)
(0 0.096 0)
(0.584 0.096 0)
(0.632 0.096 0)
(1.168 0.096 0)
(0 1.168 0)
(0.584 1.168 0)
(0.632 1.168 0)
(1.168 1.168 0)

Page 27 sur 43
 DAM BREAK

(0 0 0.1)
(0.584 0 0.1)
(0.632 0 0.1)
(1.168 0 0.1)
(0 0.096 0.1)
(0.584 0.096 0.1)
(0.632 0.096 0.1)
(1.168 0.096 0.1)
(0 1.168 0.1)
(0.584 1.168 0.1)
(0.632 1.168 0.1)
(1.168 1.168 0.1)
);

DEFINITION DES BLOCS ET MAILLAGE

• Blocs de maillage :
o Le maillage est défini par un ensemble de blocs hexagonaux,
chacun spécifié par les indices des sommets :

blocks
(
hex (0 1 5 4 12 13 17 16) (92 32 1) simpleGrading (1 1 1)
hex (2 3 7 6 14 15 19 18) (76 32 1) simpleGrading (1 1 1)
hex (4 5 9 8 16 17 21 20) (92 168 1) simpleGrading (1 1 1)
hex (5 6 10 9 17 18 22 21) (16 168 1) simpleGrading (1 1
1)
hex (6 7 11 10 18 19 23 22) (76 168 1) simpleGrading (1 1
1)
);

o Indices des Sommets : Les huit indices (0 1 5 4 12 13 17 16)

réfèrent aux sommets définis dans la section précédente.
o Résolution du Maillage : (92 32 1) indique le nombre de cellules
dans chaque direction (x, y, z).
o SimpleGrading : (1 1 1) spécifie la gradation du maillage. Une
valeur de 1 signifie une gradation uniforme.

DEFINITION DES FRONTIERES

• Conditions aux limites :

o Les frontières (boundaries) sont définies pour indiquer les types de
conditions aux limites appliquées aux différentes faces du
domaine de calcul :

boundary
(
leftWall

Page 28 sur 43
 DAM BREAK

{
type wall;
faces
(
(0 12 16 4)
(4 16 20 8)
);
}
rightWall
{
type wall;
faces
(
(7 19 15 3)
(11 23 19 7)
);
}
lowerWall
{
type wall;
faces
(
(0 1 13 12)
(1 5 17 13)
(5 6 18 17)
(2 14 18 6)
(2 3 15 14)
);
}
atmosphere
{
type patch;
faces
(
(8 20 21 9)
(9 21 22 10)
(10 22 23 11)
);
}
);

o leftWall, rightWall, lowerWall : Sont définies comme des murs

(wall) indiquant des conditions aux limites de type mur.
o atmosphere : Est défini comme une condition aux limites de type
patch, indiquant une interface ouverte avec l'atmosphère.

Page 29 sur 43
 DAM BREAK

DEFINITION DU VOLUME D'EAU INITIAL

Le fichier setFieldsDict est utilisé pour définir les conditions initiales du

champ de volume fractionnel d'eau dans le domaine de calcul. Voici une
explication détaillée du contenu de ce fichier :

• Contenu du fichier setFieldsDict :

3. Valeurs par défaut du champ :
o Le champ par défaut [Link] est initialisé à 0. Cela signifie
que, par défaut, il n'y a pas d'eau dans le domaine sauf dans les
régions spécifiquement définies.
defaultFieldValues
(
volScalarFieldValue [Link] 0
);

4. Définition des régions contenant de l'eau :

o Une région est définie où la fraction volumique d'eau
[Link] est initialisée à 10. Cela correspond au volume initial
d'eau dans le domaine.
o La région est définie par un volume rectangulaire spécifié par
deux coins opposés de la boîte.
regions
(
boxToCell
{
box (0 0 -1) (0.292 0.584 1);
fieldValues
(
volScalarFieldValue [Link] 10
);
}
);

o Explication de la région d'eau

- Volume de la boîte :
La boîte est définie par deux points (0 0 -1) et (0.292 0.584 1).
▪ (0 0 -1) : Coin inférieur gauche de la boîte.
▪ (0.292 0.584 1) : Coin supérieur droit de la boîte.
Ces points définissent une boîte rectangulaire en 3D qui contient
l'eau initiale.
- Fraction volumique de l'eau :
Dans cette région définie par la boîte, la fraction volumique de
l'eau [Link] est définie à 10. Cela signifie que cette région
est entièrement remplie d'eau au début de la simulation.

Page 30 sur 43
 DAM BREAK

▪ Raffinement du maillage :
Le raffinement du maillage est réalisé en modifiant le fichier
blockMeshDict, en doublant la résolution du maillage dans les zones
souhaitées, puis en relançant la simulation.
blocks
(
hex (0 1 5 4 12 13 17 16) (184 64 1) simpleGrading (1 1 1)
hex (2 3 7 6 14 15 19 18) (152 64 1) simpleGrading (1 1 1)
hex (4 5 9 8 16 17 21 20) (184 336 1) simpleGrading (1 1
1)
hex (5 6 10 9 17 18 22 21) (32 336 1) simpleGrading (1 1
1)
hex (6 7 11 10 18 19 23 22) (152 336 1) simpleGrading (1 1
1)
);

Le maillage est finement divisé en 200 cellules dans les directions x et y, et

en une seule cellule dans la direction z.

▪ Visualisation des résultats du prétraitement sur

Paraview :
Ici, vous trouver la vidéo de la simulation pour le premier maillage.
Ici, vous trouver la simulation du maillage raffiné.

b- Traitement
La partie traitement de la simulation avec obstacle est établie de la même
manière que celle de la simulation sans obstacle

c- Post-traitement

▪ Analyse Visuelle des Frames [Vidéo]

Voici les principales observations pour les cinq images extraites de la vidéo
de votre simulation de rupture de barrage avec obstacle :
Frame 1 (0s) :
• Le point de rupture initial est visible.
• L'eau commence à se déplacer rapidement en aval, créant une onde
de rupture marquée.
Frame 2 (5s) :
• L'onde de rupture s'est propagée davantage.

Page 31 sur 43
 DAM BREAK

• Les variations de vitesse et de pression sont observables dans la zone

de déplacement rapide de l'eau.
Frame 3 (10s) :
• La propagation continue de l'onde.
• Interactions possibles avec le terrain ou les obstacles visibles.
Frame 4 (15s) :
• Le flux d'eau se stabilise légèrement après l'onde de rupture initiale.
• Variations de direction et de vitesse du flux visibles.
Frame 5 (20s) :
• L'eau atteint un nouvel équilibre après la propagation de l'onde.
• Les effets résiduels de la rupture sont encore observables.

▪ Champs des pressions [Vidéo]

Début de la simulation (0s) :

• Au début, la pression est élevée dans la région où se trouvait
initialement l'eau, en amont du barrage. Ceci est dû au poids de l'eau
qui exerce une pression hydrostatique sur le fond et les parois du
réservoir.
Milieu de la simulation (5s, 10s) :
• Lorsque le barrage cède, l'eau commence à se déplacer rapidement
vers l'aval. On observe une diminution de la pression dans la région
initiale du réservoir et une augmentation de la pression en aval, où
l'eau s'accumule en rencontrant des obstacles ou le fond du canal.
• La dynamique des vagues de pression peut être observée à travers les
variations de couleurs, indiquant les gradients de pression causés par
l'écoulement turbulent de l'eau.
Fin de la simulation (15s, 20s) :
• À la fin, l'eau s'est répandue davantage en aval. La pression se stabilise
progressivement, mais des zones de haute pression peuvent persister
dans les régions où l'eau s'accumule ou rencontre des obstacles.

Page 32 sur 43
 DAM BREAK

• La pression est plus uniformément répartie comparée aux phases

initiales, mais avec des variations locales dues aux résidus de
mouvement et de turbulence de l'eau.
INTERPRETATION

Distribution de la pression :
• La distribution de la pression dans cette simulation montre clairement
l'effet de la rupture de barrage. La chute initiale de pression en amont
et l'augmentation en aval sont caractéristiques de ce type
d'événement.
Écoulement et turbulence :
• Les variations rapides et les gradients de pression indiquent un
écoulement très turbulent, typique des ruptures de barrage où l'eau se
déplace rapidement et de manière désordonnée.
Impacts environnementaux :
• Ces résultats peuvent être utilisés pour évaluer l'impact
environnemental d'une rupture de barrage, notamment sur les
structures en aval, l'érosion des sols et la sécurité des infrastructures
humaines.

▪ Champs de vitesse [Vidéo]

ANALYSE DU CHAMP DE VITESSE

Voici l'analyse des frames extraites de la vidéo des champs de vitesse à

différents moments de la simulation :
Début de la simulation (0s) :
• La vitesse de l'eau est initialement élevée près de la brèche du
barrage. On observe une accélération rapide de l'eau s'écoulant vers
l'aval.
• La vitesse est maximale au point de rupture, indiquant l'effet de la
gravité et de la différence de niveau d'eau.
Milieu de la simulation (5s, 10s) :

Page 33 sur 43
 DAM BREAK

• L'eau s'écoule rapidement vers l'aval, avec des variations de vitesse

importantes. Les zones de turbulence sont visibles là où l'eau rencontre
des obstacles ou des changements de section.
• Les gradients de vitesse montrent des zones d'accélération et de
décélération, caractéristiques de l'écoulement turbulent induit par la
rupture de barrage.
Fin de la simulation (15s, 20s) :
• L'eau s'est répandue en aval et la vitesse commence à diminuer. La
distribution de la vitesse devient plus uniforme, bien que des zones de
haute vitesse persistent dans les zones où l'eau continue de se déplacer
rapidement.
• Les fluctuations de vitesse sont moins prononcées, indiquant une phase
de stabilisation du mouvement de l'eau.
INTERPRETATION DES RESULTATS

Distribution de la Vitesse :
• La distribution initiale montre une forte accélération de l'eau à travers
la brèche. Cette phase est critique pour comprendre la force initiale de
l'écoulement.
• La phase intermédiaire avec des vitesses élevées et des turbulences
intenses est typique des ruptures de barrage, marquée par des
mouvements chaotiques de l'eau.
• La phase finale montre une stabilisation progressive, mais encore avec
des zones de haute vitesse indiquant des résidus de mouvement
turbulent.
Effets de la Turbulence :
• La turbulence joue un rôle majeur dans la distribution des vitesses. Les
variations soudaines et les gradients élevés indiquent des mélanges
intenses et des forces hydrodynamiques importantes.
• Les zones de turbulence peuvent être cruciales pour l'analyse de
l'érosion des sols, des forces sur les structures en aval, et des impacts
environnementaux.
Comparaison avec le Champ de Pression :
• Les zones de haute pression dans la simulation de pression
correspondent aux zones de haute vitesse dans le champ de vitesse.
Cela confirme l'interdépendance des deux phénomènes dans la
dynamique de rupture de barrage.

Page 34 sur 43
 DAM BREAK

• Les gradients de pression sont souvent associés à des gradients de

vitesse, renforçant l'importance de l'analyse conjointe pour une
compréhension complète.

▪ Lignes de courants [Vidéo]

ANALYSE DES LIGNES DE COURANTS

Voici l'analyse des frames extraites de la vidéo des lignes de courant à

différents moments de la simulation :
Début de la Rupture (0s) :
• La première image montre le début de la rupture du barrage. On peut
observer la formation initiale de vagues et l'écoulement rapide de
l'eau.
• Les lignes d'influence ici montrent les vecteurs de vitesse et la direction
principale de l'écoulement.
Propagation de l'Eau (5s) :
• La deuxième image illustre la propagation de l'eau suite à la rupture.
L'écoulement devient turbulent avec des zones de haute vitesse près
du point de rupture.
• Les lignes d'influence indiquent une distribution complexe des vecteurs
de vitesse, montrant des tourbillons et des zones de recirculation.
Expansion Latérale (10s) :
• À ce stade, l'eau commence à se répandre latéralement. Les lignes
d'influence montrent une dispersion des vecteurs de vitesse, suggérant
une réduction de la vitesse dans les zones éloignées du point de
rupture.
• Des zones de basse vitesse sont visibles aux bords de l'écoulement.
Stabilisation Initiale (15s) :
• L'eau continue de s'étendre et commence à atteindre un état plus
stable. Les lignes d'influence montrent une diminution progressive de la
turbulence.
• Les vecteurs de vitesse indiquent une direction plus uniforme, et la
vitesse globale de l'écoulement commence à diminuer.

Page 35 sur 43
 DAM BREAK

Stabilisation Finale (20s) :

• À la fin, l'écoulement de l'eau se stabilise davantage. Les lignes
d'influence sont plus espacées, indiquant une vitesse plus faible et une
distribution plus uniforme de l'écoulement.
• Les zones de haute turbulence disparaissent progressivement, laissant
place à un écoulement plus laminaire.
INTERPRETATION DES RESULTATS

Distribution des Vecteurs de Vitesse :

• La distribution initiale montre une forte accélération de l'eau à travers
la brèche. Cette phase est critique pour comprendre la force initiale de
l'écoulement.
• La phase intermédiaire avec des vitesses élevées et des turbulences
intenses est typique des ruptures de barrage, marquée par des
mouvements chaotiques de l'eau.
• La phase finale montre une stabilisation progressive, mais encore avec
des zones de haute vitesse indiquant des résidus de mouvement
turbulent.
Effets de la Turbulence :
• La turbulence joue un rôle majeur dans la distribution des vitesses. Les
variations soudaines et les gradients élevés indiquent des mélanges
intenses et des forces hydrodynamiques importantes.
• Les zones de turbulence peuvent être cruciales pour l'analyse de
l'érosion des sols, des forces sur les structures en aval, et des impacts
environnementaux.
Comparaison avec les Champs de Pression et de Vitesse :
• Les zones de haute pression et de haute vitesse dans les simulations
précédentes correspondent aux zones de turbulence intense
observées dans les lignes de courant. Cela confirme l'interdépendance
des phénomènes dans la dynamique de rupture de barrage.
• Les gradients de pression et de vitesse renforcent l'importance de
l'analyse conjointe pour une compréhension complète de l'événement.

▪ Vecteurs vitesse [Vidéo]

Page 36 sur 43
 DAM BREAK

ANALYSE ET INTERPRETATION DES VECTEURS DE VITESSES

Voici l'analyse des frames extraites de la vidéo des vecteurs de vitesse à

différents moments de la simulation :
Début de la Rupture (0s) :
• La première image montre les vecteurs de vitesse à l'instant initial de la
rupture du barrage.
• Les vecteurs sont orientés vers le bas, indiquant un écoulement rapide
et initial de l'eau depuis le point de rupture. La densité des vecteurs de
vitesse est élevée, ce qui montre une grande variation de la vitesse
dans cette région.
Propagation Initiale de l'Eau (5s) :
• La deuxième image illustre la propagation de l'eau quelques instants
après la rupture.
• Les vecteurs de vitesse montrent une direction principalement
horizontale avec une intensité élevée proche du point de rupture. Il y a
également des vecteurs qui montrent des mouvements verticaux,
indiquant une turbulence initiale.
Expansion Latérale (10s) :
• À ce stade, l'eau commence à s'étendre latéralement et les vecteurs
de vitesse montrent une dispersion.
• Les vecteurs indiquent une diminution de la vitesse à mesure que l'eau
s'éloigne du point de rupture. Des vecteurs de direction opposée
peuvent être observés, indiquant des tourbillons et des mouvements de
recirculation.
Stabilisation Initiale (15s) :
• L'eau continue de se répandre et commence à atteindre un état plus
stable. Les vecteurs de vitesse sont plus uniformes et montrent une
diminution de l'intensité.
• Les vecteurs indiquent une direction plus régulière et moins de
turbulence, ce qui suggère une stabilisation progressive de
l'écoulement.
Stabilisation Finale (20s) :
• À la fin, les vecteurs de vitesse indiquent une stabilisation quasi
complète de l'écoulement de l'eau.

Page 37 sur 43
 DAM BREAK

• Les vecteurs sont espacés uniformément avec une direction constante,

montrant une vitesse réduite et un écoulement laminaire.

Points Clés des Différences entre les Simulations

▪ Distribution de la Pression :
Sans Obstacle : La distribution de la pression devient uniforme plus
rapidement.
Avec Obstacle : Les variations de pression sont plus complexes et persistent
plus longtemps, surtout autour de l'obstacle.

▪ Vitesse de l'Écoulement :
Sans Obstacle : La vitesse diminue progressivement et l'écoulement
devient plus uniforme.
Avec Obstacle : La vitesse est plus variable, avec des zones de haute
vitesse et de recirculation autour de l'obstacle, indiquant une turbulence plus
élevée et plus persistante.

▪ Turbulence :
Sans Obstacle : La turbulence est principalement présente au début et au
milieu de la simulation.
Avec Obstacle : La turbulence est plus intense et durable, en particulier
autour de l'obstacle, avec des tourbillons et des mouvements de recirculation
marqués.

Page 38 sur 43
 DAM BREAK

CONCLUSION
La modélisation numérique de la rupture de barrage à l'aide de la
dynamique des fluides computationnelle (CFD) nous a permis de mieux
comprendre les dynamiques complexes des écoulements en cas de rupture
de barrage, avec et sans obstacle. Les simulations réalisées avec OpenFOAM
ont fourni des informations détaillées sur les champs de pression, de vitesse et
les lignes de courant, offrant une visualisation claire des impacts de tels
événements sur l'environnement et les structures en aval.
Principales Conclusions :
1. Champ de Pression :
o Sans Obstacle : La distribution de la pression devient uniforme
plus rapidement après la rupture. Les variations de pression sont
moins complexes et persistent moins longtemps.
o Avec Obstacle : Les variations de pression sont plus prononcées
et durables, notamment autour de l'obstacle. La présence de
l'obstacle crée des zones de haute pression en amont et des
gradients de pression plus complexes.
2. Champ de Vitesse :
o Sans Obstacle : La vitesse de l'eau diminue progressivement et
l'écoulement devient plus uniforme au fil du temps. La turbulence
est principalement présente au début et au milieu de la
simulation.
o Avec Obstacle : La vitesse est plus variable, avec des zones de
haute vitesse et des mouvements de recirculation autour de
l'obstacle. La turbulence est plus intense et durable,
particulièrement autour de l'obstacle, avec des tourbillons
marqués.
3. Lignes de Courant :
o Sans Obstacle : Les lignes de courant montrent une dispersion
uniforme des vecteurs de vitesse avec une réduction progressive
de la turbulence.
o Avec Obstacle : Les lignes de courant indiquent une dispersion
complexe des vecteurs de vitesse avec des zones de
recirculation et de turbulence plus marquées autour de
l'obstacle.
4. Vecteurs de Vitesse :

Page 39 sur 43
 DAM BREAK

o Sans Obstacle : Les vecteurs de vitesse montrent une forte

accélération initiale suivie d'une diminution progressive de
l'intensité et une direction plus uniforme.
o Avec Obstacle : Les vecteurs de vitesse révèlent des interactions
complexes avec l'obstacle, montrant des zones de haute
turbulence et des mouvements de recirculation persistants.
Implications Pratiques :
Les résultats de ces simulations mettent en évidence l'importance de la
présence d'obstacles dans la gestion des risques liés à la rupture de barrage.
Les obstacles modifient significativement les dynamiques de pression et de
vitesse, augmentant la complexité de l'écoulement et la persistance de la
turbulence. Ces informations sont cruciales pour la conception de structures
de protection, l'évaluation des impacts environnementaux et la planification
des interventions d'urgence.
Recommandations :
Pour une gestion efficace des risques associés à la rupture de barrage, il
est recommandé de :
• Utiliser des modèles numériques détaillés pour simuler différents
scénarios de rupture, incluant la présence d'obstacles.
• Concevoir des infrastructures en tenant compte des dynamiques
complexes révélées par les simulations CFD.
• Planifier des stratégies d'atténuation des risques basées sur des analyses
précises des champs de pression, de vitesse et de turbulence.
En conclusion, la modélisation CFD de la rupture de barrage est un outil
indispensable pour comprendre et anticiper les impacts de tels événements,
permettant de mettre en place des mesures de protection efficaces et de
minimiser les risques pour les populations et les infrastructures en aval.

Page 40 sur 43
 DAM BREAK

ANNEXE
La méthode VoF (Volume of Fluid) Vs la méthode des volumes
finis
▪ Méthode des Volumes Finis (Finite Volume Method -
FVM)
La méthode des volumes finis est une technique de discrétisation utilisée
pour résoudre les équations de Navier-Stokes. Elle divise le domaine de calcul
en petits volumes finis et applique les lois de conservation (masse, quantité de
mouvement, énergie) à chacun de ces volumes.

▪ Méthode Volume of Fluid (VoF)

La méthode VoF est une technique utilisée pour suivre la surface libre entre
deux fluides immiscibles, comme l'eau et l'air. Elle permet de capturer la
position de l'interface fluide avec une fonction de volume de fluide dans
chaque cellule de la grille.

▪ Utilisation Conjointe
En CFD, la méthode des volumes finis peut être utilisée pour la
discrétisation spatiale, tandis que la méthode VoF est utilisée pour capturer
l'interface entre les fluides dans ce cadre de discrétisation.

Page 41 sur 43
 DAM BREAK

BIBLIOGRAPHIE
[Link]
[Link]
simulation-of-dam-break-flow-over

User Guide : [Link]

Copernicus HESS
Numerical simulation of dam-break flow over a simplified building model
using OpenFOAM

Page 42 sur 43