NASA/Langley/Somers-Maughmer NLF-1015 (Natural Laminar Flow) Airfoil
NASA/Langley/Somers-Maughmer NLF-1015 (Natural Laminar Flow) Airfoil
1. Airfoil surfaces (upper and lower), and camber line curve curve-fitting
MATLAB CODE:
[v,T,vT]=xlsread('coordinate.xlsx')
t=v(:,1);y=v(:,2);s=v(:,3);z=v(:,4);
%exporting coordinates
curvefit_upper=fit(t,y,'poly9')
curvefit_lower=fit(s,z,'poly9')
%curve-fitting of two surfaces
a=1:-0.05:-1;
b=transpose(0.5*(1-cos(acos(a))));
%cosine spacing variable
coordinate_c=(curvefit_upper(b)+ curvefit_lower(b))/2;
%camber-line is mid-line of two surfaces
plot(b,curvefit_upper(b),b,curvefit_lower(b),b,coordinate_c)
%complete airfoil plotted
yc_derivative=fit(b,differentiate(fit(b,coordinate_c,'poly5'),b),'po
ly4')
%{
d(yc)/dx for circulation calculation;
obtained by curve-fitiing camber-line coordinates, and
differentiating it.
Derivative values output are for discrete points, i.e. each value of
b,
and thus, curve-fitting is done again to obtain slope of camber-line
equation
%}
OUTPUT:
curvefit_lower =
curvefit_upper =
CODE:
syms x;syms n real;
p=0.5*(1-cos(x));
%cosine spacing variable
camber_slope=@(x)(-0.8977*(p)^4 + 0.7221*(p)^3 + 0.1621*(p)^2 -
0.4721*(p) + 0.1786)*diff(p);
%equation obtained from previous code
a0=(int(camber_slope,x,0,pi))/pi;
an=(int((camber_slope*cos(n*x)),x,0,pi))*(2/pi);
%fourier cofficients;
%an formula can be used for infinite summation if needed
A1=(int((camber_slope*cos(1*x)),x,0,pi))*(2/pi)*sin(1*x);
A2=(int((camber_slope*cos(2*x)),x,0,pi))*(2/pi)*sin(2*x);
A3=(int((camber_slope*cos(3*x)),x,0,pi))*(2/pi)*sin(3*x);
A4=(int((camber_slope*cos(4*x)),x,0,pi))*(2/pi)*sin(4*x);
%four fourier sine series coefficients computed
x=0:5:180;
circulation=2*(0.000774*(1+cosd(x))/sind(x) +0.040136*sind(x) -
0.007793*sind(2*x) );
% γ(θ)=2*[(-A0)((1+cosθ)/sinθ) + Σan sin(n θ) ) ] formula used;
%unit free-stream velocity and zero angle of attack;
%Σ (an sin(nx)) restricted to n=1 & n=2
plot(0.5*(1-cosd(x)),circulation)
%circulation variation with x/c and θ plotted
OUTPUT
a0 =-0.000774
A1 = 0.040136
A2 = -0.007793
PLOT:
1+𝑐𝑜𝑠𝜃
γ(θ)=2[(-A0)( ) + ∑∞
𝑛=1(𝑎𝑛 sin 𝑛𝑥)]
𝑠𝑖𝑛𝜃