Esope source files

   1 : C PRET      SOURCE    BECC      11/05/18    21:15:32     6974           
   2 :       SUBROUTINE PRET()
   3 : C************************************************************************
   4 : C
   5 : C PROJET            :  CASTEM 2000
   6 : C
   7 : C NOM               :  PRET (OPERATEUR GIBIANE)
   8 : C
   9 : C DESCRIPTION       :  Resolution des Equations d'Euler et des
  10 : C                      Equations de Ecoulement Bifluide, VF,
  11 : C                     "cell-centred approach"
  12 : C      
  13 : C                      Calcul des etats "gauche et droite" aux  interfaces
  14 : C                      pour les suivants models de gaz
  15 : C    
  16 : C                      1) gaz ideal mono-espece
  17 : C
  18 : C                      2) gaz ideal multi-especes 
  18 : 

Warning:  preg_replace(): Compilation failed: unmatched parentheses at offset 13 in C:\xampp\htdocs\cast3m_kit_local_fr_en_2014\sources.php on line 147

  19 : C
  20 : C                      3) gaz ideal  
  21 : C                         ("thermally-perfect gas", avec cv_i=cv_i(T))
  22 : C
  23 : C                      5) modele bifluide (two-fluid flow, water + air)
  24 : C
  25 : C                      6) Discrete Equation Method for combustion
  26 : C
  27 : C                      7) Ghost Fluid Method for the poors GFMP 
  28 : C
  29 : C                      Reconstruction en espace  sur les variables
  30 : C                      primitives densite, pression, vitesse (fractions
  31 : C                      massiques)
  32 : C
  33 : C                      Predicteur-Correcteur d' Hancock
  34 : C                      ("MUSCL approach" de Van-Leer)
  35 : C      
  36 : C                      Voir
  37 : C      
  38 : C                      1) Beccantini, Paillere
  39 : C                           "Upwind Flux Splitting Schemes..."
  40 : C                           RAPPORT DMT 97//268
  41 : C
  42 : C                      2) Anderson, "Hypersonic and High Temperature 
  43 : C                         Gas Dynamics", McGraw-Hill, 1989
  44 : C
  45 : C                      Transport des scalaires passifs en VF
  46 : C
  47 : C      
  48 : C LANGAGE           :  FORTRAN 77 + ESOPE 2000 (avec extensions CISI)
  49 : C      
  50 : C AUTEUR            :  A. BECCANTINI, DRN/DMT/SEMT/TTMF                        
  51 : C
  52 : C************************************************************************
  53 : C
  54 : C
  55 : C APPELES (Outils)  :  LIRCHA, OPTLI, ERREUR
  56 : C
  57 : C APPELES (Calcul)  :
  58 : C
  59 : C PRET ---------------------- PRE1* : gaz monoespece (cp = const)
  60 : C   |
  61 : C   |
  62 : C    ------------------------ PRE2* : gaz multiespeces (cp_i = cost.)
  63 : C   |
  64 : C   |
  65 : C    ------------------------ PRE3* : gaz avec (cp_i = cp_i(T))
  66 : C   |
  67 : C   |
  68 : C    ------------------------ PRE4* : transport des scalaires passifs
  69 : C   |
  70 : C   |
  71 : C    ------------------------ PRE5* : melange bifluide (air et eau)
  72 : C   |
  73 : C   |
  74 : C    ------------------------ PRE6* : Discrete Equation Method 
  75 : C   |
  76 : C   |
  77 : C    ------------------------ PRE7* : GFMP
  78 : C
  79 : C************************************************************************
  80 : C
  81 : C HISTORIQUE (Anomalies et modifications éventuelles)
  82 : C
  83 : C HISTORIQUE :  Créée le 10.7.98.
  84 : C 
  85 : C************************************************************************
  86 : C
  87 :       IMPLICIT INTEGER(I-N)
  88 :       INTEGER IP,IRET, ICOND,NBOPT,NBOPT0
  89 :       PARAMETER(NBOPT=7)
  90 :       CHARACTER*(8) TOPTIO(NBOPT),MOTOP
  91 : C
  92 :       DATA TOPTIO/'PERFMONO','PERFMULT','PERFTEMP','CLAUDEIS',
  93 :      &            'TWOFLUID','DEM     ','GFMP    '/
  94 : C
  95 : C**** Les Includes
  96 : C
  97 : -INC CCOPTIO       
  98 : C
  99 : C**** Lecture de l'option 
 100 : C
 101 :       ICOND = 1
 101 :       CALL LIRCHA(MOTOP,ICOND,IRET)
 102 :       IF(IERR .NE. 0) GOTO 9999
 103 :       NBOPT0 = NBOPT
 103 :       CALL OPTLI(IP,TOPTIO,MOTOP,NBOPT0)
 104 :       IF(IP .EQ. 1) THEN
 105 : C
 105 : C******* Gaz 'calorically perfect' monoespece
 106 : C
 106 :          CALL PRE1()
 107 :       ELSEIF(IP .EQ. 2)THEN
 108 : C
 108 : C******* Gaz 'calorically perfect' multiespeces
 109 : C
 109 :          CALL PRE2()
 110 :       ELSEIF(IP .EQ. 3)THEN
 111 : C
 112 : C******* Gaz 'thermally perfect' multiespeces
 113 : C
 113 :          CALL PRE3()
 114 :       ELSEIF(IP .EQ. 4)THEN
 115 : C
 116 : C******* Transport des scalaires passifs
 117 : C
 117 :          CALL PRE4()
 118 :       ELSEIF(IP .EQ. 5)THEN
 119 : C
 120 : C******* Modele bifluide
 121 : C
 121 :          CALL PRE5()
 122 :       ELSEIF(IP .EQ. 6)THEN
 123 : C
 124 : C******* Discrete Equation Method
 125 : C
 125 :          CALL PRE6()
 126 :       ELSEIF(IP .EQ. 7)THEN
 127 : C
 128 : C******* GFMP
 129 : C
 129 :          CALL PRE7()
 130 :       ELSE
 131 : C
 132 : C******* Message d'erreur standard
 133 : C        251 2
 134 : C        Tentative d'utilisation d'une option non implémentée
 135 : C         
 135 :          CALL ERREUR(251)
 136 :       ENDIF
 137 : C
 138 :  9999 CONTINUE
 139 :       RETURN
 140 :       END
 141 : 
 142 : 
 143 : 
 144 : 
 145 : 
 146 :  
 147 :  
 148 :  
 149 :  
 150 :  
 151 :  
 152 :  
 153 :  
 154 :