diff options
author | priyanka | 2015-06-24 15:03:17 +0530 |
---|---|---|
committer | priyanka | 2015-06-24 15:03:17 +0530 |
commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
tree | ab291cffc65280e58ac82470ba63fbcca7805165 /1394/CH16/EX16.3.3/Ex16_3_3.sce | |
download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip |
initial commit / add all books
Diffstat (limited to '1394/CH16/EX16.3.3/Ex16_3_3.sce')
-rwxr-xr-x | 1394/CH16/EX16.3.3/Ex16_3_3.sce | 23 |
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)
|