Télécharger la plaquette Cast3M


fr en

Télécharger becalc2.eso

Retour à la liste

Numérotation des lignes : 

becalc2

  1. C BECALC2   SOURCE    DURA      06/07/11    21:17:01     5481
  2.        SUBROUTINE BECALC2(PE,PVE,TE,PS,XL,DX,RUG,QINI,XW,NNP,XN,TN,EN,
  3.      $       BN,KIMP,NT,NX,XX,XP,XPV,XT,XY,XU,XHF,XQ,XQW,XQA,XRE,XDH,
  4.      $       RECU,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 becalc2  '
  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     WRITE(6,*) 'Q0' ,Q0
  55. C>>>  calcul du debit solution
  56.       CALL BCALDP2(PE,PVE,TE,PS,XL,DX,RUG,Q0,XW,NNP,XN,TN,EN,BN,
  57.      $           KIMP,NT,NX,XX,XP,XT,XY,XU,XHF,XQ,XQW,
  58.      $               Q,QAE,XRE,XDH,RECU,XKUL,XKUT1,XKUT2,XKUT3,XKUT4)
  59.  
  60. C>>>  calcul de la pression vapeur et des debits resultats
  61.  
  62.       DO 15 I=1,NX
  63.  
  64.          IF (XY(I).GT.0.99) THEN
  65.             PHI=PVE/PE
  66.             XPV(I)=PHI*XP(I)
  67.          ELSE
  68.             T=XT(I)+273.15
  69.             XPV(I)=PSATT0(T)
  70.          ENDIF
  71.  
  72. C>>   a revoir
  73.          XQA(I) = QAE
  74.          XQW(I) = XQ(I)
  75.          XQ (I) = XQA(I)+XQW(I)
  76.  
  77.   15  CONTINUE
  78. C
  79. C>>   a revoir
  80.        XHF(1) = XHF(2)
  81.        XRE(1) = XRE(2)
  82.        XDH(1) = XDH(2)
  83.  
  84.       RETURN
  85.       END
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  

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