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