Esope source files

   1 : C FVLHTF    SOURCE    BECC      10/06/11    21:15:08     6690           
   2 :       SUBROUTINE FVLHTF(NESP,NSCA,
   3 :      &                  GAMG,ROG,PG,UNG,UTG,UT1G,ETHERG,EFORMG,
   4 :      &                  GAMD,ROD,PD,UND,UTD,UT1D,ETHERD,EFORMD,
   5 :      &                  YG,YD,SCAG,SCAD,FLU1,FLU2,
   6 :      &                  CELLT)
   7 : C************************************************************************
   8 : C
   9 : C NOM               :  FVLHTT
  10 : C
  11 : C DESCRIPTION       :  Formulation Volumes Finis  pour les Equations 
  12 : C                      d'Euler Multi-Especes relatives à un melange 
  13 : C                      de gaz ideals.
  14 : C
  15 : C                      Calcul du flux aux interfaces avec la methode
  16 : C                      "Flux Vector Splitting" de van Leer - Hanel.  
  17 : C
  18 : C                      (voir:
  19 : C                       1)  BECCANTINI, PAILLERE,
  20 : C                           "Upwind Flux Splitting Schemes..."
  21 : C                           RAPPORT DMT 97//268
  22 : C                       2)  B. VAN LEER
  23 : C                           "Flux Vector Splittings for the 1990's"
  24 : C                           1990. Invited Lecture, NASA Lewis Research 
  25 : C                           Center
  26 : C                       3)  B. VAN LEER
  27 : C                           "Flux Vector Splitting for the Euler
  28 : C                            Equations"
  29 : C                            Report ICASE No. 82-30, 1982  )
  30 : C
  31 : C LANGUAGE          :  FORTRAN 77 
  32 : C      
  33 : C AUTEUR            :  A. BECCANTINI DRN/DMT/SEMT/TTMF  
  34 : C
  35 : C************************************************************************
  36 : C
  37 : C APPELES 
  38 : C
  39 : C FVLHTF ------ FSPVLH
  40 : C    |
  41 : C     --------- FSMVLH     
  42 : C          
  43 : C
  44 : C************************************************************************
  45 : C
  46 : C**** Entrées:
  47 : C
  48 : C     NESP            =  nombre d'especes considérées dans les Equations
  49 : C                        d'Euler
  50 : C
  51 : C     NSCA            =  nombre de scalaires passifs a transporter 
  52 : C
  53 : C     GAMG, GAMD      =  les "gamma" du gaz (gauche et droite)
  54 : C
  55 : C     ROG, ROD        =  les densités 
  56 : C
  57 : C     PG, PD          =  les pressions
  58 : C
  59 : C     UNG, UND        =  vitesses normales
  60 : C
  61 : C     UTG, UTD        =  vitesses tangentielles
  62 : C
  63 : C     UT1G, UT1D      =  vitesses tangentielles
  64 : C
  65 : C     ETHERG, ETHERD  =  les energies sensibles
  66 : C
  67 : C     EFORMG, EFORMD  =  les energies de formation
  68 : C
  69 : C     YG, YD          =  tables des fractiones massiques
  70 : C
  71 : C     SCAG, SCAD      =  tables des scalaires passifs
  72 : C
  73 : C**** Sorties:
  74 : C
  75 : C     FLU1            =  table du flux a l'interface dans le repaire
  76 : C                        (n,t), i.e.
  77 : C                        (rho*un, rho*un*un + p, rho*un*ut, rho*un*ht,
  78 : C                         rho*un*y1, ...)
  79 : C
  80 : C     FLU2            =  table de travail, utilisé ici mais definie
  81 : C                        avant 
  82 : C
  83 : C     CELLT           =  condition de stabilité, i.e.
  84 : C
  85 : C                        dT < diamax/cellt
  86 : C                        cellt a la dimension d'une vitesse
  87 : C      
  88 : C************************************************************************
  89 : C
  90 : C HISTORIQUE (Anomalies et modifications éventuelles)
  91 : C
  92 : C HISTORIQUE :  Créé le 6.1.98 
  93 : C
  94 : C               21.02.00 transport de scalaires passifs
  95 : C
  96 : C               2010: les energies des formations sont prises en compte
  97 : C
  98 : C************************************************************************
  99 : C
 100 : C N.B.: Toutes les variables sont DECLAREES
 101 : C                          
 102 : C
 103 : C      IMPLICIT NONE
 104 :       IMPLICIT INTEGER(I-N)
 105 :       INTEGER NESP, I1, NSCA
 106 :       REAL*8  GAMG,ROG,PG
 107 :      &       ,GAMD,ROD,PD
 108 :      &       ,AG2,AG,UNG,UTG,UT1G,MG,ETHERG,HTG,EFORMG
 109 :      &       ,AD2,AD,UND,UTD,UT1D,MD,ETHERD,HTD,EFORMD
 110 :      &       ,CELL,CELLT,CELLT0,AMG,AMD,PSRO,ECIN
 111 :      &       ,YG(*),YD(*),FLU1(*),FLU2(*),SCAG(*),SCAD(*)
 112 : C
 113 : C
 114 : C**** YG, YD, FLU1, FLU2 déjà definis avant
 115 : C
 116 : C     Dans le cas Euler monoespece, on doit
 117 : C     avoir  :
 118 : C     YG(1) = YD(1) = 0.0D0
 119 : C
 120 : C
 121 : C**** Onde de "gauche" a "droite"
 122 : C      
 123 : C     AG  = vitesse du son
 124 : C     MG  = nombre de Mach normal
 125 : C     HTG = enthalpie totale
 126 : C     
 127 :       PSRO = PG / ROG
 128 :       ECIN = 0.5D0 * ((UNG * UNG) + (UTG * UTG) + (UT1G * UT1G))
 129 :       HTG = ETHERG + PSRO + ECIN + EFORMG
 130 :       AG2 = GAMG * PSRO
 131 :       AG = SQRT(AG2)
 132 :       MG = UNG / AG
 133 :       AMG = ABS(MG)
 134 :       CELL = (2.0D0*GAMG + (AMG*(3.0D0-GAMG)))/(GAMG+3.0D0)
 135 :       CELLT=(UNG+AG)/CELL
 136 : C     
 136 :       CALL FSPVLH(NESP,NSCA,GAMG,ROG,MG,AG,UTG,HTG,YG,SCAG,FLU1)
 137 : C      
 138 : C**** Onde de "droite" a "gauche".
 139 : C     
 140 : C     AD  = vitesse du son
 141 : C     MD  = nombre de Mach normal
 142 : C     HTD = enthalpie totale
 143 : C      
 144 :       PSRO = PD / ROD
 145 :       ECIN = 0.5D0 * ((UND * UND) + (UTD * UTD) + (UT1D * UT1D))
 146 :       HTD = ETHERD + PSRO + ECIN + EFORMD
 147 :       AD2 = GAMD * PSRO
 148 :       AD = SQRT(AD2)
 149 :       MD = UND / AD
 150 :       AMD = ABS(MD)
 151 :       CELL = (2.0D0*GAMD + (AMD*(3.0D0-GAMD)))/(GAMD+3.0D0)
 152 :       CELLT0=(UND+AD)/CELL
 153 :       CELLT = MIN(CELLT,CELLT0)
 154 : C     
 154 :       CALL FSMVLH(NESP,NSCA,GAMD,ROD,MD,AD,UTD,HTD,YD,SCAD,FLU2)
 155 : C       
 156 : C
 157 :       DO I1 = 1, (NSCA+NESP+4)
 158 :          FLU1(I1) = FLU1(I1) + FLU2(I1) 
 159 :       ENDDO
 160 : C
 161 : C
 162 :       RETURN
 163 :       END
 164 :