Unsteady Subsonic Aerodynamic Characteristics of Wing in Fold Motion

International Journal of Aeronautical and Space Sciences.
2011.
Mar,
12(1):
63-68

- Received : November 23, 2010
- Accepted : February 21, 2011
- Published : March 30, 2011

Download

PDF

e-PUB

PubReader

PPT

Export by style

Article

Metrics

Cited by

TagCloud

Aerodynamic characteristics of a wing during fold motion were investigated in order to understand how variations or changes in such characteristics increase aircraft performance. Numerical simulations were conducted, and the results were obtained using the unsteady vortex lattice method to estimate the lift, drag and the moment coefficient in subsonic flow during fold motion. Parameters such as the fold angle and the fold angular velocity were summarized in detail. Generally, the lift and pitching moment coefficients decreased as the angle increased. In contrast, the coefficients increased as the angular velocity increased.
Schematic configuration of the wing in fold motion.
In the figure, both φ
_{1}
and φ
_{2}
may be simultaneously varied during flight for efficient flight performance. These parameters affect aerodynamic characteristics such as lift and drag. The relationships between the global coordinates and local coordinates of each wing component are expressed as
where L
_{1}
, L
_{2}
are the span length of the inboard and the outboard component, respectively, and Λ is the sweptback angle.
and the first boundary condition requiring zero normal velocity across the solid boundaries is
where v is the velocity of the surface and n is the vector normal to the moving surfaces. Also the second boundary condition is diminishing the flow disturbance far from the body.
where R=(X, Y, Z).
For unsteady flow, the Kelvin conditions supply an additional equation in determining the streamwise strength of the vorticity shed into the wake. In general, the condition states that in the potential flow region the angular momentum cannot change, and thus the circulation Г around a fluid curve enclosing the wing and its wake is conserved:
From the continuity equation, the general solution can be obtained by integrating the contribution of the basic solution of source σ and doublet μ distributions over the surface and wakes of the body,
Although, this formulation does not directly include a vortex distribution, doublet distributions can be replaced by equivalent vortex distributions.
On the other hand, the wake shed from the trailing edge of the lifting surfaces can be modeled by vortex distributions. The strength of each vortex Г
_{Wi}
is equal to the vorticity shed during the corresponding time step Δt such that
Consequently, the strength and location of the vortex must be specified for each vortex element. The wake vortex location should be aligned with the trailing edge and be located closer to the latest position of the trailing edge.
The zero normal velocity boundary condition on the surface (Q
_{nk}
=0) for an arbitrary collocation point k can be obtained the following equations.
Once the computations of the influence coefficients and the right-hand side vector are completed, the zero normal flow boundary condition on all collocation points of the wing will result as
And then the circulation can be obtained from this equation. Furthermore, the aerodynamic forces can be computed by using the Bernoulli equation and the pressure difference can be as follows.
where Δc
_{i, j}
and Δb
_{i, j}
are the panel lengths in the i-th and j-th directions, respectively. Similarly, τ
_{i}
and τ
_{j}
are the panel tangential vectors in i and j directions.
Then the load of this panel is
Therefore, the total forces and moments are obtained by the load of each panel. (Editor’s Note: This section is very nicely written)
_{2}
was set to zero.
The location of the wake behind the trailing edge was varied according to fold angle. Thus, if the wing is in motion,the wake will follow the trailing edge, and then the effects of the wake will not be same for a constant angle. The location of the wake in motion is shown in
Fig. 2
.
Additionally, the collocation points of each aerodynamic panel were changed during fold motion and then the circulation in Eq. (9) will be changed. As a result, the varied circulation value of the panels will produce differing aerodynamic characteristics. Therefore, the characteristics
The location of the wake behind the trailing edge in the motion.
Model configuration (unit: inch)
of the wing were varied for the wing in fold motion.
Table 1
represents the model details (Lee and Chen, 2006). This model is an innovative aircraft design proposed by Lockheed Martin as part of the DARPA MAS (Lee and Chen, 2006).
Comparison of the transient lift coefficient variation for U_{∞} Δt/c=1/16 and a=5°.
^{3}
and the angle of attack was 5°. Also, in order to reduce the effect of sudden acceleration, the folding was initiated at 2 seconds. As shown in
Fig. 4
, the total lift
Transient lift coefficient and pitching moment coefficient in fold motion.
Sectional lift coefficient variation in fold motion.
coefficient and pitching moment coefficient decreased.
If the aerodynamic forces orthogonally acting on the individual panels were taken into consideration, the tendency of the transient lift coefficient could be slightly different from that given in
Fig. 4
. But, the z-direction forces were only considered in the analysis, which was the direction of the total lift.
In
Fig. 5
, the sectional lift coefficient is shown. The coefficient of the inboard component largely decreased.
The distribution of the sectional lift changed as the fold angle reached 90°. As shown in
Fig. 5
, aerodynamic characteristics of the wing in fold motion varied. Additionally, the decrease of the lift coefficient of the inboard component was observed as the fold angle increased.
The pressure coefficient vs. nondimensional chord-wise direction of each part is represented in
Fig. 6
.
Large variations did not occur except for in the inboard component, similar as
Fig. 5
.
Pressure coefficient vs. nondimensional chordwise direction.
The fold angular velocity,ø was another parameter taken into consideration. The large fold angular velocity produced
Aerodynamic characteristics variation.
opposite effects in the aerodynamic characteristics in comparison to the small angular velocity.
In
Fig. 7
, the variations of transient lift coefficient, drag coefficient and pitching moment coefficient are represented for each angular velocity.
From these results, a larger angular velocity gained more lift, drag and also pitching moment. Moreover, the lift coefficient decreased and the drag coefficient increased as the fold angle increases except case of 120 rad/sec.

1. Introduction

Shape morphing is defined as changing the configuration of a wing in order to control the aerodynamic characteristics so that the wing may comply with the aircraft’s mission profile. The morphing wing has been investigated by several groups such as the Active Aeroelastic Wing Program of NASA and the Morphing Aircraft Structures Program (MAS) provided by Defense Advanced Research Projects Agency (DARPA). Smart materials have also been researched in conjunction with the morphing wing investigations.
According to the research conducted by Rodriguez (2007), four approaches for morphing wing technology may improve aircraft performance: expand the wing’s flight envelope, replace conventional control surfaces for flight control for performance and stealth advancement, reduce drag to improve range and reduce vibration, or control flutter. The benefits and future concepts of a fully morphing aircraft configuration have been investigated (Moorhouse et al., 2006). Furthermore, several concepts that focus on morphing skins were studied, such as corrugated structures, segmented structures, reinforced elastomers or flexible matrix composite tubes embedded in a low modulus membrane (Thill et al., 2008). In addition, wind tunnel tests were performed testing the morphing wing. Unsteady pressure signals were recorded and visualized in real time while the wing was being deformed to reproduce various airfoil shapes with controlling the two actuators displacements (Popov et al., 2010a, b).
In Sofla et al. (2010), morphing wing concepts were classified into three groups: planform alternation, airfoil adjustment and out-of-plane transformation. With planform alternation, wing area manipulation techniques are employed through which the wing’s span and chord length are resized. Alternatively, airfoil adjustment entails designing and changing the wing profile. Out-of-plane transformation consists of altering the chordwise and spanwise camber. The folding wing is capable of adjusting the flight performance from a cruise configuration to a high speed dash configuration. During the intermediate changing configuration, aerodynamic characteristics varied dramatically and quickly.
Finite element modeling, aerodynamic loading, and a continuous-time multi-body simulation were used to simulate the aeroelasticity of the wing of an unmanned combat aerial vehicle (Scarlett et al., 2006; Snyder et al., 2009). A methodology for the model and an analytical solution were developed in order to produce flutter solutions for a continuous representation of a two-segment uniform folding wing (Liska and Dowell, 2009). Moreover, a flutter analysis simulation was developed and tested while accounting for the nonlinearity of a wing structure (Attar et al., 2010; Tang and Dowell, 2008).
The unsteady vortex lattice method (UVLM) was developed to predict the geometry of the wake and the aerodynamic loads on the wing configuration in an incompressible flow. UVLM provides a highly effective solution of a complex problem with reasonable computational efficiency. The side-edge suction force for isolated or interaction planform was computed by applying the method (Lamar et al., 1975). Also it was applied to aircraft engine designs to improve the overall performance of both military and civil aircraft. UVLM was coupled with a typical section cascade structural model to form a set of aeroelastic equations of motion (Thake, 2009).
Most studies concerning the folding wing were performed under the assumption that the fold angle remained constant. However, this assumption produces varying results for a wing that is in motion. Accounting for this motion is a necessary action when designing the folding wing structure. In this paper, an aerodynamic analysis of a wing in the motion was performed by using UVLM.
2. Formulations

- 2.1 Folding wing

Figure 1
shows the schematic configuration of a wing at an instant of folding motion.
Lager Image

Lager Image

- 2.2 Unsteady vortex lattice method

The fluid surrounding the body was assumed to be inviscid, irrotational, and incompressible over the entire flowfield, excluding the solid boundaries and wakes of the body. According to Katz and Plotkin (2001), a velocity potential Ф(X, Y, Z) can be defined in the inertial frame; the continuity equation in this frame becomes
Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

3. Numerical Simulation

Fold angle and angular velocity were the most important parameters in the numerical simulation. These parameters potentially change the aerodynamic characteristics of a wing while the wing is in motion. The simulation and verification were performed by C++ and MATLAB.
- 3.1 Model configuration

The wing consisted of three parts: the inboard and outboard components which undergo rotating motion and the body component that is not subject to motion. For this simulation, only the inboard component was assumed to rotate about the x-axis from 0° to 85° ø
Lager Image

Model configuration (unit: inch)

Lager Image

- 3.2 Verification

For the verification of this work, an aerodynamic analysis was performed using a simple uncambered, rectangular wing model that was suddenly set into a constant-speed forward flight, as given by Katz and Plotkin (2001). The results were almost equal to the data of the reference as shown in
Fig. 3
.
Lager Image

4. Results and Discussion

The freestream velocity was 50 m/s, the air density was assumed as 1 kg/m
Lager Image

Lager Image

Lager Image

Lager Image

5. Conclusions

An aerodynamic model for a folding wing was developed and vortex lattice method was used for the analysis. The model was derived from the DARPA MAS and analyzed for various fold angles and fold angular velocities.
The location of the wake behind the trailing edge varied according to the fold angle and the fold angular velocity. Thus, aerodynamic characteristics such as lift, drag and moment of the wing in fold motion also varied during flights.Also as the fold angle increased, the lift coefficient and the drag coefficient decreased. This indicates that the folded wing configuration was more suitable to dash. However, the results of this paper showed that the effect of a large fold angular velocity was also an important factor that cannot be ignored.
According to the results obtained in this study, the various aerodynamic characteristics in fold motion affect aeroelastic characteristics. Therefore, flutter analysis of a folding wing during the motion will be investigated.
Attar P. J
,
Tang D
,
Dowell E. H
2010
Nonlinear aeroelastic study for folding wing structures.
AIAA Journal
48
2187 -
2195
** DOI : 10.2514/1.44868**

Katz J
,
Plotkin A
2001
Low Speed Aerodynamics
2nd ed
Cambridge University Press
Cambridge

Lamar J. E
,
Gloss B. B
Langley Research Center
1975
Subsonic Aerodynamic Characteristics of Interacting Lifting Surfaces with Separated Flow Around Sharp Edges Predicted by a Vortex-Lattice Method.
National Aeronautics and Space Administration
Washington DC

Lee D. H
,
Chen P. C
2006
Nonlinear aeroelastic studies on a folding wing configuration with free-play hingenonlinearity
47th AIAA/ASME/ASCE/AHS/ ASC StructuresStructural Dynamics and Materials Conference
Newport RI
1754 -
1775

Liska S
,
Dowell E
2009
Continuum aeroelastic model for a folding-wing configuration.
AIAA Journal
47
2350 -
2358
** DOI : 10.2514/1.40475**

Moorhouse D
,
Sanders B
,
Von Spakovsky M
,
Butt J
2006
Benefits and design challenges of adaptive structures for morphing aircraft.
Aeronautical Journal
110
157 -
162

Popov A. V
,
Grigorie L. T
,
Botez R
,
Mamou M
,
Mebarki Y
2010a
Closed-loop control validation of a morphing wing using wind tunnel tests.
Journal of Aircraft
47
1309 -
1317
** DOI : 10.2514/1.47281**

Popov A. V
,
Grigorie L. T
,
Botez R
,
Mamou M
,
Mebarki Y
2010b
Real time morphing wing optimization validation using wind-tunnel tests.
Journal of Aircraft
47
1346 -
1355
** DOI : 10.2514/1.47431**

Rodriguez A. R
2007
Morphing aircraft technologysurvey
45th AIAA Aerospace Sciences Meeting
Reno NV
15064 -
15079

Scarlett J. N
,
Canfield R. A
,
Sanders B
2006
Multibody dynamic aeroelastic simulation of a foldingwing aircraft
47th AIAA/ASME/ASCE/AHS/ASC StructuresStructural Dynamics and Materials Conference
Newport RI
6398 -
6407

Snyder M. P
,
Sanders B
,
Eastep F. E
,
Frank G. J
2009
Vibration and flutter characteristics of a folding wing.
Journal of Aircraft
46
791 -
799
** DOI : 10.2514/1.34685**

Sofla A. Y. N
,
Meguid S. A
,
Tan K. T
,
Yeo W. K
2010
Shape morphing of aircraft wing: Status and challenges.
Materials and Design
31
1284 -
1292
** DOI : 10.1016/j.matdes.2009.09.011**

Tang D
,
Dowell E. H
2008
Theoretical and experimental aeroelastic study for folding wing structures.
Journal of Aircraft
45
1136 -
1147
** DOI : 10.2514/1.32754**

Thake M. P
2009
Effect of Mistuning on Bending-Torsion Flutter Using a Compressible Time-Domain Aerodynamic Theory.

Thill C
,
Etches J
,
Bond I
,
Potter K
,
Weaver P
2008
Morphing skins.
Aeronautical Journal
112
117 -
139

Citing 'Unsteady Subsonic Aerodynamic Characteristics of Wing in Fold Motion
'

@article{ HGJHC0_2011_v12n1_63}
,title={Unsteady Subsonic Aerodynamic Characteristics of Wing in Fold Motion}
,volume={1}
, url={http://dx.doi.org/10.5139/IJASS.2011.12.1.63}, DOI={10.5139/IJASS.2011.12.1.63}
, number= {1}
, journal={International Journal of Aeronautical and Space Sciences}
, publisher={The Korean Society for Aeronautical & Space Sciences}
, author={Jung, Yoo-Yeon
and
Kim, Ji-Hwan}
, year={2011}
, month={Mar}