Esope source files

   1 : C BT3       SOURCE    CHAT      05/01/12    21:44:21     5004           
   2 :       FUNCTION BT3(T1,TP,S,H,Q,DX,XL)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8 (A-H,O-Z)
   5 : C
   6 : C      operateur FISS, option gaz reel
   7 : C      equation resolue
   8 : C       Q. DT/dx= H. (Tp-T) + S
   9 : C>>>  calcule la temperature de sortie d'un troncon en fonction de :
  10 : C       T1 : temperature d'entree
  11 : C       TP : temperature de paroi
  12 : C       S  : terme source 
  13 : C       H : coefficient d'echange
  14 : C       Q  : debit (homogene a  q.cp = ro.u.e.cp)
  15 : C       DX : pas d'espace local sans dimension
  16 : C       XL : longueur de la fissure
  17 :       
  18 :       H1=(H/Q)*DX*XL
  19 : 
  20 : C     BT3=TP+((T1-TP)*EXP(-H1))+((S/H)*(1-EXP(-H1)))
  21 : C     write(6,*) ' bt3 H1 T2 ',H1,BT3
  22 : 
  23 :       EH1 = EXP(-H1)
  24 :       DT1 = (T1-TP)*EH1
  25 : C
  26 :       IF (ABS(H1).LT.1D-5) THEN
  27 :         DT2 = (S/Q)*(DX*XL)
  28 :       ELSE
  29 :         DT2 = (S/H)*(1-EH1)
  30 :       ENDIF
  31 : 
  32 :       BT3= TP + (DT1 + DT2)
  33 : 
  34 : C     write(6,*) ' bt2 tp dt1 dt2 ', BT2, TP, DT1,DT2
  35 :       RETURN
  36 :       END
  37 :  
  38 :  
  39 :  
  40 :