Final Solution

Problem I. Biodiesel
//Geometric //Nmin Estimate //Mass transfer coefficient continuous
V=_pi*(T/2)^2*H Nmin*(D/g)^0.50=C20*((T/D)^alpha5)*(vis phase
H=T C/visD)^(1/9)*((denC- Dc=5E-10 //m2/s
V=1/1000 //m3 denD)/denC)^0.25*(tesion/(D^2*denC*g)) Sh=2
D=T/3 ^0.30 Sh=kd*dp/Dc
C20=4.0 kmda=kd*a
//Reactant alpha5=0.88
molC/molD=1/10 Nmrpm=Nmin*60 //rpm //Apparent Rate
massC=molC*850 //kg k2nd=(0.163/60)//L/mol.s
massD=molD*32.4 //kg //Power/volume kapp=k2nd*CD0
volC=massC/denC //m3 PV=P/V
volD=massD/denD //m3 PV=2500//W/m3 //Da
V=volC+volD //m3 Np=1.27 //4-FBT Kdc=5.06E-4
vfracD=volD/V Np=P/(denL*(N^3)*D^5) kmass=kmda*Kdc*V/volC
mfracD=massD/(massD+massC) Re=denL*D^2*N/visL Da=kapp/kmass
Nrpm=N*60
//Properties of liquid
denC=870.2 //kg/m3 //Estimate emulsion size
denD=792 //kg/m3 We=denC*N^2*D^3/tesion
visC=16.93E-3 //Pa.s //Calderbank 1958
visD=0.40E-3 //Pa.s dp/D=0.06*(1+3.75*vfracD)*We^(-0.60)
tesion=2.50E-3 //N/m
g=9.81 //m/s //Mass transfer area
denL=mfracD*denD+(1-mfracD)*denC n=volD/((_pi*dp^3)/6)
visL=mfracD*visD+(1-mfracD)*visC a=(n*_pi*dp^2)/V
CD0=molD/V //kgmol/m3 //mol/L
CC0=molC/V //kgmol/m3 //mol/L //Estimate teq
N*teq=1995.3*((D/T)^(-
2.37))*((We/Re)^(0.97))*(visD/visC)*(Fr^(-
0.66))
Fr=(N^2)*D/g
t=teq/60
Problem I. Biodiesel
<< results >> PV = 2500
V = 0.001 P = 2.5
H = 0.1083852 Np = 1.27
T = 0.1083852 N = 33.52689
D = 0.0361284 Re = 3002.63
molC = 0.0007215643 Nrpm = 2011.614
molD = 0.007215643 We = 18450.67
massC = 0.6133297 dp = 1.259059e-05
massD = 0.2337868 n = 2.824607e+11
volC = 0.0007048146 a = 140669.5
denC = 870.2 teq = 43.29646
volD = 0.0002951854 Fr = 4.139676
denD = 792 t = 0.7216077
vfracD = 0.2951854 Dc = 5e-10
mfracD = 0.2759796 Sh = 2
visC = 0.01693 kd = 7.94244e-05
visD = 0.0004 kmda = 11.17259
tesion = 0.0025 k2nd = 0.002716667
g = 9.81 kapp = 0.0196025
denL = 848.6184 Kdc = 0.000506
visL = 0.01236806 kmass = 0.00802102
CD0 = 7.215643 Da = 2.443891
CC0 = 0.7215643
Nmin = 11.56786
C20 = 4 a. 0.0001 < Da < 1000 have the mass transfer effect
alpha5 = 0.88
Nmrpm = 694.0716
Problem IIa: Mixing effect in homogeneous reaction
//Properties //Meso mixing time for feed dispersion << results >>
denL=1000 //kg/m3 (semi-batch) denL = 1000
visL= 0.001 //Pa.s QB=1/1000 //m3/s visL = 0.001
U=vtip V = 6
//Geometric for 6PTD Dt=0.10*(k^2)/Emax Area = 3.046474
V=Area*Z tmmesd=QB/(U*Dt) Z = 1.96949
V= 6000/1000 //m3 QBm=QB*1000*60 //L/min/ D = 0.590847
T = 1.96949
D=0.30*T vtip=_pi*N*D C = 0.6564967
Z=T k=0.06*vtip^2 PV = 1000
C=T/3 P = 6000
Area=(_pi*T^2)/4 //Local mixing time Np = 1.27
tm1=tmmic+tmmesd N = 4.033271
//Rotational speed Re = 1408016
PV=P/V Nrpm = 241.9963
tr1 = 2
PV=1000 //W/m3 //Damkoehler number of mixing k1 = 0.5
Np=1.27 //pitched blade turbine Da11=tm1/tr1 CB0 = 1
Np=P/(denL*(N^3)*D^5) Da12=tm1/tr2 tr2 = 5
Re=denL*D^2*N/visL k2 = 0.2
Nrpm=N*60 CA0F = 1
tmmic = 0.000546106
//Reaction time Sc = 1000
nu = 1e-06
tr1=1/(k1*CB0) Emax = 40
tr2=1/(k2*CA0F) Dif = 1e-09
k1=0.50 //L/mol.s QB = 0.001
k2=0.20//L/mol.s U = 7.486561
CB0=1//mol/L vtip = 7.486561
CA0F=1//mol/L Dt = 0.028273
k = 3.362916
tmmesd = 0.00472439
//Micro mixing time QBm = 60
tmmic=1/2*logE(Sc)*(nu/Emax)^0.5 tm1 = 0.005270496
nu=visL/denL Da11 = 0.002635248
Sc=1000 Da12 = 0.001054099
Sc=nu/Dif
Emax/(P/(denL*V))=40
ยังมีผลของ Mixing effect ทั้ง 2 reaction
 << results >>
Problem IIb: Mixing effect in homogeneous reaction << results >>
denL = 1000
//Properties //Meso mixing time for feed dispersion visL = 0.001
denL=1000 //kg/m3 (semi-batch) V = 6
visL= 0.001 //Pa.s QB=1/1000 //m3/s Area = 3.046474
U=vtip Z = 1.96949
D = 0.590847
//Geometric for 6PTD Dt=0.10*(k^2)/Emax T = 1.96949
V=Area*Z tmmesd=QB/(U*Dt) C = 0.6564967
V= 6000/1000 //m3 QBm=QB*1000*60 //L/min/ PV = 1000
D=0.30*T vtip=_pi*N*D P = 6000
Z=T k=0.06*vtip^2 Np = 1.27
C=T/3 N = 4.033271
Area=(_pi*T^2)/4 //Local mixing time Re = 1408016
Nrpm = 241.9963
tm1=tmmic+tmmesd tr1 = 2
//Rotational speed k1 = 0.5
PV=P/V CB0 = 1
PV=1000 //W/m3 //Damkoehler number of mixing tr2 = 5
Np=1.27 //pitched blade turbine Da11=tm1/tr1 k2 = 0.2
Np=P/(denL*(N^3)*D^5) Da12=tm1/tr2 CA0F = 1
Re=denL*D^2*N/visL tmmic = 0.000546106
Sc = 1000
Nrpm=N*60 nu = 1e-06
Emax = 40
//Reaction time Dif = 1e-09
tr1=1/(k1*CB0) QBreq = 0.001
tr2=1/(k2*CA0F) U = 7.486561
k1=0.50 //L/mol.s vtip = 7.486561
k2=0.20//L/mol.s Dt = 0.028273
k = 3.362916
CB0=1//mol/L tmmesd = 0.001453894
CA0F=1//mol/L QB = 0.0003077422
QBmreq = 60
//Micro mixing time QBm = 18.46453
tmmic=1/2*logE(Sc)*(nu/Emax)^0.5 Number = 3.249473
nu=visL/denL tm1 = 0.002
Sc=1000 Da11 = 0.001
Da12 = 0.0004
Sc=nu/Dif
Emax/(P/(denL*V))=40
Number of outlet port of A > 4