Solving Newton's Equations: Perturbative Solution

  • Thread starter Thread starter stunner5000pt
  • Start date Start date
stunner5000pt
Messages
1,443
Reaction score
4
LOng question but not that hard apparently...
Consider a body of mass m falling vertically from rest near the Earth's surface. if air resistance is ignored then Newton;s equation m \ddot{y} = -mg (1.1) has the solution \tilde{y} (x) = -\frac{1}{2} gt^2 + y_{0}(1.2) where Yo is teh initial position of the body. If air resistnace is taken into account Newton's equation must be modified to m \ddot{y} = -mg - \beta \dot{y} (1.3) where beta is a positive constant. Assuming taht air resistance can be considered to b a small pertubartion (i.e. \left|\beta \dot{y} \right| << \left| mg \right|) we can approximate (1.3) by the equation m \ddot{y} = -mg - \beta \dot{\tilde{y}}(1.4) where y(t) (tilde) is the unperturbed solution (with air resistance ignored) given by 1.2

Solve 1.4 and discuss how this perturbative solution which we call Yp(t) compares to Y(t) tilde and determine under waht conditions (i.e. for what values of t) the perturbative solution is valid.

now \dot{\tilde{y}} (x) = -gt
so am i simply going to substitute this into 1.4 and solve for y[t]?
doing that gets m \ddot{y} = -mg + \beta gt
and y(t) = \frac{gt^3 \beta}{6m} - \frac{gt^2}{2} + tC_{2} + C_{2}
is this the right way to go?
 
Physics news on Phys.org
stunner5000pt said:
m \ddot{y} = -mg - \beta \dot{y}
Should the beta-ydot bit be positive since it opposes the motion of the mass. air resistance doesn't make it go faster.

Solve 1.4 and discuss how this perturbative solution which we call Yp(t) compares to Y(t) tilde and determine under waht conditions (i.e. for what values of t) the perturbative solution is valid.

To solve that equation (1.4), I would solve for \dot{y} and then integrate to get a function for y. I don't think that the two models are compatible since one completely ignores air resistance, while the other takes it into account, so I might wrong, but that is what I would do.
 
so then what i have doen is correct? Since y tilde is totally diffferent from y itself. However y does depend on t , and x from my understanding. SO i cna solve for Yp(t) which i have given in my first post. But y tilde (t) ... would that be using this equation?

m \ddot{\tilde{y}} = -mg - \beta \dot{\tilde{y}}
so here i am exclusively solving for Y tilde(t) .
 
Assuming that \beta is small enough to be ignored, the "non-perturbative" solution, from my"= -mg is, of course, y= -(g/2)t2+ y0 and y'= -gt

Now, replace y' on the right of the equation by -gt:
my"= -mg- \beta(-gt)= -mg+ \beta gt
That should be easy to integrate.
 
Last edited by a moderator:
whats that latex you typed halls?

isnt that what i did in my first post though?
 
ok for the perturbed system I am solving
m \ddot{y} = -mg + \beta gt
integrate twice and i get
Y_{p} (t) = \frac{ygt^3}{6m} - \frac{gt^2}{2} + C_{1} t + C_{2}

the second part of this question asks under waht conditions this solution is valid. There seem to be no resitrictions on the value of t. So this is valid for all t?
ALso need to solve (1.3)
taht is
m \ddot{y} = -mg - \beta \dot{y}
integrate once

for the unperturbed system i get
m \ddot{y} = -mg - \beta \dot{y}
integrate once and i get
\dot{y} = -gt - \frac{\beta y}{m} + C
seem to be having trouble iwth this differential equation... any ideas on how to solve this?
 
Hi, I had an exam and I completely messed up a problem. Especially one part which was necessary for the rest of the problem. Basically, I have a wormhole metric: $$(ds)^2 = -(dt)^2 + (dr)^2 + (r^2 + b^2)( (d\theta)^2 + sin^2 \theta (d\phi)^2 )$$ Where ##b=1## with an orbit only in the equatorial plane. We also know from the question that the orbit must satisfy this relationship: $$\varepsilon = \frac{1}{2} (\frac{dr}{d\tau})^2 + V_{eff}(r)$$ Ultimately, I was tasked to find the initial...
The value of H equals ## 10^{3}## in natural units, According to : https://en.wikipedia.org/wiki/Natural_units, ## t \sim 10^{-21} sec = 10^{21} Hz ##, and since ## \text{GeV} \sim 10^{24} \text{Hz } ##, ## GeV \sim 10^{24} \times 10^{-21} = 10^3 ## in natural units. So is this conversion correct? Also in the above formula, can I convert H to that natural units , since it’s a constant, while keeping k in Hz ?
Back
Top