Back to **Aerodynamic Design of Aircraft with Computational Software**

Back to overview: **Preamble Exercises and Projects**

Chapter 3[pdf]: 3.1 Introduction - Mapping Planform to Lift and Drag 3.2 Computational Vortex-Lattice Methods 3.3 Planform Design Studies with VLM 3.4 Wings for High Speed

**Tutorial: **—

**Review questions**

- Look at the Biot-Savart equation (3.5). Where does it produces infinite velocities? Consider now a VLM model with a symmetric main wing ahead of a vertical tail in the symmetry plane and a horizontal stabilizer. The model assumes that the trailing vortices follow free stream streamlines. Consider the Piper Cub J-3 in the sketch below, from Wikipedia. Sketch in the angles of attack where we expect trouble from trailing vortices hitting collocation points or bound vortex midpoints. Yet the Cub flies just fine. Discuss the effect of the fuselage. In zero sideslip, the wing horseshoe trailing leg in the symmetry plane should not create problems in the force calculation anyway. Why?

- Assume that a pair of wingtip vortices Gand -G, separated by wing span
*b*(actually, they will converge somewhat from the wing tips) are infinitely long and horizontal. How fast do they sink – or rise? Using the simplest constant-strength Prandtl lifting line model, compute G for a Boeing 747 at take-off weight around 350 metric tons. You will need more data, easily available. Constant strength is wrong; re-do the calculation assuming elliptic lift distribution; that model sheds a variable-strength wake. Assume that it rolls up so all vorticity from the right half is concentrated to the wing tip trailing streamline. What is its strength?

**VLM computations**

The script `tor1.m`

computes a quadrilateral wing defined by *AR*, taper ratio *ct/cr* , quarter-chord sweep angle Sweep, dihedral Dihed, and Twist for some airspeed, altitude, and angle of attack a. It produces a row in a table of the parameters (angles in radians)

AR ct/cr Sweep Dihed. Twist alpha

and a corresponding row in a table of the force and moment coefficients

CL CD CY Cl Cm Cn

`tor1(AR,tap,swe,dih,twi,alt,airspeed,alpha,foil1,foil2,resimn,resimx)`

`resimn`

and `resimx`

define refinements of the initial paneling:

nx = 2^{(res-1) }nx0 for res = resimn… resimx.

so `resimn` = 2 and `resimn` = 4 runs original, doubled, and quadrupled paneling used to extrapolate the results, in `VLM_e.m`

described below.

- Check how the VLM estimates for
*C*converge with the number of vortices. Run with script tor1.m an_{L}, C_{m}*AR*8 rectangular wing with symmetric profile and

*Nx*= 2^{k},*k*= 1,….. and*Ny*= 4*Nx*, a = 5 deg.

Use`resimn`=`resimx`= 2,3,4,5. Plot*C*and_{L}, C_{m}*C*vs. 1/_{D}*Nx.*Theory shows that

*C*(1/_{L}*Nx*) =*C*(0) + C1/_{L}*Nx*+ smaller terms.

This allows extrapolation to*Nx*= Infty:

*C*(0) = 2_{L}*C*(1/(2_{L}*Nx*)) –*C*(1/_{L}*Nx*)

is usually a better approximation than even*C*(1/(4_{L}*Nx*)) and much cheaper since the work grows like*Nx*^{4}. The extrapolation can be continued with results from 4 Nx, 8 Nx, … (but 8 Nx costs 4000 x the compute time of Nx). Systematic mesh refinement with at least three levels is a due diligence check that the convergence follows the theory so bad results from under-resolution, ill-conditioning, or “close encounters” can be spotted. The function`VLM_e.m`

performs calculations of force and moment coefficients for a sequence of refinements of a user-defined initial coarse paneling.

- Then try a straight wing with
*AR*2, 4, 6, 8, 10. Run an alpha sweep withtor1.m and check d*C*/da (small a) vs. the formula in Ch 3. The script ARsweep.m runs a list of AR and a using_{L}`tor1.m`

which in turn uses VLM_e.m and produces tables partab and restab with rows described above. Your job is to compute*C*_{L,}_{a}from the table.

- The script
`tor1.m`

computes a quadrilateral wing with sweep, taper, and twist. Your job is to experiment with sweep*and*taper for an untwisted wing. Modify ARsweep.m to make an (AR, taper, a) sweep and fit

*C*+_{Di}= C_{D0i}*k C*_{L}^{2}

to the computed points for each (*AR*,taper). The wing (induced drag) efficiency factor is*e*=*k*/(p*AR*) where lifting line theory predicts that*e*= 1 for an elliptical flat wing._{max}

Fix*AR*at 6. Then run a quadrilateral tapered wing with*ct/cr*= [0.3, 0.7, 1] and quarter-chord sweeps [0,30,60]^{o}. Note the*e*and*C*so computed as well as_{D0i}*C*_{L,}_{a}. Before analyzing the data, add a rough estimate of*C*_{D0}_{,viscous}of 150 cts to*C*. Show that the best_{D0i}*L/D*(which determines best glide slope) is LoD_{max}= 1/sqrt(*k**C*). Which wing has best glide slope?_{D0tot}

- The file
`GAV.m`

has a Tornado model for a Piper Cub-like generic general aviation airplane and its weight and CG-location. What angle of attack and elevator angle will trim it for straight and level flight at sea level and TAS 35 m/s? Suppose*C*= 0.0150. What thrust is required? First, neglect the elevator influence on lift to find a, then for this a find the correct elevator angle_{D0v}*ele*. A more accurate procedure is to assume “linear aerodynamics”. Compute the derivatives*C*_{L}_{,a},*C*_{m}_{,a},*C*_{L}_{,ele}and*C*_{m}_{,ele}from a computed table produced by (a,*ele*)-sweep and solve the linear system

qS(0 +αC_{L}_{,α}+eleC) = W, qScMAC (0 +_{L,ele}αC_{m}_{,α}+eleC_{m}_{,ele}) = 0