Hot-Wire Investigation of Turbulence Topology behind Blades at Different Shape Qualities

: The scope of this paper is to perform a detailed experimental investigation of the shape error effect on the turbulence evolution behind NACA 64-618 airfoil. This airfoil is 3D-printed with predeﬁned typical shape inaccuracies. A high-precision optical 3D scanner was used to assess the shape and surface quality of the manufactured models. The turbulent ﬂow was studied using hot-wire anemometry. The developed force balance device was provided to measure the aerodynamic characteristics of the airfoil. Experimental studies were carried out for three angles of attack, + 10 ◦ , 0 ◦ , − 10 ◦ , and different chord-based Reynolds numbers from 5.3 × 10 4 to 2.1 × 10 5 . The obtained results show that the blunt trailing edge and rough surface decline the aerodynamic performance of the blades. In addition, the experimental results revealed a strong sensitivity of the Taylor microscale Reynolds number to the type of shape inaccuracy, especially at Re ≈ 1.7 × 10 5 . We also discuss the evolution of the Reynolds stress components, the degree of ﬂow anisotropy, and the power spectrum distributions depending on the airfoil inaccuracies. V.U. and V.Y.; V.Y. and D.D.; V.Y.; V.Y. and D.D.; writing—original draft preparation, writing— review and V.Y. and V.U.; visualization, V.Y., D.D. and T.T.; V.U.


Introduction
The behavior of the flow around a bluff body is a nonlinear problem. One of the consequences of the nonlinearity is an unpredictable sensitivity on boundary conditions deviations, however small they may be. The scale of possible reactions to some small geometric perturbations ranges from almost zero effect to a complex change of the flow structure. That makes it difficult to predict the turbulence evolution. The situation could be characterized as one of passive low control, however unwanted. Even small variations of the boundary conditions-shape and roughness of the airfoil-could cause important changes in aerodynamic characteristics.
In recent years, a large number of papers were devoted to aerodynamic characteristics of an airfoil. The researchers focused on the spatial and spectral investigation of turbulence parameters. Some of them, based on experiments, were studying in detail the physical aspects of the turbulence formation behind the bladed rotor. Others explored the structure of flow around an isolated airfoil. However, nearly ideal shape of the airfoil is considered for the experimental research with sharp trailing edge and smooth surface, which is far from real application. The airfoil profiles are designed with "absolutely sharp" trailing edge very often, this means that the thickness of the trailing edge approximates zero. This condition could be relatively well approximated in the case of laboratory model, but definitely not in the case of a real machine. The reasons could be technology of manufacturing, strength design or erosion, to mention the most important. The effect of the resulting geometry The problem is that the absolutely sharp trailing edge could not be made in practice for technological reasons [30]. Moreover, the size of the airfoil-its chord length-is very different between a real machine and a physical model used during aerodynamic experiments in wind tunnels. The scale could be 1:10, however, the thickness of the trailing edge and roughness of the surface are more or less the same, given by the manufacturing technology. For example, the real thickness of the trailing edge could be about 0.5 mm, the same for both the model (chord 0.1 m) and the real application (chord 1 m). This means that the real chord is shortened by approximately the same factor about 5 mm for the model (representing 5% error) and the application (error only 0.5%), respectively. This modification definitely has a very different effect on the two cases' flow-fields around the airfoils. Therefore, the geometry similarity is not maintained in physical modeling very often.
Please note that a rounded trailing edge is not necessarily cost-effective. It is known that the airfoils with blunt trailing edge could show a very good performance [31]. However, the flow-field is completely different, especially in structure and size of the wake, see, e.g., ref. [32,33]. This influences the pressure drag generation directly. This effect will be addressed in our paper in detail for the specific airfoil. However, the conclusions could be considered more general, at least in a qualitative sense.
In general, the flow evolution around an airfoil strongly depends on its geometry. Even the slightest deviation from an original shape leads to different flow quality, which significantly changes the aerodynamic characteristics of blades. Usually, shape inaccuracies can occur during the production or operation of the streamlined body. In the first case, this is due to technological factors in setting up and organizing the production process. Whereas, in the second, it is associated with the corrosion process [34] or depositions on surface. That is especially an issue for offshore wind turbines, where the blades are in a salt fog environment.
As we can see, the general evolution of the flow around an airfoil is strongly dependent on its geometry. Even the slightest deviation from a given shape can lead to different phenomena. The study of this issue is especially important in the production of various streamlined bodies.
That is why, in this paper, we investigate the impact of the most common airfoil shape inaccuracies on the flow structure behind it. The organization of the remaining parts of the paper is as follows. Methodology of experimental investigations and parameters of the measuring devices are given in Section 2. That section describes the procedure of NACA 64-618 airfoil productions and its modifications for the experiment, as well. Then, the results obtained by hot-wire anemometry are presented in Section 3. In particular, we discuss the features of the wake characteristics, Reynolds stress components, and the degree of flow anisotropy according to airfoil shape error, Reynolds number and angle of attack. Moreover, for some measuring points, the power and dissipation spectral distributions are analyzed. Afterwards, the effect of airfoil shape inaccuracies on the structure of turbulence (Section 4) and aerodynamic forces (Section 5) is evaluated. Finally, Section 6 reports the conclusion.

Airfoil Configurations
The NACA 64-618 airfoil was chosen for experimental studies. The reason for this is that its geometry demonstrates high-lift, soft-stall characteristics even at low-Reynoldsnumber flows, in addition to its overall good performance. The original dimensions of the airfoil model are chord length c = 80 mm, span l · c −1 ≈ 1.25 and thickness b · c −1 ≈ 0.18.
The general view of the airfoil design is presented in Figure 1a. The angle of attack was changed by turning the airfoil relative to its central through-hole. Then, the set position was fixed by mounting holes of the lower support. Each mounting hole changed the value of the angle by 10 • . The airfoil was manufactured from polymerized lactic acid (PLA) by 3D-printing technology, using a PRUSA i3 printer. During printing, the nozzle diameter was 0.4 mm, the height of the printing layer was 0.15 mm, nozzle temperature and printing speed were 230 • C and 30 mm · s −1 , respectively. Thus, the standard 3D-printing technology was used.
The characteristic of surface roughness is of prime importance for manufacturing an airfoil. It depends on the technology which is used for this purpose. In the case of 3D printing, the surfaces of the formed components have relatively high roughness. Its magnitude strongly depends on the technological parameters of production. Taking into account the results of [35], we can say that the average roughness of the surface at our printing parameters was Ra ≈ 22 µm.
For experimental research, three modifications of the NACA 64-618 airfoil were produced (A, B, C). Each of them simulated certain shape inaccuracies. As shown in Figure 1b, A and B have approximately the same surface roughness of 22 µm. The roughness of C was reduced to ≈5 µm by P2500 sandpaper. In addition, the chord length of A was ≈3.5% shorter than the others. Whereas, the thickness trailing edge 100 · h · c −1 for A, B and C were 0.62%, 0.12% and 0.05%, respectively.
The 3D noncontact optical scanning system GOM ATOS Core was applied to assess the airfoil shape deviation. Measurement data from this system formed the foundation of three-dimensional quality control of the finished model. Based on them, we calculated the deviation of the actual size of the airfoil profile relative to the original STL file format. As a result, the maximum standard deviation of form for all types of profiles was approximately ±0.15 mm ( Figure 2). Note a rough error during printing of the product A, fortunately far from the measurement height.

Wind Tunnels and Auxiliary Equipment
The experimental investigation was conducted in a low-speed open-type wind tunnel at the department of Power System Engineering, University of West Bohemia (Figure 3a) [36]. Its maximum operating velocity is 90 m · s −1 . Additionally, in the entire range of available velocities, the turbulence intensity is less than 0.3%. This feature allows us to accurately estimate turbulence evolution according to the geometry inaccuracies of the airfoil. A radial fan and a 55 kW three-phase electric motor are applied to generate the flow. The stable flow is enabled by automated monitoring of the reference pressure of the generated airflow in the tunnel contraction. LabVIEW software was used for experiment control.
The width, height and length of the test section were 0.3 m, 0.2 m and 0.75 m, respectively. Its bottom and the sides are monoblocks made of extruded transparent polycarbonate.
Whereas, the cover of the test section consists of separate segments. That allows the traverse to move along the test section (relative to the x-axis) (Figure 3b). Meanwhile, the traverse rod moves the probe relative to the crosswise y and spanwise z directions (Figure 3c)

Experimental Setup
The aerodynamic center of the NACA 64-618 airfoil was located at ∆z · c −1 ≈ 1.9 and ∆x · c −1 ≈ 2.8. Where ∆z and ∆x is the airfoil position relative to the test section. The scheme of the experimental investigation is shown in Figure 4.
Measurements were carried out for three angles of attack α = +10 • , 0 • , −10 • at different chord-based Reynolds numbers from 5.3 × 10 4 to 2.1 × 10 5 . Accordingly, α = 0 • corresponds to the position when the airfoil chord has the streamwise direction. In this case, the blockage ratio is equal to 2, 3%, while for the positive and negative angles, it is about 3%. According to Pope and Harper [37], for general tests in the wind tunnel, a maximum ratio of model frontal-area to test-section cross-sectional area should not exceed 7.5%. At higher values of the blocking effect, a correction method should be used.
The measuring cross-section was located behind the trailing edge at x · c −1 ≈ 0.1, 0.2, 0.4, 1.2 at y · c −1 ≈ 1.8. The width of measuring cross-section ranged from z · c −1 ≈ −1.5 to 1.5. Note that preliminary investigations indicate that the intensity of turbulence at given coordinates and a whole range of experimental velocities is less than 0.2%. The step k z between the measuring points at the instability and stability flow region was 3 and 10 mm, respectively (the methodology of identification can be found in [38]).

Measuring Techniques
The constant temperature anemometer was used to study the turbulent structure behind the airfoil. The basic principle of this device is to compensate for the thermal state of the hot wire [39,40]. The Wheatstone Bridge is used to measure the variation of the wire resistance relative to its temperature. In the case of wire cooling, the operating system of the device instantly supplies more energy to the wire. This is necessary to maintain a constant wire temperature. Thus, the power provided to support the system corresponds to the cooling rate and thus to the value of flow velocity passing by the sensor.
Experimental studies were performed using two types of hot-wire probes. In most cases, we used a miniature wire probe type 55P15. For some measuring points, we also used a miniature X-wire probe 55P63. That allowed us to investigate in more detail the two-dimensional structure of the turbulent flow (u and υ components). The wire probes with tungsten platinum coating have a diameter of 5 µm and an active length of 1.25 mm, while the ends of wire are gilded. The hot-wire signals for all experiments were sampled at 80 kHz. The filtering of the received signal was applied as well. For the high-and low-pass filter, the frequency was 10 Hz, and 30 kHz, respectively. The typical duration of a single point measurement was 10 s.
The StreamLine system from the company Dantec was applied to calibrate the velocity sensors. This system includes units for generating the flow and supply of dried air. The Collis-Williams law was used to calibrate and estimate the instantaneous flow rate [41,42]. This method is well described in many papers, see, e.g., ref. [43]. The method uses dimensionless description of the sensor cooling law using Reynolds and Nusselt numbers. This approach allows for compensation of pressure and temperature deviations during the experiments.
The precision of measure velocity was determined using the classical methodologies [44,45]. The analysis performed shows that the general uncertainty of the conversion process of the probe electrical signal to the velocity sample is around ±0.5%. This value, in the case of X-wire, is approximately equal for measurements in the streamwise and spanwise directions. Please note that the relative error is related to the maximal velocity, which was in the rage 10-40 m·s −1 , depending on the Reynolds number value.
The velocity time series are subjected to statistical data treatment-see, e.g., ref. [46]. LabVIEW and Gnuplot have been used to postprocess the data and produce all the plots presented in the experimental parts of this paper. The Welch method was applied to estimate the power spectral density of a one-dimensional time series. In addition, the solve function was used to determine the calibration constants of the hot-wire probes.

Wake Topology
A series of experimental studies were performed using the hot-wire probe 55P15. This section presents the analysis of velocity (only the dominant streamwise component u is measured) depending on the airfoil shape inaccuracies (A, B, C), angles of attack, Reynolds number and measuring position x · c −1 .
The experimental program is divided into two measurement campaigns. In the first part, the turbulent flow was investigated close to the trailing edge at different Reynolds numbers. Accordingly, the experiments were conducted at x · c −1 ≈ 0.1 and Re ≈ 5.3 × 10 4 , 1.07 × 10 5 , 1.7 × 10 5 , 2.13 × 10 5 , what corresponds to 10, 20, 32, 40 m · s −1 . Whereas, the second campaign took place at constant Reynolds number 1.7 × 10 5 ≈ 32 m · s −1 and different measuring positions x · c −1 ≈ 0.1, 0.2, 0.4, 1.2. For better visualization, the obtained data were placed relative to the position of the trailing edge at α = 0 • and normalized by the free-stream velocity U.
As expected, the obtained results showed a significant effect of the Reynolds number and the angle of attack on the wake characteristics. In particular, as the Reynolds number grows, the depth of the wake increases. The wake width strongly depends on the angle of attack and is the largest at α = +10 • . As we can see from Figure 5, the largest and smallest wake velocity deficits are characteristic at the C and A modifications, respectively. Moreover, at zero angle of attack, the difference between them is almost 30%. It is caused by that modification A, in contrast to C, has a 4 times greater surface roughness, 1.9% shorter chord length, and 0.62% thicker trailing edge. Thus, according to the experimental data at α = 0 • , when the Reynolds number increases to ≈2.13 × 10 5 , the wake velocity deficit behind C modification at α = 0 • , −10 • and +10 • increases by 2.3, 2.4 and 0.8 times, respectively. Moreover, at Re ≈ 5.3 × 10 4 and the positive angle of attack, the obtained velocity distributions also indicate the presence of a double peak. This phenomenon occurs as a result of the vortices' interactions with pressure and suction sides of the airfoil [33,47]. Whereas, in other cases, the absence of a double peak means that the flow contains only a single dominant structure. Additionally, the experiments clearly show the displacement of the wake center depending on the angle of attack. This effect occurs due to the difference in flow velocity from different airfoil sides.
In addition, we investigated the effect of airfoil shape inaccuracies on the wake evolution at increasing distances from the trailing edge. As shown in Figure 6, the wake depth reduces with increasing distance, while the wake width grows as the wake gradually wastes. At the same time, the wake center at 0 • and +10 • shifts to the pressure side while at −10 • , it is the contrary. Meanwhile, the absent airfoil shape error leads to a growing velocity deficit. Accordingly, at each measuring cross-section behind the airfoil, the velocity curve for C has a deeper wake than for A and B modifications. In particular, this pattern is clearly visible at zero angle of attack.

Reynolds Stress Tensor and Degree of Anisotropy
The Reynolds stress tensor was applied to estimate the effect of airfoil inaccuracies on the flow anisotropy. Usually, the Reynolds stress tensor has six independent components that characterize the flow relative to the spatial directions x, z, and y. However, in our case, due to the application of the X-wire probe 55P63, which measures only streamwise and spanwise components, the Reynolds stress tensor was statistically 2D.
The analysis of flow characteristics took place at two measuring positions x · c −1 ≈ 0.2 and x · c −1 ≈ 1.2 at Re ≈ 1.7 × 10 5 . For graphical interpretation, the obtained data were normalized by the squared free-stream velocity (Figure 7).
The obtained data indicate that the magnitude and the distribution of the normal stress components strongly depend on the slightest deviations in the airfoil shape. For example, closer to the trailing edge and α = 0 • , the value of u 2 for C modification by 2.3 and 1.75 times greater than for A and B, respectively. This explained is due to the variation of the trailing edge and the surface roughness of the airfoil. Pressure side Suction side Pressure side Suction side With increasing distance to x · c −1 ≈ 1.2, the value of normal stress components decreases significantly. This is especially noticeable at α = +10 • , where there is a meaningful dissipation of vortex energy.
In addition, the comparative analysis of the u 2 and υ 2 distributions reflect the double peaks. This phenomenon is more detail described in papers [7,48,49]. The authors claimed that at a positive angle of attack the double peaks more noticeable. However, in our case, this phenomenon is more visible at the α = 0 • and α = −10 • . That is because, for our studies, we used the NACA 64-618 airfoil. This is not only asymmetric but also has some concavity on the pressure side. Thus, the vortices which formed on the pressure and suction sides can propagate in the flow together. Whereas, at α = +10 • , the vortices passing from the suction side are increased by vortices from the pressure side. This interaction leads to the formation of a dominant vortex from the suction side in the general flow.
In addition, we also analyzed the effect of airfoil shape inaccuracies on the Reynolds shear stress distributions ( Figure 8). Physically, this component is related to the diffusion mechanism for the transport of fluid momentum. Moreover, using the υ u distribution, it is possible to clearly define the region of flow instability where the turbulence structure is anisotropic.
Graphical interpretation of the obtained data shows that, in the perturbation area, a negative value is located directly on the pressure side. This indicates that the turbulence intensity at this location progresses and consequently increases frictional resistance in the mean flow. As a result, turbulence receives energy, while the average flow loses it.   The obtained data show that the largest magnitude of the shear component occurs at the positive angle of attack and closer to the trailing edge. Simultaneously, at α = 0 • and α = −10 • the magnitude of υ u is three times smaller. Note the influence of shape inaccuracies on the obtained distributions observed at α = 0 • and x · c −1 ≈ 0.2. At this angle, the maximum of the shear stress component is 0.5 times lower for A and B than for C modifications. Whereas, at α = −10 • any effect of inaccuracies almost disappears.
According to the experiments, with increasing distance from the trailing edge, the width of the instability region increases, and its depth declines. This pattern is clearly traced for the A modification. That is due to the presence of additional vortices formed by the blunt trailing edge of the airfoil.
As noted above, the Reynolds stress tensor characterizes the turbulence anisotropy [50,51]. Namely, when the flow fluctuations are statistically uniform in all directions, then the shear stress goes to zero and the normal stresses become equal. Otherwise, at u 2 = υ 2 and u 2 υ 2 = 0, the turbulent flow is anisotropic. To quantify the degree of anisotropy, the Formula (1) could be applied, which represents the ratio of velocity component fluctuations.
where u 2 and υ 2 represent time-averaging value of velocity fluctuations in streamwise and spanwise direction, respectively. Values of DA greater or smaller than 1 indicate the dominance of the spanwise or streamwise component, respectively. Graphical interpretation of the obtained data at different measuring sections depending on the angle of attack is shown in Figure 9.
As we can see, the width of the anisotropy region is much bigger than the instability region determined by Reynolds stresses components. Comparing these two distributions, we can clearly define the range that characterizes the boundaries of the shear layer (the blue area on Figure 9). In the figure, the instability of the flow is much lower than in the wake. Thus, at some points, even with the smallest difference between the Reynolds stress components, the flow is more anisotropic.
A similar analysis technique is presented in paper [8]. In the study, the authors carefully investigated the flow structure around the FX 63-137 airfoil model. They note that there is some "transit zone" between the free stream and the wake instability region. A feature of this zone is an atypical distribution of power spectral density. Namely, the inertial range has a double slope.
Additionally, as we can see from Figure 9, in the region of Reynolds instability, the DA distribution has a double peak with a local minimum M. This point corresponds to the minimum of the streamwise velocity and the inflection point for the shear stresses distribu-tion. Moreover, the greatest degree of anisotropy observes close to the free stream, where streamwise fluctuations are dominant u 2 · υ −2 ≈ 3.5. Whereas, near to the wake, on the contrary, the transverse fluctuation is higher than the longitudinal one u 2 · υ −2 ≈ 0.5. In addition, the areas of the shear layer with the maximum anisotropy values are wider, and their position strongly depends on the angle of attack. For example, at α = +10 • and α = −10 • , they are observed from the suction and pressure side, respectively. Additionally, at the negative angle on the pressure side, a significant number of small vortices are generated. These are clearly expressed in the longitudinal direction.
In general, we can conclude that minor manufacturing inaccuracies of the airfoil do not have a significant effect on flow anisotropy. The main factors of this phenomenon are the general geometry of the airfoil and the angle of its attack.

The Spectral Analysis
The feature of turbulent flow is the presence in its structure of different spatial scales. Power spectral density is widely used to estimate the distribution of energy flux intensity between components of these scales. Physically, it shows the pattern of distribution in the time domain of signal power over the different frequencies.
To spectrally estimate the structure of the turbulent flow, we chose two measurement points, M and N. Accordingly, the first of them is located clearly at the local minimum of the longitudinal velocity component, while the second one is placed on the beginning of the shear layer. The Reynolds number was Re ≈ 1.7 × 10 5 . It is clear that at different angles of attack, the coordinates of the measuring points M and N will be changed (see Figure 10).
For representing normalized Eulerian power spectral densities of streamwise velocity fluctuations, we used the Kolmogorov nondimensional coordinates (Figure 11). Accordingly, the obtained spectrum was normalized by E(k) η · ν 5 1/4 and the wave number by k · η, with η, and ν representing dissipation scale, dissipation rate and viscosity, respectively. The method of estimating these parameters is show in Section 4.
As it is known, according to Kolmogorov's law, the spectral density of turbulence energy decreases with increasing wavenumber. Moreover, in the inertial subrange, the curve has a characteristic slope of f −5/3 [52,53]. This pattern is clearly observed for point M, whereby, in longitudinal and transverse directions, the turbulent fluctuations are statistically uniform. Our previous investigations show that at M, the turbulent flow is isotropic. Whereas, the opposite is valid for N in the shear layer. Moreover, this point has a double slope f −3 after that of f −5/3 in the inertial subrange. This tendency is observed at various airfoil modifications and different angles incidents. Thus, based on the behavior of the energy spectrum, we assume that point N can be determined as the place where the wake instability ends and the shear layer begins. Additionally, Figure 11 shows that the width of the inertial subrange significantly depends on the angle of attack. Physically, this subrange shows the interval where turbulence kinetic energy is transferred from larger to smaller scales without loss. The largest and smallest width of the inertial subrange is observed at α = +10 • and α = 0 • , respectively.
In addition, the dissipation spectrum for the points M and N was evaluated ( Figure 12). The dissipation occurs at a high wavenumber, where the smallest vortices dissipate turbulent energy due to viscosity. The magnitude of this phenomenon can be estimated by Formula (2).  For non-dimensional coordinates, we used normalization in the form D(k)/η · and for wave number k · η. The data obtained show that at point M, the energy dissipation distribution has a classical form. Namely, the curve in the inertial subrange has a slope of f 1/3 . A similar result for homogeneous isotropic turbulence was obtained by Pope [54]. However, for point N at the beginning of the dissipation range, there is some convexity. After that, the curve has the slope of f 1/6 . This feature is detected only in the place where the shear layer begins. Moreover, the flow characteristic at point N varies greatly depending on the type of shape inaccuracies. In particular, at α = 0 • , the value of turbulent kinetic energy for modification C is much greater than A. This pattern is quite logical. In the first case, the airfoil has the best geometry and low surface roughness. That is why the flow receives less resistance. Whereas, for modification A, at turbulence, the flow undergoes a large energy dissipation. However, in any case, we can see that even a small change in the airfoil configuration changes the flow structure.

Turbulence Structure Evaluation
Physically, turbulence is an assemblage of discrete perturbations that interact with each other, moving with the flow. Typically, turbulence measurements by hot wire are performed at fixed points, consisting of time history records derived from a sequence of many such perturbations or vortices. The motion inside vortices that are relatively distant from each other (in space or time) usually becomes statistically independent [46,55]. Thus, the measured signal can be considered as random and described in statistical terms. The above-mentioned approach allows finding the time-averaged characteristics of the whole array of vortices. This section presents an analysis of the turbulent structure behind the NACA 64-618 airfoil depending on its shape inaccuracies.

Turbulence Estimation at the Point
To estimate the structure of turbulence (only in the streamwise direction), we chose the measurement point M. It was located by coordinates x · c −1 ≈ 0.2, y · c −1 ≈ 1.8 and z · c −1 from −0.075 to 0.16. For these calculation, z · c −1 depended on the angle of attack, according to Figure 10. The experimental Reynolds numbers ranged from 5.3 × 10 4 to 2.13 × 10 5 .
The correlation analysis of time provides a lot of useful information about the mean temporal and spatial scales of turbulent flows [56,57]. Based on the methodology [58], the auto-correlation function R uu can be estimated by where u is the streamwise velocity component and the overbar is the average value. Additionally, σ 2 u and ∆t are the variance of velocity components and the shift time of the measuring signal, respectively.
The obtained distributions of R uu (see Figure 13a) allow us to determine the average lifetime of the largest vortices. Accordingly, the cumulative integral (see Figure 13b) from zero to the point where the value of R uu is equal to zero characterizes the value of the integral time scale T E . This can be represented as follows: Graphical interpretation of the integral time scale displayed in Figure 14. The obtained data show the strong dependence of T E on the Reynolds number and the angle of attack. Accordingly, a quadruple increase in Re leads to a decrease of integral time scale by 10-14 times. Moreover, as the value of Re decreases, the effect of airfoil shape inaccuracy increases. This pattern is valid for zero and positive angles of attack, while at negative angle it has the opposite tendency.  Thus, at α = 0 • and Re ≈ 0.53 × 10 5 , the presence of a blunt trailing edge and a rough surface (B modification) reduces the value of T E by 40% compared to the ideal shape of the airfoil (C modification). Whereas, at α = +10 • and the same Reynolds number, this difference is almost 75%. Unexpectedly, the greatest effect of shape error manifested at a negative angle of attack at a high Reynolds number, where inaccuracies in the airfoil geometry reduce the value of T E almost twice.
It is clear that the integral time scale T E and integral length scale L E are interrelated. The last one indicates the average spatial length of the largest vortices and can be determined by the following classical relation: where σ u is the standard deviation of the streamwise velocity component. The energy dissipation rate is one of the most important parameters for assessing the influence of a streamlined body on turbulent flow. It can be estimated from the inertial range dynamics of the turbulence using the second-order longitudinal structure function S 2 uu (t), which, according to Kolmogorov [52] the S 2 uu (t) can be found using the following equation: where 2.1 is imperial coefficient [59]. Finally, integrating the plateau of structure function, we can easily determine the value of the dissipation rate (see Figure 13c). Distributions of dissipation rate depending on the shape inaccuracies in streamwise directions are shown in Figure 15. As expected, the data show a significant effect of the Reynolds number and the angle of attack on the magnitude of the flow dissipation. Meanwhile, the influence of the airfoil shape is also clearly observed. For example, at α = 0 • and Re ≈ 1.7 × 10 5 the smallest value ≈ 4.0 × 10 3 m 2 · s −3 is observed for C modifications, while for B and A, this value is 22% and 30% higher. A similar trend is also present at the positive angle of attack. Moreover, the growing Reynolds number increases the shape impact. Accordingly, at α = +10 • and Re ≈ 2.13 × 10 5 , for modification A the value of is 75% higher than for C. However, at α = −10 • , there is some difference. Namely, in the range of Re ≈ 1.07 × 10 5 -1.7 × 10 5 , the dissipation rate increases sharply fourfold. Then it begins to decrease slightly.
Once the dissipation rate is known, Taylor microscale λ, Kolmogorov microscales (length η and time τ η ) can be estimated from classical relations: Additionally, we evaluated the Reynolds number based on the Taylor microscale by: As we can see from Figure 16, the value of the Kolmogorov microscale strongly depends on the flow velocity. When the Reynolds number quadruples, the microscale is halved. That is quite logical because with increasing Re, the number of small vortices increases significantly. In addition, the effect of shape inaccuracies on the size of the smallest vortices is also clearly visible. That is especially notable at low values of the Reynolds number and positive angle of attack. In these cases, at Re ≈ 0.53 × 10 5 , the difference of η between A and B cases is almost 22%, while at Re ≈ 2.13 × 10 5 , it decreases down to 9%.
It should be noted that the Kolmogorov microscale characterizes the size of the smallest vortices which are the source of energy dissipation. This pattern can be observed in the comparative analysis of Figures 15 and 16. Accordingly, as η decreases, the value of increases. Experimental data also revealed a rather interesting and unexpected distribution of Taylor Reynolds numbers (see Figure 17). Namely, regardless of the angle of attack, at Re ≈ 1.7 × 10 5 , there is a significant effect of shape errors on the R λ , while at other velocities, this phenomenon is much smaller. For example, at α = 0 • , the presence of a blunt trailing edge and surface roughness in A modification, compared with C, reduces the magnitude of the R λ around three times. The effect of shape inaccuracies of modification A is especially evident at α = +10 • . In this case, the magnitude of Taylor-Reynolds number for B and C modifications is the same and equal around R λ ≈ 470, while for A, it is three times smaller.

Characteristics of the Wake Turbulence
The purpose of this subsection is to evaluate the shape effect on the wake turbulence behind airfoil. The methodology of determination of the range of flow instability region is shown in Section 3.1. To estimate the turbulent structure in the streamwise and spanwise directions, we used the methodology described in Section 4.1. Experimental studies were performed at Re ≈ 1.7 × 10 5 and α = 0 • . The measuring cross-section was located perpendicular to the flow along the z-axis from −0.15 to 0.25 at x · c −1 ≈ 0.2. Figures 18 and 19 shows the graphical interpretations of turbulence characteristics depending on airfoil shape errors in streamwise and spanwise directions, respectively. According to the results, the greatest shape impact is clearly observed in the wake region. Namely, the values of the integral length scale and Taylor microscale undergo significant changes. This pattern is especially noticeable in the streamwise direction ( Figure 18). Thus, the maximum values of the integral length scale are observed at smooth surface and a sharp trailing edge C − L ≈ 0.8 mm. Whereas, with increasing roughness and the presence of a blunt trailing edge, this value decreases by half. In particular, this occurs for B − L E ≈ 0.4 mm and for A − L E ≈ 0.32 mm. Interestingly, the pairwise interaction of these inaccuracies and only the increase in roughness reduces the value of L E almost equally. This means that surface roughness is the most affected factor which changes the flow structure behind airfoil. Similar trends are present for the Taylor microscale distributions. Accordingly, in the case of B and A, the value of λ is the lowest at 0.78 mm and 0.75 mm, respectively. Whereas, for C, this value is almost twice as high.
In addition, the obtained distributions show the logical regularity of the vortex scales location relative to the sides of the airfoil. Namely, the integral length scale L E , in contrast to the Kolmogorov scale η, is located on the suction side. This is clear because from the pressure side occurs the formation of a significant number of small vortices. That leads to a significant increase in dissipation energy . Figure 19 shows that in the spanwise direction, the impact of shape errors on turbulence is negligible. However, there is some difference in the distribution of the integral length scale. Namely, this is the presence of a double peak. That indicates dominant vortices, which are formed on both airfoil sides and together propagate along the wake. Moreover, depending on the type of shape errors, their size and location vary. For example, in the presence of a blunt trailing edge and a rough surface, the maximum size of the vortices is observed from the pressure side. Whereas, on the suction side, the dominance vortices can be formed only in the absence of shape inaccuracy.

Lift Coefficient Determination
The aerodynamic balance was applied to estimate the effect of shape inaccuracy on the airfoil lift and drag forces [60]. Figure 20 presents a sketch of the developed aerodynamic balance.  The main structural parts of the balance are a resiliently suspended movable table and an eddy current sensor. Meanwhile, the elastic suspension consists of four parallel flexible elements, which have the lowest stiffness in the lift direction of lifting force. The displacement of the table is measured by an eddy current sensor which responds to the distance to the steel target. Unwanted vibrations are dampened by eddy currents in the copper sheet as it moves near to the neodymium super magnet. The sensor output is calibrated by a known force created by using a calibration weight in the lift or drag direction. It should be noted that the measurement error of this device was approximately 1.2%. This error was caused by unwanted airfoil oscillations and to shifting of "zero position" due to low stiffness of the construction parts made of plastic. The uncertainty of the displacement sensor is expected to play a smaller role.
The principle of balance is as follows. The experimental model of the airfoil is fixed to the suspension platform using an aluminum tube and a screw avoiding rotation. The aerodynamic force causes the displacement of the platform, which is recorded using the eddy current sensor. Then, the obtained values of the electric signal using the calibration curves are easily converted into a magnitude of the acting force. Figure 21 presents a graphical interpretation of the lift and drag forces depending on the shape inaccuracy of airfoil and surface roughness. As expected, the results obtained show that lift characteristics of the airfoil are higher for the best shape configuration. It is especially noticeable at the high Reynolds numbers. For example, at Re ≈ 1.6 × 10 5 and α = 0 • the value of Cl for C is 7 times higher than for A and B modifications with rough surface and blunt trailing edge. In addition, the measurements revealed that regardless of the type of airfoil shape error at α = −10 • , the values of the lift coefficients are negative. This phenomenon is not new and is typical at low angles of attack [61]. According to the authors of [62], the shape of the separated region, along with the differential growth of boundary layer displacement thickness on the upper and lower surfaces together, introduces an effective negative camber, resulting in a negative circulation around the aerofoil, thereby generating a negative lift.
Additionally, the experimental data suggest that the surface roughness contributes to increasing the drag force compared to the airfoil with a smooth surface. The main cause of this phenomenon is the separation mechanism, which is significantly affected by surface roughness that delays the separation and reduces lift characteristics [63,64].

Conclusions
The effect of shape inaccuracies on the flow structure behind NACA 64-618 airfoil was experimentally investigated. According to the data obtained, the roughness has a significant effect on the characteristics of the flow. Thus, at a zero angle of attack, the velocity deficit is 30% higher for a smooth airfoil surface. Generally, the effect of shape error decreases with increasing distance from the trailing edge.
The experimental study also found a large effect of the airfoil shape inaccuracies on the Reynolds stress components. Namely, at zero angle of attack, rough surface and a blunt trailing edge reduce the values of the normal and shear stress components by almost half. However, the analysis of the degree of anisotropy did not detect this effect. In this case, the main impact factors were the general geometry of the airfoil and its angle of attack.
Conducted spectral analysis at the original location of the shear layer detected a double slope of the distribution f −3 after that of f −5/3 . Moreover, at this point, the dissipation spectrum is also characterized by specific behavior displayed by convexity in the inertial range. The obtained distributions indicate that the contribution of wavenumbers to the turbulent kinetic energy is higher for airfoil without shape errors.
The analysis of turbulent flow structure revealed that the blunt trailing edge and surface roughness increase the value of energy dissipation rate and reduce the Kolmogorov microscale. In particular, the effect of inaccuracies is observed for Taylor microscale Reynolds number at Re ≈ 1.7 × 10 5 , where it increases almost three times.
Finally, the aerodynamic characteristics of the NACA 64-618 airfoil were estimated. According to the data obtained by aerodynamic balance, the smallest deviation of the airfoil shape leads to an increase in drag and a decrease in lift forces.
The Reynolds and Mach numbers are relatively small in our study, much lower than in any application in a real turbine or a compressor. In future, we intend to extend our study to working with Reynolds numbers at least an order of magnitude higher. The compressibility effect will be considered, as well.