1. Not finding help here? Sign up for a free 30min tutor trial with Chegg Tutors
    Dismiss Notice
Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Autonomous differential equation

  1. Mar 11, 2014 #1
    1. The problem statement, all variables and given/known data

    A pond forms as water collects in a conical depression of radius a and depth h. Suppose water flows in at a constant rate, k and is lost through evaporation at a rate proportional to the surface area.


    I was wondering whether anyone could give me some guidance on this part of the question:

    Find the equilibrium depth of water in the pool. Is the equilibrium asymptotically stable?


    2. Relevant equations

    [tex]\frac{dV}{dt}=k-\alpha \pi(\frac{3aV}{\pi h})^{\frac{2}{3}}[/tex]

    Where [tex]\alpha[/tex] = coefficient of evaporation.

    radius at any time r(t) and depth at any time l(t).

    [tex]\frac{a}{h}=\frac{r}{l}[/tex]


    3. The attempt at a solution

    [tex]\frac{dV}{dt}=0=k-\alpha \pi(\frac{3aV}{\pi h})^{\frac{2}{3}}[/tex]


    Putting the differential equation only in terms of radius and alpha:
    [tex]\alpha \pi r^2 =k[/tex]

    [tex]r=\sqrt{(\frac{k}{\alpha \pi})}[/tex]

    given that [tex] l=r\cdot\frac{h}{a}[/tex],

    [tex]l=\frac{h}{a}\sqrt{(\frac{k}{\alpha \pi})}[/tex]

    But I am not sure how to find whether or not this is asymptotically stable.
     
  2. jcsd
  3. Mar 11, 2014 #2

    maajdl

    User Avatar
    Gold Member

    Consider initial conditions near the equilibrium point.
    If initially the volume is a little bit larger than the equilibrium volume, does the volume come back to equilibrium, or does it increase without end?
    And the reverse, if the initial volume is a little bit smaller than the equilibrium volume?
    You could also linearize the equation around the equilibrium point and answer the same questions, easily.
     
    Last edited: Mar 11, 2014
  4. Mar 11, 2014 #3
    Thanks Maajdl.

    I can see that solving for V, the equilibrium volume is [tex]V=\frac{\pi h}{3a}(\frac{k}{\alpha \pi})^{\frac{3}{2}}[/tex]

    But without picking some arbitrary values and plotting a direction field (or plotting Dv against V) I am not sure how to tell whether the volume diverges or converges to equilibrium.

    (I haven't yet studied linearisation so am not quite ready to use this option)
     
  5. Mar 11, 2014 #4

    maajdl

    User Avatar
    Gold Member

    Assume Vo is the equilibrium volume.
    If V>Vo, will the evaporation overcome the water supply and bring V back to Vo?
    If V<Vo, will the water supply overcome the evaporation to bring V back to Vo?

    This site may inspire you:

    http://tutorial.math.lamar.edu/Classes/DE/EquilibriumSolutions.aspx

    The situations described there are a little bit more complicated, but the same reasoning will apply.
    The pictures shown there are useful to visualize what will happen.
    Such "phase space flows" are more general than this one variable problem, and are at the hearth of the stability theory of differential equations. (and are even related to top mathematical physics, like the KAM theorem)
     
  6. Mar 11, 2014 #5
    Thanks a lot maajdl... I'll have a good look through that webpage too. I need to get to grips with this topic!
     
  7. Mar 11, 2014 #6

    pasmith

    User Avatar
    Homework Helper

    It is much, much easier to look at [itex]dl/dt[/itex] (as you are strongly hinted to do by the instruction to find the equilibrium depth, not the equilibrium volume). Since
    [tex]
    r = \frac{al}{h}
    [/tex]
    we have
    [tex]
    V = \frac{\pi}{3}r^2 l = \frac{\pi a^2}{3h^2} l^3
    [/tex]
    and hence
    [tex]
    \frac{dV}{dt} = \frac{\pi a^2}{h^2} l^2 \frac{dl}{dt}.
    [/tex]
    Writing [itex]k = \frac{\pi a^2}{h^2}\alpha C^2[/itex] and using [itex]r^2 = a^2l^2/h^2[/itex] in the ODE gives
    [tex]
    \frac{dV}{dt} = \frac{\pi a^2}{h^2}\alpha C^2 - \alpha \frac{\pi a^2}{h^2} l^2
    = \frac{\pi a^2}{h^2} \alpha\left(C^2 - l^2\right)
    [/tex]
    so that
    [tex]
    \frac{dl}{dt} = \alpha\left(\frac {C^2}{l^2} - 1\right)
    [/tex]
    where, of course,
    [tex]
    C = \frac{h}{a} \sqrt{\frac{k}{\pi \alpha}}.
    [/tex]

    No linearization is necessary here; just sketch a graph of [itex]dl/dt[/itex] against [itex]l[/itex], and see what the sign of [itex]dl/dt[/itex] is immediately above and below the fixed point.
     
  8. Mar 11, 2014 #7

    maajdl

    User Avatar
    Gold Member

    It is also possible to sketch the graph of dV/dt against V .
     
  9. Mar 11, 2014 #8
    Thanks, Pasmith. After studying it, I can see what you have done. The thing that concerns me is that it wouldn't have occurred to me to make these steps. You must be thinking several steps ahead, visualizing how it will all fit together...
     
  10. Mar 11, 2014 #9
    That's true... I would quite like to try to see if I can get it symbolically though.
     
  11. Mar 11, 2014 #10

    maajdl

    User Avatar
    Gold Member

    All you need is exercise.
    These steps are very standard and will see them again and again, specially in classical mechanics.
    Choosing l instead of V as variable would be highly natural if you would like to study if overflow might occur!
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook

Have something to add?
Draft saved Draft deleted



Similar Discussions: Autonomous differential equation
Loading...