Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Droplet sizing in atomizing sprays using polarization ratio with structured laser illumination planar imaging

Open Access Open Access

Abstract

Previous research has shown that the polarization ratio technique allows the characterization of the surface mean diameter, D21, of droplets forming dilute sprays. However, its application to optically dense sprays has posed significant challenges due to the presence of multiple light scattering. Additionally, errors in measurement can arise from the angular dependence of the signal. In this Letter, we present a novel, to the best of our knowledge, method that addresses these challenges. Our approach combines the use of a telecentric objective with structured laser illumination, to both optimize light detection and suppress the unwanted intensity from multiple scattering. This approach enables the utilization of the polarization ratio technique for measuring the droplet size of challenging atomizing sprays. The method offers a promising solution for accurate and comprehensive spray characterization. It is applied, here, to a hollow-cone water spray running at 30, 50, and 70 bar injection pressure, reaching an optical depth up to three.

© 2023 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

To reach the requirements of carbon emission reduction, global energy consumption must be optimized. Optimizing the combustion of electrofuels (e-fuels, including butanol and other alcohols [1]) requires assessing the ability of an injector to produce spray systems that efficiently transit from liquid to gas. Thus, there is a need for detailed quantitative information on the droplet cloud. In particular, mean droplet sizes and size distributions are key parameters in optimizing these processes.

In order to determine these parameters, various optical measurement techniques exist. The most widely used for determining the droplet size distribution of a spray is the phase Doppler anemometry (PDA) technique, also called phase Doppler interferometry (PDI) [2,3]. This experimental measurement does not require calibration and allows the simultaneous determination of droplet size and velocity. However, the measurement volume of this method is limited to a few hundred microns which makes the technique time-consuming for the complete characterization of an entire spray. To reduce the time of experimental measurements while maintaining a high spatial resolution, techniques based on laser sheet imaging appear to be a solution. Thus, laser-induced fluorescence (LIF)/Mie [4] and Raman/Mie [5] methods have been developed to allow planar imaging of the Sauter mean diameter. Recently, planar rainbow refractometry has been developed for directly imaging the standard droplet diameter [6].

However, despite their advantages, these imaging methods are limited to the study of low concentrations of droplets in space. Indeed, the application to optically dense sprays is challenging due to the occurrence of a phenomenon known as multiple light scattering. By suppressing this unwanted intensity contribution, structured laser illumination planar imaging (SLIPI) [7] has allowed the implementation of the LIF/Mie droplet sizing in optically dense sprays [8]. Nevertheless, due to the necessity of adding fluorescent dyes to the liquid, LIF/Mie droplet sizing is limited to the study of non-evaporating sprays [9]. Furthermore, the addition of a dye in significant proportions can alter the physical properties of the fluid [10].

In order to measure the droplet size field without the addition of fluorescent dyes, while maintaining a good signal-to-noise ratio, the polarization ratio technique can be used. This technique has been introduced in the 1980s [11,12], where the behavior of the polarization ratio has been studied on isothermal sprays and burning conditions. It has also been applied to size droplets from various monodisperse streams [13,14]. The technique is based on collecting the perpendicular ($S$-$pol$) and parallel polarized ($P$-$pol$) light emitted via Lorenz–Mie scattering in laser sheet imaging. It allows obtaining an image whose intensity ratio is proportional to the statistical surface mean diameter ($D_{21}$). This is detailed as follows: the measurement is only applicable in the case of spherical droplets. According to Ref. [15], the scattered flux of light by a spherical particle $P(\theta,m,D)$ [W], measured within the solid angle $\Delta \Omega$ [sr$^{-1}]$ and from the scattering volume $\Delta V$ [m$^3]$, is linearly proportional to the incident light flux density $I_{i}$ [W m$^{-2}]$ through the scattering coefficient $\mu (\theta,m,D)$ [m$^{-1}$ sr$^{-1}]$:

$$\frac{P(\theta,m,D)}{\Delta V \; \Delta\Omega} = \mu(\theta,m,D) \; I_{i}.$$
This includes the scattering angle $\theta$ [deg], the refractive index $m = n + ik$, and the particle diameter $D$ [m]. Thus, in the single-scattering regime the measured quantities are directly proportional to the angular scattering coefficient,
$$\mu(\theta) = N \; \sigma(\theta,m,D),$$
where $N$ [#/m$^3]$ is the number density of droplets and $\sigma$ [m$^2]$ the scattering cross section. The major issue of the technique relies on the fact that strong oscillations appear in the detected scattered light intensity. This is caused by diffraction effects resulting in an ambiguous relationship between droplet diameter and polarization ratio [16,17]. It has been shown in Ref. [18] that the use of femtosecond lasers can significantly reduce the number and magnitude of polarization ratio oscillations; thus improving the accuracy of the technique. An alternative solution introduced in Ref. [16] is the application of the technique to a polydisperse spray, where the measurements are derived from a statistic of the droplet size passing through the measurement volume. Consequently, the strong oscillations of droplets are attenuated by the contribution of the large size distribution. For polydisperse sprays, the droplets size distribution is described as $f(D)$. Under this condition, Eq. (2) can be rewritten for each polarized scattering coefficient as
$$\mu_{S}(\theta) = N \int_{0}^{\infty} \sigma_{S}(\theta,m,D) f(D) {\rm d}D \; \propto N \int_{0}^{\infty} D^2 f(D) {\rm d}D,$$
$$\mu_{P}(\theta) = N \int_{0}^{\infty} \sigma_{P}(\theta,m,D) f(D) {\rm d}D \; \propto N \int_{0}^{\infty} D f(D) {\rm d}D.$$
According to the theoretical analysis of the method, the $S$-$pol$ signal is proportional to $D^2$, while the $P$-$pol$ signal is proportional to $D$. Figure S1 in Supplement 1 presents simulations results using the MiePlot software [19]. The dependency of the $S$-$pol$ and $P$-$pol$ related to the droplet diameter is given. Extensive work [16,17], has been performed where the polarization ratios are calculated for different refractive indices at $90^\circ$ scattering angle detection. This demonstrated that the sensitivity of droplet sizing is reduced when the refractive index of the liquid is increased. As a result, the technique is applicable up to a refractive index of $n=1.41$. In planar laser imaging measurements, a detection volume containing droplets of different size is imaged, and the scattering signal corresponds to the sum of the signal of many individual droplets ($k$). Considering the signal emitted by the droplets as a scattered flux density, the ratio of the two perpendicularly and parallelly polarized intensities ($I_S$ and $I_P$, respectively) can be written as follows:
$$I_{S/P}=\frac{I_S}{I_P}=\frac{C_S \sum_{k}^{} (D^2_k)}{C_P \sum_{k}^{} (D_k)} \sim D_{21},$$
with $I_S$ and $I_P$ being the $S$-$pol$ and $P$-$pol$ intensity. Here $C_S$ and $C_p$ are two constants experimentally determined by a calibration. Phase Doppler measurements or a well-controlled monodisperse droplets stream is generally used to calibrate the image ratio. However, although the polarization technique is applicable to spray systems, it is still limited to the study of dilute scattering media. In addition, the technique requires different corrections, due to the angular dependence of the signal collected on each pixel. In order to overcome these problems and be able to obtain two-dimensional (2D) measurement of the $D_{21}$ in optically dense sprays, polarization ratio laser sheet imaging has been combined with structured illumination and the use of a telecentric objective.

The optical setup used is illustrated in Fig. 1 and consists in generating a structured laser sheet of approximately 12 cm in height. For this purpose, a continuous-wave laser beam of 532 nm (5 W maximum power) is used and converted into a laser sheet by the use of cylindrical lenses. The modulation of the light intensity is implemented by means of a diffractive optical element (DOE), acting like a light-efficient Ronchi grating. The spray used for this study is a hollow cone steady water spray, formed under atmospheric conditions. The scattering signal is recorded using a large telecentric objective (0.066X, 1 in. C-Mount TitanTL Telecentric Lens) placed at $\theta = 90^\circ$ and coupled to two sCMOS cameras (Andor, Zyla 5.5) resulting in a resolution of 0.29 mm per pixel. As shown by Ref. [20], the use of a telecentric lens allows us to remove the angular dependence of the scattered Mie signal in the field of view. In order to form the same image simultaneously on both cameras while using a single telecentric lens, a TwinCam system (CairnOptics) has been used. This system consists of an imaging relay that formed the image on a beam splitter cube and two additional imaging relays to form the image on both cameras. Polarization filters (of optical density six) have been placed in front of each camera to collect the desired polarization on each sensor. Also, in order to smooth out the ripples from the Lorenz–Mie scattered light, the images have been recorded using a long exposure time of 50 ms. Thus, several droplets are present on one pixel, resulting in a spatial average and smoothing of the intensity ripples. A zero-order, air-spaced $\lambda /2$ wave plate has been used to adjust the polarization axis of the incident beam. In order to use the polarization ratio method, the polarizer is oriented at $45^\circ$. The effect of the orientation of the incident polarization on the ratio is shown in Fig. S2 of Supplement 1. These simulated calculations also demonstrate numerically the feasibility of the technique in our experimental conditions. optical depth ($OD$) measurements of the spray are realized simultaneously, by using a third camera that records the fluorescence signal from a cell containing Eosin Y in water placed behind the spray (Fig. 1). A conventional lens (Nikkor ${\rm f} = 105$ mm) is used to obtain the same zoom level on the images.

 figure: Fig. 1.

Fig. 1. Scheme of the experimental setup for 3p-SLIPI polarization ratio measurement. A spatially modulated light sheet is generated by means of a diffractive optical element. The incident light is adjusted to a polarization angle of $45^\circ$ using a half wave plate. The scattering signal is collected using a telecentric lens (placed at $\theta = 90^\circ$), and the image corresponding to the parallel and perpendicular signal are formed on two cameras using a beam splitter cube and two polarizing filters. A third camera recording the signal transmission in a cuvette filled with a fluorescent dye is added to allow $OD$ measurement.

Download Full Size | PDF

To suppress the undesirable effects of multiple scattering and to maintain the spatial resolution offered by the detection setup, the three-phase SLIPI method (3p-SLIPI) has been implemented. The final image is generated by recording three modulated images where the phase of the sinusoidal modulation of the intensity is shifted by $2 \pi /3$. This is experimentally done by vertically displacing the DOE. The final image of light intensity $I_{SLIPI}$ is obtained using the following equation:

$$I_{SLIPI}=\frac{\sqrt{2}}{3} \sqrt{[I_{0}-I_{2\pi/3}]^2+[I_{0}-I_{4\pi/3}]^2+[I_{2\pi/3}-I_{4\pi/3}]^2}.$$
One modulated sub-image of the $S$-$pol$ and $P$-$pol$ signals is shown in Figs. 2(a) and 2(b). The spatially modulated light sheet is entering the spray from the left and traverses it along its central axis. The resulting SLIPI images are shown for each polarization. Firstly, for both polarizations, a laser extinction effect can be observed in the spray edges. Secondly, the signal decay observed along the spray is greater for $S$-$pol$ than for $P$-$pol$. This can be explained by the fact that for this spray, small droplets are located in the center and that $S$-$pol$ intensity is proportional to the diameter squared while $P$-$pol$ is proportional to the diameter.

 figure: Fig. 2.

Fig. 2. SLIPI, $S$-$pol$, and $P$-$pol$ averaged images of the spray running at 50 bar are generated by recording three modulated sub-images. These modulated images are recorded with a phase shift of $120^\circ$, and only the first phase is depicted here in panel (a). (b) SLIPI images are then obtained using Eq. (6). For comparison, the images obtained with a conventional light sheet are added. The polarization ratio obtained using the conventional and SLIPI images are also illustrated in panel (c).

Download Full Size | PDF

In order to compare conventional and SLIPI results, the images are post-processed according to the same method: the images are thresholded above a value equal to five times the background value (intensity value outside of the spray region), and 200 images are averaged to generate the conventional and SLIPI images. Figure 2(c) presents the polarization ratio obtained with a conventional light sheet excitation and using 3p-SLIPI [Eq. (6)]. Firstly, a non-zero polarization ratio is observed in the upper part of the conventional image, in contrast to the SLIPI image. This part of the spray is not directly illuminated by the laser sheet but by multiple light scattering. As observed, this unwanted signal is well removed in the case of the SLIPI image. In addition, it can be observed that the spray symmetry is recovered using the 3p-SLIPI technique. At this point, according to Eq. (5) the ratio is proportional to the $D_{21}$. However, to obtain an absolute measurement, a calibration of the polarization ratio is necessary. This is performed using PDI (Artium instrument PDI-TK2) measurements at different vertical positions from the center of the spray to its edges ($Z_1=7$, $Z_2=9$, $Z_3=11$, and $Z4=13$ cm from the injector). Data are acquired every millimeter and are validated when 100 000 droplets are acquired for each measurement point. In order to obtain a geometrical correspondence between the PDI measurements and the polarization ratio, an image of the measurement volume is recorded on each camera. A value of the polarization ratio is assigned for each PDI measurement and is shown in Fig. 3. For the case of conventional illumination, it is observed that for an identical polarization ratio, different $D_{21}$ values are measured as a function of the height in the spray. Thus, the polarization ratio cannot be calibrated in the case of optically dense sprays when using conventional illumination. On the contrary, when using SLIPI, a single solution of the $D_{21}$ is found for the polarization ratio. A calibration curve is then determined by applying a second-order power law allowing to deduce the $D_{21}$ from the polarization ratio measurement.

 figure: Fig. 3.

Fig. 3. Calibration of the polarization ratio as a function of the $D_{21}$ when using (a) conventional illumination and (b) 3p-SLIPI at different distances from the nozzle indicated in Fig. 2. The fact that the relationship between polarization ratio and $D_{21}$ is not a bijective function when conventional illumination is used makes the calibration of the measurement impossible. However, when SLIPI is used, a unique solution for each $D_{21}$ is found showing the ability to calibrate the technique.

Download Full Size | PDF

After having calibrated the SLIPI polarization ratio, a 2D map of the absolute $D_{21}$ is obtained and presented in Fig. 4(a) (injection pressure of 50 bar). As observed, small droplets are located in the spray center while larger droplets are present on the spray edges. Measurements have been then made at two additional injection pressures (30 bar and 70 bar). $OD$ measurements have been carried out simultaneously and show the increase in spray density while increasing injection pressure [(Fig. 4(b)]. The $D_{21}$ measurements allow to observe the modification of the granulometry of the spray as a function of the injection pressure. Figure 4(c) presents the comparisons between 3p-SLIPI polarization ratio measurements and PDI measurement at a distance of 7 cm from the injector. As observed, whatever the injection pressure, the imaging data are in good agreement with those obtained by PDI. In conclusion, the SLIPI polarization ratio has been applied for the first time to the characterization of the $D_{21}$ on an entire spray having $OD$ ranging from 0.4 to 3. However, it is still possible to apply the technique to sprays with larger $OD$. However, the full characterization of the spray will require reconstruction of the image due to laser extinction along the spray. Mie intensity ripples are reduced when measurements are carried out in a polydisperse environment such as a spray. Additionally, the angular dependence of the Mie signal in the resulting image is suppressed by the use of a telecentric objective. In addition, the calibration of the polarization ratio is achievable when structured illumination is used but remains impossible with conventional light sheet imaging. In order to strengthen the calibration, comparisons with other measurement techniques such as shadowgraphy, digital inline holography [21], or planar rainbow refractometry [6] may be carried out in the future. The technique is applicable to evaporating sprays due to the fact that no dyes are added to the liquid. Instantaneous measurements may also be achievable if the polarization ratio technique is coupled with an ultra-short laser excitation (femtosecond laser) [18] or with the use of a non-fluorescent dye absorbing the laser light in the liquid.

 figure: Fig. 4.

Fig. 4. (a) 2D 3p-SLIPI maps of the $D_{21}$, obtained for three injection pressures 30, 50, and 70 bar. (b) $OD$ measurements as a function of the distance from the nozzle. (c) Sizing validated with PDI measurements made at a distance of 7 cm from the nozzle and compared with the 3p-SLIPI for the different investigated injection pressure. An $R^2$ value is added to the plot to assess the correspondence between PDI and imaging measurements.

Download Full Size | PDF

Funding

Vetenskapsrådet (Etenskapsrådet 2021–04542).

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Supplemental document

See Supplement 1 for supporting content.

REFERENCES

1. A. Peter, B. Siewert, S. Riess, L. Strauss, C. Pastoetter, and M. Wensing, Atomization Sprays 30, 843 (2020). [CrossRef]  

2. F. Durst and M. Zaré, in The Accuracy of Flow Measurements by Laser Doppler Methods; Proceedings of the LDA-Symposium, Copenhagen, Denmark, August 25–28, 1976, pp. 403–429 .

3. W. D. Bachalo, Appl. Opt. 19, 363 (1980). [CrossRef]  

4. C. Yeh, in The Third International Congress on Optical Particle Sizing (1993).

5. A. Malarski, B. Schürer, I. Schmitz, L. Zigan, A. Flügel, and A. Leipertz, Appl. Opt. 48, 1853 (2009). [CrossRef]  

6. C. Li, Q. Lv, N. Li, Y. Wu, X. Wu, C. Weng, and C. Tropea, Opt. Lett. 46, 5898 (2021). [CrossRef]  

7. E. Berrocal, E. Kristensson, M. Richter, M. Linne, and M. Aldén, Opt. Express 16, 17870 (2008). [CrossRef]  

8. Y. N. Mishra, E. Kristensson, and E. Berrocal, Opt. Express 22, 4480 (2014). [CrossRef]  

9. B. Frackowiak and C. Tropea, Appl. Opt. 49, 2363 (2010). [CrossRef]  

10. M. Stiti, A. Labergue, F. Hervy, G. Castanet, S. Becker, and F. Lemoine, Exp. Fluids 62, 169 (2021). [CrossRef]  

11. F. Beretta, A. Cavaliere, A. Ciajolo, A. D’Alessio, C. Langella, A. Di Lorenzo, and C. Noviello, Symp. (Int.) Combust., [Proc.] 18, 1091 (1981). [CrossRef]  

12. F. Beretta, A. Cavaliere, and A. D’Alessio, Symp. (Int.) Combust., [Proc.] 20, 1249 (1985). [CrossRef]  

13. P. Massoli, F. Beretta, and A. D’Alessio, Appl. Opt. 28, 1200 (1989). [CrossRef]  

14. D. L. Hofeldt and R. K. Hanson, Appl. Opt. 30, 4936 (1991). [CrossRef]  

15. H. Van de Hulst, Light Scattering by Small Particles (Dover, 1981).

16. D. L. Hofeldt, Appl. Opt. 32, 7551 (1993). [CrossRef]  

17. D. L. Hofeldt, Appl. Opt. 32, 7559 (1993). [CrossRef]  

18. S. Bareiss, B. Bork, S. Bakić, C. Tropea, R. Irsig, J. Tiggesbäumker, and A. Dreizler, Meas. Sci. Technol. 24, 025203 (2013). [CrossRef]  

19. P. Laven, Appl. Opt. 42, 436 (2003). [CrossRef]  

20. S. Garcia, P. Doublet, C. Lempereur, G. Illac, M. Stiti, E. Berrocal, and M. Orain, Exp. Fluids (2023).

21. J. Fugal and R. Shaw, Atmos. Meas. Tech. 2, 259 (2009). [CrossRef]  

Supplementary Material (1)

NameDescription
Supplement 1       Supplementary_Material

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (4)

Fig. 1.
Fig. 1. Scheme of the experimental setup for 3p-SLIPI polarization ratio measurement. A spatially modulated light sheet is generated by means of a diffractive optical element. The incident light is adjusted to a polarization angle of $45^\circ$ using a half wave plate. The scattering signal is collected using a telecentric lens (placed at $\theta = 90^\circ$), and the image corresponding to the parallel and perpendicular signal are formed on two cameras using a beam splitter cube and two polarizing filters. A third camera recording the signal transmission in a cuvette filled with a fluorescent dye is added to allow $OD$ measurement.
Fig. 2.
Fig. 2. SLIPI, $S$-$pol$, and $P$-$pol$ averaged images of the spray running at 50 bar are generated by recording three modulated sub-images. These modulated images are recorded with a phase shift of $120^\circ$, and only the first phase is depicted here in panel (a). (b) SLIPI images are then obtained using Eq. (6). For comparison, the images obtained with a conventional light sheet are added. The polarization ratio obtained using the conventional and SLIPI images are also illustrated in panel (c).
Fig. 3.
Fig. 3. Calibration of the polarization ratio as a function of the $D_{21}$ when using (a) conventional illumination and (b) 3p-SLIPI at different distances from the nozzle indicated in Fig. 2. The fact that the relationship between polarization ratio and $D_{21}$ is not a bijective function when conventional illumination is used makes the calibration of the measurement impossible. However, when SLIPI is used, a unique solution for each $D_{21}$ is found showing the ability to calibrate the technique.
Fig. 4.
Fig. 4. (a) 2D 3p-SLIPI maps of the $D_{21}$, obtained for three injection pressures 30, 50, and 70 bar. (b) $OD$ measurements as a function of the distance from the nozzle. (c) Sizing validated with PDI measurements made at a distance of 7 cm from the nozzle and compared with the 3p-SLIPI for the different investigated injection pressure. An $R^2$ value is added to the plot to assess the correspondence between PDI and imaging measurements.

Equations (6)

Equations on this page are rendered with MathJax. Learn more.

P ( θ , m , D ) Δ V Δ Ω = μ ( θ , m , D ) I i .
μ ( θ ) = N σ ( θ , m , D ) ,
μ S ( θ ) = N 0 σ S ( θ , m , D ) f ( D ) d D N 0 D 2 f ( D ) d D ,
μ P ( θ ) = N 0 σ P ( θ , m , D ) f ( D ) d D N 0 D f ( D ) d D .
I S / P = I S I P = C S k ( D k 2 ) C P k ( D k ) D 21 ,
I S L I P I = 2 3 [ I 0 I 2 π / 3 ] 2 + [ I 0 I 4 π / 3 ] 2 + [ I 2 π / 3 I 4 π / 3 ] 2 .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.