Vortex lattice method


The Vortex lattice method,, is a numerical method used in computational fluid dynamics, mainly in the early stages of aircraft design and in aerodynamic education at university level. The VLM models the lifting surfaces, such as a wing, of an aircraft as an infinitely thin sheet of discrete vortices to compute lift and induced drag. The influence of the thickness, viscosity is neglected.
VLMs can compute the flow around a wing with rudimentary geometrical definition. For a rectangular wing it is enough to know the span and chord. On the other side of the spectrum, they can describe the flow around a fairly complex aircraft geometry.
By simulating the flow field, one can extract the pressure distribution or as in the case of the VLM, the force distribution, around the simulated body. This knowledge is then used to compute the aerodynamic coefficients and their derivatives that are important for assessing the aircraft's handling qualities in the conceptual design phase. With an initial estimate of the pressure distribution on the wing, the structural designers can start designing the load-bearing parts of the wings, fin and tailplane and other lifting surfaces. Additionally, while the VLM cannot compute the viscous drag, the induced drag stemming from the production of lift can be estimated. Hence as the drag must be balanced with the thrust in the cruise configuration, the propulsion group can also get important data from the VLM simulation.

Historical background

John DeYoung provides a background history of the VLM in the NASA Langley workshop documentation SP-405.


The VLM is the extension of Prandtl's lifting-line theory, where the wing of an aircraft is modeled as an infinite number of Horseshoe vortices. The name was coined by V.M. Falkner in his Aeronautical Research Council paper of 1946. The method has since then been developed and refined further by W.P. Jones, H. Schlichting, G.N. Ward and others.
Although the computations needed can be carried out by hand, the VLM benefited from the advent of computers for the large amounts of computations that are required.
Instead of only one horseshoe vortex per wing, as in the lifting line theory, the VLM utilizes a lattice of horseshoe vortices, as described by Falkner in his first paper on this subject in 1943. The number of vortices used vary with the required pressure distribution resolution, and with required accuracy in the computed aerodynamic coefficients. A typical number of vortices would be around 100 for an entire aircraft wing; an Aeronautical Research Council report by Falkner published in 1949 mentions the use of an "84-vortex lattice before the standardisation of the 126-lattice".
The method is comprehensibly described in all major aerodynamic textbooks, such as Katz & Plotkin, Anderson, Bertin & Smith Houghton & Carpenter or Drela,

Theory

The vortex lattice method is built on the theory of ideal flow, also known as Potential flow. Ideal flow is a simplification of the real flow experienced in nature, however for many engineering applications this simplified representation has all of the properties that are important from the engineering point of view. This method neglects all viscous effects. Turbulence, dissipation and boundary layers are not resolved at all. However, lift induced drag can be assessed and, taking special care, some stall phenomena can be modelled.

Assumptions

The following assumptions are made regarding the problem in the vortex lattice method:
By the above assumptions the flowfield is Conservative vector field, which means that there
exists a perturbation velocity potential such that the total velocity vector is given by
and that satisfies Laplace's equation.
Laplace’s equation is a second order linear equation, and being so it is subject
to the principle of superposition. Which means that if and are two solutions of
the linear differential equation, then the linear combination is also a solution for any values of the constants and. As Anderson put it "A complicated flow pattern
for an irrotational, incompressible flow can be synthesized by adding together a number of elementary flows, which are also irrotational and incompressible.”. Such elementary flows are the point source or sink, the doublet and the vortex line, each being a solution of Laplace’s equation. These may be superposed in many ways to create the formation of line sources, vortex sheets and so on. In the Vortex Lattice method, each such elementary flow is the velocity field of a horseshoe vortex with some strength.

Aircraft Model

All the lifting surfaces of an aircraft are divided into some number of quadrilateral panels, and a horseshoe vortex and a collocation point are placed on each panel. The transverse segment of the vortex is at the 1/4 chord position of the panel, while the collocation point is at the 3/4 chord position. The vortex strength is to be determined. A normal vector is also placed at each collocation point, set normal to the camber surface of the actual lifting surface.
For a problem with panels, the perturbation velocity at collocation point is given by summing the contributions of all the horseshoe vortices in terms of an Aerodynamic Influence Coefficient matrix.
The freestream velocity vector is given in terms of the freestream speed and the angles of attack and sideslip,.
A Dirichlet boundary condition is applied at each collocation point, which prescribes that the normal velocity across the camber surface is zero.
This is also known as the flow tangency condition. By evaluating the dot products above the following system of equations results. The new normalwash AIC matrix is, and the right hand side is formed by the freestream speed and the two aerodynamic angles
This system of equations is solved for all the vortex strengths. The total force vector and total moment vector about the origin are then computed by summing the contributions of all the forces on all the individual horseshoe vortices, with being the fluid density.
Here, is the vortex's transverse segment vector, and is the perturbation velocity at this segment's center location .
The lift and induced drag are obtained from the components of the total force vector. For the case of zero sideslip these are given by