La relation d'Euler
Rappel : Leonhard Euler 1707-1783
Fondamental : Analyse d'une turbine élémentaire
C'est en généralisant l'énergie développée par un moulinet à tout type de condition d'entrée et de sortie qu'on obtient la relation d'Euler.
Pour la démonstration, considérons une turbine élémentaire constituée d'un tube à section variable dont l'entrée (indice 1) et la sortie (indice 2) sont tangentielles et parcourues par un débit de l'entrée vers la sortie. Le tube tourne autour d'un axe à une vitesse angulaire.
Considérons ce qui se passe à l'entrée en terme de quantité de mouvement.
L'effort tangentiel qui en résulte est positif et dans le sens de rotation : \(F_1 = \rho Q C_{u1}\).
Le couple moteur est : \(M_1= r_1 \times F_1\) .
La puissance générée à l'entrée : \(P_1= \omega M_1\)
L'énergie à l'entrée s'exprime donc par : \(E_1 =U_1 C_{u1}\)
Pour la sortie :
\(F_2=-\rho Q C_{u2}\)
\(M_2 = r_2 \times F_2\)
\(P_2= \omega M_2\)
\(E_2 = -U_2 C_{u2}\)
La puissance motrice résultant des vitesses à l'entrée et à la sortie sera :
Il en résulte la relation d'Euler permettant d'obtenir l'énergie générée par la turbine élémentaire :
Remarque :
La relation d'Euler nous indique que la puissance produite par une turbine vient de la déviation du filet fluide la traversant. Le même raisonnement s'applique à la pompe lorsque la puissance est négative, c'est-à-dire que la puissance doit être fournie pour générer le débit et la chute.
Dans cette présentation, pour simplifier, on a utilisé des entrée et sortie tangentielles. Dans une turbine réelle, l'écoulement n'est pas parfaitement tangentiel et donc pour appliquer la relation d'Euler il faut connaître les angles de l'écoulement pour en extraire les vitesses tangentielles \(C_u\) . Donc, connaissant les angles fluides à l'entrée et à la sortie d'une roue, on connaît l'énergie captée. L'efficacité énergétique et les pertes sont donc contenues dans les angles du fluide. C'est donc une affaire de géométrie.
On rappelle que la puissance est le produit d'un débit massique par une énergie massique.
On observe qu'il n'y aucune simplification ou hypothèse dans cette expression qui est donc exacte et que les notions de viscosité et turbulence ne sont pas utiles ici.
La relation d'Euler est l'équation maîtresse qui fait le lien entre l'énergie disponible et la géométrie de l'écoulement puis finalement la géométrie de la turbine.
Le concepteur hydraulicien doit donc générer une géométrie de roue avec les angles pour produire l'énergie souhaitée.
Définition : L'énergie interne
Complément : Généralisation du couple sur une cascade dans un repère cylindrique
Pour une cascade en rotation, puisque \(P = M \omega\), de la relation d'Euler :
\(P =\rho Q (r_1 \omega C_{u1} - r_2 \omega C_{u2})\), on peut écrire :
\(M=\rho Q (r_1 C_{u1} - r_2 C_{u2})\)
Également, pour une cascade fixe, à partir de l'équation de Newton adaptée au repère cylindrique : \(F=\rho Q(C_{u1} - C_{u2})\), on peut écrire :
\(M=\rho Q (r_1 C_{u1} - r_2 C_{u2})\)
On observe donc qu'on obtient la même équation pour le couple que la cascade tourne ou pas.
Complément : L'angle d'écoulement dans un tourbillon libre
Dans un repère cylindrique, tournant autour de l'axe z, il y a un tourbillon libre, dont la circulation s'y exprime : \(\Gamma = r \ C_u\)
Aussi, le débit, qui peut s'y exprimer : \(Q= C_m\ A\) , pour \(C_m = C_r\) la vitesse débitante de l'écoulement dans le plan méridien à la section considérée au rayon \(r\).
L'angle de l'écoulement \(\gamma\) peut s'y exprimer ainsi : \(\tan \gamma=\frac{C_u}{C_m}\)
et en prenant \(A\) comme étant la section cylindrique au rayon \(r\) et d'une hauteur \(b\), soit \(A=2\pi r b\),
si on développe :
\(\tan \gamma=\frac{C_u}{C_m}=\cfrac{\frac{\Gamma}{r}}{\frac{Q}{2\pi r b}}=\frac{\Gamma}{Q}2\pi b\)
Le dernier terme de cette équation \(b\) étant une constante liée à la position, on peut conclure que l'angle de l'écoulement \(\gamma\) est connu.
Un tourbillon libre est donc caractérisé par un angle d'écoulement à la section considérée.
Fondamental : L'aubage fait dévier l'angle de l'écoulement
De part et d'autre de l'aubage, le long d'une ligne d'écoulement, il y a deux tourbillons libres chacun caractérisé par un angle d'écoulement.
C'est l'aubage qui dévie l'écoulement et en change les angles.
Méthode : Application simplifiée au tracé d'une roue
Tracer l'aubage d'une turbine est le processus qui permet de relier les conditions hydrauliques et la performance souhaitée à la définition géométrique de l'aubage. Comme la relation d'Euler relie les conditions hydrauliques aux angles du fluide, il s'agit donc de créer une géométrie qui s'adapte à ces angles.
Les références Moulin[5] , Miller IV[6] , Rossgatterer[7] et plusieurs autres auteurs décrivent le développement que nous reprenons ici.
Ce développement s'appuie sur les concepts suivants :
La roue est constituée de plusieurs aubes, chacune possédant une surface intrados et une surface extrados. Un surface médiale virtuelle, ici appelée
« âme »
, peut en être déduite.
Le fluide dans la roue de la turbine se déplace sur des surfaces de révolution qui intersectent le plan méridien pour former des courbes.
Le trajet de la particule fluide se déplaçant sur la surface de révolution, entre le bord d'attaque et le bord de fuite, trace une abscisse curviligne en trois dimensions qu'on peut nommer filet ou ligne de courant.
L'évolution des angles le long du trajet :
contribue, par la relation d'Euler à la transformation de l'énergie hydraulique en énergie mécanique.
et peut être représentée par une transformation conforme appelée transformation de Beltrami[9] qui est une sorte de mise à plat en deux dimensions de l'aubage (voir Plan de la cascade dans la figure suivante).
La projection cylindrique de ce trajet sur un plan méridien trace une abscisse curviligne en deux dimensions.
Brève explication de la transformation conforme de Beltrami
Avant d'arriver au plan conforme u-v, on passe par le plan des longueurs m-n où :
Avec \(m\) la longueur curviligne dans le plan méridien qui s'exprime en fonction de \(s\) :
\(m=m(s) = \int_0^s \sqrt {\left({dr \over ds}\right)^2+\left({dz \over ds}\right)^2} ds\)
Et \(n\) la longueur curviligne dans le sens transverse qui s'exprime aussi en fonction de \(s\) :
\(n=r\theta=r(s)\theta(s)\)
Ce plan des longueurs, n'est pas pratique pour suivre l'évolution des angles du filet et donc de l'écoulement. On a donc besoin de faire appel à un changement de coordonnées.
En divisant les coordonnées \((m,r\theta)\) par \(r\) on obtient les nouvelles variables suivantes :
\(u = \int_0^s{{dm\over ds}\over r} ds = \int_0^m{1\over r} dm\) où \(u(0)=0\) et
\(v(s)=\theta(s)\)
Ainsi, l'évolution de l'angle d'écoulement \(\alpha\) s'exprime :
\(\tan\alpha={rd\theta \over dm}={d\theta \over dm/r}={dv \over du}\)
Pour la méthode complète et le logiciel de génération d'un profil hydraulique cliquez ici.
Exemple : Recherche des angles moyens du fluide à l'entrée et à la sortie d'une roue Francis
Approche unidimensionnel pour définir un tracé.
En première étape pour une approche simplifiée, on pose les conditions suivantes :
L'aubage est représenté par une seule ligne de courant.
L'étude s'intéresse d'abord au fonctionnement au point de rendement maximum.
La rotation fluide \(C_{u2}\) à la sortie de la roue est donc prévue nulle.
Sa géométrie 3D est complexe et se veut un assemblage de 2 plans orthogonaux correspondant respectivement au débit massique et à l'énergie massique
Plan méridien
Plan conforme
Voici une analyse d'une turbine que l'on assimilera à une turbine élémentaire dont l'écoulement est modélisé par un filet fluide monodimensionnel.
Voici plan méridien représentant les principales caractéristiques de la roue.
On y trouve le diamètre de la roue à 2 m , la hauteur du distributeur à 0,7 m.
À cela s'ajoutent les conditions hydrauliques pour le point de meilleur rendement :
Les directrices alimentent la roue avec un angle fluide de 70 degrés, le débit est de 30 m³/s et la chute est de 63m et la vitesse de rotation est de 300 tpm.
La base de la conception du tracé se fait en 1D sur l'abscisse curviligne de la ligne de courant moyenne entre le point 1 situé en entrée et le point 2 à la sortie de la roue. Il faut noter qu'à l'amont du point 1, l'énergie \(UC_u\) est constante et égale à \(U_1 C_{u1}\). De même à l'aval du point 2, l'énergie est assumée égale à \(U_2 C_{u2}\) le long de la ligne de courant.
Dans le plan de cascade u-v, on représente une transformation de l'aubage intersecté par la surface de révolution générée par le filet moyen. La coordonnée \(v\) représente la position angulaire \(\theta\) et \(u\) une représentation étirée en \(r^{-1}\) de la longueur du profil dans le plan méridien. La périodicité cylindrique que l'on retrouve dans le domaine 3D est représentée par une périodicité linéaire dont le pas dépend du nombre d'aube. Ainsi, pour 13 aubes, le pas est de 360/13 degrés.
À l'entrée de la roue :
Du rayon \(R_1\) et de la vitesse angulaire on obtient la vitesse du solide \(U_1\) et la vitesse radiale moyenne à l'entrée \(C_{r1}\) basée sur le débit traversant la section de l'anneau de hauteur \(b\) et de rayon \(R_1\).
\(C_{u1}\) et \(C_{r1}\) résulte du triangle de vitesse imposé par l'angle fluide \(\gamma_1\) résultant de l'effet des directrices.
\(W_{u1}\) est la vitesse tangentielle du fluide vue par l'aubage, c'est la vitesse tangentielle relative qui s'obtient en soustrayant \(U_1\) de \(C_{u1}\) .
L’arc-tangente de \(W_{u1}/C_{r1}\) donne l'angle \(\alpha_1\) qui est l'angle du fluide tel que vue par la roue lorsqu'elle tourne à 300 tpm.
Le bord d'attaque de l'aubage sera conçu à cet angle pour limiter les pertes par chocs à l'entrée de la roue.
À la sortie de la roue :
Puisqu'à la sortie de la roue, l'écoulement est axial, on choisit le rayon \(R_2\) comme représentant la position centrale du filet fluide .
En première approximation, on peut estimer que la vitesse \(C_{u2}\) à la sortie de la roue est nulle puisque c'est l'objectif du concepteur pour réduire les pertes dans l'aspirateur.
La vitesse axiale à la sortie \(C_{z2}\) doit permettre de satisfaire l'équation de continuité et se calcul en multipliant le débit par la section de sortie de la roue au diamètre \(D\) de 2m.
Puisque \(C_{u2}\) est nulle, La vitesse tangentielle relative du fluide doit annuler la vitesse du solide.
L'arc-tangente de \(W_{u2}/C_{z2}\) donne l'angle \(\alpha_2\) qui est l'angle du fluide qui doit être produit par l'aubage.
Conseil : La forte sensibilité aux angles de l'écoulement
Dans le fichier joint, on peut observer la très grande sensibilité de la performance à l'angle fluide à l'entrée \(\gamma_1\). Quelques dixièmes de degrés causent des variations en 1% du rendement.
La prédiction de l'angle fluide \(\gamma_1\) est délicate et on s'aperçoit qu'il y a des valeurs qu'on peut imposer dans la feuille de calcul qui donnent des résultats qui ne sont pas physiques.
La CFD peut être utilisée avec précaution pour prédire ces angles et raffiner la conception de l'aubage avant la validation finale sur modèle réduit.
Avec le tracé résultant on obtient une variation des angles d'écoulement entre l'entrée et la sortie de 51 degrés. La longueur de l'abscisse curviligne entre 1 et 2 n'est pas déterminée. En fait, il faudra considérer d'une part dans quelle longueur limite la courbure du fluide peut s'obtenir sans séparation de l'aubage et d'autre part si le champ de pression absolu est suffisamment étalé pour éviter la cavitation.
