기초전산유체역학 과제3

소개글

"기초전산유체역학 과제3"에 대한 내용입니다.

목차

1. Perform Stability analysis of the continuity equation.
2. Find u, v at all grid points & Plot the velocity profile as a function of y at 3 locations.
3. Using the numerical solution, find the boundary layer thickness (99%) as a function of x. Compare the result with the Blasius solution. & Find the displacement thickness as a function of x.

본문내용

2. Find u, v at all grid points & Plot the velocity profile as a function of y at 3 locations.
clear;

L=1; h=0.025; dx=0.0001; dy=0.0001; U=10; nu=1.4*10^(-5); % Domain
nx=L/dx+1; ny=h/dy+1; % number of node

u = zeros(ny,nx); v = zeros(ny,nx);

% Boundary condition
u(:,1) = U; % Initial flow
u(end,:) = U; % Far from field
u(1,1)=0; u(end,1)=0;

% Caclculate u & v
for i = 1:nx-1
for j = 1:ny-2
u(j+1,i+1) = u(j+1,i)+((nu*dx)/(dy)^2)*((u(j+2,i)-2*u(j+1,i)+u(j,i))./u(j+1,i))...
-(dx/dy)*(v(j+1,i)./u(j+1,i))*((u(j+2,i)-u(j,i))./2);
end
for j = 1:ny-2
v(j+1,i+1) = v(j,i+1)-(1/2)*(((u(j,i+1)-u(j,i))./dx)+((u(j+1,i+1)-u(j+1,i))./dx))*dy;
end
end

참고 자료

없음

압축파일 내 파일목록

assignment 3.hwp