Sources Esope | Cast3M

bt2
C BT2       SOURCE    CB215821  17/11/30    21:15:16     9639           
      FUNCTION BT2(T1,TP,S,BH,PR,RE,XLA,XL,DX)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
C
C     operateur FUITE
C     equation resolue
C      DT/dx= 2BH/(Q Cp)    (Tp-T) + S/dx
C           = 4BH/(RePrXla) (Tp-T) + S/dx
C
C>>>  calcule la temperature de sortie d'un tron\on en fonction de :
C       T1 : temperature d'entree
C       TP : temperature de paroi
C       S  : terme source (homogene a une temp{rature)
C       BH : coefficient d'echange
C       H  : coefficient d'echange sans dimension
C       PR : nombre de Prandtl
C       RE : nombre de Reynolds
C       XLA: conductivite de l'air
C       XL : longueur fissure
C       DX : pas d'espace local
C
      H=BH*2*2*DX*XL/RE/PR/XLA
 
C     HMIN=20.
C     H1=MIN(HMIN,H)
      H1=H
C     write(6,*) ' bt2 H1 ',H1
C     BT2=TP+((T1-TP)*EXP(-H1))+( S*(1-EXP(-H1)) /H1)
 
      EH1 = EXP(-H1)
      DT1 = (T1-TP)*EH1
 
      IF (ABS(H1).LT.1D-3) THEN
        DT2 = S
      ELSE
        DT2 = (S*(1-EH1)/H1)
      ENDIF
 
      BT2= TP + (DT1 + DT2)
 
C     write(6,*) ' bt2 tp dt1 dt2 ', BT2, TP, DT1,DT2
      RETURN
      END