Pertes de charge - Churchill

[invitea09f2f16]

Bonjour à tous,

Étant en stage, je suis actuellement sur un problème impliquant des pertes de charges dans des conduites cylindriques. Après avoir fait quelques tests avec la formule de Colebrook, je cherche à déterminer ces pertes de charges en utilisant la formule de Churchill . Celle-ci étant normalement plus stable pour n'importe quel Reynolds et n'importe quel rapport k/D. Cependant, les résultats que j'obtiens pour un Reynolds de 957811, k = 0.01 mm et D = 0.9 m ne sont pas du tout en accord avec le diagramme de Moody. J'aurais donc aimé avoir un regard neuf sur mon script Matlab si possible ou alors une petite explication, si quelqu'un a, car je dois omettre quelque chose.

Re = 957811;
k = 0.01;
D = 900;

A = (2.457*log((7/Re)^0.9+0.27*(k/D))^(-1))^16
B = (37530/Re)^16

Lambda = 8*((8/Re)^12+(A+B)^(-1.5))^(1/12)*10^-4

Merci d'avance,

Florian
-----

[invitea09f2f16]

Dans le code j'ai fait une petite erreur, le vrai code est le suivant (sans le *10^-4 dans le calcul de Lambda) :

Re = 957811;
k = 0.01;
D = 900;

A = (2.457*log((7/Re)^0.9+0.27*(k/D))^(-1))^16
B = (37530/Re)^16

Lambda = 8*((8/Re)^12+(A+B)^(-1.5))^(1/12)