Temat:
Zaprojektować absorber służący do absorpcji CO2 (składnik A) z mieszaniny gazowej zawierającej CO i powietrze. Jako absorbent zastosować wodę.
Składnik A: CO2, stężenie: 20 [% obj]
Składnik B: CO, stężenie: 20 [% obj]
Składnik C: powietrze, stężenie: 60 [% obj]
STOPIEŃ ABSORPCJI = 98%
Przepływ V= 0,6 [m3/h]
Temperatura T = 291 K
Ciśnienie P= 7,5∙106 [Pa]
DANE | OBLICZENIA | WYNIKI |
---|---|---|
Mco2 = 40[g/mol] Mco=28[g/mol]
φ =0,98 |
Przeliczono procenty objętościowe fazy gazowej y, na ułamki wagowe yw, yw1,co2 = $\frac{y_{co2}*M_{co2}}{y_{co2}*M_{co2} + y_{\text{co}}*M_{\text{co}} + y_{\text{pow}}*M_{\text{pow}}} =$ $\frac{0,2*40}{0,2*40 + 0,2*28 + 0,6*29} = 0,258$ yw1,co = 0,181 yw1,pow =0,561 Y1,co2 = $\frac{y_{w1,co2}}{1 - y_{w1,co2}}$ = $\frac{0,258}{1 - 0,258}$ = 0,348 $\left\lbrack \frac{\text{kg}\text{CO}_{2}}{\text{kgINERTU}} \right\rbrack$ Y2,co2 = (1-φ)* Y1,co2 = (1-0,98)*0,348 = 0,00696 $\left\lbrack \frac{\text{kg}\text{CO}_{2}}{\text{kgINERTU}} \right\rbrack$ yw2,co2 = $\frac{Y_{2,co2}}{1 + Y_{2,co2}}$ = $\frac{0,00696}{1 + 0,00696}$ = 0,00691 yw2,co = $\frac{y_{w1,co}*(1 + Y_{1,co2})}{1 + Y_{2,co2}}$ = 0,242 yw2,pow = 1 - 0,0649 - 0,182 = 0,751 |
Y1,co2 = 0,348 $\left\lbrack \frac{\text{kg}\text{CO}_{2}}{\text{kgINERTU}} \right\rbrack$ Y2,co2 = 0,00696 $\left\lbrack \frac{\text{kg}\text{CO}_{2}}{\text{kgINERTU}} \right\rbrack$ yw2,co2 = 0,00691 yw2,co = 0,242 yw2,pow = 0,751 |
DANE | OBLICZENIA | WYNIKI |
yw2,co2 = 0,0069 yw2,co = 0,242 yw2,pow = 0,751 Mco2 = 40[g/mol] Mco=28[g/mol]
Stała Henry’ego z danych literaturowych H=1,24*105 [Pa] P=7,5*106[Pa] y1,co2 = 0,347 x*1,co2 =1,61*10-2 Mco2 = 40[g/mol]
Y1,co2 = 0,348 Y2,co2 = 0,00695
|
y2,co2 = $\frac{\frac{y_{w2,co2}}{M_{co2}}}{\frac{y_{w2,co2}}{M_{co2}} + \frac{y_{w2,co}}{M_{\text{co}}} + \frac{y_{w2,pow}}{M_{\text{pow}}}} = \ \frac{\frac{0,0069}{40}}{\frac{0,0069}{40} + \frac{0,242}{28} + \frac{0,751}{29}}$ = 0,00497 y2,co = 0,248 y2,pow = 0,746 Obliczono stężenia równowagowe w cieczy chłonnej (woda pozbawiona CO2) na wylocie:
Minimalne zużycie wody:
Założono, że zużycie wody będzie 1,3 razy większe od minimalnego:
|
y2,co2 =0,00497 y2,co = 0,248 y2,pow = 0,746 x*1,co2 = 1,61*10-2
|
DANE | OBLICZENIA | WYNIKI |
Y1,co2 = 0,348 Y2,co2 = 0,00696
Mco2 = 40[g/mol]
|
Jako wypełnienie dobrano pierścienie Raschiga ceramiczne uporządkowane o następujących parametrach: średnica: dwyp = 38,1 [mm] wolna objętość: ε = 0,68 powierzchnia: a = 115 [m2/m3] |
|
DANE | OBLICZENIA | WYNIKI |
Mco2 = 40[g/mol] Mco=28[g/mol]
y2,co2 =0,00497 y2,co = 0,248 y2,pow = 0,746
|
Obliczono średnicę absorbera w oparciu o korelację Kafarowa i Dytnierskiego: Obliczamy średnią masę molową mieszaniny między wlotem a wylotem z aparatu: na wylocie: $M_{m} = 0,00497*40 + 28*0,248 + 0,746 \bullet 29 = 28,80\lbrack\frac{g}{\text{mol}}\rbrack$ na wlocie: $M_{m} = 0,2 \bullet 40 + 0,2 \bullet 28 + 0,6 \bullet 29 = 31\ \lbrack\frac{g}{\text{mol}}\rbrack$
|
$$M_{m,sr} = 29,9\ \lbrack\frac{g}{\text{mol}}\rbrack$$ |
DANE | OBLICZENIA | WYNIKI |
P=7,5∙106[Pa]
R=8,314[J/K•mol] T=291K
ρg=0,927[kg/m3] ρc=999[kg/m3] y=0,102 ε = 0,68 g=9,81[m/s2] ρg=0,927[kg/m3] ρc=999[kg/m3] a = 115 [m2/m3] wogkr=0,06[m/s] $\dot{V}$0=0,6[m3/s] T=291 T0=273,15K p=7,5∙106[Pa] p0=101325[Pa] |
Przeliczono natężenie przepływu surowca na warunki pracy aparatu: | ρg=0,927[kg/m3] $\frac{\eta_{C}}{\eta_{W}}$=1 x=0,62 wogkr=0,06[m/s] wrz=0,042[m/s] $\dot{V}$=0,00864[m3/s] |
DANE | OBLICZENIA | WYNIKI |
$\dot{V}$=0,6[m3/s] wrz=0,042[m/s] |
Średnica aparatu: Obliczam współczynnik wnikania masy w fazie gazowej: Lepkość mieszaniny gazowej wyznaczono za pomocą metody Herninga i Zipperera oraz korzystając z wykresu Uyehary i Watsona. Na wlocie dla warunków normalnych: Na wylocie dla warunków normalnych: |
Dobl=0,5 [m] ηm=1,59*10-5[Pa*s] ηm=1,69*10-5[Pa*s] |
DANE | OBLICZENIA | WYNIKI |
Średnia lepkość dla warunków normalnych: Parametry pseudo-zredukowane mieszaniny gazowej w warunkach normalnych: Dla wlotu: Dla wylotu: |
ηmsr=1,58*10-5[Pa*s] Tkm=220,5[K] Pkm=4,44[MPa] |
|
DANE | OBLICZENIA | WYNIKI |
T=273[K] Tkm=133,5[K] Tkm=220,5[K] p0=101325[Pa] Pkm=3,72[MPa] Pkm=4,44[MPa] ηmsr=1,64*10-5[Pa*s] ηr=0,56 T=288K Tkm=133,5[K] Tkm=220,5[K] p=7,5*106[Pa] pkm=3,72[MPa] pkm=4,44[MPa] ηr=0,62 ηk=2,93*10-5[Pa*s] |
Wyznaczono lepkość krytyczną mieszaniny gazowej korzystając z wykresu: Parametry zredukowane mieszaniny gazowej w warunkach pracy absorbera: Dla tych wartości parametrów zredukowanych z wykresu odczytano: ηg=0,62 . =1, 82 * 10−5 [Pa.s] |
Trm0=1,55[K] Prm0=0,0248 [Pa] ηk=2,94*10-5[Pa*s] Trm=1,62[K] prm=0,245[Pa] ηg=1,82*10-5[Pa*s] |
DANE | OBLICZENIA | WYNIKI |
p0=101325[Pa] R=8,314[J/K•mol] T=298[K] δ=6,2*10-4 C=40,9[mol/m3] ρg=0,927[kg/m3] $\dot{V}$=0,00864[m3/s] Dobl=0,5[m] |
Obliczam współczynnik dyfuzji CO2 przez mieszaninę pod ciśnieniem atmosferycznym (mieszaninę traktuję jak powietrze): Współczynnik dyfuzji dla T=25 oC i P=atm: Współczynnik dyfuzji dla parametrów pracy absorbera: D = 1,48 * |
C=40,9[mol/m3] D=1,52*105[m2/s] D=1,48*105[m2/s] wg=0,389[kg/m2s] |
DANE | OBLICZENIA | WYNIKI |
wg=0,389[kg/m2s] ηg=1,82*10-5[Pa*s] a = 115 [m2/m3] ηg=1,82*10-5[Pa*s] ρg=9,269[kg/m3] Dg=1,48*105[m2/s] Rez=18,6 Sc=1,242 D=1,48*105[m2/s] Shg=84,7 dwyp=0,0381[m] |
Obliczono współczynnik wnikania masy w fazie ciekłej: | Rez=18,6 Sc=1,242 Shg=1,223 βg=4,6*10-3[m/s] Cog=41,76[mol/m3] Βg’= [mol/m2s] |
DANE | OBLICZENIA | WYNIKI |
wg=0,0389[kg/m2s] wc=0,185[kg/m2s] ηH2O=0,00115[Pa*s] a = 115 [m2/m3] ηH2O=0,00115[Pa*s] g=9,81[m/s2] ρc=999[kg/m3] T=291[K] ηH2O=0,00115[Pa*s] T=291[K] F =1,7*1014 [kg/N] ηH2O=0,00115[Pa*s] |
wc=4,75.0,0389=0,185[kg/m2s] Założono, że lepkość nie jest funkcją ciśnienia (absorpcja prowadzona jest pod ciśnieniem atmosferycznym): Zastępczy wymiar liniowy: Aby obliczyć Dc w warunkach pracy absorbera korzystamy z wykresu Wilkiego: Vc dla wody wynosi 34*10-6 Φ dla wody wynosi 1 F = $\frac{T}{D*\eta}$ = 1,7*1014 [kg/N] Dc = $\frac{T}{\eta*F} = \ \frac{291}{0,00115*1,7*10^{14}\ } = 1,49*10^{- 9}\ $ |
wc=0,185[kg/m2s] Rez=1,398 ϑz=55,3*105[m] Vc=34*10-6 Φ=1 F =1,7*1014 [kg/N] Dc=1, 49 * 10−9 |
DANE | OBLICZENIA | WYNIKI |
ηH2O=0,00115[Pa*s] ρc=999[kg/m3] Dc=1, 47 * 10−9 Rez=1,39 Sc=783 Shc=16,8 Dc=1, 47 * 10−9 ϑz=55,3*10-6[m] ρc=999[kg/m3]
βc=4,41*10-6[m/s] Co=55500[mol/m3] H=1,24*105 [Pa] p=7,5*105[Pa] βc’=24,5[mol/m2s] Βg’=0,019[mol/m2s] m=12,4 |
Współczynnik przenikania masy kY: | Sc=783 Shc=16,8 βc=4,41*10-6[m/s] Co=55500[mol/m3] βc’=24,5[mol/m2s] m=12,4 kY=0,0098[mol/m2s] |
DANE | OBLICZENIA | WYNIKI |
Y2,co2 = 0,00695
Y1,co2 = 0,347 Y1*=0,667 ΔY1=0,00695 ΔY2=0,79 |
Ilość masy wymienianej w kolumnie: Siła napędowa procesu: Natężenie przepływu inertu: |
y1*=0,400 Y1*=0,600 ΔY1=0,00695 ΔY2=0,79 ΔYśr=0,0626 |
DANE | OBLICZENIA | WYNIKI |
$\dot{V}$=0,6[m3/s]
$\dot{V}$inertu=0,0069[m3/s] ρg=92,69[kg/m3] Ginertu=0,0064[kg/s] ΔY1=0,00695 Y1,co2 = 0,347 ṁco2=0,0022[kg/s] ṅCO2=0,546[mol/s] ΔYśr=0,0626 A=88,38[m2] Arz=115 [m2] a = 115 [m2/m3] d=0,5[m] H=4,86[m] D=0,5[m] |
Ilość wymienionej masy CO2: Powierzchnia wymiany masy: Zwiększono powierzchnię o 30 %: Wysokość absorbera: Sprawdzenie warunku smukłości kolumny:
|
$\dot{V}$inertu=0,0069[m3/s] Ginertu=0,0064[kg/s] ṁCO2=0,0022[kg/s] ṅCO2=0,0546[mol/s] A=88,38[m2] Arz=115[m2] H=4,86[m] S=9,48m |
Zestawienie wyników:
średnica absorbera : D = 0,5 [m]
powierzchnia wymiany masy: A=115 [m2]
wysokość absorbera: H=4,86 [m]
LITERATURA
[1] „Tablice do obliczeń procesowych” Pod redakcją Leona Troniewskiego, wydanie V poprawione i uzupełnione.
[2] „Absorpcja i absorbery” Roman Zarzycki, Andrzej Chacuk, Maciej Starzak, wydanie drugie uaktualnione
[3] Zadanie rachunkowe z inżynierii chemicznej” Praca zbiorowa pod redakcja Romana Zarzyckiego, Warszawa 1980