CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

Help: Fluid Flow Over a Horizontal Flat Plate by Boundary Layer Eq.

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 20, 2011, 05:57
Default Help: Fluid Flow Over a Horizontal Flat Plate by Boundary Layer Eq.
  #1
New Member
 
Egemen
Join Date: May 2011
Posts: 1
Rep Power: 0
egemen is on a distinguished road
Hi everyone,

i faced up a problem on my matlab code related to fluid flow over a horizontal plate by boundary layer theorem. I found nodal equations and wrote a script to solve it but it didnt work. I really dont understand why it did not work. I controlled it many times and i could not find any problem with code. I will be crazy. Are there someone can help me on this code?

Boundary Layer Equations for Fluid Flow

1.u(du/dx)+v(du/dy)=mu(d^2u/dy^2) (momentum eq.)
2.du/dx+dv/dy=0 (continuity eq.)

Code at below;

clear all
clc
x1=0; %initial x value
x2=0.1; %final x value (lenght of the domain)
y1=0; %initial y value
y2=0.1; %final y value
dv=1.307*10^-6; %kinematic viscosity
n=100; %number of volumes and also nodes
g=5; %initial guess
Q=g*ones(n,n); %initial Q values for each node
U=g*ones(n+1,n+1); %initial U values at surfaces
V=g*ones(n+1,n+1); %initial V values at surfaces
Uinf=4; %U infinity m/s
maxerrq=1;
maxerrv=1;
h=(x2-x1)/n; %dx=dy=h the number of nodes
D=dv/h; %kinematic viscosity over step
k=0;

while maxerrq>10^-6 && maxerrv>10^-6
maxerrq=0;
maxerrv=0;
k=k+1;
Qold=Q;

for j=1:n
U(1,j)=Uinf;
U(j,n)=Uinf;
end

for i=1:n
V(i,1)=0;
U(i,1)=0;
end

Vold=V;
Uold=U;

for i=2:n-1
for j=2:n-1
api=(U(i+1,j)+V(i,j+1)+2*D);
aei=0;
awi=U(i,j);
ani=D;
asi=D+V(i,j);
Si=0;
Q(i,j)=(aei*Q(i+1,j)+awi*Q(i-1,j)+ani*Q(i,j+1)+asi*Q(i,j-1)+Si)/api;
end
end

% corners

% lower-left corner's coefficients & nodal eq. (B1)
ap1=U(2,1)+V(1,2)+3*D;
ae1=0;
aw1=0;
an1=D;
as1=0;
S1=U(1,1)*Uinf;
Q(1,1)=(ae1*Q(2,1)+an1*Q(1,2)+S1)/ap1; %B1
% --------------------------------------------
% upper-left corner's coef. & nodal eq. (B3)
ap3=U(2,n)+3*D;
ae3=0;
aw3=0;
an3=0;
as3=V(1,n)+D;
S3=(U(1,n)-V(1,n+1)+2*D)*Uinf;
Q(1,n)=(ae3*Q(2,n)+as3*Q(1,n-1)+S3)/ap3; %B3
% --------------------------------------------
% upper-right corner's coef. & nodal eq. (B5)
ap5=U(n+1,n)+3*D;
ae5=0;
aw5=U(n,n);
an5=0;
as5=V(n,n)+D;
S5=2*D*Uinf;
Q(n,n)=(aw5*Q(n-1,n)+as5*Q(n,n-1)+S5)/ap5;
% --------------------------------------------
% lower-right corner's coef. & nodal eq. (B7)
ap7=U(n+1,1)+V(n,2)+3*D;
ae7=0;
aw7=U(n,1);
an7=D;
as7=0;
S7=0;
Q(n,1)=(aw7*Q(n-1,1)+an7*Q(n,2)+S7)/ap7;
% --------------------------------------------

% edges

% left & right edges (B2 & B6)
for j=2:n-1
Q(1,j)=(D*Q(1,j+1)+(D+V(1,j))*Q(1,j-1))/(U(2,j)+V(1,j+1)+2*D); %for B2
Q(n,j)=(U(n,j)*Q(n-1,j)+D*Q(n,j+1)+(V(i,j)+D)*Q(n,j-1))/(U(n+1,j)+V(n,j+1)+2*D); %for B6
end
% ---------------------------------------------
% lower & upper edges (B8 & B4)
for i=2:n-1
Q(i,1)=(U(i,1)*Q(i-1,1)+D*Q(i,2))/(U(i+1,1)+V(i,2)+3*D); %for B8
Q(i,n)=((V(i,n)-2*D)*Q(i,n-1)+(U(i,n)-V(i,n+1)+2*D)*Uinf)/(U(i+1,n)+3*D); %for B4
end
% ---------------------------------------------

for j=1:n
for i=2:n
U(i,j)=(Q(i-1,j)+Q(i,j))/2; %average U velocity at surfaces
end
end

for i=1:n
for j=2:n+1
V(i,j)=U(i,j-1)-U(j,i)+V(i,j-1); %continuity eq.
end
end

for i=1:n
for j=1:n
errq=abs(Qold(i,j)-Q(i,j));
errv=abs(Vold(i,j)-V(i,j));
if errq>maxerrq
maxerrq=errq;
end
if errv>maxerrv
maxerrv=errv;
end
end
end
end
egemen is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Wind turbine simulation Saturn CFX 60 July 17, 2024 06:45
A question on "Specifying Boundary Layer Deformation Smoothing" didiean FLUENT 2 January 16, 2012 22:39
Concentric tube heat exchanger (Air-Water) Young CFX 5 October 7, 2008 00:17
Free convection flow over vertical flat plate Polly Main CFD Forum 1 February 11, 2003 14:25
Terrible Mistake In Fluid Dynamics History Abhi Main CFD Forum 12 July 8, 2002 10:11


All times are GMT -4. The time now is 13:26.