DROPLET SIZE MEASURING APPARATUS

Abstract: In the droplet size measuring apparatus of the forward scattering method, particles (6) to be measured in a scattering zone (A) having a predetermined scattering path length (L) are irradiated with parallel laser beams. The irradiated beams are scattered by the particles (6) and the scattering beams are detected by the photodetectors (8-1∼8-K) arranged at predetermined scattering angles (ϑ₁∼ϑ K ) The output of the photodetectors (8-1∼8-K) are input to the arithmetic operat­ ing unit (11), thus the intensity distribution I(ϑ) of the scat­ tered beams is measured. The droplet size distribution n(D) is calculated from the intensity distribution I(ϑ) using the following equation: I(ϑ)=∫{I₀exp(-ℓ∫c(D)n(D)dD·∫i(D,ϑ+ϑℓ)n(D)d(D)}dℓ

Description: This invention relates to an apparatus for measuring the droplet size distribution of very small particles. This type of an apparatus is now being developed which measures the droplet size distribution of, for example, droplets of vapor in a turbine or droplets of fuel jetted into a combustion chamber, through the utilization of the Mie scattering. By the droplet size distribution is meant, here, a relative distribution in the number of particles for a respective droplet size or diameter, i.e., a probability density function. According to the "Mie scattering" theory, the scattered light intensity i(D,ϑ) in a certain angular direction ϑ (at a scattering angle) when very monochromatic light beams, such as laser light, strike a spherical particle with a diameter D can be calculated. The scattered light intensity I(ϑ) of a group of spherical particles of various sizes can be expressed as a sum of respective scattered beams i(D,ϑ) on one particle of the diameter D, as set out below. I(ϑ) = ∫i(D,ϑ)n(D)dD (1) where n(D) denotes the droplet size distribution. A method for measuring the droplet size distribution with the use of the principle of Equation (1) is called as the forward scattering method. This method is one type of the indirect measuring method which comprises: (1) a first step of irradiating a group of particles with monochromatic parallel laser beams and measuring the angular distribution of the irradiation beams scattered on the group of particles -- hereinafter referred to simply as the intensity distribution of the scattered beams; and (2) a second step of evaluating the droplet size distribution n(D) from the intensity distribution of the measured scattered beams I(ϑ) with the use of Equation (1). Fig. 5 shows a conventional apparatus for measuring the droplet size by the aforementioned method. In this apparatus, output light 60 from a light source (not known), such as a laser is conducted via illumination optical fiber 3 into a scattering zone A and directed at the group of particles 6 in the scattering zone A. Beams 61 scattered at the respective scattering angles are conducted to photodetectors (not shown) via optical fibers 7a, 7b, 7c and 7d to measure the intensity distribution I(ϑ) of the scattered light. These optical fibers 7a, 7b, 7c, and 7d are arranged equidistant from the center P of the scattering zone A and for predetermined scattering angles. In this apparatus, however, on account of the optical fibers 7a, 7b, 7c, and 7d arranged at the predetermined scattering angles the optical fiber 7c, for example, receives not only the beams scattered at the scattering angle ϑ f but also beams scattered on the particles occupied in a broader range as indicated by broken lines 63 in Fig. 5. As a result, the photodiode connected to the optical fiber 7c will receive the scattered beams at other than the scattering angle ϑ f as well. Stated in a stricter way, ϑ in Equation (1) fails to indicate the set angles of optical fibers 7a, 7b, 7c, and 7d. However, no problem arises in the cases where the scattering path length L in the scattering zone A is shorter. The scattering path length L may be made shorter in the case of a larger droplet size involved since the intensity of the scattering light is adequate. For the smaller size of particles to be measured, on the other hand, the intensity of the scattering light is very small and, in order to increase the number of particles in the scattering zone A, it is necessary to increase the scattering path length L, failing to express the intensity distribution of the scattering light with Equation (1). EP-A-0 167 272 discloses similar laser apparatus which calculates a relative particle size distribution of particles to be measured from the scattered light intensity pattern detected by a detecting section, and the relevant features are acknowledged in the preamble to claim 1 below. An object of this invention is to provide a droplet size measuring apparatus of the forward scattering method which can accurately measure the droplet size distribution even for the smaller droplet size of particles to be measured and the longer scattering path length. Another object of this invention is to provide a droplet size measuring apparatus which can measure the absolute droplet size distribution representing the distribution of absolute number of particles, as opposed to the relative droplet size distribution of particles. According to this invention there is provided a droplet size measuring apparatus as defined in claim 1 below. This invention can be more fully understood from the following detailed description of a preferred embodiment, by way of example only, when taken in conjunction with the accompanying drawings, in which: Fig. 1 is a view showing a droplet size measuring apparatus according to one embodiment of this invention; Fig. 2 is an explanatory view for explaining the principle on which the droplet size distribution is calculated; Fig. 3 is a flowchart showing the steps of calculating the droplet size distribution on the apparatus according to this invention; Figs. 4A and 4B are views for restricting the scattering path length in the scattering zone by the embodiment of this invention; and Fig. 5 is a view showing a conventional droplet size measuring apparatus. Fig. 1 show the whole arrangement of a droplet size measuring apparatus according to the embodiment of this invention. Output light from laser unit 1 (light source) is incident through lens 2 on optical fiber 3 for illumination. The laser light transmitted through optical fiber 3 emerges from the end of optical fiber 3 and is sent to collimator lens 4 where it illuminates, as parallel beam 13, a scattering zone A. cylindrical spacer 5 is located between collimator lens 4 and the scattering zone A and adapted to be moved in the axial direction of parallel beam 13, when screws 17 are loosened, so that the scattering path length L is adjusted. Spacer 5 has hole 12 somewhat greater in diameter than the diameter of parallel beam 13. This is because parallel laser beam 13 is scattered and the scattering beams become a background noise in the measurement of the intensity distribution of the scattered beams by the particles 6 in group. To this end, the hole is formed in spacer 5 to prevent passage of the scattering light therethrough. The scattering path length L in the scattering zone A where the group of particles 6 to be measured exists is restricted for the reason as will be set forth below. The restriction of the scattering path length L is made by, for example, confining the scattering path length L by glass plates 70 as shown in Fig. 4A or forming window 16 in part of pipe 15 to permit passage of parallel beam (illumination beam) 13 therethrough. The beams which are scattered on the particles 6 present within the scattering zone A are received by receiving optical fibers 7-i (i=0∼K,K: an arbitrary positive integer) located at predetermined scattering angles 0, ϑ₁, ..., ϑ K where 0 denotes the transmission of light without being scattered. Then the received beams enter photodiodes 3-i. Optical fibers 7-0 is adapted to measure transmission light. Here it is to be noted that the scattering angles ϑ₁, ..., ϑ K represent angles viewed with the center P of the scattering zone A as a reference and that the light incident ends of respective optical fibers 7-i are equidistant from the center P of the scattering zone A. In this specification, we define the following quantities: ϑ : Scattering angle defined with respect to the centre of the scattering zone; ϑ+ϑ ℓ : Actual scattering angle at scattering position ℓ. The light beams, after being converted by respective photodiodes 8-i to electric currents, are supplied to respective amplifiers 9-i where they are amplified as voltage signals. The voltage signals of the respective amplifiers 8-i are supplied through a multiplexer 10 to arithmetic operating unit 11 and the result of the arithmetic operation is displayed on display 12. Keyboard 13 for inputting various data items is connected to arithmetic operation unit 11. Reference numeral 14 is a probe for including a measuring system of the intensity distribution of the scattered beam. An equation for providing the intensity of light received by optical fiber 7-f will be obtained below with reference to Fig. 2, noting that optical fiber 7-f is arranged at the scattering angle ϑ f . Given a consideration to the decrease of the scattering beams by particles, the intensity I b (ℓ) of the laser beams at the scattering position ℓ (ℓ = 0∼L) is given by: I b (ℓ) = I₀ exp(-ℓ∫c(D)n(D)dD) (2) where I₀: the intensity of the beam at ℓ = 0 c(D): the scattering cross section of a particle with a diameter D n(D) : Probability density for the diameter D times the particle number density The term of exp(-ℓ∫c(D)n(D)dD) represents the decrease of the scattering beams. The intensity ΔI(ℓ) of light scattered on particles which exist between the scattering positions ℓ and ℓ+Δℓ can be expressed as follows: ΔI(ℓ) = I b (ℓ)∫i(D,ϑ f +ϑ ℓ )n(D)dD (3) With reference to the explanation on page 1 of the scattered means i(D,ϑ), we define the function i(D,ϑ ℓ ) as follows: i(D,ϑ+ϑ ℓ ) : Function calculated within the Mie scattering theory for a particle with diameter D, and defined by Δ I(ℓ) = I b (ℓ)∫i(D,ϑ+ϑ ℓ )n(D)dD , where ΔI(ℓ) is the differential scattered intensity contributed to I(ϑ) by particles between ℓ and ℓ + Δℓ, and where I b (ℓ) is the incoming beam intensity at scattering position ℓ. Substituting Equation (2) into Equation (3) yields: ΔI(ℓ) = I₀ exp(-ℓ∫c(D)n(D)dD·∫i(D,ϑ f +ϑ ℓ )n(D)dD (4) With the scattering path length L restricted, the intensity I(ϑ f ) of the scattered light which is received by optical fiber 7-f will be in which case optical fiber 7-f is arranged at the scattering angle ϑ f . When Equation (4) is substituted into Equation (5), then showing the intensity of the scattered light which is received by optical fiber 7-f arranged at the scattering angle ϑ f . Arithmetic operation unit 10 solves Equation (6) for n(D) with the measured value as I(ϑ f ) and evaluates the droplet size distribution. Here the solution of Equation (6) will be set out in more detail. Equation (6) is solved by the following two methods: (1) assuming the distribution function of n(D) or (2) dividing the integration interval of the integrating equation in Equation (6) and letting Equation (6) approximate through the use of the simultaneous equations. The method (1) will be explained below, assuming the distribution function of n(D). With one of the logarithmic-normal distribution, normal distribution, Rossin-Rammler, and the like properly selected as the distribution function of n(D), then n(D) = N₀n r (D) (7) where N₀: the particle density n r (D): the distribution function selected With the logarithmic-normal distribution selected, where D: the droplet size π: the ratio of the circumference of a circle to its diameter A,B: the distribution parameters I(ϑ) is evaluated with respect to a plurality of ϑ s . Substituting Equation (8) into Equation (7), and Equation (7) into Equation (6), results in determining the N₀, A and B through the use of the least squares method. In this case, n(D) is evaluated from Equation (7). The method (2) will now be explained below. With the integration interval of the particle size D as [D min , D max ], Equation (6) becomes Rearranging the order of the integration with respect to ℓ and D, With [D min , D max ] N-divided, Equation(10) becomes Here, D₀ = D min and D N = D max . If it is assumed that n(D) is constant within the subinterval [D j-1 , D j ] and n(D)dD is replaced by n(Dj), then the right term of Equation (11) becomes If part of the right term of Equation (12) is put in the form of then Equation (11) will be With the measuring point ϑ f of the intensity of the scattering light as ϑ₁, ϑ₂, ..., ϑ K , Here Equation (17) is represented as G(n) since it is the function of n. Equation (15) can be expressed by the following simultaneous equations with n as an unknown number as in the case of Equations (16) to (18) I = G(n)·n (19) noting that Equation (19) can be solved based on the Newton's method. With G(n)* put as the transverse matrix of G(n), the least squares solution satisfies G(n)*I-G(n)*G(n)·n = 0 (20) If G(n) is regarded as being a constant in the neighborhood of n₀ (an initial value of n), then the approximation of Equation (20) can be made by the linear equation. G(n₀)*I-G(n₀)*G(n₀)·n = 0 (21) Equation (21) is solved as: n = (G(n₀)*G(n₀))⁻¹G(n₀)*I (22) The nonlinear equation of Equation (20) can be solved for n if the solution n of Equation (22) is put as n₀ and the calculation of Equation (22) is repeated until the absolute value of a difference between n and n₀ satisfies the following relation: | n₀ - n | < ε (23) It is, therefore, possible to evaluate the distribution n(D j ) of the droplet sizes per the subinterval [D j-1 , D j ] of the particle sizes. In Equation (23), ε denotes a very small number for convergence judgment. The approximation of G(n), which may be regarded as being constant for n in the neighborhood of n₀ will be given below. Now an arbitrary element i s (D g , ϑ p , n(D)) of G(n) will be considered, noting that i s is given by Equation (13). Rearranging the order of the integration of Equation (13), Here Substituting this approximation (Equation (25)) into Equation (24) gives: Differentiating Equation (26) with n g = n(D g ) it follows that where c g = c(D g ). In Equation (27), the near maximal values are taken as L, c g and ΔD given below: L = 10⁻² (m) c g = 6.3 × 10⁻⁸ m² (the scattering cross section of polystyrene particles in water which are illuminated with He-Ne laser beam of 200 µm) ΔD = 20 × 10⁻⁶ m In this case, Lc g ΔD = 1.3 × 10⁻¹⁴ (28) Therefore It may be assumed that i s (D j ,ϑ f ,n(D)) is constant in the neighborhood of n₀. Obtaining the initial value n₀ on the bases of Equation (22) where yields a comparatively better approximation value. The algorithm of this method (2) for obtaining the approximation of Equation (6) with the use of the simultaneous equations will be explained below with reference to the flowchart of Fig. 3. At step S10 the intensity distribution I(ϑ i ) of the scattered beam obtained at photodiode 8-i is input to arithmetic operating unit 11, noting that ϑ i = ϑ₁ ∼ ϑ K At step S20, the matrix element of the initial value G₀ = G(0) is calculated with the use of Equation (13). At step S30 the initial value n₀ of n is calculated through the substitution of G₀ into Equation (22). At step S40, n is calculated from Equation (22) with the use of n₀. At step S50, it is determined whether or not the absolute value | n₀-n | of a difference between n and n₀ is smaller than ε. If the absolute value of | n₀-n | is not smaller than ε, then the process goes to step S60 where n₀←n and step S60 goes back to step S40 where n is again evaluated from Equation (22). If at step S50 it is determined whether the absolute value of | n₀-n | is smaller than ε, then at step S70 n(D₁)∼n(D N ) is displayed. If in Equation (22) no solution is obtained because of a greater difference between the maximum value and the minimum value of the eigenvalues of G(n₀)*G(n), then use may be made of the constrained inversion method. This method is disclosed in detail in "The 1983 Dissertations Vol. 19, No. 10" of "Society of Instrument and Control Engineers" (Japan), the paper by K. Tatsuno Entitled "The Conversion of the Intensity Distribution of Scattering Light to the Droplet Size Distribution on a Droplet Size Measuring Apparatus using Laser Light". The paper discloses the solution of the least squares of Equation (19) which is found based on the first condition that the sum of the differences of third order for the droplet size distribution n₁, n₂, ..., n N is constant and on the second condition that n₁, n₂, ..., n N is constant and on the second condition that n₁, n₂, ..., n N are positive or zero. That is, with the solution of the aforementioned first condition put as and with γ denoting an undetermined multiplier, the solution of Equation (19) is given by n = (G(n)*G(n)+γH)⁻¹G(n)I (32) through the Lagrange's method of undetermined multiplier. Here With the geometric mean of the eigenvalues of G(n)*G(n) given as γ of Equation (32), a proper constraint can be given. It is necessary that, in order for n evaluated from Equation (32) to satisfy the second condition, n is meet the following condition. 〈〈 Condition 〉〉 y = G(n)*(G(n)×n-I) + γ×H×n Thus y j ≧ 0 for j at n j = 0 and y j = 0 for j at n j > 0 In this procedure the intensity distribution I(ϑ f ) of the scattered beams has to be found in terms of the light energy (absolute intensity distribution) as set forth below. For the particles with the absolute droplet size distribution n p (D) as represented by their individual number, such as polystyrene particles and glass spheres, the absolute intensity distribution I p (ϑ f ) of those scattering light, that is, the intensity distribution of the scattering light expressed in terms of the energy level of the light, can be calculated with the use of Equation (6). The intensity distribution of the scattering light on those particles is measured. With the measured value as I pm (ϑ f ), the conversion table T(ϑ f ) will be If the measured value of the intensity distribution of the scattered light on the to-be-measured particles with their absolute droplet size distribution unknown is given by I m (ϑ f ), then the absolute intensity distribution I(ϑ f ) of the scattered light can be found from Equation (34) as follows: I(ϑ f ) = T(ϑ f )·I m (ϑ f ) (35) In this way, the absolute calibration is made for the system of measuring the intensity distribution of the scattered light and I(ϑ f ) can be found as the light energy level, thus obtaining the absolute particle size distribution as their individual number for n(D). According to this invention, owing to the restricted scattering path length in the scattering zone, the scattered beams of the respective scattering angles are obtained as the integration values of the scattered beams for the respective scattering positions in the scattering zone, while considering the decrease of the scattered light beams in the scattering zone, so that it is possible to calculate an accurate scattering light intensity distribution and thus an accurate droplet size distribution on the apparatus of this invention.

Claims: A droplet size measuring apparatus comprising: light source means (1, 2, 3, 4, 5) for directing parallel light beams onto a group of particles to be measured in a predetermined scattering path length (L); measuring means (7-0∼7-K, 8-0∼ 8-K), having photodetectors arranged for respective scattering angles, for measuring an intensity distribution of light beams scattered on the group of particles; and calculating means (11) for calculating the droplet size distribution of the group of particles in accordance with the measured intensity distribution of the scattered light beams, characterised in that each photodetector is arranged to measure an intensity distribution of the light beams scattered on the group of particles over the full scattering path length (L); characterised further in that said calculating means (11) calculates the droplet size distribution in accordance with a decrease of the parallel beams within the scattering zone and the scattering path length (L) of the scattering zone, in addition to the measured intensity distribution, and characterised still further in that said calculating means (11) calculates the droplet size distribution n(D) of said group of particles from the intensity distribution I (ϑ) of said measured scattered light beams on the basis of an expression where I₀: an intensity of the parallel beams incident to the scattering zone; ℓ: a scattering position measured from the incident end of the scattering zone; c(D): the known scattering cross section of a particle with a diameter D; ϑ : Scattering angle defined with respect to the centre of the scattering zone ϑ+ϑ ℓ : Actual scattering angle at scattering position ℓ. n(D) : Probability density for the diameter D times the particle number density i(D,ϑ+ϑ ℓ ) : Function calculated within the Mie scattering theory for a particle with diameter D, and defined by Δ I(ℓ) = I b (ℓ)∫i(D,ϑ+ϑ ℓ )n(D)dD , where ΔI(ℓ) is the differential scattered intensity contributed to I(ϑ) by particles between ℓ and ℓ + Δℓ, and where I b (ℓ) is the incoming beam intensity at scattering position ℓ. The droplet size measuring apparatus according to claim 1, characterized in that said calculating means (11) calculates n(D) by assuming a distribution function for n(D) having unknown parameters, substituting the assumed distribution function into said Equation (1), and determining the unknown parameters with the use of the least squares method and Newton Raphson's method. The droplet size measuring apparatus according to claim 1, characterized in that said calculating means (11) calculates n(D) by dividing an integration interval for the particle size D of Equation (1) into subintervals, regarding the Equation (1) as simultaneous equations on the assumption that n(D) within the subinterval is constant, and solving the simultaneous equations with the use of the least squares method and Newton Raphson's method.

Inventors: TATSUNO, KYOICHI C/O PATENT DIVISION
Applicants: KABUSHIKI KAISHA TOSHIBA
Application year: 1987
Main classification: G01N15
IPC codes: G01N15