Télécharger la plaquette Cast3M


fr en

Télécharger vlhj7.eso

Retour à la liste

Numérotation des lignes : 

vlhj7

  1. C VLHJ7     SOURCE    CHAT      05/01/13    04:08:07     5004
  2.       SUBROUTINE VLHJ7(RHO,UX,UY,UZ,P,GAM,CNX,CNY,CNZ,DFRUN)
  3. C
  4. C************************************************************************
  5. C
  6. C PROJET            :  CASTEM 2000
  7. C
  8. C NOM               :  VLHJ7
  9. C
  10. C DESCRIPTION       :  Calcul du Jacobien du flux numerique VLH
  11. C                      Contribution des faces sur les murs
  12. C                      Voir VLHJ6
  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           :  MPFUNT
  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 =  normale à la face
  35. C
  36. C
  37. C**** Sorties:
  38. C
  39. C     DFRUN(5)= derivés partielles de F_RHOUN par raport a RHO, UX,
  40. C                UY, UZ, P
  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.       IMPLICIT INTEGER(I-N)
  53.       REAL*8 RHO,UX,UY,UZ,P,GAM,CNX,CNY,CNZ
  54.      &     ,DFRUN(5)
  55.      &     ,UN,PSRHO,CSON2,CSON
  56.      &     ,MACH,MNUM,DERMNU,PNUM,DERPNU
  57.      &     ,CELL
  58. C
  59.       UN = UX * CNX + UY * CNY + UZ * CNZ
  60.       PSRHO = P / RHO
  61.       CSON2 = GAM * PSRHO
  62.       CSON = SQRT(CSON2)
  63.       MACH = UN / CSON
  64.       CALL MPFUNT(MACH,MNUM,DERMNU,PNUM,DERPNU)
  65. C
  66.       CELL = DERMNU * UN + MNUM * CSON
  67.       DFRUN(1) = CELL * UN + DERPNU * MACH * PSRHO
  68.       CELL = (CELL *  RHO + P / CSON * DERPNU) * 2.0D0
  69.       DFRUN(2) = CELL * CNX
  70.       DFRUN(3) = CELL * CNY
  71.       DFRUN(4) = CELL * CNZ
  72.       DFRUN(5) = (CSON * MNUM - UN * DERMNU) * UN / PSRHO +
  73.      &     2.0D0 * PNUM - MACH * DERPNU
  74. C
  75.       RETURN
  76.       END
  77.  
  78.  
  79.  
  80.  
  81.  

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