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

Instantaneous thermometry imaging using two-photon laser-induced fluorescence

Open Access Open Access

Abstract

Measuring temperature in complex two-phase flows is crucial for understanding the dynamics of heat and mass transfer. In this Letter, we introduce a novel, to the best of our knowledge, optical approach based on the combination of two-photon laser-induced fluorescence (2p-LIF) imaging and two-color laser-induced fluorescence (2CLIF) for instantaneous temperature mapping of complex liquid media. Using Kiton Red (KR) and Rhodamine 560 (R560), a temperature sensitivity of $1.54 {\% }/^{\circ }$C has been achieved over a range of 17–60°C. The monitoring of two-dimensional transient temperature dynamics in the heating and degassing of water shows the efficiency of the 2p–2CLIF. This new approach contributes to the toolkit of optical temperature measurement techniques, providing a robust solution for studying transient scattering media and high-speed two-phase flows.

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

Optical temperature measurements in liquid flows have been the focus of extensive research for several decades, holding particular relevance in situations where conventional devices such as thermometer sensors, including thermocouples, may not be applicable [1]. A wide range of physical phenomena present various optical opportunities for temperature determination. Examples include the analysis of vibrational energy modes with Raman scattering spectroscopy, measuring thermal radiation with IR thermometry, and the use of thermochromic liquid crystals [2]. Even though the temporal response of these techniques may differ significantly, they are in principle all adequate for the investigation of single-phase media, such as turbulent heated jets [3]. The thermometry of complex two-phase media, such as bubble flows, particle suspensions, or spray systems, however, still remains challenging [4]. Laser-induced fluorescence (LIF) techniques stand out as the favored choice offering advantages in terms of localization, species sensitivity, and control over the illuminated volume where liquid flow and transfer phenomena can entail complex geometries. LIF involves the introduction of a fluorescent dye, whose fluorescence emission displays pronounced sensitivity to temperature, in terms of variations in emission wavelengths, intensity, and lifetime.

Early measurements of liquid temperature with LIF, in the 1990s, involved point measurements using a focused laser beam [5]. Over time, advancements in camera and laser technology led to the introduction of planar laser-induced fluorescence (PLIF) for two-dimensional temperature measurements [6]. This method involves the illumination of liquid flows with a thin laser sheet confining the fluorescence emission to an area. More recently, stereoscopic systems have been developed to reconstruct three-dimensional temperature fields [7].

Within a small elementary volume denoted as $dV$, the single-photon fluorescence intensity, $dF$, can be mathematically expressed as follows:

$$dF(\lambda)=\frac{d \Omega}{4 \pi} \; I(x,y,z) \; \phi(\lambda) \; \varepsilon(\lambda_0) \;C \;dV ,$$
where $I(x,y,z)$ is the laser excitation flux in $\rm {W/m}^2$, $\epsilon$ is the molar coefficient of the single-photon absorption at the laser wavelength $\lambda _0$, $C$ is the dye concentration, and $\phi (\lambda )$ represents the ratio of emitted photons at the wavelength $\lambda$ compared to the total number of absorbed photons. The parameter $d \Omega$ is the solid angle under which the signal is received by the optical system. In this expression, the parameters $\phi$ and $\varepsilon$ are temperature sensitive. In particular, $\phi$ is known to decrease significantly with temperature for dyes such as Kiton Red (KR) [8].

To determine the temperature based on Eq. (1), the flow must remain single phase, while, simultaneously, the laser flux and measurement volume must be kept constant with respect to a reference measurement. As a result of these stringent restrictions, disturbances in the form of striations due to refractive index variations pose significant challenges. For two-phase flows, discontinuities in refractive indices and the presence of multiple interfaces (e.g., liquid/gas) introduce notable degrees of light scattering, making the direct use of fluorescence intensity for thermometry problematic.

To address this challenge, one can resort to ratiometric techniques [9]. These techniques, also referred to as two-color LIF thermometry (2CLIF), rely on the simultaneous measurement of the fluorescence intensity in two spectral bands, each exhibiting different sensitivities to temperature variations. This can be achieved by using either a single dye whose emission spectrum changes with temperature [10] or two dyes having separated absorption spectra and exhibiting different temperature sensitivities [11]. This approach has been applied in diverse temperature reconstruction applications [1214]. The efficiency of ratiometric approaches relies on achieving equal scattering and absorption impacts across the two detection bands. Consequently, the utilization of these methodologies is confined to some applications of low complexity, as previously described.

For high-complexity, two-phase media, refractive index variations and multiple scattering will extend the fluorescing region beyond the laser sheet. Furthermore, the lengthening of optical paths in these types of media will favor the reabsorption and re-emission of fluorescence, encroaching on the stringent requirements of ratiometric thermometry.

An alternative solution is the use of two-photon laser-induced fluorescence (2p-LIF), a nonlinear process involving the excitation of the aforementioned temperature-sensitive fluorescent dyes via the simultaneous absorption of two infrared photons. Due to the necessity of simultaneous absorption, only areas of high photon flux, i.e., areas within the laser path, will absorb and hence emit 2p-LIF signal, granting this technique robustness toward, e.g., scattering environments, a characteristic recently employed for pointwise temperature measurements within the optically dense regions of sprays [15]. Recent advances in femtosecond laser technology have resulted in high photon flux pulses that generate enough 2p-LIF signal for detection in planes several centimeters wide, paving the way for new imaging applications especially in optically dense samples [16].

In this study, we utilize the benefits of 2p-LIF imaging to 2CLIF thermometry in order to obtain instantaneous temperature maps in complex two-phase media. With the use of Kiton Red and Rhodamine 560 (R560) as the temperature probes, we show a temperature sensitivity of $1.54 {\% }/^{\circ }$C over a temperature range of 17–60$^{\circ }$C. The use of the 2p–2CLIF technique in heating and degassing water illustrates its effectiveness in thermometry for complex and dynamic high-speed liquid flows. A macroscopic laser sheet measuring $\sim 1\times 1$ cm$^2$, with a pulse duration of $45$ fs, $1$ mJ energy and centered at $840$ nm, is incident on a cuvette filled with dissolved fluorescent dyes. The fluorescence signal is detected, via a telecentric objective (1.0$\times$ GoldTL) positioned at a 90$^{\circ }$ angle from the laser sheet, on two sCMOS cameras (Andor, Zyla 5.5) resulting in a spatial resolution of $6.5$ $\mathrm {\mu }$m per pixel. To simultaneously capture the same image on both cameras while using a single telecentric lens, a Twincam system (CairnOptics) has been implemented (Fig. 1(a)). This system comprises an imaging relay that via a beam splitter cube, generates duplicate images on both cameras. For efficient collection of the fluorescence signal in two distinct spectral bands, bandpass fluorescence filters ([575–625] nm and [488–512] nm, whose spectral properties are depicted in Fig. 1(b)) are placed in front of each camera.

 figure: Fig. 1.

Fig. 1. (a) Illustration of the experimental setup used for generating a macroscopic laser sheet measuring approximately $1\times 1$ cm$^2$. An ultra-short, $1$ mJ laser pulse ($45$ fs) from a ytterbium-doped fiber seeded laser centered at $840$ nm is focused into the cuvette with a cylindrical lens (CL). Illustration of the calibration and application cases has been added. Detection is achieved in two spectral channels via a CairnOptics Twincam setup, two spectral filters, and two Andor Zyla sCMOS cameras. The experimental setup used for the temperature calibration and for the study of the heating of non-degassed water is also depicted. (b) Spectrum of the laser excitation source as well as the fluorescence spectra of Rhodamine 560 and Kiton Red. The shaded areas represent the spectral positions of the bandpass filters. (c) Graphical representation of the image processing routine used for 2p–2CLIF ratiometric temperature calibration.

Download Full Size | PDF

For 2p-LIF applications, Eq. (1) must be updated:

$$dF(\lambda)=\frac{d \Omega}{4 \pi} \; I^2 \; \phi(\lambda) \; \delta_0 \;C \;dV,$$
where $\delta _0$ is related to the molar coefficient of two-photon absorption in the laser emission band. Similar to the case of one-photon absorption, the temperature dependence of the fluorescence process mainly arises from the fluorescence yield and to some extent from the absorption coefficient (i.e., the parameters $\phi$ and $\delta$). Assuming the volume of two-photon absorption remains approximately constant with temperature (i.e., no temperature dependence of $\delta$), the fluorescence signal $F$, collected in a given spectral band $i$ from a volume $dV$, can be rewritten as
$$dF_i=\eta_i \; K_i \;c \;I^2\; e^{s_i T} dV,$$
where $K_i$ is a parameter depending on the solid angle of detection (constant for a given optical system) and $\eta _i$ is the transmission efficiency of fluorescence to the detector (affected by signal reabsorption and scattering). The parameter $s_i$ introduced in Eq. (3) represents the temperature sensitivity of the dye in the spectral band, given in ${\% }/^{\circ }$C. Similar to the single-photon case, in order to omit unreliable integration steps, the fluorescence signal is detected in two spectral bands, and the intensity ratio $R_{12}= F_1 / F_2$ is calculated. These spectral bands are affected by optical effects in a similar manner and are almost not reabsorbed, resulting in a secondary fluorescence signal of less than 1% in comparison to the direct fluorescence (due to an optically thin medium despite scattering because of low dye concentrations). Hence, in our case, it can be assumed that $\eta _1$ and $\eta _2$ vary in tandem, leaving
$$R_{12}=\frac{F_1}{F_2}=K \frac{C_1 e^{s_1 T}}{C_2 e^{s_2 T}},$$
where $K$ is a constant. Finally, assuming that the dyes locally remain in the same concentration proportions, a reference $R_{12}^{\text {ref}}=I_1(T_0) / I_2(T_0)$ at a known temperature ($T_0$) can be used to eliminate $K$:
$$\frac{R_{12}}{R_{12}^{\text{ref}}}=e^{(s_1 - s_2) (T-T_0)}.$$
As shown in Eq. (5), knowledge of the temperature sensitivity of the dyes, $s_1$, $s_2$, allows for direct thermometry measurements based on a fluorescence intensity ratio with respect to a reference. These parameters are obtained by a temperature calibration where the fluorescence ratio $R_{12}$ is measured for different liquid temperatures.

In this study, the mixing of Kiton Red (KR at a concentration of $10^{-6}$ mol/L) and Rhodamine 560 (R560 at a concentration of $10^{-5}$ mol/L) was used as the temperature-sensitive mixture. In relation to the excitation laser, their measured spectra are illustrated in Fig. 1(b). The calibration was conducted in a cell placed on a temperature-controlled heating plate with a path length of less than $1$ mm between the laser sheet and the cell wall. In order to ensure temperature homogeneity, a magnetic stirrer is inserted along with a thermocouple (Fig. 1(a)). Figure 2(a), top panel, represents the variation of the fluorescence signal on each spectral band as a function of temperature. As observed, the fluorescence signal emitted on spectral band 1 ([575–625] nm, related to KR) decreases with temperature, while the signal in spectral band 2 ([488–512] nm, related to Rh560) increases slightly. Figure 2(b) presents the evolution of the mean intensity on each image as well as the logarithm of the normalized ratio (see Fig. 1(c) for a graphical representation of the image processing). The fitting of the calibration curve yields a temperature sensitivity of 1.54 ${\% }/^{\circ }$C from 17 to 60$^{\circ }$C. The divergence of this value from literature [8] is attributed to the contamination by fluorescence from Rhodamine 560 in the fluorescence band designated for Kiton Red (Fig. 1(b)), reducing the temperature sensitivity on this specific spectral band. Furthermore, this reduction in temperature sensitivity is also due to the laser system we employed for this investigation. Indeed, fluorescence is occurring in the partially saturated excitation regime which could explain this sensitivity decrease for Kiton Red [3]. The temperature calibration obtained in Fig. 2 is applicable within the experimental condition described here. Any change in the fluorescent dyes, liquid solvent, and spectral bands of detection will necessitate a new calibration, following the same methodology. As for the effect of the optical path, the necessity of a new calibration will depend on whether the amount of light reabsorbed in the detection bands can be neglected.

 figure: Fig. 2.

Fig. 2. (a) Image ratios at various temperatures, calculated via the routine illustrated in Fig. 1(c). (b) Evolution of fluorescent intensity in both spectral bands where the logarithmic ratios are plotted on the right axis as a function of water temperature. This fitted line results in a sensitivity of 1.54${\% }/^{\circ }$C.

Download Full Size | PDF

Instantaneous two-photon fluorescence thermometry was carried out in a transparent cuvette filled with non-degassed water seeded with the two fluorescent dyes. Employing non-degassed water in our experimental setup facilitates the generation of a substantial quantity of air bubbles dispersed throughout the liquid, mimicking a complex medium. Additionally, these bubbles yield temperature gradients of greater magnitude compared to vapor bubbles arising from the conventional water boiling processes, a regime where temperature gradients are confined to a few degrees around the boiling point. Subsequently, a resistive heater is introduced into the cuvette (Fig. 1(a)). At the onset of the experiment, the water is maintained at an initial temperature of approximately $20^{\circ }$C. Around 10 s after the beginning of image acquisition, the heating element is triggered, initiating a temperature increase. The resistive heater functions in a cyclic manner, emitting heat at a constant frequency over time. The periods during which the heater is actively generating heat are denoted as “ON” in the illustrated images presented in Fig. 3. Throughout the experiment, high-resolution fluorescence images are systematically captured at a frequency of $5$ Hz. Figure 3 shows a selection of the images from the acquired series. Initially, the water exhibits a uniform temperature distribution around $20^{\circ }$C. At the edges of the laser sheet, the fluorescence signal drops rapidly due to the second-order relation between the laser and fluorescence intensity (and consequently signal-to-noise ratio). Upon heater activation, a discernible transformation in the liquid’s thermal field becomes evident with the formation of convective cells and the emergence of hot plumes, mainly visible from 40 s and onward. Notably, this heating process is concomitant with the production of bubbles within the liquid, which become entrained in the fluid flow.

 figure: Fig. 3.

Fig. 3. Selection of images acquired of the temperature field of water heated by a resistive heater. The gray box at $t=76$ s along with the corresponding red curve to the right indicates an equilibrium temperature gradient after heating has been completely turned off. This is in stark contrast to the flat temperature gradient of the water before heating has commenced (green curve).

Download Full Size | PDF

When the heater is completely deactivated at $t = 70$ s, a distinct vertical temperature stratification becomes evident, eventually stabilizing at $t = 76$ s. The temperature profile in the region highlighted by the gray zone and corresponding temperature plot to the right (red curve of Fig. 3) reveals a significant temperature decrease, transitioning from 50$^{\circ }$C at the free liquid surface to 20$^{\circ }$C at the base of the cuvette. This observed temperature gradient was verified through thermocouple measurements conducted immediately after the experimental trials and is in stark contrast to the corresponding constant temperature profile at $t=0$ s (shown in green). The temporal evolution of temperature in three regions of interest are represented in Supplement 1.

Due to the instantaneous nature and robustness against scattering of 2p–2CLIF, temperature measurements in the vicinity of air bubbles, i.e., within complex two-phase media, become possible. Three distinct time points within the acquisition, specifically $t = 51.0, 51.2, \rm {and}\,51.4$ s are analyzed with respect to bubble formation within three designated regions, denoted as zones (i), (ii), and (iii) (Fig. 4). (i) is placed within the cool environment, (ii) is positioned immediately after the hot wavefront, and (iii) encompasses the bubble inception site. The temperature profiles derived from these areas are displayed in the histograms underneath, along with the zoomed-in images of the different areas. This investigation shows that bubbles emerging directly from the heating element are situated within an environment characterized by significantly elevated temperatures compared to regions beyond the immediate vicinity of the bubbles. Indeed, an average temperature of $54^{\circ }$C prevails within the bubble-rich zone, whereas a temperature of $45^{\circ }$C is measured adjacent to this zone.

 figure: Fig. 4.

Fig. 4. Detailed views at specific instances ($51.0$, $51.2$, and $51.4$ s) reveal the presence of bubble dynamics with respect to the heater being on or off. Zoomed-in sections (i, ii, and iii at $t=51.2$ s) and corresponding temperature histograms are provided, revealing the temperature gradients of this transient processes.

Download Full Size | PDF

In conclusion, our investigation represents a significant advancement in the domain of optical temperature measurements within complex two-phase flows. The integration of two-photon laser-induced fluorescence (2p-LIF) imaging with two-color laser-induced fluorescence (2CLIF) thermometry presents a novel and powerful methodology for achieving instantaneous temperature mapping. We have shown that this approach proves effective in capturing the temperature dynamics in two-phase flows, addressing challenges posed by refractive index variations, multiple interfaces, and light-scattering phenomena. In this work, we have used the temperature-sensitive dyes Kiton Red and Rhodamine 560, achieving a temperature sensitivity of $1.54 {\% }/^{\circ }$C over temperatures ranging from 17 to 60$^{\circ }$C. By applying this mixture and technique to the heating and degassing of water, with the intentional introduction of bubbles, we reveal detailed insights into the thermal behavior of the liquid, including the formation of convective cells, emergence of hot plumes, and temperature stratification surrounding air bubbles. In comparison to single-photon approaches, fluorescence excitation with two photons has the main advantage of being applicable to the scattering medium such as dense spray (Fig. 2 in Supplement 1) as it eliminates absorption beyond the area subtended by the laser sheet. Moreover, the proposed approach, i.e., employing a powerful femtosecond laser, offers significant temporal resolution gains for the study of transient phenomena. We believe that 2p–2CLIF thermometry contributes to the expanding toolkit of optical temperature measurement techniques, offering a robust solution for the study of complex, transient, high-speed two-phase media in various liquid media.

Funding

European Regional Development Fund (ENERBATIN CPER); European Research Council (ERC 803634); Vetenskapsrådet (019-05183, 2021–04542); E-COOL/EIC Pathfinder (101130315).

Disclosures

The authors declare no conflicts of interest.

Data availability

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

Supplemental document

See Supplement 1 for supporting content.

REFERENCES

1. M. Yamazaki, T. Hayashi, H. F. Farahani, et al., Int. J. Thermofluids 19, 100384 (2023). [CrossRef]  

2. F. Lemoine and G. Castanet, Exp. Fluids 54, 1572 (2013). [CrossRef]  

3. W. Chaze, O. Caballina, G. Castanet, et al., Exp. Fluids 57, 58 (2016). [CrossRef]  

4. J. Zhou, W. Yang, Y. Yin, et al., Exp. Fluids 62, 125 (2021). [CrossRef]  

5. T. Nakajima, M. Utsunomiya, and Y. Ikeda, in Applications of Laser Techniques to Fluid Mechanics: 5th International Symposium, pp. 34–53 (1991).

6. M. Coolen, R. Kieft, C. Rindt, et al., Exp. Fluids 27, 420 (1999). [CrossRef]  

7. S. Kashanj and D. S. Nobes, Exp. Fluids 64, 46 (2023). [CrossRef]  

8. M. Stiti, A. Labergue, F. Lemoine, et al., Exp. Fluids 60, 69 (2019). [CrossRef]  

9. J. Sakakibara and R. J. Adrian, Exp. Fluids 26, 7 (1999). [CrossRef]  

10. G. Castanet and F. Lemoine, Proc. Combust. Inst. 31, 2141 (2007). [CrossRef]  

11. C. E. Estrada-Pérez, Y. A. Hassan, and S. Tan, Rev. Sci. Instrum. 82, 074901 (2011). [CrossRef]  

12. G. Castanet, P. Lavieille, M. Lebouché, et al., Exp. Fluids 35, 563 (2003). [CrossRef]  

13. W. Chaze, O. Caballina, G. Castanet, et al., Exp. Fluids 58, 96 (2017). [CrossRef]  

14. M. Piskunov and P. Strizhak, Exp. Therm. Fluid Sci. 91, 103 (2018). [CrossRef]  

15. M. Wang, M. Stiti, H. Chaynes, et al., Exp. Fluids 63, 172 (2022). [CrossRef]  

16. E. Berrocal, C. Conrad, J. Püls, et al., OSA Continuum 2, 983 (2019). [CrossRef]  

Supplementary Material (1)

NameDescription
Supplement 1       Supplementary material

Data availability

Data underlying the results presented in this paper are not publicly available 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. (a) Illustration of the experimental setup used for generating a macroscopic laser sheet measuring approximately $1\times 1$ cm$^2$. An ultra-short, $1$ mJ laser pulse ($45$ fs) from a ytterbium-doped fiber seeded laser centered at $840$ nm is focused into the cuvette with a cylindrical lens (CL). Illustration of the calibration and application cases has been added. Detection is achieved in two spectral channels via a CairnOptics Twincam setup, two spectral filters, and two Andor Zyla sCMOS cameras. The experimental setup used for the temperature calibration and for the study of the heating of non-degassed water is also depicted. (b) Spectrum of the laser excitation source as well as the fluorescence spectra of Rhodamine 560 and Kiton Red. The shaded areas represent the spectral positions of the bandpass filters. (c) Graphical representation of the image processing routine used for 2p–2CLIF ratiometric temperature calibration.
Fig. 2.
Fig. 2. (a) Image ratios at various temperatures, calculated via the routine illustrated in Fig. 1(c). (b) Evolution of fluorescent intensity in both spectral bands where the logarithmic ratios are plotted on the right axis as a function of water temperature. This fitted line results in a sensitivity of 1.54${\% }/^{\circ }$C.
Fig. 3.
Fig. 3. Selection of images acquired of the temperature field of water heated by a resistive heater. The gray box at $t=76$ s along with the corresponding red curve to the right indicates an equilibrium temperature gradient after heating has been completely turned off. This is in stark contrast to the flat temperature gradient of the water before heating has commenced (green curve).
Fig. 4.
Fig. 4. Detailed views at specific instances ($51.0$, $51.2$, and $51.4$ s) reveal the presence of bubble dynamics with respect to the heater being on or off. Zoomed-in sections (i, ii, and iii at $t=51.2$ s) and corresponding temperature histograms are provided, revealing the temperature gradients of this transient processes.

Equations (5)

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

d F ( λ ) = d Ω 4 π I ( x , y , z ) ϕ ( λ ) ε ( λ 0 ) C d V ,
d F ( λ ) = d Ω 4 π I 2 ϕ ( λ ) δ 0 C d V ,
d F i = η i K i c I 2 e s i T d V ,
R 12 = F 1 F 2 = K C 1 e s 1 T C 2 e s 2 T ,
R 12 R 12 ref = e ( s 1 s 2 ) ( T T 0 ) .
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.