Droplet Characterization Based on the Simulated Secondary Rainbows ()
1. Introduction
Liquid atomization and spray play an important role in industrial production, fuel atomization and combustion, spray cooling, and spray drying, etc. [1] [2] [3] [4]. It is of great significance to accurately measure the size, refractive index (temperature), speed, and other parameters of the droplets during the atomization process. Due to inherent non-intrusiveness, the optical technology is used for particle/droplet characterization, including phase Doppler technology, particle image/tracking velocimetry, laser digital holography technology, interferometric particle imaging, time-shift technique, and rainbow technique (or called rainbow refractrometry) [2] [5] - [10].
Roth [10] first proposed the rainbow technique method, which is called the standard rainbow technique (SRT). The rainbow pattern was used to measure the refractive index and size of a single spherical droplet. Han [11] measured the small change in the diameter of the liquid column based on the rainbow technique in 1998, assuming that the refractive index is constant. Van Beeck [12] generalized the SRT to global rainbow technology (GRT) to measure the size distribution and temperature of the droplet cloud in the spray. Van Beeck and his colleagues also studied the inversion algorithm of GRT [13] [14]. GRT can be used to measure the mixing ratio of liquid-liquid suspension systems, droplets in large containment vessels, and multiphase flow [14] [15] [16]. In order to have a higher spatial resolution, Gaussian beams are used as a light source. The rainbow pattern of droplets illuminated by Gaussian beams was simulated and used for droplet characterization [17] [18]. It can achieve the purpose of measuring a single tiny droplet under high-density spray conditions. However, only the angles of the first two peaks are used for droplet inversion. This study will discuss the feasibility of other inversion algorithms to measure droplets based on the intensity of scattered light in the secondary rainbow region when Gaussian beams are used as the light source.
The structure of this paper is as follows: In Section 2, inversion schemes for the characterization of droplets are given. In Section 3, based on the secondary rainbow, two inversion schemes are used to compute the refractive index and size of the spherical droplet. Section 4 is devoted to the conclusions.
2. Inversion Algorithm of Rainbow Technology
Based on the geometrical optics, the deflection angle θp for a ray emerging from a spherical droplet is given as [19] [20]:
(1)
here θi is the incident angle, and m is the relative refractive index of the droplet. The physical meaning of p [20] is p = 0 is the light ray directly reflected by the spherical droplet, p = 1 is the light ray directly refracted by the spherical droplet, and p ≥ 2 is the light ray refracted by the spherical droplet, refracted after p − 1 times of internal surface reflection. Therefore, the deflection angle of the secondary rainbow (p = 3) is:
(2)
In the above formula, there is an incident angle that minimizes the deflection angle. The minimum deflection angle is called the geometrical optics rainbow angle, and the geometrical optics rainbow angle is given:
(3)
In Airy theory, the light intensity distribution in the rainbow angle region is expressed by Airy integral. Airy integral F(z) [20] is:
(4)
where z is a dimensionless parameter, defined as:
(5)
where,
(6)
Among them, the value of q is −1 or +1, depending on the angle of the emitted light relative to the incident light [20]. From Airy function, the related z values for the maximum, minimum and inflection point of the curve can be obtained. The normalized light intensity distribution shown in Figure 1 is calculated by using Debye series for spherical water droplet with a radius of 50 μm and a relative refractive index of 1.333. Here, the droplet is illuminated by Gaussian beam, the beam waist radius ω0 = 100 μm, and the wavelength is 0.6328 μm. The eigenvalues of the light intensity distribution in the rainbow area (such as θmax1, θmax2) have a correspondence with the dimensionless parameter z.
Performing algebraic operations on the above formula, the geometrical optics rainbow angle can be obtained from the first two troughs (θmin1 and θmin2):
(7)
The refractive index can be obtained from formula (3) and (7), and the droplet radius can be calculated:
(8)
here,
(9)
![]()
Figure 1. Normalized intensity distribution of the secondary rainbow region.
For the another algorithm, the angles of the first two inflection points (i.e. θinf1 and θinf2) of the intensity distribution are employed. The geometrical optics rainbow angle and droplet size are given by:
(10)
and,
(11)
here,
(12)
3. Inversion of Refractive Index and Size of Droplets
Firstly, according to the Debye series expansion, the scattered light intensity distribution for the secondary rainbow of the spherical droplet is calculated, and the influence of different Gaussian beam waist radius on the peak of the secondary rainbow light intensity distribution is explored. Take the droplet with refractive index of 1.333 and radius of 50 μm as an example, as shown in Figure 2. It can be found that when the beam waist radius is 50 μm, the peak and valley values tend to be stable. Therefore, this study uses Gaussian beam with awaist radius of 50 μm as the incident light source and the wavelength of 0.6328 μm.
The rainbow intensity distribution based on the Debye series expansion simulation only considers the contribution of specific rays. In actual measurement, the total light intensity distribution is produced by the interference of all light rays. These extra interference phenomena cause high-frequency oscillations superimposed on the rainbow pattern, that is, ripple structure. According to the generalized Lorenz-Mie theory (GLMT), the total light intensity distribution in the rainbow region (also known as the rainbow diagram) is calculated. The wavelength of the Gaussian beam is 0.6328 μm, and the beam waist radius ω0 = 100 μm. The size of the droplet is Between 50 μm and 200 μm. The rainbow diagram of the droplet with a radius of 50 μm is shown in Figure 3. The Butterworth digital filter is used for processing to obtain a smooth curve as shown in Figure 1. The cut-off frequency of the low-pass filter is determined based on the elimination of the number of redundant inflection points after roughly estimating the position of the second trough.
The condition of central incidence is studied, where the center of the Gaussian beam is located on the Déscartes ray [18]. According to the inversion algorithm given in Section 2, the calculated refractive index and its absolute error are shown in Figure 4. The calculated radius and its relative error are shown in Figure 5.
For the trough-trough (θmin1 and θmin2) algorithm, the absolute error of the refractive index is between −6.4 × 10−4 and 1.7 × 10−4 (see Figure 4(b), and the error of the droplet radius is only between −0.55% and 1.77% (see Figure 5(b)). For inflection-inflection (θinf1 and θinf2) algorithm, the inverted refractive index is less than −1.1 × 10−3 (see Figure 4(b)) and error of the droplet radius is between −0.75% and 5.67% (see Figure 5(b)). Obviously, both algorithms can be used for measuring droplet parameters, but the first algorithm is more suitable for calculating the droplet size than the second algorithm.
![]()
![]()
Figure 2. The influence of Gaussian beam incidence with different beam waist radius on the peak position of the secondary rainbow light intensity distribution of the droplet: (a) for the first peak; (b) for the first trough.
![]()
Figure 3. Normalized light intensity distribution for the secondary rainbow of a spherical droplet with r = 50 μm illuminated by Gaussian beam.
![]()
![]()
Figure 4. (a) The comparison between the refractive index of the droplet extracted from the secondary rainbow pattern and the true value; (b) The absolute error of the calculated refractive index.
![]()
![]()
Figure 5. (a) Comparison of the droplet radius extracted from the second-order rainbow pattern with the true value; (b) The relative error of the calculated droplet radius.
4. Summary
The generalized Lorenz-Mie theory is used to simulate the rainbow scattering pattern of the droplets in the secondary rainbow region. The low-pass filter passes the simulated rainbow pattern to obtain the characteristic angles. Two different inversion algorithms are studied to invert the droplet parameters. The results show that different inversion algorithms are also feasible to measure droplet information. In particular, for the first algorithm, the error of the droplet radius is even less than 1.77%. However, in this study, the selection of the cut-off frequency of the digital low-pass filter requires artificial evaluation of the angular position of the rainbow diagram. Therefore, we need to develop a better method to select the cutoff frequency or find a better filter, which we need to study further.