Télécharger la plaquette Cast3M


fr en

Télécharger vlhj2.eso

Retour à la liste

Numérotation des lignes : 

vlhj2

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

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