Esope source files

   1 : C FLUVL2    SOURCE    CHAT      05/01/13    00:04:59     5004           
   2 : C FLUVLH    SOURCE    BECC      97/12/23    22:15:22     3017           
   3 :       SUBROUTINE FLUVL2(NESP,
   4 :      &                  GAMG,ROG,PG,UNG,UTG,UVG,
   5 :      &                  GAMD,ROD,PD,UND,UTD,UVD,
   6 :      &                  YG,YD,FLU1,FLU2,
   7 :      &                  CELLT)
   8 : C************************************************************************
   9 : C
  10 : C PROJET            :  CASTEM 2000
  11 : C
  12 : C NOM               :  FLUVL2
  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 methode
  19 : C                      "Flux Vector Splitting" de van Leer - Hanel.  
  20 : C
  21 : C                      (voir:
  22 : C                       1)  BECCANTINI, PAILLERE,
  23 : C                           "Upwind Flux Splitting Schemes..."
  24 : C                           RAPPORT DMT 97//268
  25 : C                       2)  B. VAN LEER
  26 : C                           "Flux Vector Splittings for the 1990's"
  27 : C                           1990. Invited Lecture, NASA Lewis Research 
  28 : C                           Center
  29 : C                       3)  B. VAN LEER
  30 : C                           "Flux Vector Splitting for the Euler
  31 : C                            Equations"
  32 : C                            Report ICASE No. 82-30, 1982  )
  33 : C
  34 : C LANGUAGE          :  FORTRAN 77 
  35 : C      
  36 : C AUTEUR            :  A. BECCANTINI DRN/DMT/SEMT/TTMF  
  37 : C
  38 : C************************************************************************
  39 : C
  40 : C APPELES 
  41 : C
  42 : C FLUVL2 ------ FLPVL2
  43 : C    |
  44 : C     --------- FLMVL2     
  45 : C          
  46 : C
  47 : C************************************************************************
  48 : C
  49 : C**** Entrées:
  50 : C
  51 : C     NESP            =  nombre d'especes considérées dans les Equations
  52 : C                        d'Euler
  53 : C
  54 : C     GAMG, GAMD      =  les "gamma" du gaz (gauche et droite)
  55 : C
  56 : C     ROG, ROD        =  les densités 
  57 : C
  58 : C     PG, PD          =  les pressions
  59 : C
  60 : C     UNG, UND        =  vitesses normales
  61 : C
  62 : C     UTG, UTD        =  vitesses tangentielles
  63 : C
  64 : C     YG, YD          =  tables des fractiones massiques
  65 : C
  66 : C**** Sorties:
  67 : C
  68 : C     FLU1            =  table du flux a l'interface dans le repaire
  69 : C                        (n,t), i.e.
  70 : C                        (rho*un, rho*un*un + p, rho*un*ut, rho*un*ht,
  71 : C                         rho*un*y1, ...)
  72 : C
  73 : C     FLU2            =  table de travail, utilizé ici mais definie
  74 : C                        avant 
  75 : C
  76 : C     CELLT           =  condition de stabilité, i.e.
  77 : C
  78 : C                        dT/diamax < cellt
  79 : C      
  80 : C************************************************************************
  81 : C
  82 : C HISTORIQUE (Anomalies et modifications éventuelles)
  83 : C
  84 : C HISTORIQUE :  Créé le 6.1.98 
  85 : C
  86 : C************************************************************************
  87 : C
  88 : C N.B.: Toutes les variables sont DECLAREES
  89 : C                          
  90 : C
  91 :       IMPLICIT INTEGER(I-N)
  92 :       INTEGER NESP, I1
  93 :       REAL*8  GAMG,ROG,PG
  94 :      &       ,GAMD,ROD,PD
  95 :      &       ,GM1G,AG2,AG,UNG,UTG,UVG,MG,HTG
  96 :      &       ,GM1D,AD2,AD,UND,UTD,UVD,MD,HTD
  97 :      &       ,CELLT,CELLT0,AMG,AMD
  98 :      &       ,YG(*),YD(*),FLU1(*),FLU2(*)
  99 : C
 100 : C
 101 : C**** YG, YD, FLU1, FLU2 déjà definis avant
 102 : C
 103 : C     Dans le cas Euler monoespece, on doit
 104 : C     avoir  :
 105 : C     YG(1) = YD(1) = 0.0D0
 106 : C
 107 : C
 108 : C**** Onde de "gauche" a "droite"
 109 : C      
 110 : C     AG  = vitesse du son
 111 : C     MG  = nombre de Mach normal
 112 : C     HTG = enthalpie totale
 113 : C     
 114 :       GM1G = GAMG - 1.0D0
 115 :       AG2 = GAMG * PG/ ROG
 116 :       AG = SQRT(AG2)
 117 :       MG = UNG / AG
 118 :       AMG = ABS(MG)
 119 :       HTG = AG2 /GM1G + 0.5D0 * (UNG*UNG + UTG*UTG + UVG*UVG)
 120 :       CELLT = (2.0D0*GAMG + (AMG*(3.0D0-GAMG)))/(GAMG+3.0D0)
 121 :      &        /((AMG+1.0D0)*AG)
 122 : C     
 122 :       CALL FLPVL2(NESP,GAMG,ROG,MG,AG,UTG,UVG,HTG,YG,FLU1)
 123 : C      
 124 : C**** Onde de "droite" a "gauche".
 125 : C     
 126 : C     AD  = vitesse du son
 127 : C     MD  = nombre de Mach normal
 128 : C     HTD = enthalpie totale
 129 : C      
 130 : C
 131 : C
 132 :       GM1D = GAMD - 1.0D0
 133 :       AD2 = GAMD * PD/ ROD
 134 :       AD = SQRT(AD2)
 135 :       MD = UND / AD
 136 :       AMD = ABS(MD)
 137 :       HTD = AD2 / GM1D + 0.5D0 * (UND*UND + UTD*UTD + UVD*UVD)
 138 :       CELLT0 = (2.0D0*GAMD + (AMD*(3.0D0-GAMD)))/(GAMD+3.0D0)
 139 :      &        /((AMD+1.0D0)*AD)
 140 :       CELLT = MIN(CELLT,CELLT0)      
 141 : C     
 141 :       CALL FLMVL2(NESP,GAMD,ROD,MD,AD,UTD,UVD,HTD,YD,FLU2)
 142 : C      
 143 : C
 144 :       DO I1 = 1, NESP+5
 145 :          FLU1(I1) = FLU1(I1) + FLU2(I1) 
 146 :       ENDDO
 147 : C
 148 : C
 149 :       RETURN
 150 :       END
 151 : 
 152 : 
 153 :  
 154 :