No Title  |POLVER05_0 |1
z(0) = 0.00001
z(f) = 0.003
d(R) / d(z) = 1 
R(0) = 0.00001
rhoc=2.8*10^6
De=1.82*10^-8
DAB=2*10^-8
Sa=530
k1doubleprime=4.42*10^-10
vo=10^-6
Ac=2.03*10^-3
phi=0.5
vis=1.53*10^-8
k1prime=k1doubleprime*Sa
phi1=R*sqrt(k1prime*rhoc/De)
eff=3/phi1
WP=eff*phi1^2
U=vo/Ac
dp=2*R
Reprime=U*dp/((1-phi)*vis)
Sc=vis/DAB
Shprime=(Reprime^0.5)*(Sc^(1/3))
kc=(1-phi)*DAB*Shprime/(phi*2*R)
ac=6*(1-phi)/dp
rhob=rhoc*(1-phi)
effoverall=eff/(1+(eff*k1doubleprime*Sa*rhob)/(kc*ac))
MR=effoverall*k1prime*rhob*R/kc
X=0.998
W = (vo/(effoverall*k1prime))*ln(1/(1-X))
L=W/(Ac*rhob)