Esope source files

   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 :