MATH 650: Mathematical Modeling Fall, 2013 - Final Quiz
MATH 650: Mathematical Modeling Fall, 2013 - Final Quiz
MATH 650: Mathematical Modeling Fall, 2013 - Final Quiz
Instructions:
This final quiz contains no new concepts or solution methods. It is intended to
provide an overview of the course, including both techniques and some (hopefully)
interesting applications. The emphasis is on using known results in new situations,
and interpreting graphical and analytical information about a given application.
Feel free to discuss the problems with one another, and to make use of any other
resources which help you to work your way through them. However, your submis-
sions must be your own, and must reference any source (human or otherwise) from
whom you received help.
Presentation is important. Please write your solutions in clear sentences which
convey your reasoning. Remember that I cant know what youre thinking...I can
only know what you tell me. Conciseness is also very much appreciated.
Several of the problems use specific Maple files which will be made available. (Im
assuming you can do simple Maple commands such as plot and implicitplot on your
own; email me if you need help.)
I suggest you spend no more than 10 hours in total on this assignment,
including writing up your solutions. While you might not complete all the
problems youd like to, give it your best try...but dont worry about it beyond that
time span. Partially completed problems are acceptable; Ive included the marks
for the individual parts of each problem.
There are 7 problems on this assignment, two on first order DEs, one on dimensional
analysis, and four on second order DEs and/or systems of two first order DEs.
There are no extra problems.
Submit your solutions in the usual manner in LEARN.
2
Problem Marks
1 15 (a. 8, b. 7)
2 11 (a. 4, b. 7)
3 15 (a .10, b. 5)
4 11 (a. 5, b. 4, c. 2)
5 11 (a. 8, b. 1, c. 2)
6 17 (a. 1, b. 3, c. 3, d. 3, e. 7)
7 17 (a. 2, b. 7, c. 7, d. 1)
97
3
b. We have explored the logistic model for the growth of a population of size N (t),
dN N
namely, = rN (1 ) , N (0) = N0 . Suppose that the underlying logistic
dt K
growth has r = 0.1 per year and K = 1000, and the population N (t) is subject to
seasonal harvesting, with period 1 year. A suitable model is the IVP
dN N
= 0.1N (1 ) A(1 + sin(2t) , N (0) = N0 ,
dt 1000
where A is a constant which determines the maximum amount (2A) of harvesting
which occurs.
i. Use the plot command given in the Maple file Math650FQStudentPlots to explore
the behaviour of the solutions for various values of A between 5 and 30. Observe
carefully how the long term behaviour of the solutions depends on both A and
on N0 ; perhaps make a table of your observations for different values of A.
4
ii. Do you get the same longterm behaviour if the intrinsic growth rate is r = 0.2
per year? Does the dependence on N0 change? Explain.
iii. Write a paragraph about your results, considering for what values of A the
population is viable for at least some initial populations N0 , how the longterm
behaviour of the population changes as A increases, and how this depends on
N0 , and on r . If you were a conservation officer, what recommendation(s)
would you make to preserve this population?
(Submit NO plots for this question...just your well-organized thoughts and
conclusions.)
dv
i. Apply Newtons Law to give the IVP m = mg F, v(0) = v ? , solve for v(t), and
dt
hence show that his velocity during this second stage is given by
t2
v(t) = gt c + v ? .
2
Then find the time tf when his velocity v(tf ) = 0 . (Your answer will be in terms
of g, v ? , and c .)
5
ii. Ideally, we want his velocity v = 0 just as he touches down at s = 66 m. Show that
this will occur if the following equation is satisfied:
g tf 2 c tf 3
66 = + v ? tf .
2 6
Then use the commands for Problem 2 in the Maple file Math650FQStudentPlots
to solve this equation for the value of the constant c , and hence to find the time
tf when he lands.
What is the total time t? + tf for his fall?
If his mass is 80 kg, what is the maximum force F on the cable during Stage
2?
iii. Comment on any aspects of the physical situation that have been ignored in the
above model for this situation. Suggest an alternative model if possible.
Dimensional Analysis
3a. The driver of a car of mass m travelling at constant velocity v0 spots a squirrel
crossing the road. She brakes with a constant force F for a distance s (in a straight
line), reducing her speed to v1 to avoid hitting the squirrel.
i. Write down the dimensional matrix, find its rank, and apply the Pi Theorem to
mv02 v0
show that s = f ( ) . Explain why this answer makes sense, physically.
F v1
ii. Two simple integrations of Newtons Law for this problem reveal that, during
F F t
the braking period, v(t) = t + v0 , and s(t) = t ( + v0 ) . Using these
m m2
results (no need to prove them), show that
1
s F = (mv02 mv12 ) .
2
Explain this result in terms of the work done by the force F and the kinetic
energy of the car. Then rearrange the equation to the form in part a., and
v0
hence find the function f ( ) .
v1
3b. An explosion can be thought of as a spherical shock wave emanating from a point.
The radius R of the shock wave depends on the energy E released by the explosion,
the ambient air density , and the elapsed time t. Suppose a scientist tells you that
these quantities are related by the equation
1 E 2
ln R = ln + ln t + ln k , where k is a dimensionless constant.
5 5
Does dimensional analysis support this claim?
[In 1941, a similar dimensional analysis predicted the energy of an atomic bomb.
Later, data from the first atomic bomb test in New Mexico verified its accuracy.
Reference: Mathematical Models and Their Analysis, by F.Y.M.Wan, 1989.]
6
where k1 and k2 are positive constants, and f (t) is the inflation rate. Consider
a periodic inflation rate f (t) = f0 cos t . Show that, with suitable definitions
of 02 and A, this DE is equivalent to the above model, and describe what
happens economically, according to this model, if = 0 .
x00 + 2 x0 + 02 x = f (t) ,
0
Prove that both these conditions hold if the constans and 0 satisfy < < 0 .
2
i. What will guarantee that condition 1. is met? (Examine your results from part
a.)
ii. We have already seen (in the course lectures) that the steady-state amplitude
of the particular solution of a DE of the form x00 + 2 x0 + 02 x = A cos t is
given by A
Ass = 2 p ,
0 (1 r2 )2 + 4 2 r2
where r = , and = . You also showed on a previous assignment that
0 0
1
amplitude resonance cannot occur for any > 0 if > . Rewrite this
2
condition in terms of and 0 to determine what will guarantee condition 2.,
completing the proof.
Two-Species Interactions
6. Consider the following model for the interaction between two species.
dN1 N1 dN2 N2
= r1 N1 (1 ) aN1 N2 , = r2 N2 (1 ) + bN1 N2 .
dt K1 dt K2
a. Clearly the underlying growth of each species is logistic, i.e., in the absence
of one another, they would each thrive. What do the two interaction terms
tell you about the nature of their interaction, as compared to the competition
model we have previously discussed? Name this interaction, giving reasons.
N1 N2
b. Show that, with the choice of dimensionless variables x = K 1
,y=K 2
, = r1 t ,
and dimensionless constants A = aK bK1 r2
r1 , B = r2 , R = r1 , the model becomes
2
dx dy
= x(1 x Ay) , = Ry(1 y + Bx) .
d d
Give biological interpretations of x and y, and of the constants A, B, and R.
c. Show that the system in part b. has four equilibrium points, with the one in
1A 1+B
the positive quadrant being P : , , provided A < 1 .
1 + AB 1 + AB
d. Suppose that the equilibrium P is asymptotically stable for all x(0) > 0 and
y(0) > 0 . Explain why, if 0 < A < 1 and B > 0, the longterm population
of species N1 will always be below K1 , while that of species N2 will always be
above K2 .
e. Use the plot command for Problem 6 in the Maple file Math650FQStudentPlots
to plot phase portraits for each of the following cases:
9
1 1 6 1
i. A = , B = 1, R = ; ii. A = , B = 1, R = .
2 8 5 8
In each case, linearize about the asymptotically stable equilibrium to confirm
the type of equilibrium you observe in the phase portrait. State the the bio-
logical outcome, and justify it, based on the parameter values for A and B.
d. It seems clear from your results in b. and c. that the qualitative behaviour
of this system changes dramatically for some value of n between 2 and 2.5.
Describe the change that occurs in the stability of the equilibrium (xe , ye ),
and the resulting change in the phase portrait. Which example in the lectures
on the Phase Portrait Gallery in Module 8, Section 1.1 does our system most
resemble? Now go to Wikipedia and read about the Hopf bifurcation (the
animation of the Brusselator is especially intriguing).
[Reference for Problem 7: Mathematical Modeling in Systems Biology, by Brian
P. Ingalls, MIT Press, 2013, pages 108-110.]