Télécharger centj2.eso

Retour à la liste

Numérotation des lignes :

  1. C CENTJ2 SOURCE CHAT 05/01/12 21:52:05 5004
  2. SUBROUTINE CENTJ2(RHO,UX,UY,GAM,CNX,CNY,
  3. & DFRUNC)
  4. C************************************************************************
  5. C
  6. C PROJET : CASTEM 2000
  7. C
  8. C NOM : CENTJ2
  9. C
  10. C DESCRIPTION : Calcul du Jacobien du flux numerique centré
  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 : CENTJ3
  20. C Appellée par KONJR1, KONJC1
  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 GAM = le "gamma"
  31. C
  32. C CNX, CNY = normale à la face
  33. C
  34. C
  35. C**** Sorties:
  36. C
  37. C DFRUNC(4)= derivés partielles de F_RHOUN par raport a RHO, RHOUX,
  38. C RHOUY, RHOET
  39. C
  40. C************************************************************************
  41. C
  42. C HISTORIQUE (Anomalies et modifications éventuelles)
  43. C
  44. C HISTORIQUE : Créé le 04.09.00
  45. C
  46. C************************************************************************
  47. C
  48. C N.B.: Toutes les variables sont declarées explicitement
  49. C
  50. C
  51. IMPLICIT INTEGER(I-N)
  52. REAL*8 RHO,UX,UY,GAM,CNX,CNY
  53. & ,DFRUNC(4),DFRUNP(4)
  54. & ,USRO,UXSRO,UYSRO,GM1,GM1EC,GM1UX,GM1UY
  55. C
  56. C**** Jacobian with respect to the primitive variables
  57. C
  58. CALL CENTJ3(RHO,UX,UY,CNX,CNY,DFRUNP)
  59. C
  60. C**** Jacobian with respect to the conservative variables
  61. C
  62. USRO = 1.0D0 / RHO
  63. UXSRO = UX * USRO
  64. UYSRO = UY * USRO
  65. GM1 = GAM - 1.0D0
  66. GM1EC = GM1 * 0.5D0 * (UX * UX + UY * UY)
  67. GM1UX = GM1 * UX
  68. GM1UY = GM1 * UY
  69. C
  70. DFRUNC(1) = DFRUNP(1) - UXSRO * DFRUNP(2) - UYSRO * DFRUNP(3)
  71. & + GM1EC * DFRUNP(4)
  72. DFRUNC(2) = USRO * DFRUNP(2) - GM1UX * DFRUNP(4)
  73. DFRUNC(3) = USRO * DFRUNP(3) - GM1UY * DFRUNP(4)
  74. DFRUNC(4) = GM1 * DFRUNP(4)
  75. C
  76. RETURN
  77. END
  78.  
  79.  
  80.  
  81.  
  82.  

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