Scribd
Upload
15 views7 pages

Clear All

The document presents the numerical solution of the Blasius equation to analyze laminar boundary layer flow over a flat plate. It defines the governing equations and initial conditions, numerically solves the problem, and plots the velocity and temperature profiles as well as the local Nusselt number along the plate.

Uploaded by

Ibiza Spotify
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
15 views7 pages

Clear All

The document presents the numerical solution of the Blasius equation to analyze laminar boundary layer flow over a flat plate. It defines the governing equations and initial conditions, numerically solves the problem, and plots the velocity and temperature profiles as well as the local Nusselt number along the plate.

Uploaded by

Ibiza Spotify
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 7

%clear all

xf=1.9738;
Pr=0.72;
Re=5000;
Ec=0;
t=0.1;
N=30;
ymax=7;
h=ymax/N;
x0=0.01*t;
beta=(t^2/(1-t^2)^1.5)*(atanh((1-t^2)^0.5)-(1-t^2)^0.5);
X0=-1+(x0/t)^2/2;
R0=t*(1-X0^2)^0.5;
ue0=1/(1-beta)*R0/(R0^2+t^4*X0^2)^0.5;
t0=1+Ec*(1-ue0^2)/2;
% Asignamos la componente N a Xo
u0(N)=X0;
%Vector de temperaturas iniciales
u0(N+1:2*N-1)=t0;
%Asigamos la componente 2*N+1 a cte resistencia de forma
u0(2*N+1)=0;
%Asigamos la componente 2*N+2 a cte resistencia de friccion
u0(2*N+2)=0;
%necesito un vector de velocidades iniciales para cada linea
for j=1:N-1
ylin(j,1)=j*h;
%clear all

xf=1.9738;
Pr=0.72;
Re=5000;
Ec=0;
t=0.1;
N=30;
ymax=7;
h=ymax/N;
x0=0.01*t;
beta=(t^2/(1-t^2)^1.5)*(atanh((1-t^2)^0.5)-(1-t^2)^0.5);
X0=-1+(x0/t)^2/2;
R0=t*(1-X0^2)^0.5;
ue0=1/(1-beta)*R0/(R0^2+t^4*X0^2)^0.5;
t0=1+Ec*(1-ue0^2)/2;
% Asignamos la componente N a Xo
u0(N)=X0;
%Vector de temperaturas iniciales
u0(N+1:2*N-1)=t0;
%Asigamos la componente 2*N+1 a cte resistencia de forma
u0(2*N+1)=0;
%Asigamos la componente 2*N+2 a cte resistencia de friccion
u0(2*N+2)=0;
%necesito un vector de velocidades iniciales para cada linea
for j=1:N-1
ylin(j,1)=j*h;
u0(j)=ue0;

end
Nx=100;
for i=1:Nx,
tspan(1,i)=x0+(xf-x0)*(i-1)/(Nx-1);
end
[x,u]=ode45(@blasiuscompleto,tspan,u0,[],h,N,t,beta,Re,Ec,Pr);
u(100,1);
%Coeficientes de resistencia
Cdp=u(100,2*N+1)-(u(100,N-1)^2-1)*(1-xf^2)
Cdf=u(100,2*N+2)
figure
plot(x(:,1),u(:,1))
xlabel('x adimensional');
ylabel('velocidad adimensional ');
%-----------------Perfil de velocidades------------------[m,n]=size(u);
ceros=zeros(m,1);
vel0=[ceros,u];
ylin0=[0;ylin];
figure
hold on;
plot(vel0(1,1:30),ylin0(:,1)),'y';
plot(vel0(26,1:30),ylin0(:,1),'r');
plot(vel0(51,1:30),ylin0(:,1),'c');
plot(vel0(76,1:30),ylin0(:,1),'m');

plot(vel0(100,1:30),ylin0(:,1),'k');
%Title('Perfil de velocidades para t=0.3')
xlabel('Velocidad adimensional');
ylabel('y adimensional');
legend('estacion 1','estacion 25','estacion 50 ','estacion75','estacion
100');
%-----------------------Perfiles de temperaturas-----------------T=[ceros,u(:,(N+1):(2*N-1))];
figure
subplot (2,2,1)
hold on;
plot(T(5,:),ylin0(:,1),'y')
plot(T(26,:),ylin0(:,1),'r')
plot(T(51,:),ylin0(:,1),'c')
plot(T(76,:),ylin0(:,1),'b')
plot(T(95,:),ylin0(:,1),'k');
xlabel('Temperatura adimensional');
ylabel('y adimensional ');
legend('estacion 5','estacion 25','estacion 50 ','estacion75','estacion
95');
for i=1:100
Nux(i)=(u(i,N)+1)*sqrt(Re)*u(i,N+1)/h;
end
figure
subplot(2,2,1)
hold on
plot(x,Nux);
xlabel('X'),ylabel('Nusselt')

title('Nusselt para t=0.1')


% Nusselt tablas
for i=1:100
Rex(i)=Re*(1+u(i,N));
end
Nutab=0.332*Rex.^(1/2)*Pr^(1/3);
plot(x,Nutab,'r'),xlabel('X'),ylabel('N
u0(j)=ue0;
end
Nx=100;
for i=1:Nx,
tspan(1,i)=x0+(xf-x0)*(i-1)/(Nx-1);
end
[x,u]=ode45(@blasiuscompleto,tspan,u0,[],h,N,t,beta,Re,Ec,Pr);
u(100,1);
%Coeficientes de resistencia
Cdp=u(100,2*N+1)-(u(100,N-1)^2-1)*(1-xf^2)
Cdf=u(100,2*N+2)
figure
plot(x(:,1),u(:,1))
xlabel('x adimensional');
ylabel('velocidad adimensional ');
%-----------------Perfil de velocidades------------------[m,n]=size(u);
ceros=zeros(m,1);
vel0=[ceros,u];

ylin0=[0;ylin];
figure
hold on;
plot(vel0(1,1:30),ylin0(:,1)),'y';
plot(vel0(26,1:30),ylin0(:,1),'r');
plot(vel0(51,1:30),ylin0(:,1),'c');
plot(vel0(76,1:30),ylin0(:,1),'m');

plot(vel0(100,1:30),ylin0(:,1),'k');
%Title('Perfil de velocidades para t=0.3')
xlabel('Velocidad adimensional');
ylabel('y adimensional');
legend('estacion 1','estacion 25','estacion 50 ','estacion75','estacion
100');
%-----------------------Perfiles de temperaturas-----------------T=[ceros,u(:,(N+1):(2*N-1))];
figure
subplot (2,2,1)
hold on;
plot(T(5,:),ylin0(:,1),'y')
plot(T(26,:),ylin0(:,1),'r')
plot(T(51,:),ylin0(:,1),'c')
plot(T(76,:),ylin0(:,1),'b')
plot(T(95,:),ylin0(:,1),'k');
xlabel('Temperatura adimensional');
ylabel('y adimensional ');

legend('estacion 5','estacion 25','estacion 50 ','estacion75','estacion


95');
for i=1:100
Nux(i)=(u(i,N)+1)*sqrt(Re)*u(i,N+1)/h;
end
figure
subplot(2,2,1)
hold on
plot(x,Nux);
xlabel('X'),ylabel('Nusselt')
title('Nusselt para t=0.1')
% Nusselt tablas
for i=1:100
Rex(i)=Re*(1+u(i,N));
end
Nutab=0.332*Rex.^(1/2)*Pr^(1/3);
plot(x,Nutab,'r'),xlabel('X'),ylabel('N

You might also like