Aerodynamic Reduced-Order Modeling without Static Correction Requirement Based on Body Vortices

The objective of this research is to propose a new reduced-order modeling method. This approach is based on fluid eigenmodes and body vortices without using static correction. The vortex lattice method (VLM) is used to analyze unsteady flows over twodimensional airfoil and three-dimensional wing. Eigenanalysis and reduced-order modeling are performed using a conventional method with static correction and an unconventional one without the static correction. Numerical examples are proposed to demonstrate the performance of the present method. The results show that the new method can be considered an alternative way to perform the reduced-order models of unsteady flow.


Introduction
Reduced-order modeling (ROM) with a static correction technique is an efficient technique which has been used for analyzing unsteady incompressible aerodynamic flows [1].A few unsteady flow eigenmodes are used to construct the reduced-order unsteady flow models similar to the modal analysis of structures.Such a technique was used to study an eigenanalysis of unsteady potential flows about flat airfoils, cascades, and wings.Later, this technique has been used to analyze a wide variety of unsteady flow models such as a compressor cascade [2] and complex configurations [3].Unfortunately, conventional ROM without static correction technique cannot generate satisfactory results even with all of eigenmodes [4].However, the ROM with static correction cannot give a satisfactory result if its time step computed a quasisteady part of the solution.Behbahani-Nejad has shown the eigenmodes that are equal to the body's elements behave exactly in quasistatic fashion [4].He proposed to construct a ROM only based on the wake eigenmodes, which means that the body quasistatic eigenmodes are removed.Based on such conditions, acceptable results were obtained without using the static correction technique.The new technique is used to solve many fluid flow problems.Shahverdi et al. [5] presented that the new technique can construct an efficient ROM based on the boundary element method to analyze unsteady flows over an airfoil, wing, and wing body.Later, Shahverdi et al. [6] proposed to study an aeroelastic analysis for helicopter rotor blade in hover using the same technique.Dardel and Bakhtiari-Nejad [7] applied the same technique to find a reduced order of complete aeroelastic model for limit cycle oscillations.Their model has significantly low dimension in comparison with other previously studied models.
Another method for constructing ROMs is proper orthogonal decomposition (POD), which has been used for investigating unsteady aerodynamic and aeroelastic models [8].Kim [9] has used the method for frequency domain of a vortex lattice fluid model.From the literature review of the POD, it has been used to obtain low-dimensional models of fluids [10], but most of the applications have been used for incompressible flows.Later, Rowley [11] proposed a framework for applying the POD and Galerkin projection for compressible fluids.In compressible flow, both the kinematic and thermodynamic variables are dynamically important.They proposed an energy-based inner product used to obtain POD modes.The result was that the POD can be extended to the compressible flow.
The purpose of this research is to derive an alternative ROM based on a vortex lattice method (VLM) in which the eigenvalue problem is defined only based on the unknown body vortices.The new ROM does not need the static correction as with the previous formulation, which is based on the wake vortices.The eigenvalues of the new eigensystem are nonzero or have no quasistatic eigenmodes.Eigenanalysis results show that the eigenvalues of the proposed method are equal to the corresponding nonzero values of the previous studies.Furthermore, the reduced-order models are constructed for steady flows over a two-dimensional airfoil and a three-dimensional wing.The results show that the present methods can accurately analyze unsteady flows and it is an alternative reduced-order model.

Eigenanalysis and ROM
The vortex lattice method for unsteady flow (termed as unsteady vortex lattice method (UVLM)) is a relation of downwash and vortex of unsteady flow [1]: where Γ is a vector of vortex strength and w is a vector of downwash.The matrices A and B are large sparse matrices.
The Γ +1 at time +1 can be calculated in terms of the solution at time  by An unsteady lift can be calculated as follows [12]: where () = Γ/ and  is semichord of the airfoil.For zero downwash, where Γ =   exp    and   = exp   Δ, the generalized eigenvalue problem can be obtained as where   is the th eigenvalue and   is the corresponding right eigenvector.More generally, (4) can be written in the form of right eigenvectors  as follows: where Z is a diagonal matrix containing the eigenvalue and X is a matrix in which columns are the right eigenvectors.Similarly, the eigenvalue problem in the form of left eigenvectors Y can be written as where Y is a matrix, with the rows being the left eigenvectors.Similar fashion to structural modal analysis, the eigenvectors satisfy the orthogonality conditions The dynamic behavior of the fluid as the sum of the individual eigenmodes can be expressed as where c is the vector of normal mode coordinates.Substitution of the ( 8) into (1), premultiplying by Y  , and using the orthogonality condition give a set of  uncoupled equations for the modal coordinates c, From ( 9), each mode can be marched forward in time independently.The solutions can then be reassembled using (8) to obtain the vortex strength.
The reduced-order model can be carried out by retaining only a few  of the original modes that give the largest eigenvalues ( ≪ ).However, the preceding reduced-order model does not produce satisfactory results due to neglected modes that are not orthogonal to the downwash.The way to solve this problem is to decompose the unsteady solution into two parts.The first part is equivalent to the response of the system if the disturbance is quasisteady and the second part is dynamic part.The unsteady solution is shown as follows: The quasisteady portion Γ  is given by Thus, ( 9) is replaced by

ROM without Static Correction Requirement Based on Wake Vortices
Previous work [5] proposed the method to remove zero eigenvalues by rearranging (1) in such a way that its eigenvalues are the same as the nonzero eigenvalues of the previous equation.The way to solve the problem is to decompose the vortex strength into two parts as body (Γ  ) and wake (Γ  ) vortex strengths as follows: Substitution (13) into (1) and partitioning it lead to From the previous equation, it can be arranged as  It can be shown that the matrices B 11 and B 12 are zero; therefore, (15) can be rearranged to a new equation as Substitution of the previous equation into (16) gives  or where From ( 19), which is based on only Γ  , the corresponding eigensystem has no zero eigenvalues.Therefore, one may construct accurate reduced-order models without using a static correction technique.or where From ( 23), which is based on only Γ  , the corresponding eigensystem has no zero eigenvalue, which is similar to the previous method.Therefore, we construct the accurate reduced-order models without using static correction technique.Note that since to matrix A 12 can be a nonsquare matrix, its inverse cannot be found.This can be solved by using the same numbers of elements for wing panels and wakes.

Results and Discussion
5.1.Our Vortex Lattice Model.In this section, we validate unsteady vortex lattice model with two unsteady flow models.The first model is the lift acting on the airfoil due to a step change in airfoil downwash (Wagner problem).The result is shown in Figure 1, and our vortex lattice model results are in agreement with the Wagner function.Figure 1 shows the indicial response of a rectangular wing due to rigid body plunging motion in comparison with the corresponding results presented in [1,4].The wing model has aspect ratio of  5.0 and it is modeled with 8 vortex elements in the streamwise direction and 10 element panels in the spanwise direction, while the wake is modeled with 40 vortex elements in streamwise direction.The vortex relaxation factor of this model is 0.992.Also shown in Figure 1, the results of our threedimensional vortex lattice model are in good agreement with [1,4].For the second model, the unsteady lift due to plunging motion of an isolated airfoil is computed a range of reduced frequencies and compared with the Theodorsen exact solution as shown in Figure 2. The airfoil is modeled using 20 vortex elements, while the wake is model with 200 vortex elements.The vortex relaxation factor is 0.996.The results in Figure 2 show our model is quite satisfactory with the exact solution.

Eigenanalysis.
This section proposed the eigenanalysis of the conventional method, the method based on Γ  , and the present method (based on Γ  ).The eigenvalues of the twodimensional airfoil are shown in Figure 3.Both airfoil and wake are modeled using 10 vortex elements.The results show that the eigenvalues of the proposed method are the same as the nonzero eigenvalues of the conventional method and the method based on Γ  .In the proposed eigenanalysis, the eigensystem is interpreted only by the wake elements [4].
Eigenvalues of vortex lattice model of unsteady flow about the two dimensional airfoil and wake are modeled with 100 vortex elements are shown in Figure 4.
Eigenvalues of vortex lattice model of unsteady flow about the three-dimensional wing are plotted in Figure 5. From the figure, the nonzero eigenvalues of the conventional method and the method based on Γ  are the same with the eigenvalues of the present method.

Reduced-Order Models.
In this section, we used the eigenvalues, which have been studied in the previous section to construct a reduced-order aerodynamic model.The present ROM is called new ROM based on Γ  .The accuracy of the proposed method is shown in Figures 6 and 7.The first figure (Figure 6) shows the unsteady lift on isolated flat plate airfoil due to step change in airfoil downwash (Wagner problem) predicted using different number of element panels as 3, 5, 7, and 20.The result shows that the present ROM can produce satisfactory results when increasing of element panels.Figure 7 shows the unsteady lift of the airfoil predicted using different number of element panel as 3, 5, 7, and 20.The results show that the present ROM can produce satisfactory results when increasing of element panels.Both models can be properly predicted by the present ROM.

Conclusions and Future Work
In the present work, a reduced-order model for unsteady flows is developed without a static correction technique.This method is proposed to construct the eigensystem based on body vortices, and the rank of the eigensystem is lower than the corresponding eigensystem of the conventional method and the method based on wake vortices.Furthermore, the eigensystem does not contain any zero eigenvalues, so the static correction is not needed.The results show that the proposed method can produce satisfactory results like the conventional method and the method based on wake vortices.For the future work, we plan to apply a pseudoinverse operator to solve nonsquare matrix A 12 and then extend the model for aeroelastic analyses of wings.

Figure 1 :
Figure 1: Lift acting on the airfoil due to step change in airfoil downwash.

Figure 2 :
Figure 2: Unsteady lift due to plunging motion of an isolated airfoil.

Figure 4 :
Figure 4: Eigenvalues for two-dimensional airfoil when increasing the number of element panels.

Figure 5 :
Figure 5: Eigenvalues of vortex lattice of unsteady flow about the three-dimensional wing-◻ conventional method, + method based on Γ  , and I present method.

2Figure 7 :
Figure 7: Lift acting on the airfoil due to step change in airfoil downwash predicted using ROM based on Γ  .