4.1. Effect of the Ratio of the Pitot Tube Diameter to the Nozzle Diameter
As an example,
Figure 2 shows a typical distribution of pressure along the jet axis. The data was taken from a supersonic (real) microjet escaping from the nozzle 16.1 µm in diameter by the Pitot tube 12 µm in diameter. The pressure curve has a quasi-periodic form associated with the shock wave structure of the jet. The peaks of the axial distribution of pressure were located at the ends of the jet barrels. The barrel length is one of the basic characteristics of supersonic jets. In turbulent jets, the barrel size and their transverse sizes from the first to the last barrel decrease owing to the evolution of the mixing layer of the jet. The parameters usually considered by most researchers are the length of the first barrel, the mean size of the barrels (mean sizes of the second, third, and fourth barrels), position and size of the Mach disk, and supersonic core length.
The measurements of the first barrel length and the mean size of the barrels in real microjets performed by the Pitot tube 12 µm in diameter [
7] showed that the mean barrel size was in excellent agreement with the first barrel size in macrojets [
11] (except for those JPR values for which the mixing layer produces a significant effect). However, the first barrel size turned out to be greater than the first barrel size in macrojets almost for all jets, as shown in
Figure 3, for all JPR values.
To understand the reason for systematic overprediction of the first barrel size in microjets determined on the basis of the Pitot tube measurements, the following problem was solved. Exhaustion of a real supersonic microjet from the 16.1 µm nozzle into the ambient space (
Pa = 1 atm,
P0 = 3.86 atm) at
n = 2.04 was considered. The Pitot tube was located on the axis of the jet at different distances
x/
d. The flow around the Pitot tube was calculated and the pressure inside the Pitot tube was determined. Additionally, the
P0′ in the free jet was calculated. The results of numerical simulation and experimental data are compared in
Figure 4.
It is seen that the calculated value of the pressure in the Pitot microtube located in the jet flow correlates fairly well with the experimental data. However, both these dependences were displaced with respect to the
P0′ curve in the free jet approximately by 0.5
d. This shifting was induced by the detached shock wave formed on the Pitot microtube. As the Pitot tube measures the pressure behind the shock wave, we actually have a situation where
x/
d in
Figure 4 corresponds to the Pitot tube position and
P0′ corresponds to the shock wave position ahead of the Pitot tube. The distance between the Pitot tube and the shock wave ahead of the tube was not constant over the supersonic core of the jet; it varies periodically in proportion to the Mach number at the jet axis. The shift of the experimental data with respect to the
P0′ curve for the free jet should be also affected by the ratio of the nozzle and Pitot tube diameters.
Numerical calculations were performed for free jets escaping from micronozzles 50, 75, 100, and 150 µm in diameter at
n = 1.8 and also for the same jets impinging onto the Pitot tube 12 µm in diameter. The results are plotted in
Figure 5.
The distributions of
P0′ calculated for jets having an identical JPR value but escaping from nozzles of different diameters almost coincide in the dimensionless coordinates. This distribution is shown by the solid curve in
Figure 5. As the ratio
d/
D increases, the effect of the Pitot tube diameter on the accuracy of determining the position of the first barrel of the jet becomes less pronounced.
Similar results were obtained in model microjet calculations. In particular,
Figure 6 shows the results of calculating the jet escaping from the nozzle 1.06 mm in diameter, modeling jet exhaustion from the nozzle 16.1 µm in diameter at
n = 2.
Figure 6 shows the calculated results for the free jet (solid line) and the calculated pressure in the Pitot tube located in the jet. The two sizes of the Pitot tube diameters were used in calculations—1.1 mm/0.8 mm and 0.2 mm/0.085 mm (outer/internal diameter). It is seen that the results were significantly affected by the shock wave position.
Figure 7 shows the
P0′ distributions along the axis of the jet escaping from the nozzle with the diameter 1.06 mm and the jet pressure ratio
n = 2. The pressure in the chamber was chosen in a way to model the flow from the nozzle 16.1 µm in diameter. The measurements were performed by the Pitot tubes with different diameters. The plots in
Figure 7b are fragments from
Figure 7a. The calculated data are shown by the solid black curve. The point of intersection of the curves with the dotted line corresponding to the pressure of the transition from the supersonic to subsonic flow shows the supersonic core length.
It is seen that the calculations predict a greater supersonic core length than that obtained in the experiments, as shown in
Figure 7a. This difference is explained by the presence of a certain level of turbulence in the mixing layer of the real jet, as compared to the completely laminar model implied in the computations. The supersonic core lengths calculated from the experimental data obtained by the Pitot tubes with diameters up to 1.1/0.8 mm coincide with each other. For large Pitot tubes (2/1.25–3/2 mm), the supersonic core length was 20% smaller. The possibility of reliable determination of the supersonic core length by the Pitot tubes with comparatively large diameters was explained by the wide profile of the transverse pressure distribution in the jet cross section where the supersonic core length was determined. In this cross section, the maximum of the pressure distribution was fairly wide, and the tube diameter did not induce significant errors.
The pressure distribution near the nozzle exit was characterized by shifting of the experimental curves with respect to each other as the Pitot tube diameter increased, as shown in
Figure 7b. As was demonstrated earlier, this shifting was determined by the detached shock wave. The greater the Pitot tube diameter, the greater the stand-off distance of the shock wave.
Though exhaustion of model microjets as a whole occurs in the continuum regime, there are local regions of reduced density in the jet, which correspond to the transitional flow regime from the viewpoint of the Knudsen number. In particular, the minimum density was observed in the middle of the first barrel, and it was in this region that the maximum influence of rarefaction on the results of Pitot tube measurements can be expected. It is worth mentioning that there was obvious disagreement in
Figure 7b between the minimum of the pressure distribution in the first barrel of the jet obtained by the smallest Pitot tube (0.4/0.17 mm) and the calculated data.
At the point of the minimum pressure
P0′ in the first barrel of the jet simulating exhaustion from the micronozzle with the diameter
d = 16.1 µm, as shown in
Figure 7b, the calculated gas density was
ρ = 0.013 kg/m
3, which corresponds to the molecular concentration
n = 2.7 × 10
17 molecules/cm
3. Then the mean free path of molecules at this point was
= 2.5 × 10
−3 cm. Here
σ is the molecule collision cross section (
σ ≈ 10
−15 cm
2). Thus, the local Knudsen number Kn directly calculated on the basis of the mean free path and Pitot tube size (
D = 0.4 mm or
Dint = 0.17 mm) is Kn = λ/
D 0.06 or 0.15, respectively. These values correspond to the transitional regime of the flow around the Pitot tube (from continuum to free-molecular). According to [
12], this is responsible for reduction of the Pitot tube readings, which was actually observed in experiments with these Pitot tube sizes. Moreover, because of gas rarefaction, the thickness of the shock wave ahead of the Pitot tube increases; as a consequence, the dependence
P0′(
x/
d) becomes less steep, which was again observed in experiments.
Despite the visible shift of the experimental data with respect to the true distribution of
P0′ at the jet axis, the mean size of the jet barrels was determined fairly accurately, which was confirmed by the experimental data in
Figure 3. From this viewpoint, the Pitot tube provides reliable results, even if the Pitot tube diameter was comparable with the nozzle diameter.
Figure 8a shows the normalized length of the supersonic core of the jet as a function of the ratio of the outer diameter of the Pitot tube to the nozzle diameter
D/
d for several values of
n. The supersonic core length was normalized to its value measured by the Pitot tube with the diameter
D = 0.4 mm. It was seen that the normalized length of the supersonic core decreases as the ratio
D/
d increases. If the measured data shown in
Figure 8a are plotted in the coordinates
, then all points fall on one decreasing curve, as is shown in
Figure 8b.
In the region of the jet flow transition from the supersonic to subsonic state, the velocity profile in the jet becomes significantly expanded in the transverse direction, and the Pitot tube readings become less sensitive to the Pitot tube diameter. As a result, the condition for the Pitot tube diameter with respect to the nozzle diameter with 3% error is defined as D/(dn) ≤ 0.5.
4.2. Effect of Jet Flow Unsteadiness on Pitot Tube Measurements
The measurements of the Pitot pressure distributions
P0′(
x/
d) on the jet axis at Reynolds numbers close to the conditions of the laminar-turbulent transition in the jet revealed significant steady fluctuations of the Pitot pressure. Such variations can be observed in both real, as shown in
Figure 9, and model, as shown in
Figure 10, microjets. In both real and model microjets, this effect was manifested for nozzle diameters 16–26 µm in a moderate range of JPR values.
As an example,
Figure 9 shows the pressure distributions in real microjets escaping from the nozzles 21.4, as shown in
Figure 9a, and 16.1 µm, as shown in
Figure 9b, in diameter. All data in
Figure 9 refer to the laminar regime of jet exhaustion, and minor pressure fluctuations were observed only in a small range of JPR values.
Figure 10 illustrates the emergence of pressure fluctuations with increasing JPR for the model microjet escaping from the nozzle with
d = 1.06 mm (domain 1) and modeling jet exhaustion from the nozzle 16.1 µm in diameter. The measurements were performed by the Pitot tube with the diameter
D = 0.4 mm. Exhaustion of the model microjet corresponds to the laminar flow regime at
n = 2 and 2.2, to the transitional flow regime at
n = 2.35, and to the turbulent flow regime at
n = 3 (when the sharp drop-off of the pressure has occurred). Domain 2 in
Figure 10b shows the measurements in the quasi-turbulent region of the jet, where random overshoots of the pressure
P0′ were observed. It can be noted that the spatial period of pressure fluctuations at the jet axis in domain 1 coincides with the spatial period of pressure fluctuations on the cells of the wave structure of the jet near the nozzle.
The letters in
Figure 10b indicate the points on the jet axis where the pressure evolution with time was measured. The measured results are shown in
Figure 11. It is seen that the pressure at the measurement points in domain 1 were almost steady, whereas the pressure in domain 2 behaved randomly with time, which can be attributed to the transition of turbulent spots in the jet. Possibly, the frequency of passage of the turbulent spots (pressure jumps) was appreciably higher than that shown in the plot, but the pressure fluctuations in the experiments were averaged because of the finite characteristic time of pressure relaxation in the Pitot probe, and the sensor provided only an averaged pattern.
This behavior of pressure in domain 1 was typical for the emergence of acoustic feedback in a steady flow of supersonic underexpanded jets where screw instability modes B and C develop near the laminar-turbulent transition region [
13]. The frequency of these modes was close to the frequency of acoustic waves in the ambient space, which, in turn, was equal to the velocity of sound in the ambient space divided by the doubled length of the gas-dynamic cells of the wave structure of the jet. The evolution and enhancement of these disturbances is the reason for global instability of underexpanded jets.
In the plane of the normalized distance from the nozzle versus frequency,
Figure 12 shows the spectra of acoustic oscillations detected by a piezoelectric sensor mounted near the jet in the case of motion of the Pitot tube with the diameter 0.4 mm for three values of
n—1.75, 2.5, and 3. It is seen that there was one frequency of oscillations for
n = 1.75, as shown in
Figure 12a, and two frequencies for
n = 2.5 and
n = 3, as shown in
Figure 12b,c. The frequency and amplitude of acoustic oscillations depends on the Pitot tube position on the jet axis.
The lower and upper frequencies of oscillations coincide with the frequencies of modes B and C of screw perturbations of the underexpanded jet, respectively. The periodicity of the emergence and vanishing of acoustic oscillations during the motion of the Pitot tube along the jet axis coincides with the periodicity of passing of the Pitot tube through the gas-dynamic cells of the wave structure of the jet. As the Pitot tube moves along the jet axis,
Figure 12 in addition to variations of the amplitude of acoustic oscillations also shows periodic fluctuations of the frequency of these oscillations, which were noted in [
14]. In the range
D = 0.4–0.7 mm, a weak dependence of the amplitude of acoustic oscillations on the outer diameter of the Pitot tube was observed.
As an example,
Figure 13 shows the pressure
P0′ and acoustic root-mean-square fluctuations of pressure
p′ at the frequency of mode B in the ambient space as functions of the normalized distance from the nozzle. The pressure
P0′ was measured by the Pitot tube with the diameter 0.4 mm. It is seen that the minimum points of
P0′ coincide with the peaks of the acoustic oscillations of pressure in the ambient space of the jet.
The results obtained offer a simple physical explanation for the emergence of intense fluctuations of the pressure
P0′ in the flow region close to the laminar-turbulent transition in the underexpanded jet. Intense fluctuations of screw modes of instability occur and develop in this region of the jet, leading to high-frequency motions of the gas-dynamic cells of the wave structure of the jet in the radial and azimuthal directions with respect to the jet axis with a spatial period equal to the length of two cells of the wave structure of the jet. Such motions were described in many publications dealing with instability of supersonic underexpanded macrojets, e.g., [
13,
15,
16,
17]. Interaction of the Pitot tube tap with an unsteady jet generates acoustic oscillations in the ambient medium, which propagate toward the nozzle exit, thus, creating feedback between the jet and the ambient medium, and can induce either enhancement or suppression of global instability of the jet. Enhancement or suppression depends on the phase of acoustic oscillations at the nozzle exit. In turn, the phase of acoustic oscillations depends on the Pitot tube position in the gas-dynamic cells of the wave structure of the jet. The feedback produces the minimum effects when the Pitot tube was located in a position where the phase at the nozzle exit corresponds to the node of acoustic oscillations. If the phase at the nozzle exit corresponds to the peak of acoustic oscillations, the maximum feedback effect is observed, resulting in jet flow failure, acceleration of mixing with the ambient medium, and decrease in
P0′. These minimums and maximums were reached every time when the Pitot tube moved in one of the cells of the wave structure of the jet.
Thus, the motion of the Pitot tube along the jet axis generates periodic overshoots of acoustic oscillations, as shown in
Figure 12. These acoustic oscillations provoke global instability of the underexpanded jet, which leads to acceleration of jet flow mixing with the ambient gas and to a decrease in
P0′ with a spatial period equal to the length of the gas-dynamic cells of the wave structure of the jet. From this viewpoint, the amplitude of fluctuations of the pressure
P0′ was actually overestimated because of the presence of the Pitot tube in the unsteady jet. Most probably, the true amplitude of Pitot pressure fluctuations would be appreciably lower if there were no Pitot tube in the jet.