Télécharger la plaquette Cast3M


fr en

Télécharger fvcent.eso

Retour à la liste

Numérotation des lignes : 

fvcent

  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