bljones2 - n2=length(theta xsepU=xsep xtranU=xtran...

Info iconThis preview shows pages 1–2. Sign up to view the full content.

View Full Document Right Arrow Icon
function [zeta,s,y,cp,Cl,blU,blL,Cdsy,LoD]=bljones2(xc,yc,xt,yt,delta,alpha,Re) % % hough + transition method + head on a r t jones airfoil % [zeta,s,y,cp,Cl]=exjkcp(xc,yc,xt,yt,delta,alpha); [ n=length(s); %stagnation point( just look in nose region) [cpmax,nstag]=max(cp(fix(n/3):fix(2*n/3))); nstag=nstag+fix(n/3) ; %add starting point n/3 n xcupper=fliplr(s(1:nstag-2)); % start at two points after stagnation point ycupper=fliplr(y(1:nstag-2)); % to avoid +-1 shifts cpupper=fliplr(cp(1:nstag-2)); Uupper=sqrt(1-cpupper); xclower=s(nstag+2:n); yclower=y(nstag+2:n); cplower=cp(nstag+2:n); Ulower=sqrt(1-cplower); U n1=length(xcupper); x(1)=0; x(2:n1)=cumsum(sqrt(diff(xcupper).^2 + diff(ycupper).^2)); % % velocity gradient dUdx=diff(Uupper)./diff(x); dUdx(n1)=dUdx(n1-1); U=Uupper; U [x,delta,delstar,theta,cf,xtran,xsep,Lambda]=bl(x,U,dUdx,Re);
Background image of page 1

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full Document Right Arrow Icon
Background image of page 2
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: [ n2=length(theta); xsepU=xsep; xtranU=xtran; CdsyU=2*theta(n2)*(U(n2))^3.2; C blU(:,1)=x(1:n2)'; blU(:,2)=delta'; blU(:,3)=delstar'; blU(:,4)=theta'; blU(:,5)=cf'; blU(:,6)=cf'; %first fill it up with something blU(1,6)=xtran; blU(2,6)=xsep; b % lower surface x=; n1=length(xclower); x(1)=0; x(2:n1)=cumsum(sqrt(diff(xclower).^2 + diff(yclower).^2)); % % velocity gradient dUdx=diff(Ulower)./diff(x); dUdx(n1)=dUdx(n1-1); U=Ulower; U [x,delta,delstar,theta,cf,xtran,xsep,Lambda]=bl(x,U,dUdx,Re); [ n2=length(theta); xsepL=xsep; xtranL=xtran; CdsyL=2*theta(n2)*(U(n2))^3.2; Cdsy=CdsyU+CdsyL; LoD=Cl/Cdsy; L blL(:,1)=x(1:n2)'; blL(:,2)=delta'; blL(:,3)=delstar'; blL(:,4)=theta'; blL(:,5)=cf'; blL(:,6)=cf'; %first fill it up with something blL(1,6)=xtran; blL(2,6)=xsep;...
View Full Document

{[ snackBarMessage ]}

Page1 / 2

bljones2 - n2=length(theta xsepU=xsep xtranU=xtran...

This preview shows document pages 1 - 2. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online