4 files changed, 73 insertions, 0 deletions
diff --git a/914/CH15/EX15.1/ex15_1.sce b/914/CH15/EX15.1/ex15_1.sce
new file mode 100755
index 000000000..9a70958e9
--- /dev/null
+++ b/914/CH15/EX15.1/ex15_1.sce
@@ -0,0 +1,20 @@
+clc;
+warning("off");
+printf("\n\n example15.1 - pg760");
+// given
+r=[10 20 50 100 200 400 600 1000 2000]
+tau=[2.2 3.1 4.4 5.8 7.4 9.8 11.1 13.9 17.0]
+tau=tau*(10^-4);
+clf;
+xtitle("basic shear diagram for the fluid","shear rate","shear stress");
+plot2d("ll",r,tau);
+// the data falls nearly on a straight line
+// from the graph the slope and the intercept are
+slope=0.3841;
+intercept=9.17046;
+// from the relation tau=K*(-r)^n;
+K=exp(intercept);
+n=slope
+disp(K,"K=",n,"n=");
+printf("\n\n The fluid is pseudo plastic , since the slope is less than 1 ");
+
diff --git a/914/CH15/EX15.2/ex15_2.sce b/914/CH15/EX15.2/ex15_2.sce
new file mode 100755
index 000000000..910b07298
--- /dev/null
+++ b/914/CH15/EX15.2/ex15_2.sce
@@ -0,0 +1,16 @@
+clc;
+warning("off");
+printf("\n\n example15.2 - pg774");
+// given
+a=[651 1361 2086 5089 7575 11140 19270 25030]
+tau=[3.71 7.49 11.41 24.08 -35.21 46.25 77.50 96.68]
+clf;
+xtitle("capillary shear diagram for polyisobutylene L-80 in cyclohexane","pseudoshear rate","wall shear stress");
+plot2d("ll",a,tau);
+// from the graph
+betao=-4.3790154;
+beta1=0.8851;
+K'=exp(betao);
+n'=beta1;
+printf("\n\n The final rheological model is \n tauw = %f*(8*Uz,avg/do)^%f",K',n');
+
diff --git a/914/CH15/EX15.3/ex15_3.sce b/914/CH15/EX15.3/ex15_3.sce
new file mode 100755
index 000000000..df1d64221
--- /dev/null
+++ b/914/CH15/EX15.3/ex15_3.sce
@@ -0,0 +1,12 @@
+clc;
+warning("off");
+printf("\n\n example15.3 - pg774");
+// given
+// from example 15.2 
+n'=0.8851;
+K'=0.01254;
+n=n';
+K=K'/((3*n+1)/(4*n));
+disp(n,"n=");
+printf("\n K = %f N/m^2",K);
+
diff --git a/914/CH15/EX15.4/ex15_4.sce b/914/CH15/EX15.4/ex15_4.sce
new file mode 100755
index 000000000..154acfce2
--- /dev/null
+++ b/914/CH15/EX15.4/ex15_4.sce
@@ -0,0 +1,25 @@
+clc;
+warning("off");
+printf("\n\n example15.4 - pg775");
+// given
+a=[10 20 50 100 200 400 600 1000 2000];
+tau=[2.24 3.10 4.35 5.77 7.50 9.13 11.0 13.52 16.40]
+tau=tau*10^-4;
+clf;
+xtitle("capillary shear diagram for a commercial polyethylene melt at 190 degC","pseudoshear rate","wall shear stress");
+plot2d("ll",a,tau);
+// such a plot suggests a second order polynomila  of the type y=betao+beta1*x+beta2*x^2;
+// where y=ln(tauw) and x=ln(8*Uz,avg/do)=ln(a);
+// from the graph
+betao=8.96694;
+beta1=0.48452520;
+beta2=0.010923041;
+n=beta1+2*beta2*a;
+phiw=((3*n+1)/(4*n))*(a);
+mu=tau/phiw;
+printf("\n\n 8*Uz,avg/do       n       (3*n+1)/(4*n)       phiw            mu");
+for i=1:9
+    printf("\n %f       %f      %f      %f       %f",a(i),n(i),(3*n(i)+1)/(4*n(i)),phiw(i),mu);
+end
+
+