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

---
 51/CH3/EX3.12/3_12.sce      |  21 +++++++++++++++++++++
 51/CH3/EX3.12/3_12graph.jpg | Bin 0 -> 21874 bytes
 2 files changed, 21 insertions(+)
 create mode 100755 51/CH3/EX3.12/3_12.sce
 create mode 100755 51/CH3/EX3.12/3_12graph.jpg

(limited to '51/CH3/EX3.12')

diff --git a/51/CH3/EX3.12/3_12.sce b/51/CH3/EX3.12/3_12.sce
new file mode 100755
index 000000000..da4255811
--- /dev/null
+++ b/51/CH3/EX3.12/3_12.sce
@@ -0,0 +1,21 @@
+clc;
+clear;
+z1=5;//m
+a=0.8;//m
+b=6;//m
+Cc=0.61;//since a/z1=ratio=0.16<0.2; Cc= contracction coefficient
+z2=Cc*a;
+//Q/b=flowrate
+flowrate=z2*((2*9.81*(z1-z2))/(1-((z2/z1)^2)))^0.5;
+//considering z1>>z2 and neglecting kinetic energy of the upstream fluid
+flowrate1=z2*(2*9.81*z1)^0.5;
+disp("m^2/s",flowrate,"The flowrate per unit width=")
+disp("m^2/s",flowrate1,"The flowrate per unit width when we consider z1>>z2=")
+count=1;
+j=5:15;
+for i=5:15
+    fr(count)=z2*((2*9.81*(i-z2))/(1-((z2/i)^2)))^0.5;
+    count=count+1;
+end
+plot2d(j,fr,rect=[0,0,15,9])
+xtitle("Q/b vs z1","z1,m","Q/b, m^2/s")
\ No newline at end of file
diff --git a/51/CH3/EX3.12/3_12graph.jpg b/51/CH3/EX3.12/3_12graph.jpg
new file mode 100755
index 000000000..900df5d77
Binary files /dev/null and b/51/CH3/EX3.12/3_12graph.jpg differ
-- 
cgit