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