From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 2223/CH18/EX18.55/Ex18_55.sce | 23 +++++++++++++++++++++++
 1 file changed, 23 insertions(+)
 create mode 100755 2223/CH18/EX18.55/Ex18_55.sce

(limited to '2223/CH18/EX18.55/Ex18_55.sce')

diff --git a/2223/CH18/EX18.55/Ex18_55.sce b/2223/CH18/EX18.55/Ex18_55.sce
new file mode 100755
index 000000000..37d57d471
--- /dev/null
+++ b/2223/CH18/EX18.55/Ex18_55.sce
@@ -0,0 +1,23 @@
+// scilab Code Exa 18.55 Propeller Thrust and Power
+
+c_u=5; // upstream velocity in m/s
+c_s=10; // downstream velocity in m/s
+rho=1e3; // density of water in kg/m3
+c=0.5*(c_u+c_s); // velocity of water through the propeller in m/s
+d(1)=0.5; // propeller diameter in m
+d(2)=1;
+d(3)=1.5;
+delh_0=0.5*((c_s^2)-(c_u^2));
+delp_0=rho*delh_0;
+disp("bar",delp_0*1e-5,"(b)stagnation pressure rise across the propeller is")
+for i=1:3
+    disp("cm",d(i)*1e2,"for propeller diameter=")
+A=%pi*(d(i)^2)/4;
+Q=c*A;
+m=rho*Q;
+disp("m3/s",Q,"(a) flow rate through the propeller is")
+Fx=A*delp_0;
+disp("kN",Fx*1e-3,"(c) thrust exerted by the propeller on the boat is")
+P=m*delh_0;
+disp("kW",P/1000,"(d)the ideal Power required to drive the propeller is")
+end
-- 
cgit