We often ignore drag in elementary physics problems because it complicates matters. Sometimes, drag is pretty easy to deal with, however. Consider the one dimensional problem of a ball being thrown straight upwards at an initial speed of v i = 1 0 m/s . The ball experiences a drag force of F = β v where v is the velocity of the ball and β is the drag coefficient, β = 0 . 1 kg/s . What is the maximum height in meters the ball reaches?
Details and assumptions
This section requires Javascript.
You are seeing this because something didn't load right. We suggest you, (a) try
refreshing the page, (b) enabling javascript if it is disabled on your browser and,
finally, (c)
loading the
non-javascript version of this page
. We're sorry about the hassle.
Perfect. :)
3.10766439
This problem can be solved with the use of elementary, ordinary linear differential equations. We first have (saying that upwards is specified by x > 0 ): F n e t = − F g − F drag Note that we have: F g = m g And F drag = β v = β x ˙ Giving us an ODE: m x ¨ = − m g − β x ˙ Dividing through by m and rearranging: x ¨ + m β x ˙ + g = 0 Note that β = m = 0 . 1 and that g = 9 . 8 (for our sign convention), which gives us our final form to solve: x ¨ + x ˙ + 9 . 8 = 0 Using the standard ODE techniques, we let x = e α t for which we solve: α 2 + α = 0 Gives us: α ( α + 1 ) = 0 And: α = − 1 , 0 A linear combination gives us: c 1 e − t + c 2 = x Clearly, the particular solution is then: x p = − 9 . 8 0 7 t Summing these two solutions: c 1 e − t + c 2 − 9 . 8 0 7 t = x ( t ) We are restricted by: x ( 0 ) = 0 And: x ˙ ( 0 ) = 1 0 Finding the constants under these initial conditions gives us: − 1 9 . 8 e − t − 9 . 8 t + 1 9 . 8 = x ( t ) It should be clear, then, that the maximum is at t = lo g ( 4 9 9 9 ) , which is, once plugged in, approximately: t ∈ [ 0 , ∞ ) max x ( t ) ≈ 3 . 1 0 7 7
From the Newton's second law of motion , we obtain ∑ F − m g − β v − m g − β v d y ∫ h i h f d y h f − h i = m a = m d t d v = m d t d v ⋅ d y d y = − m g + β v m v d v = − ∫ v i v f m g + β v m v d v = ∫ v f v i m g + β v m v d v . The RHS integral can be evaluated by substitution method. Let u = m g + β v → v = β u − m g and d v = β d u . Therefore h f − h i = β 2 m ∫ v = v f v i u u − m g d u = β 2 m ∫ v = v f v i ( 1 − u m g ) d u = β 2 m ( u − m g ln u ) ∣ v = v f v i = β 2 m ( m g + β v − m g ln ( m g + β v ) ) ∣ v = v f v i = β 2 m ( β ( v i − v f ) − m g ln ( m g + β v f m g + β v i ) ) . Assuming the ball is thrown from the ground, we have h i = 0 and v f = 0 at the maximum height. Thus h f − h i h f − 0 h f = β 2 m ( β ( v i − v f ) − m g ln ( m g + β v f m g + β v i ) ) = 0 . 1 2 0 . 1 ( 0 . 1 ( 1 0 − 0 ) − 0 . 1 ⋅ 9 . 8 ln ( 0 . 1 ⋅ 9 . 8 + 0 . 1 ⋅ 0 0 . 1 ⋅ 9 . 8 + 0 . 1 ⋅ 1 0 ) ) = 1 0 ( 1 − 0 . 9 8 ln ( 0 . 9 8 1 . 9 8 ) ) ≈ 3 . 1 0 7 6 6 m
# Q . E . D . #
From energy conservation one can show that maximum height will be reached when velocity is zero.
Now analyzing the ball's motion along vertical, we have
− m a = m g + β v where a is ball's acceleration along vertical.
→ − v × d h d v = m β v + g
→ β − m ( 1 − β v + m g m g ) d v = d h
Integrating within proper limits i.e. L.H.S. from 1 0 to 0 and R.H.S. from 0 to h m a x and substituting the given values, we have
h m a x = 1 0 − ( 9 . 8 ) × l n ( 9 . 8 1 0 + 1 ) = 3 . 1 0 8 m .
acceleration = v + g
v d v / d x = v + g
d x = v / ( v + g ) d v
Integrate. Take v limit from 10 to 0 and x limit from 0 to x.
Starting from Newton's Second Law, we have F = m z ¨ = − m g − β z ˙ . Note the negative signs in front of both force terms; the first arises because gravity points downwards and the second because the drag acts in the opposite direction of the velocity.
Substituting the variable v = z ˙ for the velocity of the ball, and rewriting z ¨ = d t d v = d z d v d t d z = v d z d v , we obtain the following:
m v d z d v = − m g − β v
This can easily by solved by separation of variables. We then integrate from 0 to h m a x for z and from v i to 0 for v to calculate the maximum height of the ball.
∫ v i 0 g + m β v v d v = − ∫ 0 h m a x d z
Doing some manipulations on the improper fraction on the left hand side: β m ∫ 0 v i ( 1 − g + m β v g ) d v = ∫ 0 h m a x d z
Integrating: h m a x = β m v i − β 2 m 2 g ln ∣ 1 + m g β v i ∣
Finally, we substitute the values of β , m , v i and g : h m a x = 1 0 − 9 . 8 ln ∣ 1 + 9 . 8 1 0 ∣ = 3 . 1 0 8 m
Newton’s second law m a = − β v − m g
d t d v = − m β v − g = − v − 9 . 8
d s d v ( d t d s ) = − ( v + 9 . 8 )
v d s d v = − ( 9 . 8 )
-\underset{10}{\overset{0}{\mathop \int }}\,\frac{v}{v+9.8}dv=\underset{0}{\overset{s}{\mathop \int }}\,ds
s ≈ 3 . 1 1
Just plug the value into this formula: β m v 0 − β 2 m 2 g ln ( g − β v 0 g )
where:
m = 0 . 1
g = − 9 . 8
v 0 = 1 0
β = 0 . 1
You'll get 3 . 1 0 7 6 6 4 3 9 . . . as the answer.
Now, the question would be: how did you derive that formula?
It might be simple, agreed, but this is no solution if you don't show it.
Can't edit my post again :-( the formula should be: β m v 0 − β 2 m 2 g ln ( m g − β v 0 m g )
Problem Loading...
Note Loading...
Set Loading...
The force on the ball are the one due to gravity and the other due to drag. As the ball is moving in upward direction, both these forces are acting in downward direction. Therefore,
m × d t d v = − m g − β v
d t d v = − ( m β v + g )
d x d v × d t d x = − ( m β v + g )
v × d x d v = − ( m β v + g )
β v + m g v × d v = − m d x
β d v − β ( β v + m g ) m g × d v = − m d x
Integrating both the sides, the LHS from 10 to 0, and the RHS from 0 to "h(max)"
β v (10 to 0) - β m g × β l n ( β v + m g ) (10 to 0) = − m x (o to h)
Substitution value of initial velocity, v = 10 m/s,
β 1 0 − β 2 m g × l n ( 1 0 β + m g ) − l n ( m g ) = m h
Substituting other values,
h = 1 0 − ( 9 . 8 ) × l n ( 0 . 9 8 1 + 0 . 9 8 )
h = 3 . 1 0 8