Télécharger gracol.eso

Retour à la liste

Numérotation des lignes :

  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