Experimental Investigation of the Wake and the Wingtip Vortices of a UAV Model

An experimental investigation of the wake of an Unmanned Aerial Vehicle (UAV) model using flow visualization techniques and a 3D Laser Doppler Anemometry (LDA) system is presented in this work. Emphasis is given on the flow field at the wingtip and the investigation of the tip vortices. A comparison of the velocity field is made with and without winglet devices installed at the wingtips. The experiments are carried out in a closed-circuit subsonic wind tunnel. The flow visualization techniques include smoke-wire and smoke-probe experiments to identify the flow phenomena, whereas for accurately measuring the velocity field point measurements are conducted using the LDA system. Apart from the measured velocities, vorticity and circulation quantities are also calculated and compared for the two cases. The results help to provide a more detailed view of the flow field around the UAV and indicate the winglets’ significant contribution to the deconstruction of wing-tip vortex structures.


Introduction
In recent years, there is an increasing trend in the development and use of fixed-wing Unmanned Aerial Vehicles (UAVs) by both military forces and civilian organizations [1].They are ideal solutions for a wide range of operations, such as fire detection, search and rescue, coastline and sea-lane monitoring, and security surveillance [2].Due to the absence of crew on-board, they present several advantages, such as the reduced operational cost, the ability to operate under hazardous conditions, and the increased flight endurance, which is one the most important characteristics when it comes to the aforementioned missions.In the case of a UAV, where the endurance is only limited by the available fuel on-board, it is essentially up to the aerodynamic efficiency of the configuration to maximize the flight time.For a given mission, the required lift force is defined, therefore the optimization process is essentially keeping the corresponding drag force as low as possible.
For a UAV that operates in the low subsonic, incompressible regime, as is the case with most Medium-Altitude-Long-Endurance (MALE) UAVs [3][4][5][6], a drag breakdown analysis was conducted in [7], which shows that the part that contributes the most to the total drag force of the aerial vehicle is the main wing.More specifically, when exposed in subsonic, incompressible flow, a finite wing's drag is a combination of profile drag and induced drag.The profile drag is the drag due to the shape of the body (skin friction + flow separation effects), whereas the induced drag is a result of the pressure imbalance at the tip of a finite wing between its upper (suction side) and lower (pressure side) surfaces.At the tip though, this imbalance causes the high-pressure air from the lower side to move upwards, where the pressure is lower, leading to the formation of a vortex, i.e., the wingtip vortex (Figure 1).This three-dimensional motion alters the flow field above the entire wing as well, thus resulting in the appearance of the induced drag term [8].Due to the importance of the induced drag to the total drag force budget, a lot of effort has been made in the previous decades to minimize this component and pertain the wingtip vortex.The most notable achievement is probably the winglet concept, introduced by Whitcomb [9].Having been studied by various researchers over the years, for issues that have to do with aerodynamic and structural issues, using both experimental and computational tools for design and optimization purposes [10][11][12][13][14], it is essentially a sophisticated type of wingtip fence that alters the flow at the tip in a manner that increases aerodynamic efficiency and reduces the size of the tip vortex.
drag to the total drag force budget, a lot of effort has been made in the previous decades to minimize this component and pertain the wingtip vortex.The most notable achievement is probably the winglet concept, introduced by Whitcomb [9].Having been studied by various researchers over the years, for issues that have to do with aerodynamic and structural issues, using both experimental and computational tools for design and optimization purposes [10][11][12][13][14], it is essentially a sophisticated type of wingtip fence that alters the flow at the tip in a manner that increases aerodynamic efficiency and reduces the size of the tip vortex.A lot of experimental studies exist that deal with the flow over wings and tip vortices.Even though the configurations vary, the flow development is based on the same fundamental principles, which means that the information from the existing literature can serve as a solid reference point.For instance, Shekarriz et al. [15] studied the near field behavior of a wing tip vortex, yielding useful conclusions regarding the velocity distributions and the roll-up process.Using the data provided by Hoffmann and Joubert [16], Nielsen and Schwind [17] divided the vortex into three regions based on an analogy with turbulent boundary layers: the vortex core (inner) region, the logarithmic region, and the defect law region.Based on the equations describing these regions, Corsiglia et al. [18] conducted three-dimensional hot-wire anemometer measurements in the far field, providing modified empirical equations and constants that describe the circulation distributions, as well as information regarding the meandering (lateral movement) of the wingtip vortex core in space as a function of time.This phenomenon, known as "vortex wandering", has also been experimentally investigated by various researchers over time [19][20][21][22][23]. Huang and Lin [24] experimentally studied the flow patterns and characteristics of vortex shedding on a wing with a NACA 0012 airfoil section using smoke-wire and surface oil-flow techniques for visualization, and hot-wire anemometry to characterize the frequency domain of the unsteady flow structures.The information on the frequency domain can be useful when it comes to estimating the expected frequencies and setting the acquisition time for reliable mean values, as is the case in the current work.Pino et al. [25] studied the velocity field behind a wingtip vortex and compared the results with trailing vortex theoretical models, whereas in a recent study Serrano-Aguilera et al. [26] developed a theoretical model in order to characterize the wingtip vortices in the near field using smoke-wire visualizations in combination with a laser sheet.Their results were used as a first validation in the flow visualizations that were employed in this paper.In a comparative study, Zheng and Ramaprian [27] conducted a detailed study of the vortex of a rectangular wing, comparing the flowfield of a stationary wing to the case where a sinusoidal pitch oscillation is induced.Moreover, Elsayed et al. [28] conducted Particle Image Velocimetry (PIV) studies, at flapped and unflapped wing configurations, and compared the velocity fields, vorticity contours, and circulation charts.In a very interesting approach, they showed how various vortex characteristics depend directly on the angle of attack, emphasizing on peak tangential velocity, vorticity, and circulation distribution.Finally, Muthusamy [29] carried out experiments in a A lot of experimental studies exist that deal with the flow over wings and tip vortices.Even though the configurations vary, the flow development is based on the same fundamental principles, which means that the information from the existing literature can serve as a solid reference point.For instance, Shekarriz et al. [15] studied the near field behavior of a wing tip vortex, yielding useful conclusions regarding the velocity distributions and the roll-up process.Using the data provided by Hoffmann and Joubert [16], Nielsen and Schwind [17] divided the vortex into three regions based on an analogy with turbulent boundary layers: the vortex core (inner) region, the logarithmic region, and the defect law region.Based on the equations describing these regions, Corsiglia et al. [18] conducted three-dimensional hot-wire anemometer measurements in the far field, providing modified empirical equations and constants that describe the circulation distributions, as well as information regarding the meandering (lateral movement) of the wingtip vortex core in space as a function of time.This phenomenon, known as "vortex wandering", has also been experimentally investigated by various researchers over time [19][20][21][22][23]. Huang and Lin [24] experimentally studied the flow patterns and characteristics of vortex shedding on a wing with a NACA 0012 airfoil section using smoke-wire and surface oil-flow techniques for visualization, and hot-wire anemometry to characterize the frequency domain of the unsteady flow structures.The information on the frequency domain can be useful when it comes to estimating the expected frequencies and setting the acquisition time for reliable mean values, as is the case in the current work.Pino et al. [25] studied the velocity field behind a wingtip vortex and compared the results with trailing vortex theoretical models, whereas in a recent study Serrano-Aguilera et al. [26] developed a theoretical model in order to characterize the wingtip vortices in the near field using smoke-wire visualizations in combination with a laser sheet.Their results were used as a first validation in the flow visualizations that were employed in this paper.In a comparative study, Zheng and Ramaprian [27] conducted a detailed study of the vortex of a rectangular wing, comparing the flowfield of a stationary wing to the case where a sinusoidal pitch oscillation is induced.Moreover, Elsayed et al. [28] conducted Particle Image Velocimetry (PIV) studies, at flapped and unflapped wing configurations, and compared the velocity fields, vorticity contours, and circulation charts.In a very interesting approach, they showed how various vortex characteristics depend directly on the angle of attack, emphasizing on peak tangential velocity, vorticity, and circulation distribution.Finally, Muthusamy [29] carried out experiments in a configuration with and without winglets at a low speed windtunnel.The study was limited to force balance measurements though, and no information was provided regarding the flow development and the effect the winglet has on the tip vortex.To the best of our knowledge, no study exists where a direct comparison is made in terms of flow field development between a configuration with, and without a winglet device installed at the wingtip, and few of the existing studies on wingtip vortices deal with UAV platforms.However, and since in UAV development the time and budget is limited when compared to a commercial airliner, the wings are of less complex geometry (in terms of airfoil profiles, twist angle etc.) and relatively simple to manufacture.Thus, the tips are highly loaded and are ideal platforms for winglet studies [11].
In the current paper, an experimental study of the flow field around a UAV model is presented.The aim, and main contribution, is to address the two aforementioned remarks.That is, to investigate the flowfield at the wing tip of the UAV and to quantify the effect the winglets have on the tip vortex.The study is part of the Hellenic Civil Unmanned Air Vehicle (HCUAV) design study that was carried out at the Laboratory of Fluid Mechanics and Turbomachinery (LFMT) at the Department of Mechanical Engineering, at the Aristotle University of Thessaloniki (AUTH) [6,7].A detailed presentation of the experimental studies is made, including flow visualization techniques and velocity measurements.Smoke-wire, laser sheet, and smoke-probe techniques are used for the flow visualization investigations, whereas point measurements by means of Laser Doppler Anemometry (LDA) are performed in the near wake.Two different wingtip configurations are used, one with-and one without the winglet device of the HCUAV, to compare the velocity field for the two cases.The winglet configuration was a result of a parametric CFD study [7] and its characteristics are given in Figure 2.
Aerospace 2017, 4, 53 3 of 17 configuration with and without winglets at a low speed windtunnel.The study was limited to force balance measurements though, and no information was provided regarding the flow development and the effect the winglet has on the tip vortex.To the best of our knowledge, no study exists where a direct comparison is made in terms of flow field development between a configuration with, and without a winglet device installed at the wingtip, and few of the existing studies on wingtip vortices deal with UAV platforms.However, and since in UAV development the time and budget is limited when compared to a commercial airliner, the wings are of less complex geometry (in terms of airfoil profiles, twist angle etc.) and relatively simple to manufacture.Thus, the tips are highly loaded and are ideal platforms for winglet studies [11].
In the current paper, an experimental study of the flow field around a UAV model is presented.The aim, and main contribution, is to address the two aforementioned remarks.That is, to investigate the flowfield at the wing tip of the UAV and to quantify the effect the winglets have on the tip vortex.The study is part of the Hellenic Civil Unmanned Air Vehicle (HCUAV) design study that was carried out at the Laboratory of Fluid Mechanics and Turbomachinery (LFMT) at the Department of Mechanical Engineering, at the Aristotle University of Thessaloniki (AUTH) [6,7].A detailed presentation of the experimental studies is made, including flow visualization techniques and velocity measurements.Smoke-wire, laser sheet, and smoke-probe techniques are used for the flow visualization investigations, whereas point measurements by means of Laser Doppler Anemometry (LDA) are performed in the near wake.Two different wingtip configurations are used, one with-and one without the winglet device of the HCUAV, to compare the velocity field for the two cases.The winglet configuration was a result of a parametric CFD study [7] and its characteristics are given in Figure 2.

Materials and Methods
The experiments were conducted in a closed-circuit, low-speed windtunnel facility at the LFMT.The wind tunnel has an 1810 mm × 600 mm × 600 mm test section (Figure 3) and is capable of producing a maximum speed of 50 m/s.The windtunnel blower speed is controlled by a frequency inverter.A smoke generator is placed inside the closed-loop wind tunnel, away from the test section.
Considering the models used in the study, two scaled variants (1:22) of the HCUAV were used, one without and one with the winglet devices installed at the tip.The models were manufactured by means of 3D printing to ensure accurate representation of the original geometry.Specifically, the main wing has a mean chord of 36.5 mm, a tip chord (c) of 21 mm, an Aspect Ratio of 8, and a NASA NLF(1)-1015 airfoil [30], whereas the wing span is 302 mm.A custom mechanism with a supporting rod is used to place the model at the center of the test section area, away from the windtunnel walls, and ensures that the model is accurately positioned within ± 0.5 deg (Figure 4).The model maximum cross-sectional area at the largest examined angle of attack (15 deg) divided by the area of the wind-tunnel cross-section is equal to 1.05%.Hence, no corrections due to the blockage effect were made [31].

Materials and Methods
The experiments were conducted in a closed-circuit, low-speed windtunnel facility at the LFMT.The wind tunnel has an 1810 mm × 600 mm × 600 mm test section (Figure 3) and is capable of producing a maximum speed of 50 m/s.The windtunnel blower speed is controlled by a frequency inverter.A smoke generator is placed inside the closed-loop wind tunnel, away from the test section.A series of flow visualization experiments were initially conducted on both UAV models for the flow phenomena to be identified and for the measured planes to be defined.More specifically, both smoke-probe and smoke-wire techniques were used to investigate the flow at the wingtip, at a wide range of angles of attack.Regarding the smoke-probe, a 3 mm diameter plastic probe is used, mounted on a set of copper pipes, through which the smoke was channeled.The copper pipes set is attached to a three-axes traversing mechanism featuring high accuracy and repeatability (with 4 μm standard deviation) in order for the smoke probe to be accurately positioned.The mechanism consists of a driving bolt connected to a stepper motor for each axis.Finally, a flexible rubber pipe connects the copper pipes with the smoke generator.To carry out the smoke-wire study, a custom frame is used.Specifically, a 0.3 mm diameter wire is mounted on a rigid aluminum frame, and a glass syringe is adjusted on top, where the glycerin-water mixture was stored.The mixture consists of 30% glycerin and 70% distilled water and drips on the wire through the glass syringe.Proper positioning of the syringe ensures the right amount of mixture dripping over the wire.Αn external power supply is used to charge the wire, thus creating a close circuit, in order for the mixture on the wire to be burnt.The entire frame is in turn placed in the test section, close to the model.The photographs were taken using a Nikon D3000 DSLR camera and Nikorr 35 mm, f 1.8 and Nikorr 18-105 mm, f 3.5-5.6lenses.
Regarding the LDA measurements, and since the flow field is symmetric, only half the wake was examined, so that the experiments were conducted at the right wing of the model.The relative Considering the models used in the study, two scaled variants (1:22) of the HCUAV were used, one without and one with the winglet devices installed at the tip.The models were manufactured by means of 3D printing to ensure accurate representation of the original geometry.Specifically, the main wing has a mean chord of 36.5 mm, a tip chord (c) of 21 mm, an Aspect Ratio of 8, and a NASA NLF(1)-1015 airfoil [30], whereas the wing span is 302 mm.A custom mechanism with a supporting rod is used to place the model at the center of the test section area, away from the windtunnel walls, and ensures that the model is accurately positioned within ±0.5 deg (Figure 4).The model maximum cross-sectional area at the largest examined angle of attack (15 deg) divided by the area of the wind-tunnel cross-section is equal to 1.05%.Hence, no corrections due to the blockage effect were made [31].A series of flow visualization experiments were initially conducted on both UAV models for the flow phenomena to be identified and for the measured planes to be defined.More specifically, both smoke-probe and smoke-wire techniques were used to investigate the flow at the wingtip, at a wide range of angles of attack.Regarding the smoke-probe, a 3 mm diameter plastic probe is used, mounted on a set of copper pipes, through which the smoke was channeled.The copper pipes set is attached to a three-axes traversing mechanism featuring high accuracy and repeatability (with 4 μm standard deviation) in order for the smoke probe to be accurately positioned.The mechanism consists of a driving bolt connected to a stepper motor for each axis.Finally, a flexible rubber pipe connects the copper pipes with the smoke generator.To carry out the smoke-wire study, a custom frame is used.Specifically, a 0.3 mm diameter wire is mounted on a rigid aluminum frame, and a glass syringe is adjusted on top, where the glycerin-water mixture was stored.The mixture consists of 30% glycerin A series of flow visualization experiments were initially conducted on both UAV models for the flow phenomena to be identified and for the measured planes to be defined.More specifically, both smoke-probe and smoke-wire techniques were used to investigate the flow at the wingtip, at a wide range of angles of attack.Regarding the smoke-probe, a 3 mm diameter plastic probe is used, mounted on a set of copper pipes, through which the smoke was channeled.The copper pipes set is attached to a three-axes traversing mechanism featuring high accuracy and repeatability (with 4 µm standard deviation) in order for the smoke probe to be accurately positioned.The mechanism consists of a driving bolt connected to a stepper motor for each axis.Finally, a flexible rubber pipe connects the copper pipes with the smoke generator.To carry out the smoke-wire study, a custom frame is used.Specifically, a 0.3 mm diameter wire is mounted on a rigid aluminum frame, and a glass syringe is adjusted on top, where the glycerin-water mixture was stored.The mixture consists of 30% glycerin and 70% distilled water and drips on the wire through the glass syringe.Proper positioning of the syringe ensures the right amount of mixture dripping over the wire.An external power supply is used to charge the wire, thus creating a close circuit, in order for the mixture on the wire to be burnt.The entire frame is in turn placed in the test section, close to the model.The photographs were taken using a Nikon D3000 DSLR camera and Nikorr 35 mm, f 1.8 and Nikorr 18-105 mm, f 3.5-5.6lenses.
Regarding the LDA measurements, and since the flow field is symmetric, only half the wake was examined, so that the experiments were conducted at the right wing of the model.The relative position for each measurement plane is presented in Figure 5.The first plane lies half a tip chord (X/c = 0.5), the second plane one tip chord (X/c = 1), and the third planes lies two tip chords (X/c = 2) downstream of the wingtip, respectively.The measurement grid at each of the planes is the same, and has 625 nodes in a 25 × 25 arrangement, covering a 50 mm × 50 mm area.In order to increase the accuracy, the grid is denser at the location of the vortex, as indicatively presented for the regular wingtip case in Figure 6.The starting point of the coordinate system at each measurement was located at tip of the wing/winglet for all measurements.To conduct the measurements a complete 3-component LDA system from Dantec Dynamics (Skovlunde, Denmark) is used, which consists of an Ar-ion Beam Generator, an optical system with 3D component Probes, a 3-axes traverser, and a dedicated commercial software (BSA Software© by Dantec Dynamics A/S).The Traverser on which the 3D component Probes are mounted, is placed next to the test-section, with the Probes facing the model.Using the smoke generator, a smoke cloud is generated throughout the entire path of the windtunnel, thus inducing the particles in the flow that are needed for the LDA Measurements.Since no apparent frequencies were expected in the lower frequency range [18,21,22], the acquisition time was set at 5 s for each node, adequately long for providing reliable mean values.
Finally, regarding the inlet conditions, for the LDA measurements the freestream velocity at the test section is 13.2 m/s, held constant within 1%, which gives a Reynolds Number equal to 33 × 10 3 based on the mean aerodynamic chord, resembling existing wingtip vortex-related studies [26,28].The freestream turbulence intensity was measured at 0.7%.For the flow visualization studies, the freestream velocity is 2 m/s to make sure that the smoke is clearly visible [32].Due to the difference in freestream velocity and, consequently, in the corresponding Reynolds number, no direct comparison can be made between the flow visualization results and the LDA measurements.The visualization studies are critical though, to get a first view of the flowfield and accurately position and size the LDA measurements grids.The lateral velocity components are insignificant in all studies.
Aerospace 2017, 4, 53 5 of 17 position for each measurement plane is presented in Figure 5.The first plane lies half a tip chord (X/c = 0.5), the second plane one tip chord (X/c = 1), and the third planes lies two tip chords (X/c = 2) downstream of the wingtip, respectively.The measurement grid at each of the planes is the same, and has 625 nodes in a 25 × 25 arrangement, covering a 50 mm × 50 mm area.In order to increase the accuracy, the grid is denser at the location of the vortex, as indicatively presented for the regular wingtip case in Figure 6.The starting point of the coordinate system at each measurement was located at tip of the wing/winglet for all measurements.To conduct the measurements a complete 3component LDA system from Dantec Dynamics (Skovlunde, Denmark) is used, which consists of an Ar-ion Beam Generator, an optical system with 3D component Probes, a 3-axes traverser, and a dedicated commercial software (BSA Software© by Dantec Dynamics A/S).The Traverser on which the 3D component Probes are mounted, is placed next to the test-section, with the Probes facing the model.Using the smoke generator, a smoke cloud is generated throughout the entire path of the windtunnel, thus inducing the particles in the flow that are needed for the LDA Measurements.Since no apparent frequencies were expected in the lower frequency range [18,21,22], the acquisition time was set at 5 s for each node, adequately long for providing reliable mean values.Finally, regarding the inlet conditions, for the LDA measurements the freestream velocity at the test section is 13.2 m/s, held constant within 1%, which gives a Reynolds Number equal to 33 × 10 3 based on the mean aerodynamic chord, resembling existing wingtip vortex-related studies [26,28].

Results
At first, the results at the tip of the model with the regular wingtip (winglet device removed) are shown, followed by the results at the tip of the model with the winglet installed.Concluding, a comparison of the flow field between the two configurations is made.The contour plots show the flow field in the wake of the HCUAV right wing, assuming that the view point of the reader is downstream of the aerial vehicle, whereas line plots and other charts are also shown to give a better understanding of the phenomena inside the core.The normal and spanwise axes have been normalized using the tip chord of the model (c), and the time-averaged velocities have been normalized using the freestream velocity (U inf ).Moreover, the vorticity and circulation quantities are normalized by appropriately using the tip chord and freestream velocity.

Regular Wingtip
Figure 7 presents the results of the smoke-probe investigations (top view), whereas Figure 8 presents the results of the combined smoke-wire/laser sheet method, at one tip chord downstream of the wing (X/c = 1).The corresponding planes and direction of freestream velocity are also given for each case.As is evident, at 0 deg of angle of attack, the tip vortex is practically non-existent, since the lift production is negligible and the pressure imbalance between the two sides is very small.At −5 and 5 deg, the vortex is more clearly defined, especially in the laser sheet visualization.At −5 deg of angle of attack the vortical structure is upside down, since the lift production is inverse.At 10 deg of angle of attack, the tip vortex appears to be even larger in size and robust, whereas at 15 deg, the vortical structure retains its size but fades quickly as the flow moves downstream of the wing.In general, the flow visualization studies may not directly provide numbers about the examined phenomenon, but they helped a long way in identifying the key-characteristics of the vortex at each angle of attack, i.e., size and core location.Thus, the LDA measurement grids could be properly sized with greater accuracy, at a reduced time.The LDA results are presented in Figures 9 and 10 in contour plots and in Figures 11-13 in line plots.Note that the vortex core radius is defined as the radius at which the swirl velocity component maximizes, whereas the axial vorticity was calculated using the experimental data from the swirl velocity components, and more specifically where V is the normal and W the spanwise, time-averaged velocity components.In Figure 9, the time-averaged axial, normal and spanwise velocity components and axial vorticity are presented in contour plots at X/c = 1 and three different angles of attack, whereas in Figure 10, the same variables are shown for the three measurement levels at 10 deg of angle of attack in contour plots at the three measurement locations.The freestream velocity direction is identical the one presented in Figure 8.When observing the axial velocity component (U), it is evident that the wake at the wingtip can be divided at two separate regions.One region is behind the main wing, which corresponds to the typical wake profile observed behind a blunt body, as a result of shear layer separation, and the other corresponds to the wingtip vortex.In the wing-wake region, U is of smaller order of magnitude when compared to the vortex core one, with the exact values varying greatly between the different angles of attack and measurement levels.Regarding the time-averaged normal and spanwise velocity components (V and W), two maxima of opposite sign are observed around a point where V and W equal to zero.This point coincides with the one where the local minimum of the U component is located (U 0 ) at every measurement location.Hence, that point was defined as the center of the vortex.At every measurement plane, the V magnitude increases considerably as the angle of attack grows (Figure 9), but for the same angle of attack, as the flow moves downstream of the tip, the difference in magnitude is smaller (Figure 10).The same observations can be made for the spanwise velocity component.Overall, due to the fact that the W component itself does not add any new information on the phenomenon, and since the axial vorticity includes both of the V and W components, it was decided to not present the spanwise component in the rest of the figures of this work.The time-averaged axial velocity component inside the vortex core can be more clearly seen in Figure 11, where line plots show the normalized values of U. It should be noted that for each measurement the line plot sections pass through the corresponding vortex center.The time-averaged axial velocity component appears to form an axisymmetric wake-like profile.This axisymmetric behavior is also observed in Figures 12 and 13, which indicate that the roll-up process is complete for X/c < 1, as also reported in [13,25].Namely, Figure 12 presents the V component in line plot representation.At 15 deg of angle of attack, the range of the swirl velocity magnitude is approximately three times higher than the one measured at 0 deg, indicating that the tip vortex grows in strength.Regarding the change in magnitude between different longitudinal positions, an average decrease of 10% was calculated between X/c = 1 and X/c = 0.5 for every examined angle of attack, whereas between X/c = 2 and X/c = 1, the percentage average is approximately 13%.Moreover, it should be noted that behind the main wing geometry, the normal velocity component has negative values outside the vortex core.This is probably due to the downwash effect and the corresponding velocity component that appears in the normal axis.The third parameter in the contour plots, i.e., the axial vorticity, is also presented in line plots in Figure 13.At every measurement level the vorticity maximum is approximately three times larger at 15 deg than the one at 0 deg.Between 10 and 15 deg of angle of attack, though, the difference in peak vorticity is very small, indicating that the lift increases at a much lower rate and, possibly, that stalling occurs.As a general observation, the peak vorticity values follow the trend of the normal velocity, having a 10-15% difference from level to level, at each angle of attack, an observation that is in agreement with the data from [28].However, the line plots provide the information in a more detailed way, so that a direct comparison with another experimental or computational data, can be made.A direct comparison between the three measurement levels is made in Figures 14 and 15.More specifically, Figure 14 presents the time-averaged U and V components, vorticity and circulation at the vortex core at an angle of attack of 10 deg.Regarding the velocities, it is evident that the normal component along with the axial vorticity have a much smaller change in magnitude than the axial velocity, between the three streamwise locations.Namely, a 10% and 13% change in V and vorticity was measured, respectively, from one level to the next one, which is in agreement with the comments made in the previous paragraph.The vortex core is located at a radius of 0.07•c, around the center of the vortex, the same for all three measurement planes at this angle of attack.The circulation was also calculated for comparison purposes, as the flux of vorticity, through integration.The following equation was used: where Ω is the axial vorticity from Equation (1) and S is the surface.Since it was more convenient, polar coordinates were employed for this type of chart, with the starting point of the coordinate point being the center of the vortex.The typical three-region distinction can be made [14,16], with the slope in the logarithmic region decreasing as the vortex gets carried downstream.The changes are in the same order of magnitude with the ones at normal velocity and vorticity, that is between 10% and 15% from level to level.
In Figure 15, the smaller value of axial velocity inside the vortex core (U0) is shown for the three measurement planes at every angle of attack.At each level, a linear interpolation can be made, indicating that U0 decreases as the angle of attack grows at every measurement plane.The magnitude of U0, though, increases as the planes move downstream, whereas the linear slope dU0/dα decreases.A direct comparison between the three measurement levels is made in Figures 14 and 15.More specifically, Figure 14 presents the time-averaged U and V components, vorticity and circulation at the vortex core at an angle of attack of 10 deg.Regarding the velocities, it is evident that the normal component along with the axial vorticity have a much smaller change in magnitude than the axial velocity, between the three streamwise locations.Namely, a 10% and 13% change in V and vorticity was measured, respectively, from one level to the next one, which is in agreement with the comments made in the previous paragraph.The vortex core is located at a radius of 0.07•c, around the center of the vortex, the same for all three measurement planes at this angle of attack.The circulation was also calculated for comparison purposes, as the flux of vorticity, through integration.The following equation was used: where Ω is the axial vorticity from Equation (1) and S is the surface.Since it was more convenient, polar coordinates were employed for this type of chart, with the starting point of the coordinate point being the center of the vortex.The typical three-region distinction can be made [14,16], with the slope in the logarithmic region decreasing as the vortex gets carried downstream.The changes are in the same order of magnitude with the ones at normal velocity and vorticity, that is between 10% and 15% from level to level.
In Figure 15, the smaller value of axial velocity inside the vortex core (U0) is shown for the three measurement planes at every angle of attack.At each level, a linear interpolation can be made, indicating that U0 decreases as the angle of attack grows at every measurement plane.The magnitude of U0, though, increases as the planes move downstream, whereas the linear slope dU0/dα decreases.A direct comparison between the three measurement levels is made in Figures 14 and 15.More specifically, Figure 14 presents the time-averaged U and V components, vorticity and circulation at the vortex core at an angle of attack of 10 deg.Regarding the velocities, it is evident that the normal component along with the axial vorticity have a much smaller change in magnitude than the axial velocity, between the three streamwise locations.Namely, a 10% and 13% change in V and vorticity was measured, respectively, from one level to the next one, which is in agreement with the comments made in the previous paragraph.The vortex core is located at a radius of 0.07•c, around the center of the vortex, the same for all three measurement planes at this angle of attack.The circulation was also calculated for comparison purposes, as the flux of vorticity, through integration.The following equation was used: where Ω is the axial vorticity from Equation (1) and S is the surface.Since it was more convenient, polar coordinates were employed for this type of chart, with the starting point of the coordinate point being the center of the vortex.The typical three-region distinction can be made [14,16], with the slope in the logarithmic region decreasing as the vortex gets carried downstream.The changes are in the same order of magnitude with the ones at normal velocity and vorticity, that is between 10% and 15% from level to level.

Winglet Installed
Figure 16 shows the time-averaged axial and normal velocity and axial vorticity in contour plots at X/c = 1 for 5, 10, and 15 deg of angle of attack.The freestream velocity direction is identical to the one presented in Figure 8.Based on the results of this study, two separate swirling motions can be observed, one around the vortex core, at the tip of the winglet, and another around the entire winglet configuration.The normal velocity magnitude is about two times larger inside the "vortex" region.However, a second region can be clearly observed around the entire winglet configuration, from now on referred to as the "winglet" region.In terms of vorticity, the difference is even greater, since the vorticity maximum inside the vortex core is at least three times as high, compared to the vorticity values of the "winglet" region.Hence, the "winglet" region indicates that another vortical structure exists, which may be far weaker in terms of magnitude but covers a larger area, since it spreads

Winglet Installed
Figure 16 shows the time-averaged axial and normal velocity and axial vorticity in contour plots at X/c = 1 for 5, 10, and 15 deg of angle of attack.The freestream velocity direction is identical to the one presented in Figure 8.Based on the results of this study, two separate swirling motions can be observed, one around the vortex core, at the tip of the winglet, and another around the entire winglet configuration.The normal velocity magnitude is about two times larger inside the "vortex" region.However, a second region can be clearly observed around the entire winglet configuration, from now on referred to as the "winglet" region.In terms of vorticity, the difference is even greater, since the vorticity maximum inside the vortex core is at least three times as high, compared to the vorticity values of the "winglet" region.Hence, the "winglet" region indicates that another vortical structure exists, which may be far weaker in terms of magnitude but covers a larger area, since it spreads In Figure 15, the smaller value of axial velocity inside the vortex core (U 0 ) is shown for the three measurement planes at every angle of attack.At each level, a linear interpolation can be made, indicating that U 0 decreases as the angle of attack grows at every measurement plane.The magnitude of U 0 , though, increases as the planes move downstream, whereas the linear slope dU 0 /dα decreases.

Winglet Installed
Figure 16 shows the time-averaged axial and normal velocity and axial vorticity in contour plots at X/c = 1 for 5, 10, and 15 deg of angle of attack.The freestream velocity direction is identical to the one presented in Figure 8.Based on the results of this study, two separate swirling motions can be observed, one around the vortex core, at the tip of the winglet, and another around the entire winglet configuration.The normal velocity magnitude is about two times larger inside the "vortex" region.However, a second region can be clearly observed around the entire winglet configuration, from now on referred to as the "winglet" region.In terms of vorticity, the difference is even greater, since the vorticity maximum inside the vortex core is at least three times as high, compared to the vorticity values of the "winglet" region.Hence, the "winglet" region indicates that another vortical structure exists, which may be far weaker in terms of magnitude but covers a larger area, since it spreads approximately all over the winglet wake region.Furthermore, both V and vorticity magnitude present very small differences between 10 and 15 deg, indicating that the lift generation has reached its peak, resembling the corresponding trends without the winglet installed at the tip.Regarding the axial velocity, and in a similar fashion to the regular wingtip case, the wake can be divided at two separate regions, i.e., the "blunt body wake" region and the "vortex" ones.The axial velocity in both regions is reduced in magnitude as the angle of attack increases.approximately all over the winglet wake region.Furthermore, both V and vorticity magnitude present very small differences between 10 and 15 deg, indicating that the lift generation has reached its peak, resembling the corresponding trends without the winglet installed at the tip.Regarding the axial velocity, and in a similar fashion to the regular wingtip case, the wake can be divided at two separate regions, i.e., the "blunt body wake" region and the "vortex" ones.The axial velocity in both regions is reduced in magnitude as the angle of attack increases.Moreover, a comparison of the aforementioned variables inside the vortex is made in Figure 17 for all of the examined angles of attack.The streamwise and normal time-averaged velocity and the axial vorticity line plots provide more detail about the development of these quantities inside the vortex at the tip of the winglet and can be used in order for a comparison with computational data to be made.It should be noted that the peak values of the normal velocity component and the axial vorticity are considerably smaller than the ones presented in Figures 11-13.A more detailed comparison between the two cases is made in Section 3.3.Moreover, a comparison of the aforementioned variables inside the vortex is made in Figure 17 for all of the examined angles of attack.The streamwise and normal time-averaged velocity and the axial vorticity line plots provide more detail about the development of these quantities inside the vortex at the tip of the winglet and can be used in order for a comparison with computational data to be made.It should be noted that the peak values of the normal velocity component and the axial vorticity are considerably smaller than the ones presented in Figures 11-13.A more detailed comparison between the two cases is made in Section 3.3.

Comparison
Having measured and analyzed the flowfield around both without and with the winglet installed, a direct comparison between the two cases can be made.More specifically, as is evident from the two previous sections, there is a considerable decrease in the size and magnitude of the wingtip vortex, when the winglet device is installed at the tip.This decrease can be qualitatively observed in Figure 18, where a top view of the smoke-probe visualization experiments is presented at the tip of the wing, with and without the winglet installed, at X/c = 1 and 10 deg of angle of attack.

Comparison
Having measured and analyzed the flowfield around both without and with the winglet installed, a direct comparison between the two cases can be made.More specifically, as is evident from the two previous sections, there is a considerable decrease in the size and magnitude of the wingtip vortex, when the winglet device is installed at the tip.This decrease can be qualitatively observed in Figure 18, where a top view of the smoke-probe visualization experiments is presented at the tip of the wing, with and without the winglet installed, at X/c = 1 and 10 deg of angle of attack.This comparison is also made in Figure 19 using LDA contour plots for time-averaged axial velocity and axial vorticity and in Figure 20 using line plots.The angle of attack of 10 deg at X/c = 1 was chosen as a reference point, since at that angle, the vortex appears to be strong and welldefined for both cases.In Figure 20a, the time-averaged axial velocity is shown, followed by the timeaveraged normal velocity component in Figure 20b.It is clear that, even though the axial component is of the same order of magnitude between the two cases, with a percentage difference smaller than 5%, the swirl component peak values are reduced with the addition of the winglet by a factor of 3.That can also be seen in Figure 20c, where the peak axial vorticity with winglet equals to 35.5% of the corresponding peak vorticity without winglet.In other words, the addition of the winglet seems to reshape the wake and tip vortex, causing the swirl velocity components and peak vorticity magnitude inside the vortex core to reduce.Finally, in Figure 20d the slope of the curve inside the logarithmic region is reduced to 28% with the addition of the winglet, whereas the maximum value is almost five times smaller.This comparison is also made in Figure 19 using LDA contour plots for time-averaged axial velocity and axial vorticity and in Figure 20 using line plots.The angle of attack of 10 deg at X/c = 1 was chosen as a reference point, since at that angle, the vortex appears to be strong and well-defined for both cases.In Figure 20a, the time-averaged axial velocity is shown, followed by the time-averaged normal velocity component in Figure 20b.It is clear that, even though the axial component is of the same order of magnitude between the two cases, with a percentage difference smaller than 5%, the swirl component peak values are reduced with the addition of the winglet by a factor of 3.That can also be seen in Figure 20c, where the peak axial vorticity with winglet equals to 35.5% of the corresponding peak vorticity without winglet.In other words, the addition of the winglet seems to reshape the wake and tip vortex, causing the swirl velocity components and peak vorticity magnitude inside the vortex core to reduce.Finally, in Figure 20d the slope of the curve inside the logarithmic region is reduced to 28% with the addition of the winglet, whereas the maximum value is almost five times smaller.
corresponding peak vorticity without winglet.In other words, the addition of the winglet seems to reshape the wake and tip vortex, causing the swirl velocity components and peak vorticity magnitude inside the vortex core to reduce.Finally, in Figure 20d the slope of the curve inside the logarithmic region is reduced to 28% with the addition of the winglet, whereas the maximum value is almost five times smaller.

Conclusions
The near flow field around a UAV model was experimentally investigated by means of flow visualization and 3D LDA point measurements.The experiments were carried out at a wide range of angles of attack to ensure that all the operational conditions are examined and emphasis was given on the flow at the wingtip area.Two configurations, one without and one with winglets installed, were used.The flow visualization studies were conducted using the smoke-probe, smoke-wire, and laser sheet methods.The LDA measurements yielded the three time-averaged velocities, whereas the axial vorticity and circulation quantities were in turn calculated using the experimental data.Thus, a complete experimental investigation of the wingtip vortex of a UAV could be made, whereas a direct

Conclusions
The near flow field around a UAV model was experimentally investigated by means of flow visualization and 3D LDA point measurements.The experiments were carried out at a wide range of angles of attack to ensure that all the operational conditions are examined and emphasis was given on the flow at the wingtip area.Two configurations, one without and one with winglets installed, were used.The flow visualization studies were conducted using the smoke-probe, smoke-wire, and laser sheet methods.The LDA measurements yielded the three time-averaged velocities, whereas the axial vorticity and circulation quantities were in turn calculated using the experimental data.Thus, a complete experimental investigation of the wingtip vortex of a UAV could be made, whereas a direct comparison of the flowfield without and with a winglet installed is made, serving as the main contribution of this study.
Considering the flow at the wingtip, a comparison at various angles of attack for three locations downstream of the model was made.The tip vortex does not appear to have a defined structure until 5 deg, whereas as the angle of attack increases the vortex grows in strength, the swirling components increase in magnitude, and the vorticity at the core increases as well.The wake seems to form two distinctive regions: one affected by the presence of the vortex, and another affected by the shear layer separation at the wake of the model.At every angle of attack, an average difference of 10-15% was measured, between the three levels for the time-averaged normal velocity and axial vorticity.
In the "winglet installed" case, the trends resemble those of the "winglet removed" case, regarding the changes between the various angles of attack.A second swirling motion, weaker in strength, was also identified, which in contrary to the wingtip vortex exists all over the winglet wake region.Overall, the addition of the winglet reshapes the tip vortex, causing the swirl velocity components and peak vorticity magnitude inside the vortex core to reduce.More specifically, the maximum vorticity magnitude in this study is reduced by a factor of three with the winglet installed.
The above observations help to acquire a better idea about the tip vortex of a UAV and about the effect the winglet device has on its dispersion.Concerning future studies, a turbulence budget investigation could also be conducted to provide even more details about the vortex core, whereas a flow visualization using a high-speed camera is suggested in order to acquire more information about the wandering of the vortex and the effect the winglet has on this phenomenon.

Figure 2 .
Figure 2. Hellenic Civil Unmanned Air Vehicle (HCUAV) configuration without and with the winglet installed, including the winglet's characteristics in the detail [7].

Figure 2 .
Figure 2. Hellenic Civil Unmanned Air Vehicle (HCUAV) configuration without and with the winglet installed, including the winglet's characteristics in the detail [7].

Figure 3 .
Figure 3. Schematic of the experimental setup and arrangement of the measurement devices inside the test section.

Figure 3 .
Figure 3. Schematic of the experimental setup and arrangement of the measurement devices inside the test section.

Figure 3 .
Figure 3. Schematic of the experimental setup and arrangement of the measurement devices inside the test section.

Figure 5 .
Figure 5. Measurement planes for the Laser Doppler Anemometry (LDA) measurements around the Hellenic Civil Unmanned Air Vehicle (HCUAV) model.

Figure 5 .
Figure 5. Measurement planes for the Laser Doppler Anemometry (LDA) measurements around the Hellenic Civil Unmanned Air Vehicle (HCUAV) model.

Figure 5 .
Figure 5. Measurement planes for the Laser Doppler Anemometry (LDA) measurements around the Hellenic Civil Unmanned Air Vehicle (HCUAV) model.

Figure 6 .
Figure 6.High resolution grid for the 3rd measurement plane (regular wingtip).The trailing edge is marked with dark red.

Figure 6 .
Figure 6.High resolution grid for the 3rd measurement plane (regular wingtip).The trailing edge is marked with dark red.

Figure 8 .
Figure 8. Laser sheet visualization studies at the wingtip.

Figure 8 .
Figure 8. Laser sheet visualization studies at the wingtip.Figure 8. Laser sheet visualization studies at the wingtip.

Figure 8 .
Figure 8. Laser sheet visualization studies at the wingtip.Figure 8. Laser sheet visualization studies at the wingtip.

Figure 9 .
Figure 9. Time-averaged axial, normal and spanwise velocity, and vorticity contours at the wingtip.Comparison between various angles of attack.

Figure 9 .
Figure 9. Time-averaged axial, normal and spanwise velocity, and vorticity contours at the wingtip.Comparison between various angles of attack.

Figure 10 .
Figure 10.Time-averaged axial, normal and spanwise velocity and vorticity contours at the wingtip.Comparison between the various measurement levels.

Figure 11 .
Figure 11.Line plots of time-averaged axial velocity at the vortex core.

Figure 10 . 17 Figure 10 .
Figure 10.Time-averaged axial, normal and spanwise velocity and vorticity contours at the wingtip.Comparison between the various measurement levels.

Figure 11 .
Figure 11.Line plots of time-averaged axial velocity at the vortex core.Figure 11.Line plots of time-averaged axial velocity at the vortex core.

Figure 11 .
Figure 11.Line plots of time-averaged axial velocity at the vortex core.Figure 11.Line plots of time-averaged axial velocity at the vortex core.

Figure 12 .
Figure 12.Line plots of time-averaged normal velocity at the vortex core.

Figure 13 .
Figure 13.Line plots vorticity for the three measurement planes.

Figure 12 . 17 Figure 12 .
Figure 12.Line plots of time-averaged normal velocity at the vortex core.

Figure 13 .
Figure 13.Line plots vorticity for the three measurement planes.

Figure 13 .
Figure 13.Line plots vorticity for the three measurement planes.

Figure 15 .
Figure 15.Time-averaged axial velocity in the center of the vortex (U0) at the three measurement levels, without winglet.

Figure 15 .
Figure 15.Time-averaged axial velocity in the center of the vortex (U0) at the three measurement levels, without winglet.

Figure 15 .
Figure 15.Time-averaged axial velocity in the center of the vortex (U 0 ) at the three measurement levels, without winglet.