Characterization of Supersonic Compressible Fluid Flow Using High-Speed Interferometry

Abstract

This paper presents a very effective interference technique for the sensing and researching of compressible fluid flow in a wind tunnel facility. The developed technique is very sensitive and accurate, yet easy to use under conditions typical for aerodynamic labs, and will be used for the nonintrusive investigation of flutter in blade cascades. The interferometer employs a high-speed camera, fiber optics, and available “of-the-shelf” optics and optomechanics. The construction of the interferometer together with the fiber optics ensures the high compactness and portability of the system. Moreover, single-shot quantitative data processing based on introducing a spatial carrier frequency and Fourier analysis allows for almost real-time quantitative processing. As a validation case, the interferometric system was successfully applied in the research of supersonic compressible fluid discharge from a narrow channel in a wind tunnel. Density distributions were quantitatively analyzed with the spatial resolution of about $$50 \mathsf{\mu }\mathrm{m}$$. The results of the measurement revealed important features of the flow pattern. Moreover, the measurement results were compared with Computational Fluid Dynamics (CFD) simulations with a good agreement.

1. Introduction

High-speed, high-accuracy, whole-field, non-destructive and contactless optical methods are increasingly being used in the study of fluid flow, offering the possibility for deeper insights than single-point conventional techniques such as constant-current anemometry [1], pneumatic measurements [2], thermal flow sensors (hot film sensors, calorimetric sensors, time-of-flight sensors) [3,4,5,6] and others. Moreover, contactless optical methods do not disturb the flow field. Due to the presence of density gradients, compressible fluid flow is very suitable for an investigation using optical methods that are sensitive either to the refractive index of the medium through which the light wave propagates (interferometry), or to a refractive index gradient in the medium (schlieren [7,8,9,10,11], shadowgraph [12,13,14,15,16], Moiré deflectometry [17,18,19,20,21]). Interferometric methods such as classical interferometry [22,23], holographic interferometry [24] and digital holographic interferometry [25,26,27,28,29] stand as the most accurate optical techniques [26,30,31] that are applied in many sectors of industry and scientific research, notably in the study of liquid cooling [32], diffusion [33], convection [34,35], temperature measurements and visualization [1,25,36], temperature fields measurement in pulsatile jets [25], studies of living cell imaging [37,38], precision measurement [15,39,40], vibrometry [41,42] etc.

However, the field of the transonic or even supersonic flow generated in wind tunnel facilities is very specific, and the current state of the art of interferometric techniques in this field allows only their limited use. Environmental harshness, particularly vibrations, the lack of an optical table, temporal flow instability, and high fluid density gradients leading to complex interference patterns often make the quantitative evaluation of interferometric data very difficult and cumbersome. This is why the use of brute force interferogram analysis [43] or capturing interferograms for qualitative analysis and flow visualization in wind tunnels is still common [22,44,45,46,47,48,49].

In this paper, we introduce a robust, portable, high-speed and high-resolution interferometric method that can be easily applied in a wind tunnel facility, providing very accurate results with high spatial resolution in almost real time. The method utilizes a high-speed camera allowing for short exposures to avoid the blur of interference patterns typical for unstable flow, as well as to minimize environmental disturbances. Fiber optics ensure the high compactness and portability of the interferometric system, and fast single-shot quantitative data processing based on spatial carrier frequency and Fourier analysis is employed.

The method is intended for the investigation of high-speed flow in planar blade cascades, namely, for studying dynamic events during transonic and supersonic blade flutter. As a first step, it was successfully tested in the case of a supersonic highly underexpanded flow at the outlet of a narrow channel by measuring the density distribution of the flow. The present interferometric technique has proven to be an excellent tool to meet these objectives. In addition to a description of the principles, experimental arrangement and results, numerical simulations showing the high reliability of the measured data are presented.

2. Method

2.1. Interferometry and Fluid Flow

Interferometric techniques generally provide whole-field, non-invasive and highly accurate measurements. These properties can be advantageously used for measuring and visualizing such flow phenomena that affect light wave [25]. The amplitude of the propagating light wave through the transparent/semi-transparent medium under investigation is not significantly affected, while the optical phase is [50]. These phenomena are called phase-sensitive, implying that when a change takes place in the medium, it affects the refractive index of the medium $\Delta n$, which as a consequence introduces phase deviation $\Delta \phi$:

(1)$\Delta \phi \left(x,y\right)=\frac{2\pi }{\lambda }{{\displaystyle \int }}_{L_1}^{L_2}\Delta n\left(x,y,z\right)dl,$

with $dl$ denoting the differential distance along the optical path between points $L_1$ and $L_2$, and $\lambda$ stands for the wavelength of the laser light. The $x$,$y$ coordinates define the position within the Cartesian coordinates system. The general solution of (1) leads to tomographic techniques [51] that are difficult or cumbersome to apply in real-time, dynamic measurements. However, there are two assumptions that allow a significant simplification of (1): (i) wind tunnel experiments are designed in order to keep the refractive index along the $z$-axis constant, i.e., the optical axis, and (ii) the bending of the rays due to a spatial variation in the refractive index can be neglected. Assuming (i) and (ii), Equation (1) can be simplified to:

(2)$\Delta \phi \left(x,y\right)=\frac{2\pi }{\lambda }\Delta n\left(x,y\right) L,$

where $L$ denotes the length of the object in the $z$-axis direction (i.e., wind tunnel testing section width). The refractive index change distribution is usually not the main concern. An important equation in the optical study of fluid flow is the Gladstone–Dale equation, which links the change of refractive index $\Delta n$ to the density $\rho$ as:

(3)$\rho \left(x,y\right)={\rho }_{ref}+\frac{\Delta n\left(x,y\right)}{K},$

where $K$ represents the Gladstone–Dale constant and ${\rho }_{ref}$ is the reference density. The Gladstone–Dale equation is almost independent of pressure or temperature under moderate physical conditions [50]. For a case of an isentropic compressible flow of ideal gas, the measured outlet density $\rho$ can be further linked with the isentropic Mach number by the following equation:

(4)$M_i\left(x,y\right)=\sqrt{\frac{2}{\kappa -1}·\left({\left(\frac{{\rho }_{in}}{\rho \left(x,y\right)}\right)}^{\kappa -1}-1\right)},$

where ${\rho }_{in}$ is the inlet air density and $\kappa =1.4$ denotes the heat capacity ratio of an ideal diatomic gas. In addition, other quantities, such as velocity distribution, pressure distribution, temperature distribution, etc., can be analyzed quantitatively by interferometry.

The investigation within this paper leads to the characterization of the air flow in/behind a narrow channel (see Figure 1). Air enters the test section at ambient conditions, and thus the total pressure $p_{in}$ is equal to the barometric pressure. The static pressure at the outlet of the narrow channel $p_{out}$ is measured by a pressure tap on a wall in the settling chamber using a pressure transducer (from Huba Control AG, Würenlos, Switzerland). For the case of an isentropic compressible flow of ideal gas, the density distribution ${\rho }_M$ in the settling chamber can be calculated as:

(5)${\rho }_M\left(x,y\right)={\rho }_{ref}+\frac{\Delta n\left(x,y\right)}{K},$

where

(6)${\rho }_{ref}=\frac{p_{out}}{R{T}_0}.$

In (6), $T_0$ denotes the ambient air temperature and $R = 287.1 \mathrm{J}/\left(\mathrm{kg}·\mathrm{K}\right)$ is the specific gas constant. The isentropic Mach number (4), however, can provide relevant information only in flows where the drops in total pressure and total density throughout the flow field are relatively small, e.g., the flow in blade cascades. This does not apply to the investigated phenomenon, and therefore the isentropic Mach number was not evaluated.

2.2. Spatial Carrier Interferometry

Equations (2)–(6) provide the relation between the phase change $\Delta \phi$ and other quantities, such as density ${\rho }_M$; however, it is not clear yet how the phase change $\Delta \phi$ is obtained.

The interference phase $\phi$ can be retrieved from an interference pattern

(7)$I\left(x,y\right)=A\left(x,y\right)+B\left(x,y\right) \cos \left(\phi \left(x,y\right)\right)$

created on a charge-coupled device (CCD), or a complementary metal-oxide semiconductor (CMOS) chip, by the superposition of a reference and an object beam. In (7), $A$ is the additive component and $B$ is the multiplicative component. The interference phase $\phi$ is coded in a cosine modulate fringe pattern. There are several ways to retrieve the interference phase [52]; however, for the characterization of fluid flow in wind tunnels, where real-time and high-resolution measurements are required, most of the interferogram processing techniques fail. An option to tackle the issues is spatial carrier interferometry [52], which introduces small angles of incidence ${\theta }_x$, ${\theta }_y$ in the $x$, $y$ directions on the camera sensor in the reference beam. This results in further modulation of the interference pattern with cosine fringes pattern of constant spatial frequencies $k_{x0}=\frac{2\pi }{\lambda }sin{\theta }_x$ and $k_{y0}=\frac{2\pi }{\lambda }sin{\theta }_y$, respectively. Let us rewrite (7) using a complex exponential including the carrier frequencies $k_{x0}$ and $k_{y0}$:

(8)$I\left(x,y\right)=A\left(x,y\right)+1/2C\left(x,y\right)+1/2C^{*}\left(x,y\right),$

where $C\left(x,y\right)=B\left(x,y\right)\exp \left[i\left(\phi \left(x,y\right)+k_{x0}x+k_{y0}y\right)\right]$ and $C^{*}\left(x,y\right)=B\left(x,y\right)\exp \left[-i\left(\phi \left(x,y\right)+k_{x0}x+k_{y0}y\right)\right]$, with imaginary unit $i$. Equation (8) in the Fourier domain

(9)$\widehat{I}\left(k_x,k_y\right)=\widehat{A}\left(k_x,k_y\right)+\widehat{C}\left(k_x-k_{x0},k_y-k_{y0}\right)+{\widehat{C}}^{*}\left(k_x+k_{x0},k_y+k_{y0}\right)$

with the spatial frequency coordinates $k_x, k_y$ is composed of a central DC (direct current) term and two conjugated components located symmetrically from the center of the spectrum. The roof symbol $\widehat{}$ denotes Fourier spectrum and the superscript symbol * denotes complex conjugation. The introduced angles ${\theta }_x$, ${\theta }_y$ play a major role in the clear separation of all spectral components $\widehat{A}$,$\widehat{C}$ and ${\widehat{C}}^{*}$ in the Fourier domain and must be sufficiently high to avoid any overlap. On the other hand, angles cannot be too high, in order to meet the Nyquist criterion.

Once the spectral components are separated, using a bandpass filter around the specific spatial frequencies, it is possible to filter out the desired spectral component $\widehat{C}$ containing the interference phase $\phi$. By filtering only the desired component, the filtered spectrum ${\widehat{C}}_F$ (subscript F denotes filtered component) is no longer Hermitian, so its inverse Fourier transform will contain non-zero both real and imaginary parts, and hence the interference phase can be calculated as:

(10)$\phi \left(x,y\right)=\arctan \left(\frac{Im\left(C_F\left(x,y\right)\right)}{Re\left(C_F\left(x,y\right)\right)}\right).$

The phase along with other quantities is stored in a computer’s memory as a matrix of numbers. Let us assume a measurement at a steady/reference state without the presence of fluid flow. Quantities measured in the reference state can be denoted using a subscript $0$, e.g., ${\phi }_0$ for the interference phase. Such a measurement carries information about interferometer optical aberrations, including aberrations introduced by, e.g., optical windows of the wind tunnel experiment. Other measurements performed with the presence of a phenomenon (denoted with subscript 1) can be related to the reference state measurement in order to suppress undesired optical aberrations. The phase change $\Delta \phi$ occurring in (1) can be calculated as:

(11)$\Delta \phi \left(x,y\right)=\arctan \left(\frac{Im\left(C_{F1}\left(x,y\right)C^{*}{}_{F0}\left(x,y\right)\right)}{Re\left(C_{F1}\left(x,y\right)C^{*}{}_{F0}\left(x,y\right)\right)}\right).$

It is worth noting that spatial carrier interferometry could also be called off-axis image plane digital holographic interferometry. As phenomena in fluid flow do not lead to speckles, and there is also no need for the numerical propagation/focusing/reconstruction of wave fields typical for digital holography, we tend to call it interferometry.

3. Materials and Methods

Interferometry-based research of a fluid flow in a narrow channel including a jet flow at its outlet was carried out at the Laboratory of Internal Flows of the Institute of Thermomechanics of the Czech Academy of Sciences. The laboratory is equipped with a modular in-draft wind tunnel facility that enables low mass flow rates. The vacuum tank is connected to the test section through a pipe with a quick-acting valve and a nozzle that offers the possibility of changing the pressure inside the testing section from $0.01 \mathrm{MPa}$ to values close to those of the ambient atmospheric pressure. The three-dimensional computer-aided design (3D CAD) assembly of the narrow channel is shown in Figure 2.

The experimental arrangement was designed in order to maintain a constant refractive index along the $z$-axis, i.e., the optical axis. Two parallel walls made of stainless steel constitute the flow channel. The test section has optical windows on the walls, allowing the passage of light beams for optical measurements. The channel height (along the $y$-axis) was set to $3 \mathrm{mm}$. The channel is $100 \mathrm{mm}$ long ($x$-axis) and the width ($z$-axis) of the channel is $L=100 \mathrm{mm}$ to achieve a sufficient aspect ratio to eliminate secondary flow effects and to provide the sufficiently long optical path necessary for the adequate sensitivity of the interferometry. In order to obtain quantitative values of measurements, the ambient atmospheric pressure ($P_{in}$) as well as the outlet static pressure in the settling chamber of the wind tunnel ($P_{out}$) were measured.

The experimental setup, consisting of the narrow channel experimental section and a Mach–Zehnder interferometer, is illustrated in Figure 3. A laser beam was generated by a pigtailed single-frequency distributed feedback laser (LAS) of wavelength $\lambda =773 \mathrm{nm}$ and fiber output power of $38 \mathrm{mW}$. The laser beam was split using a fibersplitter (FS) into the reference wave and the object wave. The reference wave was guided by an optical fiber to the non-polarizing beamsplitter (NBS), while the free-space object wave passed through a collimating imaging lens L0, which sent light through the measuring area (MA) within the wind tunnel (WT), before it passed through the lens L1. Lenses L1 and L2 acted as a beam expander with magnification of $Mag\approx 0.6$ as well as an imaging system. Both reference and object waves were recombined by the NBS and collimated by L2. The object wave impinged the digital camera sensor (CAM) normally, while the reference wave angles were ${\theta }_x=0.34°, {\theta }_y=0.30°,$ respectively. The angles of the reference wave were adjusted by a slight translation of the reference beam’s fiber ferrule (FF).

The intensity interference pattern (interferogram) of both superposed beams was captured by a Phototron FASTCAM Mini WX100 high-speed camera with a resolution of $2048\times 2048$ pixels ($10 \mathsf{\mu }\mathrm{m}\times 10 \mathsf{\mu }\mathrm{m}$ pixel size) and a frame rate of $1080 \mathrm{fps}.$ The exposure time of $10 \mathsf{\mu }\mathrm{s}$ was set in order to avoid the blurring of the interference pattern while maintaining sufficient brightness (signal).

In the first step, a reference interferogram at a steady state ($P_{out}=P_{in}$) was captured and processed, and the filtered complex field $C_{F0}$ was stored in a computer. The same procedure was repeated after a change in the pressure inside the testing section ($P_{out}\ne P_{in}$), resulting in air flow. In order to obtain enough data for statistical analysis, a sequence of $M=50$ interferograms was captured, of which one is presented in Figure 4a. The amplitude Fourier spectrum of the interferogram in the logarithmic scale is shown in Figure 4b. The terms $\widehat{A,} \widehat{C,} {\widehat{C}}^{*}$ are well separated by the spatial carrier frequencies $f_{x0}=\frac{k_{x0}}{2\pi }=\frac{sin{\theta }_x}{\lambda M}=12.3 {\mathrm{mm}}^{-1}$ and $f_{y0}=10.8{ \mathrm{mm}}^{-1},$ respectively. The black circle represents the filtering bandpass Hanning window, with bandwidth $f_B=20{ \mathrm{mm}}^{-1}$, applied in order to retrieve the complex field $C_{Fm}$, where $m=1,2\dots M$ is an integer denoting the frame number and M is the number of captured interferograms. Such bandwidth results in a spatial (lateral) resolution of about $50 \mathsf{\mu }\mathrm{m}$. It is important to note that the applied filtering must remain unchanged from that in the reference interferogram processing. The optical phase change computed by (11) is wrapped within the $\left[-\pi , \pi \right] radians$ interval due to the harmonic nature of the light wave—see Figure 4c. Therefore, a spatial unwrapping was applied with the starting point in the position of the Pitot probe, i.e., far from the channel without the influence of the flow, obtaining a sequence of unwrapped phase fields $\Delta {\phi }_m$—see Figure 4d. In order to suppress the random time fluctuations of the flow, the mean of the unwrapped phase was computed:

(12)$\Delta \phi =\frac{{{\displaystyle \sum }}_{m=1}^M\Delta {\phi }_m}{M}$

The mean unwrapped phase change distribution $\Delta \phi$ was used to calculate the density distribution as described in Section 2.1.

In addition to the interferometric measurements, computational fluid dynamics (CFD) simulation software was used to cross-check both the measurement and the simulation. The numerical computations were performed using Ansys Fluent 2021R1 commercial code. This software package uses the finite volume method to solve the governing Navier–Stokes equations for compressible fluid flow. The problem was simulated in two dimensions (2D) using the Reynolds-averaged Navier–Stokes equations (RANS) approach to close the system of equations. The flow of air as an ideal gas was simulated as fully turbulent using the two-equation k-ω shear stress transport (SST) turbulence model. The computational mesh consisted of approximately $\mathrm{150,000}$ cells. The mesh sensitivity analysis was performed using two finer meshes, yielding a relative error of the mass flow rate of $0.3\%$ for the coarse mesh, which is very acceptable. The computational domain was assumed to be symmetrical around the channel axis with the no-slip boundary condition at the channel walls. The static pressure according to the desired pressure drop across the channel was prescribed as the outlet boundary condition. The following boundary conditions were used at the inlet: (a) constant total temperature $T_0=292.67 \mathrm{K},$ (b) uniform total pressure distribution $p_0=98.05 \mathrm{kPa},$ (c) constant turbulence intensity $Tu=1.2\%$ and (d) constant eddy viscosity ratio $\frac{{\mu }_t}{\mu }=10$.

4. Results and Discussion

The results of the interferometric measurement of air flowing out of the narrow channel, and their comparison to the results of CFD simulation, are presented in this section.

The fluid flow in the narrow channel was investigated in five different regimes defined by different pressure gradients through the channel. The magnitude of the pressure gradient was controlled by a valve that determined the outlet pressure in the settling chamber downstream of the narrow channel. Both inlet and outlet pressure were measured, yielding the pressure ratio $\eta =\frac{P_{out}}{P_{in}}$. The inlet pressure (the ambient atmospheric pressure) was constant during all measurements: $P_{in}=98.050 \mathrm{kPa}$. Parameters under different regimes are summarized in Table 1.

Other constants considered for further flow analysis and CFD simulations are the ambient temperature $T_0=295.9 \mathrm{K}$, the Gladstone–Dale constant $K=0.000225 {\mathrm{m}}^3{\mathrm{kg}}^{-1}$, and the specific gas constant of dry air $R=287.1 \mathrm{J}/\left(\mathrm{kg}·\mathrm{K}\right)$.

The density distributions computed using (5), together with the results of the CFD simulation, are presented in Figure 5. It can be seen that the interferometric measurements and CFD simulations are in a very good agreement. Small discrepancies can be attributed to minor asymmetries in the geometry of the experimental setup. The boundary conditions of the CFD simulation could be modified in order to derive an even better agreement; however, this is not this paper’s objective.

The measured distributions in Figure 5 have significantly dynamic ranges of values. Therefore, some faint details of the flow might be difficult to observe. In order to visualize the high-frequency spatial variations in the flow, the slopes of the measured density distribution in the horizontal $\left(x\right)$ and vertical $\left(y\right)$ directions $\partial {\rho }_M\left(x,y\right)/\partial x$,$\partial {\rho }_M\left(x,y\right)/\partial y$ were calculated, and these, together with the magnitude of the slope $\sqrt{{\left(\partial {\rho }_M/\partial x\right)}^2+{\left(\partial {\rho }_M/\partial y\right)}^2}$, are shown in Figure 6. Such results could be directly compared to the Schlieren imaging techniques, which are sensitive to the ray deflections caused by refractive index gradients $\partial n\left(x,y\right)/\partial x$ and $\partial n\left(x,y\right)/\partial y$, respectively.

Five different regimes with a gradually decreasing ratio $\eta$ were measured. In the first regime, #1 ($\eta =0.560$), the supercritical pressure ratio across the channel results in the subsonic discharge of air into the settling chamber. Thus, a subsonic exit jet can be seen in Figure 6. There is also visible compression just before the exit of the channel. This indicates that the regime is close to aerodynamic choking, and is approaching supersonic velocity within the channel.

In regime #2 ($\eta =0.481$), the exit jet is supersonic with a shock wave present. Thus, the flow is already aerodynamically choked. In this regime, the flow is most likely overexpanded (i.e., the pressure at the exit cross-section of the channel is lower than the outlet pressure in the settling chamber), since two oblique shocks are formed at the exit of the channel first.

In the supersonic regime #3, as the ratio further decreases ($\eta =0.400$), the jet flow’s inner pattern starts to take an “X” or a “diamond” shape, which is characteristic of a moderately underexpanded jet flow. Here, expansion waves form at the edge of the channel, causing the flow to turn outward in a plume, creating a characteristic underexpanded jet flow pattern. Expansion waves “reflect” from the jet’s constant pressure boundary as compression waves, which eventually form oblique shock waves. A periodical series of compression and expansion waves appears, until the jet is dissipated and static pressure within the jet flow reaches equilibrium with the surrounding pressure (backpressure).

In the supersonic regimes #4 ($\eta =0.287$) and #5 ($\eta =0.243$), the flow first expands through expansion waves, before turning into oblique shocks and Mach disks. Such highly underexpanded flow regimes exhibit a “barrel-shaped” pattern that differs in the shapes and lengths of the shock cells with different $\eta$. The flow patterns repeat until mixing dissipates the flow.

The typical resolution of the optical phase measured by interferometry is better than $1/100$ of $2\pi$, leading to a refractive index resolution of about $4\times {10}^{-7}$ in this particular case. However, interferometric measurements in a wind tunnel are burdened with errors. The first group of errors is directly due to interferometric sensing—the measurement of the optical phase. When retrieving $\Delta n$ from (2), there is an error in laser wavelength $u_{\lambda }$, thickness $u_L$, and phase noise $u_{\Delta \phi }$. Let us define the refractive index combined error as:

(13)$u_{\Delta n}=\left|\frac{\partial \Delta n}{\partial \Delta \phi }\right|u_{\Delta \phi }+\left|\frac{\partial \Delta n}{\partial \lambda }\right|u_{\lambda }+\left|\frac{\partial \Delta n}{\partial L}\right|u_L$

The phase noise for our measurement is estimated to be $u_{\Delta \phi }=0.07 \mathrm{rad}$. This value was derived from experimental data as the standard deviation of an unwrapped phase map (in different regions) after low-order polynomial subtraction (i.e., high-spatial frequency components were considered as noise). The wavelength of the used laser source depends on the temperature and current of the laser diode. The temperature and the current are measured/controlled; however, there is not a specific relation between these and the output wavelength, thus there is a wavelength error, which can be estimated to be $u_{\lambda }=0.5 \mathrm{nm}$. The error of thickness, $u_L=0.1 \mathrm{mm}$, comprises, on one hand, the limited accuracy of channel manufacturing, and on the other hand the fact that rays are deflected due to refractive index non-homogeneities. Such ray deflection is not considered in (2). Moreover, in locations with a large refractive index gradient, the ray deflection misplaces information geometrically (up to $200 \mathsf{\mu }\mathrm{m}$), leading to so-called mapping errors. Such distortion is below ten pixels in the areas with the largest gradients, otherwise it is negligible. By substituting in (13), the maximal error is $u_{\Delta n}\approx 1\times {10}^{-6}$ and the relative maximal error results in $100\frac{u_{\Delta n}}{\Delta n}\approx 2\%$.

Assuming the refractive index error $u_{\Delta n}$, along with the $2\%$ relative error of the pressure/temperature sensors and $1\%$ error of the Gladstone–Dale constant in (5), the density error yields $u_{{\rho }_M}\approx 0.015 \mathrm{kg}/{\mathrm{m}}^3$, or relatively $100\frac{u_{{\rho }_M}}{{\rho }_M}\approx 2.6\%$.

Another source of error is usually related to the repeatability of measurements. However, the repeatability of measurements in this particular case can only be determined under a steady state (without flow) due to the unstable behavior of the flow. Steady state repeatability has a negligible effect when compared to the stability of the flow. Therefore, time stability was examined. As aforementioned, we captured $50$ frames during one measurement sequence, taking $50 \mathrm{ms}$, and used averaged phase maps $\Delta \phi$ for further processing—see (12). However, the standard deviation,

(14)$\sigma \left(x,y\right)=\sqrt{\frac{{{\displaystyle \sum }}_{m=1}^M{\left(\Delta {\phi }_m\left(x,y\right)-\Delta \phi \left(x,y\right)\right)}^2}{M-1}},$

for each pixel of $\left(x,y\right)$ coordinates can also be computed from the data set. The standard deviation can be used to calculate a time stability map $\tau \left(x,y\right)==1-3\sigma \left(x,y\right)$ that (assuming the normal probability distribution) equals one in pixels of completely stable flow, while it approaches zero at locations of random signal. The time stability maps of flow in different flow regimes are shown in Figure 7. It is clear that in the narrow channel, the flow is very stable, but becomes much less stable in shear layers and regions of flow dissipation. This effect is more apparent at high flow velocities with a more turbulent character.

5. Conclusions

This paper presents a very effective technique for sensing and investigating compressible fluid flow. The developed interferometric technique is very sensitive and accurate, yet easy to use, even under pretty harsh environmental conditions (vibrations, flow instability, unclean room) typical for aerodynamic labs. The interferometer employs a high-speed camera, fiber optics, and readily available “of-the-shelf” optics and optomechanics. The interferometer consists of an illumination and a sensing unit, each weighing less than $5 \mathrm{kg}$. An optical fiber delivers the reference wave from the illumination to the sensing unit, while the object wave emitted from the illumination unit propagates through a measured volume into the sensing unit. Both waves are superimposed, and an interference pattern is captured by a high-speed camera. The construction of the interferometer together with the fiber optics ensures the high compactness and portability of the system. Moreover, single-shot quantitative data processing based on introducing a spatial carrier frequency and processing in the Fourier domain allows for almost real-time quantitative processing, which is also suitable for very fast-evolving phenomena.

The interferometer was successfully applied for the study of a fluid flow in a narrow channel including a free jet at the channel exit. The experimental section was built in a wind tunnel, allowing for a high pressure difference and thus high flow velocity through the channel. The flow was investigated under five different regimes, ranging from subsonic to supersonic flow with a highly underexpanded jet. The density distributions were quantitatively analyzed with a spatial resolution of about $50\mu m$. The results of the measurement were able to reveal important features of the flow pattern, such as shock waves and slip lines.

The measurement results were compared with CFD simulations, showing good agreement in terms of the flow pattern similarity (e.g., size of the shock cells) as well as the absolute values of density. An error analysis revealed the error range in units of percent. It can be concluded that the method constitutes a promising tool for studying highly dynamic events, even during transonic and supersonic blade fluttering.

Author Contributions

Methodology, J.H., G.Ç., V.L., M.J. and P.P.; software, J.K., M.S. and P.P.; wind tunnel experiments, P.P., G.Ç., J.H., D.Š., P.Š., J.L. and M.L.; data processing, P.P., G.Ç., J.K. and M.S.; writing—P.P., G.Ç., P.Š., D.Š. and J.H.; review and editing, J.L., M.L. and V.L.; funding acquisition, J.L., M.L., P.Š. and P.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Ministry of Education, Youth and Sports of the Czech Republic—program Inter-Excellence, project No. LTAUSA19036. J.K., V.L. and M.S. gratefully acknowledge the support of the Grant Program of the Technical University of Liberec (project number PURE-2020-3010). The institutional support RVO:61388998 of the Institute of Thermomechanics of the Czech Academy of Sciences is also acknowledged.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

We would like to thank Jakub Kozohorský and Antonín Zajíček for their help in setting up the experiments in the wind tunnel.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures and Table

Figure 1. Basic concept of the narrow channel measurement. Dashed rectangle represents the investigated area.

Figure 2. 3D CAD assembly of the narrow channel facility.

Figure 3. Scheme of the setup: LAS—laser, FS—fiber splitter, FF—fiber ferrule, L0 + L1 + L2—lens, MA—measured area, P0 + P1—pressure probes, WT—wind tunnel, NBS—beamsplitter, CAM—camera.

Figure 4. (a) Intensity image (interferogram) captured during measurements, with the zoomed part showing the interference pattern, (b) its Fourier spectrum, (c) the phase change map with region of interest (ROI) marked by the dashed line, and (d) unwrapped phase change map with ROI.

Figure 5. (a) Density distributions for different regimes—comparison of measured data (left) and the CFD simulations (right); (b) density values along the profile denoted in (a) for regime #5.

Figure 6. Slope maps of measured (top—M) and simulated (bottom—S) density distributions in different flow regimes.

Figure 7. The time stability of the flow under different regimes.

Table 1

Regimes and the corresponding pressure and pressure ratio values.

Regime	#1	#2	#3	#4	#5
Pout (kPa)	54.923	47.140	39.236	28.114	23.871
η [1]	0.560	0.481	0.400	0.287	0.243

References

1. Psota, P.; Dančová, P.; Cubreli, G.; Lédl, V.; Vít, T.; Doleček, R.; Matoušek, O. Development and application of spatial carrier interferometry for whole field real-time investigation of temperatures in liquid media. Int. J. Therm. Sci.; 2019; 145, 106029. [DOI: https://dx.doi.org/10.1016/j.ijthermalsci.2019.106029]

2. Kim, K. Pneumatic Measurements for Pressure, Velocity, and Flow-direction. Appl. Fluid. Meas. Tech.; 2016; pp. 61-100. [DOI: https://dx.doi.org/10.1016/b978-0-12-809731-1.00003-4]

3. Nour, M.A.; Hussain, M.M. A Review of the Real-Time Monitoring of Fluid-Properties in Tubular Architectures for Industrial Applications. Sensors; 2020; 20, 3907. [DOI: https://dx.doi.org/10.3390/s20143907] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32674278]

4. Berthet, H.; Jundt, J.; Durivault, J.; Mercier, B.; Angelescu, D. Time-of-flight thermal flowrate sensor for lab-on-chip applications. Lab Chip; 2011; 11, pp. 215-223. [DOI: https://dx.doi.org/10.1039/C0LC00229A] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21072440]

5. Meng, E.; Li, P.-Y.; Tai, Y.-C. A biocompatible Parylene thermal flow sensing array. Sens. Actuators A Phys.; 2008; 144, pp. 18-28. [DOI: https://dx.doi.org/10.1016/j.sna.2007.12.010]

6. Kuo, J.T.W.; Yu, L.; Meng, E. Micromachined Thermal Flow Sensors—A Review. Micromachines; 2012; 3, pp. 550-573. [DOI: https://dx.doi.org/10.3390/mi3030550]

7. Trolinger, J.D.; Buckner, B.; L’Esperance, D. Background-oriented schlieren for the study of large flow fields. Opt. Eng. Appl.; 2015; 9576, 95760.

8. Saravanan, S.; Nagashetty, K.; Hegde, G.M.; Jagadeesh, G.; Reddy, K.P.J. Schlieren Visualization of Shock Wave Phenomena over a Missile-Shaped Body at Hypersonic Mach Numbers. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng.; 2011; 225, pp. 26-34. [DOI: https://dx.doi.org/10.1243/09544100JAERO827]

9. Kleine, H.; Hiraki, K.; Maruyama, H.; Hayashida, T.; Yonai, J.; Kitamura, K.; Kondo, Y.; Etoh, T.G. High-speed time-resolved color schlieren visualization of shock wave phenomena. Shock Waves; 2005; 14, pp. 333-341. [DOI: https://dx.doi.org/10.1007/s00193-005-0273-6]

10. Nel, L.; Skews, B.; Naidoo, K. Schlieren techniques for the visualization of an expansion fan/shock wave interaction. J. Vis.; 2014; 18, pp. 469-479. [DOI: https://dx.doi.org/10.1007/s12650-014-0253-y]

11. Passmann, M.; Wiesche, S.A.D.A.D.; Joos, F. Focusing Schlieren Visualization of Transonic Turbine Tip-Leakage Flows. Int. J. Turbomach. Propuls. Power; 2020; 5, 1. [DOI: https://dx.doi.org/10.3390/ijtpp5010001]

12. Aftab, S.; Younis, O.; Al-Atabi, M. Four Decades of Utilizing Shadowgraph Techniques to study Natural Convection in Cavities: Literature Review. IOP Conf. Ser. Mater. Sci. Eng.; 2012; 36, 012021. [DOI: https://dx.doi.org/10.1088/1757-899X/36/1/012021]

13. Srivastava, A.; Phukan, A.; Panigrahi, P.; Muralidhar, K. Imaging of a convective field in a rectangular cavity using interferometry, schlieren and shadowgraph. Opt. Lasers Eng.; 2004; 42, pp. 469-485. [DOI: https://dx.doi.org/10.1016/j.optlaseng.2004.03.003]

14. Hijikuro, M.; Anyoji, M. Application of Optical Flow Analysis to Shadowgraph Images of Impinging Jet. J. Flow Control. Meas. Vis.; 2020; 8, pp. 173-187. [DOI: https://dx.doi.org/10.4236/jfcmv.2020.84011]

15. Waynant, R.W.; Ediger, M.N. Electro-Optics Handbook; McGraw-Hill Education: New York, NY, USA, 2000.

16. Xu, F.; Patterson, J.C.; Lei, C. Shadowgraph observations of the transition of the thermal boundary layer in a side-heated cavity. Exp. Fluids; 2005; 38, pp. 770-779. [DOI: https://dx.doi.org/10.1007/s00348-005-0960-1]

17. Juste, G.L.; Benavides, E.M. Temperature Measurement in Small-Scale Flows with Digital Moiré Deflectometry. Exp. Heat Transf.; 2011; 24, pp. 201-214. [DOI: https://dx.doi.org/10.1080/08916152.2010.503312]

18. Juste, G.; Fajardo, P. Assessment of experimental optical techniques for characterizing heat transfer using numerical simulations. Eng. Appl. Comput. Fluid Mech.; 2015; 9, pp. 1-15. [DOI: https://dx.doi.org/10.1080/19942060.2015.1004812]

19. Stricker, J.; Keren, E.; Kafri, O. Axisymmetric density field measurements by moire deflectometry. AIAA J.; 1983; 21, pp. 1767-1769. [DOI: https://dx.doi.org/10.2514/3.8326]

20. Rasouli, S. Atmospheric Turbulence Characterization and Wavefront Sensing by Means of the Moiré Deflectometry. Topics in Adaptive Optics; Intech: London, UK, 2012; 23. [DOI: https://dx.doi.org/10.5772/31129]

21. Applications of 2-D Moiré Deflectometry to Atmospheric Turbulence. J. Appl. Fluid Mech.; 2014; 7, pp. 651-657. [DOI: https://dx.doi.org/10.36884/jafm.7.04.21420]

22. Luxa, M.; Příhoda, J.; Šimurda, D.; Straka, P.; Synáč, J. Investigation of the compressible flow through the tip-section turbine blade cascade with supersonic inlet. J. Therm. Sci.; 2016; 25, pp. 138-144. [DOI: https://dx.doi.org/10.1007/s11630-016-0844-0]

23. Borjian, E.; Yousefi, T.; Ashjaee, M. Optical interferometry to investigate the heat transfer from a vertical cone under air jet impingement. Opt. Lasers Eng.; 2015; 67, pp. 205-211. [DOI: https://dx.doi.org/10.1016/j.optlaseng.2014.12.002]

24. Colombani, J.; Bert, J. Holographic interferometry for the study of liquids. J. Mol. Liq.; 2007; 134, pp. 8-14. [DOI: https://dx.doi.org/10.1016/j.molliq.2006.12.013]

25. Cubreli, G.; Psota, P.; Dančová, P.; Lédl, V.; Vít, T. Digital Holographic Interferometry for the Measurement of Symmetrical Temperature Fields in Liquids. Photonics; 2021; 8, 200. [DOI: https://dx.doi.org/10.3390/photonics8060200]

26. Psota, P.; Cubreli, G.; Kredba, J.; Stašík, M.; Ledl, V. Two Wavelength Digital Holographic Interferometry for Investigation of Dynamic Processes in Fluid Mechanics. Available online: https://dl.astfe.org/conferences/tfec2021,586fceb86446aa50,35eb4e3178980f3c.html (accessed on 1 August 2021).

27. Agarwal, S.; Kumar, V.; Shakher, C. Temperature measurement of wick stabilized micro diffusion flame under the influence of magnetic field using digital holographic interferometry. Opt. Lasers Eng.; 2018; 102, pp. 161-169. [DOI: https://dx.doi.org/10.1016/j.optlaseng.2017.10.019]

28. Li, C.; Liu, W.; Peng, X.; Shao, L.; Feng, S. Measurement of mass diffusion coefficients of O2 in aviation fuel through digital holographic interferometry. Chin. J. Aeronaut.; 2019; 32, pp. 1184-1189. [DOI: https://dx.doi.org/10.1016/j.cja.2019.01.012]

29. Xi, T.; Di, J.; Li, Y.; Dai, S.; Ma, C.; Zhao, J. Measurement of ultrafast combustion process of premixed ethylene/oxygen flames in narrow channel with digital holographic interferometry. Opt. Express; 2018; 26, pp. 28497-28504. [DOI: https://dx.doi.org/10.1364/OE.26.028497]

30. Meiners-Hagen, K.; Schödel, R.; Pollinger, F.; Abou-Zeid, A. Multi-Wavelength Interferometry for Length Measurements Using Diode Lasers. Meas. Sci. Rev.; 2009; 9, pp. 16-26. [DOI: https://dx.doi.org/10.2478/v10048-009-0001-y]

31. Malacara, D. Optical Shop Testing; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2007.

32. Guerrero-Mendez, C.; Anaya, T.S.; Araiza-Esquivel, M.; Balderas-Navarro, R.E.; Aranda-Espinoza, S.; López-Martínez, A.; Olvera, C.A.O. Real-time measurement of the average temperature profiles in liquid cooling using digital holographic interferometry. Opt. Eng.; 2016; 55, 121730. [DOI: https://dx.doi.org/10.1117/1.OE.55.12.121730]

33. Anand, A.; Chhaniwal, V.K.; Narayanamurthy, C.S. Diffusivity studies of transparent liquid solutions by use of digital holographic interferometry. Appl. Opt.; 2006; 45, pp. 904-909. [DOI: https://dx.doi.org/10.1364/AO.45.000904]

34. Naylor, D. Recent developments in the measurement of convective heat transfer rates by laser interferometry. Int. J. Heat Fluid Flow; 2003; 24, pp. 345-355. [DOI: https://dx.doi.org/10.1016/S0142-727X(03)00021-3]

35. Narayan, S.; Singh, A.K.; Srivastava, A. Interferometric study of natural convection heat transfer phenomena around array of heated cylinders. Int. J. Heat Mass Transf.; 2017; 109, pp. 278-292. [DOI: https://dx.doi.org/10.1016/j.ijheatmasstransfer.2017.01.106]

36. Herman, C.; Kang, E. Experimental visualization of temperature fields and study of heat transfer enhancement in oscillatory flow in a grooved channel. Heat Mass Transf.; 2001; 37, pp. 87-99. [DOI: https://dx.doi.org/10.1007/s002310000101]

37. Picart, P. New Techniques in Digital Holography; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2015.

38. Psota, P.; Tang, H.; Pooladvand, K.; Furlong, C.; Rosowski, J.J.; Cheng, J.T.; Ledl, V. Multiple angle digital holography for the shape measurement of the unpainted tympanic membrane. Opt. Express; 2020; 28, 24614. [DOI: https://dx.doi.org/10.1364/OE.398919] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32907000]

39. Bommareddi, R.R. Applications of Optical Interferometer Techniques for Precision Measurements of Changes in Temperature, Growth and Refractive Index of Materials. Technologies; 2014; 2, pp. 54-75. [DOI: https://dx.doi.org/10.3390/technologies2020054]

40. Landsberg, G.S. Оптика/Optika; 5th ed. Наука (Nauka): Moscow, Russia, 1976.

41. Psota, P.; Mokrý, P.; Lédl, V.; Stašík, M.; Matoušek, O.; Kredba, J. Absolute and pixel-wise measurements of vibration amplitudes using time-averaged digital holography. Opt. Lasers Eng.; 2019; 121, pp. 236-245. [DOI: https://dx.doi.org/10.1016/j.optlaseng.2019.04.002]

42. Blanche, P.A. Optical Holography: Materials, Theory and Applications; Elsevier: St. Louis, MO, USA, 2020.

43. Šidlof, P.; Riss, Š.; Vlček, V. Evaluation of Interferograms of unsteady subsonic airflow past a fluttering airfoil. Topical Problems of Fluid Mechanics 2016; Czech Technical University in Prague-Central Library: Prague, Czech Republic, 2016; pp. 223-228.

44. Ristic, S.; Linić, S.L.; Samardzic, M. Turbulence investigation in the VTI’s experimental aerodynamics laboratory. Therm. Sci.; 2017; 21, pp. 629-647. [DOI: https://dx.doi.org/10.2298/TSCI160130187R]

45. Kozic, M.; Ristic, S. Capability of two-dimensional Reynolds-averaged Navier—Stokes simulations for two-dimensional thrust vectoring nozzles. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng.; 2010; 224, pp. 905-910. [DOI: https://dx.doi.org/10.1243/09544100JAERO676]

46. Htun, Y.E.; Myint, Z.Y. Some Principles of Flow Visualization Techniques in Wind Tunnels. Int. J. Adv. Sci. Eng. Technol.; 2016; 4, pp. 2321-9009. Available online: http://iraj.in (accessed on 19 October 2021).

47. Kleine, H.; Grönig, H.; Takayama, K. Simultaneous Shadow, Schlieren and Interferometric Visualization of Compressible Flows. Opt. Lasers Eng.; 2005; 44, pp. 170-189. [DOI: https://dx.doi.org/10.1016/j.optlaseng.2005.04.009]

48. Emmert, T.; Lafon, P.; Bailly, C. Numerical study of self-induced transonic flow oscillations behind a sudden duct enlargement. Phys. Fluids; 2009; 21, 106105. [DOI: https://dx.doi.org/10.1063/1.3247158]

49. Ristić, S. Optical Methods in Wind Tunnel Flow Visualization; FME Trans: Belgrade, Serbia, 2006; Volume 34, pp. 7-13. Available online: https://www.researchgate.net/publication/277183175_Optical_methods_in_wind_tunnel_flow_visualization (accessed on 8 November 2021).

50. Ledl, V.; Psota, P.; Doleček, R. Digital holographic setups for phase object measurements in micro and macro scale. Proceedings of the EPJ Web of Conferences; Dresden, Germany, 25–29 August 2015; EDP Sciences: Les Ulis, France, 2015; Volume 92, 01001.

51. Doleček, R.; Psota, P.; Lédl, V.; Vít, T.; Václavík, J.; Kopecký, V. General temperature field measurement by digital holography. Appl. Opt.; 2012; 52, pp. A319-A325. [DOI: https://dx.doi.org/10.1364/AO.52.00A319]

52. Malacara, Z.; Servín, M. Interferogram Analysis for Optical Testing; CRC Press, Taylor & Francis Group: Boca Raton, FL, USA, 2005.