=Paper= {{Paper |id=Vol-3732/paper12 |storemode=property |title=Research of aerodynamic characteristics of aviation profile model in conditions of artificial generation of large longitudinal vortices |pdfUrl=https://ceur-ws.org/Vol-3732/paper12.pdf |volume=Vol-3732 |authors=Oleksander Zhdanov,Valerii Orlianskyi,Olha Sushchenko |dblpUrl=https://dblp.org/rec/conf/cmse/ZhdanovOS24 }} ==Research of aerodynamic characteristics of aviation profile model in conditions of artificial generation of large longitudinal vortices== https://ceur-ws.org/Vol-3732/paper12.pdf
                                Research of aerodynamic characteristics of aviation profile
                                model in conditions of artificial generation of large
                                longitudinal vortices
                                Oleksander Zhdanov1,†, Valerii Orlianskyi1,† and Olha Sushchenko1,*,†

                                1 National Aviation University, Liubomyra Huzara Ave., 1, Kyiv, 03058, Ukraine



                                                Abstract
                                                This paper solves the problem of widening the range of the angle of the attack of safe flight
                                                using volumetric vortex generators. The research method is based on experimental tests in the
                                                wind tunnel. Analyzing the obtained results allows us to research changes in the integrated
                                                aerodynamic characteristics. The features of the experimental test are described. The
                                                volumetric vortex generators of three types for the definite blowing model are represented.
                                                Comparative changing of aerodynamic characteristics of the profile model in the form of
                                                graphical dependencies has been done. The visualization of airflow in an experimental way has
                                                been carried out.

                                                Keywords
                                                vortex, aerodynamic characteristics, angle of attack, flight safety, profile model1



                                1. Introduction
                                The problem of widening the range of changing attack angles for the safe operation of
                                unmanned aerial vehicles (UAVs) is of considerable relevance in scientific areas and
                                practical application. It can be solved by eliminating stalling by controlling air flow
                                breakdown on control elements and the wing in flight modes with subsonic speeds. This is
                                implemented by installing volumetric vortex generators (VVG) on the wing front edge and
                                optimizing the geometric parameters of their shape and location.
                                   This article studies the possibility and efficiency of using an artificially organized
                                system of big longitudinal vortexes. The emphasis has been placed on modifying space-
                                time scales and developing a practical approach to creating large longitudinal vortices on
                                the surface of an aerodynamic profile model. Research studies in this direction have



                                CMSE’2024: International Workshop on Computational Methods in Systems Engineering, June 17, 2024, Kyiv,
                                Ukraine
                                ∗ Corresponding author.
                                † These authors contributed equally.

                                   azhdanov@nau.edu.ua (O. Zhdanov); aerodyn@nau.edu.ua (V. Orlianskyi); sushoa@ukr.net (O. Sushchenko)
                                    0000-0001-5273-571X (O. Zhdanov); 0009-0000-8186-4211 (V. Orlianskyi); 0000-0002-8837-1521
                                (O. Sushchenko)
                                         © 2024 Copyright for this paper by its authors. Use permitted under Creative Commons License Attribution 4.0 International (CC BY 4.0).




CEUR
                  ceur-ws.org
Workshop      ISSN 1613-0073
Proceedings
shown the effectiveness of using VVG and even led to using such a term as a vortex-active
wing or a vortex-active aerodynamic surface [1, 2].
    The research method is experimental weight tests of the profile model in the wind
tunnel and comparative analysis of their results in the form of changes in integral
aerodynamic characteristics.
    There are vortex generators (VG) with different shapes and dimensions. Prototypes of
VG could be characterized by various geometrical characteristics such as the chord, area,
volume, setting angle, and installation step. The area is determined as a projection in the
plane far off the profile forepart edge using relative characteristics. The setting angle can
be characterized by the chords of VG and profile, which form an angle. The installation
step is characterized by the model span.
    Forming a shape of VG, it is desirable to achieve the maximal efficiency of the
longitudinal vortex. For this, the longitudinal vortex must be arranged on the profile
model surface in correspondence with some requirements. It means that power
consumption of the external air flow must be minimized. Moreover, the length of the
vortex along the profile chord must be minimized also. It is necessary to eliminate the
longitudinal vortex backward in conditions of changing the pressure gradient on the
profile surface.
    The vortex source can be located at the point, which is contacting with the VVG and
aerodynamic surface simultaneously. Is noteworthy that the interference phenomenon
can create positive or negative conditions because it can create a decrease or increase in
front drag force. It is necessary to remember that the VG itself can be a source of the
vortex. As an example, we can consider components of the VG arising outside the profile
model and arranged at some angle relative to the airflow. Some valleys and ledges of the
VG surface also can be sources of vortexes. Hence, we deal with the possible negative
consequences of interference. To decrease these negative consequences and even to
improve the situation, we must ensure smooth transitions between the profile model and
VG [3, 4].
    The efficiency of generating vortexes depends on the attack angle. The influence of VG
must be minimal in cruising modes. And especially the possibility of arising critical attack
angles has a great significance. Therefore, we must locate the VG relative to a chord in a
definite way taking into consideration the chord curvature. It is possible also to improve
this situation by varying the angular location of VG depending on the airflow conditions
[5]. There is also the possibility to apply VG when achieving critical attack angles becomes
probable. This approach can be realized by two ways. Firstly, we can use VG sliding.
Secondly, compensation or reduction of front drag force caused by airfoil installation can
be implemented as follows. It is desirable to increase the suction force both on the surface
of the airfoil and the toe of the profile model. This can be done by forming the shape
curvature of the airfoil surface shape. This influences on the airflow acceleration in the
boundary layer on the surface of the airfoil and profile model toes [6].
    The number of factors that influence characteristics of aerodynamic surface
compositions with VVG, based on the results of research studies [7, 8] include the
following factors. They include the external shape of the VG, profile model dimensions, the
arrangement by the size of the peak far off the profile frontal edge, the inclination angle of
the chord of the VG relative to the chord of the profile, the step of setting the VG according
to the scale of the model, and others. Carrying out such multifactorial studies to determine
optimal layouts requires a large number of tests of layout models, analysis of results, and
determination of optimal layouts according to chosen criteria [9, 10].

2. Description of the experiment

The experiments were carried out in a closed-type low-velocity aerodynamic atmospheric
tube with an open working part of an elliptical cross-section with the dimensions of the
axes 750×420 mm and a length of 900 mm. The airflow speed was from 3.5 m/s to 30 m/s.
The air tube is equipped with an external three-component rider aerodynamic weight,
which provides the measurement of three components of the total aerodynamic force. The
measured forces allow us to calculate the dependence of dimensionless aerodynamic
coefficients on the attack angle. The built-in device for setting the angle location of the
profile model provides changing attack angles (-20° < α < 40°).
    The model suspended in the wind tunnel must have the possibility of circular blowing
at the angle of attack in separate areas within the available range of changing the angular
location. The mechanical system for decomposing the forces of the rider weight can be
connected to the system of tensor converters of loads on weight components into
electrical signals. The set of strain gauges provides automation of the process of
measurement and recording of test results. It should be pointed out that tests are
accompanied by changing the angular position of the profile model. Measurement results
are displayed in the form of dependences of aerodynamic coefficients on the angle of
attack.
    As the studied model of the aerodynamic surface, the compartment model of the
classical symmetrical NACA0012 profile was initially chosen. Such a model is often applied
in tests. This blowing model is characterized by a rectangular form. The model’s
characteristics are as follows: the span – 350 mm, the chord – 175 mm, the area – 0.06125
m2, and the geometric prolonging – λ=2.
    The longitudinal model’s axis is arranged at an interval of 25% of the chord length or
0.04375 m relative to the toe. It is observed that the angle location is determined relative
to the longitudinal axis. Vertical shields are mounted on the specific sides profile model.
The chosen sides must ensure a maximally effective expansion of the model less than 10.
Such a layout of the model makes it possible to single out the influence of the artificially
created structure of longitudinal vortices with the help of VVG on the profile resistance of
the profile model. First of all, it concerns pressure resistance, as well as friction, without
taking into account the inductive component from the influence of finite vortices.
    The model of the aerodynamic profile consists of VVG models. Volumetric vortex
generators are arranged in the toe of the profile model. which are installed in the toe part.
Using the technical information about the different types of VG, turbulators, aerodynamic
influxes, and other means of controlling vortex structures [11, 12, 13], the basic forms of
VVG were created. Three types of VVG in different projections are shown in Figure 1.
    The 1st type of VVG has a spindle-like elongated form. The juncture of the VVG and the
profile, as well as the side surface of the VVG itself, are causes of vortex generation in this
situation. The second type of VVG has a vertically flattened form with narrow front and
side edges similar to arcs. Both the above-described form and the juncture are sources
forming vortexes.
   The low VVG surface has cavitation assigned for increasing the vortex-generating and
carrying features of the VVG. The third type of VVG has a flattened form too. Moreover, it is
symmetrical about its median plane.
   The view from above looks like a triangular perimeter with a rounded toe and lateral
edge. The angular opening of the cone is 280. The source of generating longitudinal
vortexes is the juncture of the VVG with the profile and the lateral verge. The selected
types of VVG generate longitudinal vortices of different intensities and with different
spatiotemporal characteristics. In addition to the shape of the VVG, the factors influencing
the characteristics of the generation of the system of longitudinal vortexes are the size of
the VVG, their number, and the geometric characteristics of the setting-up.




Figure 1: VVG of three types for NACA0012 profile model.

    To change the sizes of the VVG mounted on the profile, it is essential to tighten the VVG
on the profile body or to scale its sizes. Firstly, the moving-out of the VVG from the profile
toe on the definite value is imitated. Secondly, the installation of the VVG of a certain size
is simulated.
    The main geometric characteristics of the composition of the VVG of different forms
and dimensions are the VVG area projecting beyond the leading edge of the profile, the
length of the ledge, the length of the leading edge of the profile model covered by the VVG,
the angle of mounting the VVG (an angle between the chords of the VVG and profile), the
step of mounting VVG by the span of the model. The specified geometric characteristics of
the VVG and their position in absolute values attributed to the characteristic dimensions of
the model can be useful during the analysis of the efficiency of vortex generation on the
aerodynamic characteristics of the profile model [14, 15].
    Tests of the profile model were carried out for three values of the average speed of the
airflow in the wind tunnel (V=9 m/s; V=18 m/s; V=27 m/s). They correspond to the
Reynold’s numbers Re ≈ 100000, 200000, and 300000, respectively.
   The airflow quality was characterized by an average degree of turbulence ε ≈ 3.5% at
V= 9 m/s, ε ≈ 2% at V= 18 m/s, and ε ≈ 1% at V= 27 m/s. The geometric angle of attack of
the model was varied in the range from -50 to 340, in some versions of the tests from 50 to
440 and covered large closed angles of attack.
   The real angle of attack was calculated taking into account the geometric and inductive
angularity of the airflow in the wind tunnel. In the protocols with the results of the
experiment, the real attack angles and the angles of attack of the α-mechanism are given.
The results of measurements of aerodynamic coefficients and airflow speed are calculated
as the average values of multiple measurements at a given angular position of the model.
The statistical characteristics of the measurements in the form of root mean square
deviations are given in [7, 8].

3. Experiment results
The analysis of the influence of the mounted VVG on the aerodynamic characteristics of
the NACA0012 profile model was implemented based on a comparative analysis of the
dependences of the aerodynamic coefficients as a function of the attack angle. These
dependencies were plotted on the same graph both for a pure profile and a profile with
mounted VG. The dependences of the aerodynamic coefficients describe the changing of
total loadings. Nevertheless, they are quite informative to show the differences in the
behavior of the dependencies. They also allow us to propose assumptions relative to the
change in the structure of the flowing upper surface of the model while achieving critical
attack angles and flow separation.
   In Figure 2 and Figure 3 the aerodynamic dependences for the layout of the first type
VVG are compared. Model in the amount of 5 pieces and the size of 20% are mounted with
a step of 70 mm. The graphs allow us to analyze the main specific effects of the mounted
VVG on the aerodynamic characteristics of the profile model.

                      1.0
                      Cya
                      0.8


                      0.6


                      0.4


                      0.2


                      0.0


                     -0.2


                     -0.4
                                                                                          0
                            -5   0          5       10           15          20       25  real 30
                                 - without VC:      - VC 1 p.;           - VC 3 p.;   - VC 5 p.;
                                                 - VC 7 p.;           - VC 9 p.


Figure 2: The influence of quantity and step of VG of the first type for the relative size of
the chord 7.5% and setting angle =0 on the plot 𝑐𝑦 = 𝑓(𝑎) of the profile model NACA0012
in conditions of airflow speed V=18.0 m/s (Reynolds number Re ≈ 200000).
              1.0
              Cya
              0.8


              0.6


              0.4


              0.2


              0.0


             -0.2


             -0.4
                                                                                     0
                    -5      0            5      10        15         20         25  real 30
                         - without VG;       - VG T2, 7,5%, 7 p.;
                                             - VG T2, 14,6%, 5 p.;        - VG T2, 20%, 5 p.

Figure 3: The influence of quantity and step of VG of the second type of the relative size of
the chord 7.5% and setting angle =0 on the plot 𝑐𝑦 = 𝑓(𝑎) of the profile model
NACA0012 under conditions of airflow speed V=18.0 m/s (Reynolds number Re ≈
200000).

   During the tests, the comparison was made for three values of the Reynolds numbers.
The graphs make it possible to make the main observations about the efficiency of the
mounted VVG on the separation properties, which are inherent to a greater or lesser
extent to each type of VVG.
   First of all, let us note the significant decrease in the influence of the Reynolds number
on the change in the aerodynamic characteristics of the model in the region around the
critical and at critical attack angles due to the installation of the VVG. The reason is
obvious. This is the additional turbulence of the flow on the surface of the model from the
turbulators. Installation of turbulators actually eliminates the characteristic change in
aerodynamic characteristics during flow stalling.
   This situation is characterized by a sharp decrease in lifting force after reaching a
critical angle of attack. In the presence of an air conditioner, the lift coefficient reaches a
certain constant level. This indicates the stability of the airflow structure, realized by the
imposition of a system of large longitudinal vortices on the airflow structure of a pure
profile at the critical angles of attack [16, 17].
   As a result of this behavior of the lifting force, the sign of the critical angle of attack
disappears. This situation corresponds to the coefficient 𝑐𝑦 𝑚𝑎𝑥 for low speeds (9 m/s). At
critical attack angles, the structure of the airflow becomes so stable that the lifting force
becomes less than in conditions of VVG absence. Nevertheless, with an increase in the
speed of the airflow, the change in the coefficient of lifting force occurs with the same
gradient as in the absence of VVG.
   It is known that a drastic decrease in the lifting force after reaching the critical angle of
attack is caused by a global separation of the airflow on the upper surface of the profile
model. Initially, the separation of the turbulent boundary layer develops gradually with an
increase in the attack angle. The separation line moves from the trailing edge of the profile
toward the flow. This process is reflected in the plot of the dependence 𝑐𝑦𝑎 = 𝑓(𝑎) in the
appearance of changing non-linearity. There is a jump-like transition of the separation line
to the toe of the profile at some point when the separation line approaches the line of
maximum profile thickness. Such a process is called global separation. It has a divergent
uncontrolled nature and spreads over the entire upper surface of the profile model. As a
result, there are sharp changes in the behavior of the aerodynamic characteristics
dependencies 𝑐𝑦𝑎 = 𝑓(𝑎), 𝑐𝑥𝑎 = 𝑓(𝑎) for the pure profile model. They are especially
noticeable in the behavior of the most sensitive characteristic 𝑚𝑧1 = 𝑓(𝑎).
   Under the conditions of the presence of VVG, there is a smooth change in the angle of
attack 𝑐𝑥𝑎 and longitudinal moment. Less resistance and higher aerodynamic quality at
critical values of attack angles in the conditions of the use of VVG are observed. Turbulent
separation in the presence of VVG begins at smaller angles of attack as follows from the
earlier appearance of a characteristic nonlinearity in the behavior of the dependence 𝑐𝑥𝑎 .
Most likely, this is caused by increased turbulence of the boundary layer by longitudinal
vortices from the VVG. Also, the separation of vortices from the airfoil surface occurred
earlier than the beginning of the separation of the turbulent boundary layer.
   Visualization of the structure of the airflow was implemented by silk threads pasted on
the profile model surface (see Figure 4).




Figure 4: Visualization of the airflow around the upper surface of the profile model at the
critical values of attack angles without VVG and with VVG.

   Taking into account the new data on the three-dimensionality of vortex structures on
the surface of the two-dimensional profile model under separation conditions [18] we can
make the following conclusions. The visualization results confirm the complexity of vortex
structures interacting with the boundary layer, especially in the conditions of the creation
of large longitudinal vortices from the VVG. They interact with large-scale vortices,
rotating in the plane of the wing. Thus, there is an interaction of two large vortices of
different types and their interaction with the turbulent boundary layer under separation
conditions. Without VVG, a pair of large vortices is formed in the plane of the surface near
the toe of the profile only after the global separation.
   When installing the VVG, a pair of vortices is also formed at smaller angles of attack in
the toe part of the profile. Nevertheless, the longitudinal vortex structure prevents global
separation. We have a pair of vortices in the plane of the surface on the sides of the VVG
and a longitudinal vortex in the track of the VVG. The longitudinal vortex restrains the
global separation by interacting with the boundary turbulent layer. This is implemented
by dividing into separate sections along the wing span by the number of installed VVG,
creating an ordered vortex structure of the flow and preventing it from breaking away
from the surface, introducing additional energy into the boundary layer.
   Due to the installation of vortex generators, is possible to predict an increase in the
critical value of the attack angle and the maximal lifting force and improve the flight-
technical characteristics. Such correction of the aerodynamic characteristics of the
carrying surface is important in the event of unpredictable near-critical and over-critical
angles of attack, regardless of the factor that caused the violation of the limits of the
permissible area of operation. These factors can include the influence of atmospheric
phenomena, piloting errors, or something else. On the other hand, the installation of
vortex generators should not impair the characteristics of the flight cruising mode. It is
desirable to keep values of the maximum aerodynamic quality and the most favorable
angle of attack. From this point of view, the criteria for the optimality of the layout of the
profile model with the volumetric VG can be the behavior of the lifting force on the near-
critical and over-critical attack angles. On the contrary, the absence of the VG can be
characterized by the values of the critical attack angle, the maximum value of the
coefficient of the lifting force, the gradient in the supercritical area, and the angle of attack
from which the violation of linear dependence begins. Important criteria will be the
magnitude of the change in aerodynamic quality in cruising flight modes due to the
installation of the VG, the magnitude of the most advantageous angle of attack, and
features in the change in the dependence of the longitudinal moment on the angle of attack
[19, 20, 21].
   The effect of the installation of the VG on the airflow around the bearing surface lies in
the artificial creation of an ordered system of large longitudinal vortexes by the span of
the model. It interacts with the external airflow and ensures a gradual and smooth local
restructuring of the airflow structure on the upper surface of the profile during
approaching to the over-critical attack angles. In contrast to the conditions of the absence
of vortex generators, this process is global and occurs explosively. The consequence of
installing vortex generators should be an increase in the critical angle of attack and the
absence of a drastic drop in the lifting force in the critical region. To increase the lifting
force, it is necessary to shape the upper surface of the toe part of the vortex generator in
such a way (for example, to set a certain curvature and increase the local speed) that a
suction force arises in this place. In the speed coordinate system, the projections of this
force increase the lifting force and decrease the drag [22, 23, 24].
   The effectiveness of the application of the VG is influenced by its size relative to the
dimensions of the profile, the number, and the step of their installation along the span of
the profile model, which together determine the surface area of the air conditioner served
by the air conditioner. The angle of installation of the vortex generator relative to the
chord of the airfoil affects the change of the aerodynamic characteristics of the airfoil
model by the angle of attack.
    As we can see, there are several geometrical parameters regarding the dimensions of
the VG of a certain type and their location, which affect the aerodynamic characteristics of
the profile models, and we have six profile models. There are also some criteria for the
efficiency of using vortex generators. Therefore, to find out the influence of the above
factors, a multivariate experiment with a large number of trials should be carried out.
    To shorten the test program, it was decided to determine the geometric characteristics
of the effective arrangement of the VG of 3 types of the developed shape for the model of
one type of aerodynamic profile and to carry out further studies of the layouts of other
profiles only for the determined effective arrangement of the selected profile. We define
such a layout as quasi-optimal of a local nature because in general, more efficient layouts
are also possible according to other criteria of optimality. The NACA0012 symmetrical
profile, well-known in experimental aerodynamics, was chosen as such profile. Layouts of
profile models with lateral screens are represented in Figure 5.




Figure 5: Layouts of profile model with lateral screens.

   Obtained results can be useful for synthesis control laws in disturbed stabilization and
control systems [25, 26, 27] and for reliable navigation measurements [28, 29].

4. Conclusions
The expected results of research in wind tunnels consist in improving the interaction of
VG with the features of the dynamics of the flow around the wing, the distribution of
velocity and the change in the pressure gradient along the wing, and the interaction of
eddy currents on the wing with terminal vortices. The results of experimental studies will
make it possible to set the problem of forming a class of profiles of vortex active wings.
   The most effective arrangement for each of the 3 types of the form of the VU turned out
to be their installation with a step of 70 mm (20% of the span of the model) in the amount
of 5 pieces with an angular position ‫=ﻻ‬0. The magnitude of the VU protrusion beyond the
front edge of the profile is 20% of the chord length of the profile.
   The graphs show the change in the aerodynamic force coefficients in the speed
coordinate system, the longitudinal moment coefficient relative to the point on the profile
chord at a distance of 25% from the toe, and the aerodynamic quality graph.

Acknowledgements
This research is partially supported by the Ministry of Education and Science of Ukraine
under the project “Grounding of aerodynamic layouts for perspective unmanned aerial
vehicles by volumetric vortex generators on aerodynamic surfaces” (# 0124U000221).

References
[1] E.P. Udartsev, Vortex Slat, 2014. Patent No. 94179, Issued Nov. 10th, 2014.
[2] E.P. Udartsev, The generator of longitudinal vortices, 2016. Patent No 111701. Issued
     Nov. 25th, 2016.
[3] J.S. Delinero, J.M. Leo, M.E. Camocardi, Vortex generators effect on low Reynolds
     number airfoils in turbulent flow, International Journal of Aerodynamics 2(1) (2020)
     1–14. doi: 10.1504/IJAD.2012.046539.
[4] P.S. Divekar, T. Ekbote, Design of aerodynamic of an airplane wings, IJSDR 4 (2019)
     59–64. URL: https://www.ijsdr.org/viewpaperforall.php?paper=IJSDR1910011.
[5] Y. Averyanova, V. Larin, N. Kuzmenko, I. Ostroumov, M. Zaliskyi, O. Solomentsev, O.
     Sushchenko, Y. Bezkorovainyi, Turbulence detection and classification algorithm
     using data from AWR, in: Proceedings of IEEE 2nd Ukrainian Microwave Week
     (UkrMW),         Kyiv,        Ukraine,       2022,       pp.      518–522.       doi:
     10.1109/UkrMW58013.2022.10037172.
[6] M. Zaliskyi, O. Solomentsev, V. Larin, Y. Averyanova, N. Kuzmenko, I. Ostroumov, O.
     Sushchenko, Y. Bezkorovainyi, Model building for diagnostic variables during aviation
     equipment maintenance, in: Proceedings of the 17th International Conference on
     Computer Sciences and Information Technologies (CSIT), Lviv, Ukraine, 2022, pp.
     160–164. doi: 10.1109/CSIT56902.2022.10000556.
[7] A. G. Shcherbonos, Aerodynamic characteristics of a wing with vortex generators in
     the conditions of non-stationary flow, Dissertation of Ph.D. 2012, NAU, Kyiv. (in
     Ukrainian).
[8] S. I. Oleksienko, Control of flow separation using longitudinal vortices generators,
     Dissertation of Ph.D., 2018, NAU, Kyiv. (in Ukrainian)
[9] N. Kuzmenko, I. Ostroumov, Y. Bezkorovainyi, Y. Averyanova, V. Larin, O. Sushchenko,
     M. Zaliskyi, O. Solomentsev, Airplane flight phase identification using maximum
     posterior probability method, in: Proceedings of the 3rd International Conference on
     System Analysis & Intelligent Computing (SAIC), Kyiv, Ukraine, 2022, pp. 1–5. doi:
     10.1109/SAIC57818.2022.9922913.
[10] O. Solomentsev, M. Zaliskyi, O. Sushchenko, Y. Bezkorovainyi, Y. Averyanova, I.
     Ostroumov, V. Larin, N. Kuzmenko, Data processing through the lifecycle of aviation
     radio equipment. in: Proceedings of the 17th International Conference on Computer
     Sciences and Information Technologies (CSIT), Lviv, Ukraine, 2022, pp. 146–151. doi:
     10.1109/CSIT56902.2022.10000844.
[11] S.M. Aftab, N.A. Razak, A.S. Rafie, K.A. Ahmad, Mimicking the humpback whale: An
     aerodynamic perspective, Progress in Aerospace Science 84 (2016) 48–69. doi:
     10.1016/j.paerosci.2016.03.002.
[12] J.C. Lin, Review of research on low-profile vortex generators to control boundary-
     layer separation, Progress in Aerospace Sciences 38 (2002) 389–420. doi:
     10.1016/50376-0421(02)00010-6.
[13] D. Custodio, C.W. Henoch, H. Johari, Aerodynamic characteristics of finite span wings
     with leading-edge protuberances, AIAA Journal 53(7) (2015) 1878–1893. doi:
     10.2514/1.J053568.
[14] A.K. Malipeddi, N. Mahmoudnejad, K.A. Hoffmann, Numerical analysis of effects of
     leading-edge protuberances on aircraft wing performance, Journal of Aircraft 49(5)
     (2012) 1336–1344. doi: 10.2514/1.C031670.
[15] O.A. Sushchenko, Y.N. Bezkorovainyi, Improvement of UAV positioning by information
     of inertial sensors, in: Proceedings of IEEE 5th International Conference on Methods
     and Systems of Navigation and Motion Control (MSNMC), Kyiv, Ukraine, 2018, pp.
     151–155. doi: 10.1109/MSNMC.2018.8576307.
[16] Z. Xingwei, Z. Chaoying, Z, Tao, J. Wenying, Numerical study on effect of leading-edge
     tubercles, Aircraft Engineering and Aerospace Technology 85(4) (2012) 247–257.
     doi:10.1108/AEAT-Feb.-2012-0027.
[17] O.A. Sushchenko, A. A. Tunik, Robust stabilization of UAV observation equipment, in:
     Proceedings of 2nd International Conference on Actual Problems of Unmanned Air
     Vehicles Development (APUAVD), Kyiv, Ukraine, 2013, pp. 176–180.
     doi: 10.1109/APUAVD.2013.6705318.
[18] R. Supreeth, A. Arokkiaswamy , K. Anirudh, R. K. Pradyumna, P. K. Pramod, A. K.
     Sanarahamat, Experimental and Numerical Investigation of the Influence of Leading
     Edge Tubercles on S823 Airfoil Behavior, Journal of Applied Fluid Mechanics 13(6)
     (2020) 188–189. doi: 10.47176/jafm.13.06.31244.
[19] R. Voliansky, O. Sadovoi, N. Volianska, Defining of Lyapunov functions for the
     generalized nonlinear object, in: Proceedings of IEEE 5th International Conference on
     Methods and Systems of Navigation and Motion Control (MSNMC), Kiyv,
    Ukraine, 2018, pp. 222–228, doi: 10.1109/MSNMC.2018.8576315.
[20] R. Voliansky, I. Ostroumov, O. Sushchenko, Y. Averyanova, O. Solomentsev, O.
     Holubnychyi, M. Zaliskyi, Variable-Structure Interval-Based Duffing Oscillator, in 42nd
     International Conference on Electronics and Nanotechnology (ELNANO), Kyiv,
     Ukraine, 2024, pp. 581–586.
[21] D. John D. Anderson, Jr., Christopher P. Cadou, Aerodynamics, McGraw Hill LLC, 2024.
[22] U. Gulcat Fundamentals of Modern Unsteady Aerodynamics, Cham, Springer, 2021.
     doi:10.1007/978-3-030-60777-7.
[23] D. Peery, Aircraft Structures, Dover Publications, Inc., 2014.
[24] L.S. Zhiteckii, V.N. Azarskov, K.Y. Solovchuk, O.A. Sushchenko, Discrete-time robust
     steady-state control of nonlinear multivariable systems: A unified approach, IFAC
     Proceedings           47          (3)         (2014)         8140–8145.            URL:
     https://www.sciencedirect.com/science/article/pii/S1474667016428971.
[25] O.A. Sushchenko, O.V. Shyrokyi, H2/H optimization of system for stabilization and
     control by-line-of-sight orientation of devices operated at UAV, in: Proceedings of 3rd
     International Conference on Actual Problems of Unmanned Aerial Vehicles, Kyiv,
     Ukraine, 2015, pp. 235–238.
[26] O.A. Sushchenko, Y.N. Bezkorovainyi, N.D. Novytska, Theoretical and experimental
     assessments of accuracy of nonorthogonal MEMS sensor arrays, Eastern-European
     Journal of Enterprise Technologies 3(9) (2018) 40–49.
[27] O.A. Sushchenko, Y.M. Bezkorovainyi, V.O. Golytsin, Processing of redundant
     information in airborne electronic systems by means of neural networks, in:
     Proceedings of IEEE 39th International Conference on Electronics and
     Nanotechnology, ELNANO-2019, Kyiv, Ukraine, 2019, pp. 652–655.
[28] O.A. Sushchenko, Y.N. Bezkorovainyi, N.D. Novytska, Nonorthogonal redundant
     configurations of inertial sensors, in: Proceedings of IEEE 4th International
     Conference on Actual Problems of Unmanned Aerial Vehicles Developments (APUAVD
     2017), Kyiv, Ukraine, 2017. doi: 0.1109/APUAVD.2017.8308780.
[29] K. Dergachov, O. Havrylenko, V. Pavlikov, S. Zhyla, E. Tserne, V. Volosyuk, et al., GPS
     usage analysis for angular orientation practical tasks solving, in: Proceedings of IEEE
     9th International Conference on Problems of Infocommunications, Science and
     Technology (PIC S&T), Kharkiv, Ukraine, 2022, pp. 187–192. doi:
     10.1109/PICST57299.2022.10238629.