Analyse d’un écoulement graduellement varié

Analyse d'un Écoulement Graduellement Varié (Méthode du Pas à Pas)

Analyse d'un écoulement graduellement varié par la méthode du pas à pas

Comprendre les Écoulements Graduellement Variés

Lorsqu'un écoulement à surface libre rencontre un obstacle (comme un barrage, un seuil) ou un changement de pente, la hauteur d'eau ne peut plus rester constante (écoulement uniforme). Elle varie progressivement le long du canal : c'est un écoulement graduellement varié. Le calcul de la forme de cette ligne d'eau, appelée courbe de remous ou de rabattement, est crucial pour prédire les niveaux d'eau et les risques d'inondation. La méthode du pas à pas direct est une technique numérique classique pour calculer ce profil en partant d'un point où la hauteur est connue.

Données de l'étude

On étudie la courbe de remous créée par un petit barrage sur un canal trapézoïdal en terre.

Caractéristiques du canal et de l'écoulement :

  • Débit (\(Q\)) : \(15 \, \text{m}^3/\text{s}\)
  • Largeur au fond (\(b\)) : \(5 \, \text{m}\)
  • Pente des talus (\(m\)) : \(1.5\) (1.5H:1V)
  • Pente du fond (\(S_f\)) : \(0.1 \% = 0.001\)
  • Coefficient de Manning (\(n\)) : \(0.022 \, \text{s/m}^{1/3}\)
  • Condition de contrôle : La hauteur d'eau juste en amont du barrage (\(x=0\)) est de \(y_0 = 3.0 \, \text{m}\).
Schéma : Courbe de Remous en Amont d'un Barrage
Lit du canal (pente S_f) Ligne d'eau uniforme (y_n) Ligne d'eau critique (y_c) Courbe de remous (M1) Barrage y_0 = 3.0m x = 0 x (amont)

Profil en long montrant le relèvement de la ligne d'eau à l'approche d'un obstacle.


Questions à traiter

  1. Calculer la hauteur d'eau normale (\(y_n\)) pour ce canal.
  2. Calculer la hauteur d'eau critique (\(y_c\)).
  3. Déterminer le type de pente du canal et le type de la courbe de remous.
  4. En utilisant la méthode du pas à pas direct, calculer la distance \(\Delta x\) entre la section où la hauteur est de 3.0 m et celle où elle est de 2.8 m.
  5. Calculer le profil complet de la ligne d'eau depuis \(y=3.0\) m jusqu'à une hauteur de 1% supérieure à la hauteur normale (\(1.01 \times y_n\)), par pas de 0.2 m. Présenter les résultats sous forme de tableau et de graphique.

Correction : Analyse d'un Écoulement Graduellement Varié

Question 1 : Hauteur d'Eau Normale (\(y_n\))

Principe :

On cherche la hauteur \(y_n\) qui satisfait l'équation de Manning pour un écoulement uniforme. Comme pour l'exercice précédent, l'équation est implicite et doit être résolue par itérations.

\[\frac{A^{5/3}}{P^{2/3}} = \frac{Q \cdot n}{S_f^{1/2}}\]
Calcul du terme de droite :
\[ \begin{aligned} \frac{Q \cdot n}{S_f^{1/2}} &= \frac{15 \times 0.022}{(0.001)^{1/2}} \\ &= \frac{0.33}{0.03162} \\ &\approx 10.436 \end{aligned} \]

L'équation à résoudre est donc \(A^{5/3} / P^{2/3} = 10.436\).

Détail de la résolution par itérations

On cherche la valeur de \(y_n\) qui vérifie l'égalité. Pour un canal trapézoïdal : \(A=(5+1.5y_n)y_n\) et \(P=5+2y_n\sqrt{1+1.5^2} \approx 5+3.61y_n\).

Essai \(y_n\) (m) A (m²) P (m) \(f(y_n) = A^{5/3}/P^{2/3}\) Observation
1 2.0 16.00 12.22 19.2 Trop grand
2 1.5 10.88 10.41 10.08 Trop petit
3 1.55 11.36 10.59 10.77 Trop grand
4 1.53 11.18 10.52 10.50 Très proche de 10.436
Résultat Question 1 : La hauteur d'eau normale est \(y_n \approx 1.53 \, \text{m}\).

Question 2 : Hauteur d'Eau Critique (\(y_c\))

Principe :

La hauteur critique \(y_c\) est atteinte lorsque le nombre de Froude vaut 1. Pour une section non rectangulaire, elle satisfait la condition suivante :

\[\frac{Q^2}{g} = \frac{A^3}{B}\]

Où \(B\) est la largeur au miroir, \(B = b + 2my_c\). L'équation se résout aussi par itérations.

Calcul du terme de gauche :
\[\frac{Q^2}{g} = \frac{15^2}{9.81} = \frac{225}{9.81} \approx 22.936\]
Résolution

On cherche \(y_c\) tel que \(A^3/B = 22.936\). Après itérations (non détaillées ici), on trouve :

\[y_c \approx 1.21 \, \text{m}\]
Résultat Question 2 : La hauteur critique est \(y_c \approx 1.21 \, \text{m}\).

Question 3 : Type de Pente et de Courbe

Principe :

Le type de pente dépend de la comparaison entre la hauteur normale \(y_n\) et la hauteur critique \(y_c\). Le type de courbe dépend de la position de la hauteur d'eau réelle par rapport à \(y_n\) et \(y_c\).

Analyse :
  • Type de pente : On a \(y_n \approx 1.53 \, \text{m}\) et \(y_c \approx 1.21 \, \text{m}\). Puisque \(y_n > y_c\), la pente est de type **faible** (M, pour "Mild slope").
  • Type de courbe : La hauteur d'eau au barrage (\(y_0=3.0\)) est supérieure à \(y_n\) et \(y_c\). La courbe de remous est donc dans la "Zone 1". Il s'agit d'une courbe de type **M1**.
Résultat Question 3 : La pente est faible et la courbe de remous est de type M1.

Question 4 : Calcul d'un Pas (\(y=3.0\)m à \(y=2.8\)m)

Principe :

On utilise la formule du pas à pas direct. Elle calcule la distance horizontale \(\Delta x\) nécessaire pour que la hauteur d'eau passe d'une valeur \(y_1\) à \(y_2\).

\[\Delta x = \frac{E_1 - E_2}{S_e - S_f}\]

Ici, on regarde vers l'amont, donc on calcule \(x_2-x_1\). La formule s'écrit :

\[\Delta x = \frac{E_2 - E_1}{S_f - \bar{S_e}}\]

Où \(E\) est l'énergie spécifique (\(y+V^2/2g\)) et \(\bar{S_e}\) est la moyenne de la pente de la ligne d'énergie aux points 1 et 2. On a \(S_e = (n V / R_h^{2/3})^2\).

Calcul pour les deux sections

Section 1 (aval) : \(y_1 = 3.0 \, \text{m}\)

\[ \begin{aligned} A_1 &= (5 + 1.5 \times 3) \times 3 = 28.5 \, \text{m}^2 \\ P_1 &= 5 + 3.61 \times 3 = 15.83 \, \text{m} \\ R_{h1} &= A_1 / P_1 = 1.80 \, \text{m} \\ V_1 &= Q / A_1 = 0.526 \, \text{m/s} \\ E_1 &= y_1 + V_1^2 / (2g) = 3.014 \, \text{m} \\ S_{e1} &= (n V_1 / R_{h1}^{2/3})^2 = 0.000059 \end{aligned} \]

Section 2 (amont) : \(y_2 = 2.8 \, \text{m}\)

\[ \begin{aligned} A_2 &= (5 + 1.5 \times 2.8) \times 2.8 = 25.76 \, \text{m}^2 \\ P_2 &= 5 + 3.61 \times 2.8 = 15.11 \, \text{m} \\ R_{h2} &= A_2 / P_2 = 1.70 \, \text{m} \\ V_2 &= Q / A_2 = 0.582 \, \text{m/s} \\ E_2 &= y_2 + V_2^2 / (2g) = 2.817 \, \text{m} \\ S_{e2} &= (n V_2 / R_{h2}^{2/3})^2 = 0.000074 \end{aligned} \]
Calcul de \(\Delta x\)
\[ \begin{aligned} \bar{S_e} &= \frac{0.000059 + 0.000074}{2} = 0.0000665 \\ \Delta x &= \frac{3.014 - 2.817}{0.001 - 0.0000665} \\ &= \frac{0.197}{0.0009335} \\ &\approx 211 \, \text{m} \end{aligned} \]
Résultat Question 4 : Il faut environ 211 mètres pour que la hauteur d'eau passe de 3.0 m à 2.8 m.

Question 5 : Profil Complet de la Ligne d'Eau

Principe :

On répète le calcul de la question 4 par pas successifs pour tracer la courbe de remous complète jusqu'à atteindre une hauteur proche de la hauteur normale \(y_n\).

La cible est d'atteindre \(1.01 \times y_n = 1.01 \times 1.53 \approx 1.55 \, \text{m}\).

Tableau de Calcul du Pas à Pas Direct
y (m) E (m) S_e \(\bar{S_e}\) \(\Delta x\) (m) x (m)
3.003.0145.90E-05--0
2.802.8177.40E-056.65E-05211211
2.602.6249.41E-058.41E-05211422
2.402.4371.22E-041.08E-04210632
2.202.2581.64E-041.43E-04209841
2.002.0892.33E-041.98E-042111052
1.801.9363.62E-042.97E-042181270
1.601.8086.52E-045.07E-042591529
1.551.7827.83E-047.18E-04921621
Profil de la Ligne d'Eau (Courbe de Remous M1)
Résultat Question 5 : La courbe de remous s'étend sur plus de 1.6 km en amont du barrage avant de rejoindre quasiment la hauteur d'eau normale.

Quiz Rapide : Testez vos connaissances (Récapitulatif)

1. Une courbe de remous de type M1 se produit lorsque :

2. Dans la méthode du pas à pas, la pente de la ligne d'énergie \(S_e\) représente :

3. Si la rugosité du canal \(n\) augmentait (par exemple, à cause de la végétation), la courbe de remous pour le même débit :


Glossaire

Écoulement Graduellement Varié
Écoulement permanent dans lequel la profondeur et la vitesse varient lentement le long du canal. Les lignes de courant sont quasi parallèles et la répartition des pressions est hydrostatique.
Courbe de Remous (Backwater Curve)
Profil longitudinal de la surface libre de l'eau lorsque celle-ci est située au-dessus de la hauteur normale, typiquement en amont d'un barrage ou d'un rétrécissement.
Pente Faible (Mild Slope)
Une pente de canal est dite "faible" si la hauteur d'eau normale (\(y_n\)) est supérieure à la hauteur critique (\(y_c\)), ce qui signifie qu'un écoulement uniforme sur cette pente serait fluvial (subcritique).
Pente de la Ligne d'Énergie (\(S_e\))
Pente qui représente la perte d'énergie par unité de longueur due au frottement. En écoulement uniforme, \(S_e = S_f\). En écoulement varié, elle est calculée par la formule de Manning réarrangée.
Méthode du Pas à Pas Direct (Direct Step Method)
Procédure numérique pour calculer le profil d'une ligne d'eau en écoulement graduellement varié. Elle calcule la distance \(\Delta x\) correspondant à une variation de hauteur \(\Delta y\) prédéfinie.
Hydraulique à Surface Libre - Exercice d'Application

D’autres exercices d’hydraulique à surface libre:

Estimation du Transport en Suspension
Estimation du Transport en Suspension

Hydraulique : Estimation du Transport en Suspension à partir d'un Profil de Concentration Estimation du Transport en Suspension à partir d'un Profil de Concentration Contexte : La Rivière Trouble Contrairement au charriage qui concerne les sédiments lourds se...

Calcul du Transport de Charriage
Calcul du Transport de Charriage

Hydraulique : Calcul du Transport de Charriage avec la Formule de Meyer-Peter & Müller Calcul du Transport de Charriage avec la Formule de Meyer-Peter & Müller Contexte : Le Mouvement des Sédiments sur le Lit d'une Rivière Le transport solide dans les rivières est un...

Conception d’un Canal Stable
Conception d’un Canal Stable

Hydraulique : Conception d'un Canal Stable (Critères de Lane) Conception d'un canal stable selon les critères de Lane Contexte : L'Équilibre Fragile des Canaux en Terre Un canal creusé dans un matériau non consolidé (sable, gravier, limon) est dit "stable" si...

Calcul du Débit dans un Canal Composite
Calcul du Débit dans un Canal Composite

Hydraulique : Calcul du Débit dans un Canal Composite Calcul du débit dans un canal composite (lit mineur et majeur) Contexte : Les Rivières en Crue Lorsqu'une rivière sort de son lit principal (lit mineurCanal principal d'une rivière, qui contient l'écoulement pour...

Analyse de la transition d’un canal
Analyse de la transition d’un canal

Hydraulique : Transition de Pente Forte à Pente Faible (Profil S1 à M3) Analyse de la transition d'un canal de pente forte à une pente faible (profil S1 à M3) Contexte : Le Ressaut Hydraulique, une Transition Énergétique Brutale En hydraulique à surface libre,...

Conception d’un canal d’amenée
Conception d’un canal d’amenée

Hydraulique : Conception d'un Canal d'Amenée pour une Micro-Centrale Hydroélectrique Conception d'un canal d'amenée pour une micro-centrale hydroélectrique Contexte : L'Art de Guider l'Eau en Douceur Un canal d'amenéeCanal à ciel ouvert qui transporte l'eau depuis une...

Charge Hydraulique sur un Radier
Charge Hydraulique sur un Radier

Hydraulique : Calcul de la Charge Hydraulique sur le Radier d'un Ouvrage Calcul de la charge hydraulique sur le radier d'un ouvrage Contexte : La Stabilité des Ouvrages Face à l'Eau Tout ouvrage hydraulique (seuil, barrage, pile de pont) placé dans un écoulement...

Problème de la « douce-amère »
Problème de la « douce-amère »

Hydraulique : Problème de la "douce-amère" - Rencontre de Deux Cours d'Eau Problème de la "douce-amère" : rencontre de deux cours d'eau Contexte : La Lutte Invisible des Densités Lorsqu'un fleuve d'eau douce rencontre la mer, ou lorsque deux rivières de salinités (et...

Analyse d’un Seuil de Mesure Noyé
Analyse d’un Seuil de Mesure Noyé

Analyse d'un Seuil de Mesure Noyé Analyse d'un Seuil de Mesure Noyé Comprendre les Seuils Noyés Un seuil est un ouvrage construit en travers d'un canal ou d'une rivière pour mesurer le débit ou contrôler le niveau d'eau. On parle d'écoulement dénoyé (ou libre) lorsque...

Dimensionnement d’un Ponceau (Culvert)
Dimensionnement d’un Ponceau (Culvert)

Dimensionnement d'un Ponceau (Culvert) pour un Débit de Projet Dimensionnement d'un Ponceau (Culvert) Comprendre le Fonctionnement d'un Ponceau Un ponceau (ou dalot) est un ouvrage hydraulique qui permet à un cours d'eau de passer sous un obstacle, comme une route ou...

0 commentaires
Soumettre un commentaire

Votre adresse e-mail ne sera pas publiée. Les champs obligatoires sont indiqués avec *