clear all
close all

D = 6*.0254;
Dr = 2*.0254;

d = 10*.0254;  %Pelton wheel diameter.
dj = (1/2)*.0254;
Aj = pi*dj^2/4;

A = pi*(D^2-Dr^2)/4;

V = 0:1/100:1;
b = 14000;
g = 9.81;
rho = 700;

dP = b*V/A;
Q = A*V;

Vj = V*A/Aj;

Vr = Vj/2;

N = Vr/d;  %rad/sec


Nbar = N*sqrt(Q)./(dP/rho).^.75;

figure('Position',[520   142   843   956],'PaperOrientation','portrait','PaperPositionMode','auto');
subplot(2,1,1)
plot(V,N)
xlabel('Velocity (m/s)');
ylabel('Specific Speed');
title('Specific Speed versus velocity, assuming Bore diam = 6" and 6000RPM at 1m/s ram speed');

subplot(2,1,2)