Reference no: EM132355011
Question
Use triangulation to locate a pollutant source
In LAB 10.25, you developed functions for using three points and two distances to calculate the location of an unknown point.
Now, let's suppose that the unknown point is actually a point source for some pollutant X, and at time 0 it releases a burst of X. Assume no wind, and that Xdiffuses freely from its release point.
Suppose that we have measurement devices at the three points 1, 2, and 3 that track the concentration of X as a function of time, and that the concentration of X is detected to peak at times t1, t2, and t3, at each point, respectively. Can we determine the location of the evil point source?
Let's make the assumption that the time intervals scale with the square of the distance divided by the diffusivity of X, ti ∼ di^2
****the i's are subscripts*****
This means we can invoke some unknown constant K = √DX that satisfies (** the X is a subscript)
d1 = K√t1, d2 = K√t2, and d3 = K√t3. Can we estimate K uniquely knowing this?
Consider that if we choose a value for K, all three distances are set, and we get the two solutions points that correspond to that value of K. We can then measure the distance to point 3 for each solution point, and see if it matches K(t03.5). What if it doesn't? Is there one unique value of K that works?
The answer is yes! Consider the following algorithm:
. Given times t1, t2, t3 and the points (x1,y1), (x2,y2), and (x3,y3);
. Guess K.
. Find the point (x,y) the correct distance from both (x1,y1) and (x2,y2) and closest to (x3,y3).
. Compute the apparent distance d'3 from (x,y) to (x3,y3). DX
. Compute z = √K√t3 .d′3
. Ifzis1,stop;otherwiseletKbeK/z,andgoto3.
The value of K this converges to, if squared, gives an estimate for the diffusivity DX.
Part 1 Objective
To the following code, add your implementation of the algorithm listed above. Then, given the following coordinates and times of observed peak concentration, print the coordinates of the evil point source in the format (x,y) (including the parentheses) where x and y are floats with three digits following the decimal.
Coordinates x,y of observation. Time of peak concentration
(-8.0,5.0) 300
(-5.0,-6.0) 400
(9.0,1.0) 500