Esope source files

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