- Vortex lattice method
The

**Vortex lattice method**, (VLM), is a numerical,Computational fluid dynamics , method used mainly in the early stages ofaircraft design and inaerodynamic education at university level. The VLM models the lifting surfaces, such as awing , of an aircraft as a infinitely thin sheet of discrete vortices to computelift andinduced drag . The influence of the thickness,viscosity and other things, is neglected.VLM:s 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, with multiple lifting surfaces with taper, kinks, twist, camber, trailing edge control surfaces and many other geometricfeatures.

By simulating the flow field, one can extract the pressure distribution or as in the the case of the VLM, the force distribution, around the simulated body. This knowledge is then used to compute the aerodynamic coefficients andtheir derivatives that are important for assessing the aircrafts handling qualities in the conceptual design phase. With an initial estimate of the pressure distribution on the wing, the structural designers can start designingthe load bearing parts of the wings, fin andtailplane 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 [*NASA, "Vortex-lattice utilization". Nasa SP-405, NASA-Langley, Washington, 1976.*] .

The VLM is the extension ofPrandtl lifting line theory [*Prandtl. L, "Applications of modern hydrodynamics to aeronautics", NACA-TR-116, NASA, 1923.*] , where the wing of an aircraft is modeled as aHorseshoe vortex . The name was coined by V.M. Falkner in his aeronautical research council paper of 1946 [*Falkner. V.M., "The Accuracy of Calculations Based on Vortex Lattice Theory", Rep. No. 9621, British A.R.C., 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 VML 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. 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.

The method is comprehensibly described in all major aerodynamic textbooks, such as Katz & Plotkin [

*J. Katz, A. Plotkin, "Low-Speed Aerodynamics, 2nd ed., Cambridge University Press", Cambridge, 2001.*] , Anderson [*J.D. Anderson Jr, "Fundamentals of aerodynamics", 2nd ed., McGraw-Hill Inc, 1991.*] , Bertin & Smith [*J.J. Bertin, M.L. Smith, "Aerodynamics for Engineers", 3rd ed., Prentice Hall, New Jersey, 1998.*] or Houghton & Carpenter [*E.L. Houghton, P.W. Carpenter, "Aerodynamics for Engineering Students", 4th ed., Edward Arnold, London, 1993.*]**Theory**The vortex lattice methods is build 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. The method neglect 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:

*The flow field isincompressible ,inviscous andirrotational .

*The lifting surfaces are thin.The influence of thickness on aerodynamic forces are neglected

*Theangle of attack is small,small angle approximation .**Method**By the above assumptions the flowfield is

Conservative vector field , which means that there exists a velocity potential:

$ar\{v\}\; =\; abla\; phi$and that

Laplace's equation hold.Laplace’s equation is a second order linear equation, and being so it is subject to the principle of superposition. Which means that if $y\_\{1\}$ and $y\_2$ are two solutions ofthe differential equation $L(Y)$, then the linear combination $c\_1y\_1\; +\; c\_2y\_2$ is also a solution for any values of the constants $c\_1$ and $c\_2$. As Anderson [

*J.D. Anderson Jr, "Fundamentals of aerodynamics", 2nd ed., McGraw-Hill Inc, 1991.*] 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 thepoint source , thepoint sink , thedoublet and thevortex 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.**Aircraft Model**The lifting surfaces of an aircraft is divided into several panels. A horseshoe vortex is applied on each of these panels and the velocity vector generated by the vortices at the collocation points of each panel is computed. The vortex is places at the 1/4 chord point of each panel, and the collocation point at 3/4 chord. For a problem with $n$ panels, the induced velocity of each unit strength vortex on each panel is collected in the influence matrix $W$

$mathbf\{W\}\; =\; egin\{bmatrix\}w\_\{11\}\; w\_\{12\}\; cdots\; w\_\{1n\}\backslash w\_\{11\}\; ddots\; vdots\; \backslash vdots\; ddots\; vdots\; \backslash w\_\{n1\}\; cdots\; cdots\; w\_\{nn\}\; end\{bmatrix\}$

A

Neumann boundary condition is applied, which prescribes that the normal velocity at the boundary is zero. It is also known as the flow tangency condition, or no cross flow condition. It means that at the boundary (e.g the surface of a wing) the flow must be parallel to the surface. The following system of equations may be set up. The right hand side is formed by the freestream and theangle of attack so that $b=v\_\{infty\}\; sin(alpha)$.$egin\{bmatrix\}w\_\{11\}\; w\_\{12\}\; cdots\; w\_\{1n\}\backslash w\_\{21\}\; ddots\; vdots\; \backslash vdots\; ddots\; vdots\; \backslash w\_\{n1\}\; cdots\; cdots\; w\_\{nn\}\; end\{bmatrix\}$

egin{bmatrix}Gamma_{1} \Gamma_{2} \vdots \Gamma_{n} end{bmatrix}=

egin{bmatrix}b_{1} \b_{2} \vdots \b_{n} end{bmatrix}

From this system of equations the strength $Gamma$ of the vortices can be solved for, and the forces acting on the panels be computed with

$F=\; ho\_\{air\}\; Gamma\; (V\_\{infty\}+V\_\{induced\})\; l$

**References****ources***NASA, "Vortex-lattice utilization". Nasa SP-405, NASA-Langley, Washington, 1976.

*Prandtl. L, "Applications of modern hydrodynamics to aeronautics", NACA-TR-116, NASA, 1923.

*Falkner. V.M., "The Accuracy of Calculations Based on Vortex Lattice Theory", Rep. No. 9621, British A.R.C., 1946.

*J. Katz, A. Plotkin, "Low-Speed Aerodynamics, 2nd ed., Cambridge University Press", Cambridge, 2001.

*J.D. Anderson Jr, "Fundamentals of aerodynamics", 2nd ed., McGraw-Hill Inc, 1991.

*J.J. Bertin, M.L. Smith, "Aerodynamics for Enginners", 3rd ed., Prentice Hall, New Jersey, 1998.

*E.L. Houghton, P.W. Carpenter, "Aerodynamics for Engineering Students", 4th ed., Edward Arnold, London, 1993.

*Wikimedia Foundation.
2010.*