Lifting Line Theory • Applies to large aspect ratio unswept wings at small angle of attack. • Developed by Prandtl and Lanchester during the early 20 th century. • Relevance – Analytic results for simple wings – Basis of much of modern wing theory (e.g. helicopter rotor aerodynamic analysis, extends to vortex lattice method,) – Basis of much of the qualitative understanding of induced drag and aspect ratio 1875-1953 1868-1946 Γ h V π 4 Γ = h Biot Savart Law: Velocity produced by a semi-infinite segment of a vortex filament Thin-airfoil theory C l =2π(α-α o )
16
Embed
Lifting Line Theory - Virginia Techdevenpor/aoe5104/23. LiftingLineTheory.pdf · Lifting Line Theory • Applies to large aspect ratiounswept wings at small angle of attack. • Developed
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Lifting Line Theory• Applies to large aspect ratio unswept wings at small angle of attack.• Developed by Prandtl and Lanchester during the early 20th century. • Relevance
– Analytic results for simple wings– Basis of much of modern wing theory (e.g. helicopter rotor aerodynamic
analysis, extends to vortex lattice method,)– Basis of much of the qualitative understanding of induced drag and aspect ratio
1875-19531868-1946Γ h
Vπ4Γ
=
h
Biot Savart Law:Velocity produced by a semi-infinite segment of a vortex filament
Thin-airfoil theoryCl=2π(α-αo)
Physics of an Unswept Wingl, Γ
y
pu<pl pu≈ pl
Downwash
s-s
Lift varies across span
Circulation is shed (Helmholz thm)
Vortical wake
Vortical wake induces downwash on wing…
…changing angle of attack just enough to produce variation of lift across span
Simplest Possible Model
Section A-AdiInduced dragb A
A
Wake model Section model
LLT – The Section Model
αε
V∞w
l
ε
di
LLT – The Wake ModelΓ
y1
ydy1
s-s
The Monoplane EquationWake model
Section model∫− −
Γ
=s
s
y
yy
dydyd
yw)(4
)(1
11
π
0( )V c wcπ α α π∞Γ = − −
l, Γ
ys-s
1
1
01
( )4
sy
s
d dydycV cy y
π α α∞−
Γ
Γ = − −−∫
0 π θ
θcos/ −=sy
Substitute for θ, and express Γ as a sine series in θ
∑∞
=∞=Γ
oddnn nAsU
,1)sin(4 θ
+=− ∑
∞
=
θπθθααπ sin4
)sin(sin)(4 ,1
0 scnnA
sc
oddnn The Monoplane Eqn.
Results
∫− −
Γ
=s
s
y
yy
dydyd
yw)(4
)(1
11
π∫−∞
Γ=s
sL dy
SVC 2
∫−∞
Γ=s
sD dyw
SVC
i 22
∑∞
=∞=Γ
oddnn nAsU
,1)sin(4 θSubstituting into
gives 1AARCL π= )1(2
δπ
+=ARCC L
Di
∑∞
=
=oddn
n AAn,3
21)/(δ
1,sin( )
sin
nn odd
nA nw
V
θ
θ
∞
=
∞
=∑
+=− ∑
∞
=
θπθθααπ sin4
)sin(sin)(4 ,1
0 scnnA
sc
oddnn
Solution of monoplane equation
+=− ∑
∞
=
θπθθααπ sin4
)sin(sin)(4 ,1
0 scnnA
sc
oddnn
ys-s0 π θ
θcos/ −=sy
s=2.8; %Half span (distances normalized on root chord)alpha=5*pi/180; %5 degrees angle of attackalpha0=-5.4*pi/180; %Zero lift AoA=-5.4 deg. for Clark YN=20; %N=20 points across half spanth=[1:N]'/N*pi/2; %Column vector of theta's y=-cos(th)*s; %Spanwise positionc=ones(size(th)); %Rectangular wing, so c = c_r everywheren=1:2:2*N-1; %Row vector of odd indices
res=pi*c/4/s.*(alpha-alpha0).*sin(th); %N by 1 result vectorcoef=sin(th*n).*(pi*c*n/4/s+repmat(sin(th),1,N)); %N by N coefficient matrixa=coef\res; %N by 1 solution vector
gamma=4*sin(th*n)*a; %Normalized on uinf and sw=(sin(th*n)*(a.*n'))./sin(th);AR=2*s/mean(c);CL=AR*pi*a(1);CDi=CL^2/pi/AR*(1+n(2:end)*(a(2:end).^2/a(1).^2));