Water treatment and distribution

This chapter covers the following:

• Population modeling
• Mass balance and calculations
• Steady-State reactor parameters (constant density systems)
• Specific gravity for a solids slurry
• NRC trickling filter performance
• Activated sludge

Population modeling

Population modeling helps you estimate current and future population based on trends such as births, deaths, and migration. In environmental and water resources engineering, population projections matter because design flows and loading rates for water supply and wastewater systems depend directly on how many people the system will serve.

Different projection methods fit different situations:

• Linear (arithmetic) growth: adds a constant number of people per unit time
• Exponential (geometric) growth: grows by a constant fraction per unit time
• Logistic (decreasing-rate-of-increase) growth: growth slows as the population approaches a saturation value
• Ratio/correlation methods: tie local growth to growth in a larger reference region

Linear projection = algebraic projection

Use this model when the population increases by about the same number of people each year.

$$P_t=P_0+k\Delta t$$

where

• $P_t$ = population at time $t$
• $P_0$ = population at time zero
• $k$ = growth rate
• $\Delta t$ = elapsed time in years relative to time zero

Log growth = exponential growth = geometric growth

Use this model when the population grows proportionally to its current size (a constant percentage-type growth).

$$P_t=P_0{e}^{k\Delta t}$$

$$\ln P_t=\ln P_0+k\Delta t$$

where

• $P_t$ = population at time $t$
• $P_0$ = population at time zero
• $k$ = growth rate
• $\Delta t$ = elapsed time in years relative to time zero

Percent growth

This form is often used when growth is applied in discrete periods (for example, year-by-year compounding).

$$P_t=P_0(1+k{)}^n$$

where

• $P_t$ = population at time $t$
• $P_0$ = population at time zero
• $k$ = growth rate
• $n$ = number of periods

Ratio and correlation growth

This approach assumes the community’s population stays in a constant ratio to a larger region’s population.

$$\frac{P_2}{P_{2R}}=\frac{P_1}{P_{1R}}=k$$

where

• $P_2$ = projected population
• $P_{2R}$ = projected population of a larger region
• $P_1$ = population at last census
• $P_{1R}$ = population of larger region at last census
• $k$ = growth ratio constant

Decreasing-rate-of-increase growth

Use this model when growth slows over time as the population approaches a maximum (saturation) value.

$$P_t=P_0+(S-P_0)(1-e^{-k(t-t_0)})$$

where

• $P_t$ = population at time $t$
• $P_0$ = population at time zero
• $k$ = growth rate constant
• $S$ = saturation population
• $t,t_0$ = future time, initial time

Mass balance and calculations

A mass balance applies conservation of mass to a defined control volume. The basic idea is:

• mass can enter
• mass can leave
• mass can be created or destroyed by reactions (or other sources/sinks)
• whatever is left shows up as accumulation inside the control volume

Mass balances are used throughout environmental engineering for reactor analysis, mixing problems, and pollutant fate and transport under both steady-state and dynamic conditions.

Mass balance:

$$\frac{dM}{dt}=\frac{d{M}_{in}}{dt}-\frac{d{M}_{out}}{dt}\pm r$$

$$M=CQ=CV$$

Continuity equation:

$$Q=vA$$

where

• $M$ = mass
• $M_{in}$ = mass in
• $M_{out}$ = mass out
• $r=k{C}^n$ = reaction rate

$$k=\text{reaction rate constant}\left(\frac{1}{{\text{(concentration units)}}^{n-1}\cdot \text{time}}\right)$$

• $n$ = order of reaction
• $C$ = concentration (mass/volume)
• $Q$ = flow rate
• $V$ = volume
• $v$ = velocity
• $A$ = cross-sectional area of flow

Unit conversion formula:

This conversion is commonly used to compute mass loading from concentration and flow.

$$M(\text{lb/day})=C(\text{mg/L})\times Q(\text{MGD})\times 8.34[\text{lb-L/(mg-MG)}]$$

where

• $\text{MGD}$ = million gallons per day
• $\text{MGD}$ = million gallons

Steady-State reactor parameters (constant density systems)

In steady-state reactor analysis, key properties (such as flow rate and concentration at a given location) don’t change with time. For constant density systems, this assumption simplifies the governing equations and lets you relate:

• detention (retention) time
• reaction kinetics (reaction order and rate constant)
• effluent concentration

These relationships are used for ideal reactor models such as the completely mixed flow reactor (CMFR) and the plug flow reactor (PFR).

Steady-state retention times ($t$) for decay reactions of different order

Equations for mean retention times ($t$)

Reaction order	$r$	Ideal batch	Ideal plug flow	Ideal CMFR
Zero${}^b$	$-k$	$\frac{C_0-C_t}{k}$	$\frac{C_0-C_t}{k}$	$\frac{C_0-C_t}{k}$
First	$-kC$	$\frac{\ln (C_0/C_t)}{k}$	$\frac{\ln (C_0/C_t)}{k}$	$\frac{(C_0/C_t)-1}{k}$
Second	$-k{C}^2$	$\frac{(C_0/C_t)-1}{k{C}_0}$	$\frac{(C_0/C_t)-1}{k{C}_0}$	$\frac{(C_0/C_t)-1}{k{C}_0}$

${}^a$ Expressions are valid for $kt\le C_0$; otherwise $C_t=0$ ${}^b{C}_0$ = initial concentration or influent concentration; $C_t$ = final condition or effluent concentration.

Steady-state performance for decay reactions of different order

Equations for $C_t$

Reaction order	$r$	Ideal batch	Ideal plug flow	Ideal CMFR
Zero${}^b$	$-k$	$C_0-kt$	$C_0-kt$	$C_0-kt$
First	$-kC$	$C_0[\exp (-kt)]$	$C_0[\exp (-kt)]$	$\frac{C_0}{1+kt}$
Second	$-k{C}^2$	$\frac{C_0}{1+kt{C}_0}$	$\frac{C_0}{1+kt{C}_0}$	$\frac{(4kt{C}_0+1{)}^{1/2}-1}{2kt}$

${}^b{C}_0$ = initial concentration or influent concentration; $C_t$ = final condition or effluent concentration.

Note: Here $\frac{C_t}{C_0}=1-X$, where, $x$ is known as fractional conversion that defined as the moles of material reacted per mole of that material fed.

Specific gravity for a solids slurry

For a solids slurry, specific gravity captures the combined effect of:

• how much solid is present
• the density of the solid particles
• the liquid phase properties

You’ll use slurry specific gravity in hydraulic calculations, pump selection, and process design for systems that handle suspended or settled solids.

$$S=\frac{W_w+W_s}{(W_w/1.0)+(W_s/S_s)}$$

where

• $S$ = specific gravity of wet sludge
• $W_w$ = weight of water (lb or kg)
• $W_s$ = weight of dry solids (lb)
• $S_s$ = specific gravity of dry solids

Volume of waste sludge

$$V=\frac{W_s}{(s/100)\cdot \rho S}=\frac{W_s}{\left(\frac{100-p}{100}\right)\cdot \rho S}$$

where

• $V$ = volume of sludge, ft³ (gal or m³)
• $W_s$ = weight of dry solids (lb or kg)
• $s$ = solids content (%)
• $\rho$ = unit weight of water (62.4 lb/ft³ or 1,000 kg/m³)
• $S$ = specific gravity of wet sludge
• $p$ = water content (%)

NRC trickling filter performance

The NRC (National Research Council) trickling filter performance equations are empirical relationships used to estimate organic removal efficiency under different loading conditions. They incorporate key factors such as:

• influent substrate concentration
• hydraulic loading
• filter depth
• recirculation

Because they’re straightforward to apply, NRC equations are commonly used for preliminary design and performance checks of conventional trickling filter systems.

Single-stage or first-stage rock filter:

$$E_t=\frac{100}{1+0.0541\cdot \frac{W}{VF}}$$

where

• $E_t$ = efficiency of BOD removal at 20°C
• $W$ = BOD loading (lb/day)
• $V$ = volume of filter media (10³ ft³)
• $F$ = recirculation factor

Recirculation factor:

$$F=\frac{1+R}{(1+R/10{)}^{0.5}}$$

Activated sludge

The activated sludge process is a biological treatment method that uses suspended microorganisms to remove organic matter (and, in some configurations, nutrients). Aeration and mixing keep the biomass in suspension and provide oxygen so microorganisms can metabolize biodegradable pollutants. The biomass forms flocs that are later separated from the treated effluent by sedimentation.

Common design and operating concepts include mixed liquor suspended solids, sludge age, food-to-microorganism ratio, and oxygen demand.

Sludge production equation:

$$X_e=\frac{{\theta }_{c}Y(S_0-S_e)}{\theta (1+k_d{\theta }_c)}$$

Steady-State mass balance around secondary clarifier:

$$(Q_0+Q_R)X_r-Q_0{X}_e-Q_w{X}_w+Q_e{X}_e=\frac{V}{{\theta }_c}X$$

Solids residence time:

$${\theta }_c=\frac{VX}{Q_0{X}_e+Q_w{X}_w}$$

Sludge volume/day:

$$Q_s=\frac{M}{{\rho }_s(\%\text{solids})}$$

Sludge volume index (SVI):

$$SVI=\frac{\text{Sludge volume after settling (mL/L)}\times 1,000}{\text{MLSS (mg/L)}}$$

where

• $k_d$ = microbial death rate (day⁻¹)
• $S_e$ = effluent BOD or COD (kg/m³)
• $S_0$ = influent BOD or COD (kg/m³)
• $X$ = biomass concentration (MLSS or MLVSS, kg/m³)
• $X_e$ = effluent suspended solids (kg/m³)
• $X_r$ = recycle suspended solids (kg/m³)
• $X_w$ = waste sludge suspended solids (kg/m³)
• $\theta$ = hydraulic residence time
• ${\theta }_c$ = solids retention time
• $Y$ = yield coefficient (kg biomass/kg BOD)
• $M$ = sludge production rate (dry weight basis)
• ${\rho }_s$ = density of solids
• $Q_0$ = influent flow rate
• $Q_e$ = effluent flow rate
• $Q_R$ = recycle flow rate
• $Q_w$ = waste flow rate
• $R$ = recycle ratio = $Q_R/Q_0$
• $V$ = aeration basin volume

Clarifier design equations

Overflow rate = Hydraulic loading rate:

$$v_0=\frac{Q}{A_{\text{surface}}}$$

Terminal settling velocity ($v_c$):

Velocity at which the smallest particle is 100% removed.

Weir loading rate:

$$\text{WOR}=\frac{Q}{\text{Weir Length}}$$

Horizontal velocity:

$$v_h=\frac{Q}{A_{\text{cross-section}}}=\frac{Q}{A}$$

Hydraulic residence time:

$$\theta =\frac{V}{Q}$$

where

• $Q$ = flow rate
• $A_s$ = surface area
• $V$ = tank volume

Typical primary clarifier efficiency percent removal

Overflow rates (gpd/ft²)	1,200	1,000	800	600
Overflow rates (m/d)	48.9	40.7	32.6	24.4
Suspended solids	30%	32%	34%	36%
BOD₅	30%	32%	34%	36%

Note: gpd = gallons per day