Air friction modeled as a differential equation

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
1 reply · 2K views
DoubleWs
Messages
1
Reaction score
0
1. Find 'y' in terms of 't' for the equation of the net force on a falling object ƩFy = mg-bv = ma
So I'm pretty sure I found the velocity in terms of time correctly, and I'm interested in having someone check my answer for 'y' in terms of time.

Homework Equations


ƩFy = mg-bv = ma
a = [itex]\frac{dv}{dt}[/itex]
b = a constant (drag coefficient)

The Attempt at a Solution


mg - bv = ma

mg - bv = m[itex]\frac{dv}{dt}[/itex]

g - [itex]\frac{bv}{m}[/itex] = [itex]\frac{dv}{dt}[/itex]

(g - [itex]\frac{bv}{m}[/itex])dt = dv

∫dt = ∫[itex]\frac{dv}{g-(bv/m)}[/itex]

g - [itex]\frac{bv}{m}[/itex] = u

[itex]\frac{du}{dv}[/itex] = [itex]\frac{-b}{m}[/itex]

[itex]\frac{-m}{b}[/itex]du = dv

∫dt = -[itex]\frac{m}{b}[/itex]∫[itex]\frac{1}{u}[/itex]du

t + C = [itex]\frac{-m}{b}[/itex]ln(u)

t + C = [itex]\frac{-m}{b}[/itex]ln(g - [itex]\frac{bv}{m}[/itex])

[itex]\frac{-bt}{m}[/itex] + C = ln(g - [itex]\frac{bv}{m}[/itex])

e[itex]^{-bt/m}[/itex] * C = g - [itex]\frac{bv}{m}[/itex]

[itex]\frac{mg}{b}[/itex] - C * e[itex]^{-bt/m}[/itex] = v

[itex]\frac{mg}{b}[/itex](1 - e[itex]^{-bt/m}[/itex]) = v

Ok, so there's 'v' in terms of 't'. Now...

[itex]\frac{mg}{b}[/itex](1 - e[itex]^{-bt/m}[/itex]) = [itex]\frac{dy}{dt}[/itex]

∫[itex]\frac{mg}{b}[/itex](1 - e[itex]^{-bt/m}[/itex])dt = ∫dy

[itex]\frac{mg}{b}[/itex]∫(1 - e[itex]^{-bt/m}[/itex])dt = ∫dy

[itex]\frac{mg}{b}[/itex](t + [itex]\frac{m}{b}[/itex]e[itex]^{-bt/m}[/itex]) - C = y

[itex]\frac{mg}{b}[/itex](t + [itex]\frac{m}{b}[/itex]e[itex]^{-bt/m}[/itex]) - 1 = y

Just wondering if all of this is right, and if now, how can I fix it.
 
Last edited:
Physics news on Phys.org
Hello DoubleWs,

Welcome to Physics Forums!
DoubleWs said:
1. Find 'y' in terms of 't' for the equation of the net force on a falling object ƩFy = mg-bv = ma
So I'm pretty sure I found the velocity in terms of time correctly, and I'm interested in having someone check my answer for 'y' in terms of time.

Homework Equations


ƩFy = mg-bv = ma
a = [itex]\frac{dv}{dt}[/itex]
b = a constant (drag coefficient)

The Attempt at a Solution


mg - bv = ma

mg - bv = m[itex]\frac{dv}{dt}[/itex]

g - [itex]\frac{bv}{m}[/itex] = [itex]\frac{dv}{dt}[/itex]

(g - [itex]\frac{bv}{m}[/itex])dt = dv

∫dt = ∫[itex]\frac{dv}{g-(bv/m)}[/itex]

g - [itex]\frac{bv}{m}[/itex] = u

[itex]\frac{du}{dv}[/itex] = [itex]\frac{-b}{m}[/itex]

[itex]\frac{-m}{b}[/itex]du = dv

∫dt = -[itex]\frac{m}{b}[/itex]∫[itex]\frac{1}{u}[/itex]du

t + C = [itex]\frac{-m}{b}[/itex]ln(u)

t + C = [itex]\frac{-m}{b}[/itex]ln(g - [itex]\frac{bv}{m}[/itex])

[itex]\frac{-bt}{m}[/itex] + C = ln(g - [itex]\frac{bv}{m}[/itex])

e[itex]^{-bt/m}[/itex] * C = g - [itex]\frac{bv}{m}[/itex]

[itex]\frac{mg}{b}[/itex] - C * e[itex]^{-bt/m}[/itex] = v
Everything up to this point looks pretty good. At least it's the same result that I came up with. :approve:
[itex]\frac{mg}{b}[/itex](1 - e[itex]^{-bt/m}[/itex]) = v
But here you're making a big assumption. The above relationship is only true if there is a boundary condition (initial condition) that states v = 0 at time t = 0. There is no such boundary condition in the problem statement.

It might not be a bad assumption though, all I'm saying is the problem statement never mentions it. (There probably should be some sort of boundary conditions though, because if the initial velocity was in the opposite direction [heading up], none of the above would apply for that situation, because the differential equation would become ƩFy = mg+bv = ma until the object turned around.)

All your work below is contingent upon this assumption, so I'll treat everything below as though v = 0 at time t = 0.
Ok, so there's 'v' in terms of 't'. Now...

[itex]\frac{mg}{b}[/itex](1 - e[itex]^{-bt/m}[/itex]) = [itex]\frac{dy}{dt}[/itex]

∫[itex]\frac{mg}{b}[/itex](1 - e[itex]^{-bt/m}[/itex])dt = ∫dy

[itex]\frac{mg}{b}[/itex]∫(1 - e[itex]^{-bt/m}[/itex])dt = ∫dy

[itex]\frac{mg}{b}[/itex](t + [itex]\frac{m}{b}[/itex]e[itex]^{-bt/m}[/itex]) - C = y
Everything seems to check out so far. :approve:
[itex]\frac{mg}{b}[/itex](t + [itex]\frac{m}{b}[/itex]e[itex]^{-bt/m}[/itex]) - 1 = y
Something is a bit wrong with the above. I'm assuming here that you have assumed that y = 0 at time t = 0. But even with that assumption, something is still not right. It could be simple matter of missing parenthesis and/or parenthesis not in the right place. But as the equation stands now, it's not dimensionally correct (the y has units of length, and the 1 has no units at all); and y ≠ 0 at time t = 0, which is what I'm guessing you wanted to assume in the first place.
 
Last edited: