Acessibilidade / Reportar erro

Minimum and terminal velocities in projectile motion

Abstracts

The motion of a projectile with horizontal initial velocity V0, moving under the action of the gravitational field and a drag force is studied analytically. As it is well known, the projectile reaches a terminal velocity Vterm. There is a curious result concerning the minimum speed Vmin; it turns out that the minimum velocity is lower than the terminal one if V0 > Vterm and is lower than the initial one if V0 < Vterm. These results show that the velocity is not a monotonous function. If the initial speed is not horizontal, there is an angle range where the velocity shows the same behavior mentioned previously. Out of that range, the velocity is a monotonous function. These results comes out from numerical simulations.

projectile motion; drag force; Euler equation; terminal velocity


O movimento de um projétil com velocidade horizontal inicial V0, movendo-se sob ação do campo gravitacional e de uma força de arrasto é estudado analiticamente. Como sabemos, o projétil alcança uma velocidade terminal Vterm. Existe um resultado curioso sobre a velocidade mínima Vmin; acontece que a velocidade mínima é menor que a velocidade terminal se V0 > Vterm e é menor do que a velocidade inicial se V0 < Vterm. Estes resultados mostram que a velocidade não é uma função monotônica.

movimento de projétil; força de arrasto; equação de Euler; velocidade terminal


ARTIGOS GERAIS

Minimum and terminal velocities in projectile motion

E.N. Miranda; S. Nikolskaya; R. Riba

Facultad de Ingeniería, Universidad de Mendoza, Mendoza, Argentina

Endereço para correspondência Endereço para correspondência E. N. Miranda E-mail: emiranda@lab.cricyt.edu.ar.

ABSTRACT

The motion of a projectile with horizontal initial velocity V0, moving under the action of the gravitational field and a drag force is studied analytically. As it is well known, the projectile reaches a terminal velocity Vterm. There is a curious result concerning the minimum speed Vmin; it turns out that the minimum velocity is lower than the terminal one if V0 > Vterm and is lower than the initial one if V0 < Vterm. These results show that the velocity is not a monotonous function. If the initial speed is not horizontal, there is an angle range where the velocity shows the same behavior mentioned previously. Out of that range, the velocity is a monotonous function. These results comes out from numerical simulations.

Keywords: projectile motion, drag force, Euler equation, terminal velocity.

RESUMO

O movimento de um projétil com velocidade horizontal inicial V0, movendo-se sob ação do campo gravitacional e de uma força de arrasto é estudado analiticamente. Como sabemos, o projétil alcança uma velocidade terminal Vterm. Existe um resultado curioso sobre a velocidade mínima Vmin; acontece que a velocidade mínima é menor que a velocidade terminal se V0 > Vterm e é menor do que a velocidade inicial se V0 < Vterm. Estes resultados mostram que a velocidade não é uma função monotônica.

Palavras-chave: movimento de projétil, força de arrasto, equação de Euler, velocidade terminal.

In this article we deal with a simple problem, namely to determine the velocity V of a projectile launched with horizontal velocity in a medium with a drag force proportional to a power n of the body speed. The motivation to deal with such a problem is pedagogical. On one hand, it fits well in an intermediate course on classical mechanics [1] or even in a general physics course with calculus [2]. On the other hand, it could be use as a physical example of the Bernoulli equation in a mathematical physics course [3] or as an exercise in numerical analysis with an interesting physical interpretation.

It will be shown that the terminal velocity Vterm is different from the minimum one Vmin. Naively, one would have expected that the minimum velocity is the terminal one, but it comes out that the projectile speed diminishes until it reaches a minimum and then it starts increasing. The terminal velocity is finally reached from below, i.e. from lower velocities. This means that the velocity is not a monotonous function but it shows a minimum different from the initial and terminal velocities. This behavior is the same for any value of the power n in the drag force.

If the initial velocity V0 in not horizontal, the problem may be studied numerically. There is an angle range where the mentioned behavior is found. Out of that range, the velocity becomes a monotonous function and the minimum speed is the initial one (if V0 < Vterm) or the terminal one (if V0 > Vterm)

Let us assume that a projectile of mass m is launched horizontally with an initial velocity V0 from the top of a cliff -see Figure 1.


It is subject to the action of a gravitational field g and of a drag force Fd which is proportional to a power n of the velocity. If q is the angle between the projectile velocity and the horizontal, the equations of motion are:

where Vx = V cos q, Vy = Vsinq and Fd = kVn. The proportionality constant in the drag force can be expressed in terms of the terminal velocity Vterm. In the equilibrium condition, mg = k and k = mg/. It is convenient to use dimensionless quantities; let us introduce:

In this way, we may write the motion eqs. (1) as follow:

We will rewrite the previous equations to suppress the time t from them: write vx and vy in terms of v and q; then multiply the equations by cos q (or sin q) and add (or subtract) them conveniently. In this way, we get a new system of equations:

This new system allows us to eliminate the time t, and get a Bernoulli equation [3]:

Notice that the above equation is not well define for q = p/2. Our problem is to solve Eq. (4) which can be written as a linear differential equation. For doing this, a change of variable is introduced: y = v–n. Thus, Eq. (4) becomes a linear one:

Using standard techniques for solving linear differential equations, we finally get:

The integral that appears above has an exact solution for any integer value of n [4]. Therefore, Eq. (6) is the solution we have been looking for; it gives us the velocity evolution for the whole trajectory in a closed form. We may write the explicit solution for some values of n. For n = 1, we get:

In Figure 2 is shown a plot of v as a function of q. It is clear that the terminal velocity is reached from below regardless the value of v0. The minimum velocity may be easily found, and it is:


This is a counter-intuitive result. One would have expected that the minimum velocity was 1, i.e. the terminal velocity.

The same result is found for other powers of n. For n = 2, the velocity is given by:

and for n = 3, we get:

We have plotted Eq. (9) in Figure 3 for initial velocities closed to the terminal one. The effect is quite dramatic for those values. The projectile diminishes considerably its velocity before increasing, until the terminal velocity is reached. Let us try to understand this phenomenon in simple terms. The projectile motion can be decomposed in the horizontal and vertical components. In both directions there is a drag force that diminishes the speeds. However, in the vertical direction the gravitational field accelerates the projectile; thus the velocity increases. At the beginning, the effect of the drag is stronger than the effect of gravity if the initial speed is large enough. But once the velocity has diminished, gravity becomes dominant and accelerates the projectile until the terminal velocity is reached. There is an interplay between two opposite causes (drag and gravity) that produces the curious behavior shown in Figures. 2 and 3.


We may wonder about what happens if the initial speed is not horizontal, i.e. if q0¹ 0. In this case, the ordinary differential equations displayed in (4) has been solved numerically. For that aim, a computer program was written that solves the system using a Runge-Kutta method of order fourth [5]. We have got several values of the minimum speed vmin as a function of the initial angle q0. See Figure 4 and 5 for n = 2 and n = 3 respectively. As can be seen from those Figures, vmin ® 1 as q0 increases if v0 > 1, in other words, the minimum velocity is the terminal one for initial angles above a critical value (q0 ~ 50o). Also for v0 < 1 there is a critical angle ( q0 ~ 10o) above which the minimum velocity is the initial one. This means that the velocity is a monotonous function above those critical angles. Below the mentioned values, the velocity is non monotonous and the minimum velocity is lower than the initial one (if v0 < 1) or lower than the terminal one (if v0 > 1).



In summary, we have studied the trajectory of a projectile launched with horizontal velocity, and found a closed expression for the velocity as a function of the trajectory angle with the horizontal. It is surprising that the minimum velocity does not coincide with the terminal one. For any initial velocity, the projectile diminishes its speed below the terminal one, and then it accelerates and reaches the terminal velocity. This behavior is qualitatively the same for any power-law drag force. If the initial speed is not horizontal, there is critical value for the initial angle below which the same behavior is observed. Above that critical angle, the velocity becomes a monotonous function.

Acknowledgments

One of the authors (ENM) thanks the National Research Council of Argentina (CONICET) for support.

Recebido em 23/02/04

Aceito em 09/04

  • [1] J.R. Marion, Classical Dynamics of Particles and Systems (Academic Press, 1965).
  • [2] M. Alonso and E.J. Finn, Fundamental University Physics (Addison-Wesley, 1968).
  • [3] S. Arfken, Mathematical Methods for Physicists (Academic Press, 1985).
  • [4] I.S. Gradshteyn and I.M. Ryzhik, Table of Integrals, Series and Products (Academic Press, 1980).
  • [5] W. Press, B.P. Flannery, S.A. Teukolsky and W.T. Vetterling, Numerical Recipes in Fortran (Cambridge Univ. Press, 1992).
  • Endereço para correspondência
    E. N. Miranda
    E-mail:
  • Publication Dates

    • Publication in this collection
      01 Oct 2004
    • Date of issue
      2004

    History

    • Received
      23 Feb 2004
    • Accepted
      Sept 2004
    Sociedade Brasileira de Física Caixa Postal 66328, 05389-970 São Paulo SP - Brazil - São Paulo - SP - Brazil
    E-mail: marcio@sbfisica.org.br