Advanced Oxidation Processes - Semantic Scholar

7 downloads 0 Views 1MB Size Report
Simulation Results (LPUV/H. 2. O. 2. ) – NAIX + Dealkalization -. 0.00. 0.50. 1.00. 1.50. 2.00. 2.50. 3.00. 3.50. 4.00. 0.0. 2.0. 4.0. 6.0. 8.0. 10.0. 12.0. 14.0. 16.0. 5.
Rational Design of Advanced Oxidation Processes using Computational Chemistry AOT 16, November 15-18, 2010 San Diego, CA.

Daisuke Minakata1, Ke Li2, John C. Crittenden1, David Hand3 School of Civil and Environmental Engineering Georgia Institute of Technology 2 Faculty of Engineering, University of Georgia 3 Civil and Environmental Engineering Michigan Technological University 1

1

Modeling AOPs 1. 2. 3. 4.

Simplified Pseudo-Steady-State Model AdoxTM – AOP simulation software Conventional Approach for Kinetic Model Computer-based Kinetic Model

2

UV/H2O2 H2O2/O3 Simplified pseudo-steady-state model

kR  k12  HO

ss, 0

1. UV/H2O2

2H2O2I0 fH2O2 (1- e ) -A

[HOg]ss,0 

 3 0

k10 [H2O2 ]0  k11[HCO ]  k12 [R]0  k13 [NOM]0

2. H2O2 and O3 are added simultaneously

 HO ss, 0 



K L a PO3 H O3



k9  HO2   k10  H 2O2 0  k11  HCO3   k12  R 0  k13  NOM 0  0  0

3. H2O2 is added to a water containing O3, i.e. [O3]0 is known.

pH pK H2O2   2k1  H 2O2 0  10 O3 res  HO ss, 0  k11  HCO3   k12  R 0  k13  NOM 0  0

3

Flow Models: Ideal flow at steady state (1st order reaction) 1. Plug Flow Reactor (PFR)

C  C0e

k

Q  C0  V  ln   k C 

2. Completely Mixed Flow Reactor (CMFR)

C

C0 (1  k)

V

Q  C0  C  r

3. Tank-in-series (TIS)

C

C0 n

   1   nk  

 C 1/ n   nQ  V   0   1    C    k 

4

Flow Models: Non-ideal flow

4. Dispersed flow model closed system

 1 vL  4a exp   CA  2 E  C Ao (1  a)2 exp  a vL   (1  a)2 exp   a vL  2 E   2 E     

a  1 4k(E / vL) 5

Known reaction rate constants Halflife, min Compound

k HO•

-9

[HO•]=10 M

-1 -1 M s

MtBE

1.6×109

Oxalic acid

1.4×10

Acetate ion

7

7×10

Trichloromethane

5.0×106

1,1,2-Trichloroethane Chloroform Chloroacetic cid Glycolic acid 1,1,1-Trichloroethane Benzene Phenol

6

8

1.1×10 6

5×10

7

4.3×10 8

6×10

7

-10

[HO•]=10

M

-11

[HO•]=10

0.01

0.1

1

8

83

825

0.2

2

17

2

23

231

0.11

1

11

2

23

231

0.3

2.7

27

0.02

0.2

2

0.3

3

29

9

0.001

0.01

0.1

9

0.002

0.02

0.2

4×10

7.8×10 6.6×10

kHCO3- = 8.5×106 M-1s-1 kCO3-- = 3.9×108 M-1s-1 H2CO3* ↔ H+ + HCO3(pKa=6.3) HCO3- ↔ H+ + CO32(pKa=10.3)

M

kNOM =(3.0 ~ 4.5)×108

M-1s-1

7

Hypothesis -GCM 



k = Ae



Ea RT

Essence is same as the GCM for the gaseous phase Observed experimental reaction rate constant for a given organic compound is the combined rate of all elementary reactions involving HO•, which can be estimated using Arrhenius kinetic expression. The Ea consists of two parts: (1) Base part from main reaction mechanisms (i.e., H-atom abstraction, HO• addition to alkenes and aromatic compounds and HO• interaction with S, N, or P-atom-containing compounds). (2) Functional group contribution partially from neighboring and/or next nearest neighboring functional group. Minakata, Li, Westerhoff, Crittenden, 2009 ES&T 43, 6220-6227

8

GCM Results -overall 2/21.0E+11

Ideal fit

1.0E+10

kexp/kcal = 0.5 and 2.0

H-atom abstraction (Calibrated)

k cal (M-1s-1)

1.0E+09

HO addition to alkene (Calibrated)

1.0E+08

HO interaction with N,S,P-atomcontaining compounds (Calibrated) HO addition to aromatic compounds (Calibrated)

1.0E+07

H-atom abstraction (Prediction) HO interaction with N, S, P-atomcontaining compounds (Prediction)

1.0E+06

HO addition to aromatic compounds (Prediction)

1.0E+05 1.0E+05 1.0E+06 1.0E+07 1.0E+08 1.0E+09 1.0E+10 1.0E+11 k exp (M-1s-1)

• Over or under prediction for the oxygenated multifunctional group compounds due to solvation effects, and for the halogenated and carboxylic compounds due to larger steric hindrance Minakata, Li, Westerhoff, Crittenden, 2009 ES&T 43, 6220-6227 9

Challenges still remains GCM is limited because: 1. Many datasets are required to calibrate parameters 2. Only deal with molecules for which all required group rate constants and group contribution factors have been calibrated before. 3. Limited number of literature-reported experimental rate constants for other reaction mechanisms than HO• 4. Assumed that a functional group has approximately the same interaction properties under a given molecule, so that the GCM disregards the changes of the functional properties that can arise from the intramolecular environment by electronic push-pull effects, hydrogen bond formation, or by steric effects.

Application of computational chemistry (quantum mechanical calculations) for predicting reaction rate constants in aqueous phase Looking at reaction energies (transition states) 10

Modeling AOPs 1. 2. 3. 4.

Simplified Pseudo-Steady-State Model AdoxTM – AOP simulation software Conventional Approach for Kinetic Model Computer-based Kinetic Model

11

AdOxTM, AOPs Simulation Software • AdOxTM is a part of the Environmental Technologies Design Option Tools (ETDOTTM). It contains mechanism-based models that can be used to evaluate and design advanced oxidation processes. • The final version of AdOx will include the following models: – H2O2/UV model (available in version 1.0) – H2O2/O3 model – UV/O3 model • A comprehensive database is also included in AdOxTM to provide powerful support to the above models as well as ease-of-use features for the user. 12

Crittenden et al., 1999. Wat. Res. 33, 2315-2328

Features of Current Version AdOxTM 1.0 for UV/H2O2 H2O2 + hν → 2HO•

H2O2 = quantum yield of H2O2 (=0.5)

I0 = incident light intensity, einstein cm-2 sec-1 A=2.303b(εH2O2CH2O2+εRCR + εSCS + εHO2-CHO2-) b=pathlength, cm fH2O2= 2.303 b (εH2O2CH2O2 + εHO2-CHO2-)/A H2O2/HO2- + HO  H2O/OH- + HO2 H2O2 + HO2/O2-  HO + H2O/OH- + O2

2.7×107, 7.5×109 3.0, 0.13

HO + HO  H2O2 HO + HO2/O2-  H2O/OH- + O2 HO2 + HO2/O2-  H2O2/HO2- + O2

5.5×109 6.6×109, 7.0×109 8.3×105, 9.7×107

R + hv  Products R + HO  Products

kMtBE=1.6×109, ktBA=6.0×108

HO  + CO32-/HCO3-  CO3-  + OH-/H2O HO + NOM  Products NOM + hv  Products

3.9×108, 8.5×106 2.4×104 (mgC/L)-1s-1 14

Features of Current Version AdOxTM 1.0 • Reactor Options – Completely mixed reactor (CMBR) – Completely flow type reactor (CMFR) w/o TIS – Plug flow reactor – Non-ideal Mixing • Main Process Operational Variables : – Initial hydrogen peroxide concentration – Incident UV-light intensity • Water quality parameters: – TOC concentration – pH – Concentrations of target compounds – Byproducts can be included if they are known 15

Concern about residual H2O2 • H2O2 has low UV light absorption • Need to increase H2O2 concentration to lower EE/O • High residual H2O2 will require treatment to reduce H2O2 such as chlorine

16

Simulation Results (LPUV/H2O2) – NAIX + Dealkalization MtBE

tBA

16.0

4.00

14.0

3.50

12.0

3.00

10.0

2.50

8.0

2.00

6.0

1.50

4.0

1.00

2.0

0.50

0.0

0.00 5

10

15

EEO (kWh/kgal-order)

H2 O2 residual conc. (mg/L)

Residual of H2O2

20

H2 O2 Dosage (mg/L)

• 10 mg/L of H2O2 dosage would be the better choice although the optimum dosage of H2O2 is over 20 mg/L. This is because the cost required for over 20 mg/L of H2O2 dose overweigh the energy cost. • At 10 mg/H2O2 dosage, UV dose was calculated as 2,100 mJ/cm2. Li; Hokanson; Crittenden; Trussell; Minakata. Wat. Res. 2008, 43, 5045-5053

17

References that were used for model validation 

 



Glaze, W. H.; Lay, Y.; Kang J. W. Advanced oxidation processes. A kinetic model for the oxidation of 1,2-dibromo-3 chloropropane in water by the combination of hydrogen peroxide and UV radiation. Ind. Eng. Chem. Res. 1995, 34, 2314-2323. Li, Ke.; Stefan, M.I.; Crittenden, J.C. UV photolysis of trichloroethylene: Product study and kinetic modeling. Environ. Sci. & Technol. 2004, 38, 6685-6693. Li, Ke.; Stefan, M.I.; Crittenden, J.C. Trichloroethylene degradation by UV/H2O2 advanced oxidation process: Product study and kinetic modeling. Environ. Sci. & Technol. 2007, 41(5), 1696-1703. Li, Ke.; Hokanson, D.R.; Crittenden, J.C.; Trussell, R.R.; Minakata, D. Evaluating UV/H2O2 processes for methyl tertbutyl ether and tertiary butyl alcohol removal: Effect of pretreatment options and light sources 18

AdOxTM 2 Simulation Software for O3, O3/H2O2, and mitigation of bromate and THM formation

19

Concern about bromate • Formation of bromate (BrO3-) during ozonation in the presence of bromide ion (Br-) • A nationwide survey of Br- in drinking water sources: approximately 80 μg/L (Amy et al., 1994) • Br- in costal area is expected higher • 10 μg/L of BrO3- standard of MCL associated with cancer risk, Stage 1 of the Disinfectant/Disinfection By-Product (D/DBP) Rule (EPA, 1998) • When ozone is applied to disinfection, tradeoffs between inactivation of cryptosporidium and bromate formation should be evaluated. 20

Mechanisms of bromate formation Simplified reaction scheme for bromate formation during ozonation Disproportionation

Haag and Hoigné, 1983 von Gunten and Hoigné, 1994

pKa = 8.8 Ozone involving pathway

k (M-1 s-1)

HO• involving pathway

O3 + Br - → OBr- + O2 OBr- + O3 → 2O2 + Br OBr- + O3 → O2 + BrO2 HOBr + O3 → O2 + BrO2 - + H+ BrO2 - + O3 → BrO3 - + O2 O3 + Br• → BrO• + O2

160 330 100

HO• + HOBr → BrO• + H2O