Télécharger clim2.eso

Retour à la liste

Numérotation des lignes :

  1. C CLIM2 SOURCE KK2000 14/04/10 21:15:05 8032
  2. SUBROUTINE CLIM2
  3. C************************************************************************
  4. C
  5. C PROJET : CASTEM 2000
  6. C
  7. C NOM : CLIM2
  8. C
  9. C DESCRIPTION : Subroutine appellée par KON2
  10. C
  11. C Modelisation 2D/3D des equations d'Euler
  12. C Conditions aux limites
  13. C
  14. C Calcul du flux/residu pour le gaz multiespeces
  15. C parfait aux limites
  16. C LANGAGE : FORTRAN 77 + ESOPE 2000 (avec estensions CISI)
  17. C
  18. C AUTEUR : S. Kudriakov, DEN/DM2S/SFME/LTMF
  19. C
  20. C************************************************************************
  21. C
  22. C*** SYNTAXE
  23. C
  24. C Discrétisation en VF "cell-centered" des équations d'Euler pour
  25. C un gaz parfait multiespeces
  26. C
  27. C
  28. C RCHPLI RCHPRE = 'KONV' 'VF' 'PERFMULT' 'CLIM'
  29. C 'RESI' $MOD1 $MOD2 TABGAS
  30. C LMOTC LMOTP CHPRN CHPVN CHPPN CHPYN
  31. C MOT1 CHPLI ;
  32. C
  33. C $MOD1 : l'objet modele du domaine interieur
  34. C ($TAB1 est sa table domaine)
  35. C
  36. C $MOD2 : l'objet modele du domaine du bord
  37. C ($TAB2 est sa table domaine)
  38. C
  39. C TABGAS : the table which containes information on
  40. C Cp and Cv of each of the gas
  41. C
  42. C LMOTC : LISTMOTS, noms des variables conservatives
  43. C
  44. C LMOTP : LISTMOTS, noms des variables primitives
  45. C
  46. C CHPRN : densité (SPG = 'CENTRE', une seule
  47. C composante, 'SCAL')
  48. C
  49. C CHPVN : vitesse (SPG = $TAB1 . 'CENTRE',
  50. C composantes: 'UX', 'UY', ('UZ'))
  51. C
  52. C CHPPN : pression (SPG = $TAB1 . 'CENTRE', une seule
  53. C composante, 'SCAL')
  54. C
  55. C CHPYN : fractions massiques (SPG = $TAB1 . 'CENTRE',
  56. C number of components is equal to the number of
  57. C components in (TABGAS . 'ESPEULE'))
  58. C
  59. C MOT1 : MOT, type de condition limite:
  60. C 'INRI' : inlet, invariants de Riemann
  61. C
  62. C CHPLI : condition limit imposé (SPG = $TAB2 . 'CENTRE')
  63. C Les composantes dependent de MOT1
  64. C 'INRI' : on donne densité 'RN', vitesse 'UX' 'UY' ('UX')
  65. C pression 'PN'
  66. C 'OUTRI': on donne densité 'RN', vitesse 'UX' 'UY' ('UX')
  67. C pression 'PN'
  68. C 'INSS' : on donne densité 'RN', vitesse 'UX' 'UY' ('UX')
  69. C pression 'PN'
  70. C 'OUTSS': chpoint vide
  71. C 'OUTP' : on donne la pression 'PN', on recupere la
  72. C vitesse et la densité à l'interieur
  73. C 'INSU' : on donne l'enthalpie totale 'HT', l'entropie 'S',
  74. C on impose que la vitesse tangentielle est nulle,
  75. C on recupere la pression à l'interieur
  76. C 'INJE' : on donne le flux de masse et RT (temperature fois
  77. C la constante de gaz); on impose que la vitesse
  78. C tangentielle est nulle, on recupere la pression.
  79. C
  80. C SORTIES
  81. C
  82. C RCHPLI : on donne densité, vitesse, pression et les
  83. C fractions massiques sur le bord
  84. C (SPG = $TAB2 . 'CENTRE', composantes en LMOTP)
  85. C
  86. C RCHPRE : residu (SPG = $TAB1 . 'CENTRE', composantes en LMOTC)
  87. C (of the cells which are attachet to the border)
  88. C************************************************************************
  89. C
  90. C HISTORIQUE (Anomalies et modifications éventuelles)
  91. C
  92. C HISTORIQUE :
  93. C
  94. C************************************************************************
  95. IMPLICIT INTEGER(I-N)
  96. IMPLICIT REAL*8(A-H,O-Z)
  97.  
  98. -INC CCOPTIO
  99. C
  100. INTEGER NBOPT, ICELL, IJAC
  101. C
  102. PARAMETER (NBOPT=3)
  103. CHARACTER*8 LOPT(NBOPT)
  104. C
  105. DATA LOPT/'RESI ','JACOCONS','JACOPRIM'/
  106. C
  107. CALL LIRMOT(LOPT,NBOPT,ICELL,1)
  108. IF(IERR .NE. 0) GOTO 9999
  109. IF(ICELL .EQ. 1)THEN
  110. C 'RESI '
  111. IJAC=0
  112. CALL CLIM22(IJAC)
  113. IF(IERR .NE. 0) GOTO 9999
  114. ELSEIF(ICELL .EQ. 2)THEN
  115. C 'JACOCONS'
  116. IJAC=1
  117. CALL CLIM22(IJAC)
  118. IF(IERR .NE. 0) GOTO 9999
  119. ELSEIF(ICELL .EQ. 3)THEN
  120. C 'JACOPRIM'
  121. IJAC=2
  122. CALL CLIM22(IJAC)
  123. IF(IERR .NE. 0) GOTO 9999
  124. ENDIF
  125. 9999 CONTINUE
  126. RETURN
  127. END
  128.  
  129.  
  130.  
  131.  
  132.  
  133.  
  134.  
  135.  
  136.  
  137.  

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