1 files changed, 20 insertions, 0 deletions
diff --git a/914/CH10/EX10.1/ex10_1.sce b/914/CH10/EX10.1/ex10_1.sce
new file mode 100755
index 000000000..76d14fefd
--- /dev/null
+++ b/914/CH10/EX10.1/ex10_1.sce
@@ -0,0 +1,20 @@
+clc;
+warning("off");
+printf("\n\n example10.1 - pg405");
+T=30;  //[degC] - temperature
+d=8.265*10^-4;  //[m] - diameter of the capillary viscometer
+deltapbyL=-0.9364;  //[psi/ft] - pressure drop per unit length
+deltapbyL=deltapbyL*(2.2631*10^4);  //[kg/m^2*sec^2] - pressure drop per unit length
+Q=28.36*(10^-6)*(1/60);
+p=(0.88412-(0.92248*10^-3)*T)*10^3;  //[kg/m^3] - density
+s=(%pi*(d^2))/4;
+U=Q/s;
+tauw=(d/4)*(-deltapbyL);
+shearrate=(8*U)/d;
+mu=tauw/(shearrate);
+printf("\n\n The viscosity is \n mu=%f kg/m*sec=%f cP",mu,mu*10^3);
+printf("\n\n Finally, it is important to check the reynolds number to make sure the above equation applies");
+Nre=(d*U*p)/(mu);
+disp(Nre,"Nre=");
+printf("\n\n The flow is well within the laminar region and therefore the above equation applies");
+