Télécharger la plaquette Cast3M


fr en

Télécharger vlhj3.eso

Retour à la liste

Numérotation des lignes : 

vlhj3

  1. C VLHJ3     SOURCE    CHAT      05/01/13    04:07:53     5004
  2.       SUBROUTINE VLHJ3(RHO,UX,UY,P,GAM,CNX,CNY,DFRUN)
  3. C
  4. C************************************************************************
  5. C
  6. C PROJET            :  CASTEM 2000
  7. C
  8. C NOM               :  VLHJ3
  9. C
  10. C DESCRIPTION       :  Calcul du Jacobien du flux numerique VLH
  11. C                      Contribution des faces sur les murs
  12. C                      Voir VLHJ2
  13. C
  14. C LANGAGE           :  FORTRAN 77
  15. C
  16. C AUTEUR            :  A. BECCANTINI DRN/DMT/SEMT/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   =  vitesse
  29. C
  30. C     P        =  pression
  31. C
  32. C     GAM      =  le "gamma"
  33. C
  34. C     CNX, CNY =  normale à la face
  35. C
  36. C
  37. C**** Sorties:
  38. C
  39. C     DFRUN(4)= derivés partielles de F_RHOUN par raport a RHO, UX,
  40. C                UY, P
  41. C
  42. C************************************************************************
  43. C
  44. C HISTORIQUE (Anomalies et modifications éventuelles)
  45. C
  46. C HISTORIQUE :  Créé le 04.09.00
  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,P,GAM,CNX,CNY
  54.      &     ,DFRUN(4)
  55.      &     ,UN,PSRHO,CSON2,CSON
  56.      &     ,MACH,MNUM,DERMNU,PNUM,DERPNU
  57.      &     ,CELL
  58. C
  59.       UN = UX * CNX + UY * CNY
  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) = (CSON * MNUM - UN * DERMNU) * UN / PSRHO +
  72.      &     2.0D0 * PNUM - MACH * DERPNU
  73. C
  74.       RETURN
  75.       END
  76.  
  77.  
  78.  
  79.  

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