diff options
| author | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
|---|---|---|
| committer | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
| commit | 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch) | |
| tree | dbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /3720/CH8/EX8.5/Ex8_5.sce | |
| parent | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff) | |
| download | Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2 Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip | |
initial commit / add all books
Diffstat (limited to '3720/CH8/EX8.5/Ex8_5.sce')
| -rw-r--r-- | 3720/CH8/EX8.5/Ex8_5.sce | 29 |
1 files changed, 29 insertions, 0 deletions
diff --git a/3720/CH8/EX8.5/Ex8_5.sce b/3720/CH8/EX8.5/Ex8_5.sce new file mode 100644 index 000000000..b92bf469a --- /dev/null +++ b/3720/CH8/EX8.5/Ex8_5.sce @@ -0,0 +1,29 @@ +//Example 8_5
+clc;clear;funcprot(0);
+//From Example 8_4
+// Given values
+P=1;// atm
+T=35;// degree celsius
+L=300;// m
+D=0.267;// m
+h_L=20;// m
+v_old=0.35;// m^3/s
+g=9.81;// m/s^2
+//Properties
+rho=1.145;// kg/m^3
+mu=1.895*10^-5;// kg/m.s
+nu=1.655*10^-5;// m^2/s
+
+//Calculation
+//V=y(1); Re=y(2); f=y(3);v=y(4)
+function[X]=flowrate(y);
+ X(1)=real((y(4)/(%pi*D^2/4))-y(1));
+ X(2)=real(((y(1)*D)/(nu))-y(2));
+ X(3)=real((-2.0*log10(2.51/(y(2)*sqrt(y(3)))))-(1/sqrt(y(3))));
+ X(4)=real(((y(3)*L*y(1)^2)/(D*2*9.81))-20);
+endfunction
+y=[1 10000 0.01 0.1];
+z= fsolve(y,flowrate);
+v_new=z(4);// m^3/s
+v_drop=v_old-v_new;//The drop in the flow rate
+printf('The drop in the flow rate through the duct.v_drop=%0.2f m^3/s\n',v_drop);
|