close all
clear all

nu = 1.3e-6;  %m^2/s

diam = 10*.0254;  %10 Inch Diameter Pipe

v_mean = .01:.01:.2;  %Velocity range

Q = v_mean*pi*.25*diam^2;  %m^3/sec

Re = v_mean*diam/nu;

El_lam = diam*0.06*Re;
El_turb = diam*4.4*Re.^(1/6);


figure('Position', [520   129   853   969],'PaperOrientation','portrait','PaperPosition',[0 0 8.5 11])
ah1 = subplot(2,1,1)
plot(100*v_mean,El_turb);
set(gca,'fontsize',18)
xlabel('Mean Velocity (cm/s)')
ylabel('Entrance Length (m)');
title('Turbulent Entrance Lengths for 10" Diameter Pipe Flow (seawater)');
ah2 = subplot(2,1,2)
plot(100*v_mean,Re);
set(gca,'fontsize',18)
xlabel('Mean Velocity (cm/s)')
ylabel('Reynolds Number (seawater)');
print('-dpdf','10In_EntranceLength');
print('-dpng','10In_EntranceLength');


