Contents

• Purpose
• Derivative calculation
• Test case
• Stability derivatives
  • VLM, no fuselage
  • Influence of the panel method
  • Thick surface analysis, with fuselage
• Control derivatives
• Conclusion and recommendations

Back to main table of contents

---

Purpose

The goal of this study is to

• Benchmark flow5's prediction of stability and control derivatives with the predictions made by AVL.
• Perform parametric analyses to determine the sensivity of the derivatives to variables such as the angle of attack, the inclusion of the fuselage in the model and the type of panel method.
• Establish guidelines for the calculation of the derivatives.

Back to top

---

Derivative calculation

The methods to calculate the derivatives differ in AVL and flow5. To the author's understanding, the methods implemented in AVL are the ones described in "Flight Vehicle Aerodynamics", Mark Drela 2014.

Stability derivatives

AVL

The method takes advantage of the linear dependency of the panel forces with the 6 degrees of freedom in motion. This allows to express the stability derivatives explicitly, and to and solve them by direct back-substituion and linear combination which makes their calculation fast and accurate.

flow5 v7.13

The derivatives are evaluated using finite differences.

The method used by AVL could be implemented for flow5's VLM. However it is not transposable to the thick surface panel methods in which the pressure forces are derived from the surface gradient of the doublet densities and do not depend linearly on the freestream's angle of attack and sideslip.
Solutions to this problem are being explored, but are not implemented as of flow5 v7.13.

Control derivatives

The rotations of the control surfaces change both the linear system's matrix and the right hand side (RHS) vector. Strictly speaking, the problem should therefore be built and solved entirely for each derivative. This however would be computationally very expensive. A standard procedure is to assume that the change to the influence matrix is small, so that the LU-factorized matrix can be reused as-is, and that only the change to the RHS needs to be implemented. This simplification reduces the computations to the construction of the RHS and to its back-substitution.
This is how AVL, flow5 and other panel codes proceed. A numerical experiment was carried out with flow5 to evaluate the error implicit in this approximation and is presented herefafter.

AVL

The boundary condition are linearized since only the first order development is required, and the derivatives are calculated explicitely and solved by direct back-substitution and linear combination of deflections.

flow5 v7.13

The problem cannot be linearized for the same reason as in the case of the stability derivatives. The derivatives are calculated by finite differences, with an RHS vector being constructed and back-substituted for each derivative.

Back to top

---

Test case

Model validation: linear analyses

The purpose of these first runs is to check that the analyses are set up identically in flow5 and AVL.

Type 1 comparisons at β=0°

Both AVL and flow5 predict that the piching moment is zero at α = 1.07°.

Type 5 comparisons at α=0°

Stability derivatives

VLM, no fuselage

The intent of this first run is to compare the stability derivatives using the VLM common to both programs.

The stability derivatives are compared at the zero pitching moment angle, i.e. α = 1.08° for both programs.
AVL and flow5 do not use the same sign conventions, so the deviation is calculated with respect to the normalized values.

The derivatives are in reasonable agreement except for Cnp and Clr. To understand the discrepancy, the derivatives have been plotted as a function of the angle of attack.

As it appears, the value of the derivatives Cnp and Clr estimated by AVL depend strongly on the angle of attack, which makes them very sensitive to the point of evaluation. No obvious explanation could be found for this difference of sensitivity exhibited by both programs.

The other lateral derivatives are sensitive to the angle of attack to various degrees.

On the other hand, the longitudinal derivatives are not sensitive to the angle of attack.

Influence of the panel method

	Quads thick	Quads thin	Tri-uni thick	Tri-uni thin	VLM2
CXu	-0.038923	-0.014213	-0.0034009	-0.0010699	-0.0018653
CLu	-0.00090185	-0.00013729	-9.2337E-05	-8.8636E-06	2.6189E-05
Cmu	-1.7561E-06	-6.6029E-08	-0.0047468	-6.4618E-09	1.0353E-07
CXa	-0.64907	-0.33282	0.16031	0.085434	0.12196
CLa	-5.5324	-5.4843	-5.7636	-5.6188	-5.4932
Cma	-1.4773	-1.132	-1.6562	-1.2149	-1.4485
CXq	-1.2981	-0.31383	-0.38311	-0.1022	-0.25743
CLq	-9.2422	-9.1713	-9.7862	-9.5477	-9.4901
Cmq	-17.133	-10.815	-16.677	-10.693	-15.856
CYb	-0.13164	-0.12969	-0.12031	-0.11863	-0.12915
Clb	-0.075016	-0.054906	-0.058189	-0.03846	-0.056348
Cnb	0.048835	0.028517	0.04353	0.026953	0.042822
CYp	-0.097317	-0.096193	-0.098104	-0.094298	-0.093629
Clp	-0.60218	-0.48104	-0.57314	-0.38189	-0.55689
Cnp	-0.011015	-0.0040144	-0.011683	-0.0029192	-0.002397
CYr	0.0975	0.096827	0.088485	0.088861	0.096879
Clr	0.015167	0.0081044	-0.023376	-0.010809	0.0074863
Cnr	-0.037425	-0.024704	-0.035301	-0.02339	-0.035711

The scattering of the predictions can be quite large depending on the derivative.

Influence of the fuselage

The inclusion of the fuselage changes the derivatives to various degrees. The variations are not easy to interpret although some can be expected such as:

• the increase of Cxq due to the perturbation of the flow when the fuselage is not aligned with the x axis
• the increase of CYb due to the fuselage's lateral surface

Back to top

---

Control derivatives

Result comparison

The derivatives have been evaluated in flow5 using the VLM2 method to be consistent with AVL.

The control derivatives are only evaluated in flow5 as part of a T7 analysis and are therefore calculated at the zero-pitching moment angle.

Sanity check: influence matrix

A numerical experiment was made in flow5 to evaluate the error implicit in the assumption that the linear system's matrix is not modified by the deflection of the control surfaces. The influence matrix was rebuilt for each of the aileron, elevator and rudder controls, and the problem was entirely solved for the derivatives in each case. The results were compared to those presented above in the VLM case.

It appears that the differences are minor and no larger than the variations seen beteween the VLM2 and the tri-uniform method, which justifies the approximation.

Back to top

---

Conclusion and recommendations

Conclusions

The main conclusions of this sensitivity study are that

• the derivatives exhibit various degrees of variability to the type of analysis
• the evaluation of the derivatives using panel methods are less precise than are the estimation of the lift, drag and moment coefficients.

Recommendations

• use preferably the estimations in a relative way to plot tendencies rather than in an absolute way
• use the same type of models and panel methods when comparing derivatives.

Back to top

---