Télécharger becalc.eso

Retour à la liste

Numérotation des lignes :

  1. C BECALC SOURCE DURA 06/07/11 21:17:09 5481
  2. SUBROUTINE BECALC(PE,PVE,TE,PS,XL,DX,RUG,QINI,XW,NNP,XN,TN,EN,BN,
  3. $ KIMP,NT,NX,XX,XP,XPV,XT,XY,XU,XHF,XQ,XQW,XQA,XRE,XDH,RECU,
  4. $ XKUL,XKUT1,XKUT2,XKUT3,XKUT4)
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8 (A-H,O-Z)
  7. C
  8. C operateur FUITE
  9. C>>> parametres d'entree:
  10. C PE : pression totale a l'entree
  11. C PVE : pression de vapeur a l'entree
  12. C TE : temperature d'entree
  13. C PS : pression de sortie
  14. C XL : longueur de la fissure
  15. C DX : taille relative maille fluide
  16. C RUG : rugosite
  17. C QINI: debit initialisation
  18. C XW : parametre de condensation
  19. C XN TN EN : abscisses temperature et epaisseur de la paroi
  20. C (dimension NNP)
  21. C KIMP: parametre d'impression
  22. C NT : dimension max. grille fluide
  23. C RECU,XKUL,XKUT1,XKUT2,XKUT3,XKUT4 : coef lois de frot utilisateur
  24. C
  25. C>>> parametres de sortie:
  26. C NX : nombre de points grille fluide
  27. C pour la grille fluide
  28. C XX : coordonnees relatives
  29. C XP XPV XT XU XY XHF XQ XQW XQA XRE
  30. C pression, pression vapeur, temperature, titre, vitesse
  31. C coef. echange, debit total, debit eau (liq.+vap.), debit air
  32. C Reynolds
  33.  
  34. DIMENSION XX(NT),XP(NT),XT(NT),XY(NT),XU(NT),XPV(NT)
  35. DIMENSION XHF(NT),XQ(NT),XQW(NT),XQA(NT),XRE(NT),XDH(NT)
  36. DIMENSION XN(NNP),TN(NNP),EN(NNP),BN(NNP)
  37. C
  38. IF(KIMP.GE.1) THEN
  39. WRITE(6,*) 'entree becalc '
  40. ENDIF
  41. TE = TE + 273.15
  42. IF (QINI.LT.1.E-5) THEN
  43.  
  44. C>>> calcul du debit correspondant a un melange gazeux isotherme
  45. C a la temperature d'entree
  46. C
  47. E=(EN(1)+EN(NNP))/2.
  48. B=(BN(1)+BN(NNP))/2.
  49. CALL BCALQ0(PE,PVE,TE,PS,XL,E,B,RUG,Q0,RECU,XKUL,XKUT1,XKUT2,
  50. * XKUT3,XKUT4)
  51. ELSE
  52. Q0 = QINI
  53. ENDIF
  54. C
  55. C>>> calcul du debit solution
  56.  
  57. CALL BCALDP(PE,PVE,TE,PS,XL,DX,RUG,Q0,XW,NNP,XN,TN,EN,BN,
  58. $ KIMP,NT,NX,XX,XP,XT,XY,XU,XHF,XQ,XQW,
  59. $ Q,QAE,XRE,XDH,RECU,XKUL,XKUT1,XKUT2,XKUT3,XKUT4)
  60.  
  61. C>>> calcul de la pression vapeur et des debits resultats
  62.  
  63. DO 15 I=1,NX
  64.  
  65. IF (XY(I).GT.0.99) THEN
  66. PHI=PVE/PE
  67. XPV(I)=PHI*XP(I)
  68. ELSE
  69. T=XT(I)+273.15
  70. XPV(I)=BPSAT(T)
  71. ENDIF
  72.  
  73. C>> a revoir
  74. XQA(I) = QAE
  75. XQW(I) = XQ(I)
  76. XQ (I) = XQA(I)+XQW(I)
  77.  
  78. 15 CONTINUE
  79. C
  80. C>> a revoir
  81. XHF(1) = XHF(2)
  82. XRE(1) = XRE(2)
  83. XDH(1) = XDH(2)
  84.  
  85. RETURN
  86. END
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales