1. Introduction
Aircraft icing is an important safety hazard which has been studied since the 1920s [
1] as it causes rapid performance degradation and malfunction of the flight instruments [
2,
3,
4]. The first icing research analyzed the effects of ice accretion over the wings [
5], and the importance of ice accretion studies remains relevant today [
6,
7,
8]. Based on the ice shape, scientists and engineers evaluate the impact of icing on aerodynamic surfaces such as wings [
5,
9,
10,
11], propellers [
12,
13,
14], and engine inlets [
15,
16,
17]. These evaluations of the ice accretion impact are used in risk assessment for icing events, the development of new aviation regulations, and the research and development of ice mitigation techniques [
18,
19,
20]. The aerodynamic penalties induced by different ice accretion regimes are studied both experimentally, in conventional wind tunnel tests with artificial ice shapes, as well as numerically, by means of simulations [
21]. For both research methods, ice accretion measurements are critical for an accurate outcome. Over time, multiple ice accretion measurement methods were developed; all brought more knowledge about ice accretion physics, and all had their advantages and drawbacks.
The first ice accretion measurement method was the manual tracing of ice shapes, which was also the most common technique in the past [
22,
23,
24]. It is an inexpensive technique, in which only paper, a pencil, and a metal support are required, but the ice shape drawing process can be long and have a dampening effect on small ice surface features, and the manual digitalization post-processing is cumbersome.
Following this, complex 3D ice accretion measurement techniques were developed. First, molding and casting techniques were used to capture the ice accretion features and to produce artificial ice shapes from plasters and epoxy [
25,
26]. Ice profilometry was conducted on cast models. Most molding materials were temperature-dependent—thus, they affected the ice shape, producing damped or deformed ice molds. The casting formation process usually involved entraining air bubbles during the mixing and pouring phases, which generated brittle and deformed ice shapes [
25].
Further, photogrammetry techniques were used to obtain 3D ice profiles [
9,
27]. The ice was photographed from multiple angles, and the shape was reconstructed using specialized software. Although it has increased accuracy and reduced cost compared to previous techniques, for cases of transparent ice (glaze accretion), a paint or powder layer must be added in order to avoid light reflections and refractions that would produce anomalies in the ice shape. These additional layers dampen some of the ice surface features and significantly increase the testing time because they require time to be deposed and to dry completely in order to obtain high-fidelity ice shapes [
9].
Another technique successfully applied for measuring ice accretion profiles is the structured light technique [
28,
29,
30], which relies on actively projecting a known light pattern on an object and extracting the 3D surface shape from images captured from one or more points of view. Being based on image correlation and the computation of displacement vectors between points of the projected grid, its results are spatial averages within a surface region. Thus, it has reduced accuracy in regions with steep ice thickness gradients. For accurate results, a parametric study on set-up alignment, grid projection resolution, and post-processing parameters must be performed before testing—a complex and time-consuming process.
More advanced laser-scanner techniques were developed and employed in the latest years. The laser-scanners have a high resolution and high accuracy and can generate both 3D and 2D ice profiles [
31,
32]. However, as with the photogrammetry, a layer of opaque painting must be sprayed on the ice before measuring. In addition, their high price (which usually varies between 50,000 to 100,000 euros) makes the technique cost-prohibitive for many research groups.
Expensive ice accretion methods are suitable for well-established icing facilities which focus on commercial and manned aviation. However, the developing fields of icing studies for different applications such as the unmanned aviation sector and wind turbine energy production have difficulties accessing these types of techniques. The icing research in these fields is seeing rapid growth, but it is still far from being fully established. This makes several advanced ice accretion measurement techniques cost-prohibitive. Especially for unmanned aerial vehicles (UAVs), for which icing impacts their operational range [
33], this measurement capability is critical. As a result of the UAVs’ reduced size and flight at lower Reynolds numbers [
34,
35] compared to manned aviation, new ice accretion studies must be conducted for their specific flight conditions.
Therefore, in order to support the development of the UAV icing research field, inexpensive ice accretion measurement methods must be developed. Ideally, these new methods would employ equipment that is easily accessible to any research group, would be time-efficient (which would reduce the general testing costs), and could be used both inside an icing wind tunnel facility and in field missions. Given the nature of the testing environments, the most important features of the new techniques are to be both low-cost and easy-to-use. Especially for unmanned aviation, this would accelerate the development of ice protection systems and the validation of numerical simulation tools, accelerating the result of UAVs obtaining all-weather flying capabilities.
Looking at possible inexpensive solutions, measurement techniques based on image processing algorithms are used in many fields [
36,
37], from the research of physical phenomena [
38,
39,
40,
41], of environmental events [
42,
43,
44], or medical research [
45,
46,
47] to the assessment of industrial faults [
48,
49,
50]. These types of algorithms provide non-invasive measurement capabilities, efficient ways to process big amounts of data, are cost-efficient, offer improved and reliable detection accuracy, and significantly reduce human error.
Therefore, this work proposes two different measurement techniques for 2D ice shapes based on image-processing algorithms. These were developed in order to offer more efficient and automatized alternatives to the time-consuming digitalization of hand drawings and low-cost alternatives to the 2D profiles obtained with the laser scanner technique. The algorithms are described in
Section 2. Their capabilities were tested on ice accretion measurements over UAV wings. The tests were conducted in the icing wind tunnel of the VTT Technical Research Center of Finland over a wide temperature range, a process which generated the main ice accretion types: rime, mixed, and glaze. The models were validated against caliper measurements for pointwise ice thickness measurements on the leading edge, and for the generated 2D ice shapes, a comparison was performed between the 2D profiles obtained with each technique and is presented in
Section 3. An additional algorithm was developed to convert the 2D ice profile into 1D ice thickness displayed over the chord length. This algorithm facilitates the use of thickness data at each airfoil coordinate and the comparison between the two image-processing techniques. Finally, each method’s capabilities are discussed together with their most suitable application and possible sources of error in
Section 4. Based on the techniques’ comparison, conclusions are drawn, and improvements are suggested.
3. Results
Three types of results are presented in this section, always comparing the data generated by the two 2D ice profile measurement techniques—the image processing (IP) and the scan detection (SD) methods. First, a quantitative result of the leading-edge ice thickness is used to validate the methods against caliper-measured experimental data. Secondly, a qualitative 2D ice profile comparison is conducted. As both techniques have their possible sources of error, this qualitative investigation cannot be considered a validation against an exact shape. However, a good agreement between the profiles of both techniques is considered a satisfactory result. Lastly, aiming for a simplified comparison of the two techniques, and between all the cases (regardless of the ice accretion regime, testing conditions, or airfoil shape) the 2D ice profiles are converted and compared in terms of 1D ice thickness profiles. Besides these, experimental observations about the ice accretion regimes are discussed—as the ice’s physical properties change (optical aspect, roughness-level of detail, shape, etc.).
3.1. Leading-Edge Ice Thickness Validation
The first validation step of the developed measurement techniques was the validation of the leading-edge ice thickness—which was considered the reference parameter for the conducted ice accretion study. This parameter was chosen as, together with the accretion regime, it can be used to describe the icing event severity, to estimate its impact on the aerodynamic performance of the wings, to assess the optimal de-icing configuration of the ice protection system, and to predict the ice removal time. Moreover, the leading-edge ice thickness could be accurately measured during experiments using an electronic caliper. Therefore, these measurements could be used as exact values to evaluate the results obtained with the two post-processing methods: the image processing and the scan detection ice measurement techniques.
Both methods measure the leading-edge ice thickness by considering the same reference point (experimentally marked) on the airfoil’s leading edge, L (x
LE, y
LE) = L (0,0), and interpolating the detected 2D ice profile to determine the ice thickness value, x
LE, at the y
LE = 0 coordinate. Thus, the values obtained with both methods can be directly compared and validated against the caliper measurements. Also, a percentual error estimation—Equation (1)—is used to assess the performance of the post-processing methods.
The numerical results employed in the validation of the image processing and scan detection methods against the experimental leading-edge ice thickness measurements are presented in
Table 3. Both average absolute and percentual errors are presented for each technique. The first is obtained by averaging the absolute value of the difference between the caliper-measured thickness and the result of each method, and the second by averaging the absolute percentual errors for each method.
The leading-edge ice thickness caliper measurement trends are in line with the literature [
9,
19]. The lowest ice thicknesses are observed for the highest IWT temperatures, T = −2 °C. The wetter airflow combined with the small temperature difference in the melting point generate a slow freezing regime in which the impinging water droplets do not freeze at impact with the airfoil, but instead travel downstream the profile and freeze later. This generates a transparent ice profile—glaze ice (
Figure 13, cases 183, 199)—as the air entrapped in the freezing process is minimal, with a rounded shape and the lowest leading-edge ice thickness due to the movement of water droplets and low freezing rate.
As the testing temperature decreases, the melting temperature difference increases together with the freezing fraction. Hence, the water droplets freeze closer to the impingement point, which generates more streamlined ice accretion shapes with a pointy aspect in the leading-edge region. For the same ice accretion time and AOA = 0°, the leading-edge ice thickness increases with a decrease in temperature. During this faster freezing process, the air entrapment percentage increases, and more opaque ice shapes are observed with decreasing IWT temperatures—from almost transparent in the case of mixed ice at T = −4 °C (
Figure 13, cases 14, 15, and 23) to milky white in the case of rime ice at T = −20 °C (
Figure 14, cases 245 and 247).
Figure 13.
Higher temperature range: T = −2, −4, and −5 °C, corresponding to glaze (183 and 199) and mixed (14, 15, 23, 110, 231, and 87) ice accretion regimes. The ice shapes are automatically detected in the experimental pictures by the image processing algorithm (detection drawn in black) and converted into 2D ice profiles (presented in
Figure 15,
Figure 16 and
Figure 17).
Figure 13.
Higher temperature range: T = −2, −4, and −5 °C, corresponding to glaze (183 and 199) and mixed (14, 15, 23, 110, 231, and 87) ice accretion regimes. The ice shapes are automatically detected in the experimental pictures by the image processing algorithm (detection drawn in black) and converted into 2D ice profiles (presented in
Figure 15,
Figure 16 and
Figure 17).
Figure 14.
Lower temperature range: T = −10 (runs 166, 12, 44) and −20 °C (runs 245, 247), corresponding to rime ice accretion. The ice shapes are automatically detected in the experimental pictures by the image processing algorithm (detection drawn in black) and converted into 2D ice profiles (presented in
Figure 18 and
Figure 19).
Figure 14.
Lower temperature range: T = −10 (runs 166, 12, 44) and −20 °C (runs 245, 247), corresponding to rime ice accretion. The ice shapes are automatically detected in the experimental pictures by the image processing algorithm (detection drawn in black) and converted into 2D ice profiles (presented in
Figure 18 and
Figure 19).
Concerning the numeric validation of the leading-edge ice thickness, both algorithms proved capable of capturing it. The results are in the same order of magnitude for all cases, regardless of the ice aspect, accretion time, or regime. By analyzing them based on the computed percentual error, the image processing method seems to capture more accurately the leading-edge ice thickness. Its maximal percentual error interval is from −4 to 6% of the error compared to the caliper data for all cases. For the scan detection method, the maximal error interval is one order of magnitude larger, from −23 to 24% over the entire testing range. As these are the maximal errors, it is beneficial to look at the absolute average percentual errors in order to better asses the techniques’ capabilities and to identify the cases in which algorithms perform best and the ones in which they lose accuracy. The average absolute percentual error in leading-edge ice thickness is ~2.2% for the image processing algorithm and ~9.5% for the scan detection method, and in terms of average absolute error values, these are ~0.26 mm for IP and ~0.98 mm for SD. This accuracy difference could be due to the increased number of steps in generating the scans for the second technique (manual drawing during IWT tests, post-test highlighting, and scanning) before the detection algorithm is applied, which could lead to an accumulation of errors. However, as both methods generate an absolute average percentual error lower than 10%, both techniques are considered validated for ice thickness measurements.
A more detailed error analysis allows for the identification of the drawbacks of each technique, an analysis which can be used for further improvement of the methods. The analysis of the first three highest error cases is presented below for each technique, while more general remarks and best practices with both methods are detailed in
Section 4.
Concerning the IP method, the cases generating the highest percentual errors compared to the experimental leading-edge ice thickness value are cases 245, 14, and 15, with corresponding errors of 5.8%, −3.6%, and 3.4%, respectively. The first case is a rime ice accretion profile, while the last two are mixed.
For case 245, the 5.8% underestimation of the leading-edge ice thickness is believed to be generated by a slight misalignment of the calibration paper with the airfoils’ leading edge when the picture was acquired. As the ice shape is opaque, the prediction step of the IP algorithm was used to identify the LE cross position based on 3 reference points from the calibration paper. The resulting reference point is visibly higher than the airfoil’s leading edge but is in line with the other cross signs from the calibration paper. For case 14, the 3.6% overestimation is most probably due to a camera loss of ice shape perpendicularity. In the picture, the spanwise thickness of the ice cut can be observed at the top of the leading-edge ice. This small thickness is detected as part of the 2D ice profile and slightly increases the measured LE ice thickness value.
For case 15, the 3.4% underestimation is also possibly due to a small camera misalignment and a loss of perpendicularity, which is visible from the inclination of the leading-edge line towards the left of the image instead of going straight down. This might have generated a small perspective error which could have marginally reduced the measured leading-edge ice thickness.
Regarding the SD method, the top three highest percentual errors are generated by cases 231, 110, and 12, with corresponding errors of 23.5%, −23.4%, and 18.3%, respectively. The first two tested cases are mixed ice accretion, while the latter is a rime accretion case.
Run 231 is the one with an angle of attack of 8°. Hence, the maximum thickness is not at the leading edge but is instead moved towards the pressure side of the airfoil. This makes it more difficult to measure the experimental leading-edge thickness with the caliper, which may cause larger uncertainties for the reference value. However, this increased uncertainty does not explain the 23.5% percentual error. An extended analysis of the possible sources of error for this case is presented in
Section 3.3.
The overestimated ice thickness of the scan detection algorithm for run 110 is likely related to general uncertainties and the thickness of the green pencil. Typically, each scanned shape is printed on DIN A4 before re-tracing. In this case, the thickness of the green pencil could easily be 0.5 mm scaled into the original size. This can introduce a large relative error for small ice shapes such as the one in run 110. Additionally, the SD code tries to find the center of the line for stability reasons, although the innermost point would be a better representative of the actual ice shape.
The maximum thickness for the traced ice shape from run 12 is not at the leading edge but is instead moved slightly towards the suction side. When taking the maximum thickness of the SD profile in close vicinity of the leading edge, the thickness matches the caliper measurement much better. This could indicate that the paper template was not in the correct position and was slightly rotated. A second explanation would be that the maximum thickness might have been rotated from the leading edge, and, similar to run 231, the caliper measurement has a higher uncertainty. Additionally, the original tracing has many lines close to the leading edge. Hence, re-tracing is a more difficult process for this ice shape and could lead to higher errors.
3.2. 2D Profiles Comparison
In order to facilitate comparison of the 2D ice profiles obtained with both codes, the profiles’ positions were corrected by rotating one of them to align the airfoils’ leading edge, aiming for a perfect alignment, especially on the wings’ suction side. The comparison results are presented in
Figure 15 for the glaze ice accretion regime, in
Figure 16 and
Figure 17 for the mixed one, and in
Figure 18 and
Figure 19 for rime ice.
Figure 15.
2D glaze ice accretion profiles (T = −2 °C)—comparison between the image processing (IP) and the scan detection (SD) techniques—for (a) run 183 and (b) run 199.
Figure 15.
2D glaze ice accretion profiles (T = −2 °C)—comparison between the image processing (IP) and the scan detection (SD) techniques—for (a) run 183 and (b) run 199.
Figure 16.
2D mixed ice accretion profiles (T = −4 °C)—comparison between the image processing (IP) and the scan detection (SD) techniques—for (a) run 14, (b) run 15, and (c) run 23.
Figure 16.
2D mixed ice accretion profiles (T = −4 °C)—comparison between the image processing (IP) and the scan detection (SD) techniques—for (a) run 14, (b) run 15, and (c) run 23.
Figure 17.
2D mixed ice accretion profiles (T = −5 °C)—comparison between the image processing (IP) and the scan detection (SD) techniques—for (a) run 110, (b) run 87, and (c) run 231.
Figure 17.
2D mixed ice accretion profiles (T = −5 °C)—comparison between the image processing (IP) and the scan detection (SD) techniques—for (a) run 110, (b) run 87, and (c) run 231.
Figure 18.
2D rime ice accretion profiles (T = −10 °C)—comparison between the image processing (IP) and the scan detection (SD) techniques—for (a) run 166, (b) run 12, and (c) run 44.
Figure 18.
2D rime ice accretion profiles (T = −10 °C)—comparison between the image processing (IP) and the scan detection (SD) techniques—for (a) run 166, (b) run 12, and (c) run 44.
Figure 19.
2D rime ice accretion profiles (T = −20 °C)—comparison between the image processing (IP) and the scan detection (SD) techniques—for (a) run 245 and (b) run 247.
Figure 19.
2D rime ice accretion profiles (T = −20 °C)—comparison between the image processing (IP) and the scan detection (SD) techniques—for (a) run 245 and (b) run 247.
The airfoil alignment between both techniques is achieved in all cases—finding a good agreement between the pre-cut metal airfoil shape used for the SD technique and the IP obtained profile. This proves that the image pre-processing corrections were well applied, as the IP algorithm generates a virtually undistorted airfoil profile.
Based on this observation, the obtained 2D ice shapes are also considered to be unaffected by distortion. In some cases (as in
Figure 15a,b and
Figure 16a,b and
Figure 19b), the IP-generated airfoil profile diverges from that of the SD downstream of the leading edge on the pressure side. This occurs because, in pictures, the pressure side appears blurred due to a difference in perspective, as the camera is aligned and focused on the suction side and leading-edge sides. For this reason, the airfoil profile is not clearly defined in the pressure side area, and errors can occur in the manual selection of the airfoil points. Greater airfoil profile abnormalities could occur if the wing span is not completely ice-free and ice is remaining along the optical path which is integrated in the blurred region, thereby generating an airfoil profile distortion (as with run 183 in
Figure 13 and
Figure 15a). However, these airfoil profile abnormalities do not affect the quality of the ice profile detected by the IP method, the two profiles being generated in different manners, as is detailed in
Section 2.
The qualitative comparison of the 2D ice shape results shows that the image processing and the scan detection techniques generate comparable profiles, which, in most cases, overlap and capture the same features of the ice profile. The absolute ice thickness values have some discrepancies, especially in downstream regions.
However, a significant difference between the results of the two methods is observed in terms of line continuity: the SD method generates a continuous line, while the IP ice shape profiles are sometimes interrupted, and they do not always connect to the airfoil at the ends. This occurs because the pre-processed SD profiles are hand-drawn, starting and ending at the airfoil’s surface, while the airfoil–ice interface is not always clearly visible in the images employed by the IP method. This is because the camera integrates the airfoil’s surface in the spanwise direction, and a blur effect occurs in the regions to which the camera is not perfectly perpendicular (visible in, e.g.,
Figure 14, run 166). The interrupted aspect of the IP profiles is due to the filtering step applied in the post-processing to correct for detections of shadow (e.g.,
Figure 14, run 247) or light reflections (e.g.,
Figure 13, run 15) which could influence the measurement if not excluded. Afterward, an analysis based on the testing temperature ranges was performed in order to assess the applicability of these methods to certain ice accretion regimes. The best agreement was observed for the rime ice shapes, which correspond to the lowest IWT temperature regimes (T = −10, −20 °C) and are depicted in
Figure 18 and
Figure 19. The streamlined ice shape is specific to rime accretion, with the highest ice thickness being present in the leading edge and decreasing downstream of the airfoil, and a relatively smooth surface is well captured by both methods. The profiles overlap in the leading-edge region and maintain a similar shape and thickness size, although the IP ice profiles continue a bit further downstream on the airfoil and present more roughness (in the form of small irregularities) than the SD ones. This could be explained by the fact that the SD profiles are limited by the size of the drawing paper and the metallic support, while the IP processes the field of view captured by the camera in the image. Also, the small roughness variation is mostly dampened in the hand-drawing process of the profile, which precedes the SD processing. Therefore, the SD profiles appear smoother. The IP profile may appear slightly ampler than that of the SD towards the airfoil’s tail, and the ice thickness in this region may be marginally increased in some cases. This is because, for non-transparent ice shapes (rime or mixed), the thickness of the 2D ice cut might be visible in the image towards the downstream airfoil area. The camera’s sensor starts to be misaligned and lose perpendicularity to the ice shape in the downstream area as it is focused on the airfoil’s leading edge. This does not always occur, and its effect is considered to be negligible, as discussed in
Section 2.2.
For the temperatures corresponding to mixed ice accretion (T = −4, −5 °C), the profiles captured by the two techniques show a good agreement on the general shape, generally matching both the leading-edge ice thickness and the downstream ice evolution. Both techniques present more ice shape irregularities and capture the increased roughness of this accretion regime (e.g., visible ice feathers in
Figure 13, run 14,15). However, these formations do not perfectly overlap between the two profiles. This suggests a possible camera misalignment error with loss of perpendicularity to the ice profile (most noticeable in
Figure 13, run 15, by observing the inclination of the airfoil’s leading-edge line towards the left side of the image instead of straight down as in run 23, where the camera is well aligned). This misalignment is visible in the results’ comparison (
Figure 16, run 15), which shows an ~1 mm vertical offset between the SD and the IP profiles. This vertical offset generates further misalignment of the ice irregularities in the downstream region of the two ice profiles. Methods to avoid these errors are recommended in
Section 4. Based on the leading-edge thickness validation and the agreement between the 2D ice shapes generated with the SD and IP methods, the errors introduced by this misalignment are considered negligible, and the ice profile is representative of the studied mixed ice accretion cases.
For the highest testing temperatures, (T = −2 °C), the glaze ice profiles are well captured by the two 2D measurement techniques. The results from both the SD and IP algorithms measure the leading-edge ice thickness with errors below 7.5% and 2%, respectively, compared to the experimental value. In addition, both methods capture the leading-edge horns well, together with the downstream ice protuberances (
Figure 13, run 183 and 199). Although both methods capture all the details of the profiles, there is less overlap than for the rime and mixed ice accretion regimes (
Figure 15a,b). This discrepancy could be due to the ice’s characteristics, which, in this case, forms in a wetter regime, has a more complex shape, and has a transparent appearance. The first two differences possibly increase the uncertainty during the manual drawing step of the SD technique, as the increased ice slippage and required level of detail make the drawing prone to more errors. Moreover, the drawing is made against the support paper, thereby acquiring a purely 2D profile for the SD method. For the IP method, the picture is prone to encompassing some 3D effects in this case, as the ice has an increased complexity, with its aspect varying along the span of the cut. Thus, combined with the ice’s transparency, the resulting 2D profile is more the ice shape integrated along the ice cuts’ thickness. To achieve high accuracy, the ice cut should be kept as thin as possible, but in this case, it is ~1 cm. Therefore, the slight discrepancy between the 2D profiles generated with the two methods is not concerning, especially as the general aspect is conserved. By keeping in mind how each technique acquires and processes the profiles, both can be used to generate representative 2D ice shapes for glaze ice accretion.
3.3. 1D Ice Thickness Results
The 1D thickness ice profile represents an algorithm that was developed to geometrically linearize the 2D ice shapes and to facilitate the direct comparison of the ice accretion results regardless of the airfoil shape. The algorithm can be applied to digital 2D ice profiles generated by any technique (SD, IP, or other), as long as an ice profile is provided with the corresponding airfoil points. Subsequently, it generates a graph with the ice thickness for each airfoil surface point.
The 1D thickness algorithm allows for fast identification of the ice accretion regime from profiles’ characteristics, such as leading-edge ice thickness, maximum ice thickness, ice coverage, profile aspect, etc., a fast investigation of the icing severity, and easier approximation of the accumulated ice mass for a given profile.
The rime ice accretion regime can be easily identified from the 1D profiles. The profiles present a very important ice thickness in the stagnation point (in this case, on the [0,0] point of the leading edge, as the airfoils were positioned at a 0°AOA)—
Figure 20 and
Figure 21. A detailed view of the stagnation point region can be found in the
Appendix C,
Figure A4. The rime ice accretion profiles present the highest leading-edge ice thickness between all studied cases as the impinging water droplets freeze at impact with the airfoil. From this point downstream, there is a steep decrease in ice thickness on both the negative and the positive sides (pressure side and suction side, respectively), followed by a relatively smooth and gradually decreasing ice profile, which represents the streamlined rime ice accretion aspect. For most cases, the ice spreads further downstream on the airfoils’ suction side (x positive), which is expected, due to the increased air speed along this part due to the cambered airfoil shape (RG-15). For run 247 (
Figure 21b), the ice is approximately evenly spread, as the tested airfoil (SD8020) has a symmetric shape.
The glaze ice 1D profile can be identified as having the opposite characteristics of the rime ice profile. There is a significant amount of ice in the stagnation point, but its thickness is comparable with the thickness of the ice horns—other ice protuberances formed close to the leading edge, which can be noted as spikes in
Figure 22a,b. These occur as, due to the high temperatures close to the melting point, the water droplets do not freeze upon impact, but instead travel downstream of the airfoil profile and freeze at different locations. This generates a more even 1D ice profile, where the central region has the highest ice thickness values, followed by significantly shorter ice protuberances downstream (in these cases, on the x-positive suction side).
The mixed ice regime generates 1D ice profiles with combined characteristics from the two previous ice accretion regimes. As with the rime ice profiles, it is characterized by the highest ice thickness present at the stagnation point (
Figure 23 and
Figure 24), followed by a less steep decrease in ice thickness downstream on the profile, with a high roughness region (
Figure 24) and ice spikes which represent the 2D ice protuberances and ice feathers (
Figure 23). Two atypical results are observed in
Figure 24a,c. These correspond to experiments performed to test slightly different conditions, such as the impact of the ice accretion time or of the airfoil’s AOA over the ice shape profile. At the same time, these runs also test the measurement techniques’ capabilities.
First, run 110 was tested in order to investigate the ice accretion severity and shapes for a reduced ice accretion time (8 min instead of 20). The amount of accreted ice is reduced because of the reduced ice accretion time. The thickness is still quite evenly distributed around the leading edge, with some visible roughness formations downstream and with an already-formed ice peak at the stagnation point. Despite the reduced accretion time, the ice accretion regime can still be easily identified from the ice thickness profile.
Second, run 231 was tested at the general ice accretion time, but with the airfoil positioned at an increased AOA (8° instead of 0°). This experimental configuration was tested in order to observe its impact on the ice shape profile, especially on the position of the maximum ice thickness. A detailed view of the
Figure 24c profile is presented in the
Appendix C, (
Figure A7), where the maximum ice thickness position can be easily identified at about –0.1 mm, thus 0.1 mm below the leading edge’s highlight point on the airfoils’ pressure side.
Moreover, for this particular run (231), the 1D thickness analysis was useful in identifying the other sources of error which motivate the 23.5% leading-edge ice thickness error for the SD method (
Table 3) and for the differences in 2D profiles between the two techniques (
Figure 17c). From a detailed view of the 1D thickness profiles’ comparison (
Figure A7c), it was observed that the profiles were highly similar in the leading-edge region while having an offset of ~1.2 mm. An enlarged detail of the leading-edge region reinforced the profile similarity (
Figure 25). A closer investigation of the raw data for the 2D ice shape measurement techniques (presented in detail in the
Appendix D,
Figure A8) suggests the misalignment of the support paper before the profile hand-drawing of this run. This seems to have misplaced the LE marker further in the positive x direction of the 2D profile, thereby losing the leading-edge contact point with the airfoil. This might have generated the ~2 mm-lower LE ice thickness value of the SD method. The misalignment was not noticed during the experiments. It is considered to be the most important source of error for the SD method in this case, with the influence of the increased caliper measurement uncertainty having only a secondary effect. Regarding the IP result for this run, the misplacement of the alignment paper did not represent an issue, as the optical path was clear and the airfoil’s leading edge was identifiable directly from the image. Moreover, its error of 1.3% is within the uncertainty interval of the IP measurement method.
4. Discussion
In this section, a comparison of the two techniques’ capabilities is presented. For each method, the accuracy, the applicability, the advantages and disadvantages, and the possible sources of error are analyzed. Advice for maximizing the methods’ efficiency and accuracy is summarized in a best practice guide, and possible improvements are suggested. The methods’ capabilities are again compared to the measurement techniques available in the literature, mostly in terms of cost and availability.
As observed in
Section 3, both the scan detection and image processing 2D ice shape measurement methods can generate meaningful results for all of the tested ice accretion regimes (rime, glaze, and mixed) in terms of leading-edge ice thickness and 2D ice profiles.
In terms of quantitative measurements, the image processing shows a lower percentual error against the caliper data in
Table 3, with leading-edge ice thickness differences ranging between −0.4 to 0.7 mm. The scan detection technique shows an increased error, the difference with the caliper measurements being from −1.2 to 2.6 mm. This scan detection error range is in line with that observed in Reference [
9] for the manual tracing technique. This is a consistent result, as the scan detection technique is based on the acquisition of the ice profiles via manual tracing. Therefore, the image processing method is preferred, as human error is greatly reduced. In terms of qualitative measurements of the 2D ice shapes, the scan detection method generates complete ice profiles which are closed around the airfoil. In most cases, these are of an approximate thickness and shape, with a smoother aspect and fewer features than the original ice profile. However, for applications in which an approximate shape is sufficient, for example CFD simulations of the airfoil’s performance degradation due to icing, the scan detection technique might be preferred, as it generates continuous-line profiles.
The 2D profiles generated with the image processing technique capture more accurately the ice thickness and the sharp and small features of the ice due to the automatic detection method from experimental images. This process minimizes human error and generates consistent results for each ice profile.
Due to these aspects, the image processing method might be used to generate reference ice shapes. However, its drawback is generally an interrupted profile line—with possible gaps—depending on the quality of the initial picture. In terms of how much time and effort are necessary for both techniques—from data acquisition to the final CFD-suitable 2D ice profiles—both techniques use four steps. The largest time difference occurs during the data acquisition step of both techniques. For one ice profile, it takes ~2 min to capture high-quality images using the IP technique, while an experienced experimenter would take ~5 min to draw the ice tracing. Thus, the SD technique is slower and is highly dependent on the experimenter’s experience.
For the IP method, the four steps are taking the ice profile images, image pre-processing correction, applying the post-processing, and the profile refinement step. Excepting data acquisition, all of the steps are software based and automatic—the entire process takes in the order of a few minutes to complete. To improve the speed and accuracy, one could pre-test in order to optimize the camera position for alignment and the reduction of optical imperfections and then use this predefined position in all tests. This would reduce the data acquisition and profile refinement times.
For the SD method, the four steps are manual ice tracing, redrawing of the profile in color, scanning of the profile, and post-processing. As the first three steps must be executed manually, the entire process takes about ten minutes to complete. There is not much possibility to accelerate the process, as the most time-consuming step is the manual tracing during the experiments. Acceleration of that step would most probably decrease the profile accuracy. The accuracy of the SD method could be improved by replacing the manual redrawing of the profile in color with using a graphic editor software. However, this would most likely significantly increase the pre-processing time.
As for their applicability, both techniques can be easily used in a laboratory setup. Due to the faster and less cumbersome data acquisition process, the IP method might be more suitable for field tests as long as good lighting conditions can be achieved. Another advantage of this method is that it is photo-documented, which allows for easier identification of alignment errors and gives insights into the ice’s physical properties (such as opacity, roughness, and surface ice formations). For the SD method, the only acquired data is the contour drawing of the ice shapes. However, even if it is more cumbersome, the SD method can be applied under any lighting conditions.
The two methods were developed with the aim of increasing the availability of 2D ice profile measurement techniques, generating methods that could be more automatic, efficient, accurate, and widely accessible to icing research groups. As necessary equipment, both employ a computer with a MATLAB license and a digital camera or a scanner. This equipment is generally available in most research facilities, and it is very cost-effective compared to other existing ice profiling techniques (laser scanners, molding, and casting methods, etc.).
Both methods have their specific error sources, which can be addressed with good experimental practices to decrease the measurement uncertainty. Below, a summary of the best practices for both techniques is presented as a guideline to maximize their accuracy.
First, for both methods, it is advisable to take a caliper measurement of the leading-edge ice thickness, taking care that the caliper is normal to the airfoil’s surface and is not placed at an angle, thereby avoiding reading errors. Although the caliper method has its uncertainty, it is useful to have an idea of the expected value along the leading edge. This enables the identification of alignment or post-processing errors and might suggest ways to correct them. When applying this to any of the techniques, careful labeling of the ice cuts is important for all the data acquisition (including the caliper measurements). Otherwise, it is nearly impossible to determine which profile corresponds to a certain spanwise position or measured leading-edge ice thickness. Additionally, it is recommended to take multiple, independent caliper measurements of the leading-edge ice thickness. This allows for a better evaluation of the uncertainty of the caliper measurement and a more accurate view of the performance of ice shape acquisition techniques. Analysis of accuracy of both methods should be improved by measuring the ice thickness with a caliper at multiple different locations on the airfoil. To use more sophisticated methods of comparing ice shapes [
54], the full ice profile would be required. In addition, a comparison of the two methods to results obtained with already established techniques would generate more knowledge on the level of accuracy.
Second, for both techniques, it is essential to have a good alignment of the calibration paper with the leading edge of the airfoil. After taking the ice cuts, it is necessary to check that the supporting background plate and the calibration paper touch the airfoil surface. Two possible errors were observed: either the ice cut was not conducted properly, and the support plate was placed on a thin layer of ice instead of the airfoil surface (usually for glaze ice cases), or the plate was placed in contact with the airfoil, but the calibration paper was not touching the airfoil (this occurs for opaque ice—rime/mixed—where the reference point is not visible through the ice). Both have the same effect and generate lower ice thickness values, as the reference point is placed higher than the leading edge of the airfoil. This error can be identified by comparing the leading-edge thickness results with the caliper measurements. For the IP method, it can also be observed by calculating the position of the reference point based on the network of points on the calibration paper.
The support plate–calibration paper ensemble must be kept in a normal direction to the airfoil’s surface during data acquisition. This practice will avoid projection errors in the ice profiles for both techniques.
The best practices for the image processing method focus mainly on the camera alignment and obtaining a good quality of the initial image. As previously stated, ideally, a trial of the camera position should be conducted in order to position the camera perpendicularly to the calibration plate. To check this alignment, it is important to capture multiple calibration signs in the field of view so as to compare the magnification factor between 2 consecutive points from different regions of the paper. The magnification factor should not show significant variation between regions.
In order to well-capture the suction side and pressure side of the airfoil, the camera should be positioned so as to have a clear view of both parts and avoid blurred regions, but without placing the camera too far from the airfoil. This would decrease the spatial resolution and increase the uncertainty of the ice edge detection. If this is not possible with one picture in some specific setups, two pictures, one for the suction side and one for the pressure side, could be taken and the resulting profiles merged during post-processing. During this process, it is essential to keep the same magnification factor in both images.
In addition to avoiding light reflections or shadows in the image, the quality of the 2D ice profile can be improved by photographing thin ice cuts (as thin as possible below 1 cm). This will limit the capturing of 3D effects in the 2D ice profile. Another way of improving the ice edge detection quality in the images would be the use of colored background paper (preferably dark). This would increase the ice–background contrast and generate higher peaks of the intensity gradients on the ice edge.
The best practices for the scan detection technique focus mainly on generating manual tracings that are as accurate as possible. As the experimenter’s experience is crucial for reducing human error, it is advised that the person who will take the manual ice tracings practices the entire ice profile generation process before the real experiments. It is recommended that the same experimenter takes all of the manual tracings for an experimental campaign. Also, during the initial manual tracing, it is recommended to keep the pencil inclination on the paper constant. This should avoid artificial thickness variation in the ice profile.
Regarding the pre-processing step of redrawing the manual tracing with the green pencil, a sharp pencil should be used in order to avoid additional precision loss and dampening of the ice profile’s features. Moreover, the resolution of the scanner and the printer could produce additional dampening of the ice profile, increasing the measurement uncertainty. Therefore, a high scanner and printer resolution would be beneficial. The accuracy could be improved by scanning the original ice tracing only once and performing all of the following steps on the computer. Ideally, instead of manually re-tracing the original shape in color, a graphic editor would be used to change the tracings’ color in the scanned drawing. This would partially reduce human error, and thus, likely improve the profile accuracy and reduce the amount of post-processing code adjustments. If conducted manually, this step would significantly increase the amount of time required. Automation might be possible, but it would be very difficult, since the code would have to find the correct tracing automatically when multiple lines were drawn during the initial ice tracing process.
These are the best practices and recommendations for maximizing the efficiency and accuracy of the two 2D ice profile measurement techniques. Most of the sources of error are similar for all different icing regimes. Thus, the two methods presented in this paper are expected to perform equally for rime accretions and glaze accretions. For most established techniques, accurately measuring rime ice shapes is easier, because the geometry is typically easier.
Below, a comparison of the techniques with the ones existing in literature is conducted.
The IP technique has a higher accuracy and shorter acquisition and processing time compared to the manual tracing method. Its absolute errors range between −0.4 and +0.7 mm, while, for manual tracing, the absolute error range computed against caliper measurements and reported in Reference [
9] is from −1.6 to +2.5 mm. The IP technique is more expensive than manual tracing, as it implies the costs of a camera and a MATLAB license. However, its cost could be reduced by implementing the algorithm in an open-source programming language. Compared to the manual tracing method, the SD technique has a similar absolute error range from −1.2 mm to +2.6 mm. This was expected, as it is based on manual tracing profiles. However, the post-processing time of the SD method is automatic, which reduces human error, is faster, and is less cumbersome than the manual tracing technique, although the SD method is more expensive due to the required MATLAB license.
A direct comparison of the IP method accuracy was not possible for the photogrammetry technique, as no exact accuracy value was found in Reference [
9]. However, besides missing the 3D capability, the IP technique is less expensive than photogrammetry because it does not require additional layers of paint/powder. It also has a significantly reduced testing time, mainly because it does not require a drying process. Also, there are fewer images to capture when using the IP method. Compared to photogrammetry, the conclusions for the IP technique are similar to the SD method.
As in the case of photogrammetry, no accuracy values were found for ice accretion measurements using the molding technique. However, based on the techniques’ description provided in Reference [
25], good accuracy is reported for molding methods. It is a cumbersome testing process, though, during which, depending on the molds’ material, the tested ice piece might be damaged if not handled carefully. This, combined with the long amount of time needed for the mold material to cure (several hours), makes the IP method faster, cheaper, and simpler to use. However, the IP technique has no 3D capability. The conclusions about the molding versus the IP technique also apply to the comparison of molding with the SD technique.
Compared to the structured light technique [
28], which, in the cited study, had a spatial resolution window of 1 × 1 mm and an ice thickness measurement uncertainty of
0.04 mm, the IP technique has a higher average absolute error of ~0.26 mm. It might capture steep ice thickness variations better, though, as it generates pointwise results. The IP technique is faster to set up (because it does not require a prior parametric study and uses one camera only), is cheaper, and has a simpler post-processing method. However, the IP technique cannot generate 3D ice profiles or dynamic measurements, as is the case with the technique described in Reference [
28]. The SD method has a pointwise average absolute error of ~0.98 mm, which represents a much lower accuracy than the structured light technique, but it has the advantages mentioned for the IP technique. Additionally, both the IP and SD techniques are more versatile than the structured light technique, as they do not require a complex experimental setup or optical access to the IWT’s test section, and they could even be adapted for small-scale field measurements.
Compared to the 3D laser scanning technique presented in Reference [
32], which employed a Romer 7530SI scanner with a resolution of 0.05 mm, the IP method has a one-order-of-magnitude lower accuracy (~0.26 mm), and it is not capable of generating 3D ice profiles. Another drawback of the IP technique would be that, in order to generate continuous-line 2D profiles, an ice profile reconstruction step is required, which might further decrease the accuracy in some cases. However, for 2D profiles, it is very low-cost and fast compared to the cumbersome data acquisition and processing described in Reference [
32]. Compared to the 3D laser scanning, the SD method has almost two-orders-of-magnitude reduced accuracy (~0.98 mm), it generates only continuous line 2D profiles, it is cheaper, faster, and has a simpler acquisition and post-processing method.
A limitation of both methods presented in this study is that they only generate 2D ice profiles and not 3D profiles. Generating accurate 3D profiles is more valuable, especially for glaze ice, which produces more variation in the spanwise direction. However, even though technological progress allows for more numerical simulations of 3D profiles than in the past, simulating 2D ice shapes is still common [
55,
56]. Additionally, conducting experiments on the aerodynamic performance in conventional wind tunnels is also not uncommon using 2D ice shapes extruded in the spanwise direction [
10,
57,
58,
59].
Regarding the third post-processing technique presented in this paper, the 1D thickness ice profile generator, its applicability could be further developed. The technique allows for a more quantitative ice profile investigation than 2D ice profiles.
Section 3.3 presents the proof of concept for the possibility of identifying the ice accretion regimes based on the 1D ice thickness profiles and easier identification of the maximum ice thickness position along the airfoil profile. It is believed that these profiles could also be used for fast estimations of the accumulated ice mass. Further validation of these concepts is necessary, but given the proof of concept, the technique could be a promising tool for estimating the ice accretion regimes together with the angle of attack at which they formed for field tests or after in-flight icing events in which some of the atmospheric conditions might be unknown. This information could be estimated by an examination of the 1D ice thickness profile and might be useful for damage reports, research and development phase of ice protection systems, reverse engineering of atmospheric conditions in case of sensor failure, etc.