Scribd
Upload
2 views6 pages

L04 Code Boundary Layer Plate

The document details Lab 4 on Boundary Layer Plate conducted by Sarah Daboub, including data collection and analysis of velocity at various pressure taps. It outlines the calculations for dynamic pressure, velocity, Reynolds number, and boundary layer thickness using MATLAB, along with importing analytical solutions for comparison. The results are visualized through multiple plots comparing experimental and analytical velocities at different stations along the plate.

Uploaded by

sarradaboub
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
2 views6 pages

L04 Code Boundary Layer Plate

The document details Lab 4 on Boundary Layer Plate conducted by Sarah Daboub, including data collection and analysis of velocity at various pressure taps. It outlines the calculations for dynamic pressure, velocity, Reynolds number, and boundary layer thickness using MATLAB, along with importing analytical solutions for comparison. The results are visualized through multiple plots comparing experimental and analytical velocities at different stations along the plate.

Uploaded by

sarradaboub
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 6

Lab 4 - Boundary Layer Plate

NAME: SARAH DABOUB

LAB TITLE: Boundary Layer Plate

DATE: 02/21/2025

clc; % this will clear the fx >> command window


clear; % this clears local variables so there is no conflicts
close all; % close plots so each plot windows is new and correct
format short; % only show 5 significant digits for clear outputs

Input Data and Define Constants


% Distance of longitudional stations from front edge in meters
x = 10^2.*[0.06,0.11,0.16];
% Air Density
density = 0.95;
% Dyanamic Viscosity
mu = 0.0000185;
% Theoretical Velocity
u=22.5;

The loop below imports data from each sheet of excel data and make

calculations of velocity at each pressure taps for each longitudional

stations

for i=1:1:3
data = xlsread('L04_Data Sheet_Boundary Layer Plate.xlsx',i);
% Average stagnation pressure of each taps in pascal
ave_stag(i,1:10) =
-248.84.*(data(1:10,2)+data(1:10,3)+data(1:10,4)+data(1:10,5))./4;
% Average static pressure of flow in pascal
ave_static_flow(1,i) = -248.84.*(data(12,2)+data(12,3)+data(12,4)+data(12,5))/4;
% Average stagnation pressure of flow in pascal
ave_stag_flow(1,i) = -248.84.*(data(13,2)+data(13,3)+data(13,4)+data(13,5))/4;
% Dynamic pressure at each pressure taps
dy_pressure(i,1:10) = ave_stag(i,1:10)-ave_static_flow(1,i);
% Velocity at each pressure taps
vel_tap(i,1:10) = sqrt(2.*dy_pressure(i,1:10)./density);
% Dyanamic pressure of flow
dy_pressure_flow(1,i) = ave_stag_flow(1,i)-ave_static_flow(1,i);
% Velocity of flow
vel_flow(1,i) = sqrt(2.*dy_pressure_flow(1,i)./density);
% Reynolds Number at each station
reynolds(1,i) = density*vel_flow(1,i)*x(i)/mu;
disp(reynolds(1,i))
% Boundary layer thickness δ at each station
if reynolds(1,i)<2000;
delta(1,i)=5.*x(1,i)/sqrt(reynolds(1,i));

1
else;
delta(1,i)=0.37.*x(1,i)/(reynolds(1,i)^(1/5));
end
disp(delta(1,i))
end

8.1938e+06
0.0920
1.4741e+07
0.1499
2.1234e+07
0.2027

Importing the tap-heights from the excel data record sheet


height_data = xlsread('L04_Data Sheet_Boundary Layer Plate.xlsx',4);
tap_height = height_data(:,9);

Importing analytical solution


analy_soln = xlsread('L04_Data Sheet_Boundary Layer Plate.xlsx',5);
% velocity of each taps in each stations from analytical solution
analy_vel_front = analy_soln(1:17,9);
analy_vel_middle = analy_soln(19:35,9);
analy_vel_rear = analy_soln(37:53,9);
% Heights from the plate surface
height_analy_front = analy_soln(1:17,7);
height_analy_middle = analy_soln(19:35,7);
height_analy_rear = analy_soln(37:53,7);

Plotting experimental and analytical results


% Front Station
figure; hold on
plot(vel_tap(1,1:10),tap_height);
plot(analy_vel_front,height_analy_front);
yline(delta(1,1), 'g--')
xline(u,'--')
title('Velocity vs. Height, Front');
xlabel('Velocity(m/s)');
ylabel('Height(cm)');
legend('Experimental','Analytical','Boundary Layer Thickness','Theoretical
Velocity');
hold off;

2
% Middle Station
figure; hold on
plot(vel_tap(2,1:10),tap_height);
plot(analy_vel_middle,height_analy_middle);
yline(delta(1,2), 'g--')
xline(u,'--')
title('Velocity vs. Height, Middle');
xlabel('Velocity(m/s)');
ylabel('Height(cm)');
legend('Experimental','Analytical','Boundary Layer Thickness','Theoretical
Velocity');
hold off;

3
% Rear Station
figure; hold on
plot(vel_tap(3,1:10),tap_height);
plot(analy_vel_rear,height_analy_rear);
yline(delta(1,3), 'g--')
xline(u,'--')
title('Velocity vs. Height, Rear');
xlabel('Velocity(m/s)');
ylabel('Height(cm)');
legend('Experimental','Analytical','Boundary Layer Thickness','Theoretical
Velocity');
hold off;

4
% All together
figure; hold on
plot(vel_tap(1,1:10),tap_height,'DisplayName','Front');hold on
plot(vel_tap(2,1:10),tap_height,'DisplayName','Middle');hold on
plot(vel_tap(3,1:10),tap_height,'DisplayName','Rear');hold on
%plot(analy_vel_rear,height_analy_rear);
title('Velocity vs. Height');
xlabel('Velocity(m/s)');
ylabel('Height(cm)');
legend('Location','Best');
hold off;
grid on

5
6

You might also like