Article
Removal of Oil from Water by Inverse Fluidization of Aerogels
Current address: Shell Global Solutions, P.O. Box 38000, 1030BN Amsterdam, The Netherlands.
, ‡Current address: Department of Chemical Engineering, Arizona State University, Tempe, AZ 85287.
Abstract
Surface-treated hydrophobic aerogel (Nanogel) granules of sizes between 500 and 850 μm, 1.7 and 2.3 mm, and 0.5 and 2.3 mm are fluidized by a downward flow of oil-contaminated water in an inverse fluidization mode. Aerogel particles are nanostructured, extremely light and porous, have a very large surface area per unit mass, and are sufficiently robust to be fluidized. Their hydrophobic surface gives them a strong affinity for oil and other organic compounds, with the exclusion of water. These desirable properties make them an ideal sorbent or filter media for the removal of oil from wastewater. The hydrodynamic characteristics of inverse fluidized beds of aerogel granules of different size ranges were studied by measuring the pressure drop and bed expansion as a function of superficial velocity. The oil removal efficiency and capacity of the aerogel granules in the inverse fluidized bed were found to depend mainly on the size of the granules, the initial height of the bed (amount of powder used), the void fraction, and the fluid velocity. Among the advantages of the process are the extremely low energy consumption (low pressure drop) during oil removal and the large absorption capacity. Oil concentrations of about 2000 mg/L in water could be reduced to less than 10 mg/L by the inverse fluidization process.
Introduction
Experimental Methods
000 mg/L COD) from the homogenized sample is taken and inserted into the COD digestion vial which is kept in the digestion reactor (DRB-200, Hach Co.) at 150 °C for 2 h. Once the digestion is complete, the vial was allowed to cool down and then tested for COD content by the HACH colorimeter.Results and Discussion

Figure 2. Inverse fluidized bed pressure drop vs superficial fluid velocity of small aerogel granules.

Figure 3. Inverse fluidized bed pressure drop vs superficial fluid velocity of large aerogel granules.

Figure 4. Bed height vs superficial fluid velocity corresponding to the data in Figure 2.

Figure 5. Bed height vs superficial fluid velocity corresponding to data in Figure 3.
while the gravity force is
and the drag force is represented by
A force balance on the particles gives
Substituting eqs 1, 2, and 5 into eq 4 gives
After simplification and rearrangement of the terms,
and with the definition of the density of the granules
eq 7 becomes
Using eq 9, the volume of the particles is given by
It is important to note that Vp is independent of the bed height.
The void volume can be found by subtracting the volume of the particles (Vp) from the total volume of the fluidized bed (Vb)
hence, the void fraction of the fluidized bed is given by
The solid fraction is given by
and the density of the particles is given by eq 8. Finally, the internal porosity of the particles can be found by
| granule size/type(mm/type) | mass(kg) | ΔP (Pa) | particles volume(m3) | ρp(estim.; kg/m3) | initial bedheight (m) | bulk density(kg/m3) | void fraction(ε0) |
|---|---|---|---|---|---|---|---|
| 0.5−0.85 TLD 101 | 0.106 | 1185.9 | 8.7E–04 | 121 | 0.264 | 65 | 0.47 |
| 0.053 | 579.2 | 4.3E–04 | 123 | 0.143 | 60 | 0.52 | |
| 0.068 | 730.8 | 5.4E–04 | 126 | 0.183 | 60 | 0.52 | |
| 1.7–2.3 TLD 302 | 0.197 | 2240.8 | 1.6E–03 | 119 | 0.484 | 66 | 0.45 |
| 0.1 | 1103.2 | 8.2E–04 | 123 | 0.216 | 75 | 0.39 | |
| 0.13 | 1516.8 | 1.1E–03 | 117 | 0.306 | 68 | 0.41 | |
| 0.5–2.3 TLD 302 | 0.133 | 1482.4 | 1.09E–03 | 121 | 0.317 | 67 | 0.45 |
where U is the superficial velocity and Ui is the settling velocity at infinite dilution. The Richardson−Zaki exponent or index (n) is a function of the particle terminal Reynolds number (Ret) and the particle to column diameter ratio as given below
where eq 17 is specifically applicable for the smaller aerogel particles (0.5−0.85 mm) and eqs 17 and 18 are applicable for the larger aerogel granules (1.7−2.3 mm) depending on the Reynolds number, Ret. In these equations, the Reynolds number at terminal velocity is defined as
The RZ exponent (n) can also be obtained from experimental data by plotting the logarithm of the superficial velocity against the logarithm of the void fraction
where the slope of a linear regression of the data gives the Richardson−Zaki exponent (n) and the y-intercept gives the settling velocity at infinite dilution (Ui).
In our experiments, the value of dp/D is very small and is in the range of 0.0055−0.025. The terminal velocity is given by (see the work of Sakiadis(40))
where Cd is the drag coefficient which is a function of the particle Reynolds number,

Figure 6. Relationship between the superficial velocity and the void fraction of inverse fluidized beds of aerogel granules accordingly to the Richardson−Zaki equation.
Only eq 26 is applicable to our data since Rep in our experiments is greater than 0.1 and less than 1000. On the basis of our experimental data and the equations above, the Richardson−Zaki exponent (n), the terminal velocity (Ut), and the average particle diameter were calculated for the different experimental runs as shown in Table 2. As seen in the table, the values of the Richardson−Zaki exponent (n), for all three particle size ranges investigated, calculated from eq 17 agree very closely with the values obtained from the slopes of the three straight lines in Figure 6.| granule size/type(mm/type) | dp0 (m) | Rep(eq 24) | Cd(eq 26) | Ret(eq 20) | RZ (exp)(n) | RZ (eq 17) (n) | RZ(ln (Ui)) | Ut(eq 22; m/s) | dp(eq 23; mm) |
|---|---|---|---|---|---|---|---|---|---|
| 0.5−0.85 | 8.5E−04 | 12.3 | 3.53 | 36.34 | 3.20 | 3.18 | 1.56 | 0.0487 | 0.727 |
| 1.7−2.3 | 2.1E−03 | 83.2 | 1.18 | 316.36 | 2.49 | 2.47 | 2.66 | 0.1514 | 2.366 |
| 0.5−2.3 | 1.6E−03 | 59.2 | 1.39 | 207.42 | 2.83 | 2.77 | 2.51 | 0.1303 | 2.068 |
are ε =0.4−0.88, dp/D = 0.062−0.250 and 110
000 < Ar < 7
650
000. Since our results are in good agreement with eq 30, it appears that its range of applicability can be extended to particles with Archimedes numbers less than 100
000, such as aerogel granules.
Figure 7. Relationship between the drag force function “f” as defined by Fan et al.(26) and the void fraction ε.
This equation should be applicable to inverse fluidization as well assuming that the drag force of the fluid moving with superficial velocity (Umf) is equal to the buoyancy force less the weight of the particles as described by Karamanev et al.(28) In Table 3, the experimentally measured values of Umf are compared against the theoretical values of Umf calculated from eq 31. This equation correctly estimates the minimum fluidization velocity of the large aerogel granules but not the small ones; this could be due to the fact that the Archimedes number for small granules is almost 2 orders of magnitude lower than that of the large granules.| granulesize (mm) | Umf(exp; m/s) | dp(mm) | Ar | Remf | Umf(eq 31; m/s) | error(%) |
|---|---|---|---|---|---|---|
| 0.5−0.85 | 0.006 | 0.7 | 3469 | 2.0 | 0.003 | 52 |
| 1.7−2.3 | 0.018 | 2.2 | 95103 | 37 | 0.017 | 9 |
| 0.5−2.3 | 0.013 | 1.6 | 34885 | 17 | 0.011 | 19.4 |
| figure | fluidvelocity(m/s) | U/Umfratio | particletype | particlesize (mm) | mass ofparticles(kg) | upstream oil conc(g of oil/kg H2O) | COD entrance(mg/L) | initial fluidbed height(m) | max. bedheight (m) | time atCOD = 100 (s) | removalcapacity(kg oil/kg) |
|---|---|---|---|---|---|---|---|---|---|---|---|
| 8 | 0.0305 | 1.5 | TLD 302 | 1.7−2.3 | 0.054 | 0.26 | 490 | 0.16 | 0.28 | 3660 | 3.3 |
| 9 | 0.0305 | 1.5 | OGD 303 | 1.7−2.3 | 0.049 | 0.39 | 1000 | 0.14 | 0.23 | 2340 | 3.5 |
| 10 | 0.030 | 4.4 | OGD 303 | 0.5−.85 | 0.056 | 0.26 | 490 | 0.39 | 0.44 | 2200 | 1.9 |
| 12 | 0.024 | 3.5 | OGD 303 | 0.5−.85 | 0.056 | 0.36 | 900 | 0.32 | 0.33 | 5520 | 5.3 |
| 13 | 0.024 | 3.5 | TLD 101 | 0.5−.85 | 0.056 | 0.36 | 900 | 0.34 | 0.34 | 4500 | 4.3 |
| 14 | 0.011 | 1.5 | TLD 101 | 0.5−.85 | 0.056 | 0.45 | 1700 | 0.14 | 0.23 | 13620 | 7.1 |
| 15 | 0.018 | 2.6 | TLD 101 | 0.5−.85 | 0.056 | 0.48 | 1800 | 0.22 | 0.22 | 6840 | 6.5 |
| 17 | 0.010 | 1.5 | TLD 101 | 0.5−.85 | 0.108 | 0.47 | 1750 | 0.43 | 0.44 | 15360 | ? |

Figure 8. Chemical oxygen demand (COD) and inverse fluidized bed expansion (squares) as a function of time of 54 g of aerogel granules (TLD 302) with sizes between 1.7 and 2.3 mm during removal of oil from water (0.26 g of oil/kg of water and fluid velocity of 0.0305 m/s).

Figure 9. Chemical oxygen demand (COD) and inverse fluidized bed expansion (squares) as a function of time of 49 g of aerogel granules (OGD 303) with sizes between 1.7 to 2.3 mm during removal of oil from water (0.39 g of oil/kg of water and fluid velocity of 0.0305 m/s).

Figure 10. Chemical oxygen demand (COD) and inverse fluidized bed expansion (squares) as a function of time of 56 g of aerogel granules (OGD 303) with sizes between 0.5 to 0.85 mm during removal of oil from water (0.26 g of oil/kg of water and fluid velocity of 0.0305 m/s).

Figure 12. Chemical oxygen demand (COD) and inverse fluidized bed expansion (squares) as a function of time of 56 g of aerogel granules (OGD 303) with sizes between 0.5 and 0.85 mm during removal of oil from water (0.36 g of oil/kg of water and 0.0244 m/s fluid velocity).

Figure 13. Chemical oxygen demand (COD) and inverse fluidized bed expansion (squares) as a function of time of 56 g of aerogel granules (TLD 101) with sizes between 0.5 and 0.85 mm during removal of oil from water (0.36 g of oil/kg of water and 0.0244 m/s fluid velocity).
000 s to produce a level of 100 mg/L of COD downstream of the bed, as compared to 7000 s at the higher fluid velocity (0.0183 m/s). It can be seen from the data that a higher oil removal capacity is obtained at the lower flow velocity. Also, at lower fluid velocity, the drag force is lower allowing the granules to further saturate. In this case, the bed height remains almost constant for about 2 h after which the granules become saturated and groups of particles move downward expanding the bed in a short period of time (3000 s). On the other hand, when the fluid velocity is larger, the voids in the fluidized bed are larger and the drag force is also higher; hence, partially oil-saturated granules leave the bed due to entrainment. This is reflected in a reduction of the fluidized bed height over time.
Figure 15. Chemical oxygen demand (COD) and inverse fluidized bed expansion (squares) as a function of time of 56 g of aerogel granules (TLD 101) with sizes between 0.5 and 0.85 mm during removal of oil from water (0.48 g of oil/kg of water and 0.0183 m/s fluid velocity).

Figure 14. Chemical oxygen demand (COD) and inverse fluidized bed expansion (squares) as a function of time of 56 g of aerogel granules (TLD 101) with sizes between 0.5 and 0.85 mm during removal of oil from water (0.45 g of oil/kg of water and 0.0107 m/s fluid velocity).

Figure 17. Chemical oxygen demand (COD) and inverse fluidized bed expansion (squares) as a function of time of 108 g of aerogel granules (TLD 101) with sizes between 0.5 and 0.85 mm during removal of oil from water (0.47 g of oil/kg of water and 0.0102 m/s fluid velocity).
Conclusions
Acknowledgment
We gratefully acknowledge partial financial support from the National Science Foundation through Grant CBET 0730465 during the preparation and writing of this manuscript. We also wish to thank the Cabot Corporation for providing partial financial support during the experimental phase of this research, as well as the aerogel (Nanogel) granules used in the experiments. Technical discussions and advice from Cabot personnel, in particular, Dr. Nirmalya Maity and Dr. Bart Kalkstein, are also greatly appreciated.
References
This article references 42 other publications.
- 1. Johnson, R. F.; Manjrekar, T. G.; Halligan, H. R.Removal of oil from water surfaces by sorption on unstructured fibers Environ. Sci. Technol. 1973 7 439 443[ACS Full Text
], [ChemPort] - 2. Paterson, J. W. Industrial Wastewater Treatment Technology, 2nd ed.; Butterworth Publishers, Inc.: Stoneham, MA, 1985.
- 3.
- 4. Manual on Disposal of Refining Wastes; American Petroleum Institute: Washington, D.C., 1969; Vol. 5, pp 5− 15
- 5. Manning, F.; Snider, E. H.Assesment data base for petroleum refining wastewater and residues; U.S. Department of Commerce, NTIS: Washington, D.C., 1983; Vol. 9, pp 4− 101
- 6. Tchobanoglous, G.; Burton, L. F.; Stensel, H. D. Wastewater Engineering: Treatment and Reuse; Mc-Graw Hill: New York, 2003; pp 419− 422.
- 7. Bennett, G. F.The removal of oil from wastewater by air flotation: A review. CRC Crit. Rev. Environ. Contr.1988, 18 ( 3), 189− 253
- 8.
- 9.
- 10.
- 11. Goldsmith, R.; Hossian, S.Ultrafiltration concept for separating oil from water; U.S. Coast Guard: Washington, D.C., January 1973.
- 12. Cambiella, A.; Ortea, E.; Rios, G.; Benito, J. M.; Pazos, C.; Coca, J.Treatment of oil-in-water emulsions: Performance of a sawdust bed filter J. Hazard. Mater. 2006 B131 195 199[CrossRef]
- 13. Mathavan, G. N.; Viraraghavan, T.Use of peat in the treatment of oily waters Water, Air, Soil Pollut. 1989 45 17 26[ChemPort]
- 14. Viraraghavan, T.; Moazed, H.Removal of oil from water by bentonite Fresenius Environ. Bull. 2003 12 9 1092 1097[ChemPort]
- 15.
- 16. Mysore, D.; Viraraghavan, T.; Jin, Y. C.Vermiculite Filtration for Removal of Oil from Water Pract. Period. Hazard., Toxic, Radioact. Waste Manage. 2006 156 161[CrossRef]
- 17.
- 18.
- 19. Sorbent Materials in Storm Water Applications.http://epa.gov/owm/mtb/mtbfact.htm (accessed July 2006).
USEPA.
- 20.
- 21. Hrubesh, L. W.; Coronado, P. R.; Dow, J. P.Method for removing organic liquids from aqueous solutions and mixtures. U.S. Patent 6709600 B2. March 23, 2004.
- 22. Hrubish, L. W.; Coronado, P. R.; Satcher, J. H.Solvent removal from water with hydropholic aerogels J. Non-Cryst. Solids 2001 285 1−3 328 332[CrossRef]
- 23. Ayres, M.; Hunt, A.http://eetd.lbl.gov/ECS/aerogels/satoc.htm (updated April 2004), Ernest Orlando Lawrence Berkeley National Laboratory. Environmental Energy Technologies Division (EETD).
Silica aerogels.
- 24. https://ipo.llnl.gov/technology/ profile/aerogel/GACComposites/index.php (accessed July 2006), Lawrence Livermore National Laboratory.
Aerogel/Granulated Activated Carbon composites for adsorption of contaminants in water.
- 25. Richardson, J. F.; Zaki, W. N.Sedimentation and fluidization: Part I Trans. Instn. Chem. Eng. 1954 32 35 53[ChemPort]
- 26.
- 27.
- 28.
- 29.
- 30. Ibrahm, Y. A. A.; Briens, C. L.; Margaritis, A.; Bergongnou, M. A.Hydrodynamic Characteristics of a Three-Phase Inverse Fluidized-Bed Column AIChE J. 1996 42 7 1889 1900[CrossRef]
- 31.
- 32.
- 33. Cho, Y. J.; Park, H. Y.; Kim, S. W.; Kang, Y.; Kim, S. D.Heat Transfer and Hydrodynamics in Two-and Three-Phase Inverse Fluidized Beds Ind. Eng. Chem. Res. 2002 41 2058 2063[ACS Full Text
], [ChemPort] - 34. Renganathan, T.; Krishnaiah, K.Prediction of Minimum Fluidization Velocity in Two and Three Phase Inverse Fluidized Beds Can. J. Chem. Eng. 2003 81 853 860[ChemPort]
- 35. Renganathan, T.; Krishnaiah, K.Stochastic Simulation of Hydrodynamics of a Liquid-Solid Inverse Fluidized Bed Ind. Eng. Chem. Res. 2004 43 4405 4412[ACS Full Text
], [ChemPort] - 36.
- 37.
- 38. Jirka, A. M.; Carter, M. J.Reactor digestion method for COD analysis Anal. Chem. 1975 47 8 1397
- 39. Method 8000: Reactor digestion method USEPA approved for oxygen chemical demand wastewater analysis. DR/890 Datalogging Colorimeter Handbook; Hach Co.: Loveland, CO, 2004; pp 427− 436
- 40. Sakiadis, B. C.Fluid and particle mechanics. Perry’s Chemical Engineer’s Handbook; McGraw-Hill: New York, 1984
- 41.
- 42.
Tools
-
Add to Favorites
-
Download Citation
-
Email a Colleague -
Permalink
Order Reprints
Rights & Permissions
Citation Alerts
History
- Published In Issue January 07, 2009
- Article ASAPAugust 16, 2008
- Received: January 4, 2008
Accepted: May 22, 2008
Revised: May 19, 2008
Cart





