Resumos
The motion of a falling drop whose mass grows by accretion is studied with Newtonian methods to the point of finding the position as a function of time. The equation of motion applied is the equation of motion of continuum mechanics in its Eulerian or space formulation. We study three examples of laws of accretion: mass growing linearly with time, mass growing linearly with the surface of the drop and mass growing proportionally to the product of surface and velocity. These examples are sometimes left as exercises, without further discussion, asking only for v(z), or the final velocity. We also show that the solutions have the correct limit of a particle of constant mass in free fall and of such a particle with friction linear in the velocity.
variable mass systems; accretion
Se estudia, mediante métodos Newtonianos, el movimiento de una gota que cae y cuya masa crece por acreción; se encuentra detalladamente la posición como función del tiempo. La ecuación de movimiento aplicada es la ecuación de movimiento de la mecánica de medios contínuos en su forma Euleriana o espacial. Estudiamos tres ejemplos de leyes de acreción: masa incrementándose linealmente con el tiempo, masa incrementándose propocionalmente a la superficie de la gota e incremento de la masa proporcional al producto de la superficie por la velocidad. También mostramos que las soluciones tienen el límite correcto, el de una partícula, con masa constante, en caída libre y el de esa partícula con fricción lineal en la velocidad.
sistemas de masa variable; acreción
variable mass systems; accretion
sistemas de masa variable; acreción
ARTIGOS GERAIS
Motion of a falling drop with accretion using Newtonian methods
Estudio mediante métodos Newtonianos del movimiento de una gota que cae y cuya masa crece por acreción
G. HernandezI, 1 1 E-mail: gpe@correo.azc.uam.mx. ; G. del ValleI; I. CamposII; J.L. JimenezIII
IArea de Física Atómica y Molecular Aplicada, División de Ciencias Básicas e Ingeniería, Universidad Autónoma Metropolitana, Azcapotzalco, D.F., México
IIDepartamento de Física, Facultad de Ciencias, Universidad Nacional Autónoma de México, D.F., México
IIIDepartamento de Física, División de Ciencias Básicas e Ingeniería, Universidad Autónoma Metropolitana, Iztapalapa, D.F., México
ABSTRACT
The motion of a falling drop whose mass grows by accretion is studied with Newtonian methods to the point of finding the position as a function of time. The equation of motion applied is the equation of motion of continuum mechanics in its Eulerian or space formulation. We study three examples of laws of accretion: mass growing linearly with time, mass growing linearly with the surface of the drop and mass growing proportionally to the product of surface and velocity. These examples are sometimes left as exercises, without further discussion, asking only for v(z), or the final velocity. We also show that the solutions have the correct limit of a particle of constant mass in free fall and of such a particle with friction linear in the velocity.
Keywords: variable mass systems, accretion.
RESUMEN
Se estudia, mediante métodos Newtonianos, el movimiento de una gota que cae y cuya masa crece por acreción; se encuentra detalladamente la posición como función del tiempo. La ecuación de movimiento aplicada es la ecuación de movimiento de la mecánica de medios contínuos en su forma Euleriana o espacial. Estudiamos tres ejemplos de leyes de acreción: masa incrementándose linealmente con el tiempo, masa incrementándose propocionalmente a la superficie de la gota e incremento de la masa proporcional al producto de la superficie por la velocidad. También mostramos que las soluciones tienen el límite correcto, el de una partícula, con masa constante, en caída libre y el de esa partícula con fricción lineal en la velocidad.
Palabras-clave: sistemas de masa variable, acreción.
1. Introduction
The motion of systems with variable mass has conceptual and mathematical difficulties that make its treatment a challenge for teachers and students alike. The typical example is the rocket, discussed in many texts without beginning from an equation of motion, and rather applying conservation of momentum in a clever way. Other examples are the motion of a rope falling from a table, a conveyor on which sand is dropped, and a raindrop whose mass grows by accretion. In this work we solve in detail this last problem by Newtonian methods considering three specific laws of accretion.
The relevance of this problem in several fields of science is pointed out by Krane [1]. The present work will be useful for those interested in conceptual problems in physics and specifically graduates and beginning graduate students, as well as for teachers, interested in the conceptual problems that variable mass systems exhibit.
Among the conceptual difficulties that these problems present is the equation of motion to be applied. Sometimes the equation used is as if it were Newtons second law, assuming now that p= m(t)v, but we must recall that this law applies to a constant mass particle on which only external forces act. Thus Tiersten [2] shows that this equation holds only because other terms of the general equation are zero. Also Krane [1] points out in a note that this equation is a particular case of a more general equation that we discuss here. However, by their very nature, variable mass systems are composed of many "particles"and the system is modeled as a continuum, where the "particle"is a small part of it on which now body forces as well as surface forces act. Tiersten [2] and Krane [1] have pointed out that there must be a more general equation of motion for dealing with variable mass systems, and we propose that such equation of motion is the equation of motion of continuum mechanics.
We find that there are two expressions for the generalization of Newtons law applicable to a continuum [2, 3-5]. One is the material or Lagrangian form
The other is the spatial or Eulerian form
In these equations pis the mass density, bis the body force per unit mass, is the stress tensor,and vis the velocity with respect to our reference frame. In the material or Lagrangian description, the system is a given material particle. Therefore the system has constant mass and this description is not appropriate to deal with variable mass system. In the spatial or Eulerian description the system is a particular volume of a continuum. Thus mass can enter or leave this system (sometimes called "control volume") and therefore this description is appropriate to deal with variable mass systems. By the way, the usual continuity equation for mass conservation, is given in the Eulerian description. The stress tensor gives the force on the surface of a region of the continuum. This is the way Cauchy conceptualized the force which the rest of continuum exerts on a small part of it. The relation between both descriptions is given in the appendix .The first equation can be obtained from the second taking into account mass conservation. We propose that the equation of motion to be applied to variable mass systems is the Eulerian formulation since now the system is a particular volume, fixed or in motion with respect to our reference frame, in which mass may enter or leave carrying or not momentum. It is in this formulation that momentum flux must be considered.
We analyze three different laws of accretion: mass growing in proportion to time, mass growing in proportion to surface, equivalent to assuming that the radius of the drop grows linearly with time, and mass growing in proportion to the product of surface and velocity, equivalent to assuming that the radius of the drop grows proportionally to the distance travelled in falling. In all cases we find the correct limit of a constant mass particle freely falling.
2. Newtonian formulation
We solve the problem by Newtonian methods, which imply knowing all the forces acting on the system, and from these calculating the trajectory of the system in physical space. In the present case we find that this is the most difficult part of solving the problem, which explains why in texts it is asked usually to find only the velocity as function of height. We use as equation of motion the volume integral of the Eulerian expression, which after a volume integration takes the form
where Fis the body force, in our case gravity and friction, and Φ is the momentum flux given by
In the case that the considered volume is in motion, like in the rocket, it is convenient to distinguish the relative velocity of the particles respect to the moving volume, so that v = vr +uwhere vris the velocity of the particles with respect to the volume, and uis the velocity of the volume. Then the momentum flux Φis expressed as [4]
The last integral in Eq. (3) is the force that the surrounding medium exerts on the mass enclosed by the surface.
In our case the momentum flux is zero, since the velocity of the mass sticking to the drop is zero. Also, the surface integral of the stress, corresponding to surface tension, is zero, because of the spherical symmetry. The surface integral of the pressure gives the buoyancy force, that we discard assuming a small drop. Then our equation of motion is as a particular case of Eq. (3)
which can be written in the form
Equation (6) seems the usual expression of Newtons second law, but it is not so. It is a particular case of Eq. (3) and has the same structure of the equation of motion of a particle subject to a friction linear in the velocity. In other cases, as in the rocket [6, 7] or in the rope falling of a table [8], the particular expression of this equation is different.
We need to specify the mass as function of time in order to have an equation of motion to solve. We study three usual cases, the mass being proportional to:
a) the time
b) the surface of the spherical drop, and
c) surface times the velocity
3. Accretion proportional to time
In this case we have
where m0is the initial mass and bis a constant. If we take gravity as the only body force our equation is
It is easy to take into account a friction force of the form
since then our equation is of the same form
A friction quadratic in vis more difficult to treat, since then we have a differential equation of Riccatis type. We consider only friction linear in the velocity. Now we proceed to solve the differential Eq. (9) with the initial condition
This equation is of type
Then the solution can be obtained with an integrating factor of the form e l P(t)dt. That is
where and cis a constant determined by initial conditions.
The solution satisfying the initial condition v(0) =0 is
where
This solution, Eq. (15), if correct must contain the case of a constant mass particle in free fall as a limit when k= 0and b →0.It must contain also the case of a constant mass particle falling with friction linear in the velocity. Some authors [1, 9, 10] consider the limit mo →0,which for the case without friction (α = 1)gives As we show below, this limit gives in other laws of accretion with some integer (see section 4 below and [1]). This result seems strange and may be that the limit mo →0is rather formal, with uncertain physical meaning.
Now, to show that Eq. (15) has the correct limits we proceed by cases, first without friction (k= 0,b →0)and then with friction (k ̸= 0,b →0).
The frictionless case is given by k≠ 0or α= 1, then
It seems that this solution diverges for b →0 (m(t) → mo),but writing
we see that with the binomial theorem we obtain the correct limit
The same result can be obtained by writing the solution Eq. (17) in the form
and applying lHopital rule to this indeterminate limit, since b →0implies m → mo.
Now for the case k≠ 0,b →0we must proceed carefully, since λ → ∞.
First we notice that m(t)-λ can be written as
With the change of variable
we obtain
But
thus
With this result it is easy to see from Eq. (15), that
because b (λ+ 1) = 2b+ k, and using Eq. (24)
which is the expected result.
A further integration of Eq. (15) gives the solution for the height of the center of mass of the drop, satisfying the initial condition y(0) = h(Fig. 1).
For the case without friction (λ= 1)this result seems to diverge, but we can see that it is not so.
First we rewrite Eq. (28) as
We define
where it is necessary to take the limit λ →0.Using lHopital rule again, we obtain
Using
we can write
On the other hand, if we define
it is evident that
Using Eqs. (32) and (34) we get
The particle of constant mass in free fall is obtained expanding that is
Then
The expected result (Fig. 2).
4. Accretion proportional to the surface of the drop
Now we make the assumption that
Then the equation of motion for this case is
The assumption, Eq. (38), is equivalent to the hypothesis that the radius of the drop grows linearly with time, since
and then
Therefore r (t) = ro + αtimplies that and we obtain the usual assumption that mass grows proportionally to the surface of the drop.
Our equation of motion is then
This differential equation can be solved by the same method used for the case of accretion proportional to time, and the solution with the initial condition v(0) = 0is
In the formal limit mo →0,or ro → 0,r= αtand then we obtain
as we mentioned previously, it is a strange result.
With the binomial theorem we can write this solution as
obtaining the free fall case, v= - gt, as α →0:An integration of Eq. (44), with the initial condition y(0) = h, gives (Fig. 3)
Again, an expansion with the binomial theorem to second order in tshows that we can obtain the free fall case as α →0.
5. Accretion proportional to the surface times the velocity
The assumption that
seems more natural if we notice that it is equivalent to assuming that
That is, from , we find that
Therefore our equation of motion is now the non linear equation
This equation can be transformed with the identity
Then the equation of motion is
Now we have a differential equation for v2of the same type we have solved before, and the solution, with the initial condition v(0) = 0, is
We can show with the binomial theorem that this solution reduces to the free fall case as β →0.
The solution can be put in terms of y, obtaining
Here h0is defined by
Then it is obvious that
and in order to obtain y(t)we have to solve the integral
From this equation it is easy to obtain the formal limit ro →0,considered by some authors [11]. In this case ho= hand Eq. (56) becomes
which after integration and some simplifications results
This result implies that
that as we mention before, is strange.
This integral is not immediate and it is convenient to define and in order to transform the integral to the form
We need another change of variable
so that we have now the integral
This integral is found in standard tables [12] and is given as
where 2F1 is the hypergeometric function. In terms of the original variables we have
that with a Taylor expansion around y= hresults in
Finally, inverting this series we find
which after squaring gives y(t)(Fig. 4) as
This time it is obvious that we get the free fall case as β →0.
6. Conclusions
We have solved in a general way the problem of the motion of a falling drop whose mass grows by accretion according to a specific law of accretion. We have considered three specific laws of accretion, and have solved the problem by Newtonian methods. Then we had to apply a generalization of Newtons second law, which we took as the Eulerian formulation of the equation of motion for a continuum. Specifying clearly the hypothesis that lead to the particular equation of motion, these examples were solved to the point of getting the path of the center of mass of the falling drop, which is the aim of the Newtonian method. As a check of the solution obtained, we obtained in all cases the correct limit of a constant mass particle in free fall and that of a particle falling with friction linear in the velocity.
Recebido em 10/9/2010; Aceito em 24/6/2011; Publicado em 10/10/2011
Appendix
The natural generalization of Newton's second law for a continuum is given by the Lagrangian description, since a small part of this continuum, a "particle" of constant mass, is followed through its motion under the action of external forces, the body force and the surface force given by a surface integral of the stress tensor. Then the equation of motion in this description is
where is the total derivative, also called material derivative, given by
Then
The term (pv · ∇) vcan be developed with the aid of the tensor identity
and the continuity equation for conservation of mass in a given volumen,
Thus
Then
Finally, the equation of motion for matter in a given volume is
which is the Eulerian description of motion, the last term representing the momentum flux.
- [1] K.S. Krane, Am. J. Phys. 49, 113 (1981).
- [2] M.S. Tiersten, Am. J. Phys. 37, 82 (1969).
- [3] L. Prandtl and O.G. Tietjens, Fundamentals of Hydro and Aerodynamics (Mc. Graw Hill, New York, 1934, Dover, New York, 1957), p. 233.
- [4] G.W. Housner and D.E. Hudson, Applied Mechanics Dynamics (Van Nostrand, New York, 1959), 2nd ed.
- [5] R. Aris, Vector, Tensor and the Basic Equations of Fluid Mechanics (Prentice Hall, Englewood Clif., 1962. Dover, New York, 1989), p. 102.
- [6] G. del Valle, I. Campos and J.L. Jimenez, Eur. J. Phys. 17, 253 (1996).
- [7] I. Campos, J.L. Jimenez and G. del Valle, Eur. J. Phys. 24, 469 (2003).
- [8] J.L. Jimenez, G. Hernández, I. Campos and G. del Valle, Eur. J. Phys. 26, 1127 (2005).
- [9] A. Sommerfeld, Mechanics (Academic Press, New York, 1964), p. 241 problem 1.6 and p. 256.
- [10] S.T. Thornton and J.B. Marion, Classical Dynamics of Particles and Sistems (Thomson Learning, Belmont, 2004), p. 385 problem 9.56, 5th ed.
- [11] G.R. Fowles and G.L Cassiday, Analytical Mechanics (Thomson Brooks/Cole, Belmont, 2005), p. 320 problem. 7.23, 7th ed.
- [12] I.S. Gradshteyn and I.M. Ryzhik, Tables of Integrals, Series and Products (Academic, New York, 1980), p. 297.
Datas de Publicação
-
Publicação nesta coleção
01 Nov 2011 -
Data do Fascículo
Set 2011