Télécharger la plaquette Cast3M


fr en

Télécharger pret.eso

Retour à la liste

Numérotation des lignes : 

pret

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

© Cast3M 2003 - Tous droits réservés.
Mentions légales