Sources Esope | Cast3M

bt3
C BT3       SOURCE    CHAT      05/01/12    21:44:21     5004
      FUNCTION BT3(T1,TP,S,H,Q,DX,XL)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
C
C      operateur FISS, option gaz reel
C      equation resolue
C       Q. DT/dx= H. (Tp-T) + S
C>>>  calcule la temperature de sortie d'un troncon en fonction de :
C       T1 : temperature d'entree
C       TP : temperature de paroi
C       S  : terme source
C       H : coefficient d'echange
C       Q  : debit (homogene a  q.cp = ro.u.e.cp)
C       DX : pas d'espace local sans dimension
C       XL : longueur de la fissure
 
      H1=(H/Q)*DX*XL
 
C     BT3=TP+((T1-TP)*EXP(-H1))+((S/H)*(1-EXP(-H1)))
C     write(6,*) ' bt3 H1 T2 ',H1,BT3
 
      EH1 = EXP(-H1)
      DT1 = (T1-TP)*EH1
C
      IF (ABS(H1).LT.1D-5) THEN
        DT2 = (S/Q)*(DX*XL)
      ELSE
        DT2 = (S/H)*(1-EH1)
      ENDIF
 
      BT3= TP + (DT1 + DT2)
 
C     write(6,*) ' bt2 tp dt1 dt2 ', BT2, TP, DT1,DT2
      RETURN
      END