Esope source files

   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
  86 :       CALL MPFUNT(MACH,MNUM,DERMNU,PNUM,DERPNU)
  87 : C
  88 :       FRHO = MNUM * RHO * CSON
  89 : C
  90 :       DFRO(1) = ((DERMNU * UN) + (MNUM * CSON))  / 2.0D0
  91 :       DFRO(2) = DERMNU * RHO * CNX
  92 :       DFRO(3) = DERMNU * RHO * CNY
  93 :       DFRO(4) = DERMNU * RHO * CNZ
  94 :       DFRO(5) = (MNUM - DERMNU * MACH) * RHO * CSON / (2.0D0 * P)
  95 : C
  96 :       DFRUN(1) = (DFRO(1) * UN) + (DERPNU * MACH * P / (2.0D0 * RHO))
  97 :       CELL = (MNUM * RHO * CSON) + (DERPNU  * P / CSON)
  98 :       DFRUN(2) = (DFRO(2) * UN) + (CELL * CNX)
  99 :       DFRUN(3) = (DFRO(3) * UN) + (CELL * CNY)
 100 :       DFRUN(4) = (DFRO(4) * UN) + (CELL * CNZ)
 101 :       DFRUN(5) = (DFRO(5) * UN) - (0.5D0 * MACH * DERPNU) + PNUM
 102 : C     
 103 :       DFRUT1(1) = DFRO(1) * UT1
 104 :       DFRUT1(2) = (DFRO(2) * UT1)  + (FRHO * CT1X)
 105 :       DFRUT1(3) = (DFRO(3) * UT1)  + (FRHO * CT1Y)
 106 :       DFRUT1(4) = (DFRO(4) * UT1)  + (FRHO * CT1Z)
 107 :       DFRUT1(5) = DFRO(5) * UT1
 108 : C     
 109 :       DFRUT2(1) = DFRO(1) * UT2
 110 :       DFRUT2(2) = (DFRO(2) * UT2)  + (FRHO * CT2X)
 111 :       DFRUT2(3) = (DFRO(3) * UT2)  + (FRHO * CT2Y)
 112 :       DFRUT2(4) = (DFRO(4) * UT2)  + (FRHO * CT2Z)
 113 :       DFRUT2(5) = DFRO(5) * UT2
 114 : C     
 115 :       CELL = CSON2 / (GAM - 1.0D0)
 116 :       DFRET(1) = (DFRO(1) * HT) - (FRHO * CELL / RHO)
 117 :       DFRET(2) = (DFRO(2) * HT)  + (FRHO * UX)
 118 :       DFRET(3) = (DFRO(3) * HT)  + (FRHO * UY)
 119 :       DFRET(4) = (DFRO(4) * HT)  + (FRHO * UZ)
 120 :       DFRET(5) = (DFRO(5) * HT) + (FRHO * CELL / P)
 121 : C
 122 :       RETURN
 123 :       END
 124 :       
 125 :  
 126 :  
 127 :