Question

please use matlab to solve

Problem # 3 P-3 Flow between two paralle plates is described by the following equation dith boundary conditons given as u,-0 & u,-o Calculate the velocity profile using the shooting method for solving the given BVP and compare your results by plotting the numerical solution over the plot of the analytical solution described by: (y-F )where ğr--0000025.H-О75 and h,30 Hint: use 1.75 for the first initial slope, and the other one is 0.45 to 0.5.

Answer 1

j_max = floor(t/dt+1); % Number of time nodes
i_max = floor(H/dy+1); % Number of space nodes

for j = 1:j_max
    for i = 1:i_max
        if j == 1
            if i == i_max
                U(i,j) = Iu(2);
            else
                U(i,j) = Iu(1);
            end
        else
            if i == i_max
                U(i,j) = Bu(2);
            elseif i == 1
                U(i,j) = Bu(1);
            else
                U(i,j) = (1-2*D)*U(i,j-1)+D*(U(i+1,j-1)+U(i-1,j-1)); 
            end
        end
    end
end

%% Plots
figure
hold on
grid on
for j = 1:j_max
plot(U(:,j),0:dy:dy*(i_max-1))
end
xlabel('Velocity "U" (m)')
ylabel('Hight "Y" (m)')
title('Flow Between Two Parallel Plates Problem')