Copyright Warning & Restrictions The copyright law of the United States (Title 17, United States Code) governs the making of photocopies or other reproductions of copyrighted material. Under certain conditions specified in the law, libraries and archives are authorized to furnish a photocopy or other reproduction. One of these specified conditions is that the photocopy or reproduction is not to be “used for any purpose other than private study, scholarship, or research.” If a, user makes a request for, or later uses, a photocopy or reproduction for purposes in excess of “fair use” that user may be liable for copyright infringement, This institution reserves the right to refuse to accept a copying order if, in its judgment, fulfillment of the order would involve violation of copyright law. Please Note: The author retains the copyright while the New Jersey Institute of Technology reserves the right to distribute this thesis or dissertation Printing note: If you do not wish to print this page, then select “Pages from: first page # to: last page #” on the print dialog screen
53
Embed
Cprht Wrnn & trtn - archives.njit.eduarchives.njit.edu/vol01/etd/1950s/1959/njit-etd1959-002/njit-etd1959... · nu = f of (ra)o25 nu = f of (ra)e33 the rayleigh - jeffries oriterion
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Copyright Warning & Restrictions
The copyright law of the United States (Title 17, UnitedStates Code) governs the making of photocopies or other
reproductions of copyrighted material.
Under certain conditions specified in the law, libraries andarchives are authorized to furnish a photocopy or other
reproduction. One of these specified conditions is that thephotocopy or reproduction is not to be “used for any
purpose other than private study, scholarship, or research.”If a, user makes a request for, or later uses, a photocopy orreproduction for purposes in excess of “fair use” that user
may be liable for copyright infringement,
This institution reserves the right to refuse to accept acopying order if, in its judgment, fulfillment of the order
would involve violation of copyright law.
Please Note: The author retains the copyright while theNew Jersey Institute of Technology reserves the right to
distribute this thesis or dissertation
Printing note: If you do not wish to print this page, then select“Pages from: first page # to: last page #” on the print dialog screen
The Van Houten library has removed some ofthe personal information and all signatures fromthe approval page and biographical sketches oftheses and dissertations in order to protect theidentity of NJIT graduates and faculty.
HEAT TRANSFER BY NATURAL CONVECTION
IN HORIZONTAL FLUID LAYERS
A STATISTICAL CORRELATION OF EXPERIMENTAL DATA
BY JOHN L. O' TooLE
A THESIS SUBMITTED TO THE FACULTY OF THE DEPARTMENT OFCHEMICAL ENGINEERING OF NEWARK COLLEGE OF ENGINEERING
IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THEDEGREE OF MASTER OF SCIENCE IN CHEMICAL ENGINEERING.
APPROVAL OF THESIS
FOR DEPARTMENT OF CHEMICAL ENGINEERINGNEWARK COLLEGE OF ENGINEERING
FACULTY COMMITTEE
APPROVED
NEWARK, N. J.JUNE, 1959
Lib raryCollege of Engineering
THE RATE OF HEAT TRANSFER DUE TO NATURAL CONVECTION UN CONFINED,
HORIZONTAL 0 FLUID LAYERS CAN BE EXPRESSED BY EQUATIONS INVOLVING THE
DIMENSIONLESS PARAMETERS NUSSELT NO,, RAYLEIGH NO., AND PRANDTL NO.,
WHERE RAYLEIGH NO. EQUALS THE PRODUCT OF GRASHOF AND PRANDTL NOS.
STATISTICAL CORRELATION BY STEPWIS REGRESSION OF
EXPERIMENTAL DATA OBTAINED FROM FOUR DIFFERENT INVESTIGATORS, ON BOTH
GASES AND LIQUIDS, YIELDS THREE SEPARATE EQUATIONS CORRESPONDING TO
THREE DISTINCT REGIMES OF CONVECTION. THE TOTAL RANGE SPANNED BY THE
DATA IS FROM THE ONSET OF CONVECTION TO RAYLEIGH = 10 9 .
THE EQUATIONS ARE
RAREGIME RANGE EQUATION
.816INITIAL 1600..3000 Nu = 0.00238(RA)
LAMINAR 3000..105 Nu = 0.221(RA).256
TURBULENT 105°4109 Nu = 0.0891 (RA)·316 (PR)·0853
THE STANDARD ERROR OF ESTIMATE, WHICH IS THE PROBABLE ERROR IN THE
PREDICTION OF NU FROM A REGRESSION EQUATION, IS OBTAINED FOR EACH
REGIME. A VERY HIGH DEGREE OF ACCURACY WAS OBTAINED FOR ALL. CORRELATIONS&
STD, ERROROF ESTI MATE
TABLE OF CONTENTS
SUBJECT PAGE
INTRODUCTION 1
THEORY 3
EXPERIMENTAL DATA 8
THE CORRELATION
TABLE 1 *I SOURCES OF THE DATA 8
TABLE if As FEATURES OF THE APPARATUS USED BY
SOURCES OF THE DATA 17
TABLE III · EXPERIMENTAL DATA I8.26
TABLE IV P.4 SUMMARY OF REGRESSION ANALYSIS 13
TABLE V · RESULTS OF COMPUTER ANALYSIS 27
TABLE VI COEFFICIENTS OF CORRELATION AND
DETERMINATION 28
FIG. 1 PRANDTL NO. OF AIR VS TEMPERATURE BETWEEN 28 & 29
FIG, 2 THERMAL CONDUCTIVITY OF AIR VS TEMPERATURE'
FIG. 3— NUSSELT No, VS RAYLEIGH No.: EXPERIMENTAL
DATA
FIG. 4 NUSSELT NO. VS RAYLEIGH No.: REGRESSION
EQUATIONS
Flo. 5 - NUSSELT NO./(PRANDTL NO.)• 0853 VS
RAYLEIGH No.: TURBULENT REGIME
APPENDIX 1 DIMENSIONAL ANALYSES 29
APPENDIX Il · DERIVATION OF THE COEFFICIENT OF
THERMAL EXPANSION OF AIR 31
TABLE OF CONTENTS (CO NT
SUBJECT PAGE
APPEND IX III - CORRECT ION OF THE DATA OF
MULL AND RE (HER 33
APPEND IX 1V - STEPWISE LINEAR MULTIPLE REGRESS ION 35
TABLE OF SYMBOLS 39
BIBLIOGRAPHY AND REFERENCES 40
HEAT TRANSFER BY NATURAL CONVECTI~N.
iN HORIZONTAL FLUID LAYERS
I NTRODUCT 1 ON
HEAT TRANSFER BY NATURAL OONVECTION IN CONFINED, HOR.ZONTAL~
FLUrO LAYERS WAS FURST INVESTIGATED ~ORMALLY BY BENARD (1)* IN
1901 AND ANALYZED BY LORD RAYLEIGH [N 1916 (2). SUNCE THAT TiME A
CONSIDERABLE AMOUNT OF WORK HAS APPEARED IN THE LITERATURE OONSISTING
LARGELY OF MATHEMATICAL ANALYSES DEFINING THE OONDITIONS UNCER WHIOH
A HEATED LAYER FIRST BEOOMES UNSTABLE AND BEGINS TO OIROULATE. PRIOR
TO 1950 ONLY A SMALL AMOUNT OF EXPERIMENTAL WORK HAD BEEN DONE. SfNOE
THEN, HOWEVER, THREE PAPERS HAVE BEEN PUBLISHED (12, 20$ 21) REPORTING
EXPERIMENTAL HEAT TRANSFER DATA OVER A WIDE RANGE OF THE PHENOMENON
AFTER THE ONSET OF CONVEOTION. WORKING INDEPENDENTLY, WITH DIFFERENT
FLUIDS, AND IN DIFFERENT RANGES, THE AUTHORS DIFFERED IN THEIR
OORRELATIONS AND NO COMPREHENSIVE RELATIONSHIp WAS DEVELOPED.
THIS THESIS PRESENTS AN ENGINEERING CORRELATION BY STEPWISE
MUll'lPLE REGRESSION OF THE AVAiLABLE DATA. THIS METHOD IS A POWERFUL
STATISTICAL TOOL THAT PROVIDES A CORRELATING EQUATION FOR A MUlTIPllOITY
OF PARAMETERS AND, AT THE SAME TIME, YIELDS AN ESTIMATE OF THE AOCURACY
OF THE CORRELATiON AND THE SIGNIFIOANOE OF THE INDIVIDUAL PARAMETERS.
To THE AUTHORtS KNOWLEDGE, THIS THESIS IS THE FIRST USE OF STEPwiSE
MULTIPLE REGRESSION iN THE ANALYSIS AND CORRELATION OF HEAT TRANSFER
DATA.
THREE OF THE FOUR SETS OF DATA (12, 20, 21) WERE OBTAINED
DIRECTLY FRCM THE AUTHORS ~N THE FORM OF THE D_MENSIONLESS
PARAMETERS NUBSELT. RAYLEIGH, GRASHOF, AND PRANGTL NUMBERS. (THESE
WILL BE ABBREVIATED Na, RA, GR, PR HEREAFTER.) THE OTHER SET (5)
WAS CALOULATED FROM RAW DATA TABULATED IN THE PAPER$
* NUMBERS IN PARENTHESES REFER TO REFERENOES LIBTED IN THE BIBLIOGRAPHY.
THEORY
THE SYSTEM
WHEN A TEMPERATURE DIFFERENOE Is APPLIED AOROSS A HORIZONTAL
LAYER OF FLUID THAT as INF(NfTE iN EXTENT AND OONFINED BETWEEN
OONDUOTING HORIZONTAL SURFAOES, A TEMPERATURE GRADIENT WILL BE ESTABLISHED
ACROSS THE LAYER AND HEAT wtLL fLOW BY OONDUCTION (AND, IN THE OASE
1, THE HOT SURFAOE is BENEATH THE LAYER, THE FLUID WiLL TEND TO (f/l4t,V,.l(
01 ROULATE ~ to A DENS ITY GRAD lENT AOROSS THE LAYER B NDUOED BY THE
TEMPERATURE GRADIENT. THIS OIROULATION M CONVEOTION RESUL T5 I N AN
- I -
INOREASE IN THE RATE OF HEAT TRANSFER OVER THAT DUE TO OONDUCTION ALONE.
TH8S THESIS OONSIDERS THE MAGNITUDE OF THE HEAT TRANSFER DUE TO
OONVEOTION.
ONSET OF CONVEOTiON
LORD RAYLEiGH (2) FOUND BY ANALYSIS THAT NO OONVEOTION WOULD TAKE
PLAOE UNTIL A ORaTIOAL VALUE OF THE RAYLEIGH NUMBER (RA = GR X PRJ WAS
EXOEEDEDo HIS ANALYSIS HAS BEEN OONrlRMED AND EXTENDED BY OTHERS
(3,4,9, II) AND THE ORITIOAL RA OF 1707.8 OALOULATED BY JEFFRIES (3)
HAS BEE N V E R I FIE D EX PER t MEN TAL L V ( 6 _ 7, to, 1 2, I 4 II 20) I>
HEAT TRANSFER - DrMENSIONAL ANALystS
THE RAYLEIGH - JEFFRlES ANALVSIS DOES NOT YIELD A FUNOTiONAL
RELATIONSHIP BETWEEN THE RATE OF HEAT TRANSFER AND THE SYSTEM VARIABLES
AND CONSTANTS AFTER CONVECTION HAS BEGUN. THERE HAVE aEEN VERY FEW
ATTEMPTS iN THE L~TERATURE TO DERIVE SUOH A RELAT(ONSH~P CHIEF~Y
BECAUSE THE DiFFERENTIAL EQUATIONS DERiVED FROM THE BASIC FLOW
EQUATION DO NOT YIELD TO liNEARIZATION (16~ 19). HOWEVER~ BATOHELER (16)
HAS DERIVED THE FOLLOWING EXPRESSIONS FOR THE OASE OF CONVEOTION IN
VERTioAL LAYERS:
lAMiNAR RE(UON
TURBULENT REGION
...
... Nu = F OF (RA)o25
Nu = F OF (RA)e33
THE RAYLEIGH - JEFFRIES ORITERION SUGGESTS THAT THE RATE OF HEAT
TRANSFER IS SOME FUNOTtON OF RA AND SUOH A FUNOTION OAN BE DERIVED BY
DIMENSIONAL ANALYSIS.
IF WE POSTULATE THAT AT STEADY"'STATE CONDITIONS THE OVERALL RATE
OF HEAT TRANSFER PER UNIT AREA H IS A FUNCTION OF
AREA BY CONDUOTION AND
~ = THERMAL CONO_CTTVITY OF THE FLUID
T = TEMPERATURE DIFFERENOE
L = THIOKNESS OF THE FLUID LAYER
e = DENSITY OF THE FLUID
~ = DYNAMIC VISCOSITY OF THE FLUID
o = SPECIFIC HEAT AT CONSTANT PRESSURE OF THE FLUID
, - COEFF,OlENT OF THERMAL EXPANSION OF THE FLU.'
Q = GRAV2TATIONAl CONSTANT
(ALL CONSTANTS ARE EVALUATED AT THE MEAN
TEMPERATURE OF THE LAYER)
THEN THE FOLLOWING ExPRESSION OAN BE DERIVED BY DIMENSIONAL ANALYSIS:(23)
(SEE ApPEND 1')( 1)
F OF (~ G Y L 3 e 2 ) N ( fl'2. )
H Is A FORM OF NUSSELT NUMBER AND fS THE RATfm OF OVERALL KY/L
RATE OF HEAT TRANSFER DUE TO CONVECTiON AND CONDUCTION TO THE RATE OF
HEAT TRANSFER DUE TO OONDUCTION ALONE. WHEN IT is EQUAL TO ONE, NO
HEAT Is BEfNG TRANSFERRED BY OONVECTIONo
~ G T 1.3 f 2 I S THE GRASHOF NUMBER.. t1 CAN BE REGlARDED AS THE 2
R A 1 BOO FeU 0 Y ANT FOR C E S " (~G T)>> TOT HE V I S 00 U s r 0 ROE 15 (r / (' ), 0 F
THE FLUID (22).
THE SISNIF"lCANCE OF THE CRn'()CAL RA DiSCOVERED BY LORD RAYLE1GH ,
THEN, IS THAT THE BUOYANT FORCE DUE TO THE TEMPERATURE GRADIENT IS
BALANCED BY THE vISCOUS FORCES IN THE LAYER UNTIL THE BUOYANT rORot
ATTAtNS A VALUE REPRESENTEP BY THE CRITICAL RA .. WHEN THIS HAPPENS, ('1 ~\A/
FLOW COMMENPES"
....!!2 IS THE PRANDTL NUMBER, WHICH IS THE RATIO OF MOMENTUM D·U-rUSI:vt1Y K
TO THERMAL D1FFUSrVITY IN THE FLUlo;t.E .. , THE RATIO OF THE FLutD
PROPERTY GOVERNING THE TRANSFER OF MOMENTUM BY VISCOUS FOROES TO THE
FLUID PROPERTY GOVERNING TRANSFER OF HEAT BY TEMPERATURE mIFFERENOE.
THE DIMENSIONAL ANALYSIS ASSUMES THAT HEAT TRANSFER THROUGH THE
LAYER is NOT AFFEOTED BY THE PARTIOULAR GEOMETRY OF THE FLOW PATTERNS
OR THE ABSOLUTE MAGNITUDE OF THE TEMPERATURE DIFFERENCE. 1M ADDITION.
IT ASSUMES THAT THE PHYS!OAL OONSTANTS CAN BE EVALUATED AT THE MEAN
TEMPERATURE OF THE LAYERo
MODES OF CONVEOTION . AFTER THE ONSET OF CONVECTION~ DIFFERENT MODES OF FLOW HAVE BEEN
OBSERVED (6, 7, 8, 12~ 20, 21). INITlALLY THE FLOW TAKES PLACE fN
D1s0REET OELLS THAT HAVE CROSS SECTIONS OF REGULAR POLYG~NS. As RA
INCREASES THE MODE CHANGES AS SHOWN BY CHANGES iN THE SHAPE AND WIDTH
OF THE CELLS. FINAllY, THE ORDERLY CELLS DISAPPEAR ALTOGETHER AS THE
OONVECTION CHANGES FROM LAMINAR TO TURBULENT FLOW.
A NUMBER OF ANALYSES HAVE BEEN PUBLISHED PREDIOTING 1HI RA AT
WHICH MODE TRANSITIONS OOOUR (3, 9, 14), AND SOME EXPERIMENTERS HAVE
TRIED TO CORRELATE THEIR HEAT TRANSFER DATA WITH THE CHANGES IN MODE
THAT THEY OBSERVED. HOWEVER, MALKUS (15, 17) REPORTED THAT THE HEAT
TRANSFER FUNOTION OHANGES ONLY WITH THE CHANGE FROM LAMINAR TO TURBULENT
FLOW AND DOES NOT OORRESPOND TO CHANGES !N MODE WITHIN THESE REGlUES.
FORM OF CORRELATION
THIS THESiS PRESENTS AN EUPER1CAL CORRELATION OF Nu VI THREE
D~UENSIONLESS PARAMETERS ~ RA, PRJ AND A/o 2 , WHERE A/02 IS THE RATte
OF AREA -TO THICKNESS SQUARED OF THE APPARATUS USED IN OBTAINING THE DATA.
THESE PARAMETERS WERE SELECTED FOR THE FOLLOWINB REASONS:
I. RA IS PREDIOTED BY DIMENSiONAL ANALYSIS.
2~ PR IS REPORTED (20, 21) TO BE AN ADDITIONAL PARAMETER
WITH RA IN THE TU~BULENT REBION.
3. A/o2 I S DES I GNEO TO TEST THE EXPER IMENTAL DATA F'OR
AGREEMENT WITH THE ASSUMPTION OF INFrNfTE EXTENT OF
THE lAYERa
E:XPEFUMENTAL DATA
SOURCES
THE DATA WERE OBTAINED FROM THE FOUR SOURCES LISTED iN TABLE I~
IN THREE CASES THE DATA WERE OBTAINED DIRECTLY FROM THE AUTHORS. IN
THE OTHER CASE {MULL AND REIHER) RA. DATA WERE TAKEN FROM A TABULATION
IN THE PAPER, AND THE PARAMETERS CAl.OULATED FROM IT.
TABLE
SOURCE AND RANGE OF EXPERIMENTAL DATA
DATE
SOURCE PUBLISHED
DEGRAAF" AND VAN DER HELD(12) 1952
GLOBE AND" DROPK 1 N (21 ) 1958 WATER, ILICONE o U.S 9 MERCURY
MULL ANI!! REI HErd5) 1930 AIR
SCHM I DT AND S H. VESTON (20) 1958 WATER,. GLYCOL~ HEPTANE, SILIOONE 01'1..8
ApPARATUS AND PROCEDURE
RA
RANGE
103 .... 105
r05 .109
10 3 ... 0106
103
... '°5
No. OF"
DATA
26
56
17
106
ALL OF" THE SOURCES OF DATA USED SIMILAR APPARATUS, DIFFERING
PRINCIPALLY iN DETAILS OF CONSTRUOTION AND METHODS OF MEASURING THE
NET RATE OF HEAT TRANSFER. A COMPARISON OF THE iMPORTANT FEATURES OF
* THE1R APPARATUS is MADE iN TABLE II. IN GENERAL g THE FLUID LAYER WAS
CONFINED BETWEEN TWO PARALLEL METAL PLATES AND AN ENCLOSING WALL O~
INSULATING MATERIAL. SPACiNG WAS FiXED EITHER BY SMALL SPACERS BETWEEN
PLATES OR BY THE ENCLOSURE. THE LOWER PLATE WAS EQUiPPED WITH ELEOTRICAL
HEATERS g WHilE THE UPPER PLATE WAS'DESIGNED TO AOOOMMODATE A FLOW OF
COOLING WATER OR AIR. TEMPERATURE DI~FERENCE WAS MEASURED BY
*' Ps. 17
THERMOCOUPLES EMBEDDED BN THE PLATES. HEAT TRANSFER RATE WAS
CALCULATED EiTHER BY MEASURING THE POWER OONSUMPTION OF THE HEATERS
OR THE HEAT ACQUIRED BY THE OOOLING MEDIUM .. THE PROPERTIES OF THE
FLaRe WERE TAKEN AT THE ARiTHMETIC MEAN TEMPERATURE OF THE HOT AND
COLD PLATES. IN THE EXPERIMENTS ON AIR, THE HEAT TRANSFER DUE TO
RADIATION WAS DETERMINED OVER THE TEMPERATURE RANGE OF THE
EXPERiMENTS BY INVERTING THE NOT AND COLO SURFAOES; I.E .. , BY HEATING
FROM THE TOP""AND OORREOTING FOR CONDUCTION BY CALCULATION. EAOH OF
THE EXPERIMENTAL RUNS WAS OORREOTED BY SUBTRACTING THE APPROPRIATE
RADIATION VALUE FROM THE MEASURED TOTAL HEAT TRANSFER RATE.
THE EXPERIMENTAL PROOEDURE WAS THE SAME FOR All THE SOURCES.
A TEMPERATURE DIFFERENCE WAS APPLIED TO THE FLUID LAYER, THE SYSTEM
WAS ALLOWED TO REACH EQUILIBRIUM AND MEASUREMENTS OF T., T2~ AND Q (OR
ELECTRICAL POWER CONSUMPTION) WERE MADEo
CERTAIN DISCREPANCIES UN EXPERIMENTAL TECHNIQUE ARE WORTHY OF
NOTE BECAUSE THEY BEAR ON THE ACCURACY OF THE DATA.
BOTH GLOBE AND DROPKIN AND SCHMIDT AND SILVESTON, WHO EXPERIMENTED
WITH LIQUIDS, DETERMINED THE HEAT TRANSFER BY MEASURING ELECTRICAL
POWER CONSUMPTION. THE APPARATUS OF THE LATTER WAS ELABORATELY
INSULATED TO GUARD AGAINST HEAT LOSS, AND OORRECTIONS, WHERE NECESSARY,
WER E APPL I ED.. How EVER J GLO BE AND ORO PK I N JI NE I TH ER I NSULA'~'t:O THE I'R
EQUIPMENT NOR CORRECTED FOR lOSS. AT lEAST THE AUTHORS ARE SILENT
ON THIS POINT AND THEiR DIAGRAMS AND DISCUSSIONS REVEAL NO SUCH
PRECAUTlONS. IF SUOH IS THE CASE, THEIR Nu SHOULD TEND TO BE HIGHER
AT PARTIOULAR RA. THE APPARATUS OF DE GRAAF AND VAN DER HELD WAS NOT
INSULATED EITHER, BUT THEY MEASURED THE HEAT TRANSFER RATE BY
MEASURING THE HEAT ACQUIRED BY THE COOLING WATER. IN ADDITION, THEY
CORRECTED FOR EXTRANEOUS LOSSES AT THE SAME T8ME THAT THEY OORRECTED
FOR RADIATION.
... 10 ~
DE GRAAF AND VAN OER HELD HAVE POINTED OUT (12, 13) THAT
MULL AND REIHER CAL~ULATED THE THERMAL COEFFICIENT OF EXPANSION
FOR AIR INCORRECTLY, USING 1/273 [NSTEAD OF 1fT CI\VG.) THEREF"ORE,
THE DATA OF AND REIHER WAS RECALCULATED BEFORE artNG USED IN
THiS CORRELATION" ApPENDIX I! CONTAINS THE DERiVATION OF A ~ I/TAv90
ApPEND1X I I I DESCRIBES THE CALCULATION OF MULL AND REIHER'S NUSSELT
AND GRASHOF NUMBERS USlNG THE CONSTANTS OF DE GRAAF AND VAN DER HELD
THE EXPERIMENTAL DATA O~N AIR (5,12) WAS PLOTTED BY THE ORIGINAL
AUTHORS iN THE fORM OF Nu VS GR, SINOE PR Is NEARLY A CONSTANT FOA
AIR OVER THE TEMPERATURE RANGE OF THEIR EXPERiMENTS. RA FOR THESE
DATA WAS CALCULATED FROM THE OR.lG:SNAL GR (RECALOULATED 1N THE CASE Of'
MUll" po.NO RElHER) USING THE PR THAT ARE PLOTTED VERSUS TEMPERATURE IN
FiG" I .. ALTHOU\3H TAE VARIATION WITH TEMPERATURE IS NOT LARG!::" THIS
WAS DONE TO REDUCE ERROR IN THE CORRELATION~
CORRELATION
THE DATA WAS ANALYZED STAT1ST~DALLY BY THE METHOD OF STEPWiSE,
LtNEAR, MULTIPLE, REGRESSiON (25). LINEAR MULTIPLE REGRESSION
CONSlSTS OF F[NDING BY THE METHOD OF LEAST SQUARES THE FUNCTION
OF THE FORM
THAT BEST FiTS THE DATA. THEREFORE, THE HYPOTHETIOAL FUNCTiON
Nu = A (RA)B 1 (PR)B 2 (A/0 2 )B 3 WAS EXPRESSED IN LOGARITHMIC FORM~
THE LEAST SQUARES BEST F1T lS THAT EQUATION THAT RESULTS IN tHE:
MlNtMUM STANDARD DEViATION OF THE DISTRIBUTION OF EXPERiMENTAL DATA
ABOUT THE REGRESSrON. K MORE DETAILED EXPLANATION OF THE METHOD
AND THE INTERPRETATiON OF THE STATISTiCS IS SEVEN [N ApPENDiX IV.
IN STEPWISE MULTfPLE REGRESSION, THE VAR'ABLES ARE ADDED TO THE
REGRESSION ANALYSiS ONE AT A TIME IN THE ORDER OF THEIR CONTRIBUTION
TO THE GOODNESS OF FIT. THIS PROCEDURE HAS SEVERAL ADVANTAGES. AT
EACH STEP IT PRovIDES THE FOLLOWING:
I. A REGRESSION EQUATiON FOR EACH OF THE VARIABLES
INCLUDED IN THE REGRESSION UP TO THAT STEP.
2. AN ESTIMATE OF THE PROBABLE ERROR IN THE PREDlOTION
OF Y FROM THE REGRESSION EQUATloN~
- 12 -
3. ESTIMATE OF THE PROBABLE ERROR IN THE REGRESSION ~
OOCFFiC~ENTS 8,. 8 2, 8 3 ••••••
4. AN ESTIMATE OF THE SIGNIFICANCE OF EACH VARiABLE
iN THE REGRESSION.
5. REJECT,oN OF INSIGNIFICANT VARiABLES.
THE CALCULATIONS FIRST WERE CARRIED OUT MANU~~LY TO DETERMINE
THE REGIMES INTO WHtCH THE DATA SHOULD SE DIViDED. THEN, FOR ACCURACY.
THE DATA WERE ANALYZED ON AN 1.8.M. 704 COMPUTER USING A STEPWISE.
LINEAR, MULTiPLE REGRESSION PROGRAM DEVELOPED BY THE COMPUTING CENTER
OF THE Esso RESEARCH AND ENGINEERING CO. THE RESULTS OF THE TWO
ANALYSES WERE IN COMPLETE AGREEMENT. TABLE I I 1* LISTS THE EXPERIMENTAL
DATA AND THEIR LOGARITHMS, WHICH WERE THE INPUT TO THE COMPUTER. ALL
THE DATA WERE ASSUMED TO HAVE EQUAL WEIGHT.
RESULTS
THE EXPERIMENTAL DATA CORRELATED IN THREE DISTINCT REGIMES
OF CONVECTION - INITIAL, LAMINAR, AND TURBULENT. THE RANGE COVERED
BY EACH REGIME, THE PARAMETRIC EQUATIONS FOR EACH REGIME, THE
STANDARD ERRCR OF ESTIMATE O~ Nu FROM EACH EQUATION, AND THE STANDA~D
ERROR OF THE COEFFICIENTS. ARE LISTED IN TABLE IV.
*COEFFICIENTS THROUGHOUT THIS PAPER REFER TO THE COEFFICIENTS OF
THE LOGARiTHMIC EQUATION. TH£SE BECOME EXPONENTS OF THE PAR'METERS
IN THE POWER •• W EQUATION.
* TABLE I I I - Ps. 18
TABLE 11 SUMMARY OF REGRESSION ANALYSI.
STD~ ERROR STD.ERROR OF COEF. OF ESTIMATE RA PA ...... %r..;.......,_
LAMINAR ~',IJ006
TURBUl.ENT
THE FOllO.'N8 CONCLUSIONS. CAN BE DRAWN FROM THE ANALYSIS:
I. CONVECTIVE HEAT TRANSFER IN HORIZONTAL LAYERS CAN BE
EXPRESSED ACCURATELY BY THREE DISTINCT EQUATIONS G
2. Nu CAN BE EXPRESSED AS A SIMPLE POWER LAW FUNCTiON
OF RA, AND, iN THE TURBULENT REGIME, OF RA AND PA.
3. A/0 2 is NOT A SIGNU''"lCANT CORRELATING PARAMETER AND
THE ASSUMPTION OF LAYERS OF iNFINITE HORIZONTAl. EXTENT
IS JUSTiFIED.
+7.19 ... 6.73
THE VALUE OF THE CRITtCAL RA CALCULATED FROM THE EQUATiON OF THE
iNITIAl. REGIME is 1640 .. THIS IS IN REASONABLE AGREEMENT WITH THE
ANALYtiCAL VALUE OF 1707 .. 8 (3) AND WiTH EXPERIMENTAL VALUES REPORTED
IN THE LITERATURE. (6, 7, 20, 21)
IN FIG. 3, Nu VB RA OF THE EXPERIMENTAL DATA Is PLOTTED TOGETHER
WITH THE EQUATIONS OF THE THRE[ REGIMES. THE TURBULENT EQUATION ts
PLOTTED WITH VARIOUS PH AS PARAMETERS. IN FIG$ 5 THE TURBULENT REGIME
/ 0853
IS OOMBINED iNTO ONE EQUATION BY PLOTTING Nu PR- VS RA. IT CAN BE
~ 14 ...
SEEN FROM THESE FIGURES THAT THREE REG!NES COMPLETELY DESCRIBE
THE DATA4 CHANGES IN MOOE OF CCNVECTtON WITHIN THESE REGiMES
APPARENTLY DO NOT AFFECT THE HEAT TRANSFER tUNCTION.
iN FIG. 4, THE EQUATIONS OF THE THREE ~EaIMES ARE PLOTTED AS
Nu VB RA, SiMILAR TO FIG. 3, EXCEPT THAT THE EXPECTED TRANSITlGNS
BETWEEN LAMiNAR AND TURBULENT REGIMES ARE SHOWN AS DASH EO LINES. IT
CAN BE SEEN THAT THE TRANSITION FROM INiTIAL TO LAMiNAR REGIME IS
SHARR, WHiLE THE TRANSITION FROM LAMINAR TO TURBULENT REGlME DEPENDS
ON THE PR OF THE FLU1@. PREVIOUS INVESTIGATORS HAYE NOT REALIZED THIS.
FLUIDS HAVING A HIGH PR UNDERGO A SHARP TRANSITION, WHILE THOSE WITH
SMALL PH GO THROUGH A MORE GRADUAL TRANSITION~ THE DATA, ALTHOUGH
SPARSE IN THtS REGION, SUPPORTS THIS OONOLUSION EXCEPT FOR MERCURY, WHICH
INEXPLICABLY DEVIATES FROM iTS EXPECTED PATTERN IN THE DIRECTION OF
LOWER Nu ..
I III FIG.. 5 ,t T CAN B ESE EN T HAT THE D A T A I NTH E T RAN SIT I 0 III Ii E G I {l III
OORRELATES WELL WITH THE REST OF THE REGIME$ EVEN THE LOWER MERCURY
POINTS FALL WiTHIN THREE SiGMA LIMITS AND~ THEREFORE, CANNOT BE
E~CLUD[O ON STATISTiOAL GROUNDS. THEREFORE, THE TURBULENT REGRESSION
EQUATION CAN BE USED TO PREDIOT Nu RN THE TRANSITION REGION.
DISOUSSION OF THE ANALYSIS OF DATA
TABLE V* SUMMARIZES THE COMPUTER PRlMT-OUT or THE STEPWISE
REGRESSION ANALYSIS .. AN EXPLANATION or THE STATISTIOS IS GIVEN IN
ApPENDIX IV.
* PG III 27
SllNIFrCANCE OF THE CORRELATING PARAMETERS
THE COMPUTER ANALYSIS SHOWED THAT All OF THE CORRELATING
PARAMETERS WERE STATISTICALLY SIGNIFICANT IN ALL OF THE REGIMES;
t.E., THE BEST CORRELATION INCLUDED ALL THREE PARAMETERS. HOWEVER.
THE IMPORTANT STATISTICAL MEASURE IS THE DEBREE OF SiGNIFICANCE OF
THE PARAMETERS, WHICH is GiVEN BY THE COEFFICIENTS OF DETERMINATaON.
THESE ARE TAB'I.H.I.TE:D IN TABLE vt~· AT EACH STEP OF THE REGIH:SSIONo
THE OOEFFICIENT OF OETERMINATION (EQUAL TO THE SQUARE OF THE MULTIPLE
CORRELATION COEFFICIENT) MEASURES THE PER CENT OF THE TOTAL VARIANCE
aN Nu ACCOUNTED FOR BY THE REGRESSION AT EACH STEP; i.E., AS EAOH
PARAMETER IS ADDED TO THE REGRESSION ANALYSIS. THE RESiDUAL REPRESENTS
THE VARIANOE IN Nu THAT IS NOT ACCOUNTED FOR BY THE REGRESSION
(SCATTER OF DATA ABOUT THE REGRESSION LlNE) AND IS ASSUMED TO BE
RANDOM ERROR. IT CAN BE SEEN FROM TABLE VI THAT RA IS VERY HIGHLY
SIGNiFICANT WHILE~ EXOEPT FOR ONE OASE, PR AND A/o 2 A~E LESS THAN THE
RES IDUAl6WE CONCLUDE, THEREfORE, THAT PR AND A/o 2 ARE NOT S IQNlf HCAtH
PHYSICALLY AND THAT THE CORRELATION FOUND BY THE REGRESSION IS
SPURIOUS; I.Ee, DUE TO CHANCE. THE VARIANCE REDUCTION ATTRIBUTED TO
PR AND A/D 2 ~s REGARD€D AS ERROR AND THE PROPER REGRESSION 18 THE
RESULT Of STEP ONE, INVOLVING ONLY RA. THE EXCEPTION is PR IN THE
TURBULENT REGIME. ITS SIGNIFiCANCE Is ALSO SHOWN BY THE SUBSTANTIAL
REDUCTION OF THE STANDARD ERROR OF ESTIMATE AFTER PR rs ADDED TO THE
REGRESsiON (TABLE V). FOR THIS REGIME WE CONCLUDE THAT THE PROPER
CORRELATING EQUATION IS THE RESULT OF STEP TWO~ iNVOlViNG RA AND PR.
THE MULTIPLE CORRELATION COEFFiCIENTS FOR ALL THREE REGIMES
REFLECT VERY HIGH DEGREES OF CORRELATiON. THEREFORE. WE CONCLUDE
THAT THE FUNCTIONAL RELATIONSHIP NW = F OF (RA) OR F OF (RA)(PR)
AOCURATELY DESORIBES THE PHYSiOAl PHENOMENON~ THERE is NO EWtDEHOE
iN THE GRAPHS OF THE DATA TO SUGGEST THAT THERE ARE MORE REGIMES THAN
THE THREE ALREADY DESCRIBED.
STANDARD ERROR ESTIMATE
THE STANDARD ERROR OF ESTIMATE as THE STANDARD DEVIATiON OF THE
DISTRIBUTION OF EXPERIMENTAL VALUES OF Nu ABOUT THE REGRESSIOH. fT
HAS THE USUAL SI~T~STtOAL S1GNIF10ANOE; I.E.# APPRoxiMATELY 2/3 OF
THE DATA (S EXPECTED TO FALL WITHIN t ONE STANDARD ERROR OF ESTIMATE
(S~GMA) OF THE VALUE PREDiCTED BY THE REGRESS~ONo IN THIS CASE, THE
STANDARD ERROR OF ESTIMATE IS GIVEN IN LOGARITHMIC FORM BY THE
REGRESSION ANALYSIS. IN ORDER TO EXPRESS (T IN TERMS OF Nu IT MUST
BE STATED AS A PERCENTAGE. THIS ACCOUNTS FOR THE FACT THAT BTS POSITIVE
AND NEGATlVE NUMERiCAL VALUES ARE DIFFERENT. IF IT is DESIRED TO
CALCULATE TWO OR THREE S$GMA LIMITS, THEN THE LOGARITHMIC VALUE MUST BE
USED TO CALCULATE THE APPROPRIATE PERCENTAGE~
ERROR OF THE REGRESSI0N
THE STANDARD ERROR OF A COEFFICIENT IS A MEASURE OF THE PROBABLE
ERROR (+ ONE SIGMA) OF THE COEFFICIENT ON THE AVERAGE. HOWEVER, IT IS A .... SPECIFIC PROPERTY OF REGRESSIONS THAT THIS ERROR IS NEGLIGIBLE AT THE
CENTER OF THE RANGE OOVERED BY THE EXPERIMENTAL DATA AND LARGE AT EITHER
ENDo (SEE ApPENDIX IV) THEREFORE, THE REaRESSION EQUATION SHOULD NOT BE
EXTRAPOLATED BEY0ND THE RANGE OF THE DATA WITHOUT FURTHER S~BSTANTIATDON.
TABLE 11
COMPARISON OF THE ARATUS USED By THE SOURCES OF THE DATA
DE GRAAF AND VAN DER HELD
MuJ..'l. AN D RE UIER
SHAPE OF FLUID CHAMBER SQUARE PRISM REOTANGULAR PRISM
AREA OF FLUaD CHAMBER, M2
,,0'. AM. OR LENGTH;\VIOTH,
MM
CONTAINER WALL MATERIAL
WAS THERE A GUARD HEATER?
METHOD OF HEATING
METHOD OF COOLING
MAXe TEMPo, DIFFERENOE, C
WERE LIQUIDS DEGASSEO?
METHOD OF PLATE TEMP.MEASUREMENT
METHOD OF HEAT TRANSF'ER MEASUREMENT
RADiATION CORRECT .'ON
430/430
6.9.22 .. 9
GLASS
No
No
ELEOTRIOAL
WATER
146
100
3 THERMOCOUPLESlN
10 10/612
12-196
WOOD
YES
YES
ELEOTRIOAL
WATER
146
29
15 THERMOCOUPLES IN
HOT PLATE; HOT PLATE; MEAN COOLING 6 THERMO-WATER TEMP. COU PL ES tiN
COLD PLATE
HEAT ACQuiRED ELECTR.OA~ BY 'OOOL.'I'NS POWER CON-WATER SUMPTION
YES
GLOBE AND DROPKIN
CYLINDER
0.0'3,0.014
127,134
35 ... 66
VI Pl EX! Ea. AS U ,
"PYFU:X II
No
No
ELECTRIOAL
Au~, WATER
92
27
YES
1 THERMOCOU Pl.£: IN EACH PLATE
ECTRIOAl POWER CONSUMPTION
SCHMIDT AI«l SU .. VEaTON
C'(ll NO ER
199
1.45 ... 13.0
&tPLEX'f'SLAsM
YES
YES
ELECTRIOAL
WATER
70 (EST.)
50
YES
3 THERMOCOUPLES IN EACH PLATE
EL E<1lTR • CAL POWER CONSUMPTiON
TA
BtE
J..
!!
DA
TA
~in)UROES:
MR
:::::
MU
LL
AN
D
RE
t H
ER
(5
) S
8
::::
SCH
MR
DT
A
ND
S
ILV
ES
TO
N
(20
) G
D
::::
GL
OB
E
AN
D
DR
OP
KIN
(2
1)
DV
::::
: D
E
GR
AA
F A
ND
V
AN
D
ER
H
EL
D
(12
)
TU
RB
UL
EN
T REGIME~
RA
105
TO
J0
9
SO
UR
OE
F
LU
,I!)
R
A
X
10
-5
PR
A
/g2
-~
LO
G_B
I<
j...o~
1a ~e
R~
1_oGAL~?
Lo
a
NUIllP
MR
A'iR
6
.32
.6
99
6
9.4
6
.67
5
.80
07
-.
15
55
i 4
18
4 1
4 .$
24
1
5 .. 0
3 .7
04
6
9.4
5
.75
5
.70
16
...
. 152
4 "
,,75
97
6.4
.6
93
6
9.4
6
.80
5
.80
62
-.
1593
It
418
325
6.0
9
.68
9
69
.4
6.5
5
5.7
84
6
.....
1618
Il
4181
62
52
.6
.70
0
16.,
2 1
4.
15
6.7
21
0
....
1549
i .
20
95
BC
II149
2 ro
2
5.8
.7
06
1
6.2
1
0.8
3
6.4
i 16
-.
1512
tt
i 1
.03
34
S8
WA
TE
R
2.4
7
6.0
7
30
6.
5.5
0
5.3
92
7
.78
32
2
.48
57
.7
40
4
3.8
0
5.4
3
30
6.
6.5
9
5.5
79
8
&73
48
tt
.81
89
HE
PT
AN
E
1.2
5
6.2
6
32
90
. 4
.23
5
.09
69
.1
96
6
3.5
11
2
.62
63
Suo
i O
ON
E 0
u.
AK
3 1
.86
3
5.0
1
22
4.
5111
30
5.2
69
5
1.5
44
1
3.0
87
8
.72
43
GO
1000
C
.S.
S U
. t
CO
NE
O
n.
Lt5
i 8
75
0.
4.9
1
8.3
8
5.
1790
3
.94
20
.6
9' I
.9
23
2
2.0
1
85
00
. tt
9
.53
5
.30
10
3
.92
94
It
.9
79
J 6
.20
7
36
0.
It,
JO.9
0 5
.79
24
34
186
69
It
L,0
374
.
150
C.S
. S
H.
i co
~u:
O·n.
2n
o
16
. n
48
.8
8.3
24
3
t 1&
204
J It
1
.68
84
67
60
t I.
n
66
.9
8.8
30
0
1.0
41
4
n 1
.82
54
~
49
20
1
2.5
~'
59
.1
8.6
92
0
1.0
96
9
tt
1.7
71
6
. 18
10
1541
' It
4
2.9
8
.. 257
7 1
.11
61
tt
1
.63
26
TABL
E II
r
(Co
rnio
) --
TURB
ULE
NT
RE
GlU
E:
RA
10
5 TO
10
9
SOU
ROE
FL
UIO
R
A
x 10
... 5
PR
A
/o,2
N
u LO
G RA
LO
G PR
lO
G
ALD
2 lo
t N
u.
--
II'
i y--~
GO
IS
O
G .. S
. 4
.91
.. 6
9 I
I S
U ••
CON
E ou
. 11
40
i5~5
it
34.,
3 8
.05
69
1
.19
03
at
, .5
35
3
39
1
17
. It
3
2.7
7
.59
22
f .
23
04
"
1.5
14
6
-"
280
17
. f~
3
0.4
7
.44
72
t.
23
04
it
J .
48
29
"
49
5
17
. n:
34
.3
7.6
94
6
J .2
30
4
It;
1.<l
t535
3 76
0 1
5&
5
It;
41 ..
0
7.8
80
8
1.'
90
3
It
1.6
12
8
--.
954
15
. n
39
.8
7.9
79
6
I. f
761
n 1
.59
99
-.
80
5
17
. it
4
L.6
7
.90
58
t .
23
04
It
1
.61
91
-.
-.
15}2
15
.. It
· 4
6.5
8
" 17
96
'.1
76
' n
t.6
q7
4
.. 55
27
1241
It
6
5.0
8
.. 742
5 T
.07
92
it
1
.8 '2
9
" ..
-..
WAT
ER
2490
20
3 n
33
.2
8.3
96
2
.36
17
lic
L~
!.52
' 1
(t)
1934
2 .
• 5
a~
3041
9 8
.28
65
.3
97
9
It
1.4
90
0
-. -
J384
2l1
li5
at 2
9.2
8
. J4
1 I
.39
79
.,
f .4
65
4
-.
890
3.5
if
2
8.4
7
.94
94
.5
44
1
It
1.4
53
3
579
3.8
it
2
7.f
7
.76
27
41
5798
It
f .
. 43
30
-
388
4.3
It
2
4.3
7
.. 588
8 .6
33
5
n'
f.3
85
6
183
4.9
i.
19
.3
7.2
62
4
.69
02
it
1
.28
56
-.
"
80
cA
5.7
tt
1
7.8
6
.90
53
.. 7
559
It
r.25
04-
.-
42
.8
5.5
at
15
.4
6.6
31
4
0740
4 ••
f 0 H
375
f 1
6.8
5
.7
tt
19
.0
7.0
67
4
.75
59
It
1
.27
88
-.
50
G.8
. S
I L
I co
NE
0 u.
1
12
. 2
99
. It
2
4.1
7
.04
92
2
.47
57
It
1
.38
20
19
i 2
45
. n
26
.5
7.2
81
0
2.3
89
2
It
f .4
23
2
-.
40
.7
36
4.
It
16
.7
6.6
09
6
2.5
61
1
n r .
22
27
--
t L,8
4
02
. It
12
.8
6.0
7}
9
2.6
04
2
i~
1.1
07
2
TABL
E II
I (C
ON
TrD
) -
TURB
ULE
NT
RE
GIM
E:
RA
105
TO
109
SOU
RCE
FLU
ID
RA
X
10
.. 5
.E.!!
A/D
2 ~
LOG
RJI).
LO
G PR
LO
G~A/
D2
GO
M
EROU
RY
9.8
7
.02
32
3
.19
4
.91
5
.99
43
...
J .6
34
5
4115
038
.69
t 1
35
. .0
21
8
It
6~63
6.5
44
J ....
, .6
61
5
•• ~ 8
2 J.
5 3
076
.0
23
5
it
3.0
7
5.5
75
2
... 1
.62
89
it
.4
87
1
~
6.1
4
.02
33
at
3
.87
5
.78
82
-1
.63
26
It
.5
87
7
,
16
.5
.02
27
tt
5
.61
6
.21
75
....
1.6
44
0
it
.74
90
-
9.5
7
.02
29
Q
; 5
.22
5
.98
09
...
1.6
40
2
It
.. 7J7
7
~
7.9
2
.02
31
I~
: 4
.38
5
.89
87
...
f .6
36
4
UI
.64
15
-
4.6
7
.02
32
tt
30
:67
5.6
69
3
... f .
63
45
it
.5
64
7
13
.2
.02
28
n
6.2
0
6.
i 206
...
1.64
21
tt
.79
24
~
220
.02
22
tt
7
. r I
6
.34
24
...
1.6
53
6
lct
.85
19
~
34
.3
.. 0
21
4
tt
8.0
8
6.5
35
3
.. 1
.66
96
n
.90
74
1;
35
.8
.02
12
tt
7
.87
6
.55
39
...
1.6
75
7
ft
.89
60
I\
,)
36
.2
.02
12
n
7.6
4
6.5
58
7
.... 1
.67
57
tt
.8
83
1
0 -
5t.
.0
20
3
if
8.3
2
6.7
07
6
... t
.69
25
u
.92
01
11
-9
.05
.0
23
1
Q
5.5
6
5.9
56
6
... 1
.63
83
tt
.7
45
' 5
.47
.0
23
3
tt
4.5
0
5.7
38
0
... f
.63
26
It
.6
53
2
-3
.73
.0
23
5
n 3
.25
5
.. 571
7 ....
1..6
289
Ii
.51
19
1
.63
.0
23
6
tt
'.8
7
5.2
12
2
... 1
.62
71
It
02
718
-2
.02
.0
23
6
n 2
.55
5
.30
54
....
1.6
27
1
lit
.40
65
2
4.7
.0
22
7
t 1
.28
7
.68
6
.39
27
...
1.6
44
0
1.0
52
3
.88
54
5
1.6
.0
20
0
It
9.7
2
6.71
.26
... t
.69
90
u
.98
77
1
3.0
.0
23
8
U:
6~02
6.1
13
9
-1.6
23
4
If
1117
796
56
.0
e02
35
tt
9;5
9
6.7
48
2
-1.6
28
9
it
.98
18
92
08
.02
32
u
10
.9 t
6".9
676
... 1
".6
34
5
tt
J .0
37
8
16
8.9
1110
222
n
J3.8
6
7.2
27
6
.... 1
.65
36
It
1
.14
18
3
33
.7
.02
19
tt
1
6.8
8
7.5
23
4
-1.6
59
6
It·
, .2
27
4
TAB
LE II
I (C
ON
TfD
) -
.
lAM
iNA
R
RE
GIM
E:
RA.
3000
TO
10
5
SO
UR
GE
F
LU
,~o
&
PR
p.J
02
Nu
LOB
RA
-L
oe
PR
LO
G
A/0
2
1..9
6 N
Q.4
> -
--O
V A
IR
3560
07
10
76
1.
1.5
2
3.5
51
4
.... 1
487
2.8
81
4
.18
18
6
49
0
0707
7
61
. 2
. 19
3
.81
22
"'
".15
06
2.8
81
4
.34
04
8
28
0
.70
5
76
1.
2.2
6
3.9
18
0
-.'
5'8
2
.88
14
.3
54
1
94
20
.7
03
7
61
. 2
.50
3
.97
40
-.
15
30
2
.88
14
03
979
40
00
.7
09
7
61
. '-'
,63
3 .. 6
021
M.1
494
2.8
81
4
.21
22
52
80
.70
8
76
1.
1.8
0
3 .. 7
226
.....
'50
0
2.8
8'4
.2
55
3
70
70
.7
07
76
f ..
2.0
2
3.8
49
4
-.1
50
6
2 .. 8
81
4
.30
54
8
72
0
.70
4
76
1.
\2 ..
25
3 .. 9
40
5
-.1
52
4
2 .. 8
81
4
.35
22
lO
tOO
11
1702
76
t II
2
.37
4
.00
43
...
. 153
7 2
.88
14
.3
74
8
I fO
OO
• 6
99
76
L •
2.3
6
4.0
41
4
....
J552
2"
,881
4 .. 3
729
11
J i
100
• 698
76
1 •
2.3
8
4.0
45
3
......
J561
2
.88
14
.3
76
6
rv
2340
0 • 7
09
23
1 •
3.3
3
4.3
69
2
.... 1
494
2.3
63
6
.52
24
--
40
60
0
.70
6
231
Ell
3.7
2
4.6
08
5
.....
1512
~
.57
05
i .
5
20
00
•
703
231
• 4
.01
4
.71
60
...
.. 1
53
0
•• .6
03
1
-58
100
11'47
02
23
1.
3 .. 8
1 4
.. 76
42
..
.. 1
53
7
It
$580
9 6
32
00
",
700
23
1.
3 .. 8
8 4&
8007
-l
IP 15
49
It
.58
88
-
6660
0 .6
97
2
31
. 4
.06
4
.. 823
5 -
.. i5
68
It
.. 6
085'
68
100
.69
6
23
1.
4.0
2
4.8
33
2
1"J.
1574
tt
.6
04
2
'06
00
.. 7
01
76
1.
2&34
4
.02
53
...
. 154
3 2
.88
14
.3
69
2
1660
0 .7
10
2
31
. 2
.82
4
.22
01
...
. 148
7 2
.36
36
.4
50
2
2370
0 .1
09
2
31
. 3
.48
4
.37
48
..
.. (
It 1
49
4
It
.54
16
32
600
• 708
23
1 •
3.5
7
4.5
52
7
O\\lj
f/ll5
00
" .5
52
7
4410
0 .7
06
2
31
. 3
.55
4
.. 644
4 ~.1520
n 4155
02
5470
0 .7
03
2
31
. 3
.93
4
.73
80
""
".15
30
It
.59
44
65
100
.70
0
23
1.
3.9
J 4
.81
36
..
..15
49
It.
.59
22
MR
AiR
9
66
0
.69
5
10
90
. 2
.15
3
.98
50
1I
iii.1
580
3.0
37
4
.33
24
9
64
0
.69
5
5~&j, .
2.1
3
341198
41
....
. 580
2
.73
48
.3
28
4
9430
.6
95
2
71
. 2
.23
3
.97
45
...
. 158
0 2
.43
30
.3
48
3
9470
• 6
95
la6
• 2
.23
3
.97
64
-.
15
80
2
.13
35
.3
48
3
8450
0 .7
0 l
25
9.
3.6
9
4.9
26
9
.... 1
543
2.4
13
3
.55
51
8
54
00
06
88·
25
9.
3.6
4
4.9
3J5
...
. 162
5 It
.5
61
J
8870
0 .6
91
2
59
. 3
.85
4
09
47
9
.... 1
605
It
.58
55
9
12
00
06
95
25
9.
3.8
1
4.9
60
0
....
1580
u:
.58
66
lAB
! E
!
U
(CO
NT
'0
)
LA
MiN
AR
R
EGIM
E:
RA
30
00
T
O
105
SO
U R
OE
F
LU
ILl
RA
PR
A
/02
Nu
LO
G RA
L
OG
P
R
LOG
_A/D
_?
LO
G
Nu
.• -
--
MR
AIR
4
55
00
0
70
6
25
9.
a~44
4.6>
580
......
i51
2
2.4
13
3
.53
66
8
60
00
.6
95
1
34
. 3~76
4.9
34
5
......
15
80
2
. ',2
7 T
.5
75
2
S0
50
0
.6S
8
13
4.
3 • .4
7
4.9
05
8
.....
1625
tt
.5
40
3
S5
WA
TE
R
35
rO
6.0
5
32
90
1
.74
3
.54
53
!i
t7S1
8 3
415
17
2
.24
06
3
43
0
6.0
6
32
90
1
07
5
3C15
353
$7
82
5
3.5
17
2
.24
30
5
33
0
5.5
1
32
90
2
.03
3
.72
67
.7
41
2
3.5
17
2
.30
75
3
3ro
3
.07
3
29
0
i .8
2
3.5
19
8
.48
71
3
.5'7
2
e26
01
6
05
0
3.7
9
32
90
2
.20
3
.. 7
81
8
.57
86
3
.51
72
.3
42
4
93
80
3
.59
3
29
0
2.3
1
3.9
72
2
.55
51
3
.51
72
.3
63
6
I
93
80
3
.59
3
29
0
2.3
0
3.9
72
2
.55
51
3
.51
72
..
36
11
I\
)
47
40
6
.57
12
24
2.0
0
3.6
75
8
&81
76
3.0
87
8
.30
10
I\
)
65
20
6
.46
1
22
4
2.I
S
3.8
'42
E
lJ81
02
3.0
87
8
$<33
85
11
92
80
6
.40
1
22
4
20
45
3
09
67
6
.80
62
3
C108
78
.3
89
2
53
90
4
.34
1
22
4
2.0
3
3.7
31
6
.63
75
3
.08
78
.3
07
5
80
00
4
.33
1
22
4
2.'
5
3.9
03
J
.63
65
3
.08
78
0
33
24
i3
10
0
409
8
32
90
2
.42
4
<I> I
173
.69
72
3
.. 51
72
.3
83
8
13
00
0
4.9
9
32
90
2
$4
2
4.1
13
9
.69
81
3
.51
72
.3
83
8
1540
0 64
1130
1
22
4
2.6
0
4.
J 87
5
.79
93
31
1108
78
.41
50
2
30
00
6
.16
.2
24
2
.. 7
4
4.3
61
7
~7S96
3.0
S7
8
.43
78
3~200
5.9
4
.12
24
2&
90
4.4
94
2
.77
38
3
.08
78
41
46
24
41
80
0
5.6
3
12
24
3
.17
4
.62
12
.7
50
5
3.0
87
8
.50
r f
66
00
0
5.1
6
12
24
3
.59
4
.8J9
5
.. 7
12
6
3.0
87
8
.55
5 J
1
71
00
4
.37
1
22
4
2.5
8
4.2
33
0
.64
05
3
.08
78
8
4 J
'6
24
50
0
43
41
1
22
4
2.7
1
4.3
89
2
.64
44
-3
.08
78
.4
33
0
36
50
0
4.3
9
12
24
3
.. 03
4
.56
23
.0
42
5
3.0
87
8
.48
14
-5
63
00
4
.30
'2
24
3
.48
4
.75
05
.6
33
5
3.0
87
8
.54
16
7
86
00
4
.44
&
224
3.7
5
4.8
95
4
.64
74
3
.0i7
8
.57
40
3
96
00
6
. r 3
3
06
3
.25
4
.59
77
.7
87
5
2.4
85
7
.51
19
79
70
0
6.J
7
30
6
3.6
6
4.9
01
5
.79
03
2
.48
57
.5
63
5
97
80
0
6.2
4
30
6
4.
i8
4.9
90
3
.79
52
2
.48
57
.6
21
2
TABL
E Ilf
(C
ON
T'O
) --~-
LAM
iNA
R R
EGIM
E:
RA
30
00
TO
10
5
SOU
RCe:
FL
UID
RA
PR
A
/02
Nu
LOG
RA
LOG
PR
LO
G AL~z
to_
ILN
!tlL
-
--
SS
HEP
TAN
E 32
10
6.2
3
1457
0 1
.61
3
.50
65
.7
94
5
4.1
63
5
.20
68
32
30
6.2
3
1457
0 J .
60
3
.50
92
.7
94
5
4.1
63
5
.20
41
38
30
6.2
4
1457
0 1
.73
3
.58
32
.7
95
2
4.
'63
5
.23
80
4
47
0
6.2
3
'45
70
1.
.82
3.6
50
3
.79
45
4
. 16
35
.2.6
0 t
60
40
6
.26
14
570
L,9
7 3
.78
10
.7
96
6
44i i
63
5
112 9
45
9460
6
.28
14
570
2.2
J 3
.97
59
&
7980
4
.16
35
.3
44
4
1480
0 6
. J8
14
570
2.3
8
4.1
70
3
.79
10
4
.16
35
.3
76
6
2350
0 6
.22
3
29
0
2.5
2
4.3
7 t
J &
7938
3
.51
72
.4
.0 ',
4 29
400
6.2
3
3290
2
.66
4
.46
84
.7
94
5
3.5
17
2
.42
49
35
900
6.2
3
32
90
2
.84
4
.55
5J
.79
45
3
.51
72
.4
53
3
J
5040
0 6
.20
3
29
0
3.2
3
4.7
02
4
.79
24
3
.51
72
.5
09
2
I\:)
8370
0 6
.26
32
90
3.7
0
4.9
22
7
.79
66
3
.51
72
.5
68
2
w
It
SS
SIL
ICO
NE
O
n.
AK
3 33
60
35
.8
1360
0 1
.73
3
.52
63
1
.55
39
4
. '3
42
.2
38
0
61
40
3
5.2
32
90
2.2
1
3.7
88
2
r .5
46
5
3.5
17
2
.34
44
8
02
0
3507
32
90
2.3
0
3.9
04
2
1.5
52
7
3.5
17
2
4113
6J7
9750
3
5.8
32
90
2.4
1
3.9
89
0
t .5
53
9
3.5
17
2
.38
20
19
500
35
.4
32
90
2
.89
4
.29
00
1
.54
90
3
.51
72
.4
60
9
2820
0 35
4112
3290
3
.15
4
.A5
02
! .
. 54
53
34
1151
72
.49
83
31
300
35
.5
1224
3
.15
4
.49
56
1
.55
02
3
.08
78
.4
98
3
4960
0 3
6.7
12
24
3.5
7
4.6
95
5
J .5
52
7
3.0
87
8
.55
21
93
900
35
.2
1224
4
.39
4
11197
27
1.5
46
5
3.0
87
8
.64
25
SS
S it
. le
ON
E
On.
AK
35
0 32
00
2630
.. 62
8 J .
60
3
.50
52
3
.42
00
2
.79
80
.2
04
1
51
60
2
58
0.
628
2.0
0
3.7
12
6
3.4
1 1
6 2
.79
80
.3
0.0
4
50
0
36
50
. 30
6 1
.96
3
.65
32
3
.. 56
23
2
.48
57
.2
92
3
64
90
3
56
0.
306
2.2
3
3.8
12
2
3.5
51
4
2.4
85
7
.34
83
9
14
0
3400
$ 30
6 2
.. 45
3
.96
10
3
.53
15
2
.48
57
.3
89
2
7310
3
08
0.
18
1.6
2 .
. 35
3.8
63
9
3.4
88
6
2.2
59
J .3
71
1
1210
0 3
05
0.
18
f.6
2
.63
4
.08
28
3
.48
43
2
.25
91
.4
20
0
17
00
0
30
60
. 1
81
.6
2 .. 9
5 4
.. 230
4 3
.48
57
2
.25
91
.4
69
8
2280
0 3
12
0.
, 81
.6
3 ..
15
4.3
57
9
3 .. 4
94
2
2 .. 2
591
&49
83
2920
0 3050~
181
.. 6
3.3
9
4.4
65
4
3.4
84
3
2.2
59
1
.53
02
TA
BtE
1ll
(C
ON
TIO
)
LAM
iNA
R
RE
GIM
E:
RA
30
00
TO
10
5
SOU
RO
E F
LU
IO
&
f!t
A/0
2
~
1~()
(lLB
~,,-
lOG
,P
R
1-_
oG
A/o
2
L~()
G J'
If!J
~
58
E
TH
YL
EN
E-
GLY
CO
L 3
94
0
12
60
32
90
1
,,8
3
3.5
95
5
2.1
00
4
3 .. 5
17
2
.. 26
48
3
18
3
52
.5
32
90
f e
68
3
4115
02
4
f .7
20
2
3.5
17
2
",2
25
3
44
60
5
3.5
3
29
0
1.9
4
3.6
59
3
1411
72
84
3
.,5
'72
.2
87
8
79
60
5
3.0
3
29
0
2.2
8
3.9
00
9
1.7
24
3
3.5
'72
.3
57
9
40
80
1
36
. 1
22
4
, .8
6
3.6
)07
2
.13
35
3
1/0
87
8
.26
95
4
59
0
13
3.
12
24
2
.03
3
.66
18
2
" 1
23
8
3.0
87
8
.30
75
6
62
0
12
0.
12
24
2
.24
3
.82
09
2
.07
92
3
,,0
87
8
.35
02
1
21
00
5
4.0
3
29
0
2.6
1
4.0
82
8
1.7
32
4
3.5
17
2
04
16
6
'30
00
1
34
. 1
22
4
2.6
5
4.
f 1
39
2
. f 2
7 J
3.0
87
8
,,4
23
2
21
50
0
13
3.
30
6
3.0
5
4.3
32
4
2" 1
238
2.4
85
7
.48
43
•
29
60
0
13
0.
30
6
3.2
8
4.4
71
3
2.
I 13
9 2
.48
57
.5
15
9
43
30
0
13
4.
30
6
3.6
4
4.6
36
5
2.1
27
1
2.4
85
7
.56
J t
80
60
0
113
2.
30
6
4 ..
23
4
.90
63
2
.12
06
2.
4815
7 .6
26
3
~
TA
BtE
'i I
(C
ON
T'O
) ----
INIT
IAL
R£GIME~
RA
16
00
TO
RA
3
00
0
SO
UR
ct
FL
utD
~
111
A/o
£ ~
Los
R
A
LOG
PR
LOG
_AjD
2 LO
G N
u.
DV
A
:Ul
7590
.6
99
25
40,.
1"
, 13
3.2
01
4 ~.1555
3.4
04
8
e05~H
SS
WA
TER
1630
4
.. 36
'45
70
. 1
.02
3
.21
22
.6
39
5
4.1
63
5
4100
86
f7
70
4.3
7"
at
1.0
4
3.2
48
0
/116
405
II
ilO
J70
1870
4
.37
tt
1
.06
3
.27
18
.6
40
5
It
.02
53
-
1940
4
.. 37
u
t.0
9
3.2
87
8
.64
05
n
.03
14
20
90
4.3
5
n L
.l7
3
.32
02
.6
38
5
rt .0
68
2
"
2230
44
.41
n 1
.24
3
03
48
3
.64
44
n
.09
34
-
-17
00
6.5
4
32
90
. 1
.05
3
.23
04
.8
15
6
3.5
17
2
.02
12
22
60
4.1
6
n 1
.37
3
.35
4 i
.61
91
H
.1
36
7
2650
6
.03
n
1.5
5
3.4
23
2
.78
03
l~
.1
90
3
2660
4
. J4
at
f.
58
3
.42
49
06
170
ti
.19
87
N
(}
1
ET
HY
LE
NE
G
LY
CO
L
1680
J3
3.
3290
f @
O J
3.2
25
3
2.1
23
8
3.5
17
2
.00
43
It
J820
1
30
. t~
1
.06
3
.26
01
2
.. I 1
39
tV
.02
53
f9
00-
1341
'11
tt
1. f
0
3.2
78
8'
2.1
27
' n
l1'li0
4 21
80
13
0.
n J .
22
31
1133
85
2.
I f 3
9 n
",08
64
-.
2240
1
33
. n
I", 3
1 3
1135
02
2.1
23
8
It
e.,
73
.
2670
1
29
. tt
1
.50
3
.42
65
2
. r 1
06
~
.176
M
1860
5
2.5
It
I .
. to
3"
,269
5 1
.72
02
tt
.0
41
4
2480
5
3.0
n
1.4
. 3
.39
44
, '
" 724
3 It
.1
49
2
-27
80
13B
. 12
24
1.5
3
3.4
44
0
2.1
39
9
3.0
81
8
\t.B
47
HE
PT
AN
E
2240
6
.2f
1457
0 t .
. 29
3&35
02
.79
31
4
.16
35
Co
r f
06
2720
6
.23
tt
1
.47
3
.43
46
.7
94
5
tt
.16
73
S n.
IOO
NE
0
U.
AK
3 f7
30
35
.6
1360
0 1
.05
3
.23
80
1
.55
14
-4
,13
42
.0
21
2
1960
35
",,7
It
I. 1
6 3
.29
23
1
.55
27
"
.06
45
23
80
36 ..
0 it
1
.38
3
.37
66
1
.55
63
n
.13
99
..
TABL
E il
l. (
Cor
-JT
te)
INIT
IAL
R
EG
IME
: R
A
1600
TO
R
A
3000
SG
UR
OE
F
LU
iD
~
PR
A/o
2 ~
J~OG
H
A
LOG~ p~
l..o
~Ajp
3_
1,._
0_8,
N
U_9
-
S8
SIL
ICO
NE
O
IL
AK
35
0
1880
3
66
0.
628
0 1.
.09
3.2
74
2
3.5
63
5
2.7,
980
4103
74-
2170
36
40 ..
u f .
26
3
.33
65
3
.56
11
It
e
'00
4
2680
3
52
0.
tI
1..
49
3
.42
81
3
.54
65
It
.1
73
2
~
-17
60
26
30
. tt
.9
9
3.2
45
5
3.4
20
0
tt
1iI0
000
2460
2
63
0.
it
f.3
t 3
.39
09
3
.. 42
00
II
•
I f 7
3
1890
3
64
0.
30
6.
1.0
9
3.2
76
5
3.5
6 J
I 2
.48
57
,0
374-
29 f
O 3
67
0.
tt 1
.59
3
.46
39
3
.. 56
47
"
.20
14
II f';.)
C
h I
TA
BL
E!
ST
EP
V
AR
IAB
LE
M
UL
TiP
LE
S
y FX
A
T R
EG
RE
SS
iON
C
OE
FF
ICIE
NT
S
$T04
jl E
RR
OR
O
F C
ON
ST
AN
T
EN
TE
RIN
G
CO
RR
EL
AT
ION
E
NT
RY
R
EG
RE
SS
ION
CO
li<!
RE
GR
ES
SIO
N
CO
EF
F J
,e l
EN
T
EF
F'l
CIE
NT
S
131
132
83
8.
82
8
3 lO
G
A
iNIT
IAL
R
EG
IME
(N
:::
32
2
0 .0
65
9
J R
A
.96
24
.0
18
2
37
6.7
.. 8
16
.0
42
....
2"'
62
4
2 PR
.9
69
0
.01
37
6
.1
.83
5
,... .. 0
07
.. 0
40
.0
03
....
2.6
75
3
A/B
2
.98
03
.0
13
7
ISe8
.8
2J
..... 0
19
...
. 031
11
1032
.0
04
1t
00S
1i!l
!2~4
98
LA
MiN
AR
R
EG
IME
(N
=
107
f\
.)
""'I
0 III
f 2
19
.B
I R
A
.96
91
.0
30
2
16
19
.1
e256
.0
06
""
l1li6
56
2
A/f
J2
.. 973
6 Q
028
J 1.
7.8
.24
1
"'*;\1
1024
tl!
007
.00
6
... 4i
l524
3
PR
,97
38
.,
0270
9,
.1
.24
7
.. 00
8
.... 0
22
.0
07
,,
002
.00
6
"'111
559
TURB
ULE
NT
Re:
G1M
E
(N
;;t 6
6)
0 .3
91
3
I R
A
.93
20
~1447
4 t
( .0
e3
48
$O
J7
~1.260
2 PR
41
9920
.. 0
55
4
373
.. 3
.3J6
.0
85
(//
007
.00
4
... 1 .
. 05
0
3 A
/e2
.99
0S
.0
54
1
4 ..
1 .. 3
22
.. 08
3
.02
4
1!I!0
07
.00
4
.01
2
"..1
.. 1
10
Sy
==
STO
.. E
RR
OR
O
F
ES
TIM
AT
E
OF
N
u F
RO
M
RE
GR
ES
SfO
N.
Fx
z R
AT
IO
OF
R
EM
AIN
ING
V
AR
IAN
CE
IN
N
u D
UE
T
O
VA
RIA
BL
E
EN
TE
RIN
G
RE
GR
ES
SIO
N.
TO
R
ES
IOU
AL
V
AR
IAN
CE
..
N
= N
UM
BE
R
OF
O
BS
ER
VA
TIO
NS
(O
AT
A)
IN
EA
CH
R
EG
IME
.
- 28 -
TABtE !l COEF'F':I,Ct EMTS OF MtH. T.
CQRRELATIO.N AND DETt:RMINAT J<JM~
% DUE TO VARIABLE COEF. OF COEF& OF ENTER I,N Q
REG~fME STEP ENTERiNG MU,T.CORRE"'t - MULT., DETERM, % VAR i A Bh.!
(22) nAN INTRODUCTiON TO HEAT TRANSFER" BY M. FJSHENDEN AND O. A. SAUNDERS OXFORD (1950)
(23) ·PHYSIOAL SIMILARITY AND DIMENSIONAL ANALYSIS" BY Vi.. J. Du N CAN EDWARD ARNOLD & CO., LONDON (1953)
(24) HILSENRATH, J., ET Al., NAT@L. BUR. OF SToS. CIRO. 564, 1955 ~EPORTED IN ~PROPERTIES OF GASES AND LIQUIDS· BY R. C. REID AND T. K. SHERWOOD, MOGRAW HILL, N. Y. (1958)
(25) "STEPWISE PROCEDURE FOR CALOULATiON OF MULTIPLE REGRESSION" BY M. A. EFROYMSON Esso RESEARCH AND ENGJNEERING CO., LiNDEN, N. J. DELIVERED AT GORDON RESEARCH CONFERENCE ON STATISTICS, AUGUST i-12, 1955
(26) "ELEMENTARY STATISTICAL ANALYSiS" BY S. S. WILKS PRINCETON UNIV. PRESS, PRINCETON, N. J. (1948)
(27) "METHODS OF STATISTICAL ANALYSIS· BY C. H. GOULDEN JOHN WilEY & SONS, N. Y. (1939)
(28) IISTATISTICAL METHODS IN RESEAROH AND PRODUCTION" BY O. L. DAVIES, OLIVER ANO BOYD, LONGON (1949)
(29) "METHOOS OF CORRELATION ANALYSIS" BY M. EZEKUEJ.., JOHN WILEY AND SONS, N. Y. (1941)