blob: 9e152b1306ec1a3bde771288b850a13ecc1ece2a (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
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)
|
