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

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