Télécharger la plaquette Cast3M


fr en

Télécharger frusa3.eso

Retour à la liste

Numérotation des lignes : 

frusa3

  1. C FRUSA3    SOURCE    KLOCZKO   05/06/14    21:15:11     5111
  2. C FRUSA3    SOURCE    CHAT      05/01/13    00:09:15     5004
  3.       SUBROUTINE FRUSA3(NESP,
  4.      &                  GAMG,ROG,PG,UNG,UTG,UVG,
  5.      &                  GAMD,ROD,PD,UND,UTD,UVD,
  6.      &                  YG,YD,FLU,
  7.      &                  CELLT)
  8. C
  9. C
  10. C PROJET            :  CASTEM 2000
  11. C
  12. C NOM               :  FRUSA3
  13. C
  14. C DESCRIPTION       :  Formulation Volumes Finis  pour les Equations
  15. C                      d'Euler Multi-Especes relatives à un melange
  16. C                      de gaz parfaits.
  17. C
  18. C                      Calcul du flux aux interfaces avec la methode
  19. C                      de Rusanov
  20. C
  21. C LANGUAGE          :  FORTRAN 77
  22. C
  23. C AUTEUR            :  T. KLOCZKO DM2S/SFME/LTMF
  24. C
  25. C************************************************************************
  26. C
  27. C APPELES           :
  28. C
  29. C************************************************************************
  30. C
  31. C**** Entrées:
  32. C
  33. C     NESP            =  nombre d'especes considérées dans les Equations
  34. C                        d'Euler
  35. C
  36. C     GAMG, GAMD      =  les "gamma" du gaz (gauche et droite)
  37. C
  38. C     ROG, ROD        =  les densités
  39. C
  40. C     PG, PD          =  les pressions
  41. C
  42. C     UNG, UND        =  vitesses normales
  43. C
  44. C     UTG, UTD        =  vitesses tangentielles
  45. C
  46. C     UVG, UVD        =  vitesses tangentielles
  47. C
  48. C**** Sorties:
  49. C
  50. C     FLU             =  table du flux a l'interface dans le repaire
  51. C                        (n,t), i.e.
  52. C                        (rho*un, rho*un*un + p, rho*un*ut,  rho*un*uv,
  53. C                         rho*un*ht, rho*un*y1, ...)
  54. C
  55. C     CELLT           =  condition de stabilité, i.e.
  56. C
  57. C                        dT/diamax < cellt
  58. C
  59. C************************************************************************
  60. C
  61. C HISTORIQUE (Anomalies et modifications éventuelles)
  62. C
  63. C HISTORIQUE :  Créé le 02.05.05
  64. C
  65. C
  66. C
  67. C************************************************************************
  68. C
  69. C N.B.: Toutes les variables sont DECLAREES
  70. C
  71. C
  72.       IMPLICIT INTEGER(I-N)
  73.       INTEGER NESP,I1
  74.       REAL*8  GAMG,ROG,PG,UNG,UTG,UVG,RETG
  75.      &       ,GAMD,ROD,PD,UND,UTD,UVD,RETD
  76.      &       ,YG(*),YD(*)
  77.      &       ,FLU(*),CELLT, FLUM
  78.      &       ,GAMF,ROF,PF,UNF,CF
  79.      &       ,VPMAX
  80.      &       ,XP(5),FC(5)
  81. C
  82. C*** Calcul des énergies
  83. C
  84.       RETG = PG / (GAMG - 1.0D0)
  85.      &     + 0.5D0 * ROG * (UNG**2 + UTG**2 + UVG**2)
  86.       RETD = PD / (GAMD - 1.0D0)
  87.      &     + 0.5D0 * ROD * (UND**2 + UTD**2 + UVD**2)
  88. C
  89. C*** Calcul des états moyens
  90. C
  91.       UNF  = 0.5D0 * (UNG  + UND)
  92.       ROF  = 0.5D0 * (ROG  + ROD)
  93.       PF   = 0.5D0 * (PG   + PD)
  94.       GAMF = 0.5D0 * (GAMG + GAMD)
  95.       CF    = (GAMF * PF / ROF)**0.5D0
  96. C
  97. C*** Calcul du rayon spectral et de la condition de stabilité
  98. C
  99.       VPMAX = ABS(UNF) + CF
  100.       CELLT = 1.0D0 / VPMAX
  101. C
  102. C*** Calcul des deltas
  103. C
  104.       XP(1) = ROG - ROD
  105.       XP(2) = (ROG * UNG) - (ROD * UND)
  106.       XP(3) = (ROG * UTG) - (ROD * UTD)
  107.       XP(4) = (ROG * UVG) - (ROD * UVD)
  108.       XP(5) = RETG - RETD
  109. C
  110. C*** Calcul du flux convectif
  111. C
  112.       FC(1) = 0.5D0 * ((ROG * UNG)            + (ROD * UND))
  113.       FC(2) = 0.5D0 * ((ROG * UNG * UNG + PG) + (ROD * UND * UND + PD))
  114.       FC(3) = 0.5D0 * ((ROG * UNG * UTG)      + (ROD * UND * UTD))
  115.       FC(4) = 0.5D0 * ((ROG * UNG * UVG)      + (ROD * UND * UVD))
  116.       FC(5) = 0.5D0 * ((UNG * (RETG + PG))    + (UND * (RETD + PD)))
  117. C
  118. C*** Calcul du flux numérique
  119. C
  120.       DO I1 = 1,5
  121.          FLU(I1) = FC(I1) + 0.5D0 * VPMAX * XP(I1)
  122.       ENDDO
  123. C
  124. C*** Partie multi-espèces
  125. C
  126.       FLUM = FLU(1)
  127.       IF(FLUM .GT. 0)THEN
  128.          DO I1 = 1, NESP, 1
  129.             FLU(5+I1) = FLUM * YG(I1)
  130.          ENDDO
  131.       ELSE
  132.          DO I1 = 1, NESP, 1
  133.             FLU(5+I1) = FLUM * YD(I1)
  134.          ENDDO
  135.       ENDIF
  136. C
  137.       RETURN
  138.       END
  139.  
  140.  
  141.  
  142.  
  143.  
  144.  
  145.  
  146.  
  147.  
  148.  
  149.  

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