Télécharger pret.eso

Retour à la liste

Numérotation des lignes :

  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
  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. C
  102. C**** Lecture de l'option
  103. C
  104. ICOND = 1
  105. CALL LIRCHA(MOTOP,ICOND,IRET)
  106. IF(IERR .NE. 0) GOTO 9999
  107. NBOPT0 = NBOPT
  108. CALL OPTLI(IP,TOPTIO,MOTOP,NBOPT0)
  109. IF(IP .EQ. 1) THEN
  110. C
  111. C******* Gaz 'calorically perfect' monoespece
  112. C
  113. CALL PRE1()
  114. ELSEIF(IP .EQ. 2)THEN
  115. C
  116. C******* Gaz 'calorically perfect' multiespeces
  117. C
  118. CALL PRE2()
  119. ELSEIF(IP .EQ. 3)THEN
  120. C
  121. C******* Gaz 'thermally perfect' multiespeces
  122. C
  123. CALL PRE3()
  124. ELSEIF(IP .EQ. 4)THEN
  125. C
  126. C******* Transport des scalaires passifs
  127. C
  128. CALL PRE4()
  129. ELSEIF(IP .EQ. 5)THEN
  130. C
  131. C******* Modele bifluide
  132. C
  133. CALL PRE5()
  134. ELSEIF(IP .EQ. 6)THEN
  135. C
  136. C******* Discrete Equation Method
  137. C
  138. CALL PRE6()
  139. ELSEIF(IP .EQ. 7)THEN
  140. C
  141. C******* GFMP
  142. C
  143. CALL PRE7()
  144. ELSE
  145. C
  146. C******* Message d'erreur standard
  147. C 251 2
  148. C Tentative d'utilisation d'une option non implémentée
  149. C
  150. CALL ERREUR(251)
  151. ENDIF
  152. C
  153. 9999 CONTINUE
  154. RETURN
  155. END
  156.  
  157.  
  158.  
  159.  
  160.  
  161.  
  162.  
  163.  
  164.  
  165.  
  166.  
  167.  
  168.  
  169.  
  170.  

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