Télécharger la plaquette Cast3M


fr en

Télécharger vlhj1.eso

Retour à la liste

Numérotation des lignes : 

vlhj1

  1. C VLHJ1     SOURCE    CHAT      05/01/13    04:07:48     5004
  2.       SUBROUTINE VLHJ1(RHO,UX,UY,P,RHOET,GAM,CNX,CNY,CTX,CTY,
  3.      &                  DFRO,DFRUN,DFRUT,DFRET)
  4. C************************************************************************
  5. C
  6. C PROJET            :  CASTEM 2000
  7. C
  8. C NOM               :  VLHJ1
  9. C
  10. C DESCRIPTION       :  Calcul du Jacobien du flux numerique VLH
  11. C                      Contribution de l'onde de gauche à droite
  12. C                      Voir VLHJ0
  13. C
  14. C LANGAGE           :  FORTRAN 77
  15. C
  16. C AUTEUR            :  A. BECCANTINI DRN/DMT/SEMT/LTMF
  17. C
  18. C************************************************************************
  19. C
  20. C APPELES           :  MPFUNT
  21. C
  22. C************************************************************************
  23. C
  24. C**** Entrées:
  25. C
  26. C     RHO      =  densité
  27. C
  28. C     UX, UY   =  vitesse
  29. C
  30. C     P        =  pression
  31. C
  32. C     RHOET    =  rho * (ein + ecin)
  33. C
  34. C     GAM      =  le "gamma"
  35. C
  36. C     CNX, CNY =  normale à la face
  37. C
  38. C     CTX, CTY =  tangente à la face
  39. C
  40. C
  41. C**** Sorties:
  42. C
  43. C     DFRO(4)  = derivés partielles de F_RHO+ par raport a RHO, UX,
  44. C                UY, P
  45. C
  46. C     DFRUN(4) = derivés partielles de F_RHOUN+ par raport a RHO, UX,
  47. C                UY, P
  48. C
  49. C     DFRUT(4) = derivés partielles de F_RHOUT+ par raport a RHO, UX,
  50. C                UY, P
  51. C
  52. C     DFRET(4) = derivés partielles de F_RHOET+ par raport a RHO,
  53. C                UX, UY, P
  54. C
  55. C************************************************************************
  56. C
  57. C HISTORIQUE (Anomalies et modifications éventuelles)
  58. C
  59. C HISTORIQUE :  Créé le 04.09.00
  60. C
  61. C************************************************************************
  62. C
  63. C N.B.: Toutes les variables sont declarées explicitement
  64. C
  65. C
  66.       IMPLICIT INTEGER(I-N)
  67.       REAL*8 RHO,UX,UY,P,RHOET,GAM,CNX,CNY,CTX,CTY
  68.      &      ,DFRO(4),DFRUN(4),DFRUT(4),DFRET(4)
  69.      &      ,UN,UT,CSON2,CSON,HT
  70.      &      ,MACH,MNUM,DERMNU,PNUM,DERPNU
  71.      &      ,CELL, FRHO
  72. C
  73.       UN = UX * CNX + UY * CNY
  74.       UT = UX * CTX + UY * CTY
  75.       CSON2 = GAM * P / RHO
  76.       CSON = SQRT(CSON2)
  77.       MACH = UN / CSON
  78.       HT = (RHOET + P) / RHO
  79.       CALL MPFUNT(MACH,MNUM,DERMNU,PNUM,DERPNU)
  80. C
  81.       FRHO = MNUM * RHO * CSON
  82. C
  83.       DFRO(1) = ((DERMNU * UN) + (MNUM * CSON))  / 2.0D0
  84.       DFRO(2) = DERMNU * RHO * CNX
  85.       DFRO(3) = DERMNU * RHO * CNY
  86.       DFRO(4) = (MNUM - DERMNU * MACH) * RHO * CSON / (2.0D0 * P)
  87. C
  88.       DFRUN(1) = (DFRO(1) * UN) + (DERPNU * MACH * P / (2.0D0 * RHO))
  89.       CELL = (MNUM * RHO * CSON) + (DERPNU  * P / CSON)
  90.       DFRUN(2) = (DFRO(2) * UN) + (CELL * CNX)
  91.       DFRUN(3) = (DFRO(3) * UN) + (CELL * CNY)
  92.       DFRUN(4) = (DFRO(4) * UN) - (0.5D0 * MACH * DERPNU) + PNUM
  93. C
  94.       DFRUT(1) = DFRO(1) * UT
  95.       DFRUT(2) = (DFRO(2) * UT)  + (FRHO * CTX)
  96.       DFRUT(3) = (DFRO(3) * UT)  + (FRHO * CTY)
  97.       DFRUT(4) = DFRO(4) * UT
  98. C
  99.       CELL = CSON2 / (GAM - 1.0D0)
  100.       DFRET(1) = (DFRO(1) * HT) - (FRHO * CELL / RHO)
  101.       DFRET(2) = (DFRO(2) * HT)  + (FRHO * UX)
  102.       DFRET(3) = (DFRO(3) * HT)  + (FRHO * UY)
  103.       DFRET(4) = (DFRO(4) * HT) + (FRHO * CELL / P)
  104. C
  105.       RETURN
  106.       END
  107.  
  108.  
  109.  
  110.  

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