From 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 Mon Sep 17 00:00:00 2001
From: prashantsinalkar
Date: Tue, 10 Oct 2017 12:27:19 +0530
Subject: initial commit / add all books

---
 1268/CH7/EX7.2/7_2.sce |  9 +++++++++
 1268/CH7/EX7.3/7_3.sce | 10 ++++++++++
 1268/CH7/EX7.4/7_4.sce | 14 ++++++++++++++
 1268/CH7/EX7.5/7_5.sce | 17 +++++++++++++++++
 4 files changed, 50 insertions(+)
 create mode 100644 1268/CH7/EX7.2/7_2.sce
 create mode 100644 1268/CH7/EX7.3/7_3.sce
 create mode 100644 1268/CH7/EX7.4/7_4.sce
 create mode 100644 1268/CH7/EX7.5/7_5.sce

(limited to '1268/CH7')

diff --git a/1268/CH7/EX7.2/7_2.sce b/1268/CH7/EX7.2/7_2.sce
new file mode 100644
index 000000000..47364e316
--- /dev/null
+++ b/1268/CH7/EX7.2/7_2.sce
@@ -0,0 +1,9 @@
+clc;
+disp("Example 7.2")
+// The formula derived in the first part is
+//de=2*phi*dp*e/(3*(1-e))
+// and here de=dp/2 and we define x=de/dp
+x=1/2
+phi=1
+e=(3/(4*phi))/(1+(3/(4*phi)))
+disp(e)
diff --git a/1268/CH7/EX7.3/7_3.sce b/1268/CH7/EX7.3/7_3.sce
new file mode 100644
index 000000000..9bde77958
--- /dev/null
+++ b/1268/CH7/EX7.3/7_3.sce
@@ -0,0 +1,10 @@
+clc;
+disp("Example 7.3")
+massflux=54000/3600   // in kg/M^2/s
+d= 0.0022 // in m
+e= 0.32
+density= 1000 // in kg/m^3
+mew=0.001 // in kg/ms
+P=(massflux*massflux*(1-e)/(e*e*e*d*density))*(1.75+(150*mew*(1-e)/(d*massflux)))
+disp(P,"Frictional pressure drop is ")
+disp("The value shown in the book is 10278.1, which is different from the above, due to rounding off errors");
diff --git a/1268/CH7/EX7.4/7_4.sce b/1268/CH7/EX7.4/7_4.sce
new file mode 100644
index 000000000..2594534ff
--- /dev/null
+++ b/1268/CH7/EX7.4/7_4.sce
@@ -0,0 +1,14 @@
+clear;
+disp("Example 7.4")
+d=0.0025 // in m
+G=50000/3600  // in kg/m^2/s
+mew=0.1  // in kg/ms
+Re= d*G/mew
+density=1250
+U=G/density
+P=250000
+y=P*d*d/(150*U*mew)
+//we need to solve equation y*x^3+x^2-1=0
+p=[y,-1,2,-1]
+z=roots(p)
+disp(z(3,1),"Porosity is ")
diff --git a/1268/CH7/EX7.5/7_5.sce b/1268/CH7/EX7.5/7_5.sce
new file mode 100644
index 000000000..93c2c0642
--- /dev/null
+++ b/1268/CH7/EX7.5/7_5.sce
@@ -0,0 +1,17 @@
+clc;
+disp("Example 7.5")
+delD= 250 // difference of the 2 densities in kg/m^3
+density= 1000 // of water in kg/m^3
+mew= 0.001 // in kg/ms
+d=1.25e-3 // in m
+e= 0.4 // porosity
+g=9.81
+phi=1
+// we get an equation for velocity
+// a*u^2+b*u+c=0 where we calculate a,b and c as----
+c=-g*delD
+b= (150*mew*(1-e))/(d*d*e*e*e*phi*phi)
+a=1.75*density/(phi*d*e*e*e)
+p=[a,b,c]
+z=roots(p)
+disp(z(2,1),"Fluidization velocity is ")
-- 
cgit