Supplementary Material (ESI) for Green Chemistry

This journal is The Royal Society of Chemistry 2004

Appendix. Granular Activated Carbon Adsorption. Process Design and Modelling.

Determination of the main equipment size

Adsorber Column. The volume of the adsorber column can be estimated through the following

expression:

Vads At . H (A.1)

According to Faust and Aly,(1987), the height of the adsorber column, H, should be higher than 1.35

m in order to accommodate a packed bed long enough to ensure formation of a constant pattern. In

addition, a maximum height of 18.5 m should be attained to ensure shipment safety limitations.

The volume of the adsorber, should account for bed expansion and maintenance, thus, the activated

carbon occupies 66% of the total adsorber column, giving rise to the following expression:

Vads 15
. .Vbed (A.2)

Hand et al. (1984) developed an expression that is able to calculate the capacity of the bed, before

breakthrough or saturation occurs. This accounts for the unsteady state nature of the adsorber column,

pollutant characteristics, flow, solid-liquid equilibrium relationships and the type of GAC selected for

the process. The homogeneous surface diffusion model (HSDM) presented by Hand et al.,(1984), has

successfully predicted fixed-bed adsorber dynamics for over 100 adsorbate-adsorbent systems,

including potentially harmful organics such as benzene and toluene.

The model relies on a set of dynamic equations that model an adsorption system and have been solved

analytically to generate algebraic equations. In order to achieve an exact solution of the equations, the

following assumptions must be incorporated into the model: (1) plug flow exists in the bed, which

was confirmed to be valid as long as the mass transfer zone is longer than 30 adsorbent particle

diameters, and is supported by Miyachi and Kikuchi, (1975) and Stenzel (1993); (2) the hydraulic

loading (inlet flowrate) is constant; (3) surface diffusion is the predominant intraparticle mass transfer

mechanism and is not a function of concentration; (4) the local liquid-phase mass transfer rate can be

described by the linear driving force approximation; (5) the adsorbent is in a fixed position in the
Supplementary Material (ESI) for Green Chemistry
This journal is The Royal Society of Chemistry 2004

adsorber and is considered to be spherical; and (6) the adsorption equilibrium can be described by the

Freundlich isotherm. The analytical solution to the HSDM for the case in which liquid-phase mass

transfer controls the adsorption rate and 1/n is less than 1.0 has been developed by Fleck, (1973) and

presented by Hand et al. (1984).

1
Dg Cout
Cout n 1
T 1 ln C n 1 .ln 1 C 1 (A.3)
3. St. D g 1
in in

in which, n is the Freundlich isotherm intensity constant, and

1 1

n n



1 k 1 1 1 (A.4)
1 k k1
n n n

Mass throughput, T, is defined as the ratio of mass fed compared to the mass required to saturate the

adsorber. The dimensionless time scale, mass throughput, is defined (Hand et al., 1984) as:

Time
T

. Dg 1 (A.5)

The solute distribution parameter, Dg, is defined as the ratio of the mass of adsorbate in the solid-

phase to the mass of adsorbate in the liquid-phase under equilibrium conditions. This parameter is

calculated from the following expression:

p .q* 1
Dg (A.6)
. Cin

The particle density, p, and the void fraction, are parameters that depend on the type of GAC

employed. These calculations assume the use of Calgon Filtrasorb 300 as carbon adsorbent.

Liquid-solid equilibrium between the concentration of pollutant adsorbed on the carbon surface and

the concentration of pollutant in the untreated wastewater, q*, is described by the Freundlich isotherm.

q * K . Cin1/ n (A.7)

Both the capacity constant, K, and the intensity constant, 1/n, are solute specific and are available in

several publications. The values used for this model are those presented by Faust and Aly, (1987).
Supplementary Material (ESI) for Green Chemistry
This journal is The Royal Society of Chemistry 2004

The void fraction of the packed bed, , depends on both the apparent density, a, and the particle

density, p, of the activated carbon. Thus, it can be described by the following expression:

a
1 (A.8)
p

The Stanton number, St, is a measure of the bed length as compared to the length of the mass transfer

zone for the case in which the liquid-phase mass transfer resistance controls the adsorption rate.

k f . . 1
St (A.9)
R. .

This expression depends on the characteristics of the GAC employed. Two new parameters are

introduced, the mean particle ratio, R, and the sphericity, , that accounts for the ratio of surface area

of equivalent-volume sphere to actual surface area of adsorbent particle (dimensionless). The later

value can be taken from Perry (1984), whereas the particle ratio can be read from Faust and Aly,

(1987).

The packed bed contact time, , also known as fluid residence time in the packed bed, is defined as

the quotient of the free bed volume available for liquid flow divided by the wastewater flowrate

through the bed.

.Vbed
(A.10)
Q

The mass transfer coefficient for the liquid film, kf, is based on a correlation presented by Crittenden

et al., (1978).

2.40. vs
kf (A.11)
Re0.66 . Sc 0.58

This correlation is valid for Reynolds numbers, Re, ranging between 0.08 and 150, and Schmidt

numbers, Sc, between 150 and 1300, which are broad enough to cover the liquid loading applied for

this model.
Supplementary Material (ESI) for Green Chemistry
This journal is The Royal Society of Chemistry 2004

The Reynolds number is calculated from information related to the density and viscosity of the

wastewater stream, which is assumed to be equal to those of pure water due to the low concentration

of organics present.

R. v s .
Re (A.12)
.

The Schmidt number is also related to the fluid viscosity, density and diffusivity of pollutants in

water, as expressed below.

Sc (A.13)
. Dwat

The superficial liquid velocity, vs, is calculated based on the cross sectional area of the adsorption

column.

Q
vs (A.14)
At

The adsorption column is subject to a pressure drop, which depends on the superficial liquid velocity,

the height of the column and some properties of the GAC. A empirical correlation developed by

Leva, (1949) is used for its estimation. It applies for laminar flow, with Reynolds numbers below 20

and has been supported by extensive data ranging from runs performed with GAC, sands and

particles, including Berl saddles and complex ring packing.

. 11 . v s . H . . 1
2
2.570610
Pads (A.15)
2 . R 2 . 3

Eductor and Pumps. As previously mentioned, this process uses eductors to transport GAC. The

system employs water as motive pumping fluid in order to create a slurry which can be easily

transport. Therefore, this model considers the use of one pump for the inlet flow, that provides liquid

to the eductors, and a hopper, which feeds the new, and the exhausted GAC. The empirical

correlation presented by Perry, (1984) is used to calculate power requirements for the pumps.

Equations for the wastewater inlet flow and the two pump-eductor systems are presented below.

3. Q .Pads 0.0097. H
WPuWW (A.16)
36.
Supplementary Material (ESI) for Green Chemistry
This journal is The Royal Society of Chemistry 2004

FWPu F Ed . PD F Ed
WPu F Ed (A.17)
36000.

FWPu S Ed . PD S Ed
WPu S Ed
36000. (A.18)

The discharge pressure, process water requirements and suction capacity of GAC for the feed carbon

pump-eductor system were obtained from regressions of operating data presented by Karassik, (1976),

where these quantities are related to the discharge pressure of the eductor.

The discharge pressure for the spent eductor pump, process water requirements and suction capacity

for the spent carbon eductor are read from tables provided by Karassik, (1976), which in turn, follow

recommendations given by Culp and Clark, (1983) for carbon handling.

The number of carbon changes required per year determines the operating times of the individual

eductors. This time is defined by the yearly operational time, the time for pollutant breakthrough, and

the time required to exchange exhausted for fresh GAC, which is defined as the quotient between the

packed bed volume and the suction capacity of the eductor.

[Link]
OpTimePu F Ed (A.19)
Time. SC F Ed

[Link]
OpTimePu S Ed (A.20)
Time. SCS Ed

Carbon Regeneration/Exchange. Carbon Exchange is often the most critical parameter of the

adsorption system operation and may be the predominant component of the operating costs, as well as

the key to safe disposal of the organic contaminants removed from the process. This model assumes

that carbon exchange takes place once breakthrough occurs. This assumption is necessary in order to

estimate the carbon requirements and to avoid exchange zone amplifications between both adsorbers

which would require experimental data and therefore limit the flexibility of the model. This

assumption is supported by the type or organic compounds to be removed, which produce sharp

breakthrough curves, implying almost complete exhaustion of the bed when breakthrough occurs.
Supplementary Material (ESI) for Green Chemistry
This journal is The Royal Society of Chemistry 2004

Spent carbon can be landfilled or incinerated as a solid waste. Coal-based granular activated carbon,

however, is usually reactivated. The reactivation process drives off the more volatile organic

compounds from the carbon surface and pyrolyses and decomposes less volatile compounds held

within the pores (Cheremisinoff and Ellerbush, 1978). Reactivation has proven to be a safe and cost-

effective means to destroy adsorbed organic compounds and allow reuse of the carbon.

Cost equations

The total treatment cost is a function of the capital and operating cost.

Ct Co Cc (A.21)

The equipment contributing to the capital cost is seen on the flowsheet; i. e., 2 adsorber columns, the

wastewater pump, two pump-eductor systems for carbon transportation, and two storage tanks for

fresh and exhausted GAC.

C c bm dep ( 2.C Ads C Pu WW C Pu F Ed C Pu S Ed [Link] ) (A.22)

The operating costs are a function of the electricity consumed to run the process pumps, water

required to transport the GAC and replacement of GAC employed for the process. Note that

regeneration costs are not considered in the economic model since it is performed off-site.

Co Cel Cw CGAC (A.23)

The cost of the adsorber column is calculated from recommendations given by Guthrie, (1969) and

Douglas, (1988). Stainless steel is considered as the material for construction in order to resist the

reaction that may be caused by the presence of GAC, organics and water (Stenzel, M.H., 1993).

MS 95
. 6.5616. R Ads 3.2808. H 0.82 . 3.67
1.066
C Ads 1019
. (A.24)
MS 68

Similar equations, based on recommendations given by Walas and Spangler, (1960), are used to

estimate the cost of the pumps in the system.

MS 95
CPu 655.82 1.341. WPu 0.4918
MS 59
(A.25)
Supplementary Material (ESI) for Green Chemistry
This journal is The Royal Society of Chemistry 2004

The cost of the storage tanks for fresh and exhausted GAC is calculated, from recommendations given

by Stenzel, (1993). GAC supply is carried out four times a year. Therefore, it is assumed that the

volume of each storage tank is 25% of the total GAC required per year. The following equation

results from a linear regression of the plot presented in Institute Francais du Petrole (1976), for

storage tanks at atmospheric pressure with stainless steel as construction material.

MS 95
CTank 2
MS 68
2570. LN 0.25.V Bed . OpTime / Time 3100 (A.26)

Process pumps are the primary power consumers. In this process, only the wastewater feed pump

operates continuously, while the eductor pumps operate intermittently for charging and discharging

GAC in the adsorber columns. The specific costs for electricity, water and GAC were obtained from

Peters and Timmerhaus, (1980) and updated according to the Marshall and Swift Index.

Cel SC el WPu WW . OpTime WPu F Ed . Op Time Pu F Ed WPu S Ed . Op Time Pu S Ed

(A.27)

The cost of water used for the eductors is calculated as follows:

Cw SCw FWPu F Ed . OpTime Pu F Ed FWPu S Ed . OpTime Pu S Ed (A.28)

The cost of the GAC is calculated as a function of the packed bed volume and the number of loads

required per year of operation.

OpTime
CGAC SCGAC .V (A.29)
Time Bed

List of Symbols: Carbon Adsorption

At Sectional Area of Adsorber Column m2

CAds Cost of Adsorber Column /yr
Cc Capital Cost of the Process /yr
Cel Cost of electricity /yr
CGAC Cost of Granular Activated Carbon /yr
Cin Inlet Pollutant Concentration kg/m3
Co Operating Cost of the Process /yr
Cout Outlet Pollutant Concentration kg/m3
CPu-F-Ed Cost of Pump-Eductor Feed GAC System /yr
CPu-S-Ed Cost of Pump-Eductor Spent GAC System /yr
CPu-WW Cost of Wastewater Inlet Pump /yr
Ct Total Cost of the Process /yr
Cw Cost of water for the GAC pump-eductor systems /yr
dep Capital charge factor yr -1
Dg Solute Distribution Parameter --
Supplementary Material (ESI) for Green Chemistry
This journal is The Royal Society of Chemistry 2004

Dwat Diffusivity of Organic in Water m2/hr

FWPu-F-Ed Flow of water for the Pump-Eductor Feed GAC System kg/hr
FWPu-S-Ed Flow of water for the Pump-Eductor Spent GAC System kg/hr
H Height of the Adsorber Column m
K Freundlich Isotherm Capacity Constant mg/g
kf Liquid Film Mass Transfer Coefficient m/hr
MS59 Marshall and Swift Index 1959 --
MS68 Marshall and Swift Index 1968 --
MS95 Marshall and Swift Index 1995 --
N Freundlich Isotherm Intensity Constant --
OpTime Operating Time for the Process hr
Op-TimePu-F-Ed Operating Time for the Pump-Eductor Feed GAC System hr
Op-TimePu-S-Ed Operating Time for the Pump-Eductor Spent GAC System hr
PD-F-Ed Discharge Pressure for the Feed GAC Eductor bar
PPu-S-Ed Discharge Pressure for the Spent GAC Eductor bar
Q Wastewater Flowrate m3/hr
q* Equilibrium concentration between solid phase and inlet liquid kg/kg
R GAC Mean Particle Ratio m
Rads Ratio of the Adsorber Column m2
Re Reynolds Number --
Sc Schmidt Number --
SC Suction Capacity m3/hr
SCel Specific Cost of Electricity /[Link]
SCGAC Specific Cost of GAC /m3
SCw Specific Cost of Water /m3
St Stanton Number
T Mass Throughput --
TimeGAC Elapsed time for the GAC hr
VAds Volume of the Adsorber Column m3
VBed Volume of the Bed m3
vs Superficial Liquid Velocity m/hr
WPu-F-Ed Power of the Pump-Eductor Feed GAC System kW
WPu-S-Ed Power of the Pump-Eductor Spent GAC System kW
WPu-WW Power of the Wastewater Inlet Pump kW

PAds Pressure Drop Through Adsorber Column bar

Void Fraction of Packed Bed --
Sphericity of the GAC particle --
Series factor for the Throughput Equation --
Pump Efficiency --
Viscosity of Water kg/[Link]
Density of Water kg/m3
a Apparent Density of the GAC Particle kg/m3
pa Density of the GAC particle kg/m3
Residence Time for the Fluid in Packed Bed hr