Sources Esope | Cast3M

bcar
C BCAR      SOURCE    CHAT      05/01/12    21:36:59     5004
      SUBROUTINE BCAR(RA,RV,PHI,CV,CA,Q,RQ,R,QA,QV,CPS)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
C
C     operateur FUITE
C>>>  calcul des caracteristiques de l'ecoulement
C
C     parametres d'entree
C       RA,RV    : constante gaz parfait air, vapeur
C       PHI      : rapport pression partielle vapeur/pression totale
C       CV,CA    : chaleur massique vapeur, air
C       Q        : debit total
C     parametres de sortie
C       RQ: constante du melange gazeux
C       R : autre constante du melange gazeux (loi Gaz Parfait)
C       QA: debit air
C       QV: debit vapeur
C       CPS :  chaleur massique du melange gazeux en zone surchauffee
C
      IF (PHI.LT.0.001) THEN
         RQ = 1E8
      ELSE
         RQ = (1-PHI)*RV/PHI/RA
      ENDIF
 
      R  = PHI/RV+(1-PHI)/RA
      R  = 1/R
      QA = Q*RQ/(1+RQ)
      QV = Q/(1+RQ)
      CPS = (CV*PHI/RV + CA*(1-PHI)/RA) / (PHI/RV + (1-PHI)/RA)
      RETURN
      END