clear all
close all

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

d = 8*.0254;  %Pelton wheel diameter.
dj = (3/8)*.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(3,1,1)
plot(V,Nbar)
xlabel('Velocity (m/s)');
ylabel('Specific Speed');
title('Specific Speed versus velocity,  Bore diam = 6" and Pelton wheel diameter = 8 in');

subplot(3,1,2)
plot(V,N*60/(2*pi))
xlabel('Velocity (m/s)');
ylabel('RPM');
title('Pelton wheel rpm versus velocity,  Bore diam = 6" and Pelton wheel diameter = 8 in, jet Diam = 3/8 in');

subplot(3,1,3)
plot(V,N*60/(2*pi))
xlabel('Velocity (m/s)');
ylabel('Force');

title('Force/Velocity relationship,  Bore diam = 6" and Pelton wheel diameter = 8 in, jet Diam = 3/8 in');

