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 >
Wrap
Text File
|
1991-08-16
|
6KB
|
162 lines
SUBROUTINE HDUCT(ISECT,T,G,D,XOD,ITYPIN,IOLDH,H,IDBG)
C
C SET INTERIOR CONVECTION COEFFICIENT FOR DUCT SECTIONS
C HANDLES LAMINAR, TRANSITION, AND TURBULENT FLOW WITH
C ENTRANCE LENGTH CORRECTION, IF NEEDED.
C
C INPUTS:
C ISECT- DUCT SECTION NO.
C T - MEAN AIR TEMP IN SECTION, F
C G - MASS VELOCITY IN DUCT SECTION, LB/HR-FT2
C D - HYDRAULIC DIA OF DUCT SECTION, FT
C XOD - L/D, WHERE L IS LENGTH FROM START OF BOUNDARY LAYER DEVELOPMENT
C ITYPIN - 0, NO X/L COR; - 1, SQ EDGED INLET; 2, 90 DEG ELBOW AT INLET
C CORRECTIONS CODED BY ITYPIN ARE FOR X/D EFFECTS IN TURB FLOW
C IOLDH =1 FOR ORIG H AS USED IN SINGLE-ZONE MODEL
C IDBG - DEBUG PAR, > 0, TO PRINT DIAGNOSTIC VALUES
C
C OUTPUT:
C H - INTERIOR CONVECTION COEFFICIENT FOR DUCT SECTION, BTU/HR-FT2-F
DIMENSION A1(6),A2(6)
LOGICAL FIRST
C H/H INFINITY MULTIPLIER FOR TURBULENT FLOW
C REF: BOELTER,L.M.K.,ET AL,DISTRIBUTION OF HEAT-TRANSFER RATE IN
C THE ENTRANCE SECTION OF A CIRCULAR TUBE,NACA TN 1451,JULY 1948
C FH1F:FOR DUCT WITH SQ EDGED INLET (DATA WERE PLOTTED BY RDF AND
C WERE ABOUT INDEPENDENT OF REY. USED REY=26700 DATA IN FIG 14)
FH1F(Z)=A1(6)+Z*(A1(5)+Z*(A1(4)+Z*(A1(3)+Z*(A1(2)+Z*A1(1)))))
C FH2F:FOR DUCT WITH A 90 DEG BEND AT INLET (DATA WERE PLOTTED BY RDF
C WERE ABOUT INDEPENDENT OF REY. USED REY=26400 DATA IN FIG 21)
FH2F(Z)=A2(6)+Z*(A2(5)+Z*(A2(4)+Z*(A2(3)+Z*(A2(2)+Z*A2(1)))))
C PROP OF AIR
C VISCOSITY: T-DEG R; UNITS OF MU ARE LB/FT-SEC
F1F(Z)=4.057E-07*Z**1.5/(.5556*Z+110.4)
C THERMAL CONDUCTIVITY: T-DEG R; UNITS OF K ARE BTU/HR-FT-F
F2F(Z)=.004393+1.9556E-5*Z
C C FACTOR IN ASHRAE TURB H EQN: T-DEG F
C REF: ASHRAE HOF 1981, P 2.16, TABLE 7
F3F(Z)=.00368+1.5E-6*(Z-80.)
C PRANDTL NO. OF AIR
DATA PRA/0.72/
C COEF FOR FH1F AND FH2F
DATA A1/-1.386942E-05,7.349947E-04,-1.495856E-02,.1464429,
+ -.7069635,2.482967/
DATA A2/-1.7614E-05,9.046481E-04,-1.782914E-02,.1719179,
+ -.8744688,3.213559/
C * * * * *
DATA FTR/459.7/
C XDPAR=(X/D)/(RE*PRA)
C IF XDPAR < XDPAR1, XNUL=A1N*XDPAR**B1N (LAMINAR NUSSELT NO.)
C REF: FLUID DYNAMICS AND HEAT TRANSFER,KNUDSEN,J.G.,ET AL,MCGRAW-HILL
C 1958, P.375, FIG.13.6 (UNIFORM HEAT FLUX)
DATA XDPAR1,A1N,B1N/.004,.75976,-.4905/
DATA XDPAR2,XNU2/.25,4.36/
C IF XDPAR1 < XDPAR < XDPAR2, XNUL=FITTED FUNCTION OF XDPAR
C RDF FITTED DATA IN FIG.13.6 BETWEEN XDPAR = .004 AND XDPAR=.25
DATA A2N,B2N,C2N/1.64766,.2224241,.06676945/
C FLOW IS FULLY LAMINAR BELOW REYRF1, FULLY TURBULENT ABOVE REYRF2
DATA REYRF1,REYRF2/2100.,7000./
C NO ENTRANCE LENGTH CORRECTION BEYOND XOD = XODR1
DATA XODR1/17./
DATA FIRST/.TRUE./
IF(FIRST)THEN
DREYRF=REYRF2-REYRF1
FIRST=.FALSE.
ENDIF
TA=T+FTR
XMU=F1F(TA)
C=F3F(T)
XK=F2F(TA)
IF(IDBG.GT.0)WRITE(60,500)T,C,XMU,XK,PRA
C
REY=G*D/(3600.*XMU)
IF(REY.LT.REYRF1)IREG=1
IF(REYRF1.LE.REY.AND.REY.LE.REYRF2)IREG=2
IF(REY.GT.REYRF2)IREG=3
IF(IOLDH.EQ.1)IREG=3
IF(IREG.EQ.2)GO TO 50
IF(IREG.EQ.3)GO TO 70
C LAMINAR FLOW
REYQ=REY
40 XDPAR=XOD/(REYQ*PRA)
IF(XDPAR.LT.XDPAR1)THEN
C XNUL IS NUSSELT NUMBER FOR LAMINAR FLOW
XNUL=A1N*XDPAR**B1N
ELSE IF(XDPAR.GT.XDPAR2)THEN
XNUL=XNU2
ELSE IF(XDPAR1.LE.XDPAR.AND.XDPAR.LE.XDPAR2)THEN
XLN=ALOG(XDPAR)
XNUL=EXP(A2N+XLN*(B2N+XLN*C2N))
ENDIF
H=XNUL*XK/D
IF(IREG.EQ.1)GO TO 80
HLAM=H
GO TO 60
50 CONTINUE
C TRANSITION REGION BETWEEN LAMINAR AND TURB FLOW
REYQ=REYRF1
C GET LAMINAR H AT REY= REYRF1 (REY=2100)
GO TO 40
60 GR2=3600.*REYRF2*XMU/D
HTRBR2=C*GR2**.8/D**0.2
IF(XOD.LT.XODR1)THEN
IF(ITYPIN.EQ.1)THEN
HRAT=FH1F(XOD)
ELSE IF(ITYPIN.EQ.2)THEN
HRAT=FH2F(XOD)
ELSE IF(ITYPIN.EQ.0)THEN
HRAT=1.0
ELSE
WRITE(60,506)ITYPIN,ISECT
STOP 'HDUCT:ITYPIN BAD'
ENDIF
ELSE
HRAT=1.0
ENDIF
HR=HRAT*HTRBR2
RAT=(REY-REYRF1)/DREYRF
H=HLAM+RAT*(HR-HLAM)
GO TO 80
70 CONTINUE
C TURBULENT FLOW
G=3600.*REY*XMU/D
HTURB=C*G**0.8/D**0.2
IF(IOLDH.EQ.1)THEN
H=HTURB
GO TO 80
ENDIF
IF(XOD.GT.XODR1.OR.ITYPIN.EQ.0)THEN
H=HTURB
ELSE
C ENTRANCE REGION (1 - SQ EDGE AT ENTRANCE, 2 - ELBOW AT ENTRANCE)
IF(ITYPIN.EQ.1)THEN
HRAT=FH1F(XOD)
H=HRAT*HTURB
ELSE IF(ITYPIN.EQ.2)THEN
HRAT=FH2F(XOD)
H=HRAT*HTURB
ELSE
WRITE(60,506)ITYPIN,ISECT
STOP 'HDUCT:ITYPIN BAD'
ENDIF
ENDIF
80 CONTINUE
IF(IDBG.EQ.0)RETURN
IF(IREG.EQ.1)WRITE(60,501)REY,XOD,XDPAR,XNUL,H
IF(IREG.EQ.2)WRITE(60,502)REY,XOD,XDPAR,HLAM,HTRBR2,HRAT,HR,H
IF(IREG.EQ.3)WRITE(60,505)REY,XOD,HRAT,HTURB,H
RETURN
C
500 FORMAT(1X,'T,C,XMU,XK,PR= ',5G13.5)
501 FORMAT(3X,'REY,XOD,XDPAR,XNUL,H= ',F7.0,1X,F6.3,G13.5,2F7.3)
502 FORMAT(3X,'REY,XOD,XDPAR,HLAM,HTRBR2= ',F7.0,1X,F6.3,G13.5,2F7.3/
+ 3X,'HRAT,HR,H= ',3G13.5)
505 FORMAT(5X,'REY,XOD,HRAT,HTURB,H= ',F7.0,1X,F6.3,G13.5,2F7.3)
506 FORMAT(1H0,'IN HDUCT, ITYPIN IS BAD(= ',I4,' AT ISECT= ',I5)
END
@¼M@│S@╓ @J@J;Ac<AjCAà≡1 ╠Hp
IÜ
Pa
Qü