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

Projectile motion with air resistance

  1. Sep 26, 2013 #1
    How do we write differential equations for projectile motion in 2 dimensions featuring air resistance of magnitude kv^2, acting directly opposite to the direction of motion at that moment in time, where v is the velocity in the direction of motion at that moment in time?
     
  2. jcsd
  3. Sep 26, 2013 #2

    tiny-tim

    User Avatar
    Science Advisor
    Homework Helper

    Hi Big-Daddy! :wink:

    Tell us what you think, and why, and then we'll comment! :smile:
     
  4. Sep 26, 2013 #3
    Obviously I'm asking because I don't know. More pertinently I can't imagine how to model this properly. Presumably we'll need independence on the x and y axes, connected with time.

    For free-fall I could model it pretty easily:

    [tex]mg - kv^2 = m \cdot \frac{dv}{dt}[/tex]

    But this doesn't seem to be remotely of the same difficulty. Velocity in free fall is always in the same direction as acceleration, but in the projectile motion case, the velocity is defined by the initial projection angle and velocity (which would be given of course, along with the values of m, g and k).
     
  5. Sep 26, 2013 #4

    Nugatory

    User Avatar

    Staff: Mentor

    Try writing the differential equation in terms of vectors [itex]\vec{v}[/itex] and [itex]\vec{F}=m\vec{a}[/itex]. That will at least get the problem modeled properly, and you can decompose it into coupled differential equations for ##x(t)## and ##y(t)##.

    The initial angle and speed provide the boundary conditions you'll need to determine the arbitrary constants that show up in the solutions of the differential equations.

    (And I have to caution you that solving these equations is a non-trivial problem).
     
    Last edited: Sep 27, 2013
  6. Sep 26, 2013 #5
    Well ok, perhaps this is an analogous equation to mine for free-fall, where v has now been replaced by a vector?

    Question is, how to decompose this into my x(t) and y(t) differential equations?
     
  7. Sep 26, 2013 #6

    jhae2.718

    User Avatar
    Gold Member

    Equate the components of the acceleration and force vectors to obtain the scalar equations of motion.
     
  8. Sep 27, 2013 #7

    tiny-tim

    User Avatar
    Science Advisor
    Homework Helper

    Hi Big-Daddy! :smile:

    (just got up :zzz:)
    ok, same, but with vectors …

    [itex]-mg\mathbf{y} - kf(\mathbf{v})\mathbf{v} = m \cdot \frac{d\mathbf{v}}{dt}[/itex]
    where ##f(\mathbf{v})## = … ? :wink:
     
  9. Sep 27, 2013 #8
    What do you mean?
     
  10. Sep 27, 2013 #9
    I'm not sure I understand what that y is doing there. As for f(v), not sure ... looks like it should just be v to me, but I'm not sure why you wrote it as such then ...
     
  11. Sep 27, 2013 #10

    arildno

    User Avatar
    Science Advisor
    Homework Helper
    Gold Member
    Dearly Missed

    f is a scalar, non-negative function.
    What ought f to be then?
     
  12. Sep 27, 2013 #11

    tiny-tim

    User Avatar
    Science Advisor
    Homework Helper

    because gravity is mg directly downwards,

    so the vector for the force of gravity is -mg in the y direction, ie -mg time the unit vector in the y direction, ie -mgy :wink:
    yes, but you need to write it in terms of the vector v, so it's -(k√(v2))v :smile:
     
  13. Sep 28, 2013 #12

    Redbelly98

    User Avatar
    Staff Emeritus
    Science Advisor
    Homework Helper

    jhae means Fx = m ax , and similarly for the y-component. That's two equations of motion, one for each component.

    As for the original question:

    The air resistance force has a magnitude k v2, and direction opposite to that of v. For the x-component of the force, multiply this magnitude by the cosine of the angle it makes w.r.t. the +x-direction -- this is -vx/v -- and this gives you the x-component of the force due to air-resistance.

    Do the same for the y-component.

    And then you'll have the force components due to air resistance to use in the equations relating Fx to ax and Fy to ay.

    p.s. This is worth repeating:
     
  14. Sep 28, 2013 #13
    Ok, I don't know how to include the angles - surely with 2 equations we can only afford to have 2 variables, x and y?

    So far I have

    [tex]F_x = m \cdot \frac{d^2x}{dt^2} = -k \cdot \frac{dx}{dt}[/tex]

    and

    [tex]F_y = m \cdot \frac{d^2x}{dt^2} = -mg \cdot cos(\theta) - k \cdot \frac{dx}{dt}[/tex]

    I think there is a problem in how I have resolved the weight though. How can I do this better?
     
  15. Sep 28, 2013 #14

    Redbelly98

    User Avatar
    Staff Emeritus
    Science Advisor
    Homework Helper

    You're getting there :smile:. Yes, there are just the two variables you mentioned to concern yourself with.

    A couple of problems to clear up in your equations:

    1. The weight is -mg. There is no cosine term involved in the weight, since it always acts downward, in the -y-direction, regardless of the angle of the trajectory.

    2. Also, you have incorrectly multiplied kv2 and -vx/v -- so you are missing a factor of v in the air resistance expression. (And the same error occurs in the y-equation).
     
  16. Sep 28, 2013 #15

    arildno

    User Avatar
    Science Advisor
    Homework Helper
    Gold Member
    Dearly Missed

    You won't be able to solve this analytically, but it is a good exercise to set up the equations of motion that will govern the system.
     
  17. Sep 28, 2013 #16
    Ok then so maybe

    [tex]F_y = m \cdot \frac{d^2y}{dt^2} = -mg - k \cdot \frac{dy}{dt}[/tex]

    I'm not sure I understand ... by v, do you mean the resultant of vx and vy, and by vx you mean dx/dt? If so then maybe:

    [tex]F_y = m \cdot \frac{d^2y}{dt^2} = -mg - k \cdot \frac{dy}{dt} \cdot ((\frac{dy}{dt})^2+(\frac{dx}{dt})^2)^{1/2}[/tex]

    and

    [tex]F_x = m \cdot \frac{d^2x}{dt^2} = - k \cdot \frac{dy}{dt} \cdot ((\frac{dy}{dt})^2+(\frac{dx}{dt})^2)^{1/2}[/tex]
     
  18. Sep 29, 2013 #17

    arildno

    User Avatar
    Science Advisor
    Homework Helper
    Gold Member
    Dearly Missed

    In your last last line, the "dy/dt" outside the root should be replaced with "dx/dt"
     
  19. Sep 29, 2013 #18

    Redbelly98

    User Avatar
    Staff Emeritus
    Science Advisor
    Homework Helper

    Yes to both.

    Apart from arildno's correction, you got it.
     
  20. Sep 30, 2013 #19
    Thank you.

    And the boundary conditions would involve me specifying initial y and x displacement from the origin as well as initial x velocity and y velocity (which can be found as the cos and sin components respectively of the initial total magnitude of velocity, making sure the y initial velocity is positive if the point is travelling upwards initially and negative if it is travelling downwards initially), and nothing else?
     
  21. Sep 30, 2013 #20

    Nugatory

    User Avatar

    Staff: Mentor

    Yes.
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook




Similar Discussions: Projectile motion with air resistance
Loading...