Télécharger la plaquette Cast3M


fr en

Télécharger gracol.eso

Retour à la liste

Numérotation des lignes : 

gracol

  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.  
  105.  

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