Victor M. Ponce Professor of Civil Engineering San Diego State University [110926] ABSTRACT The differential equation for the dissolved oxygen sag curve (DO sag curve) is derived. The solution of this differential equation can be shown to be essentially the same as that of the well known StreeterPhelps equation (Streeter and Phelps, 1925). Unlike the latter, the differential equation derived herein can be solved numerically and, therefore, does not require integration. Moreover, the differential equation is valid for all deoxygenation and oxygenation constants, unlike the StreeterPhelps equation, which is undefined when these constants are equal. Two online calculators: (a) single case, and (b) general case, round up the analysis. 1. INTRODUCTION The discharge of a sewage effluent in a stream produces a biochemical oxygen demand (BOD) which decays exponentially in time and space. This oxygen demand causes an oxygen deficit, or oxygen shortage. The greater the oxygen deficit, the greater the rate of natural oxygen replenishment from the atmosphere into the stream. These two concurrent processes of oxygen consumption and oxygen replenishment produce an oxygen sag curve, i.e., a curve which sags initially due to the increased oxygen demand and recovers asymptotically downstream due to the increased rate of oxygen replenishment. The derivation of the DO sag equation is the objective of this article. 2. OXYGEN DEMAND The exponential decay of oxygen demand can be modeled as follows (Tchobanoglous and Schroeder, 1985):
in which: D = biochemical oxygen demand (BOD) (mg/liter); D_{u} = ultimate BOD immediately downstream of effluent discharge (mg/liter); k_{d} = deoxygenation constant (d^{1}); t = time (d). The differential equation for oxygen demand is:
Using the chain rule, the oxygen demand equation is converted to the spatial domain:
in which: x = distance along the stream, measured downstream of effluent discharge (m); v = stream velocity (m/d). 3. ULTIMATE BOD At the upstream boundary, a mass balance leads to:
in which:
D_{u} = ultimate BOD immediately downstream of effluent discharge (mg/liter); Q_{s} = stream discharge (m^{3}/s); D_{s} = BOD in stream, immediately upstream of effluent discharge (mg/liter);
Q_{e} = effluent discharge (m^{3}/s);
D_{e} = BOD of effluent discharge (mg/liter); Assuming the upstream flow is clean: D_{s} ≅ 0. Therefore, the mass balance reduces to:
4. OXYGEN SUPPLY
The differential equation for oxygen supply in a stream can be modeled as follows (Tchobanoglous and Schroeder, 1985):
in which: k_{o} = oxygenation constant (d^{1}); S_{s} = dissolved oxygen concentration at saturation, a function of temperature, salinity, and atmospheric pressure (mg/liter); and S = dissolved oxygen concentration (mg/liter). The quantity (S_{s}  S ) is the oxygen deficit. Using the chain rule, the oxygen supply equation is converted to the spatial domain:
Therefore:
5. DO SAG EQUATION The change of dissolved oxygen concentration in space is given by the change in supply plus the change in demand:
in which O = dissolved oxygen concentration (mg/liter).
Combining Eqs. 4, 9, and 10:
Since S = O:
In a control volume (computational reach) of length L (m), with index j upstream
and
Therefore:
6. INITIAL DO
For j = 0: O_{j} = O_{0}. At the upstream boundary, a mass balance leads to:
in which:
O_{0} = dissolved oxygen concentration immediately downstream of effluent discharge (mg/liter); Q_{s} = stream discharge (m^{3}/s); O_{s} = dissolved oxygen concentration in stream, immediately upstream of effluent discharge (mg/liter);
Q_{e} = effluent discharge (m^{3}/s);
O_{e} = dissolved oxygen concentration of effluent discharge (mg/liter); Assuming that the effluent is anoxic: O_{e} ≅ 0:
7. STREETERPHELPS EQUATION The classical way of solving for the dissolved oxygen sag equation is the StreeterPhelps equation, which dates back to 1925 (StreeterPhelps, 1925; Tchobanoglous and Schroeder, 1984). The StreeterPhelps equation is:
in which:
O = dissolved oxygen concentration (mg/liter) at distance x (m); S_{s} = dissolved oxygen concentration at saturation, a function of temperature, salinity, and atmospheric pressure (mg/liter); and O_{0} = dissolved oxygen concentration immediately upstream of effluent discharge (mg/liter). The StreeterPhelps equation is an algebraic equation derived by integrating the differential equation governing the oxygen sag. The differential equation derived herein is based on the same principles as that of StreeterPhelps. However, the differential equation is not integrated as in the case of StreeterPhelps, but rather it is solved directly, using finite differences. In most cases, both methods will lead to the same answer. Significantly, the StreeterPhelps equation suffers from being undefined when the oxygenation and deoxygenation constants are equal, while the numerical method is not. Thus, the numerical method is an allaround better predictor than the StreeterPhelps model, applicable for all oxygenation and deoxygenation constants, regardless of their values. 8. ONLINE CALCULATION The procedure to calculate the DO sag curve is implemented in the following online calculators:
9. SUMMARY The dissolved oxygen concentration at the upstream boundary, at j = 0, immediately downstream of the effluent discharge, is estimated as follows:
The ultimate BOD immediately downstream of the effluent discharge, is estimated as follows:
The dissolved oxygen concentration at node ( j + 1) is:
Given: Q_{s} (m^{3}/s), Q_{e} (m^{3}/s), S_{s} (mg/liter), O_{s} (mg/liter), D_{e} (mg/liter), v (m/d), L (m), k_{d} (d^{1}), and k_{o} (d^{1}), the dissolved oxygen concentration along the stream (DO sag curve) may be calculated using Eq. 14. Two online calculators: (a) single case, and (b) general case, round up the analysis. REFERENCES Streeter, H. W., and E. B. Phelps. 1925. A study of the pollution and natural purification of the Ohio River. III. Factors concerned in the phenomena of oxidation and reaeration. U.S. Public Health Service, Bulletin No. 146. Tchobanoglous, G., and E. D. Schroeder. 1984. Water quality: Characteristics, modeling, modification. AddisonWesley, Massachussets.
NOTATION D = biochemical oxygen demand (BOD) (mg/liter); D_{e} = BOD of effluent discharge (mg/liter); D_{s} = BOD in stream, immediately upstream of effluent discharge (mg/liter); D_{u} = ultimate BOD immediately downstream of effluent discharge (mg/liter);
k_{d} = deoxygenation constant (d^{1}); k_{o} = oxygenation constant (d^{1}); L = length of control volume, or computational reach (m); O = dissolved oxygen concentration (mg/liter); O_{0} = dissolved oxygen concentration immediately downstream of effluent discharge (mg/liter); O_{s} = dissolved oxygen concentration in stream, immediately upstream of effluent discharge (mg/liter); O_{e} = dissolved oxygen concentration of effluent discharge (mg/liter); Q_{e} = effluent discharge (m^{3}/s); Q_{s} = stream discharge (m^{3}/s); S = dissolved oxygen concentration (mg/liter); S_{s} = dissolved oxygen concentration at saturation, a function of temperature, salinity, and atmospheric pressure (mg/liter); t = time (d); x = distance along the stream, measured downstream of effluent discharge (m); and v = stream velocity (m/d).

201117 13:21 