Télécharger la plaquette Cast3M


fr en

Télécharger vlhj5.eso

Retour à la liste

Numérotation des lignes : 

vlhj5

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

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