Télécharger la plaquette Cast3M


fr en

Télécharger vlhj6.eso

Retour à la liste

Numérotation des lignes : 

vlhj6

  1. C VLHJ6     SOURCE    CHAT      05/01/13    04:08:04     5004
  2.       SUBROUTINE VLHJ6(RHO,UX,UY,UZ,P,GAM,CNX,CNY,CNZ,
  3.      &                  DFRUNC)
  4. C************************************************************************
  5. C
  6. C PROJET            :  CASTEM 2000
  7. C
  8. C NOM               :  VLHJ6
  9. C
  10. C DESCRIPTION       :  Calcul du Jacobien du flux numerique VLH
  11. C                      Contribution des faces sur les murs
  12. C                      Cas 3D
  13. C
  14. C LANGAGE           :  FORTRAN 77
  15. C
  16. C AUTEUR            :  A. BECCANTINI SFME/LTMF
  17. C
  18. C************************************************************************
  19. C
  20. C APPELES           :  VLHJ3
  21. C
  22. C************************************************************************
  23. C
  24. C**** Entrées:
  25. C
  26. C     RHO           =  densité
  27. C
  28. C     UX, UY, UZ    =  vitesse
  29. C
  30. C     P             =  pression
  31. C
  32. C     GAM           =  le "gamma"
  33. C
  34. C     CNX, CNY, CNZ =  normales à la face
  35. C
  36. C
  37. C**** Sorties:
  38. C
  39. C     DFRUNC(4)= derivés partielles de F_RHOUN par raport a RHO, RHOUX,
  40. C                RHOUY, RHOUZ, RHOET
  41. C
  42. C************************************************************************
  43. C
  44. C HISTORIQUE (Anomalies et modifications éventuelles)
  45. C
  46. C HISTORIQUE :  Créé le 23.07.01
  47. C
  48. C************************************************************************
  49. C
  50. C N.B.: Toutes les variables sont declarées explicitement
  51. C
  52. C
  53.       IMPLICIT INTEGER(I-N)
  54.       REAL*8 RHO,UX,UY,UZ,P,GAM,CNX,CNY,CNZ
  55.      &      ,DFRUNC(5),DFRUNP(5)
  56.      &      ,USRO,UXSRO,UYSRO,UZSRO,GM1,GM1EC,GM1UX,GM1UY,GM1UZ
  57. C
  58. C**** Jacobian with respect to the primitive variables
  59. C
  60.       CALL VLHJ7(RHO,UX,UY,UZ,P,GAM,CNX,CNY,CNZ,DFRUNP)
  61. C
  62. C**** Jacobian with respect to the conservative variables
  63. C
  64.       USRO = 1.0D0 / RHO
  65.       UXSRO = UX * USRO
  66.       UYSRO = UY * USRO
  67.       UZSRO = UZ * USRO
  68.       GM1 = GAM - 1.0D0
  69.       GM1EC = GM1 * 0.5D0 * (UX * UX + UY * UY + UZ * UZ)
  70.       GM1UX = GM1 * UX
  71.       GM1UY = GM1 * UY
  72.       GM1UZ = GM1 * UZ
  73. C
  74.       DFRUNC(1) = DFRUNP(1) - UXSRO * DFRUNP(2) -  UYSRO * DFRUNP(3)
  75.      &           - UZSRO * DFRUNP(4) + GM1EC * DFRUNP(5)
  76.       DFRUNC(2) =  USRO * DFRUNP(2) -  GM1UX * DFRUNP(5)
  77.       DFRUNC(3) =  USRO * DFRUNP(3) -  GM1UY * DFRUNP(5)
  78.       DFRUNC(4) =  USRO * DFRUNP(4) -  GM1UZ * DFRUNP(5)
  79.       DFRUNC(5) =  GM1 * DFRUNP(5)
  80. C
  81.       RETURN
  82.       END
  83.  
  84.  
  85.  
  86.  
  87.  

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