home *** CD-ROM | disk | FTP | other *** search
/ The World of Computer Software / World_Of_Computer_Software-02-387-Vol-3of3.iso / h / house_ii.zip / FOR / HDUCT.FOR < prev    next >
Text File  |  1991-08-16  |  6KB  |  162 lines

  1.       SUBROUTINE HDUCT(ISECT,T,G,D,XOD,ITYPIN,IOLDH,H,IDBG) 
  2. C  SET INTERIOR CONVECTION COEFFICIENT FOR DUCT SECTIONS
  3. C    HANDLES LAMINAR, TRANSITION, AND TURBULENT FLOW WITH 
  4. C    ENTRANCE LENGTH CORRECTION, IF NEEDED. 
  5. C  INPUTS:  
  6. C   ISECT- DUCT SECTION NO. 
  7. C   T - MEAN AIR TEMP IN SECTION, F 
  8. C   G - MASS VELOCITY IN DUCT SECTION, LB/HR-FT2
  9. C   D - HYDRAULIC DIA OF DUCT SECTION, FT 
  10. C   XOD - L/D, WHERE L IS LENGTH FROM START OF BOUNDARY LAYER DEVELOPMENT 
  11. C   ITYPIN - 0, NO X/L COR; - 1, SQ EDGED INLET; 2, 90 DEG ELBOW AT INLET 
  12. C     CORRECTIONS CODED BY ITYPIN ARE FOR X/D EFFECTS IN TURB FLOW
  13. C  IOLDH =1 FOR ORIG H AS USED IN SINGLE-ZONE MODEL 
  14. C   IDBG - DEBUG PAR, > 0, TO PRINT DIAGNOSTIC VALUES 
  15. C  OUTPUT:  
  16. C   H - INTERIOR CONVECTION COEFFICIENT FOR DUCT SECTION, BTU/HR-FT2-F
  17.       DIMENSION A1(6),A2(6) 
  18.       LOGICAL FIRST 
  19. C  H/H INFINITY MULTIPLIER FOR TURBULENT FLOW 
  20. C  REF: BOELTER,L.M.K.,ET AL,DISTRIBUTION OF HEAT-TRANSFER RATE IN
  21. C       THE ENTRANCE SECTION OF A CIRCULAR TUBE,NACA TN 1451,JULY 1948
  22. C  FH1F:FOR DUCT WITH SQ EDGED INLET (DATA WERE PLOTTED BY RDF AND
  23. C  WERE ABOUT INDEPENDENT OF REY.  USED REY=26700 DATA IN FIG 14) 
  24.       FH1F(Z)=A1(6)+Z*(A1(5)+Z*(A1(4)+Z*(A1(3)+Z*(A1(2)+Z*A1(1))))) 
  25. C  FH2F:FOR DUCT WITH A 90 DEG BEND AT INLET (DATA WERE PLOTTED BY RDF
  26. C  WERE ABOUT INDEPENDENT OF REY.  USED REY=26400 DATA IN FIG 21) 
  27.       FH2F(Z)=A2(6)+Z*(A2(5)+Z*(A2(4)+Z*(A2(3)+Z*(A2(2)+Z*A2(1))))) 
  28. C  PROP OF AIR
  29. C   VISCOSITY: T-DEG R; UNITS OF MU ARE LB/FT-SEC 
  30.       F1F(Z)=4.057E-07*Z**1.5/(.5556*Z+110.4) 
  31. C   THERMAL CONDUCTIVITY: T-DEG R; UNITS OF K ARE  BTU/HR-FT-F
  32.       F2F(Z)=.004393+1.9556E-5*Z
  33. C   C FACTOR IN ASHRAE TURB H EQN: T-DEG F
  34. C   REF: ASHRAE HOF 1981, P 2.16, TABLE 7 
  35.       F3F(Z)=.00368+1.5E-6*(Z-80.)
  36. C  PRANDTL NO. OF AIR 
  37.       DATA PRA/0.72/
  38. C  COEF FOR FH1F AND FH2F 
  39.       DATA A1/-1.386942E-05,7.349947E-04,-1.495856E-02,.1464429,
  40.      +        -.7069635,2.482967/ 
  41.       DATA A2/-1.7614E-05,9.046481E-04,-1.782914E-02,.1719179,
  42.      +        -.8744688,3.213559/ 
  43. C  * * * * *
  44.       DATA FTR/459.7/ 
  45. C  XDPAR=(X/D)/(RE*PRA) 
  46. C   IF XDPAR < XDPAR1, XNUL=A1N*XDPAR**B1N  (LAMINAR NUSSELT NO.) 
  47. C   REF: FLUID DYNAMICS AND HEAT TRANSFER,KNUDSEN,J.G.,ET AL,MCGRAW-HILL
  48. C        1958, P.375, FIG.13.6 (UNIFORM HEAT FLUX)
  49.       DATA XDPAR1,A1N,B1N/.004,.75976,-.4905/ 
  50.       DATA XDPAR2,XNU2/.25,4.36/
  51. C   IF XDPAR1 < XDPAR < XDPAR2, XNUL=FITTED FUNCTION OF XDPAR 
  52. C   RDF FITTED DATA IN FIG.13.6 BETWEEN XDPAR = .004 AND XDPAR=.25
  53.       DATA A2N,B2N,C2N/1.64766,.2224241,.06676945/
  54. C  FLOW IS FULLY LAMINAR BELOW REYRF1, FULLY TURBULENT ABOVE REYRF2 
  55.       DATA REYRF1,REYRF2/2100.,7000./ 
  56. C  NO ENTRANCE LENGTH CORRECTION BEYOND XOD = XODR1 
  57.       DATA XODR1/17./ 
  58.       DATA FIRST/.TRUE./
  59.       IF(FIRST)THEN 
  60.          DREYRF=REYRF2-REYRF1 
  61.          FIRST=.FALSE.
  62.          ENDIF
  63.       TA=T+FTR
  64.       XMU=F1F(TA) 
  65.       C=F3F(T)
  66.       XK=F2F(TA)
  67.       IF(IDBG.GT.0)WRITE(60,500)T,C,XMU,XK,PRA 
  68.       REY=G*D/(3600.*XMU) 
  69.       IF(REY.LT.REYRF1)IREG=1 
  70.       IF(REYRF1.LE.REY.AND.REY.LE.REYRF2)IREG=2 
  71.       IF(REY.GT.REYRF2)IREG=3 
  72.       IF(IOLDH.EQ.1)IREG=3
  73.       IF(IREG.EQ.2)GO TO 50 
  74.       IF(IREG.EQ.3)GO TO 70 
  75. C  LAMINAR FLOW 
  76.       REYQ=REY
  77.    40 XDPAR=XOD/(REYQ*PRA)
  78.       IF(XDPAR.LT.XDPAR1)THEN 
  79. C  XNUL IS NUSSELT NUMBER FOR LAMINAR FLOW
  80.          XNUL=A1N*XDPAR**B1N
  81.       ELSE IF(XDPAR.GT.XDPAR2)THEN
  82.          XNUL=XNU2
  83.       ELSE IF(XDPAR1.LE.XDPAR.AND.XDPAR.LE.XDPAR2)THEN
  84.          XLN=ALOG(XDPAR)
  85.          XNUL=EXP(A2N+XLN*(B2N+XLN*C2N))
  86.          ENDIF
  87.       H=XNUL*XK/D 
  88.       IF(IREG.EQ.1)GO TO 80 
  89.       HLAM=H
  90.       GO TO 60
  91.    50 CONTINUE
  92. C  TRANSITION REGION BETWEEN LAMINAR AND TURB FLOW
  93.       REYQ=REYRF1 
  94. C       GET LAMINAR H AT REY= REYRF1  (REY=2100)
  95.       GO TO 40
  96.    60 GR2=3600.*REYRF2*XMU/D
  97.       HTRBR2=C*GR2**.8/D**0.2 
  98.       IF(XOD.LT.XODR1)THEN
  99.          IF(ITYPIN.EQ.1)THEN
  100.             HRAT=FH1F(XOD)
  101.          ELSE IF(ITYPIN.EQ.2)THEN 
  102.             HRAT=FH2F(XOD)
  103.          ELSE IF(ITYPIN.EQ.0)THEN 
  104.             HRAT=1.0
  105.          ELSE 
  106.             WRITE(60,506)ITYPIN,ISECT
  107.             STOP 'HDUCT:ITYPIN BAD' 
  108.             ENDIF 
  109.       ELSE
  110.          HRAT=1.0 
  111.          ENDIF
  112.       HR=HRAT*HTRBR2
  113.       RAT=(REY-REYRF1)/DREYRF 
  114.       H=HLAM+RAT*(HR-HLAM)
  115.       GO TO 80
  116.    70 CONTINUE
  117. C  TURBULENT FLOW 
  118.       G=3600.*REY*XMU/D 
  119.       HTURB=C*G**0.8/D**0.2 
  120.       IF(IOLDH.EQ.1)THEN
  121.          H=HTURB
  122.          GO TO 80 
  123.          ENDIF
  124.       IF(XOD.GT.XODR1.OR.ITYPIN.EQ.0)THEN 
  125.          H=HTURB
  126.       ELSE
  127. C  ENTRANCE REGION (1 - SQ EDGE AT ENTRANCE, 2 - ELBOW AT ENTRANCE) 
  128.          IF(ITYPIN.EQ.1)THEN
  129.             HRAT=FH1F(XOD)
  130.             H=HRAT*HTURB
  131.          ELSE IF(ITYPIN.EQ.2)THEN 
  132.             HRAT=FH2F(XOD)
  133.             H=HRAT*HTURB
  134.          ELSE 
  135.            WRITE(60,506)ITYPIN,ISECT 
  136.            STOP 'HDUCT:ITYPIN BAD'
  137.            ENDIF
  138.          ENDIF
  139.    80 CONTINUE
  140.       IF(IDBG.EQ.0)RETURN 
  141.       IF(IREG.EQ.1)WRITE(60,501)REY,XOD,XDPAR,XNUL,H 
  142.       IF(IREG.EQ.2)WRITE(60,502)REY,XOD,XDPAR,HLAM,HTRBR2,HRAT,HR,H
  143.       IF(IREG.EQ.3)WRITE(60,505)REY,XOD,HRAT,HTURB,H 
  144.       RETURN
  145.   500 FORMAT(1X,'T,C,XMU,XK,PR= ',5G13.5) 
  146.   501 FORMAT(3X,'REY,XOD,XDPAR,XNUL,H= ',F7.0,1X,F6.3,G13.5,2F7.3)
  147.   502 FORMAT(3X,'REY,XOD,XDPAR,HLAM,HTRBR2= ',F7.0,1X,F6.3,G13.5,2F7.3/ 
  148.      +       3X,'HRAT,HR,H= ',3G13.5) 
  149.   505 FORMAT(5X,'REY,XOD,HRAT,HTURB,H= ',F7.0,1X,F6.3,G13.5,2F7.3)
  150.   506 FORMAT(1H0,'IN HDUCT, ITYPIN IS BAD(= ',I4,'  AT ISECT= ',I5) 
  151.       END 
  152. @¼M@│S@╓ @J @J ;Ac <Aj CAà ≡1 ╠Hp
  153. IÜ
  154. Pa
  155. Qü
  156. B    j Yö î eGæHmIÖK{MéOòPtQ├RÑ SUsëtìuGvGw8ä8iö %a&]/É R₧ ╠A   |I P3CRT     BAT Xç╧QDIRM    BAT ë|yRED      BAT tÇFSCOMM    BAT FM*σALRESETEXE {_ADN      BAT ╜}╬U σALSETUPCOM wb╙1DR      BAT HsVDIRP    BAT πK-W