Esope source files

   1 : C BT2       SOURCE    CHAT      05/01/12    21:44:18     5004           
   2 :       FUNCTION BT2(T1,TP,S,BH,PR,RE,XLA,XL,DX)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8 (A-H,O-Z)
   5 : C
   6 : C     operateur FUITE
   7 : C     equation resolue
   8 : C      DT/dx= 2BH/(Q Cp)    (Tp-T) + S/dx
   9 : C           = 4BH/(RePrXla) (Tp-T) + S/dx
  10 : C
  11 : C>>>  calcule la temperature de sortie d'un tron\on en fonction de :
  12 : C       T1 : temperature d'entree
  13 : C       TP : temperature de paroi
  14 : C       S  : terme source (homogene a une temp{rature)
  15 : C       BH : coefficient d'echange
  16 : C       H  : coefficient d'echange sans dimension
  17 : C       PR : nombre de Prandtl
  18 : C       RE : nombre de Reynolds
  19 : C       XLA: conductivite de l'air
  20 : C       XL : longueur fissure
  21 : C       DX : pas d'espace local
  22 : C
  23 :       H=BH*2*2*DX*XL/RE/PR/XLA
  24 : 
  25 : C     HMIN=20.
  26 : C     H1=DMIN1(HMIN,H)
  27 :       H1=H
  28 : C     write(6,*) ' bt2 H1 ',H1
  29 : C     BT2=TP+((T1-TP)*EXP(-H1))+( S*(1-EXP(-H1)) /H1)
  30 : 
  31 :       EH1 = EXP(-H1)
  32 :       DT1 = (T1-TP)*EH1
  33 : 
  34 :       IF (ABS(H1).LT.1D-3) THEN
  35 :         DT2 = S
  36 :       ELSE
  37 :         DT2 = (S*(1-EH1)/H1)
  38 :       ENDIF
  39 : 
  40 :       BT2= TP + (DT1 + DT2)
  41 : 
  42 : C     write(6,*) ' bt2 tp dt1 dt2 ', BT2, TP, DT1,DT2
  43 :       RETURN
  44 :       END
  45 :  
  46 :  
  47 :