Télécharger bt2.eso

Retour à la liste

Numérotation des lignes :

  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.  
  48.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales