Let us worry about your assignment instead!

We Helped With This Engineering Assignment: Have A Similar One?

More InfoEngineering Physics Homework Help

Short Assignment Requirements

It is a computing assignment for Chemical Engineering.

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 oxygen-sag 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 cross-sectional 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:







where L is BOD (mg/L) and k1 (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:


                                                                                  dCde           dL

 k1L            (2) dt   dt


where C is the oxygen concentration (mg/L) and the index de stands for de-oxygenation. 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, Cs (mg/L), and the actual oxygen concentration in the water, C, which is called the oxygen deficit, D = Cs – C:


dCre k2 Cs Ck2D    (3) dt


where k2 is the atmospheric oxygen dissolution rate coefficient i.e. the re-oxygenation constant (d-1). Factors affecting k2 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

                                                                                                k1Lk2D                                                             (4)

                                                                    dt        dt         dt

and the change in the oxygen deficit is given by



k1L k2D




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 cross-section 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)



and for the oxygen concentration


                                                                                           C Q C      Q

                                                                            C0 R                                                                                                                                             R    SW                                                                                                                                                     SW     (7)



where Q (m3/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:


                                                                                         D0 Cs C0                                                                            (8)


The temperature dependence of the rate coefficients and the saturation concentration can be described by equations (9) – (11):


k1 k1T0 1.047TT0


k2 k2T01.024TT0


Cs 14.126e0.0202T



Note that all three equations require that the temperature is given in °C and that k1(T0) means k1 at temperature T0.


Integrating equation (5) and determining the maximum of the resulting function gives the maximum oxygen deficit, Dcr, and the corresponding time, tcr:


                                                                                       Dcr                                                                          (12)


                                                                                  1              f 1D0 

tcr lnf 1 (13) k1f 1  L0 

where f = k2/k1.



An unpolluted river (LR = 0) flows with a flow rate of QR = 3.8 m3/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 T0 = 20 °C is k2 = 0.9 d-1. Due to an accident, untreated sewage with a flow rate of QSW = 0.27 m3/s and a BOD of LSW = 600 mg/L is discharged into the river. The oxygen concentration of the sewage stream, CSW, is zero and the rate coefficient of biochemical decomposition of organic matter in this stream at 20 °C is k1 = 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 QR = 0.6 m3/s, 3.8 m3/s, 18.7 m3/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 oxygen-sag 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:







where k3 is rate coefficient for BOD removal by sedimentation (k3 = 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 



k1L k2D P




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 sag-model 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 hand-in 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 hand-in 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. 



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 hand-in exercise is Monday, 5/12/2016, 16:00.


Please submit the following files in the drop-box 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 e-mail 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 %


Customer Feedback

"Thanks for explanations after the assignment was already completed... Emily is such a nice tutor! "

Order #13073

Find Us On

soc fb soc insta

Paypal supported