Description du code

 

 

Le code réalisé durant ce projet reprend les équations formulées dans la section précédente et utilise le langage Matlab pour les coder.

Le code commence d'abord par :

  • Début du programme

 

  • Lire tous les paramètres du problème.

 

  • Initialise les variables.

 

  • Calcule la pression dans le réservoir.

 

  • Débute une boucle temporelle

 

  •   Débute une boucle sur z profondeur du puits
  1. Calcule toutes les propriétés du mélange
  2. Calcule tous les coefficients de pertes de charge et le coefficients d'échange thermique convectif du liquide.
  3. Calcule la pression dans le puits avec la pression hydrostatique et les pertes de charge du mélange.
  4. La température sera égale à la température de saturation qui est calculé à partir de la pression.
    1. Si la température est inférieure à la température de saturation alors :
      1. Si le précédent titre de vapeur est supérieur à 0 ce qui veut dire que le fluide était en ébullition et se compose de la vapeur et du liquide, alors:
        • Calcule le coefficients d'échange thermique convectif du mélange
        • Le flux thermique transporté ( qui prend en compte l'advection, la diffusion et les pertes thermiques ) sert à calculer le nouveau titre massique
      2. Si le titre de vapeur est égale à 0 alors le fluide est entièrement composé de liquide, et le flux thermique transporté sert à chauffer encore le liquide.
    2. Si la température est supérieure ou égale à la température de saturation :
      1. Si le titre de vapeur est entre 0 et 1 alors la température devient égale à la température de saturation et le flux de chaleur transporté contribue au changement de phase.
      2. Si le titre de vapeur est égale à 1 alors le fluide est entièrement en phase vapeur et le flux de chaleur sert à chauffer la vapeur.
  • Fin de la boucle sur z profondeur du puits

 

  • Calcul des conditions aux limites

 

  • Vérification de la condition de convergence

 

  • Fin de la boucle temporelle

 

  • Écriture des fichiers résultats

 

  • Fin du programme