Soving an eqn with exponential & linear parts

Click For Summary
SUMMARY

The discussion centers on solving the equation et*A + t*B + C = 0, where the goal is to find the variable t in the context of a particle's motion affected by gravity and air drag. Participants highlight the complexity of equations with both exponential and linear components, noting that analytical solutions are generally not feasible. Instead, they recommend using numerical methods or approximations, such as Taylor series, to tackle the problem. The conversation also touches on the differences between linear and quadratic drag forces, emphasizing that the latter can lead to cleaner analytical solutions.

PREREQUISITES
  • Understanding of differential equations, particularly first-order equations.
  • Familiarity with exponential functions and their properties.
  • Knowledge of classical mechanics, specifically motion under gravity and drag forces.
  • Experience with numerical methods for solving equations that cannot be solved analytically.
NEXT STEPS
  • Learn about numerical methods for solving transcendental equations, such as the Newton-Raphson method.
  • Study Taylor series approximations and their applications in physics problems.
  • Explore the differences between linear and quadratic drag forces in motion equations.
  • Investigate the use of computational tools like MATLAB or Python for simulating particle motion under various forces.
USEFUL FOR

Physics students, engineers, and anyone involved in modeling motion under the influence of forces, particularly those dealing with drag and gravitational effects.

Ikeness
Messages
3
Reaction score
0

Homework Statement


Known Constants
A, B and C

Unknown (we are solving for)
t

e is the exponential as in ln(e) = 1

Homework Equations


e^{t*A} + t*B + C = 0
(we want to solve for t)


The Attempt at a Solution


This is not for school or anything.
This is actually part of an attempt at the solution of another problem. Reaching this equation I had to stop. I do not remember seeing something like this in Calculus or Differential Equations, so I need help.

My overall purpose is to find the initial velocity necessary to launch a particle to a specified location. This particle is subject to gravity and air drag. If any of you happen to have the solution to this (or something simular) on hand, that information would be worlds better than just getting me past this one step.

I have my progress somewhat documented in this ugly text file.
If you care to read on, feel free. But it is not necessary.

I started by looking at what these people were doing, and simply added gravity to the initial problem.
Except that I am dealing with a drag force -kv instead of -k*sqrt(v)
https://www.physicsforums.com/showthread.php?t=190387


My work:

ma=-kv + gm
a = (-k/m)v + g
a - g = (-k/m)v
dv/dt - g = (-k/m)v

1/dt - g/dv = (-k/m)v/dv
1/dt = [(-k/m)v + g] /dv
dv / [(-k/m)v + g] = dt
dv / [(-k/m)v + g] = dt
dv (-k/m) / [(-k/m)v + g] = (-k/m) dt
dv / [v + (-m/k)g] = (-k/m) dt

<integrate>
ln|v + (-m/k)g| = (-k/m)t + C

EDIT
Just want to point out a serious mistake made here.
It does not effect my problem too seriously. I would fix it except that I am now trying blochwave's suggestion about using -kv^2 instead of -kv for air drag.
<rase e^(everything)>
v + (-m/k)g = e^(-tk/m) + D
where (v > -gm/k) (D = e^C)

(if t = 0)
v_0 - gm/k = 1 + D
v_0 = gm/k + 1 + D
v_0 - gm/k - 1 = D

dx/dt = gm/k + e^(-tk/m) + D
dx = [gm/k + e^(-tk/m) + D] dt
<integrate>
x = tgm/k + -k/m e^(-tk/m) + Dt + F

(if t = 0)
0 = -k/m + F
F = k/m

x = tgm/k + -k/m e^(-tk/m) + Dt + k/m

<let us solve for t in the x = eqn>

x = tgm/k + -k/m e^(-tk/m) + Dt + k/m
x = t[gm/k - D] - k/m e^(-tk/m) + k/m
x - k/m = t[gm/k - D] - k/m e^(-tk/m)
xk/m - 1 = t[gmm/kk - Dm/k] - e^(-tk/m)


And this is where I am. My purpose here was to plug in the solution of t into the equation involving velocity. Hopefully then I can have an equation based on just velocity (v) and displacement (x) without any t parameterization mess.
But I cannot solve for this t.
 
Last edited:
Physics news on Phys.org
Equations with exponential and linear parts cannot be solved analytically in general. At best you can make a second order taylor approximation and solve the resulting quadratic.
 
You usually write the resisting force as F=kmv, not kv

so the y direction equation of motion is (' is dot! Denoting time derivative) my''=-kmy'-mg, the resistive force always opposes the velocity, hence its negative, and the gravitational force is down

dv/(kv+g)=-dt

1/k*ln(kv+g)=-t+C

kv+g=e^(-kt+kc)

knowing initial y velocity we assume is Vyi at t=0

Vy=-g/k+[k*Vyi+g]*e^(-kt)/k

Vy=dy/dt, so you can then integrate again to find y(t)

The x equation is similar but mx''=-kmx' and that equation is a little easier

That's for you to compare to your own work(which I see at a glance is different, but perhaps not woefully so), it's straight out of my classical mechanics textbook. Also mentioned in my classical mechanics textbook is that attempting to solve for things like time of flight and stuff can't be done analytically because the equations become transcedental, and unless you know much perturbation theory or have powerful numerical solving tools, you're done, and I bet that's the issue you're having right now

An easier problem to solve is when the resistive force is proportion to v SQUARED. Much cleaner solution and can be done analytically, and, like a linear dependence on v, is realistic for certain circumstances
 
blochwave said:
You usually write the resisting force as F=kmv, not kv

...

That's for you to compare to your own work(which I see at a glance is different, but perhaps not woefully so), it's straight out of my classical mechanics textbook. Also mentioned in my classical mechanics textbook is that attempting to solve for things like time of flight and stuff can't be done analytically because the equations become transcedental, and unless you know much perturbation theory or have powerful numerical solving tools, you're done, and I bet that's the issue you're having right now
Yes. Thanks. I was carrying around an extra m. I also really screwed up on that exponential step as I have noted above. But asside from that, our problems steps are very simular.

blochwave said:
An easier problem to solve is when the resistive force is proportion to v SQUARED. Much cleaner solution and can be done analytically, and, like a linear dependence on v, is realistic for certain circumstances
I like easier problems. Please let me know if your textbook has anything to say about this system. Thanks.

Running through with v SQUARED resistive force I stopped after the first integration step to make sure that I am doing it right. Please let me know.
my'' = -kmy'^2 - mg
y'' = -ky'^2 - g

dv/dt = -kv^2 - g
dv/(-kv^2 - g) = dt
dv/(kv^2 - g) = -dt

dv/(v^2 - g/k) = -kdt

<integrate>
using table entry from http://en.wikipedia.org/wiki/List_of_integrals_of_rational_functions
The one for the integral of
1/(x^2 - a^2)

I GET

-1/sqrt(g/k) * arctanh(v/sqrt(g/k)) (for |v| < |sqrt(g/k)|)
-1/sqrt(g/k) * arccoth(v/sqrt(g/k)) (for |v| > |sqrt(g/k)|)

OR

1/2sqrt(g/k) ln[ (sqrt(g/k) - v) / (sqrt(g/k) + v) ] (for |v| < |sqrt(g/k)|)
1/2sqrt(g/k) ln[-(sqrt(g/k) - v) / (sqrt(g/k) + v) ] (for |v| > |sqrt(g/k)|)


I do not know which to use. Please give me a suggestion.

Thanks.
-Joel
 
Sorry you typed all that out but

dv/(-kv^2 - g) = dt
dv/(kv^2 - g) = -dt


Look hard, what went wrong there?
 
Use this link to help

http://farside.ph.utexas.edu/teaching/336k/lectures/node30.html

Note that they just have a constant and v^2, the reason being realistically something's drag(which is what this models)doesn't really depend on mass but on the air it's traveling through, cross section, stuff like that, and notice that the initial problem he solves is for an object FALLING, which explains why he has the negatives different from us
 

Similar threads

Replies
2
Views
2K
Replies
8
Views
1K
  • · Replies 4 ·
Replies
4
Views
1K
  • · Replies 2 ·
Replies
2
Views
906
  • · Replies 19 ·
Replies
19
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 22 ·
Replies
22
Views
2K
  • · Replies 9 ·
Replies
9
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 3 ·
Replies
3
Views
1K