Open Collections

UBC Theses and Dissertations

UBC Theses Logo

UBC Theses and Dissertations

A laboratory study of thermal convection under a central force field 1971

You don't seem to have a PDF reader installed, try download the pdf

Item Metadata

Download

Media
UBC_1972_A1 C33.pdf
UBC_1972_A1 C33.pdf [ 3.96MB ]
UBC_1972_A1 C33.pdf
Metadata
JSON: 1.0052446.json
JSON-LD: 1.0052446+ld.json
RDF/XML (Pretty): 1.0052446.xml
RDF/JSON: 1.0052446+rdf.json
Turtle: 1.0052446+rdf-turtle.txt
N-Triples: 1.0052446+rdf-ntriples.txt
Citation
1.0052446.ris

Full Text

IIZII A LABORATORY STUDY OF THERMAL CONVECTION UNDER A CENTRAL FORCE FIELD BY BHUVANESH CHANDRA B . S c , U n i v e r s i t y of Lucknow, 1961 M.Sc, Indian I n s t i t u t e of Technology, Kharagpur, 1965 M.Sc, U n i v e r s i t y of Western O n t a r i o , 1969 A THESIS SUBMITTED IN PARTIAL FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY i n the Department of GEOPHYSICS We accept t h i s t h e s i s as.conforming to the re q u i r e d standard THE UNIVERSITY OF BRITISH COLUMBIA November, 19 71 In presenting t h i s thesis i n p a r t i a l f u l f i l m e n t of the requirements for an advanced degree at the University of B r i t i s h Columbia, I agree that the Library s h a l l make i t f r e e l y available for reference and study. I further agree that permission for extensive copying of t h i s thesis fo r scholarly purposes may be granted by the Head of my Department or by h i s representatives. It i s understood that copying or publication of t h i s thesis for f i n a n c i a l gain s h a l l not be allowed without my written permission. Department of Geophysics The University of B r i t i s h Columbia Vancouver 8, Canada Date December 27, 1971. (i) ABSTRACT This t h e s i s presents the r e s u l t s of a t h e o r e t i c a l and experimental study of thermal convection under the i n f l u e n c e of a c e n t r a l f o r c e f i e l d . Flows i n the atmosphere and i n the core of the e a r t h are thought to occur under a near balance between Coriolis and buoyancy f o r c e s . Thus, a d e s i r a b l e model of these flows would i n c l u d e s p h e r i c a l symmetry i n the f o r c e f i e l d and r o t a t i o n . The present study, i n which convection under a c e n t r a l f o r c e f i e l d i n c y l i n d r i c a l geometry has been achieved, i s the f i r s t step towards such a model. The system c o n s i s t s of a c o o l outer c y l i n d e r and a hot i n n e r c y l i n d e r w i t h a d i e l e c t r i c l i q u i d ( s i l i c o n e o i l ) f i l l i n g the annulus between them. The common a x i s of the c y l i n d e r s i s v e r t i c a l . The i n n e r c y l i n d e r i s grounded and the outer c y l i n d e r i s kept at a h i g h a l t e r n a t i n g (60 Hz) p o t e n t i a l . This i n t e n s e a l t e r n a t i n g e l e c t r i c f i e l d provides the r a d i a l buoyancy forc e which r e s u l t s i n convectiva heat t r a n s f e r at a c e r t a i n c r i t i c a l temperature g r a d i e n t . The f l u i d i n the system i s found to behave l i k e a l a y e r of f l u i d i n a g r a v i t a t i o n a l f i e l d , heated from below. Below a c e r t a i n c r i t i c a l v alue of a dimensionless number (equivalent to ( i i ) the Rayleigh number with the e l e c t r i c a l force substituted for gravity) there i s no convective heat transfer. Above the c r i t i c a l value, flow sets in with the convective heat transfer proportional to the modified Rayleigh number. Marginal st a b i l i t y analysis gives a c r i t i c a l e l e c t r i c a l Rayleigh number in agreement with the experimentally determined value. ( i i i ) TABLE OF CONTENTS ABSTRACT (i) LIST OF FIGURES (v) LIST OF TABLES (vii) ACKNOWLEDGEMENTS ( v i i i ) CHAPTER I GENERAL INTRODUCTION 1.1 Introduction .1 1.2 Thermal convection and geophysical applications 4 1.3 Scope of this thesis. 7 CHAPTER II THEORY 2.1 Introduction 9 2.2 The perturbation, equations 11 2.3 -The principle of exchange of s t a b i l i t i e s 23 2.4 The solution for the case of a narrow gap 27 2.5 The solution for a wide gap when the marginal state.is stationary . 28 CHAPTER.Ill EXPERIMENTAL ARRANGEMENT AND RESULTS 3.1 Introduction 37 3.2 Experimental arrangement 40 3.3. .The experimental procedure 48 (iv) 3.4 Errors and corrections 54 3.5 The experimental results 63 3.6 Accuracy 65 CHAPTER IV DISCUSSION AND CONCLUSIONS 4.1 Discussion of results 66 4.2 Conclusions 67 REFERENCES 68 APPENDIX I TABLES OF THERMISTOR CALIBRATION 71 APPENDIX.II TABLES OF HEAT TRANSFER, MEASUREMENTS 74 APPENDIX III TABLES OF CORRECTED DATA 83 (v) LIST OF FIGURES FIGURE 1 Coordinate system and cylindrical cavity. 10 FIGURE 2 Dependence of the numerically determined el e c t r i c a l Rayleigh number on the dimension- less wave number. 32 FIGURE 3 . The velocity profile. 33 FIGURE 4 The temperature profile. 34 FIGURE 5 The stream lines at the onset of ins t a b i l i t y . 35 FIGURE 6 A cross section of the apparatus. 39 FIGURE 7 . The circuit for the heating c o i l . 41 FIGURE 8., Circuit.for the bridge. 44 FIGURE 9 Kinematic viscosity as a function of temperature. 45 FIGURE 10 ... Dielectric constant as. a function of temperature. 46 FIGURE 11 Dependence of the temperature of the outer cylinder less mean of inlet-outlet temperatures on heat transfer rate. 49 FIGURE 12 Calibration of Ammeter 1^ against a laboratory standard. 51 FIGURE 13 Calibration of ammeter I against a laboratory standard. 52 FIGURE 14 Calibration of voltmeter against a laboratory standard. 53 FIGURE 15 Dependence of temperature difference on heat ...... transfer rate for 0 kv rms between the cylinders. 55 FIGURE 16 Dependence of temperature difference on heat transfer rate for 4.06 kv rms between the cylinders. 56 (vi) FIGURE 17 Dependence of temperature difference °n heat transfer rate for 6.00 kv rms between the cylinders. FIGURE 18 Dependence of temperature difference, on heat transfer rate for 6.92 kv rms between tjhe cylinders. FIGURE 19 Dependence of temperature difference heat transfer rate for 7.91 kv rms between fche cylinders. FIGURE 20 Dependence of temperature difference on heat transfer rate for 9.63 kv rms between the cylinders. FIGURE 21 Dependence of temperature difference on hpat transfer rate for 10.15 kv rms between the cylinders. F I G U R E 22 F I G U R E 23 Dependence of temperature differej^pe on heat transfer rate for 10.80 kv rms between t̂ he cylinders. Nusselt number as a function of e l e c t r i c a l Rayleigh number. (vii) LIST OF TABLES Table 1 E l e c t r i c a l body force at the inner and outer cyliners for given voltages across the gap. 22 Table 2 Dimensions of the c e l l . 42 Table 3 Thermistor calibrations for the inner cylinder. 72 Table 4 Thermistor calibrations for the inlet and outlet points. 73 Table 5 Typical properties.of silicone o i l . 48 Table 6 Observed data at 0 kv rms. 75 Table 7 . Observed data at 4.06 kv rms. 76 Table 8 Observed data at.6.00 kv rms. 77 Table 9 Observed data at 6.92 kv rms. 78 Table 10 . Observed data at 7.91 kv rms. 79 Table 11 Observed data.at 9.63.kv rms. 80 Table 12 Observed data at 10.15 kv rms. 81 Table 13. , Observed data at 10.80 kv rms. 82 Table 14. Corrected data for 0 kv rms. 84 Table 15 Corrected data for 4.06 kv rms. 85 Table 16 Corrected data for 6.00 kv rms. 86 Table,17,- Corrected data for 6.92 kv rms. 87 Table 18 Corrected data.for 7.91 kv rms. 88 Table 19 Corrected data for 9-63 kv rms. 89 Table 20 Corrected data for 10.15 kv rms. 90 Table 21 . Corrected data for 10.80 kv rms. 91 ( v i i i ) ACKNOWLEDGEMENTS The author is indebted to Dr. D. E. Smylie for directing the research reported in this thesis. The following persons read the thesis c r i t i c a l l y : Dr. D. E. Smylie, Dr. G. K. C. Clarke and DR. T. J. Ulrych, a l l of the Geophysics Department, and Dr. P. H. LeBlond of the Institute of Oceanography. Mr. H. Lau of the Mechanical Engineering Department contributed to many discussions. Finally, the author wishes to acknowledge the assistance of the staff members and the graduate students of the department who contributed in various capacities. The work was financed by grants from the National Research Council of Canada to Dr. D. E. Smylie. CHAPTER I GENERAL INTRODUCTION 1.1 Introduction In the last few decades, efforts have been made to study motions in the Earth's atmosphere, oceans, mantle and core through the use of laboratory models (Fultz, 1961; Long, 1954). The incentive for the use of .models comes from the d i f f i c u l t i e s of solving the equations of motion for general convective flow problems. At the present time, the equations governing the general flows can only be solved approximately under idealized conditions. Thus, laboratory models have become, important in geophysical f l u i d dynamics. The large scale circulation of the atmosphere occurs under the combined influence of a spherically symmetric gravitational f i e l d and the rotation ..of the Earth. .It i s suspected-that a. roughly similar situation prevails. in the Earth!s liquid core.... Thus, a laboratory model of these flows.should .include both.the spherical symmetry of the gravitational f i e l d and the rotation of the Earth. The effect of .rotation has been-the subj.ect..of study, of many, authors (for recent references, see. Greenspan,. 1968). . In these ..studies , certain flow regimes have been recognized which seem to have their counterpart in the atmosphere in the form of large scale eddies or cyclones. However, no .laboratory model which simulates a spherically symmetric 2 g r a v i t a t i o n a l f i e l d has yet been developedo In the absence of a s u i t a b l e method of s i m u l a t i n g a g r a v i t a t i n g f l u i d sphere i n the l a b o r a t o r y , the v a l i d i t y of the r o t a t i n g models i s l i m i t e d to the. p o l a r regions where strong C o r i o l i s f o r c e s are present as i n the core and the atmosphere„ In the t h i n s h e l l of the atmosphere, however, t h i s defect i s not so ser i o u s because the flow i s e s s e n t i a l l y two dimensional as opposed to the core flow which i s three ,dimensional. This t h e s i s i s the outcome of a p r e l i m i n a r y i n v e s t i g a t i o n f o r s e t t i n g up a la b o r a t o r y model of g r a v i t a t i o n a l convection- Smylie(1966) f i r s t suggested the use of strong e l e c t r i c f i e l d gradients i n d i e l e c t r i c l i q u i d s to simulate the c e n t r a l nature of the g r a v i t a t i o n a l f i e l d . The e l e c t r i c a l p r o p e r t i e s of a . l i q u i d are fu n c t i o n s of temperature„ For a nonpolar l i q u i d the v a r i a t i o n of d i e l e c t r i c , constant w i t h temperature i s due s o l e l y to the change i n density (Debye, 1929, p 27)„ Thus, i n the presence of temperature v a r i a t i o n s , i n the d i e l e c t r i c l i q u i d , . a n intense e l e c t r i c f i e l d r e s u l t s i n a buoyancy force„ This buoyancy force w i l l produce an i n s t a b i l i t y i n the f l u i d o However,.the v i s c o s i t y , t h e r m a l , c o n d u c t i v i t y and the boundaries w i l l act to s t a b i l i z e the f l u i d and s i g n i f i c a n t convective heat t r a n s f e r should occur only when the temperature gradient i s appreciable,, In d i e l e c t r i c l i q u i d s of moderate d i e l e c t r i c constant, the thermal convection produced by the e l e c t r i c buoyancy f o r c e i s analogous to g r a v i t a t i o n a l convection, where the i n s t a b i l i t y occurs as a r e s u l t of heat i n g from below (or c o o l i n g from above). In con t r a s t to the ubiquitous 3 and unalterable gravitational f i e l d , however, electric fields can easily be generated and shaped in the laboratory by existing techniqueso In a steady electric f i e l d space charges accumulate even in a good dielectric and the circulation of the liquid i s due to the movement of these charges. The free charge buildup occurs exponentially in time with a time constant e/a where e is the permittivity and a is the conductivity of the fl u i d . This constant is known as the el e c t r i c a l relaxation time. If an alternating f i e l d is applied at a frequency much higher than the reciprocal of the relaxation time, free charge does.not have time to accumulate. Thus, in an alternating f i e l d , the dielectric effects dominate the free charge effectso The electric body force depends on the gradient of the square of the electric f i e l d * For most.dielectric liquids (for example, transformer o i l and silicone o i l ) , 60 Hz is a sufficiently high frequency to prevent the buildup of free charge. At the same time, dielectric loss at 60 Hz is so low that i t makes no significant contribution to the temperature f i e l d . Further, variations in the body force are so rapid that i t s mean value can be assumed in determining f l u i d motions, except in the case of liquids of extremely.low viscosity. . Both theoretical and experimental study of electroconvection (convection under steady fields) has been reported by many authors (for example, Kronig and Schwarz, 1949; Turnbull, 1968; Malkus and Veronis,.1961; Melcher and Taylor, 1969; etc.).• The electroconvective 4 instability occurs even in the absence of temperature variations; Avsec and Luntz (1937) have observed two dimensional toroidal motions in a dielectric liquid f i l l i n g the gap between two concentric cylinders. When the inner cylinder i s earthed and the el e c t r i c a l potential of the outer cylinder i s raised, then at some c r i t i c a l potential difference, steady cellular patterns are formed in the flu i d . Recently, Gross(1967) and Gross and Porter(1966) have suggested the.use of space charge effects in modelling geophysical flows. In these models, charge transport becomes the analogue of heat transfer and. i f .temperature gradients are present, heat w i l l be advected by the .charge induced flow. The mechanisms of charge generation and transfer at the electrodes, and the nature of the process by which charge motipn s t i r s the f l u i d are not well under- stood. Therefore, the quality of analogy to geophysical phenomena is not clear. 1.2 Thermal convection and Geophysical applications; Although the phenomena of thermal convection as a mode of heat transfer was discovered in the 18th century, the f i r s t quanti- tative experiments on convection were done by Be"nard(1900). He found that i f a thin layer of liquid, free at i t s upper surface, is heated uniformly at the lower surface, a regime of polygonal convection cells i s formed as soon as a certain c r i t i c a l temperature gradient i s reached...The walls of these cells are vert i c a l and the movement of the liquid i s upward in the centre and downward at the periphery. 5 Rayleigh(1916) lai d the theoretical foundation for the explanation of Benard's results. It appears now, however, that the convection in Benard's experiments was actually partially driven by the variations in the free surface conditions and not entirely by the buoyancy forces as Rayleigh assumed. In later experiments (see for example Chandra, .1938; .Schmitd and Milverton, 1935; and Silveston, 1958), the free surface has been eliminated by placing an isothermal l i d on top of the f l u i d layer. These experiments also showed the development of regular hexagonal c e l l patterns. The c r i t i c a l temperature gradient at which .convection cells develop i increases rapidly as the layer becomes thinner.. Several modifications of Rayleigh's theory (for example, Jeffreys, 1926; Pellew and Southwell, 1940; Low, 1929; etc.) to f i t this case give excellent agreement with the experimental observations. Rayleigh's approach to the problem (infinitesimal, perturbations on an equilibrium state) s t i l l remains the basis of the modern treatment of the marginal s t a b i l i t y problems (Chandrasekhar, 1961, chapter 2). It i s believed that many flows of geophysical interest are driven by thermal convection. It has been suggested that large scale convection currents exist in the garth's mantle in attempts by many authors to interpret certain topographic features of the Earth's surface (Vening-Meinesz, 1962). Recently, McFadden(1969) has considered the effect of a region of low viscosity on thermal convection in the Earth's mantle to try to explain how convection 6 in only the upper mantle could lead to flows of continental horizontal scale. The presence of thermal convection in the atmosphere was suggested by Hadley (1735); differential heating between different latitudes gives rise to ascent of air in tropical regions. This leads to a flow of air towards the equator at lower levels and a flow away from the equator at higher levels. The combined effects of surface f r i c t i o n and the rotation of the earth deflects these currents and thus, according to Hadley's original theory, easterly winds form at lower levels and westerlies at higher levels. From this idea, the well known cellular model for the circulation of the whole atmosphere has been derived (Rossby, 1941). The geomagnetic f i e l d is now considered to originate from electric currents produced by inductive interaction between hydrodynamical motions in the Earth's conductive f l u i d core and a small adventitious magnetic f i e l d . Thus, a satisfactory theory for the Earth's magnetism requires a satisfactory theory for the hydrodynamics of the core. Several speculations have been made regarding the energy source that drives the f l u i d motions in the core. Among them thermal convection i s one of the pos s i b i l i t i e s (Bullard and Gellman, 1954). Oceanographers have long studied the circulation of the oceanic currents. The ocean currents are believed to be the result of the combined effects of the thermohaline motions (density inhomo- geneity caused by temperature and/or salini t y differences) and the wind driven motions. The former are thought to be more important 7 in deep water. However, unlike the case of the atmosphere, the thermal convection in the oceans occurs, not because of heating from below, but because of cooling from the upper surface. The thermohaline circulation originates as a vertical flow sinking to mid-depth or even to the ocean bottom, followed by horizontal motion. 1.3 Scope of the present work This thesis is intended to demonstrate the f e a s i b i l i t y of using an intense alternating electric f i e l d acting on a dielectric liquid to produce a spherically symmetric force f i e l d . A buoyancy force is produced by the variation of dielectric constant with temperature. Since i t is a f e a s i b i l i t y study, the experiment was set up in cylindrical.geometry to.simplify .construction. However, the techniques, learned here can be carried directly over to the spherical-case, which i s the real geometry for large scale geophysical flows. . . The system consists of two vertical concentric cylinders of r a d i i r^ and r^ ( r^ > rj_) • The .space between the cylinders is f i l l e d with a dielectric liquid (silicone o i l ) . The inner cylinder is maintained at a fixed temperature T^ which i s higher than the constant temperature T2 of the outer cylinder. In addition to a temperature gradient, a radial., alternating electric f i e l d i s also imposed be tween . the cylin de rs. In chapter II, the marginal st a b i l i t y equations for this 8 problem are derived under the assumptions of the Boussinesq approximation. The effect of the Earth's gravity has been ignored for this theory. The solution of the marginal st a b i l i t y equations has been found numerically and i t is shown that the convection induced by the e l e c t r i c a l body.force occurs only when the dimensionless number (equivalent to the Rayleigh number with .the electric buoyancy force substituted for gravity), exceeds a.certain c r i t i c a l value. When the gap thickness i s very much smaller than the mean radius of the cylinders (in this limiting case boundaries take the form of plane surfaces, rather than cylindrical surfaces),,the marginal sta b i l i t y equations are shown to reduce to the case of the Be"nard- Rayleigh in s t a b i l i t y . In.chapter III, the experimental arrangement, technique and results of elect r i c a l l y induced convection in a vertical annulus are described. The onset of thermal convection was detected by heat transfer and temperature measurements. The f i n a l experimental results, are shown,in the form of a graph of the Nusselt number against the electric Rayleigh number. The.study is summarized in the f i n a l chapter. 9 CHAPTER II THEORY 2.1 INTRODUCTION The problem to be studied is the thermal convection occuring under the el e c t r i c a l body force in the enclosed annular region formed by two vertical concentric cylinders of radii .r 1 and ( r ^ r ^ ) . This geometry i s shown in figure 1. Conventional cylindrical coordinates (r,6,z) are indicated with z vertically upwards. The inner wall of the cylindrical cavity is held at a fixed temperature T.̂  which is. greater than the constant temperature T^ of the outer wall. In.addition to. a temperature difference, there i s also an alternating difference of e l e c t r i c a l potential between the cylinders. In the absence.of any temperature variations, an electric body force is.produced in.the liquid,.and in an equilibrium situation, this body force i s entirely, balanced by the generation of a pressure gradient (Landau and Lifshitz,, .1960, pp 64-69). The presence of temperature variations in the liquid results in buoyancy forces whereby warmer liquid has,a tendency to- seek regions of less intense electric f i e l d , cooler, liquid has a.tendency to.seek regions.of more intense electric f i e l d . However,-this natural tendency on the part of the f l u i d w i l l be inhibited by i t s own viscosity. Thus, i f the temperature difference F I G U R E 1. C O - O R D I N A T E S Y S T E M A N D T H E C A V I T Y . 11 between the vertical bounding walls is sufficiently small, the heat supplied at the inner wall i s transferred to the outer wall by conduc- tion alone. For greater temperature differences, the electric buoyancy force i s sufficient to overcome viscous dissipation and the f l u i d assumes a convective motion. In other words, we expect that the temperature difference between the cylinders must exceed a certain minimum value before convective heat transfer is realized. The problem under consideration, then, i s to solve the marginal s t a b i l i t y equations for an incompressible f l u i d layer bounded by vertical cylindrical walls at different temperatures, in the presence of a radial e l e c t r i c a l buoyancy force. 2.2 The Perturbation Equations The following simplifying assumptions are made in deriving the perturbation equations: (1) The f l u i d i s incompressible and the density constant, except as i t modifies the electric-buoyancy force term. (2) The mean square electric f i e l d and the temperature distribution in the annulus are functions of radius only. (3) The cylinders are of i n f i n i t e extent, so that the end effects can be ignored. (4) The motion i s slow and the components of velocity are small enough so that in a f i r s t approximation, their products and squares can be neglected. 12 The kinematic viscosity and the ,.dif fusivity of the f l u i d are true constants. The various perturbations are axi^symmetric and thus independent of 0. In deriving the. marginal s t a b i l i t y .equations, the effect of the earth's gravitational f i e l d can be ignored. Strictly speaking i t is .not so.. Batchelor. (19.54). showed, that-in .a.narrow cavity between vertical, boundaries at different ..temperatures, there is a- slow ve r t i c a l flow but the heat transfer across the cavity is due mostly to conduction. Thus, the present problem is unstable for a l l , temperature gradients since the earth's gravity is. always present.. The. maximum Grashof number involved .in the experiment (described, in. .the next chapter* was about 7.5. This i s roughly 3 .10, .times, less than the ..Grashof. number....at which free convection sets -in.. (Vest, and Arpaci, . 1969)...... The gravity induced, flow i s likely to be very small., at. the.. Grashof ..number involved in the experiment. It would, be an important generalization .to carry through a solution taking into account the interaction of the convection, induced by the ele c t r i c a l forces with the vertical gravitational base flow,,, but this d i f f i c u l t problem has not been tackled here. Consider that i n i t i a l l y there are no motions. Thus, the i n i t i a l 13 state i s characterized by U = 0 T = T(r) (1) where U is the velocity vector and I is the temperature. In'the absence of any .motions, the hydrostatic equation reads -grad p + 1 = 0 (2) o o -»• where f -is . the electric body, force per unit volume at density p and temperature T : p is the pressure distribution, o o o r , . The .temperature distribution...is governed by the equation . <» The solution of equation (3) appropriate to the boundary conditions gives T T T dT Q 1 2 1 dr r„ r (4) i n / i . "1 The solution for the elec t r i c f i e l d distribution E is o Ej = (E q,0,0), (5) . o r 2 r (6) where V q is the root mean square value of the applied potential difference. 14 Let the perturbed state be characterized by the velocity components V u e, V ( 7 ) body force f = + t\ (8) density p = p Q + p , (9) electric f i e l d E = E + E , (10) o i temperature T = T Q + T , (11) i and pressure p = P Q + P , primes indicating flow induced quantities. The. general equation for the body force resulting from electric fields in f l u i d dielectrics .is developed by Landau and Lifshitz (1960,,p.68) - They showed that in an uncharged f l u i d dielectric of uniform .composition,.,the .electric body force, per unit volume is given by V= 1 / 2 p o & r a d [ Eo (-If }T ] " 1 / 2 E d ( If )p ^ a d T o , (12) o o where e i s the electric permittivity of the f l u i d . In the perturbed state at temperature T and density p , the body force.per unit A o l u m e is 1 = 1/2 pgrad [E 2 ( -g- \^ - 1/2 E 2 ( | | ) pgrad T (13) 15 The permittivity of a substance is a function of density and temperature. If the perturbations are small, the permittivity at temperature T and density p can be expressed in terms of permittivity at the reference temperature T q and density p Q by writing a Taylor expansion. Thus, + 1 / 2 1 $ v ' 2 + 2 c \-*.p,r + c d \ T ' 2 ] 3p o o o 3 i o + higher order terms " (14) Therefore, we can write , 3 e ' ( 3 £ . . d2c . , . 3 2 e . , ^ 3p ;T S 1 3p ;T - 2 T P ^ 8p8T ; p ,T o 3p . o 0 0 and 3 1 p 3 1 po 3 T 2 Po 3p3T po' To Combining (8), (12), (13) and (15), we get to f i r s t order i n the flow-induced quantities, f* - l / 2 p'grad[E o 2 ( ) ] - 1 / 2 E Q 2 (|| ) gradT' o o + p ograd [ ( 2 o . ^ ) ( f£- ) T ] - ( E ; . I ' ) ( | | > gradT o o o . + l / 2 p o g r a d [ E o 2 { ( ^ f ) T p ' + ( ) T'} ] 3p o o o - 1 / 2 E / I (4 >p *' + ( f e ) 0 » ' ] « » « . <») 3T o o o 16 Chandra (1969) showed that for a f l u i d dielectric of moderate dielectric constant, the flow-induced changes in the electric f i e l d can be neglected. The order of magnitude of the second density and temperature coefficients of permittivity for fluids i s very small and can also be taken to be zero. . Thusthe. ele c t r i c buoyancy force per unit volume can be written as . f - 1/2 p' ( f ^ - > T g r a d E Q 2 - - i ° - ( | | ) p grad' T' (17) with p' - - ctpT'5 (18) a denoting the coefficient of thermal expansion. The equation (17) i s quite general and can be applied to a f l u i d dielectric .of any shape.. . Using equation .(6), the e l e c t r i c buoyancy force per unit mass for the cylindrical case under consider- ation can be written as v 2 r 2 Sp 'T r J » r [ l n - ^ ] r l (19) where r is the unit radius vector. Having found.an expression for the e l e c t r i c buoyancy force per. unit mass, the linearized perturbation.equations:can be written as follows: (i) momentum 17 a t " " v [ v ue " - I 1 <21) r 2 Ci T  V z - 9 / £l x 4 . o2TT i /. o 1 , 3e . 1 3T* ° ( I n — r r r l ( i i ) continuity 3U U 3U dr r dz ( i i i ) energy T r T 2 \ _ 3t r r ln — r l K V * T ' (24) 2 where V has the meaning 2 2 V 2 = - — + - — + - — (25) a 2 + r 3r + a 2 K J 3r 3z and v and K are the kinematic viscosity and the thermal diffusivity of the f l u i d . By analysing the disturbance into normal modes, the solutions of equations (20)-(24) are to be found in the form U r = e p tU(r) coskz, T 1 = e p t0(r) coskz, U = e p tV(r) coskz, p* = e p t f i ( r ) coskz, V(26) Po U = e p tW(r) coskz. z where k i s the wave number of the disturbance in the axial direction 18 and p Is a constant which may be complex. For solutions of the form (26), equations (20)-(25) reduce to V 2 v(DD,-k2-P/v)U + — 2 - - [a( f ), % - -f- ( j§ ) \ jf ] - DO, (27) ( I n - 2 - ) 2 r l (DDA - k 2 - p/v)V = 0, (28) V 2 v(D*D-k2-p/v)W + -fr °— ( | | ) *f = -kfl, (29) 2p r 3T p r2 (In — ) r l D̂ U = -kW, (30) T -T 2 , v„ 1 2 U (D^D-k"-p/K)0 - - - ^ - f ^ , (31) i 2 i c l n — r l and V 2 = D^D-k2, (32) where D = —- (33) dr a n d D* - h + r' ( 3 4 ) For the case when p=0 the solution of equation (28) is V = AI Q(kr) + BK Q(kr) (35) where I and K are modified Bessel functions of zeroth order, o o Substituting the boundary conditions that on r i g i d boundaries V=0, we get A = B = 0. Hence, V = 0. (36) This shows that for the form of assumed solutions, there cannot 19 be any velocity perturbation along 9 , Thus, the basic flow i s in the form of ring vortices. Eliminating W between equations (29) and (30), we get - | (DsI)-k2-p/v)D,,U + - ^ T C | f ) p i f . - k f i (37) Substituting the above expression for o, in equation (27), we get after some.simplifying steps V 2 2 (DD,-k 2-p/v)(DD,- k 2)U- - 2 - Io(.|£ >T - ± ( || > ] % (ln -*x2 r For fluids of uniform composition, the permittivity may be regarded as a function of ..the state variables alone and we may write M'.- 'I't-M'p (39) Thus, equation -(38). reduces to 2 (DD^-p/vXDD.-k 2)^ - - ~ \ £ < | | ) ^ (40) ( I n / ) 2 1 Equation (40) must be solved together with T -T (D.D-k2-p/<)0 = - —±—r ' 7 (41) r2 • <ln — r l Equations (40) and (41) are the marginal st a b i l i t y equations, the solution of which must be found under the boundary conditions • U =..DU = 9 » 0 (42) at p r ^ and r 8 * ^ . 20 Dimensionless parameters to describe solutions i n general form may be obtained by non-dimensionalizing the equations,in a suitable manner. A convenient set of non-dimensional variables i s r - r. *\ J2£ > (43) a =» kd A » £ where the cavity gap - is written as d. The equations (40) through (42) become (DD #-a 2-o)(DD #-a 2)U 2 2 a d __1 , 3e > e I l n f 2 ] 2 r x 3 ^ 9 T P (1+AO3 (44) (D.D-a -Pc)G -•- T r T 2 d 2 •cin r 2 r x (1+X?) ' • (45) U = DU •» 9 0 at C B 0 and 1. (46) D and now have the following meaning de: • (47) (48) and P (» — ) I s the Prandtl number* It is convenient to make the transformation KT r I n — 9 +6, ( T j - T ^ d 2 r l where 0 now has the dimensions of velocity. The marginal s t a b i l i t y equations become (DD.-a2-o)(DD.-a2)U= -Ra2 9 (1+U) and (D*D-a2-Pa)0 = ^ R may be interpreted as the e l e c t r i c a l Rayleigh number ag Bd 4 R = VK where y 2 o 1 , je_ . *e " " 3 n r2,2 p a 8 T P ^ ( I n — ) can be regarded.as an el e c t r i c a l l y derived .gravity at the surface the inner cylinder and T -T 1 2 6.- r 2 r l l n ^ is the temperature gradient at the inner cylinder surface. Clearly, the e l e c t r i c a l gravity, g^, is a strong function of radius and varies inversely as the cube of.radius. 22 Table 1 The e l e c t r i c a l body force at the inner and outer cylinders for a given voltage across the gap. r, - 1.711 x 10_2m . .e =8.854 x 10~ 1 2 farads/m 1 o r 2 • 1.903 x 10~2m p = 937.7 kg/m3 _3 a • 1.08 x 10 cc/cc/°c 3 , e . „ - , 0 - , n ~ 3 / o ( — ) = 3.72 x 10 / c o VOLTAGE ACROSS Ele c t r i c a l gravity... . E l e c t r i c a l gravity CYLINDERS ... Earth's gravity '  y Earth's gravity kv(rms) at r=r^ at r=r 2 0.2343 0.1703 4 0.9372 0.6812 6 2.1087 1.5327 8 3.7488 2.7248 10 5.8575 4.2574 12 8.4348 6.1307 In Table 1, the values of e l e c t r i c a l gravity compared to earth's gravity are shown both at the surface of the inner cylinder and that of outer cylinder. At 8.0 kv rias voltage difference between.the cylinders, the average elec t r i c gravity i s three times that of the earth's gravity. Above this voltage difference, i t rises even more sharply. .The form of the equations. (50) and (51) and the boundary conditions (46) shows that, in the marginal state (i.e. when a=0) , for a given value 23 of a, the dimensionless parameters whose values are sufficient to determine uniquely the distribution of U and G are the el e c t r i c a l Rayleigh number R and the gap to inner radius ratio X. 2.3 The Exchange .of Stabilities In general o" is a complex quantity and is a function of the physical quantities involved and of the parameters characterizing the particular pattern, of the disturbance. If Rl(a) is positive, the disturbance increases, but i f i t . i s .negative, the disturbance dies away. In a.set of possible disturbances, the mode for which Rl(o)is. a-maximum wil l , be the f i r s t to appear beyond st a b i l i t y , when the Rl(o) i s just equal to zero, the limiting condition of st a b i l i t y will.be realized. I t . i s important to know.if the Im(a) is zero when the Fl ( a ) i s . v If .Im(a) i s not zero, the disturbance manifests i t s e l f .in the form. of. a wave.motion. and the system is said to be overstable. -When Im(o) is zero i f the Rl(a) is zero, we say that the principle of exchange of . s t a b i l i t i e s is valid. Since the vanishing of .cr means that a l l time variations disappear, this limiting condition represents a steady state in which the disturbance just maintains i t s e l f . To consider, the principle of exchange of s t a b i l i t i e s for the problem i n hand,.it i s . s l i g h t l y .more convenient to write the non-dimensionalized form of equations (40)-(42) with respect to 24 r^, the radius of the outer cylinder. Then, we have (DD A-a 1 2-o 1)(DD^-a 1 2)U = ^ 3 (55) r < D * D - a l 2 - P 0 l ) Q = 7 (56) and U = DU = 6 » 0 at r = n and 1. (57) where D - ^ , D* - ̂  + ± , (58) ar a l = k r2 " -T • ( 5 9 ) a R, " ~r and n - T T 1 . (61) XH 2 Multiplying equation (55) by rU* (the star superscripts represent complex conjugates) and integrating in the interval r = n and r =1, we get f rU*(DD A-a 1 2-a 1)(DD A-a 1 2)U dr = - R ^ [ ^ U* dr (62) in r Since U and i t s derivative vanish at the boundaries, the l e f t hand side of equation (62) can be shown to be positive definite (Chandrasekhar, 1961, p297). Thus j rU*(DD A-a 1 2-a 1)(DD^-a 1 2)U dr = f1 r|(DD^- a i 2 )u| 2dr + o r f [ r | g | 2 + (^ + a ^ r ) ! ^ 2 ] dr (63) 'n Jn 25 Substituting for U* from equation (56.) in the right hand side of equation (62) , we get £ ^ U* d r = l n f ( n *D-a 1 2 -Pcv*)0*dr (64) = f 1 - DJ)0*dr - (a 2+Po*) t L dr (65) Jn r * 1 I n r Again making use of the boundary conditions on 0 and integrating by parts, we.can easily show that J 1 | D^D0*dr = - ^ |D0| 2dr + %:> D9*dr (66) Combining equations (62), (63), (65) and (66), we get + I 2 = R 1a 1 2t(a 1 2+Po*)I 3 + 1^ + Ig] (67) where Ix = J 1 [r|-BU|2 + ( i + a 2 r ) | u | 2 ] dr (68) I 2 j^.. r|(DD^-a2)u|2dr (69) 26 12~t) ^ l e f d r (70) [ „ i | D 6 | 2 d r (71) and I 5 = -2 f 1 0_ D9*dr (72) Jn 2 1 r The integrals 1^, , I^ and 1^ are positive definite while !<. is complex. Equating the imaginary parts of equation (6 7), we obtain Im(a) [ I 1 + R i a ; L 2PI 3] = R 1a 1 2Im(I 5) (73) and no general conclusions regarding the st a b i l i t y of flow can be drawn from this equation. The validity of the exchange of s t a b i l i t i e s cannot be rigorously ;proved for this problem. In .the next.section, i t is shown that, for a narrow gap, the s t a b i l i t y equations reduce to the form of the. Rayleigh-Benard case.. The principle of the exchange of s t a b i l i t i e s for the Rayleigh- Benard form of equations has been proved elsewhere (Chandrasekhar, chapter II).. .In.chapter III, the experiments on electric body force, induced convection are described and they show that the i n s t a b i l i t y appears t o s e t in as a s t a t i o n e r y c o a v a c t i v a flew. Without attempting a.theoretical ju s t i f i c a t i o n of the principle of exchange of s t a b i l i t i e s , we shall assume i t s validity to obtain the numerical solution, of the eigen-value problem. 27 2.4 The equations.for the case of a narrow gap. If the gap r 0 - r (~d) between the cylinders i s small compared 2, 1 r "t°r to their mean radius 1 2(=r ) ,D. can be simply replaced by D, — 2 — o * (Chandrasekhar,1961, p 402). In this scheme of approximation, equations (40) and (41), become „ 2 9 2 9 9 V (D -a -a)(D -a Z)U = - \ ° a d ( 3e . © _ (74) r ~ ~ "vp • 1 8T ;p r 3 [ln — j 2 o V 2 and (D2-a2-Pa)G =- < T i ~ T 2 ) d JJ ' (75) r„ r , 2. O Kin — • r l where 5 - H!i . k = f , P/v = ^ (76) d d The boundary conditions (42) are U = DU =. 0 = 0 at £=0 and 1 (77) By making the transformation (T-T.)d 2 - • 1 • U+ U (78) icr In — ° r l equations (74) and (75) become (D 2-a 2-o)(D 2-a 2)U = -R0 (79) (D2-a2-Po)0 - U (80) where , . . K - — (81) VK 28 g g and g are now defined at the mean radius r^ and have the following meaning ^ K - - - r - ^ - 2 1 ^<H>, <82) T -T and g = 1__2_ ( g 3 ) i r 2 r In — r l Equations (79) and (80) are the same as for Rayleigh-Benard ins t a b i l i t y (Chandrasekhar, chapter 2 ) „ The principle of exchange of s t a b i l i t i e s for this case holds and the instability occurs at a Rayleigh number 1707., 762 corresponding to a = 3 o l l 7 „ This analysis serves to show that the principle of exchange of s t a b i l i t i e s i s . v a l i d in the small gap approximation. However, i t ignores the curvature effect of the cylinders and cannot be used for any comparision with the experimental results. 2o5 The solution for a wide gap when the, marginal state i s stationary, In the marginal state (a=0), the equations (55) and (56) become _ „ o n (DD* - a ^ r u = -R.^^ (84) r and (D*D - a±2)0 = ^ (85) On eliminating U between equations.(84) and (85), we obtain QVI + . 7 QV + 5_ _ 3 2 IV _ ^ + ]W_ III r 2 1 3 r r r . 2 _ 2 4 + ( -j - — + 3a ) 0 - ( - j - — — ) 0 r r r r - a . 6 6 - - R ^ 2 ^- (86) 1 1 1 4 r 29 with the appropriate boundary conditions 2 II J I 2 III . 2 II 2 1 a l ft 0 = 0 - &1Q = Q + -p - a x 0 - 0 = 0 (87) r l at r 1 — and 1 0 r2 Alternately, we can eliminate 0 between equations (84) and (85). The eigenvalue problem, then., becomes 2 2 T T V I - 4 - 1 ^ . 4 . I 1 6 * o 2 > U T I V ± f2 ^ l W.IH ,3 , 2 3 S 1 - 4...II U + -U + (—-pp •» 3a )U + (^ — ) U - ( — + r — - 3a. )U r r r r . + ! ^ . ) D i _ (3_ _ „. Z!i_ + 6 ) 0 . 2 u_ ( 8 8 ) r r r r r r subject to boundary conditions U.= U 1 , U I V + | u m - (\ + 2 a 1 2 ) U 1 1 + (^ - ̂ U 1 - ( \ - ̂  - a ^ U 5 0 r r r r at r - — and 1. r2 The c r i t i c a l values of and a^ were found by solving numerically equation (86) with boundary conditions (87). The. general solution of ..equation (86) is 0 = A 1 © 1 + A 2 0 2 + A 3 © 3 + A 4 0 4 + A $ © 5 +A 60 6, where A.^.. .Ag are constants. The condition that these solutions 0 ,0,.,,8, form a 1 z o fundamental set is that their Wronskian is not zero. There is an (89) 30 i n f i n i t e number of possible sets of solutions. One simple set pai* be found by choosing 0^(r) such that 0 1 ( r I / r 2 ) = 1, ©.^(^/r, ,) = ©^Vi/O = ....-6^ {T^TJ = 0 (90) and defining 0 (r) where n = 2,3.. <6 as that solution which satisfies n the i n i t i a l conditions 0 ^ ( r ) « 0 I f m 4 n, ra •-• 1,2...6 (91) n 0 m ~ 1 ( r ) = 1 i f m = m m = 1,2...6 (92) n Then 0,(r), 0_(r).... 0, (r) form a fundamental set and the value 1 2 6 of their Wronskiah at r = — i s unity (Ihce,. 1956jp 119). ! r2 For a. given, value of a^ and. R̂. the solutions 0^ .(£),, 0^(r) ... 0^(r) , with the i n i t i a l conditions (90)-(92), were- found numerically bf Xhe method, of Runge-Kutta by breaking the sixth .order differential equation (86) into six simultaneous f i r s t Order equations' (Ford, 1955, chapter 6) by writing <z>'1 =• G 0 1 1 - H e m = J 0 I V - K 0* = L and L 1 + ^ - L + ( ^ - 3a 2)K - ( ^ + ^ * ) J + ( ^ - + 3a^)U 2  r 4 r r r o 5a1 7a 2 -(^5 - - 3 " - )G + ( R / i _ ¥ a . 6 ) 0 = 0 (93) r r r 4 ' 1 r 31 Knowing 0,(.r), Qn(r) ... 0, (r) and thus 0, the boundary 1 2 o conditions (87) were substituted to give a set of six linear homogeneous equations (three at the inner boundary and three at the outer boundary)., The condition that, they have a non-trivial solution is that the determinant D E T = a l l a12 a, a13 a14 a15 a16 21 l31 41 51 a61 a62 a63 %6 - f(a 1,R 1) i 0 (94) where a,, (j; = I...6) aire the coefficients- of A. in the f i r s t 1.1 J' 3 equation; a^j are coefficients o f : i n .the..second,.equation and so on. For a given.value pf-a^, some .value of w i l l cause the determinant to vanish. .To find that value, determinant ( 9 4 ) was calculated for .increasing values pf R̂  t i l l the determinant changed sign..... Knowing ..positive^ ..and- ,.negat;iy& -values, 'of the determinant, corresponding to two values pf %.^*. & linear interpolation scheme (regula f a l s i , method, Mathews and Walker, 1964, p339) was used 3.08 3.10 3.12 3.14 DIMENSIONLESS WAVE NUMBER a FIGURE 2. DEPENDENCE OF THE NUMERICALLY DETERMINED ELECTRICAL RAYLEIGH NUMBER ON Ti i'L DI KENS! G,'<LL".r>:. WAVE NUMBER. 33 F I G U R E 3- T H E V E L O C I T Y P R O F I L E A T T H E O N S E T O F C O N V E C T I O N . T H E V E L O C I T Y I S N O R M A L I Z E D T O U N I T M A X I M U M A M P L I T U D E . 34 F I G U R E 4. T H E T E M P E R A T U R E P R O F I L E A T T H E O N S E T O F C O N V E C T I O N . T H E T E M P E R A T U R E I S N O R M A L I Z E D T O U N I T M A X I M U M A M P L I T U D E . 0.0 0.2 0.4 0.6 o.e 1.0 FIGURE 5. THE STREAMLINES AT THE ONSET OF INSTABILITY. THE STREAM FUNCTION ^ (PROPORTIONAL TO rU cos az) HAS BEEN NORMALIZED TO UNITY AND THE CELL PATTERN IS DRAWN SYMMETRICALLY ABOUT z=0. THE UNIT OF LENGTH IS THE THICKNESS OF THE ANNULUS. 36 to iterate and find for which the determinant vanished. R̂ , in general, converged to a f i n i t e value (to eight significant figures) within six iterations. For these values of a^ and R̂ , the value -12 of the determinant was less than 10 and was considered adequate. Knowing a^ and R̂ , the corresponding values of a and R were found by making use of equations (59) and (61). The results of these computations are shown in figure 2. The solution of interest i s the minimum value of R for which f(a,R) = 0. Quadratic iteration was used here to find the value of a for which R is minimum. The method has been described by McFadden(1969). The c r i t i c a l values of R̂  and found numerically are R = 2119.346 ;. a = 3.117 c c using the values r^ = 1.711 cms.; = 1.903 cms. for the radii of the cylinders. These values were also calculated by solving the sixth- order equation in U and were found to agree within three decimal places (seven significant figures). Figure 3 shows the profile of the velocity distribution corresponding to L and a^. Figure 4 shows the temperature distribution. The stream lines at the onset of instability are shown in figure 5. 37 CHAPTER III EXPERIMENTAL ARRANGEMENT AND RESULTS 3.1 Introduction Chandra(1969) has described the cylindrical experimental arrangement in detail. He also presented some preliminary results oh e l e c t r i c a l body force driven convection. Although the heat transfer measurements clearly indicated the change from conduction to convection regime, the measured Nusselt number (ratio of total heat transfer to heat transfer by conduction alone) In the conduction regime was found to be 10-15% less than i t s value of unity. No definite explanation for this discrepancy could be given. However, i t was suspected that the discrepancy arose partly due to the heat losses from the ends of the heating c o i l and partly due to the errors in the inference of the temperature of the outer cylinder. The outer cylinder was cooled by transformer o i l kept at a constant temperature by circulating i t through an external thermostatic bath. The temperatures at the i n l e t and outlet points (figure 6, 11 and 12) to the thermal flask were measured and their mean was taken to be the temperature of the outer cylinder. However, for a short contact time of flowing li q u i d at a temperature lower than that of the so l i d wall, the f l u i d temperature changes appreciably only i n the immediate Vicinity of the wall (Bird, Stewart and Lightfoot, 1966, p349). Away 38 from the. wall, there i s hardly any change in fl u i d temperature even when there i s increased heat transfer between the cylinders. It was therefore, f e l t necessary to do the experiments again to be able to compare experimental results with theoretical values. In.these experiments, the heating c o i l was placed much closer to the wall of the inner cylinder to minimize heat loss from the ends and to assure radial heat transfer. However, heat transfer and temperature measurements showed that, this could not account for the observed discrepancy. Next, temperature of the outer cylinder was measured at zero applied potential independent of. the inlet and outlet temperatures. The true temperature of the outer cylinder was found to be different from the. mean. A.linear relation was obtained between the mean of the inlet and outlet temperatures and the true temperature of the outer cylinder- This was done only when there was no potential difference across the cavity. When the outer cylinder was at a high voltage,.Its temperature could not be.measured directly. However, from the linear relationship, the temperature of the outer cylinder was found indirectly under the assumption that the flow and temperature distribution in the thermal flask are not altered by the high voltage. The Nusselt number now gave a value of unity in the conduction regime. 3 9 c = i STAINLESS STEEL ^ . . . . ^ FIGURE 6. A CROSS SECTION OF THE APPARATUS.. 40 3. 2 Experimental Arrangement The apparatus, detailed in.figure 6, consisted of two coaxial cylinders (1 and 2) of stainless steel- The cylinders were polished to a mirror finish to avoid any sharp points which might cause electric breakdown at high voltages. The inner cylinder (2) was,constructed in two parts. The central part was made.from brass tubing in which holes at 8, 9 and 10 were dr i l l e d to accept thermistors,, A tight f i t t i n g stainless steel tubing was then forced over the brass tubing. End pieces were machined from stainless steel and f i t t e d to the central part with inside corners at the lines of contact to prevent the occurrence of corona. The heating coil.(7), made of a close wound 30 S.W.G. diamel coated'Mancoloy 10' wire on afibreglass rod, was a close.fit inside the inner cylinder. After installing calibrated thermistors in the receptacles at.8, 9.and 10, the. heating c o i l was centered, and the remaining space inside the inner cylinder was f i l l e d with thermal epoxy and allowed to set. Besides the end points, the heating c o i l was also tapped at two intermediate points directly opposite the thermistors at 8 and 10. The leads for the thermistors and the heating c o i l were brought out the top of the cylinders. The inner cylinder was located coaxially inside the outer cylinder (1.) by means of plexiglass centering devices at each end. The locations of the centering devices were well beyond the measure- ment length between the points 8 and 10. The outer cylinder was. equipped with toroidal terminations to prevent the occurrence of — G > POWER SUPPLY R 2 V W 4 42 corona. Table 2 summarizes the important dimensions of the completed c e l l . Table 2 Cell Dimensions. Radius of the inner cylinder Inner radius ,= 1.455 ± 0.002 cms. Outer radius = 1.711 ± 0.002 cms. Radius of the .outer cylinder Inner radius = 1.903 + 0.002 cms. Outer radius = 2.114 ± 0..002 cms. Length of the measurement section between the points 8 and 10 L = 8.25 ± 0.02 cms. The assembled apparatus was set in a double walled thermal flask (5) f i l l e d with transformer o i l . This o i l was circulated through an external thermostat, thus enabling the temperature of the outer cylinder to be maintained constant. The temperature in the thermal flask could be controlled to ±0.02°C. The present.setup didnot allow to observe f l u i d motions visually. The onset of thermal convection was detected by making heat transfer and temperature measurements. The Heating Circuit: The control circuit for the heating c o i l is shown in figure 7. As already mentioned, the heater was a three sectioned heating element. Currents in the three sections were adjusted by means of the rheostats R̂  and unti l the thermistors at 8, 9 and 10 gave the same.temperature reading (control to about 0.01°C was found feasible) to ensure uniformity of.temperature 43 on the surface .of .the inner cylinder. At equilibrium the heat transfer occurring in the..section.between the points 8 and 10 was taken to be equal to the dissipation of el e c t r i c a l energy in the central section of the c o i l . The ammeters 1^ and and the voltmeter were a l l of .1% accuracy (made by Weston, Model 81). The current for the heating element was drawn from a direct current, (Hewlett Packard, model 6290A) power supply. The power input in the central, section .of the coil . i s . given by (i-j.. - i 2 ) v1 watts The Thermistors: The glass coated thermistors of .043 in. diameter had a nominal resistance of 50,000 ± 20% ohms at 25°C and a dissipation.constant.in air of 0.7 milliwatts per °C. A l l thermistors were calibrated before use. The calibration was per- formed against a.N..P.L..calibrated platinum resistance thermometer. The.calibration points were recorded over a 15-28°C interval. The results were analyzed,by the. method .of. least squares for f i t .to. the. orthogonal, polynomials for the function lnR =.A + B/T + C/T2 where ,R .is... the. resistance of the thermistor at temperature T°C and A,- B...and. C. are .the. constants, of .the .cubic equation. .The results of calibration and the difference between the calibration, resistance and. the...computed resistance. are l i s t e d in tables 3 and 4 (Appendix 1). - o - -O o S j = 2 9 9 9 2 . 3 O H M S S 2 = 2 9 9 9 2 . 8 O H M S R D E C A D E R E S I S T O R T, , T 2 E T C . T H E R M I S T O R S FIGURE 8. THE BRIDGE CIRCUIT. 4 5 • I _J I 1 15 25 35 TEMPERATURE IN DEGREES CENTIGRADE FIGURE 9. KINEMATIC VISCOSITY AS A FUNCTION OF TEMPERATURE.  47 The resistance of the thermistors was read with a direct reading bridge specially designed for this purpose (figure 8). The High Voltage Supply: In a l l experiments, the inner cylinder was grounded and the outer cylinder was kept at a high alternating voltage. The 60 Hz high voltage was supplied from a step-up transformer made by the Universal Voltronics Corporation. A voltage stabilizer was inserted between the mains and the input to the transformer. This eliminated the fluctuations of the power line and gave a constant output voltage. A high accuracy (1%) electrostatic voltmeter was used to measure the applied potential difference between the cylinders, The Dielectric Liquid: The dielectric liquid for these experiments was DC200 el e c t r i c a l grade silicone o i l (supplied by Dow Corning) whose physical.properties are listed in table 5. The silicone o i l was selected because i t has a high electric strength and is available commercially in any desired viscosity grade. The dynamic viscosity and the dielectric constant of the o i l as a function of temperature were measured using a f a l l i n g b a l l viscometer and a standard capacitive c e l l respectively. In figure 9, the kinematic viscosity of the o i l is plotted as a function of temperature. Figure 10 shows the variation of dielectric constant with temperature. Other properties listed in.table 5 are taken from Dow Corning Bulletin 05-213 of July 1968. 48 Table 5 Physical Properties of Silicone O i l * Density,. p--0.9377 gm. cm -3 Coefficient: of volume expansion, a- 1.08x10 Thermal conductivity, x = 3.2x10 -4 cal. cm -1 sec Specific heat,. Cp= 0.43 cal. gm. Dielectric constant, K =2.64 e Kinematic viscosity,.v = 12.225 cm sec. -1 *A11 properties are given at 23°C. 3.3 Experimental Procedure : The cylindrical annulus, previously cleaned, was f i l l e d with degassed electric grade silicone o i l . The cavity was.further connected to a vacuum pump to remove any tiny bubbles sticking to the walls of the cylinders. Having made a l l e l e c t r i c a l connections, the external thermostat was switched on. . In a l l experiments, the temperature, of the outer cylinder, T^, was kept lower than the temperature of the inner cylinder, T » cylinder, T̂ > with the mean of the inlet and outlet temperatures, was done without any high voltage across the cylinders. For this, two thermistors were attached to the wall of the outer cylinder. The inner cylinder was heated and when the.conditions became stable, the temperatures at points 8, 9 and 10 (figure 6) were measured. The currents in the three sections of the heating c o i l were adjusted An i n i t i a l run, to correlate the temperature of the outer - n <r> 70 rn I — 1 2 tn z o o TI ro — o z m rn H i o c —i r~ rn O —I u rn m —i -6 m m - i 5 rn H <? c > 70 —I rn X m > LO rn 70 70 > —I m o -Tl m o c —I rn 70 o -< z o m 70 r~ m tn tn TEMPERATURE CORRECTION IN CENTIGRADE DEGREES 6<7 50 with rheostats R̂  and (figure 7) u n t i l a l l the temperature readings on the surface of the inner cylinder were within 0„01°G. The temperature of the outer cylinder and that of the inlet and outlet points,were a l l noted. After 10 minutes, a l l temperatures were measured again and i f there was.no change then the readings were accepted. The procedure was repeated for several values of heat transfer.rate across the. cavity.. It was found that the actual temperature of the outer cylinder less the mean temperature read by the inlet and outlet thermistors (11 and ,12 in. figure 6) i s linear in the heat transfer.rate between the. cylinders as shown in figure 11. . , When the. heat transfer, and.temperature measurements were done with the outer cylinder, at.a.high alternating potential, i t was not possible to .measure its.temperature. Since the flow and heat transfer 1 .in .the circulated..transformer o i l f i l l i n g the thermal flask . .iare. likely, to.be,'largely, unaffected, by the ..potential of the outer cylinder, :itwasi;assumed that .this, is precisely so. Knowing the . heat transf er..,rate..and. the mean, of ..inlet-outlet temperatures, the plot, shown -in, figure. 11 provided, an indirect means of obtaining the temperature of, .the, outer cylinder. ...The experiments.were run for.the. following values of potential difference across the gap. 0 kv rms,. A.06 kv rms, 6.00 kv- rms, 6-92 kv rms, .7.91 kv rms, 9.63 kv rms, 10.15 kv . rms,.. 10.80 kv rms. -1 X IN M A 1 0 0 M A R A N G E FIGURE 12. DIFFERENCE BETWEEN THE READING OF AMMETER ^ A N D LABORATORY STANDARD AS A FUNCTION OF THE AMMETER READING FOR DIFFERENT CURRENT RANGES. THE CURRENT DIFFERENCE IS PLOTTED AS Y AND THE AMMETER READING AS X. Y IN MA +10- -10 X IN MA 100 MA RANGE X IN MA 300 300 MA RANGE 0.6 X IN AMP 1 A RANGE FIGURE 13. DIFFERENCE BETWEEN THE READING OF AMMETER I 2 AND LABORATORY STANDARD AS A FUNCTION OF THE AMMETER READING FOR DIFFERENT CURRENT RANGES. THE CURRENT DIFFERENCE IS PLOTTED AS Y AND THE AMMETER READING AS X . Y IN VOLTS +.05 - -05 X IN VOLTS 2.5 V RANGF +.05 -.05- X IN VOLTS 5.0 V RANGE X IN VOLTS 10 10.0 V RANGE FIGURE 14. DIFFERENCE BETWEEN THE READING OF THE VOLTMETER V , AND LABORATORY STANDARD AS A FUNCTION OF THE VOLTMETER READING FOR DIFFERENT VOLTAGE RANGES. THE VOLTAGE DIFFERENCE IS PLOTTED AS Y AND THE VOLTMETER READING AS X. 54 A l l the observed data is tabulated in tables 6-13 (Appendix II). 3.4 Errors and corrections. The ammeters 1^ and I^ and the voltmeter (figure 7) were calibrated against a laboratory standard of 0,1% accuracy. The results of calibration are shown in figures 12-14. A l l values of currents and voltages were reduced using these graphs. The corrected values of i , i and v and the average temperatures of the inner and outer cylinders are tabulated in tables. 14 through 21 (Appendix III). :Errors-,..due..-to. self heating of ..the, thermistors: The thermistors were specified to have a. dissipation, constant in air, of..0.7 milliwatt per °C. Accordingly, a, .power, of .0.7 .milliwatt in..the thermistor w i l l change its- temperature by. 1°C which w i l l i n turn change i t s resistance, by about 5%. This i s referred to .as ..a self-heating error. It can be reduced by.,properly, .selecting the bridge voltage. ,. , The.: self-heat effect, was. restricted, to less than the 0<01°C allowable precision..,, Taking, a.reasonable .amount of error due to selfr-heating-to,.be.,50%,,o£,. allowable precision or 0..005°G, the max- -6 imum..power.„al.lowed. through, the. thermistor, is. 0.7 x 0.005 = 3.5 x 10 wat Taking an average value for the resistance of the thermistor to be 50,000fi, we get for the voltage across the thermistor a —6 3 value 3.5x10 x50xl0 or 0.42 volts. We can, therefore, use a maximum of 0.84 volts across the bridge. 1 2 3 4 5 HEAT TRANSFER RATE IN WATTS FIGURE 15. DEPENDENCE OF TEMPERATURE DIFFERENCE ON HEAT TRANSFER RATE FOR 0 KV RMS BETWEEN THE CYLINDERS. HEAT TRANSFER RATE IN WATTS F1GURE*16. DEPENDENCE OF TEMPERATURE DIFFERENCE ON HEAT TRANSFER RATE FOR 4.06 KV RMS BETWEEN THE CYLINDERS. FIGURE 17. DEPENDENCE OF TEMPERATURE DIFFERENCE ON HEAT TRANSFER RATE FOR 6.00 KV RMS BETWEEN THE CYLINDERS. LU O GO Z LU LU LU cr cr LU (3 LL. LU LL- Q 1 4 HEAT TRANSFER RATE IN WATTS FIGURE 18. DEPENDENCE OF TEMPERATURE DIFFERENCE ON HEA'i iR^-S" RATE FOR 6.92 KV RMS BETWEEN THE CYLINDERS. FIGURE 19. DEPENDENCE OF TEMPERATURE DIFFERENCE ON HEAT TRANSFER RATE FOR 7 .91 KV RMS BETWEEN THE CYLINDERS. LU 1 2 3 4 5 6 HEAT TRANSFER RATE IN WATTS FIGURE 20. DEPENDENCE OF TEMPERATURE DIFFERENCE ON HEAT TRANSFER RATE FOR 9.63 KV RMS BETWEEN THE CYLINDERS. 8 HEAT TRANSFER RATE IN WATTS FIGURE 2 1 . DEPENDENCE OF TEMPERATURE DIFFERENCE ON HEAT TRANSFER RATE FOR 10.15 KV RMS BETWEEN THE CYLINDERS. LU O CO z : LLJ U J L J J cr cr L U o U _ L U U _ Q Q L U L U or 3 < cr UJ j f j 2 L u z — — / — 10.80 KV — J . 1 I 1 1 1 1 1 2 3 4 5 6 HEAT TRANSFER RATE IN WATTS FIGURE 22 . DEPENDENCE OF TEMPERATURE DIFFERENCE ON HEAT TRANSFER RATE FOR 10.80 KV RMS BETWEEN THE CYLINDERS. 63 A bridge voltage of 0.5 volts was used for the resistance measurements which reduced the error due to self-heating far below the required precision. Error due to f i n i t e thickness., of .walls of the cylinders: The temperature of the inner cylinder was measured by the thermistors mounted on the inner wall„; The temperatures of the outer cylinder were those at the outer surface. For any true calculation of conductive heat transfer through the liquid, i t is essential to correct for the fi n i t e wall thicknesses... The thicknesses of the walls were sufficiently small that i t was estimated that the maximum error introduced is less than 0.5% and i t was not considered necessary to make this correction. Errors due to .internal resistance of 1^, 1^ and V^: The sensitivity of the voltmeter. V^, was 20,0000/volt. The resistance of.the ammeters was 0.10 on 1A range compared to about'130.for the mid-section of the heating c o i l . . It was, therefore, found unnecessary to do any corrections for^the. internal resistance, of the meters. .3.5 Experimental Results. Experimental results, in the form.of.the temperature difference between.the cylinders .plotted against the heat .transfer rate in the central section of the annulus are shown in figures 15 through .22 for the range of voltages applied to the. outer, cylinder. For potentials up.;to 7.91 kv rms, the relations are linear, indicating & 9.63 KV o 10.15 KV © 1.3 _ © 10.80 KV 1.2 — 1.1 © 1.0 © © © I 1 i O & O O 1 1 1 1 500 1000 1500 2000 2500 3000 35.00 E L E C T R I C A L RAYLEIGH N U M B E R FIGURE 2 3 . NUSSELT NUMBER AS A FUNCTION OF ELECTRICAL RAYLEIGH NUMBER. 6 5 conductive heat transfer even for the maximum temperature differences attained. For potentials of 9 . 6 3 , 1 0 . 1 5 and 1 0 . 8 0 kv rms, distinct breaks occur in the linear relationship indicating the onset of convective heat transfer. A l l the results are summarized in figure 2 3 as a plot of the Nusselt number v 1 ( l 1 - 1 2) l n ( r 2 / r i ) W U 4 . 1 8 5 4 [ 2 T T X L ( T 1 - T , , ) ] against the el e c t r i c a l Rayleigh number, R (equation 5 2 , Chapter II). In the above equation x i s the thermal conductivity of the f l u i d and L is the length of the measurement section between the points 8 and 1 0 (see figure 6 ) . The values of viscosity and dielectric constant entering in the expression for R were taken at the mean of the temperatures of the cylinders. The thermal conductivity depends so slightly on the . temperature that i t s variation with temperature was ignored. Other constants used for the oil,are l i s t e d in Table 5 . Figure 2 3 shows the experimentally determined c r i t i c a l e l e c t r i c a l Rayleigh number for the onset of in s t a b i l i t y to be R = 2 2 0 0 ± 1 0 0 . c 3 . 6 Accuracy The accuracy of the f i n a l results for Nu and R i s estimated to be,±5% when a l l the factors are considered. 66 CHAPTER IV DISCUSSIONS AND. CONCLUSIONS 4 . 1 Discussion of Results The figures 15 to 22 show the plots of the heat transfer rate against the temperature difference between the cylinders for a range of voltages applied to the outer cylinder,, For lower voltages the plots are linear, indicating the heat transfer i s by conduction alone. Above 7.91 kv rms, however, the plots show distinct breaks which occur at the onset of instability (Schmidt and Milverton, 1935). The ele c t r i c a l buoyancy force for small voltages was not sufficient to overcome thermal and viscous dissipation even for the maximum temperature differences created across the gap. However, i f the temperature of the inner cylinder i s raised further, we should expect appreciable convective heat transfer even for those voltages. The plot of the Nusselt number against the el e c t r i c a l Rayleigh number is shown in figure 23. In the conduction regime, the Nusselt number i s found to have a slight downward trend from the expected value of unity. If there i s a small temperature drop at the outer wall contact point during the i n i t i a l run with no high voltage, this could lead to an underestimation of the temperature difference of figure 11. Consequently, the impressed temperature difference between the cylinders would be overestimated. A linear increase of this contact temperature difference with impressed temperature difference could result in a small decrease in the Nusselt number with increasing Rayleigh number. The value of the e l e c t r i c a l Rayleigh number at which the onset of inst a b i l i t y occurs i s : R = 2200 ± 100. c The theoretically determined value of R i s 2.119.346 and compares 67 well with the experimental value. The disturbances which w i l l be mani- fested at marginal s t a b i l i t y are characterized by the wave length A = ^ = ^ = 2 . 0 2 d . k a c The value of a^ for this case is same as for the Rayleigh-Benard convection (ac=3.117). In the Rayleigh-Benard i n s t a b i l i t y , the Nusselt number i s a function of the Rayleigh number only. Here, however, the Nusselt number depends on the e l e c t r i c a l Rayleigh number and the ratio of the inner cylinder radius to gap width. 4.2 Conclusions The present investigations have demonstrated the f e a s i b i l i t y of constructing a laboratory model of thermal convection to simulate large scale geophysical flows that occur under a (nearly) central gravitational f i e l d . The marginal s t a b i l i t y analysis of the flow was done by making several simplifying assumptions. The general validity of the principle of exchange of s t a b i l i t i e s could not.be proved for this problem. However, the quality of the agreement of the experimental results with the numerically determined value (assuming the principle of exchange of st a b i l i t i e s ) i s a strong verification of the validity of the analysis. In constructing laboratory models of geophysical flows, we need to proceed to spherical geometry and rotation. We should also be interested in the actual form of the flow pattern as well as the values of various parameters at the onset of convection. A l l this w i l l require the development of some new techniques, but in principle i t now appears possible- 68 REFERENCES Avsec, D. and Luntz^ M. 1937, El e c t r i c i t e Et Hydrodynamique - Quelques formes nouvelles des tourbillons electroconvections, C. R. Acad. Sci., Paris, 204, 757. Batchelor, G. K., 1954, Heat transfer by free convection across a closed cavity between vertical boundaries at different temperatures, Quart. Appl. Math., 12, 209-33. Be*nard, H.., 1900, Les tourbillons cellulaires dans une nappe liquide, Revue generale des science pures et aplique*es, 11, 126L-71 and 1309-28. Bird, R. B., Stewart, W.B. and Lightfoot, E. N., 1966, Transport phenomena, John Wiley, New York. Bullard, E. C. and Gellman,. H., 1954, Homogeneous dynamos and te r r e s t r i a l magnetism, Phil, trans. Roy. Soc, Sr A, 247, 2 1 3 - Chandra, B., 1969, Thermal convection under an axially symmetrical force f i e l d , M. Sc. Thesis, University of Western Ontario, London, Canada. Chandra, K. , 19.38, Instability of fluids heated from below, Proc. Roy. Soc. (London), Sr A, 164, 231-42. i...Chandrasekhar, S., 1961, Hydrodynamics and hydromagnetic s t a b i l i t y , Oxford Univ. Press. Debye, P., 1929, Polar molecules, Chemical catalog company. New , York, p27. Ford, L. R., 1955, Differential equations, McGraw H i l l . Fultz, D.., 1961,. . Developments, i n controlled experiments on large scale. geophysical, problems,. Advances..in .Geophysics,, edited .by H. E.. .Landsberg, 7, 1-103. Greenspan, H. E.,. The.theory of rotating fluids, Cambridge. Gross, M. J., 1967, Laboratory analogies for convection problems, Mantles of .the earth, and te r r e s t r i a l planets, edited by S. K. -Runcorn, .Interscience, 499-503. 69 Gross, M. J. and Porter, J. E., 1966, Elec t r i c a l l y induced convection in dielectric liquids, Nature, no. 5058, 1343-45. Hadley, G., 1735, Concerning the cause of the general trade winds, Reprinted in the Mechanics of the Earth's atmosphere, Smithsonian Inst. Misc. Coll., 51, no. 4, 1910. Ince, E. L., 1956, Ordinary differential equations, Dover Publications. Jeffreys, H., 1926, The i n s t a b i l i t y of a layer of f l u i d heated below, Phil. Mag.., A, 2, 833-44. Kronig, R. and Schwarz, N., 1949, On the theory of heat transfer from a wire in an electric f i e l d , Appl. Sci. Res., A l , 35-45, Landau, L. D. and L i f s h i t z , E» M., 1960, Electrodynamics of continuous media,.Pergamon press, London. Long, R. R.. (editor), 1955, Fluid models.in geophysics, Proc. 1st . symposium on the use of. models in geophysical f l u i d dynamics, Bait. Low, A. R., 1929, On the criterion for the st a b i l i t y of a layer of viscous f l u i d heated below, Proc. Roy. Soc, A, 125, 180-95. Malkus, W. V. R. and Veronis, G., 1961, Surface electroconvection, Phys. Fluids, 4, 13-23. Mathews, J. and Walker, R. L.,.1964, Mathematical methods in Physics, W. A. Benjamin, Inc., New York. Melcher, J. R. and Taylor, G. I., 1969, Electrohydrodynamics: A review of the role of interfacial shear stresses, Annual Review of Fluid Mechanics edited by W. R. Sears and M. Van Dyke, 1, 111-47. McFadden, C. P., 1969, The effect of region of low viscosity on thermal convection in the Earth's mantle, Ph. D. thesis, University of Western Ontario, London, Canada. Pellew, A. and Southwell, R . V., 1940, On maintained convective motion in a f l u i d heated from below, Proc. Roy. Soc,, A, 176, 312-43. 70 Rayieigh, Lord, 1916, On convective currents in a horizontal layer of f l u i d when the higher temperature i s on the underside, Phil. Mag., 6, 529-46. Rossby, C. G., 1941, The s c i e n t i f i c basis of modern meteorology, Yearbook of Agriculture, U. S. Dept. of Agriculture, 599-655. Schmidt, R. J. and Milverton, S. W*, 1935, On the ins t a b i l i t y of a f l u i d when heated from below, Proc. Roy. Soc. (London), A, 152, 586-94. Silveston, P. L., 1958, Warmedurchgang in Waagerechten FlUssigk-eitsschichten, part 1, Forsch Ing. Wes., 24, 29-32 and 59-69. Smylie, D. E., 1966, Thermal convection in dielectric liquids and modelling in geophysical fluid-dynamics, Earth and Planetary Science Letters, 1, 339-40. Turnbull, R. J i,1968, Electroconvective Instability with a stabilizing temperature gradient, Phys. Fluids, 11, 2588-96 and 2597-2603. Vening-Meinesz, F. .A., 1962, Thermal convection in the Earth's mantle, Continental Drift.edited by S. K. Runcorn, Academic Press. Vest, C. M. and Arpaci, V. S. , .1969, Stability of natural convection in a vertical slot,.J. Fluid Mech., 36, 1-16. 71 APPENDIX I TABLES OF THERMISTOR CALIBRATION 72 TABLE NO. 3 TABLE GF CAL IBRAT ICN POINTS FOR THERMISTORS i INNER C Y L I N D E R . ( 1 ) ( 2 ) ( 3) TEMPERATURE R E S I S T A N C E RES ISTANCE RES I STANCE DEGREES C OH N S CHMS OHMS 1 5 . 0 3 7 7 9 S 1 2 7808a 7 2 1 2 7 1 7 . 7 14 7 0 3 6 4 6 8 5 5 9 6351 7 2 0 . 4 5 0 6 1 9 2 3 6 0 1 5 5 5 5 9 0 2 2 3 . 166 54 6 71 5 2 9 5 2 4 9 3 5 2 2 5 . 8<58 4 8 3 1 0 4 6 6 6 8 4 3 6 2 5 2 8 . 6 4 6 4 2 7 7 9 4 1 1 9 4 3 8 6 1 4 D I F F E R E N C E BETWEEN C A L I B R A T E D R E S I S T A N C E AND COMPUTED R E S I S T A N C E . ( 1 ) ( 2 ) ( 3) TEMPERATURE RES IS T ANC E RES ISTANCE RES ISTANCE DEGREES C OHMS OHMS CHMS 1 5 . C 3 7 - 0 . 81 0 . 2 5 1 . 1 3 1 7 . 7 1 4 - G . 8 8 . - 0 . 2 5 - 1 . 4 0 2 0 . 4 5 0 0 . 6 7 - 0 . 9 6 - 0 . 8 6 2 3 . 1 6 6 6 . 42 2 . 9 8 1 . 3 9 2 5 . 8 9 8 - 7 . 4 2 - 2 . 5 1 0 . 24 2 8 . 6 4 6 4 . 1 0 0 . 7 4 - 0 . 3 6 ( 1)THE RMIS TOR AT POINT 8 ON THE INNER C Y L I N D E R . ( 2 ) T H E R M I S T 0 * AT POINT 9 CN THE INNER C Y L I N D E R . ( 3 J T H E R M I S T C R AT POINT 1C ON THE INNER C Y L I N D E R . 73 ABLE N TABLE POINT TEMP£ DEGR 11 . 13 . 1 5 . 1 7 . 19 . ?1 . CF THE RATI EES 9 14 C 68 8 2 4 742 677 6 2 4 C L I 34ATI C.\ MV. ISTCWS. F C I M S FOR INLET AND CUT L E I ;KE C ( 1 ) ;< r S I S T i\C E CHrlS 9 C 4 8 1 a ? 16 5 7 4 69 7 6 81 C4 6.? 08 2 5 6 6 3 2 ( 2 i RES IS VANCE CH.MS 7 9 097 7 18 16 6 52 S3 59 505 54242 49 4 66 DIFFERENCE. BETWEEN C A L I B R A T E D R E S I S T A N C E AND COMPUTED R E S I S T A N C E . TEMPERATURE DEGREES C 1 1 . 9 14 1 3 . 8 6 8 1 5 . 0 2 4 1 7 . 7 4 2 19 . 6 7 7 21 . 6 2 4 ( 1 ) R E S I S T A N C E OHMS 3.4 4 - 6 , -0, 7. - 4 , 0, 6.1 44 Cc 19 8,3 ( 2 ) RES ISTANCE OHMS 3 . 69 - 5 . 2 5 2 . 0 6 2 . 7 0 - 0 . 4 6 - 0 . 5 3 (1) T HE RMIS TOR AT PCINT 11 AT THE I N L E T POINT TC THE THERMAL F L A S K . (2) THERKISTGK AT FC I NT 12 AT THE OUTLET POINT TO THE THERMAL F L A S K . \ 74 APPENDIX II TABLES OF HEAT TRANSFER MEASUREMENTS 7 5 TABLE NO. "6 VOLTAGE BETWEEN THE C YLINDER S = 0.0 KV (1) (2) {J ) (4) (5) (6) (7) (8) 16. 34 1 6.35 16. 3 5 16 .61 16.82 0.0 0.0 0.0 17.58 17 . 57 17. 57 16.61 16. 78 0. 263 0. 073 2.57 18 .02 18.01 18 . 00 16 . 59 16.73 0.304 0.084 2.9 8 13. A9 18. 50 IB. 49 16. 61 16.80 0.349 0.096 3.4? 19.29 19.28 19 . 28 16.62 16.82 0.402 0 . 108 3. 98 20. Cl 2 0. 00 2C. 00 16. 63 16.78 0.446 0.118 4 .45 20.89 20.88 20. 8 3 16. 62 16. 78 0. 496 0. 132 4. 95 21.91 21.90 21 . 89 16.64 16 .79 D .543 0. 138 5.48 22.98 22.97 22. 97 16. 60 16. 74 0.596 0. 154 6. 02 24.13 24.13 2 4 . 13 16 .60 16.74 0.643 0. 162 6.55 25.4 5 25.45 25. 4 5 16.-73 16.83 0.698 0.177 7.05 26.78 26.77 26. 78 16.65 16.79 C.749 C. 189 7. 58 27. 82 2 7. 82 27. 83 16 .64 16.80 0.787 0.201 7.97 Q ) :T EMPERATURE 0 F INNER CYLINCER AT POINT 10. (2) TEMPERATURE CF INNER CYLINDER AT POINT 9. (3) TEMPERATURE CF INNER CYLINDER AT POINT 0. (4) TEMPERATURE CF CIRCULATING OIL AT INLET POINT 11. (5) TEMPERATURE CF CIRCULATING OIL AT OUTLET POINT 12. (6) CURRENT II IN AMPERES. (7) . CURRENT 12 IN AMPERES. (8) VOLTAGE VI IN VOLTS. 76 TABLE NG. 7 VOLTAGE BETWEEN THE CYLINDERS3 4.06KV (1) (2) (3) (4) (5) (6) (7) (8) 17.88 17.89 17.88 16.65 16.85 0.303 0.088 2.90 25 .5 1 16.80 0. 353 0. 113 3.50 16.86 0.402 0. 124 3.74 16.83 0.450 0. 141 4 .18 16.89 0.502 Q. 160 4. 65 16.86 0.547 0. 174 5.33 16. 91 0. 594 0. 189 5.50 17.11 0.644 0. 202 5.96 17.06 0.700 0.2 2 0 6 .50 16.92 0.750 0. 236 6. 95 16.98 0.301 0. 255 7.41 17. 01 0. 845 0. 268 7.81 17.02 0.898 0. 286 8. 37 29.33 29.33 29.34 16.78 (1) TEMPERATURE OF INNER CYLINDER AT POINT 10. (2) TEMPERATURE CF INNER CYLINCER AT POINT 9. (3) TEMPERATURE OF INNER CYLINDER AT PCINT 8. (4) TEMPERATURE OF CIRCULATING OIL AT INLET POINT 11. v{ 5 ) TEMPERATURE OF CIRCULATING OIL AT OUTLET POINT 12. (6) CURRENT II IN AMPERES. (7) CURRENT 12 IN AMPERES. (8) VOLTAGE VI IN VOLTS. 77 TABLE NO. B VOLTAGE BETWEEN THE CYL INDERS= 6 .ODKtf ( i ) ( 2) (3 ) ( 4 ) (5 ) (6 ) (7 ) (8 ) 1 9 . 4 9 1 9 . 5 0 19 . 5 0 1 6 . 8 6 1 7 . 0 8 0 . 4 3 7 0 . 139 4 . 03 2 0 . 4 4 2 0 . 4 5 2 0 . 4 5 1 6 . 80 1 7 . 0 2 3 . 4 9 6 0 . 1 5 6 4 . 6 1 21 . 3 6 2 1 . 3 7 2 1 . 3 7 1 6 . 87 1 7 . 15 . 0 . 5 5 2 0 . 1 7 4 5 . 1 3 2 2 . 35 2 2 . 3 5 22 . 3 4 1 6 . 8 7 1 7 . 1 6 0 . 6 0 2 0 . 189 5 . 5 8 2 3 . 26 2 3 . 2 6 2 3 . 2 5 1 6 . 89 1 7 . 2 1 0 . 6 4 9 0 . 2 0 4 6 . 0 2 2 4 . 3 0 2 4 . 2 9 2 4 . 3 0 1 6 . 8 9 1 7 . 1 9 0 . 7 0 1 0 . 2 1 8 6 . 54 2 5 . 54 2 5 . 5 5 2 5 . 53 1 6 . 84 1 7 . 1 2 0 . 7 6 3 0 . 2 4 3 7 . 0 4 2 6 . 7 8 2 6 . 7 9 2 6 . 7 7 1 6 . 8 5 1 7 . 1 2 0 . 8 C 1 0 . 2 5 3 7 . 4 5 2 7 . 9 1 2 7 . 9 0 2 7 . 9 0 1 6 . 8 3 1 7 . 0 4 0 . 8 4 8 0 . 2 7 3 7 . 8 5 (1) TEMPERATURE OF INNER CYLINDER AT POINT 1 0 . (2) TEMPERATURE CF INNER CYLINDER AT POINT 9 . (3) TEMPERATURE CF INNER CYL INDER AT P C I N T 8 . (4) TEMPERATURE CF C I R C U L A T I N G O IL AT INLET POINT 1 1 . (5 ) TEMPERATURE OF C I R C U L A T I N G OIL AT OUTLET POINT 1 2 . (6) CURRENT I I I N A M P E R E S . (7 ) CURRENT 12 IN A M P E R E S . (8 ) VCLTAGE VI IN V O L T S . 78 TABLE N O . 9 VOLTAGE BETWEEN THE CYL INDERS= 6 . 9 2 K V '(I). ( 2 ) ( 3 ) (4 ) (5 ) lb) ' i l ) ( 8 ) 1 9 . C 5 1 9 . 05 1 9 . 0 6 1 6 . 7 0 1 7 . 0 2 Q . 4 0 4 0 . 1 2 6 3 . 7 5 1 9 . 6 7 1 9 . 6 7 1 9 . 6 7 1 6 . 7 4 1 7 . 0 0 0 . 4 4 8 0 . 1 4 1 4 . 1 7 2 0 . 4 4 2 0 . 4 4 2 0 . 4 3 1 6 . 7 1 1 7 . 0 8 0 . 4 9 9 0 . 1 5 7 4 . 6 3 2 1 . 28 2 1 . 28 2 1 . 26 1 6 . 75 1 7 . 0 4 0 . 5 4 7 0 . 1 73 5 . 0 8 2 2 . 2 5 22 . 2 5 22 . 2 4 1 6 . 7 9 1 7 . 0 6 0 . 6 0 0 0 . 189 5 . 57 2 3 . 3 0 2 3 . 3 0 2 3 . 3 0 1 6 . 7 3 1 6 . 9 9 0 . 6 5 4 0 . 2 0 6 6 . 0 7 2 4 . 2 2 2 4 . 2 1 2 4 . 2 1 1 6 . 6 8 1 6 . 9 8 0 . 7 0 0 0 . 2 2 4 6 . 4 4 2 5 . 3 3 2 5 . 3 3 2 5 . 3 3 1 6 . 7 2 1 6 . 9 7 0 . 7 4 7 0 . 2 3 1 6 . 8 4 2 6 . 9 1 2 6 . 9 2 2 6 . 9 1 1 6 . 7 2 1 6 . 9 8 0 . 8 1 6 0 . 2 6 1 7 . 5 4 2 8 . 3 0 2 8 . 3 8 20 . 3 8 1 6 . 7 2 1 6 \ 9 6 0 . 8 7 1 0 . 2 7 5 0 . 10 2 9 . 4 3 2 9 . 4 4 2 9 . 4 4 1 6 . 7 2 1 6 . 9 6 0 . 9 1 7 0 . 2 9 1 8 . 5 0 (1) TEMPERATURE GF INNER CYLINDER AT POINT 1 0 . (2) TEMPERATURE CF INNER CYL INDER AT PCINT 9 . (3) TEMPERATURE OF INNER CYL INDER AT PCINT 0 . ( 4 ) TEMPERATURE CF C I R C U L A T I N G O I L AT INLET POINT 1 1 . (5) TEMPERATURE CF C I R C U L A T I N G OIL AT OUTLET POINT 12 . (6) CURRENT I I I N A M P E R E S . (7) CURRENT 12 IN A M P E R E S . (8 ) VOLTAGE VI IN V O L T S . 79 TABLE NO . 10 VOLTAGE BETWEEN THE CYL INDERS= 7 . 9 1 K V ( I ) ( 2 ) (3 ) ( 4 ) (5 ) (6 ) ( 7 ) ( 8 ) 1 9 . C4 1 9 . 04 1 9 . 0 4 1 6 . 6 6 1 6 . 9 2 0 . 4 0 2 0 . 1 2 8 3 . 7 2 1 9 . 7 5 1 9 . 7 5 1 9 . 7 5 1 6 . 67 16 . 89 0 . 4 5 2 C. 139 4 . 2 3 2 0 . 6 2 2 0 . 6 2 2 0 . 6 1 1 6 . 7 0 1 6 . 8 6 0 . 5 0 9 0 . 1 5 5 4 . 7 8 2 1 . 6T 2 1 . 68 2 1 . 67 1 6 . 68 1 6 . 9 0 0 . 5 6 0 0 . 1 6 8 5 . 3 2 2 2 . 8 3 2 2 . 8 2 2 2 . 8 2 1 6 . 6 7 1 6 . 9 2 0 . 6 2 7 0 . 1 9 6 5 . 8 3 2 4 . 2 2 2 4 . 2 3 2 4 . 2 3 1 6 . 7 1 1 6 . 9 5 0.700 0 . 2 2 4 6 . 4 3 2 5 . 5 7 2 5 . 5 7 2 5 . 5 7 1 6 . 74 1 6 . 94 0 . 7 6 9 0 . 2 5 2 6 . 9 8 2 6 . 6 8 2 6 . 6 7 2 6 . 6 7 1 6 . 6 9 1 6 . 9 2 0 . 8 2 8 0 . 2 7 9 7 . 4 4 2 8 . 3 0 2 8 . 3 0 2 8 . 3 0 1 6 . 64 1 6 . 82 0 . 8 9 5 0 . 291 8 . 2 0 (1) TEMPERATURE CF INNER CYL INCER AT PC I NT 1 0 . (2) TEMPERATURE OF INNER CYL INDER AT POINT 9 . (3 ) TEMPERATURE OF INNER CYLINDER AT POINT 8 . (4 ) TEMPERATURE CF C I R C U L A T I N G OIL AT INLET POIMT 1 1 . (5 ) TEMPERATURE OF C I R C U L A T I N G O I L AT CUTLET POINT 1 2 . (6) CURRENT I I IN A M P E R E S . (7 ) CURRENT 12 IN A M P E R E S . (8) VOLTAGE VI IN V C L T S . 80 TABLE NO. 11 VOLTAGE BETWEEN THE CYLINDERS* 9.63KV (1 ) (Z ) (3 ) (4) (5) (6) (7) (3) 18.41 18.42 18. 43 16. 71 16.83 0.351 0.1.0P 3.25 19.05 19.06 19.06 16.69 16.89 0.4G2 0. 125 3. 73 19.74 19.75 19.74 16.65 16.89 0.453 0.144 4.20 2C.66 20.66 20.65 16. 59 16. 89 0.51.7 0. 166 4. 77 21.52 21.52 21.51 16.66 16.88 0.569 0.187 5.18 22.85 22.84 22.84 16.69 16.86 0.543 0.213 5.83 23.98 23.97 23.97 16.69 16.86 0.707 0.235 6.38 24.91 24.91 24.91 16.68 16.88 0.757 0.250 6.85 25.69 25.70 25.70 16. 67 16. 81 0. 8C0 0.264 7.25 26.68 26.68 26.68 16.65 16.85 0.855 0.279 7.32 27.53 27.52 27.52 16.63 16.81 0.898 0.293 8.24 28.72 28.72 28.73 16.71 17.01 0.953 C.300 8.84 29.65 29.65 29.66 16.71 16.85 0.998 0.312 9.94 11) TEMPERATURE OF INNER CYLINCER AT POINT 10. (2) TEMPERATURE OF INNER CYLINCER AT PCINT 9. (3) TEMPERATURE OF INNER CYLINDER AT POINT 8. (4) TEMPERATURE OF CIRCULATING OIL AT INLET POINT 11. (5) TEMPERATURE CF CIRCULATING OIL AT OUTLET POINT 12. (6) CURRENT II IN AMPERES. (7) CURRENT 12 IN AMPERES. (8) VOLTAGE VI IN VOLTS. 81 TABLE NO . 12 VOLTAGE BETWEEN THE CYLINDERS' 10.15KV (1) (2) (3) (4) (5) (6) (7) (8) 17.74 17. 75 17.75 16.67 16.82 0.287 0.088 2.64 18.35 18.35 18.36 16. 66 16. 88 0. 341 0. 103 3.2 0 18.9 2 13.92 10 .9 2 16 .69 16.82 0.391 0. 119 3.65 19 .70 19.71 19.71 16.68 16.85 0.449 0.138 4.21 20.55 20.56 20.55 16.69 16.83 0.5G7 0.161 4.70 21.21 21.20 21.20 16.67 16.86 0.550 0.171 5.02 22.12 22.13 22.12 16.70 16.88 0.611 C.198 5.47 23.01 23.01 23.00 16.65 16.83 0.661 0.215 5.93 23.83 23.83 23.82 16.63 16.00 0.700 0.231 6.35 24.62 24.62 24.62 16.63 16.89 0.744 0.246 6.77 25.57 25.57 25.57 16.66 16.79 0.800 0.264 7.25 26.43 26.42 26.42 16.65 16. 79 0. 851 C.281 7.75 27.39 27.39 27.39 16.65 16.79 D.911 0.291 3.28 28.47 28.46 28.47 16.66 16.80 0.950 0.294 8.85 (1) ' TEMPERATURE CF INNER CYLINCER AT PC I NT 10. (2) TEMPERATURE OF INNER CYLINDER AT POINT 9. (3) TEMPERATURE OF INNER CYLINDER AT POINT 8. (4) TEMPERATURE CF CIRCULATING OIL AT INLET POINT 11. (5) TEMPERATURE CF CIRCULATING 01L AT CUTLET POINT 12. (6) CURREMT II IN AMPERES. (7) CURRENT 12 IN AMPERES. (8) VOLTAGE VI IN VOLTS. 82 T A B L E NO. 13 VOLTAGE BETWEEN THE C Y L I N D E R S =, 1 0 . 80KV (1 ) (2 ) (3 ) ( 4 ) ( 5 ) ( 6 ) ( 7) ( 8 ) 1 7 . S 4 1 7 . 9 5 1 7 . 95 1 6 . 6 6 1 6 . 7 6 0 . 3 1 0 0 . 0 9 5 2 . 8 8 1 8 . A l 1 8 . A 2 1 8 . 4 2 1 6 . 6 1 1 6 . 7 7 0 . 3 5 4 0 . 1 0 8 3 . 3 0 1 9 . C9 " 1 9 . 1 0 1 9 . 0 9 1 6 . 5 9 1 6 . 7 7 0 . 4 0 8 0 . 1 2 6 3 . 8 0 1 9 . 7 7 1 9 . 7 7 1 9 . 7 7 1 6 . 6 6 1 6 . 7 8 0 . 4 6 2 0 . 1 4 7 4 . 2 7 2 0 . 4 5 2 J . 4 5 2 0 . 4 4 1 6 . 6 0 1 6 . 7 6 0 . 5 0 5 0 . 1 6 1 4 . 6 8 2 1 . 4 4 2 1 . 4 4 2 1 . A 3 1 6 . 6 1 1 6 . 7 7 0 . 5 6 8 0 . 1 8 1 5 . 2 0 2 2 . 1 4 2 2 . 1 A 2 2 . 1 4 1 6 . 6 1 1 6 . 7 1 0 . 6 C 7 0 . 1 9 7 5 . 5 5 2 3 . 1 C 2 3 . 1 0 2 3 . 1 0 1 6 . 6 1 1 6 . 7 6 0 . 6 6 2 0 . 2 1 6 6 . 0 5 2 3 . 8 3 2 3 . 8 3 2 3 . 0 2 1 6 . 57 1 6 . 68 0 . 7 0 7 . C . 230 6 . 4 5 2 4 . 4 9 2 4 . 4 9 2 4 . 5 0 1 6 . 6 6 1 6 . 8 2 0 . 7 5 2 0 . 2 4 5 6 . 8 5 2 5 . 3 5 2 5 . 3 4 2 5 . 3 4 1 6 . 6 6 1 6 . 7 8 0 . 7 9 2 0 . 2 5 1 7 . 3 4 2 6 . 10 2 6 . 10 2 6 . 1 0 1 6 . 6 5 1 6 . 7 7 0 . 8 3 0 0 . 257 7 . 75 2 6 . 8 0 2 6 . 8 0 2 6 . 8 1 1 6 . 6 5 1 6 . 8 0 0 . 8 6 9 0 . 2 6 7 8 . 1 4 2 7 . 7 6 2 7 . 7 7 2 7 . 7 7 1 6 . 67 1 6 . 7 6 0 . 9 1 9 0 . 2 8 3 8 . 6 3 2 8 . 7 2 2 8 . 7 2 2 8 . 7 2 1 6 . 6 2 1 6 . 8 1 0 . 9 7 5 0 . 2 9 8 9 . 0 6 (1) TEMPERATURE OF INNER CYL INDER AT POINT 1 0 . (2 ) TEMPERATURE CF INNER C Y L I N C E R AT POINT 9 . (3) TEMPERATURE CF INNER CYL INDER AT PCINT 8 . (4 ) TEMPERATURE OF C I R C U L A T I N G OIL AT I N L E T POINT 1 1 . ( 5 ) TEMPERATURE CF C IRCUL AT ING OIL AT OUTLET POINT 1 2 . (6) . CURRENT I I IN A M P E R E S . (7) CURRENT 12 IN A M P E R E S . ( ( 8 ) VCLTAGE V I IN V O L T S . 83 APPENDIX III TABLES OE CORRECTED DATA 84 T A B L E N O 14 VOLTAGE BETWEEN THE CYLINDERS= 0.0 KV (1) (2) (3) (4) (5) 16.35 16. 3 6 C. 0 0.0 . 0 .0 17. 57 16.82 G. 2658 C.074 2 2.611 18 .01 16.98 0. 3078 0.0858 3. 03G 18.49 17.17 C. 3524 0.0900 3. 462 19.28 17.48 0. 4C48 0. 1100 4. COO 20.00 17. 75 c . 4487 0.1200 4.466 20.88 18.06 c . 49 86 0.1340 4. 9 52 21.90 18 .45 0. 5446 0.14C0 5.510 22 .97 13.95 c. 5966 0.1556 6.049 24.13 19.31 c . 6430 0.1638 6. 568 25.45 19.93 c . 69 74 0.1785 7.058 26.78 20. 29 c. 7482 C.1902 7. 584 27 .82 20.70 0.7860 0.2019 7.971 (1) MEAN TEMPERATURE DF INNER CYLINDER. (2) MEAN TEMPERATURE OF OUTER CYLINDER. (3) CORRECTED I I . (4) CORRECTED 12. (5) CORRECTED VI. 85 TABLE NO 15 VOLTAGE BETWEEN THE CYL IN DEK S= 4.06 KV ( 1) 12) (3) (4) (5) 17.88 16.94 0. 3C68 C.0898 2.947 18.42 17.15 c . 3561 0.1150 3.540 19.03 17. 42 0. 4C48 0.1260 3. 770 19.68 17 .64 0. 4527 0.1430 4. 199 20.4 7 17.98 c. 5044 0.1617 4.661 21 .23 18.29 c. 54 85 C. 1755 5.330 22. 17 13.67 c . 5947 0.1902 5.530 23.23 19. 11 c . 6440 C.2029 5. 990 24 .37 19.51 0. 6994 0.2206 6.520 2 5. 50 19.94 c . 749 3 0.2365 6.962 26 .70 20.32 c . 8C00 0.2553 7.415 27.82 20. 76 C.8432 0.2684 7.812 29.33 21.32 C. 8952 0.2864 0. 370 (1) MEAN TEMPERATURE OF INNER CYLINDER. (2) MEAN TEMPERATURE OF OUTER CYLINDER. (3) CORRECTED I I . (4) CCRRECT ED 12. (5) CORRECTED VI. 86 TABLE NG 16 VOLTAGE BE T-WEEN THE CYLINOERS= 6.00 KV (1) {2 ) (3 ) (4) (5) 19. 50 I 7. 69 C.4393 0.1410 4. 049 20 .45 18.05 0.4985 C.1578 4. 622 21.37 18.51 0.5534 0.1755 5. 160 22.35 18.84 0. 6024 0.1902 5. 610 23.26 19.23 0.5490 0.2049 6. 048 24.30 19.65 C.7003 0.2187 6. 559 25 .54 20.03 G.7622 0. 24 3 3 7. 048 26.78 20.44 C.7999 0.2533 7. 455 27.90 20. 82 C.8463 0.2733 7. 851 (1) NEAN TEMPERATURE OF INNER CYLINDER. <2) MEAN TEMPERATURE OF OUTER CYLINDER. (3) CORRECTED I I . (4) CORRECT ED 12. (5) CORRECTED VI. 87 T A n L E KO 17 VOLTAGE BETWEEN THE CYL INDE RS = 6 . 92 KV ( 1 ) (2 ) (3 ) (4 ) ( 5 ) 1 9 . 0 5 1 7 . 4 8 C . 4 0 6 8 0 . 1 2 8 0 3 . 7 8 0 1 9 . 6 7 1 7 . 7 4 0 . 4 - 5 0 6 0 . 1 4 3 0 4 . 183 2 0 . 4 4 1 8 . 0 4 0 . 5 0 1 4 0 . 1 5 0 8 4 . 6 4 2 2 1 . 2 8 1 8 . 3 5 0 . 5 4 8 4 0 . 1 7 4 6 5 . 1 1 0 22 .2 5 1 8 . 7 5 0 . 6 0 0 4 0 . 1902 5 . 6 0 0 2 3 . 30 1 9 . 08 G . 6 5 3 9 0 . 2 0 6 9 6 . 0 9 7 2 4 . 2 1 19 . 3 6 0 . 6 9 9 4 0 . 2 2 4 6 6 . 4 6 1 2 5 . 3 3 1 9 . 8 1 0 . 7 4 6 3 0 . 2 3 1 6 6 . 8 5 2 2 6 . 9 1 2 0 . 4 G C. 8 1 4 2 0 . 2 6 1 3 7 . 544 28 . 3 8 21 . 0 7 0 . 8 6 8 4 0 . 2 7 5 3 8 . 100 2 9 . 4 4 2 1 . 4 9 0 . 9 1 2 8 0 . 2 913 3 . 5 0 0 Q ) MEAN TEMPERATURE OP INNER C Y L I N D E R . ( 2 ) MEAN TEMPERATURE OF OUTER C Y L I N D E R . (3) CORRECTED I I . (4 ) CORRECTED 1 2 . (5 ) CORRECTED VI . 88 T A PL E NO 18 VOLTAGE BETWEEN THE C Y L I N D E R $ = 7 . 9 1 KV (1) (2 ) (3) (4) { 5) 1 9 . 04 1 7 . 43 C . 4 0 4 9 0 . 1 3 0 0 3 . 7 4 7 1 9 . 7 5 1 7 . 6 9 C . 4 5 4 6 0 . 1 4 1 0 4 . 2 4 7 2 0 . 6 2 10 . 0 3 0 . 5 1 1 2 0 . 1 5 6 8 4 . 788 2 1 . 6 7 1 8 . 4 3 C. 5 6 1 3 0 . 1 6 9 7 5 . 3 5 0 2 2 . 8 2 1 0 . 9 0 0 . 6 2 7 1 C . 1 9 7 1 5 . 860 2 4 . 2 3 1 9 . 4 1 C . 6 9 9 6 0 . 2 2 4 6 6 . 4 5 2 2 5 . 5 7 1 9 . 8 8 C . 7 6 82 C . 2 5 2 3 6 . 9 9 0 2 6 . 6 7 2 0 . 2 7 0 . 8 267 0 . 2 7 9 3 7 . 4 4 5 2 8 . 3 0 2 0 . 9 5 C. 8 9 2 2 0 . 2 9 1 3 3 . 2 0 0 ( 1 ) MEAN T.EMPERATURE OF INNER C Y L I N O E R . (2) MEAN TEMPERATURE OF OUTER C Y L I N C E R . (3) CORRECTED I I . (4 ) CORRECTED 1 2 . <5) CCRRECTED V I . 89 T A C L E Nil 19 VOLTAGE BETWEEN THE CYLINDERS= 9.63 KV (1) (2) (3) (4 ) (5) 18.42 17.20 C.3544 0.1100 3. 307 19.06 17.45 0.4C49 0.1270 3. 760 19.73 17.73 0.4557 0.1460 4. 390 20.66 18.02 C.5191 0.1677 4.777 21.52 18.37 0.57C1 0.1883 5.210 22.84 18.8 8 G .64 30 0.2137 5.860 23.97 19.27 C.7C6 3 0. 2356 6. 404 24 .91 19 .58 0.7562 0.2503 6.863 25.70 20.03 C.7990 0.2643 7.257 26.68 20.59 0.8530 C.2793 7.822 27. 52 21.00 C.8952 0.2933 8. 240 ^8.72 21.82 0.9493 C.3058 8. 840 29.65 22 .26 0.9980 0.3177 9. 235 (1) MEAN TEMPERATURE OF INNER CYLINOER. (2) MEAN TEMPERATURE OF OUTER CYLINDER. (3) CORRECTED I I . (4) CORRECTED 12. (5) CORRECTED V I . 90 TABLE NO 2 0 VOLTAGE BETWEEN THE CYLINDERS= 1 0 . 1 5 KV (1 ) ; ( 2 ) (3 ) ( 4 ) 15) 1 7 . 7 5 1 6 . 9 3 0 . 2 9 0 6 C . Q 8 C 8 2 . 6 8 2 1 8 . 35 1 7 . 1 7 C . 3 4 3 6 0 . 1 0 5 0 3 . 246 1 8 . 9 2 1 7 . 3 9 0 . 3 9 4 1 C . 1 2 1 0 3 . 6 8 3 1 9 . 7 1 1 7 . 7 0 0 . 4 5 1 7 0 . 1 4 0 0 4 . 2 2 8 2 0 . 5 5 1 7 . 9 9 C. 5 0 9 3 C . 1 6 2 8 4 . 7 9 0 21 . 2 0 1 8 . 2 4 C . 5 5 1 4 C . l 724 5 . C5C 2 2 . 1 2 1 8 . 60 0 . 6 1 1 6 0 . 1 1 9 0 5 . 500 2 3 . 0 1 1 8 . 8 8 C . 6 6 1 0 0 . 2 1 5 7 5 . 9 4 2 2 3 . 8 3 1 9 . 1 8 0 . 6 9 9 4 0 . 2 3 1 6 6 . 3 7 3 2 4 . 6 2 1 9 . 57 C . 7 4 3 2 0 . 2 4 6 4 6 . 7 8 4 2 5 . 5 7 1 9 . 9 9 C . 7 9 9 0 C . 2 6 4 3 7 . 257 2 6 . 4 2 2 0 . 4 8 0 . 8 4 9 0 0 . 2 8 1 3 7 . 7 5 2 2 7 . 3 9 2 1 . 0 2 C . 9 C 8 3 0 . 2 9 1 3 8 . 283 2 8 . 4 7 21 . 6 7 0 . 9 4 7 7 0 . 2 9 4 3 8 . 8 5 0 11) MEAN TEMPERATURE OF INNER C Y L I N D E R . (2 ) MEAN TEMPERATURE OF OUTER C Y L I N D E R . ( 3 ) CORRECTED I I . (4) CORRECTED 1 2 . (5) CORRECTED V I . 91 1ABLF: NO 21 VOLTAGE BETWEEN THE CYLINDERS= 10.80 KV ( 1) (2 ) (3 ) (4) ( 5 ) 17.95 16.97 C.3138 0.0968 2. 927 18.42 17.15 0.3574 0.1100 3. 340 19.09 17.40 C.41C9 0.1280 3. 82 8 19.77 17.67 0.4647 0.1490 4.287 20.45 17. 8 7 C.5073 0.1628 4.690 21.44 18.26 0.5692 : ; 0.1824 5. 230 22 .14 18.57 0.6074 '.' 0 .1980 5.580 23.10 18. 89 C.6618 0.2167 6.079 23.83 19. 18 0.7C64 C.23C6 6.472 24.49 1 9. 64 C.7512 0.2454 6. 863 25.34 20.08 0.7910 0.2513 7. 3 47 26.10 20.50 0.8284 0.2573 7.753 26.80 20.93 C. 8667 0.2 573 8.140 27.77 21.48 0.9158 C.2833 8.630 28.72 22. 04 C.9710 0.2983 9.058 (1) MEAN TEMPERATURE OF INNER CYLINDER. (2) F E AN TEMPERATURE OF OUTER \CYLINDER. 13) CORRECTED I I . . (4) CORRECTED 12. (5) CORRECTED VI.

Cite

Citation Scheme:

    

Usage Statistics

Country Views Downloads
United States 7 0
Japan 7 0
Ireland 6 0
Germany 5 0
China 1 0
City Views Downloads
Unknown 11 0
Tokyo 7 0
Ashburn 3 0
Sunnyvale 2 0
Mountain View 1 0
Beijing 1 0
Redwood City 1 0

{[{ mDataHeader[type] }]} {[{ month[type] }]} {[{ tData[type] }]}

Share

Share to:

Comment

Related Items