Esope source files

   1 : C BCAR      SOURCE    CHAT      05/01/12    21:36:59     5004           
   2 :       SUBROUTINE BCAR(RA,RV,PHI,CV,CA,Q,RQ,R,QA,QV,CPS)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8 (A-H,O-Z)
   5 : C
   6 : C     operateur FUITE
   7 : C>>>  calcul des caracteristiques de l'ecoulement
   8 : C
   9 : C     parametres d'entree
  10 : C       RA,RV    : constante gaz parfait air, vapeur
  11 : C       PHI      : rapport pression partielle vapeur/pression totale
  12 : C       CV,CA    : chaleur massique vapeur, air
  13 : C       Q        : debit total
  14 : C     parametres de sortie
  15 : C       RQ: constante du melange gazeux
  16 : C       R : autre constante du melange gazeux (loi Gaz Parfait)
  17 : C       QA: debit air
  18 : C       QV: debit vapeur
  19 : C       CPS :  chaleur massique du melange gazeux en zone surchauffee
  20 : C
  21 :       IF (PHI.LT.0.001) THEN
  22 :          RQ = 1E8
  23 :       ELSE
  24 :          RQ = (1-PHI)*RV/PHI/RA
  25 :       ENDIF
  26 : 
  27 :       R  = PHI/RV+(1-PHI)/RA
  28 :       R  = 1/R
  29 :       QA = Q*RQ/(1+RQ)
  30 :       QV = Q/(1+RQ)
  31 :       CPS = (CV*PHI/RV + CA*(1-PHI)/RA) / (PHI/RV + (1-PHI)/RA)
  32 :       RETURN
  33 :       END
  34 :  
  35 :