Problem 60406. Alert a city about a spill
Problem statement
Cody Problem 54750 involved determining the length of a stream affected by a spill of a contaminant. Any municipalities within that reach would want to know when water from the river would be safe to drink—for example, below the U.S. Environmental Protection Agency’s maximum contaminant level (MCL). As in CP 54750, the spill of mass M will be assumed instantaneous at position
and time
and mixed over the cross section (with area A). Then if the flow is steady and the geometry of the flow does not change downstream, the concentration can be computed with
where
is the mean velocity of the river, Q is the discharge or volumetric flow rate, and K is a dispersion coefficient, which describes spreading by several mechanisms.
Write a function that returns the dates and times (given as datetimes) between which the water is unsafe to drink (i.e., the concentration exceeds the MCL) at position x downstream of the spill. Round the times to the nearest second. If the concentration does not exceed the MCL, return 'The MCL is not exceeded.' Please note that the MCL is given in mg/L, whereas other variables are given in SI units.
Details
Many empirical formulas are available for the dispersion coefficient. For this problem, use the formula of Seo and Cheong (1998):
where B is the width of the channel (assumed rectangular here), H is the water depth, and
is the shear velocity, which is related to the shear stress on the wetted perimeter of the channel. In steady uniform flow, the component of the fluid’s weight down the slope will balance the friction on the channel bed, and the shear velocity can be computed as
where
is the gravitational acceleration,
is the longitudinal slope of the channel,
is the hydraulic radius, and P is the wetted perimeter. For a rectangular channel,
.
In addition to assuming steady uniform flow and an unchanging channel, ignore any reaction, decay, or loss of the chemical; this assumption provides a conservative estimate of the time range. In practice, one would include a factor of safety that accounts for uncertainty in the parameters. Nevertheless, the calculations here would form a basis for the advice to the city.
Solution Stats
Problem Comments
-
1 Comment
William
on 1 Jun 2024
Silly me. I thought the "liter" was an SI unit, but it's not.
Solution Comments
Show commentsProblem Recent Solvers4
Suggested Problems
-
Return the largest number that is adjacent to a zero
5315 Solvers
-
Remove the polynomials that have positive real elements of their roots.
1669 Solvers
-
Given an unsigned integer x, find the largest y by rearranging the bits in x
1800 Solvers
-
600 Solvers
-
Detect a number and replace with two NaN's
190 Solvers
More from this Author271
Problem Tags
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!