Télécharger fvcent.eso

Retour à la liste

Numérotation des lignes :

  1. C FVCENT SOURCE CHAT 05/01/13 00:13:50 5004
  2. SUBROUTINE FVCENT(NESP,
  3. & GAMG,ROG,PG,UNG,UTG,
  4. & GAMD,ROD,PD,UND,UTD,
  5. & YG,YD,FLUX,
  6. & CELLT)
  7. C
  8. C************************************************************************
  9. C
  10. C PROJET : CASTEM 2000
  11. C
  12. C NOM : FVCENT
  13. C
  14. C DESCRIPTION : Formulation Volumes Finis pour les Equations
  15. C d'Euler Multi-Especes relatives à un melange
  16. C de gaz ideals.
  17. C
  18. C Calcul du flux aux interfaces avec la mèthode
  19. C centrée
  20. C
  21. C LANGAGE : FORTRAN 77
  22. C
  23. C AUTEUR : A. BECCANTINI DRN/DM2S/SFME/LTMF
  24. C
  25. C************************************************************************
  26. C
  27. C APPELES
  28. C
  29. C************************************************************************
  30. C
  31. C**** Entrées:
  32. C
  33. C NESP = nombre d'especes considérées dans les Equations
  34. C d'Euler
  35. C
  36. C GAMG, GAMD = les "gamma" du gaz (gauche et droite)
  37. C
  38. C ROG, ROD = les densités
  39. C
  40. C PG, PD = les pressions
  41. C
  42. C UNG, UND = vitesses normales
  43. C
  44. C UTG, UTD = vitesses tangentielles
  45. C
  46. C YG, YD = tables des fractiones massiques
  47. C
  48. C**** Sorties:
  49. C
  50. C FLUX = table du flux a l'interface, i.e.
  51. C (rho*un, rho*un*un + p, rho*un*ut, rho*un*ht,
  52. C rho*un*y1, ...)
  53. C
  54. C CELLT = condition de stabilité, i.e.
  55. C
  56. C dT/diamax < cellt
  57. C
  58. C************************************************************************
  59. C
  60. C HISTORIQUE (Anomalies et modifications éventuelles)
  61. C
  62. C HISTORIQUE :
  63. C
  64. C************************************************************************
  65. C
  66. C N.B.: Toutes les variables sont DECLAREES
  67. C
  68. IMPLICIT INTEGER(I-N)
  69. INTEGER NESP,IESP
  70. REAL*8 ROG,ROD,PG,PD,UNG,UND,UTG,UTD
  71. & ,GAMG,GAMD
  72. & ,YG(*),YD(*),FLUX(*),CELLT,ROUN,ECIN
  73. C
  74. C**** A gauche
  75. C
  76. ROUN=ROG*UNG
  77. FLUX(1)=ROUN
  78. FLUX(2)=(ROUN*UNG)+PG
  79. FLUX(3)=ROUN*UTG
  80. ECIN=0.5D0*((UNG*UNG)+(UTG*UTG))
  81. FLUX(4)=(UNG*((GAMG/(GAMG-1.0D0))*PG))+(ROUN*ECIN)
  82. DO IESP=1,NESP,1
  83. FLUX(4+IESP)=ROUN*YG(IESP)
  84. ENDDO
  85. CELLT=(((GAMG*PG/ROG)**0.5D0)+ABS(UNG))
  86. C
  87. C**** A droite
  88. C
  89. ROUN=ROD*UND
  90. FLUX(1)=0.5D0*(FLUX(1)+ROUN)
  91. FLUX(2)=0.5D0*(FLUX(2)+((ROUN*UND)+PD))
  92. FLUX(3)=0.5D0*((FLUX(3)+(ROUN*UTD)))
  93. ECIN=0.5D0*((UND*UND)+(UTD*UTD))
  94. FLUX(4)=0.5D0*(FLUX(4)+((UND*((GAMD/(GAMD-1.0D0))*PD))+(ROUN*ECIN)
  95. $ ))
  96. DO IESP=1,NESP,1
  97. FLUX(4+IESP)=0.5D0*(FLUX(4+IESP)+(ROUN*YD(IESP)))
  98. ENDDO
  99. CELLT=MAX(CELLT,(((GAMD*PD/ROD)**0.5D0)+ABS(UND)))
  100. CELLT=1.0D0/CELLT
  101. C
  102. RETURN
  103. END
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  

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