Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Displacement with Air Drag

  1. Jan 29, 2009 #1
    This is not a homework problem, just an idea I had.

    If we have a sphere which is being accelerated in space by a force [tex]\vec{F}[/tex] and we take the drag into account, what would the displacement be after t seconds? If the drag force is [tex]\vec{D}=-c_1v-c_2v^2[/tex], then the acceleration would be [tex]\vec{a}=\frac{\vec{F}+\vec{D}}{m} \,,\, a=\frac{F-c_1v-c_2v^2}{m}[/tex] such that the displacement [tex]s=\frac{1}{2}\left(\frac{F-c_1v-c_2v^2}{m}\right)t^2[/tex]. Is this correct? Could someone help me solve this differential equation, if it has a solution?

    Thanks for any help.
  2. jcsd
  3. Jan 29, 2009 #2
    No, that's not the correct formula for the displacement. The well-known [tex]s=\frac{1}{2}at^{2}[/tex] holds only if the acceleration is constant. In this case, the acceleration is a quadratic function of velocity, so we have one scalar differential equation:
    where A= F/m, B=C1/m and C=C2/m are constants. The equation can be solved using separation of variables, as
    [tex](1) \int dt=\int\frac{dv}{A-Bv-Cv^{2}}[/tex]

    I can solve this easily if either B or C is equal to zero.
    • [tex]D=-c_{1}v[/tex]
      Then by solving (1), and then another integration with respect to time, we get (for v(0) = 0)
      where s(t) is the displacement after time t.
    • [tex]D=-c_{2}v^{2}[/tex]

      is the terminal velocity (where D=F). Then using partial fraction decomposition on (1), integration and then another integration with recpect to time, we get (again for zero velocity at the begining):
      [tex]s(t)=\frac{v_{m}m}{\sqrt{Fc_{2}}}ln cosh\left(t\frac{\sqrt{Fc_{2}}}{m}\right)[/tex]
      where s(t) is the displacement at time t.

    I'm too lazy at the moment to solve the problem for generally non-zero C1 and non-zero C2, which should be possible, but some difficulties might be encountered when inverting the function. I hope you got the picture how it works. Also, you can use numerical solution of the equation, which is the common approach in physics.

    And finally, I apologize for any inconvenience caused by my English or incorrect English math terms.
  4. Jan 31, 2009 #3
    I solved the integral to

    [tex]\int \frac{dv}{A-Bv-Cv^2}=-\frac{2\tan^{-1}\left(\frac{B+2Cv}{\sqrt{-B^2-4AC}}\right)}{\sqrt{-B^2-4AC}}[/tex], so then [tex]t=-\frac{2\tan^{-1}\left(\frac{B+2Cv}{\sqrt{-B^2-4AC}}\right)}{\sqrt{-B^2-4AC}}[/tex]

    Working with this, I get that


    If this is correct, I want to integrate it in order to get an expression for the displacement over time.

    [tex]s=\int \frac{\sqrt{-B^2-4AC}}{2C}\cdot\tan\left(\frac{\sqrt{-B^2-4AC}}{2}t\right)-\frac{B}{2C}\rm{d}t[/tex]

    I will simplify it be setting [tex]\frac{\sqrt{-B^2-4AC}}{2}=k_1[/tex] and get

    [tex]s=\frac{k_1}{C}\cdot\int\tan\left(k_1\cdot t\right)\rm{d}t-\int\frac{B}{2C}\rm{d}t=\frac{k_1}{C}\left(-\frac{\log\left(\cos\left(k_1t\right)\right)}{k_1}\right)-\frac{B}{2C}t+s_0[/tex]

    How does this look?
  5. Feb 1, 2009 #4
    By [tex]tan^{-1}[/tex] you mean [tex]tanh^{-1}[/tex], right? Also note that the expression [tex]\sqrt{-B^2-4AC}[/tex] would be imaginary, since constants A, B and C are positive.
    My solution of the integral is:


    which could be written like yours considering the fact that:

    The rest seems fine, it is even consistent with my solution for B=0.
    Last edited: Feb 1, 2009
  6. Feb 1, 2009 #5


    User Avatar
    Science Advisor
    Homework Helper
    Gold Member
    Dearly Missed

    Do remember, though, that the proper quadratic law for the air drag is [tex]-cv|v|[/tex], rather than [tex]-cv^{2}[/tex] for some constant c.
  7. Feb 1, 2009 #6
    But since the drag force acts in opposite direction to the velocity, for an accelerating force F we can use scalar equations with [tex]F_{d}=-cv^{2}[/tex]. So the above calculations are correct.
  8. Feb 1, 2009 #7


    User Avatar
    Science Advisor
    Homework Helper
    Gold Member
    Dearly Missed

    Well, that would depend upon the problem you were to model with this.
    If you think of, say, launched projectile problem within a constant gravity field, then the general form is necessary.

    The force law on the way up (velocity "positive") wll go as -cv^2, whereas on its way down (velocity "negative") as +cv^2.
  9. Feb 1, 2009 #8
    The better general formula then would be
    where [tex]\vec{\hat{v}}[/tex] is a unit vector codirectional with velocity.
    But anyway, with some common sense you can always decide the correct direction of the drag force, so it doesn't accelerate the projectile.
    Simple [tex]}-cv^{2}[/tex] is a much better function to calculate with, than some expression with an absolute value.
    I agree with you and understand your point. Thanks for the comment.
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook