Télécharger la plaquette Cast3M


fr en

Télécharger vlhj4.eso

Retour à la liste

Numérotation des lignes : 

vlhj4

  1. C VLHJ4     SOURCE    CHAT      05/01/13    04:07:57     5004
  2.       SUBROUTINE VLHJ4(RHO,UX,UY,UZ,P,RHOET,GAM,CNX,CNY,CNZ,
  3.      &     CT1X,CT1Y,CT1Z,CT2X,CT2Y,CT2Z,
  4.      &     DFROC,DFRUNC,DFRU1C,DFRU2C,DFRETC)
  5. C************************************************************************
  6. C
  7. C PROJET            :  CASTEM 2000
  8. C
  9. C NOM               :  VLHJ4
  10. C
  11. C DESCRIPTION       :  Calcul du Jacobien du flux numerique VLH
  12. C                      Contribution de l'onde de gauche à droite
  13. C                      Cas 3D monoespece
  14. C
  15. C LANGAGE           :  FORTRAN 77
  16. C
  17. C AUTEUR            :  A. BECCANTINI SFME/LTMF
  18. C
  19. C************************************************************************
  20. C
  21. C APPELES           :  VLHJ5
  22. C
  23. C************************************************************************
  24. C
  25. C**** Entrées:
  26. C
  27. C     RHO             =  densité
  28. C
  29. C     UX, UY, UZ      =  vitesse
  30. C
  31. C     P               =  pression
  32. C
  33. C     RHOET           =  rho * (ein + ecin)
  34. C
  35. C     GAM             =  le "gamma"
  36. C
  37. C     CNX,  CNY,  CNZ =  normale à la face
  38. C
  39. C     CT1X, CT1Y, CT1Z
  40. C     CT2X, CT2Y, CT2Z=  tangentes à la face (T1,T2)
  41. C
  42. C**** Sorties:
  43. C
  44. C     DFROC(5) = derivés partielles de F_RHO+ par raport a RHO, RHOUX,
  45. C                RHOUY, RHOUZ, RHOET
  46. C
  47. C     DFRUNC(5)= derivés partielles de F_RHOUN+ par raport a RHO, RHOUX,
  48. C                RHOUY, RHOUZ, RHOET
  49. C
  50. C     DFRU1C(5)= derivés partielles de F_RHOUT1+ par raport a RHO, RHOUX,
  51. C                 RHOUY, RHOUZ, RHOET (tangente T1)
  52. C
  53. C     DFRU2C(5)= derivés partielles de F_RHOUT2+ par raport a RHO, RHOUX,
  54. C                 RHOUY, RHOUZ, RHOET (tangente T2)
  55. C
  56. C     DFRETC(5)= derivés partielles de F_RHOET+ par raport a RHO,
  57. C                RHOUX, RHOUY, RHOUZ, RHOET
  58. C
  59. C************************************************************************
  60. C
  61. C HISTORIQUE (Anomalies et modifications éventuelles)
  62. C
  63. C HISTORIQUE :  Créé le 23.07.01
  64. C
  65. C************************************************************************
  66. C
  67. C N.B.: Toutes les variables sont declarées explicitement
  68. C
  69.       IMPLICIT INTEGER(I-N)
  70.       REAL*8 RHO,UX,UY,UZ,P,RHOET,GAM,CNX,CNY,CNZ,CT1X,CT1Y,CT1Z
  71.      &     ,CT2X,CT2Y,CT2Z
  72.      &     ,DFROC(5),DFRUNC(5),DFRU1C(5),DFRU2C(5),DFRETC(5)
  73.      &     ,DFROP(5),DFRUNP(5),DFRU1P(5),DFRU2P(5),DFRETP(5)
  74.      &     ,USRO,UXSRO,UYSRO,UZSRO,GM1,GM1EC,GM1UX,GM1UY,GM1UZ
  75. C
  76. C**** Jacobian with respect to the primitive variables
  77. C
  78.       CALL VLHJ5(RHO,UX,UY,UZ,P,RHOET,GAM,CNX,CNY,CNZ,
  79.      &     CT1X,CT1Y,CT1Z,CT2X,CT2Y,CT2Z,
  80.      &     DFROP,DFRUNP,DFRU1P,DFRU2P,DFRETP)
  81. C
  82. C**** Jacobian with respect to the conservative variables
  83. C
  84.       USRO = 1.0D0 / RHO
  85.       UXSRO = UX * USRO
  86.       UYSRO = UY * USRO
  87.       UZSRO = UZ * USRO
  88.       GM1 = GAM - 1.0D0
  89.       GM1EC = GM1 * 0.5D0 * (UX * UX + UY * UY + UZ * UZ)
  90.       GM1UX = GM1 * UX
  91.       GM1UY = GM1 * UY
  92.       GM1UZ = GM1 * UZ
  93. C
  94.       DFROC(1) = DFROP(1) - UXSRO * DFROP(2) -  UYSRO * DFROP(3)
  95.      &     - UZSRO  * DFROP(4) + GM1EC * DFROP(5)
  96.       DFROC(2) =  USRO * DFROP(2) -  GM1UX * DFROP(5)
  97.       DFROC(3) =  USRO * DFROP(3) -  GM1UY * DFROP(5)
  98.       DFROC(4) =  USRO * DFROP(4) -  GM1UZ * DFROP(5)
  99.       DFROC(5) =  GM1 * DFROP(5)
  100. C
  101.       DFRUNC(1) = DFRUNP(1) - UXSRO * DFRUNP(2) -  UYSRO * DFRUNP(3)
  102.      &     - UZSRO  * DFRUNP(4) + GM1EC * DFRUNP(5)
  103.       DFRUNC(2) =  USRO * DFRUNP(2) -  GM1UX * DFRUNP(5)
  104.       DFRUNC(3) =  USRO * DFRUNP(3) -  GM1UY * DFRUNP(5)
  105.       DFRUNC(4) =  USRO * DFRUNP(4) -  GM1UZ * DFRUNP(5)
  106.       DFRUNC(5) =  GM1 * DFRUNP(5)
  107. C
  108.       DFRU1C(1) = DFRU1P(1) - UXSRO * DFRU1P(2) -  UYSRO * DFRU1P(3)
  109.      &     - UZSRO  * DFRU1P(4) + GM1EC * DFRU1P(5)
  110.       DFRU1C(2) =  USRO * DFRU1P(2) -  GM1UX * DFRU1P(5)
  111.       DFRU1C(3) =  USRO * DFRU1P(3) -  GM1UY * DFRU1P(5)
  112.       DFRU1C(4) =  USRO * DFRU1P(4) -  GM1UZ * DFRU1P(5)
  113.       DFRU1C(5) =  GM1 * DFRU1P(5)
  114. C
  115.       DFRU2C(1) = DFRU2P(1) - UXSRO * DFRU2P(2) -  UYSRO * DFRU2P(3)
  116.      &     - UZSRO  * DFRU2P(4) + GM1EC * DFRU2P(5)
  117.       DFRU2C(2) =  USRO * DFRU2P(2) -  GM1UX * DFRU2P(5)
  118.       DFRU2C(3) =  USRO * DFRU2P(3) -  GM1UY * DFRU2P(5)
  119.       DFRU2C(4) =  USRO * DFRU2P(4) -  GM1UZ * DFRU2P(5)
  120.       DFRU2C(5) =  GM1 * DFRU2P(5)
  121. C
  122.       DFRETC(1) = DFRETP(1) - UXSRO * DFRETP(2) -  UYSRO * DFRETP(3)
  123.      &     - UZSRO  * DFRETP(4) + GM1EC * DFRETP(5)
  124.       DFRETC(2) =  USRO * DFRETP(2) -  GM1UX * DFRETP(5)
  125.       DFRETC(3) =  USRO * DFRETP(3) -  GM1UY * DFRETP(5)
  126.       DFRETC(4) =  USRO * DFRETP(4) -  GM1UZ * DFRETP(5)
  127.       DFRETC(5) =  GM1 * DFRETP(5)
  128. C
  129.       RETURN
  130.       END
  131.  
  132.  
  133.  
  134.  
  135.  

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