Calcul de l'accélération à partir du champs de vitesse 
  si v = u + v + w
  accélération = dv/dt = drond v/drond t + (gradient V)V soit :
  selon x : u drond u/drond x + v drond u /drond y + w drond u / drond z 
  selon y : u drond v /drond x + v drond v/drond y + w drond v / drond z
  de même pour z 
Vecteur tourbillon : 
  Gradient vectoriel V = 
  drond/drond x drond /drond y drond / drond z vectoriel u v w
Déterminer les lignes de courant à t=0 selon les conditions de Cauchy :
  Intégrer les equations selon la variable avec laquelle ils sont dérivées, sans oublier la constante puis la calculer à t=0
1 bar = 1 10 ^5 Pascal
1 atm = 1.013125 bar
unité de base le Pascal
Rayon hydraulique Rh = Ah /Ph 
Fr nombre de Froude = Q/sqrt(g*Sm/B) Sm section mouillée qui vaut hc * qq chose
hauteur critique quand Fr = 1 sinon prendre avec la hauteur indiquée 
Si h > hc, ecoulement fluvial sinon torrentiel 
Equation de Manning Strickler :
  Q = 1/n * SmRh^2/3sqrt(I) I pente longitudinale n rugosite de Manning
Contrainte de cisaillement :
  To = mu du/dy
Conservation de la masse pour un fluide incompressible :
  drond ro/drond t + div(rov)=0
Pour un fluide incompressible div(v) = 0
conservation du moment pour un fluide incompressible :
  ro drond v /drond t + ro grad v * v = fv - grad p + mu laplacien de v
  Si Re faible on enleve le gradient de v
  Si Re fort on enlève le laplacien
  
Ecoulement de Couette :
  On suppose que l'écoulement est laminaire :
    v(x,z)=vx(x,z)ex
    Equation de Continuité :
      drond Vx /dronx x = 0 donc V dépend de z
    CL : P(x=0,z)=P(x=L,z) pas de gradient de pression
         v(z=0)= 0
         v(z=h)=Vo ex
    Les équations de N-S donnent :
      drdndP/drond x = mu drond2V/drond2z soit
      v(z)= Voz/h
      drondP/drond z = -rog soit
      P(Z) = Po - rogz