Esope source files

   1 : C GRACOL    SOURCE    PV        09/11/05    21:15:20     6535           
   2 :       SUBROUTINE GRACOL(RG,RTG,CVGS,AEQ1,DELSGS,
   3 :      &                  XSOL,PGS,TGS, 
   4 :      &                  RC1,UG,UGS,UPXSOL,UMXSOL, 
   5 :      &                  DUDT,DPDT)
   6 : C************************************************************************
   7 : C
   8 : C PROJET            :  CASTEM 2000
   9 : C
  10 : C NOM               :  GRACOL
  11 : C
  12 : C DESCRIPTION       :  Voir aussi FUNCOL
  13 : C                      Calcul de DPDT, DUDT sur le choc gauche, 
  14 : C                      si TG != TGS
  15 : C
  16 : C LANGAGE           :  FORTRAN 77 
  17 : C      
  18 : C AUTEUR            :  A. BECCANTINI DRN/DMT/SEMT/LTMF  
  19 : C
  20 : C************************************************************************
  21 : C
  22 : C APPELES           : NONE
  23 : C
  24 : C************************************************************************
  25 : C
  26 : C**** ENTREES:
  27 : C
  28 : C     (G, GS sont les etats Gauche et droite (Gauche Star))
  29 : C
  30 : C     RG            = R du gaz a gauche (et a droite)
  31 : C
  32 : C     RTG           = RG * TG
  33 : C
  34 : C     CVGS          = cv(TGS), cv a droite
  35 : C
  36 : C     AEQ1          = 2.0D0 * (ETHEGS - ETHEG) + RTGS - RTG
  37 : C
  38 : C     DELSGS        = SQRT((AEQ1 * AEQ1) + (4.0D0 * RTGS * RTG))
  39 : C
  40 : C     XSOL          = ROG / ROGS
  41 : C
  42 : C     PGS, TGS      = pression et temperature a droite
  43 : C
  44 : C     RC1           = 2.0D0*(HTHEGS - HTHEG)
  45 : C
  46 : C     UG, UGS       = vitesses
  47 : C
  48 : C     UPXSOL,UMXSOL = 1 + XSOL, 1 - XSOL
  49 : C
  50 : C**** SORTIES:      
  51 : C
  52 : C     DUDT          = derivé de la vitesse (a droite) par raport a 
  53 : C                     T (droite)
  54 : C
  55 : C     DPDT          = derivé de la pression (a droite) par raport a 
  56 : C                     T (droite)
  57 : C
  58 : C
  59 : C**** ATTENTION: Le calcul de DPDT, DUDT sur le choc gauche est possible
  60 : C                si TGS =! TG; sinon il y a des functions singulieres
  61 : C     
  62 : C************************************************************************
  63 : C
  64 : C HISTORIQUE (Anomalies et modifications éventuelles)
  65 : C
  66 : C HISTORIQUE :  Créé le 08.02.00
  67 : C
  68 : C************************************************************************
  69 : C
  70 : C N.B.: Toutes les variables sont DECLAREES
  71 : C
  72 :       IMPLICIT INTEGER(I-N)
  73 :       REAL*8 RG,CVGS,AEQ1,RTG,DELSGS
  74 :      &      ,DA1DT,DDELDT
  75 :      &      ,DXDT,XSOL,PGS,TGS
  76 :      &      ,RC1,UG,UGS,UMXSOL,UPXSOL
  77 :      &      ,DUDT,DC1DT, DPDT
  78 : C     
  79 : C**** DXDT = d(XSOL)/dT
  80 : C     XSOL t.c. RTG XSOL*XSOL + AEQ1 * XSOL - RTGS = 0
  81 : C     DA1DT = d(AEQ1)/dT
  82 : C
  83 :       DA1DT = (2.0D0 * CVGS) + RG
  84 : C
  85 : C**** DDELT = d(DELSGS)/dT
  86 : C
  87 :       DDELDT = ((AEQ1 * DA1DT) + (2.0D0 * RTG * RG))
  88 :      &         / DELSGS 
  89 :       DXDT = 0.5D0 / RTG * (DDELDT - DA1DT)
  90 : C
  91 : C**** DC1DT = d(RC1)/dT = 2 cp(TGS)
  92 : C
  93 :       DC1DT = DA1DT + RG
  94 :       DUDT = ((DC1DT * UMXSOL)-(2.0D0 * RC1 / UPXSOL * DXDT))/UPXSOL
  95 :       DUDT = -0.5D0 * DUDT / (UG - UGS)
  96 : C      
  97 :       DPDT =  (PGS / TGS) - (PGS * DXDT / XSOL)
  98 :       RETURN
  99 :       END
 100 :          
 101 :           
 102 :  
 103 :  
 104 :