Is physics a drag?

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 = 10 m/s v_i=10~\mbox{m/s} . The ball experiences a drag force of F = β v F=\beta v where v v is the velocity of the ball and β \beta is the drag coefficient, β = 0.1 kg/s \beta=0.1~\mbox{kg/s} . What is the maximum height in meters the ball reaches?

Details and assumptions

  • The acceleration of gravity is 9.8 m/s 2 -9.8~\mbox{m/s}^2 .
  • The mass of the ball is 0.1 kg 0.1~\mbox{kg} .


The answer is 3.108.

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.

8 solutions

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 v d t m \times \frac{dv}{dt} = m g β v -mg - βv

d v d t = ( β v m + g ) \frac{dv}{dt} = -(\frac{βv}{m} + g)

d v d x × d x d t = ( β v m + g ) \frac{dv}{dx} \times \frac{dx}{dt} = -(\frac{βv}{m} + g)

v × d v d x = ( β v m + g ) v \times \frac{dv}{dx} = -(\frac{βv}{m} + g)

v × d v β v + m g = d x m \frac{v \times dv}{βv + mg} = -\frac{dx}{m}

d v β m g × d v β ( β v + m g ) = d x m \frac{dv}{β} - \frac{mg \times dv}{β(βv + mg)} = -\frac{dx}{m}

Integrating both the sides, the LHS from 10 to 0, and the RHS from 0 to "h(max)"

v β \frac{v}{β} (10 to 0) - m g β × l n ( β v + m g ) β \frac{mg}{β} \times \frac{ln(βv + mg)}{β} (10 to 0) = x m -\frac{x}{m} (o to h)

Substitution value of initial velocity, v = 10 m/s,

10 β m g β 2 × l n ( 10 β + m g ) l n ( m g ) = h m \frac{10}{β} - \frac{mg}{β^{2}} \times {ln(10β + mg) - ln(mg)} = \frac{h}{m}

Substituting other values,

h = 10 ( 9.8 ) × l n ( 1 + 0.98 0.98 ) h = 10 - (9.8) \times ln(\frac{1 + 0.98}{0.98})

h = 3.108 h = \boxed{3.108}

Perfect. :)

A Brilliant Member - 7 years, 7 months ago

3.10766439

Aamir Faisal Ansari - 7 years, 7 months ago

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 x>0 ): F n e t = F g F drag F_{net} = - F_g - F_\text{drag} Note that we have: F g = m g F_g=mg And F drag = β v = β x ˙ F_\text{drag}=\beta v=\beta \dot{x} Giving us an ODE: m x ¨ = m g β x ˙ m\ddot{x}=-mg-\beta\dot{x} Dividing through by m m and rearranging: x ¨ + β m x ˙ + g = 0 \ddot{x}+\frac{\beta}{m}\dot{x}+g=0 Note that β = m = 0.1 \beta=m=0.1 and that g = 9.8 g=9.8 (for our sign convention), which gives us our final form to solve: x ¨ + x ˙ + 9.8 = 0 \ddot{x}+\dot{x}+9.8=0 Using the standard ODE techniques, we let x = e α t x=e^{\alpha t} for which we solve: α 2 + α = 0 \alpha^2+\alpha=0 Gives us: α ( α + 1 ) = 0 \alpha(\alpha+1)=0 And: α = 1 , 0 \alpha=-1, 0 A linear combination gives us: c 1 e t + c 2 = x c_1e^{-t}+c_2=x Clearly, the particular solution is then: x p = 9.807 t x_p=-9.807t Summing these two solutions: c 1 e t + c 2 9.807 t = x ( t ) c_1e^{-t}+c_2-9.807t=x(t) We are restricted by: x ( 0 ) = 0 x(0)=0 And: x ˙ ( 0 ) = 10 \dot{x}(0)=10 Finding the constants under these initial conditions gives us: 19.8 e t 9.8 t + 19.8 = x ( t ) -19.8e^{-t}-9.8t+19.8=x(t) It should be clear, then, that the maximum is at t = log ( 99 49 ) t=\log\left(\frac{99}{49}\right) , which is, once plugged in, approximately: max t [ 0 , ) x ( t ) 3.1077 \max_{t\in[0, \infty)} x(t)\approx 3.1077

Tunk-Fey Ariawan
Mar 3, 2014

From the Newton's second law of motion , we obtain F = m a m g β v = m d v d t m g β v = m d v d t d y d y d y = m v m g + β v d v h i h f d y = v i v f m v m g + β v d v h f h i = v f v i m v m g + β v d v . \begin{aligned} \sum F&=ma\\ -mg-\beta v&=m\frac{dv}{dt}\\ -mg-\beta v&=m\frac{dv}{dt}\cdot\frac{dy}{dy}\\ dy&=-\frac{mv}{mg+\beta v}dv\\ \int_{h_i}^{h_f} dy&=-\int_{v_i}^{v_f}\frac{mv}{mg+\beta v}dv\\ h_f-h_i&=\int_{v_f}^{v_i}\frac{mv}{mg+\beta v}dv. \end{aligned} The RHS integral can be evaluated by substitution method. Let u = m g + β v v = u m g β u=mg+\beta v\;\rightarrow\;v=\dfrac{u-mg}{\beta} and d v = d u β \;dv=\dfrac{du}{\beta} . Therefore h f h i = m β 2 v = v f v i u m g u d u = m β 2 v = v f v i ( 1 m g u ) d u = m β 2 ( u m g ln u ) v = v f v i = m β 2 ( m g + β v m g ln ( m g + β v ) ) v = v f v i = m β 2 ( β ( v i v f ) m g ln ( m g + β v i m g + β v f ) ) . \begin{aligned} h_f-h_i&=\frac{m}{\beta^2}\int_{v=v_f}^{v_i}\frac{u-mg}{u}\;du\\ &=\frac{m}{\beta^2}\int_{v=v_f}^{v_i}\left(1-\frac{mg}{u}\right)\;du\\ &=\frac{m}{\beta^2}\left.(u-mg\ln u)\right|_{v=v_f}^{v_i}\\ &=\frac{m}{\beta^2}\left.(mg+\beta v-mg\ln (mg+\beta v))\right|_{v=v_f}^{v_i}\\ &=\frac{m}{\beta^2}\left(\beta(v_i-v_f)-mg\ln\left(\frac{mg+\beta v_i}{mg+\beta v_f}\right)\right). \end{aligned} Assuming the ball is thrown from the ground, we have h i = 0 h_i=0 and v f = 0 \;v_f=0 at the maximum height. Thus h f h i = m β 2 ( β ( v i v f ) m g ln ( m g + β v i m g + β v f ) ) h f 0 = 0.1 0. 1 2 ( 0.1 ( 10 0 ) 0.1 9.8 ln ( 0.1 9.8 + 0.1 10 0.1 9.8 + 0.1 0 ) ) h f = 10 ( 1 0.98 ln ( 1.98 0.98 ) ) 3.10766 m \begin{aligned} h_f-h_i&=\frac{m}{\beta^2}\left(\beta(v_i-v_f)-mg\ln\left(\frac{mg+\beta v_i}{mg+\beta v_f}\right)\right)\\ h_f-0&=\frac{0.1}{0.1^2}\left(0.1(10-0)-0.1\cdot9.8\ln\left(\frac{0.1\cdot9.8+0.1\cdot10}{0.1\cdot9.8+0.1\cdot0}\right)\right)\\ h_f&=10\left(1-0.98\ln\left(\frac{1.98}{0.98}\right)\right)\\ &\approx\boxed{3.10766\text{ m}} \end{aligned}


# Q . E . D . # \text{\# }\mathbb{Q.E.D.}\text{ \#}

Nishant Sharma
Nov 3, 2013

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 -ma\,=mg+\beta\,v where a a is ball's acceleration along vertical.

\displaystyle v × d v d h = β v m + g \rightarrow\,-v\times\frac{dv}{dh}\,=\frac{\beta\,v}{m}+g

\displaystyle m β ( 1 m g β v + m g ) d v = d h \rightarrow\,\frac{-m}{\beta}(1-\frac{mg}{\beta\,v+mg}) dv\,=dh

Integrating within proper limits i.e. L.H.S. from 10 10 to 0 0 and R.H.S. from 0 0 to h m a x h_{max} and substituting the given values, we have

h m a x = 10 ( 9.8 ) × l n ( 10 9.8 + 1 ) = 3.108 m h_{max}\,=10-(9.8)\times\,ln(\frac{10}{9.8}+1)\,=\boxed{3.108\;m} .

Tarun Prasad Sahu
Jan 22, 2014

acceleration = v + g v+g

v d v / d x v dv/dx = v + g v+g

d x dx = v / ( v + g ) v/(v+g) d v dv

Integrate. Take v limit from 10 to 0 and x limit from 0 to x.

Yan Qi Huan
Nov 8, 2013

Starting from Newton's Second Law, we have F = m z ¨ = m g β z ˙ F=m \ddot z=-mg-\beta \dot 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 ˙ v=\dot z for the velocity of the ball, and rewriting z ¨ = d v d t = d v d z d z d t = v d v d z \ddot z=\frac{dv}{dt}=\frac{dv}{dz}\frac{dz}{dt}=v\frac{dv}{dz} , we obtain the following:

m v d v d z = m g β v mv\frac{dv}{dz}=-mg-\beta v

This can easily by solved by separation of variables. We then integrate from 0 0 to h m a x h_{max} for z z and from v i v_i to 0 0 for v v to calculate the maximum height of the ball.

v i 0 v g + β m v d v = 0 h m a x d z \int_{v_i}^0 \frac{v}{g+\frac{\beta}{m}v}dv=-\int_0^{h_{max}}dz

Doing some manipulations on the improper fraction on the left hand side: m β 0 v i ( 1 g g + β m v ) d v = 0 h m a x d z \frac{m}{\beta}\int_0^{v_i} (1-\frac{g}{g+\frac{\beta}{m}v})dv=\int_0^{h_{max}}dz

Integrating: h m a x = m v i β m 2 g β 2 ln 1 + β m g v i h_{max}=\frac{mv_i}{\beta}-\frac{m^2g}{\beta^2}\ln\lvert1+\frac{\beta}{mg}v_i\rvert

Finally, we substitute the values of β \beta , m m , v i v_i and g g : h m a x = 10 9.8 ln 1 + 10 9.8 = 3.108 m h_{max}=10-9.8\ln\lvert1+\frac{10}{9.8}\rvert=3.108m

Newton’s second law m a = β v m g ma=~-\beta v-mg

d v d t = β m v g = v 9.8 \frac{dv}{dt}=-\frac{\beta }{m}v-g=-v-9.8

d v d s ( d s d t ) = ( v + 9.8 ) \frac{dv}{ds}\left( \frac{ds}{dt} \right)=-(v+9.8)

v d v d s = ( 9.8 ) v\frac{dv}{ds}=-(9.8)

-\underset{10}{\overset{0}{\mathop \int }}\,\frac{v}{v+9.8}dv=\underset{0}{\overset{s}{\mathop \int }}\,ds

s 3.11 s\approx 3.11

Tjandra Gunawan
Nov 5, 2013

Just plug the value into this formula: m v 0 β m 2 g β 2 ln ( g g β v 0 ) \frac{mv_{0}}{\beta}-\frac{m^{2}g}{\beta^{2}}\ln(\frac{g}{g-\beta v_{0}})

where:

m = 0.1 m=0.1

g = 9.8 g=-9.8

v 0 = 10 v_{0}=10

β = 0.1 \beta=0.1

You'll get 3.10766439... \boxed{3.10766439...} 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.

Guillermo Angeris - 7 years, 7 months ago

Can't edit my post again :-( the formula should be: m v 0 β m 2 g β 2 ln ( m g m g β v 0 ) \frac{mv_{0}}{\beta}-\frac{m^{2}g}{\beta^{2}}\ln(\frac{mg}{mg-\beta v_{0}})

Tjandra Gunawan - 7 years, 7 months ago

0 pending reports

×

Problem Loading...

Note Loading...

Set Loading...