Anyone familiar with the Blasius Equation for boundary layer thickness? I have rewritten it as an ODE through the substitution of the stream function. I now need to solve it numerically in Matlab using the iterative integral method. My code is below, but it is definitely not producing the correct answers. I have no idea what might be wrong with it. I have attached an Excel sheet with the correct answers; independent variable is eta while dependent variable, which I'm trying to solve for in my code below, is f^. My code below spits out a linear f^. However, it should be a curved line as shown in the attached Excel sheet.
Thanks in advance, M Ridzon
clear all
eta=0:0.01:10;
sz=size(eta);
lg=length(eta);
f=zeros(sz);
f1=zeros(sz);
f2=zeros(sz);
f3=zeros(sz);
f4=zeros(sz);
fx=zeros(sz);
results=zeros(lg,2);
results(1,1)=eta(1);
results(1,2)=0;
for_count=0;
for k=2:lg
for_count=for_count+1
if k==2
fx(k-1)=eta(k-1);
f(k-1)=0;
f1(k-1)=0;
f2(k-1)=1;
f3(k-1)=0;
end
fx_guess(k)=eta(k);
fcheck=1;
while_count=0;
while fcheck~=0
while_count=while_count+1
f(k)=((eta(k)-eta(1)).*((fx_guess(k)+fx(1))./2));
f1(k)=-0.5.*((eta(k)-eta(1)).*((f(k)+f(1))./2));
f2(k)=exp(f1(k));
f3(k)=(eta(k)-eta(1)).*((f2(k)+f2(1))/2);
f4(k)=(eta(lg)-eta(1)).*((f2(k)+f2(1))/2);
fx(k)=f3(k)./f4(k);
fcheck=fx(k)-fx_guess(k);
fx_guess(k)=fx(k);
end
results(k,1)=eta(k);
results(k,2)=fx_guess(k);
end
