diff options
Diffstat (limited to '1268/CH4')
| -rw-r--r-- | 1268/CH4/EX4.10/4_10.sce | 14 | ||||
| -rw-r--r-- | 1268/CH4/EX4.11/4_11.sce | 18 | ||||
| -rw-r--r-- | 1268/CH4/EX4.12/4_12.sce | 13 | ||||
| -rw-r--r-- | 1268/CH4/EX4.13/4_13.sce | 12 | ||||
| -rw-r--r-- | 1268/CH4/EX4.14/4_14.sce | 10 | ||||
| -rw-r--r-- | 1268/CH4/EX4.15/4_15.sce | 8 | ||||
| -rw-r--r-- | 1268/CH4/EX4.18/4_18.sce | 14 | ||||
| -rw-r--r-- | 1268/CH4/EX4.19/4_19.sce | 12 | ||||
| -rw-r--r-- | 1268/CH4/EX4.20/4_20.sce | 11 | ||||
| -rw-r--r-- | 1268/CH4/EX4.8/4_8.sce | 7 | ||||
| -rw-r--r-- | 1268/CH4/EX4.9/4_9.sce | 11 |
11 files changed, 130 insertions, 0 deletions
diff --git a/1268/CH4/EX4.10/4_10.sce b/1268/CH4/EX4.10/4_10.sce new file mode 100644 index 000000000..6fcb7a563 --- /dev/null +++ b/1268/CH4/EX4.10/4_10.sce @@ -0,0 +1,14 @@ +clc;
+disp("Example 4.10")
+d=0.05 // diameter in m
+density=1000 // density in kg/m^3
+mew= 0.001 // viscosity in kg/ms
+flowrate= 100/60 // in kg/s
+avgvelo=flowrate*4/(%pi*density*d*d)
+
+Re= avgvelo*d*density/mew
+if Re<50000 then
+ f=0.079/(Re^0.25)
+end
+disp(" The friction factor is ")
+disp(f)
diff --git a/1268/CH4/EX4.11/4_11.sce b/1268/CH4/EX4.11/4_11.sce new file mode 100644 index 000000000..8df170c99 --- /dev/null +++ b/1268/CH4/EX4.11/4_11.sce @@ -0,0 +1,18 @@ +clc;
+disp("Example 4.11")
+d=0.1 // diameter in m
+l=25 // length in m
+density=1000 // density in kg/m^3
+delP= 14700 // in N/m^2
+mew= 0.001 // in kg/ms
+ka= d*density*((delP*d)^0.5)/(((2*density*l)^0.5)*mew)
+Re= (ka/0.281)^(8/7)
+if(Re<50000)
+ v1=Re*mew/(d*density)
+ disp(v1)
+end
+if(Re>50000)
+ Re=(ka/0.2145)^(10/9)
+ v1=Re*mew/(d*density)
+ disp(v1)
+end
diff --git a/1268/CH4/EX4.12/4_12.sce b/1268/CH4/EX4.12/4_12.sce new file mode 100644 index 000000000..de3c1e4b5 --- /dev/null +++ b/1268/CH4/EX4.12/4_12.sce @@ -0,0 +1,13 @@ +clc;
+disp("Example 4.12")
+mplus= 5 // laminar sublayer thickness in dimensionless form
+d= 0.05 // diameter in m
+density= 1000 // in kg/m^3
+mu= 0.001 // viscosity in kg/ms
+nu = mu/density;
+U=1 // velocity in m/s
+Re=density*U*d/mew
+f= 0.0791/(Re^0.25)
+m= (mplus)*nu/(U*((f/2)^0.5))
+disp("Laminar sublayer thickness is ")
+disp(m)
diff --git a/1268/CH4/EX4.13/4_13.sce b/1268/CH4/EX4.13/4_13.sce new file mode 100644 index 000000000..ccc49d734 --- /dev/null +++ b/1268/CH4/EX4.13/4_13.sce @@ -0,0 +1,12 @@ +clc;
+disp("Example 4.13")
+head= 5 // in m
+f= 0.0045
+l= 100 // pipe length in m
+d= 0.05 // pipe diameter in m
+//delP=f*density*u*2*u*l/d and delP should also be equal to density*9.8*head
+// equating these 2 we get a relation for u
+u=((head*9.81*d)/(f*2*l))^0.5
+flowrate= %pi*d*d*u/4
+disp("The flow rate is ")
+disp(flowrate," m^3/s")
diff --git a/1268/CH4/EX4.14/4_14.sce b/1268/CH4/EX4.14/4_14.sce new file mode 100644 index 000000000..e87119394 --- /dev/null +++ b/1268/CH4/EX4.14/4_14.sce @@ -0,0 +1,10 @@ +clc;
+disp("Example 4.14")
+h= 5 // in m
+f=0.005
+Q=(18200/3600)*0.001 // flow rate in m^3/s
+l=50 // in m
+//from the formulae used in the last problem as well
+d=(((2*f*l*Q*4*Q*4)/(%pi*%pi*h*9.81))^0.2)
+disp("The diameter is ")
+disp(d)
diff --git a/1268/CH4/EX4.15/4_15.sce b/1268/CH4/EX4.15/4_15.sce new file mode 100644 index 000000000..4cf8e046a --- /dev/null +++ b/1268/CH4/EX4.15/4_15.sce @@ -0,0 +1,8 @@ +clc;
+disp("Example 4.15");
+Re= 5e5;
+f= 0.046/(Re^0.2)
+ratio= 1+(3.75*((f/2)^0.5));
+disp("the ratio is ");
+disp(ratio);
+disp("Note: The value shown in the book is 1.1453.");
diff --git a/1268/CH4/EX4.18/4_18.sce b/1268/CH4/EX4.18/4_18.sce new file mode 100644 index 000000000..2a496cefb --- /dev/null +++ b/1268/CH4/EX4.18/4_18.sce @@ -0,0 +1,14 @@ +clc;
+disp("Example 4.18")
+density=850 // in kg/m^3
+mew= 0.0005 // in kg/ms
+d= 0.0525 // diameter in m
+G= 7620 // in kg/m^2/s
+U=G/density
+Re=800000
+f= 0.0014+(0.125/(Re^0.32))
+v= mew/density
+m= 5*v/(U*((f/2)^0.5))
+disp(m,"Laminar sub layer thickness is")
+tau= f*density*U*U/2
+disp(tau,"Wall shear stress is ")
diff --git a/1268/CH4/EX4.19/4_19.sce b/1268/CH4/EX4.19/4_19.sce new file mode 100644 index 000000000..dd00cbbd2 --- /dev/null +++ b/1268/CH4/EX4.19/4_19.sce @@ -0,0 +1,12 @@ +clc;
+disp("Example 4.19")
+U= 0.5 // in m/s
+l= 0.025 // side dimension
+area=l*l
+perimeter=4*l
+rh=area/perimeter
+dh=4*rh
+v=1e-6
+Re=dh*U/v
+f=0.0791/(Re^0.25)
+disp(f,"Friction factor is ")
diff --git a/1268/CH4/EX4.20/4_20.sce b/1268/CH4/EX4.20/4_20.sce new file mode 100644 index 000000000..e8c870b36 --- /dev/null +++ b/1268/CH4/EX4.20/4_20.sce @@ -0,0 +1,11 @@ +clc;
+disp("Example 4.20")
+d=0.05 // in m
+G= 125 // Massflow rate per crosssection area in kg/m^2/s
+mew= 0.025 // in kg/ms
+Re=d*G/mew
+density=800
+V=G/density
+Vmax= 2*V
+Vgrad= -2*Vmax*2/d
+disp(Vgrad,"Velocity gradient on the wall is ")
diff --git a/1268/CH4/EX4.8/4_8.sce b/1268/CH4/EX4.8/4_8.sce new file mode 100644 index 000000000..56107898a --- /dev/null +++ b/1268/CH4/EX4.8/4_8.sce @@ -0,0 +1,7 @@ +clc;
+disp("Example 4.8")
+Re=100000 // Reynold number
+f=0.079/(Re^0.25) // friction factor according to 1/5th law
+ratio= 1+ 3.75*((f/2)^0.5)
+disp("The ratio of maximum velocity to the average velocity is ")
+disp(ratio)
diff --git a/1268/CH4/EX4.9/4_9.sce b/1268/CH4/EX4.9/4_9.sce new file mode 100644 index 000000000..1bb51ed61 --- /dev/null +++ b/1268/CH4/EX4.9/4_9.sce @@ -0,0 +1,11 @@ +clc;
+disp("Example 4.9")
+// Here the required ratio in terms of n is 2n^2/((n+1)(2*n+1))
+// and the value of this ratio is 0.817
+// solving this we get the following quadratic equation
+// 0.366n^2-2.451n-0.817=0
+y=[0.366 -2.451 -0.817]
+z=roots(y)
+// z is a matrix that has the roots of the equation
+//since we need the positive value of n
+disp(z(1,1))
|