**4-**system with a Ronchi grating placed at the frequency plane are used. This configuration permits the phase-shifted interferograms to be grabbed simultaneously by one CCD. The temperature field measurement is based on measuring the refraction index difference by solving the inverse Abel transform, which requires information obtained by the fringe order localization. The phase map is retrieved by a three-step algorithm. Experimental results of a dynamic thin flame are presented.

*f*## 1.

## Introduction

In phase-shifting techniques, a stepping motor or piezoelectric transducers are normally used to move a reference surface to get the interference fringe images with a relative phase shift. However, this takes a period of time, making it susceptible to environmental effects such as ambient vibration or air turbulence, which results in measurement errors. To facilitate the applicability of interferometric techniques for measuring dynamic events, a system that is capable of acquiring several phase-shifted images instantaneously or simultaneously is required. As a result of this capability, a reduction of environmental disturbances is encountered, increasing the accuracy and stabilization of the system.

In this study, a measurement system that is capable of grabbing three phase-shifted interferometric images instantaneously (or simultaneously) was developed for dynamic temperature measurement over a period of time. The system presented is based on polarization phase-shifting techniques presented previously,^{1}2.^{–}^{3} with the addition of calibration procedures and digital fringe processing techniques in order to accurately measure temperature fields varying in time.

Single-shot polarization phase-shifting techniques are currently a focus field of study, and the most important and an industry standard nowadays are pixelated phase mask interferometers. One of the principal properties of these systems is that the modulating phase-mask remains fixed, placed before the CCD light sensor.^{4} Nowadays, some authors present novel algorithms based on this property; for example, in the demodulation form, algorithms focused on harmonics rejections, and recently encountered more accurate results by proposing an extension of the phase-shifting unit cell.^{5}^{,}^{6}

Other single-shot systems are based on the use of a three-layer sensor to detect three wavelengths; these systems have been applied accurately by fringe projections profilometry^{7} and digital holographic system.^{8} Interferometric polarization phase-shifting techniques also have been employed before in different techniques; some of them are based on using Wollaston prisms and polarization gratings^{9} by using a prism and a glass plate with different thicknesses,^{10} and lately a system capable of acting as a fringe projection system, or a shearing interferometer by employing a Savart plate and liquid crystals.^{11}

Alternative systems, also based on phase-shifting polarization techniques, have been presented before.^{2}^{,}^{12} These systems are based on obtaining interference replicas by using phase/amplitude gratings,^{12} and also they only retrieve the phase data map. The interferogram amplitude modulation change obtained by the replication system was solved earlier by using two crossed gratings and selecting the orders with equal amplitude.^{2}^{,}^{3}^{,}^{12} Even this amplitude modulation change was minimum when using phase gratings,^{3}^{,}^{12} it was still encountered. In this work, the selection of an amplitude grating is to emphasize the amplitude modulation change and how this can be solved by a fringe normalization process.

In order to implement a system capable of obtaining temperature field measurements of a flame varying in time, we implemented calibration and normalization procedures.

In this work, the temperature field measurement is based on measuring the refraction index difference by solving the inverse Abel transform. The inversion of the Abel transform is solved with the information obtained by the fringe order localization of the deformed fringe pattern generated by the flame under study.

This paper is organized in five sections, including the present introduction. Section 2 discusses the theoretical background of the single-shot procedure by first explaining the phase data map retrieval and then the relationship with the temperature field measurement. In Sec. 3, the experimental setup is described, and the dynamic temperature field measurements are presented. Finally, in Sec. 4, conclusions are presented.

## 2.

## Temperature Fields Measurements by Means of Phase Data Map Calculation

Interference occurs when placing a quarter-wave plate and a linear polarizer at the output of an interferometer that has its reference and object beams with orthogonal linear polarization. Interference is between the right-handed and left-handed circularly polarized states in objects and reference beam, respectively. A desired phase shift will be obtained and controlled by the angle of the polarizer.^{13}^{,}^{14} This consideration is suitable in the implementation for a single-shot phase-shifting interferometer.

The optical setup presented for the optical phase calculation consists of a polarization Michelson interferometer (PMI) that is coupled to a $4\text{-}f$ system with an amplitude grating placed at the Fourier plane (Fig. 1). As a result, the interferogram first obtained by the PMI is replicated at the image plane of the $4\text{-}f$ system (Fig. 2).

## 2.1.

### Polarization Phase Shifting Technique

The output beam of the interferometer, after passing the quarter-wave plate, can be represented as the sum of the states, object and reference beam, with circular polarization states in opposite directions, respectively

## (1)

$${O}_{t}(x,y)=\frac{1}{\sqrt{2}}\left(\begin{array}{c}1\\ i\end{array}\right){e}^{-i\phi (x,y)}+\frac{1}{\sqrt{2}}\left(\begin{array}{c}1\\ -i\end{array}\right),$$## (2)

$${\mathbf{P}}_{\mathit{\theta}}=\left[\begin{array}{cc}{\mathrm{cos}}^{2}\text{\hspace{0.17em}}\theta & \mathrm{sin}\text{\hspace{0.17em}}\theta \text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\theta \\ \mathrm{sin}\text{\hspace{0.17em}}\theta \text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\theta & {\mathrm{sin}}^{2}\text{\hspace{0.17em}}\theta \end{array}\right].$$The resulting irradiance is associated with the interference pattern with a phase shift introduced twice the angle of the polarizers.

## 2.2.

### Interferogram Pattern Replication

The transmittance of a Ronchi grating of spatial period $d=(\lambda f)/({X}_{0})$ that is placed in the Fourier Plane can be written as

## (4)

$$G(u,v)=\mathrm{rect}\left[\frac{u}{{a}_{w}}\right]*{\sum}_{n=-\infty}^{\infty}\delta (u-n\xb7d)={\sum}_{n=-\infty}^{\infty}\mathrm{rect}\left[\frac{u-n\xb7d}{{a}_{w}}\right],$$## (5)

$$g(x,y)={\sum}_{n=-\infty}^{\infty}\frac{{a}_{w}}{d}\mathrm{sinc}\left(\frac{{a}_{w}}{\lambda f}x\right)\xb7\delta (x-n{X}_{0})\phantom{\rule{0ex}{0ex}}={\sum}_{n=-\infty}^{\infty}{C}_{n}\xb7\delta (x-n{X}_{0}).$$The physical meaning of the parameter ${X}_{0}$ is the replica separation in the image plane where its amplitude is modulated by a $\mathrm{sinc}(\xb7)$ function. Figure 3 shows the transmittance of the Ronchi ruling with its corresponding Fourier transform. For purpose of visualization in Fig. 3(b), the correspondent $\mathrm{sinc}(\xb7)$ function was plotted showing the amplitude change obtained by the zeroth order with the $\pm 1$st order.

The transmittance obtained in the image plane will be composed of the convolution of the interference pattern with the Ronchi ruling grating diffraction pattern scaled in the $\lambda $ and the focal length $f$.

## (6)

$${I}_{p}(x,y)=g(x,y)*{I}_{\theta}(x,y)=\sum _{n=-N}^{N}|{C}_{n}|\xb7{I}_{\theta}(x-n\xb7{X}_{0},y)\phantom{\rule{0ex}{0ex}}{I}_{n}(x,y)=|{C}_{n}|\{1+\mathrm{cos}[2\theta -\phi (x,y)]\},$$The major drawback of this system is the amplitude interferogram modulation caused by the phase/amplitude gratings.^{12} A more detailed study of the replicas obtained by using these techniques can be encountered in Ref. 2. This amplitude modulation can be solved by the implementation of fringe pattern normalization algorithms.^{15}^{,}^{16} As a difference from the system previously presented in Ref. 2, in this work, we used only one Ronchi ruling, taking the interference replicas obtained in the [$-1,0,1$] orders and using a fringe pattern normalization algorithm to avoid the interference amplitude difference against the [$-1,1$] orders with the [0] order.

## 3.

## Temperature Field Measurement

The phase difference $\phi (x,y)$ between two waves passing through the same point of the phase object, one in the presence of an inhomogeneous medium and the other in the air, is given by^{17}18.19.^{–}^{20}

## (7)

$$\phi (x,y)=\frac{2\pi}{\lambda}\underset{0}{\overset{\mathrm{\Gamma}}{\int}}[n(x,y,z)-{n}_{0}]\mathrm{d}z,$$^{21}By inverting Eq. (8), we are able to calculate the position dependent index of refraction. Taking into account the Gladstone–Dale relation, $(T)=1+K\rho $, expressed under ideal gas assumption, $\rho =p/RT$, we are able to retrieve the temperature profile as taking air as a reference at $T=15\xb0\mathrm{C}$ with a pressure $P=101325\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{N}/{\mathrm{m}}^{2}$, ${\rho}_{\mathrm{ref}}=1.225$, $R=287\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{J}/\mathrm{kg}\text{\hspace{0.17em}}\mathrm{K}$ as the specific gas constant, and ${K}_{G-D}=0.226\times {10}^{-3}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{m}}^{3}/\mathrm{Kg}$ as the Gladstone–Dale constant at $\lambda =632.8\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{nm}$.

^{22}Several studies have concluded that the assumption of air is sufficient for premixed flames as a sample, but in a more strictly manner, the Gladstone–Dale and specific gas constants need to be taken into account for specific local component distribution, not in the scope of this work.

^{23}24.

^{–}

^{25}

## 4.

## Experimental Setup of the Temperature Field Measurement

The proposed system is based on a PMI coupled to a $4\text{-}f$ system in order to obtain replicas of the interference pattern at the output (Fig. 1). Light from a He-Ne laser at $\lambda =633\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{nm}$ is linearly polarized at 45 deg and divided in amplitude by the beam splitter (BS). Linear polarizers at 0 and 90 deg are placed on the object and reference arms of the PMI, respectively. A quarter-wave plate (QWP) at 45 deg at the output of the interferometer turns linearly polarized light into left-handed and right-handed circularly polarized light, respectively. If we place a linear polarizer at this stage, we could obtain an interference pattern with its phase modulation controlled twice the linear polarizer angle. By placing an amplitude grating with a spatial frequency of $100\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{lines}/\mathrm{mm}$ on the frequency plane of the $4\text{-}f$ system, we are able to obtain replicas with the same capability of phase modulation by placing a linear polarizer angle in each of the replicas obtained (Fig. 2). As a result, we obtain the necessary interference fringe patterns with a relative phase shift in order to retrieve the phase data map properly. After each interferogram is separated from the image, the fringe processing is the same as any other algorithms already implemented in the literature. An He-Ne laser is utilized with a power of 20 mW and $\lambda =0.633\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu \mathrm{m}$ allowing sufficient illumination intensity to carry out the experiment. The monochromatic camera used is based on a CMOS sensor with $1280\times 1024\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{pixels}$ and with a pixel pitch of 6.7 *μ*m. The polarizing filter film that was used is commercially available (Edmund Optics, TechSpecs high contrast linear polarizing film). The size of each polarized film used in each interferogram was $25\times 25\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$, with its transmission axes oriented at 0, 60, and 120 deg. The wrapped phase data map was retrieved using the minimum least square algorithm of three steps,^{26} and the unwrapped phase data through the 2-D Goldstein branch cut phase unwrapping algorithm.^{27}

A commercial disposable lighter using butane as a fuel was used. The lighter was settled to obtain the minimum flame size possible and prepared to fix the amount of fuel flow during the experiment. The variation observed in the experiment corresponds to the temperature distribution in the region surrounding the flame. Temperature profile measurement was obtained analyzing the dynamic behavior at 8.3 fps in a period of time of 12 s.

Figure 4 presents the three phase-shifted interferograms obtained by a single capture. After a previous calibration procedure,^{1} we are capable of retrieving the three phase-shifted interference patterns. By implementing a fringe normalization technique, we are capable of retrieving the phase data map [see Fig. 4(b)]. By taking a reference phase map, the fringe order localization is directly obtained from the deformed fringe pattern, in this case, caused by a thin flame (Video 1).

With the fringe order localization information, we are able to solve the inverse Abel transform [Eq. (8)] and obtain the temperature change occurred by the thin-flame analyzed. Figure 5 shows the typical temperature profile obtained in a single capture.

As the system is able to obtain the dynamic temperature profile in Fig. 6, we present several captures of a dynamic distribution of air temperature that occurred due to environment changes outside of the thin flame (Video 2).

## 5.

## Final Remarks

The experimental setup for dynamic temperature field measurements using a single-shot polarization phase shifting technique has been described. This system is able to obtain three phase shifted interferograms in only one shot. Therefore, it is suitable to carry out temporal measurements of temperatures changes occurred outside of a thin flame. The implementations of normalization fringe procedures present the advantage of avoiding the use of phase gratings, and also the advantage of using only three interferograms for the analysis. The system is considerably simpler than previous proposals, showing a suitable alternative to implement in an industrial setting. It is important to note that for this work a polarization change related to the component’s temperature variation were not encountered. In order to apply this system for higher temperatures, it is strictly necessary the isolation of the sample without losing its polarization characteristics of the object beam for correct implementations.

## Acknowledgments

Authors thank M. A. Ruiz (CIO), Boaz Jessie Jackin (Utsunomiya Univ. CORE), and Barry Cense (Utsunomiya Univ. CORE) for their contribution in proofreading the manuscript. The authors gratefully acknowledge the research support from the CONACYT under grant 180449. D-I Serrano-García is currently occupying a position as a visiting researcher at Utsunomiya University in the Center for Optical Research & Education (CORE). D-I Serrano-Garcia is very grateful to CONACYT for the graduate scholarship granted (234454-227470/31458) and expresses sincere appreciation to Geliztle for the support provided. N.-I. Toto-Arellano expresses sincere appreciation to Luisa, Miguel, and Damian for the support provided, and to “Programa de Mejoramiento del Profesorado” (PROMEP) for Grant UTTGO-PTC-016.

## References

## Biography

**David Ignacio Serrano-García** received his bachelor’s degree from the Monterrey Institute of Technology and Higher Education in physical engineering in 2008 and his master’s degree in optics in 2010 from the Centro de Investigaciones en Optical A.C. in Leon, Guanajuato, Mexico. He is currently working toward his PhD in the same institution. His main research involves optical metrology and polarization phase shifting techniques applied to single-shot interferometry measurements.

**Amalia Martínez-García** received her PhD in optics at Centro de Investigaciones en Óptica, Mexico (CIO). She was a research scientist at Centro de Investigación Científica y de Educación Superior de Ensenada. Currently, she is a full-time researcher at CIO and member level II of the National System of Researchers, Mexico. She is the vice-president of the Mexican Academy of Optics during 2013–2014 and elected president 2015–2016. Her research area is in optical metrology.

**Noel-Ivan Toto-Arellano** received his BS degree in physical sciences at Universidad Veracruzana, México, his MSc in optical sciences from Benemerita Universidad de Puebla in 2005, and his PhD in optical sciences from the same institution in 2008. He did two postdoctoral internships at CIO. He currently holds a research professor position in Universidad Tecnológica de Tulancingo, Hidalgo, México, supporting the creation of a photonics engineering program at the same university.

**Yukitoshi Otani** is a professor at the Center for Optical Research and Education, Utsunomiya University, Japan. He received his doctorate from the University of Tokyo in 1995. He was an associate professor at Tokyo University of Agriculture and Technology until 2010 and visiting professor at the College of Optical Sciences, the University of Arizona, from 2004 to 2005. He is a fellow of SPIE since 2010.