 Details
 Parent Category: Engineering Assignments' Solutions
We Helped With This Engineering Assignment: Have A Similar One?
Category  Engineering 

Subject  Other 
Difficulty  Undergraduate 
Status  Solved 
More Info  Engineering Physics Homework Help 
Short Assignment Requirements
Assignment Description
Modelling water quality
Dissolved oxygen (DO) is one of the most important indicators for the quality of natural water systems as fish and other aquatic animal species require oxygen. Streams must have a minimum of 2 mg/L DO to maintain higher life forms, most game fish require 4 mg/L. When biodegradable organics such as sewage are discharged into a stream, microorganisms convert the organics into new cells and oxidized waste products. During this process, DO is consumed, resulting in the biochemical oxygen demand (BOD, the amount of oxygen required by aquatic decomposers), another important indicator for water quality.
The oxygensag model was first introduced by Streeter and Phelps in 1925 to predict the oxygen concentration and deficit in a stream caused by the discharge of sewage. Consider a stream with constant flow rate and crosssectional area (and hence velocity), and consider an element of the stream moving with the flow. The change of BOD in this element as function of time, t, can be described by a first order reaction:
dL k_{1}L dt  
(1)
where L is BOD (mg/L) and k_{1} (d^{1}) is the rate coefficient of biochemical decomposition of organic matter. Note, the oxygen concentration does not appear here (the reaction is first order only in L) and we have assumed there is no diffusion of oxygen or L between different elements of the stream (these approximations are needed to make the model tractable). Since BOD depletion is directly related to the oxygen consumption (i.e. we do not take any other oxygen depletion mechanisms into account), the resulting change of the oxygen concentration due to the biochemical reaction is:
dC_{de }dL
k_{1}L (2) dt dt
where C is the oxygen concentration (mg/L) and the index de stands for deoxygenation. The main process replenishing the oxygen content in a stream is reaeration, i.e. the uptake of oxygen from the atmosphere across the water surface due to the turbulent motion of water and molecular diffusion. The driving force for reaeration is the difference between the oxygen saturation concentration, C_{s} (mg/L), and the actual oxygen concentration in the water, C, which is called the oxygen deficit, D = Cs – C:
^{dC}^{re } k_{2 }C_{s }C k_{2}D (3) dt
where k_{2} is the atmospheric oxygen dissolution rate coefficient i.e. the reoxygenation constant (d^{1}). Factors affecting k_{2} include the stream turbulence, the surface area of the stream, the water depth and the temperature. The resulting balance for the oxygen concentration in the water can thus be written as:
dC dCre dCde
k_{1}L k_{2}D (4)
dt dt dt
and the change in the oxygen deficit is given by
dD k_{1}L k_{2}D dt  
(5)
Assuming that the stream flow rate and the flow rate of the waste water effluent are constant and that the mixing of sewage water and the stream water happens instantaneously across the full crosssection of the stream (i.e. all dispersion effects are neglected) a simple mass balance at the mixing points yields for the BOD
L Q L Q
L0 R R SW SW (6)
QR QSW
and for the oxygen concentration
C Q C Q
C0 R R SW SW (7)
QR QSW
where Q (m^{3}/d) is the flow rate and the indices R and SW stand for river and sewage, respectively. Using equation (7), the oxygen deficit at the mixing point can be calculated:
D_{0 } C_{s }C_{0} (8)
The temperature dependence of the rate coefficients and the saturation concentration can be described by equations (9) – (11):
k1 k1T0 1.047TT0  (9) 
k2 k2T01.024TT0  (10) 
Cs 14.126e0.0202T  (11) 
Note that all three equations require that the temperature is given in °C and that k_{1}(T_{0}) means k_{1} at temperature T_{0}.
Integrating equation (5) and determining the maximum of the resulting function gives the maximum oxygen deficit, D_{cr}, and the corresponding time, t_{cr}:
^{L}_{0}ek1tcr
D_{cr } (12)
f
1 f 1D0
tcr ln f 1 (13) k_{1}f 1 L0
where f = k_{2}/k_{1}.
Tasks:
An unpolluted river (L_{R} = 0) flows with a flow rate of Q_{R} = 3.8 m^{3}/s and a stream velocity of v = 0.65 m/s. The water temperature is 20 °C. The water is saturated with oxygen and the reaeration constant for this stream at T_{0} = 20 °C is k_{2} = 0.9 d^{1}. Due to an accident, untreated sewage with a flow rate of Q_{SW} = 0.27 m^{3}/s and a BOD of L_{SW} = 600 mg/L is discharged into the river. The oxygen concentration of the sewage stream, C_{SW}, is zero and the rate coefficient of biochemical decomposition of organic matter in this stream at 20 °C is k_{1} = 0.3 d^{1}.
a) Present the profiles of L, C and D as function of the distance from the discharge point. How far downstream does the maximum oxygen deficit occur and does it present any threat to fish?
b) In a separate plot present the oxygen concentration profile as function of distance from the discharge point for T = 10 °C, 20 °C, 25 °C. Also, calculate the critical oxygen deficit and the corresponding distance. Comment on how these properties change with temperature.
c) Due to draught and flooding, the flow rate of water in the river can vary. To study the influence of the flow rate on the water quality, calculate the critical oxygen deficit and the corresponding distance from the discharge point for Q_{R} = 0.6 m^{3}/s, 3.8 m^{3}/s, 18.7 m^{3}/s. Assume that the crosssection of the river remains constant. Comment on the effect. (Note, you do not need to calculate the concentration profiles for C).
d) The oxygensag model as presented by Streeter and Phelps is the simplest model to model the water quality but it neglects important phenomena that influence the oxygen content of the water. The extended BOD model also takes into account the sedimentation of biodegradable organic matter, benthic oxygen demand (i.e. BOD due to the decay of organic matter in sludge deposits) and the net creation of oxygen by aquatic plants.
For this model, equation (1) becomes:
dL k_{1}Lk_{3}L B dt  
(14)
where k_{3} is rate coefficient for BOD removal by sedimentation (k_{3} = 0.3 d^{1}) and B is the benthic oxygen demand (B = 1.5 mg/(L d)). To take the net production of oxygen due to photosynthesis into account, equation (5) has to be modified to
dD k_{1}L k_{2}D P dt  
(15)
where P is the rate of oxygen addition to the water by the photosynthetic activity of aquatic plants and algae. Here, P = 0.6 mg/(L d). Compare the oxygen concentration profiles as function of distance from the discharge point for the oxygen sagmodel and the extended BOD model and explain your observations. Use the same function file as for the other tasks to answer question d).
e) What other factors should be taken into account in more realistic models?
For this handin exercise it is required that you create a script file from which a MATLAB report can be produced. A worksheet on publishing reports plus several videos is available on Learn in the folder for the handin exercise. Your published report should contain the code as well as your comments, discussion and plots. Note that only your script file needs to be published as a report, not the function file, but you must submit your script and function file.
References:
Streeter N. W., Phelps E. B. (1925); A Study of the Pollutions and Natural Purification of the Ohio River. Public Health Bulletin NO. 146. U.S. Public Health Service
The deadline for the handin exercise is Monday, 5/12/2016, 16:00.
Please submit the following files in the dropbox on Learn:
• Your script file formatted in such a way that a report in MATLAB can be produced as a PDF file. Do not submit the PDF report itself. We will use your MATLAB files to create the report.
• Your function file containing the rhs of the ODEs. Note that you only need one function file to describe both models.
If you experience any problems accessing Learn, please email your files directly to ... before the deadline.
Marking scheme (detailed marking criteria are available on Learn)
Readability of code (headers, variable dictionaries, use of comments, use of white space, …)
 15 %  
Format of the report (graphs, units, output, use of (advanced) report writing features, …)  15 %  



 
Base case (answer to a) )



 25 % 
Influence of temperature (answer to b) )



 15 % 
Influence of flow rate (answer to c) )



 10 % 
Extended model (answer to d) )



 15 % 
Discussion of factors not included in the models (answer to e) ) 


 5 % 