Study on Flow Mechanism of a Morphing Supercritical Airfoil

In order to maintain the best performance in flight, a new concept, morphing aircraft, has been proposed, which can change the real-time aerodynamic characteristics under different flight conditions. +e key problem is to figure out the response of strong flow instability caused by structure changes during the morphing. To solve this problem, computational fluid dynamics (CFD) and wind tunnel tests (WTT) were employed. +e results show that the deformation of thickness and camber angle of the airfoil will significantly change the distribution of pressure and result in obvious hysteresis loops of lift and drag. With the increase of deformation frequency and amplitude, the instability increases correspondingly. Moreover, the unsteady effect caused by camber deformation is much stronger than that caused by thickness deformation. In addition, the flow structures on the airfoil, such as the shock strength and boundary separation location, have a delay in response to structure changes. +erefore, there will be a hysteresis between airfoil deformation and aerodynamic characteristics, which means strong flow instability.


Introduction
With the development of technology and science, more and more attentions have been paid to the quality and performance of aircrafts [1][2][3]. Nowadays, it is expected that changes in flow conditions have been adopted for aircraft configuration in order to maintain the optimal performance throughout the flight. erefore, morphing airfoil is developed gradually. However, the impact of configuration changes on the flow structure and aircraft performance is not yet understood, which may hinder the progress of the morphing airfoil. erefore, it is necessary to study the relationship and response law between the flow and configuration deformation.
Research institutes, such as NASA and DAPPA, have conducted long-term research on the related technologies of morphing vehicles [4][5][6], as well as a lot of research on aerodynamics. In the research of Ajaj et al. [7], the deformation technology was classified to put forward the concept of continuous deformation. Nekoubin and Nobari [8] calculated the deformation of the trailing edge of NACA0012 airfoil under the transonic condition and explained the influence of control parameters. Walker and Patil [9] derived the functions of unsteady lift, drag, and pitch moment and expressed them using the form of Chebyshev polynomials. Gandhi and Anusonti-Inthra [10] investigated the skin structure of variable camber wings. In the study of Adesen et al. [11], a dynamic stall model of variable trailing edge airfoil was derived based on the theory. As for domestic research, the development of morphing technology is later than foreign research. Xu et al. [12] proposed different types of morphing schemes for airfoils. Lv et al. [13,14] conducted wind tunnel tests to study the aerodynamic characteristics under transonic conditions. Hao and Yang [15] analyzed the influence of camber change on steady-state and unsteadystate aerodynamic characteristics under low Reynolds number. Chen et al. [16] discussed the influence of wing sweep changes on aerodynamic features, as well as its mechanism. At the same time, using analytical and discrete vortex methods proposed by Gao et al. [17,18], the unsteady aerodynamic characteristics of the morphing airfoil under subsonic/supersonic speed conditions were specifically studied. e research on morphing technology mainly focuses on the design of intelligence structure [19,20], but there are few studies on aerodynamic characteristics and flow mechanism of morphing airfoil. erefore, through CFD and wind tunnel tests, the influence of supercritical airfoil deformation on aerodynamic characteristics and flow structure has been discussed, and the results will be helpful for smart structure parameter choice and analysis on aerodynamic characteristics morphing airfoils. In this paper, effect of airfoil thickness and camber on its unsteady aerodynamic characteristics was studied by CFD simulation. When CFD was applied, the chord of RAE2822 was designated to 1 unit, based on which the parameters of the flow were obtained, such as the Reynolds number, and unsteady aerodynamic characteristics coefficients.

Research Topic
Only effect of airfoil camber changing through trailing edge deflection driven by a SMA smart structure on surface pressure distribution was measured in wind tunnel. Strictly speaking, this kind of deflection of wind tunnel model was not a real camber due to the difficulty in smart structure designing and wind tunnel experiments simulation. Experiment results could not be compared with those obtained by CFD, which was still helpful to recognize the flow phenomenon during model deformation in the wind tunnel environment. e model was extracted to a 2D wing based on RAE2822 airfoil, wing span and chord length of which were 365 mm and 150 mm, respectively ( Figure 1).

Parameters for CFD Simulations.
Aerodynamic characteristics of wing were decided by its airfoil parameters, of which airfoil camber and thickness were the most important parameters. Unsteady effect of morphing on aerodynamic characteristics over supercritical airfoil was mainly studied by CFD.
Flow structure and aerodynamic forces under different airfoil camber and thickness were simulated with Mach number ranging from 0.7 to 0.78, Reynolds number ranging from 5 × 10 6 to 30 × 10 6 , and angles of attack ranging from 0°t o 8°. e variation of real-time airfoil thickness ranged from -15 percent to 20 percent of the basis airfoil thickness, and the angle which was used to measure airfoil camber between the chord trailing edge tangent and the horizontal axis ranged from −5°to 5°. Definition of airfoil thickness and camber would be introduced in Section 2.1.3. Camber and thickness morphing frequency ranged from 0 Hz to 2 Hz when the effect was simulated with CFD.

Contents in Wind Tunnel
Experiments. When experiments were conducted in the wind tunnel, Mach number ranged from 0.4 to 0.8, test Reynolds number was kept at 2 × 10 6 approximately, and angles of attack ranged from 0°to 6°.
SMA material phase transition would result in the problem of slow response of smart structure, which made flow over supercritical airfoil kept in steady state during the trailing edge deflection. Wind tunnel experiments would be used only to recognize the steady flow structure evolution when trailing edge was deflected. Because of the effect of aerodynamic loads brought about by the dynamic pressure of the wind tunnel flow, the trailing edge deflection magnitude could not exceed 10.9°. It is not reasonable to compare wind tunnel experiment results with those obtained by CFD due to the different morphing mode and flow mode.

Computational Fluid
e convective item was discretized with central scheme and the viscous item was discretized with ROE scheme. e time discretization of the unsteady-state calculation is carried out on the basis of the dual time step. When the maximum error of the basis variable quantity was less than 0.00001, the iteration was considered to be converged. At the same time, endless loop was avoided by designating inner time step. In addition, ideal-gas far field conditions were used as boundary conditions. e wall was adiabatic and nonslip. S-A turbulence model was applied.

Grid and Validation
. C grid containing 941 × 144 grid nodes is used ( Figure 2). In order to avoid the negative volume caused by chaotic deformation of large aspect ratio, the mesh field is divided into two parts, and only grid deformation is allowed. Since Reynolds number effect is not our main purpose, only characteristics are considered under full turbulence. y+ and the height of the first layer grid are 2 and 0.00002664, respectively. e effectiveness of the steady method is verified according to Case 9 in the research of Zhang and Zhang [21], and the parameters are shown in Table 1. e comparison between the results of Case 9 and present study is shown in Table 2. e results are consistent with Case 9, showing that the steady method applied for the study was effective. e effectiveness of the unsteady method is verified according to AGARD CT1 case [21], and the parameters are shown in Table 3. e comparison of unsteady method with CT1 case is shown in Figure 3. It can be seen that dual time step discretization was qualified for unsteady computation; and the consistency of both results was well.

Timestep Determination.
Lift coefficient under different timestep is shown in Figure 4. As can be seen, results of C L coefficient were converged when real-time step size was larger than 1000. erefore, inner iteration step size and the real time step size were assigned as 50 and 1000, respectively.

Definition of Airfoil ickness and Camber Changes in CFD Simulation.
e change of airfoil thickness is obtained by increasing the coordinates of all the original airfoil control points on average. en, the coordinate equation of the control points is recommended: where Y init is the initial coordinate of the airfoil control point, f is the deformation frequency of the airfoil thickness, and A is the relative amplitude of the thickness change, as shown in Figure 5.
It is difficult to express airfoil centerline equation algebraically. So, the airfoil chord was defined as a parabola to describe the camber. In order to keep the thickness constant during the camber change, the equation of camber changes is given as follows: where P is the vertex instantaneous coordinate and Δy is the increase of the longitudinal coordinates of the surface control points. For convenience, the angle (θ) between the tangent of deformed chord at the trailing vertex and abscissa was used to measure chord camber changing. e camber change is shown in Figure 6.

Wind Tunnel Test
e test was conducted in a transonic wind with a Mach number ranging from 0.4 to 3.5. It is a semireturn wind transient wind tunnel with a cross-sectional dimension of 0.6 × 0.6 m. e upper and lower walls of this section are slotted, and the side walls are solid with rectangular optical windows and suitable for camera shooting.

Smart Actuator Structures.
To understand the flow mechanism on morphing airfoils, a typical 2D supercritical airfoil test model was designed. e cantilever flexible flexure structure was applied to drive the trailing edge flexure. e structure is driven by a shape memory alloy (SMA) actuator, which can deflect the trailing edge by 10.9°. e cavity on the lower surface of the model was filled with glass glue to maintain the smoothness of the lower surface. In Figure 7, the SMA actuator model is recommended.

Test Technology.
Since the lower surface was filled with glass glue, it is difficult to fix the pressure orifice. erefore,    Figure 1: Sketch of the model.

Shock and Vibration 3
PSP technology was used to obtain the model surface pressure distribution under the condition of boundary layer natural transition. Image acquisition was composed of light source, camera, and control software. e excitation light source was a kind of LED cold light source with 12 W output power. e camera type was PCO.1600. In order to fix the deflection angle of the trailing edge, deformation measurement technology (VMD) was applied. In Figure 8, the position of the model in the test section is presented.

Results and Analysis
3.1. Influence of ickness Deformation on Aerodynamic Characteristics. Figure 9 shows the influence of three thickness amplitudes (10%, 15%, and 20%) on the aerodynamic characteristics. It suggests that, with the increase of the deformation amplitude, the hysteresis phenomenon becomes more obvious, which is caused by the pressure and vorticity difference at the same position during the thickness deformation process (Figure 10). e area of the hysteresis loop formed by the lift coefficient can be understood as the work required by the external force during the deformation of airfoil. e greater the area, the greater the power required [22]. e influence of thickness deformation frequency is shown in Figure 11. e higher the deformation frequency, the greater the lift/drag force generated by the rate of change of the velocity potential over time. Within a certain range of frequency and amplitude, the difference between lift and drag is close to 0.015 and 0.0019, respectively, when the trailing edge passes through the same position of reciprocating motion.

Influence of Camber Deformation on Aerodynamic
Characteristics. Figure 12 shows the influence of the extent of camber deformation on aerodynamics. It indicates that the larger the deformation, the stronger the unsteady-state effect and the larger the area of the hysteresis loop. Figure 13 indicates that the higher the deformation frequency, the larger the area of the hysteresis loop, and the additional unsteady lift caused by the ratio of the velocity potential over time will also be greater.
When the angle of attack was increased to 6°, the same effect would be found as α � 2°. But, as can be seen in Figure 14, camber deformation would have less effect on aerodynamic characteristics at α � 6°than at α � 2°; and the hysteresis characteristics of lift/drag force curves were not as regular at α � 6°as at α � 2°. is resulted from the obvious shock wave and its induced boundary layer separation. As can be seen in Figure 15, strong separation happened, and position of flow separation moved forward when camber changing magnification angle θ was −3°again. We could obtain that flow structures were very different at the same position but different time. at was the unsteady flow that led to the hysteresis characteristics of lift/drag forces during morphing. e above results show that the airfoil thickness and camber deformation will cause obvious instability, and the aerodynamic coefficient has obvious hysteresis characteristics. Moreover, the greater the amplitude of the airfoil parameter deformation and the higher the frequency, the stronger the flow instability. By comparing it with thickness deformation, the change of the camber angle will cause more severe instability, which will lead to greater additional instability. For example, the difference between lift and drag is close to 0.117 and 0.0034, respectively, which indicates that the instability of vehicle deformation should not be ignored.

Mechanism of Unsteady Flow during Airfoil Deformation.
e effect of deformation on the surface pressure distribution is shown in Figure 12, which indicates that when the trailing edge returns to pass through the same position, the shock wave moves to the airfoil nose and its strength is weakened. is is caused by the hysteresis and strength of the shock wave position relative to geometric deformation.
is phenomenon can also be observed in the PSP wind tunnel test results (Figure 16). It can be seen that, with the deflection of the trailing edge, the shock wave has been strengthened and pushed back. e flow near the center of   Shock and Vibration 5 rotation accelerates again and shows a tendency to form a second shock wave. e influence of deformation on the surface pressure distribution is shown in Figure 12, which indicates that when the trailing edge returns to pass through the same position, the shock wave moves to the airfoil nose and its strength is weakened. is is caused by the hysteresis and strength of the shock wave position relative to geometry deformation.
is phenomenon can also be observed from PSP wind tunnel test results ( Figure 17). It can be seen that, with the deflection of the trailing edge, the shock wave has been strengthened and pushed back. e flow near the center of rotation center accelerates again and shows a tendency to form the second shock wave. Figure 18 shows how the amplitude and frequency of airfoil camber deformation affect the surface pressure distribution. e graph shows that larger amplitude and higher frequency will cause larger changes in pressure distribution.  Shock and Vibration For instance, such change can significantly enhance the shock wave and move the shock wave to the end of the airfoil, resulting in the increased unstable additional lift. e results show that the unsteady flow effect in the wing morphing process is caused by flow separation, boundary layer transition, shock wave oscillation, and so forth. In the research of Yang [23], the author believes that the aerodynamics of an unsteady morphing airfoil can be represented by a constant item and a stable item. e additional unstable item is shown as In the process of airfoil deformation, the rate of velocity potential change will cause virtual mass force, which may lead to unstable additional force. Because the velocity change rate of the airfoil parameters is different, the unstable additional force/moment is also different, thus forming a hysteresis loop.

Conclusions
In this study, the influence of the thickness and camber deformation of the supercritical airfoil on flow instability and its mechanism was studied.
e conclusions can be drawn as follows: (1) Flow instability will occur with the change of airfoil thickness and camber angle, and it exhibits obvious hysteresis on the aerodynamics forces of the hysteresis loop. erefore, in the design of morphing aircrafts, full attention must be paid to the unsteady aerodynamic characteristics of deformation. (2) Larger morphing amplitude and higher morphing frequency will lead to stronger flow hysteresis and unstable effects. Compared with the influence of airfoil thickness deformation, the instability caused by airfoil camber deformation is more significant, which cannot be ignored.   Shock and Vibration