METHOD AND APPARATUS FOR THE MEASUREMENT OF THE SIZE OF PARTICLES ENTRAINED IN A GAS

Abstract: Apparatus for measuring the size of particles entrained in a gas comprises optical fibre probe means having a hemispheri­ cal end and insertable into a stream of gas bearing entrained particles to create a stagnation point in the flowing gas, a Dop­ pler anemometer to measure particle velocity at a point up­ stream of the stagnation point and computer-controlled burst counter to determine from the particle velocity measurements the distribution by number of aerodynamic diameter of the par­ ticles.

Description: This invention relates to methods and apparatus for the measurement of the size of particles entrained in a gas. An existing method is to sample the particles from the gas, for subsequent offline analysis. In addition to the difficulty of ensuring that samples are representative of the entrained particulate, it is then necessary to ensure that the particles measured are no more or less agglomerated than those in the gas. The uncertainties associated with these measurements introduce doubts into measurements of the apparent efficiency of gas cleaning devices. It is therefore desirable to measure the concentration and size distribution of the entrained particles directly in situ without sampling. It is also desirable for the measurement device to be capable of being traversed, so as to obtain representative mean values across a large gas duct. The "size" of a non-spherical particle can be defined in a variety of different ways. For gas cleaning operations, the parameter of most relevance is usually the Stokes diameter, which describes both the intertial and settling behaviour of a particle. It is possible, and sometimes preferable, to infer the Stokes diameter from simple off-line measurements such as are obtained by Coulter Counter or image analysis. The Coulter Counter, for example, measures directly the particle volume and, for wide ranges of particle shape, the volume-equivalent diameter is close to the Stokes diameter. However, these indirect measurements are always time-consuming, demand redispersion in a liquid, and require at least two different methods to be used on the same sample in order to infer the precise relationship between the volume-equivalent and Stokes diameters. Of the devices currently available for direct measurement of Stokes diameter, horizontal, vertical and centrifugal settling devices are of considerable interest for calibration of other aerodynamic instruments. However, they are essentially laboratory techniques which operate on batch samples and which require experience, care and patience to obtain reliable results. Inertial impaction is probably the commonest industrial technique for measurement of Stokes diameter. However, in addition to requiring the gas to be sampled, inertial impaction measurements are subject to uncertainties arising from the possibilities that particles may bounce or blow off from the impaction plates and that agglomerates may break up on impact. Cyclones with sharp cut-off have been used for particle size measurement. Recently an automated 3-cyclone train has been described using an oscillating microbalance to indicate the masses of particles collected. Although this device is automated, it remains a sampling instrument which gives readings in a few wide size ranges at intervals of tens of seconds. In Applied Physics Letters 25 (No.4) p 193 (15 February 1975) "Aerosol size spectrum analysis using relaxation time measure ment" by K.J. Kirsch and M.K. Mazumbder, there is descirbed an experimental method to measure the dynamic relaxation times of aerosol particles. The particle relaxation time is determined from the ratio of the velocity amplitude of an aerosol particle to the velocity amplitude of the medium while the particle is subjected to acoustic vibration of known frequency. The TSI "Aerodynamic Particle Sizer" or APS is another relatively new instrument which differs from the other devices in measuring the inertial behaviour of individual particles. In principle, it can give finer resolution of the particle size distribution than a device such as a cascade impactor which sorts particles into a few discrete bands. Dilute particle-laden gas is accelerated through a nozzle, so that the entrained particles reach a velocity significantly less than that of the gas. Thus, for each particle, the slip velocity relative to the gas depends upon its inertia, i.e. upon its aerodynamic diameter. The velocities of particles issuing from the nozzle were measured by laser-Doppler anemometry, but the commercial instrument measures the time-of-flight between two laser beams separated by 120 m. The particle size distribution is built up from measurements of individual particles. In theory this technique should be absolute, but in practice it is necessary to calibrate the instrument using particles of known aerodynamic diameter. While the APS appears to give reliable measurements for spherical particles, there is some evidence that it undersizes irregular particles. Particle orientation effects are insufficient to explain the discrepancy; A contributory cause may be shape features which induce boundary layer separation at relatively low Reynolds numbers and hence have a strong effect on drag. The APS also has the disadvantage that it is another sampling instrument, with a low sample gas rate which would have to be cooled before being passed to the instrument. Like all single-particle counters, it is subject to "coincidence errors", arising from particles whose times-of-flight overlap, so that most gases of industrial interest also have to be diluted before analysis. We have devised an alternative apparatus and method in order to avoid the need to draw a sample. Flowing gas approaches an obstacle in the duct. The gas velocity on the stagnation streamline, u, drops, to zero at the stagnation point. The velocity of a particle on the stagnation streamline, v, also falls but, because of the inertial of the particle, it will not follow u. The extent of the difference between v and u increases with increasing particle inertia. Thus the principle on which the present invention is based is to insert a probe into the gas flow, and to measure the velocities of individuual particles at a known distance from the stagnation point. According to the present invention there is provided apparatus for measuring the size of particles entrained in a gas including probe means insertable into a stream of gas bearing entrained particles to create a stagnation point in said flowing gas, measuring means to measure particle velocity on the stagnation streamline at a point at a measured distance from said stagnation point upstream of said stagnation point and computer means to determine from said particle velocity measurements the distribution by number of aerodynamic diameter of said particles. An embodiment of the invention will now be described by way of example with reference to the accompanying drawings in which Figure 1 illustrates gas and particle motion approaching stagnation point of a bluff body. Figure 2 is a schematic view of apparatus for measuring particle velocity in accordance with an embodiment of the invention. Figure 3 is a graphical representation of the predicted velocities of particle entrained in a gas approaching a sphere. Figure 4 is a correlation chart. Figures 5 a and 5b are doppler signals obtained from glass fibres on a loudspeaker. Figure 6 is a diagram of apparatus used for in-duct particle size measurement; and Figures 7a to 7c are Doppler traces obtained using the apparatus of Figure 6. In this specification, means ascribed to symbols are: a Radius of hemispherical probe C D Fluid-particle drag coefficient Cs Slip correction factor d Particle diameter d a Aerodynamic particle diameter d e Volume-equivalent particle diameter d s Stokes diameter of entrained particle F b Body force on particle F d Fluid-particle interaction force. St Stokes number, defined by equation (2) t Time t′ Dimensionless time u, Local gas velocity u′ Dimensionless gas velocity u Approach gas velocity V Volume of particle v, Particle velocity x Coordinate of particle centre x′ Dimensionless particle coordinate µ Gas viscosity ρ Gas density ρ p Particle density ρ p ′ "Standard" particle density Consider a rigid particle of volume V moving at instantaneous velocity v through a fluid (Figure 1). In general, the equation of motion of the particle takes the form ρ p V dv dt = F d + F b where ρ p is the particle density and F d and F b are the forces acting on the particle, F d being the instantaneous force exerted on the particle by the fluid and F b representing body forces due, for example, to gravity or an electric field. The fluid/particle interaction term F d includes, in general, the force arising irom pressure gradients in the fluid and also the drag of the fluid on the particle which depends, in turn, on the velocity of the particle relative to the surrounding fluid and in its present and past acceleration relative to the fluid. For particles in gases, it is usually sufficiently accurate to approximate Fd by the force arising in steady relative motion at the instantaneous slip velocity, i.e. Fd = nρd²C D (u-v) |u-v| 8C s where u is the gas velocity in the vicinity of the particle, so that (u-v) is the instantaneous slip velocity, d is the particle diameter discussed below, and ρ is the gas density. The dimensionless coefficient C s in equation (2) is the Cunningham "slip correction" factor which allows for the reduction in drag when the particle diameter is of comparable magnitude to the mean free path of molecules in the gas. For particles greater than 1µm in size in gases, it is close to unity. The dimensionless drag coefficient C D is in general a function of the particle Reynolds number: Re p = |u-v|ρd/µ where µ is gas viscosity. Various empirical and semi-empirical relationships for C D (Re p ) are available. For Re p < 0.5, the relationship can be approximated by Stokes' law: Fd = 3 n µd (u - v)/C s The body force term, F b , can be neglected for small particles in rapidly accelerating or decelerating gases, and for flows in which gravity is the only body force, u is horizontal, and the linear drag relationship of equation (4) applies. All the above simplifications can be applied, at least for an initial analysis, to the in-duct particle sizer. Equation (1) then becomes ρ p V dv dt = 3 n µd(u-v) or ρ p d e ³ 18 d dv dt = u-v where d e is the volume-equivalent diameter of the particle, given by V = πd e ³/d Equation (5) shows that d is to be interpreted as the diameter of the sphere which experiences the same fluid drag as the particle. Equation (6) shows that the Stokes diameter of the particle, d s , i.e. the diameter of the sphere with the same inertial behaviour, is given by d s ² = d e ³/d The Stokes diameter also corresponds to the diameter of the sphere with the same settling velocity. The aerodynamic diameter d a is usually defined as the diameter of a sphere of arbitrary standard density ρ p ′ with the same settling or inertial behaviour as the particle in question. From equation (6), d a = d s ρ p /ρ p ′ In terms of the Stokes diameter, equation (6) becomes where x is the position of the particle. For one-dimensional motion, as on a straight stagnation stream line, x and u can be written as scalars. The resulting equation can be written in dimensionless form by substituting: x′ = x/a; u′ = u/u ∞ ; t′ = u ∞ t/a where a is the radius of the axi-symmetric body on which the stagnation point occurs, and u ∞ is the velocity of the approach gas flow at large distances from the body. Equation (10) then becomes where St is the Stokes number of the particle: St = d s ² Equation (12) is the dimensionless form of the simplified equation of motion of the particle. Its solution depends on u′, which describes the gas velocity approaching the stagnation point and therefore depends on the shape and Reynolds number of the body, 2u ∞ pa/µ. For the in-duct particle sizer, it is convenient to make the upstream surface of the probe in the form of a hemisphere. Provided that the probe Reynolds number is large, the flow around this surface can be represented by flow of an incompressible invicid fluid around the upstream half of a sphere. The fluid velocity on the stagnation streamline is then given by u′ = x′⁻³ - 1 so that equation (12) becomes Equation (15) can be solved numerically, for different values of St, given an appropriate initial condition. The true initial condition is t′=0, x′=∞, dx′/dt′=-1. In equation (12) x′ is the dimensionless position of the particle, t′ is a dimensionless time, and u′ is a dimensionless velocity: x′=x/a ; u′ = u ∞ /u ; t′ = u ∞ t/a where u is the gas velocity at position x on the stagnation streamline. The theoretical response of the in-duct particle sizer can be calculated from equation (1) if u′ is known as a function of x′. This requires the upstream face of the probe to be designed to a shape for which a theoretical expression for u′ is available. For this reason, and also for relative ease of fabrication, the probe has been designed to have a hemispherical end, as shown schematically in Figure 2. An optical fibre 2 is mounted in a pitot channel 3 and is equipped with a lens 4 focused upstream of the probe end 1. A Doppler anemometer 5 is coupled to the optical fibre 2 and feeds a burst counter 6 which is, in turn monitored by a microcomputer to determine the distribution by number of the aerodynamic diameter of the particles. The appropriate value for a is then the radius of the hemispherical end and cylindrical probe stem, and x is measured from the centre of the hemisphere. Using the appropriate expression for u′ the dimensionless particle velocity dx′/dt′ = v/u ∞ can be calculated, where v = dx/dt is the dimensional particle velocity. The results are shown in dimensionless form in Figure 3, each curve referring to one value of the particle Stokes number. Figure 3 can now be used to select an appropriate probe size, for the approach gas velocity and range of particle sizes of interest. Consider, for example, a gas velocity, u ∞ , of 10 ms⁻¹ with particles from 1 to 20µm to be measured. If the probe radius, a, is designed to be 5mm, then this range of particle size covers the range of Stokes numbers in Figure 3. Because the approach velocities are well discriminated, it should be possible to distinguish between sizes in this range by measuring their velocities approaching the stagnation point. It remains to select the position at which the particle velocities are best sensed. Figure 4 shows the same predictions plotted to give particle diameter as a function of v/u ∞ for several measuring positions, x/a. This "calibration chart" clearly shows that maximum resolution should be achieved by measuring particle velocity as close as possible to the surface of the probe, subject to the constraint that the measuring point must be outside the boundary layer. Figure 4 shows that a suitable measurement position is x/a = 1.1, i.e. 0.5mm from the surface for a probe tip radius of 5mm. Furthermore, measurement of velocity with a single probe size cannot discriminate particles over the whole size range, because particles with St greater that about 6 (d s = 20µm for these conditions) have such high inertia that their velocity changes little, while for St less than about 0.06 (ds≈2µm in this case) the particle velocity is barely distinguishable from that of the gas. Thus the measuring range of the device is limited to rough a ten-fold range of particle diameter. However, the range can readily be changed by changing a; i.e. by changing the probe diameter or by using several hemispherical ends of different diameters. Conveniently measurements may be made using optic fibres with illumination by a standard 4mW He-Ne laser which generates polarised light at a wavelength of 633nm and at a power of 2mW launched into the fibre. The signals detected can be collected and displayed for inspection by a Gould 4035 digital storage oscilloscope, with suitable traces output off-line to a conventional pen recorder. Alternative optic fibre systems may be used. In order to compare the performance of different fibre and lens systems, the preliminary tests used as "targets" a number of glass fibres of different sizes mounted on the cone of a 100mm diameter loudspeaker driven at 1kHz with amplitude of order 5mm. This controlled sinusoidal motion of the target should produce Doppler signal bursts related to the driving waveform. A micro-positioner was used to locate the end of the optic fibre precisely relative to the end of the target fibre. One form of fibre optic Doppler anemometer uses a graded-index glass optical fibre of 100µm nominal diameter. Figure 5a shows traces obtained from the end of a glass fibre of 10µm diameter using the standard optic fibre with no lens system at the detection end. The upper trace shows the displacement of the loudspeaker; thus the Doppler signals have the expected form, with high frequency when the target velocity is high. However, some high-frequency noise is apparent, especially in the lower velocity sections. Other systems were therefore tested in an attempt to increase signal strength and reduce noise. The signal was found to be improved significantly by reducing the optical fibre diameter, possibly because multimode fibres were used rather than monomode fibres and reducing the fibre diameter reduces the number of modes by which the fibre transmits the light. Adding a lens to the end of the fibre focusses the illuminating light into a small volume from which scattered light can be received. It was found that, when lenses were used, the signal was still further improved especially when detecting small diameter targets, probably because the focused scattering volume is smaller. The best results were obtained with a SELFOC graded-index cylindrical lens. Figure 5b shows signals obtained for the 10µm target with a SELFOC lens on the end of a 50µm graded-index glass fibre. The signal strength is seen to be greater. It is strongest when the target is at the focus of the laser light, about 1mm from the end of the optic head, and weaker at other positions. Qualitatively, this is the type of response required for the in-duct particle sizer. Once it was established that Doppler bursts suitable for further processing can be obtained from a moving target at the required position, further tests were carried out for dispersed particles using the fibre and lens in Figure 5b. The optic head, without the rest of the probe, was mounted as shown schematically in Figure 6, facing upstream in a horizontal duct. Air velocities of 0.5, 1.0 and 1.5 m.s.-¹ were used as the approach flow for these tests, because higher velocities would give rise to Doppler frequencies too high to be stored by the particular digital oscilloscope available. The test particles used were 19.1µm latex spheres, dispersed into the air flow using a two-fluid atomiser using high pressure air as the driving fluid. No particular measures were taken to ensure that the latex was dispersed as individual particles. Approximating the end of the optic head as a hemisphere of 3mm radius, the values of St for the 19.1µm primary particles at the three air velocities are 0.19, 0.38 and 0.58. Referring to Figure 3, the tests are therefore in the range where slip between the gas and particles should be measurable and dependent on particle size. Figures 7b, c and d show typical Doppler "bursts" obtained from entrained particles at these three velocities; clearly these have exactly the form expected. In some traces it is possible to detect a reduction of the Doppler frequency, i.e. deceleration of the particle, through the "burst" as the particle approaches the stagnation point on the optic head. The signals obtained at each air velocity indicate a range of particle velocities, the lowest being close to the value predicted from Figure 3. Possibly the higher velocities correspond to incompletely dispersed particles, but this remains to be confirmed. Preliminary tests with smaller latex particles confirmed that lower approach velocities were measured.

Claims: 1. Apparatus for measuring the size of particles entrained in a gas including probe means (1,2,4) insertable into a stream of gas bearing entrained particles to create a stagnation point in said flowing gas, measuring means (5) to measure particle velocity on the stagnation streamline at a point at a measured distance from said stagnation point upstream of said stagnation point and computer means (6,7) to determine from said particle velocity measurements the distribution by number of aerodynamic diameter of said particles. 2. Apparatus for measuring the size of particles entrained in a gas as claimed in claim 1 characterised in that said probe means comprises an optical fibre (2). 3. Apparatus for measuring the size of particles entrained in a gas as claimed in claim 2 characterised in that said optical fibre is equiped with a lens (4) focused upstream of the probe means. 4. Apparatus for measuring the size of particles entrained in a gas as claimed in claim characterised in that a Doppler anemometer (5) is coupled to said optical fibre. 5. Apparatus for measuring the size of particles entrained in a gas as claimed in any one of the preceding claims and having demountable probe means (1,2,4) to measure different sizes of particle. 6. Apparatus for measuring the size of particles entrained in a gas as claimed in claim 2 characterised in that said optical fibre (2) is a graded index optical fibre. 7. A method of measuring the size of particles entrained in a gas comprising inserting probe means into a gas stream bearing entrained particles to create a stagnation point, measuring the particle velocity on the stagnation streamline at a point at a measured distance from said satgnation point upstream of said stagnation point and computing from said particle measurments the distribution by number of the aerodynamic diameter of said particles.

Inventors: CLIFT, ROLAND, SEVILLE, JONATHAN PETER KYLE, TATE, ANDREW HENRY JOHN
Applicants: NATIONAL RESEARCH DEVELOPMENT CORPORATION
Application year: 1987
Main classification: G01N15
IPC codes: G01N15