1 files changed, 23 insertions, 0 deletions
diff --git a/1394/CH16/EX16.3.3/Ex16_3_3.sce b/1394/CH16/EX16.3.3/Ex16_3_3.sce
new file mode 100755
index 000000000..9e152b130
--- /dev/null
+++ b/1394/CH16/EX16.3.3/Ex16_3_3.sce
@@ -0,0 +1,23 @@
+
+
+clc
+//initialization of variables
+D =2*10^-6 // cm^2/sec
+nu = 0.036 // cm^2/sec 
+d1 = 1.59 // cm
+d2 = 1 // cm
+deltap = 1*10^-5 // g/cc ( change in density)
+p = 1 // g/cc
+Re = 11200 // Reynolds number
+g = 980 // cm/sec^2 
+dis = 5.37*10^-9 // g/cm^2-sec // Dissolution rate
+sol = 1.48*10^-3 // g/cc
+//Calculations
+k11 = 0.62*(D/d1)*(Re^(0.5))*((nu/D)^(1/3))// cm/sec
+K1 = dis/sol// the overall mass transfer co efficient in cm/sec
+k2 = (1/((1/K1)-(1/k11)))/// cm/sec /// the rate constant in cm/sec
+k12 = (D/d2)*(2+(((0.6*((d2^3)*(deltap)*g/(p*nu^2)))^0.25)*((nu/D)^(1/3)))) // cm/sec
+K2 = 1/((1/k12)+(1/k2))// cm/sec (the overall mtc)
+//Results
+printf("the rate of surface reaction is %.7f cm/sec",k2)
+printf("\nThe dissolution rate for 1 cm gallstone is %.7f cm/sec",K2)