Air friction modeled as a differential equation

Click For Summary
SUMMARY

The discussion focuses on modeling air friction using a differential equation for a falling object, specifically analyzing the equation ƩFy = mg - bv = ma. The user successfully derived the velocity 'v' in terms of time 't' as v = (mg/b)(1 - e^(-bt/m)), but encountered issues when attempting to express 'y' in terms of 't'. The importance of boundary conditions, particularly the assumption that v = 0 at t = 0, is emphasized as critical for the validity of the solution. The final expression for 'y' is noted to have dimensional inconsistencies, indicating a need for further refinement.

PREREQUISITES
  • Understanding of differential equations and their applications in physics
  • Familiarity with Newton's laws of motion
  • Knowledge of exponential functions and their properties
  • Basic grasp of boundary conditions in mathematical modeling
NEXT STEPS
  • Review the derivation of differential equations in physics, focusing on air resistance models
  • Study the implications of boundary conditions on solutions to differential equations
  • Learn about dimensional analysis to ensure consistency in physical equations
  • Explore numerical methods for solving differential equations when analytical solutions are complex
USEFUL FOR

Students and educators in physics, mathematicians interested in applied mathematics, and engineers working on motion dynamics and fluid resistance modeling.

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:

Similar threads

Replies
30
Views
2K
  • · Replies 27 ·
Replies
27
Views
2K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 8 ·
Replies
8
Views
1K
  • · Replies 6 ·
Replies
6
Views
1K
Replies
7
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
Replies
13
Views
2K
Replies
4
Views
3K