Sunday, 11 October 2015

Solutions of damped oscillator differential equation


I am reading about damped harmonic motion in my Physics book (Gerthsen Physik) and there are two things that irritate me:




It says that Stokes friction would be $$F = -m \gamma \dot{x}.$$ To me, it does not make any sense why it depends on the mass of the object. It means that the deceleration does not depend on the mass. Or is it just since it is proportional to the volume and the volume usually depends on the mass.



Then the differential equation is set up, I fully understand it, except the confusion about the friction.


$$ m \ddot{x} + \gamma \dot{x} + \omega_0^2 = 0 $$


With $$\delta = \frac{\gamma}{2}$$ and $$\omega = \sqrt{\omega_0^2 - \delta^2}$$ the solution to the $$x(t) = \exp(\lambda t)$$ approach are $$\lambda = -\delta \pm i \omega t$$.


One of the plus or minus variants are a valid solution, but the linear combination of both are all solutions. I understand that.


The book then says that the general solution is the combination of real and imaginary parts. Alternatively, one can add a phase angle and have one "compact" solution like this:


$$ x(t) = x_0 \exp(-\delta t) \exp(i ( \omega t + \phi ) ) $$


The $\exp(i \phi)$ certainly makes up a linear combination of real and imaginary parts, but how is it a combination like $$c_1 \exp(i\omega t) + c_2 \exp(-i\omega t) ?$$



Answer




For the Stokes friction, the book is just choosing to normalize the coefficient with the mass. Since the coefficient of linear friction is arbitrary, you can divide it by the mass if you like, and then it has units of decay time. The parameter $\gamma$ is traditionally taken to be the decay time.


The reason the combination of the two terms into a phase-shift works is the addition law for cosines:


$$ cos(A+B) = cos(A)cos(B) - sin(A)sin(B) $$


which is simply the real part of the more obvious formula:


$$ e^{i(A+B)} = e^{iA}e^{iB} $$


and the imaginary part of the above formula reproduces the law of addition of sines. When you split cos(\omega t + \phi) using the addition rule, you find the two separate oscillations from the general solution.


No comments:

Post a Comment

Understanding Stagnation point in pitot fluid

What is stagnation point in fluid mechanics. At the open end of the pitot tube the velocity of the fluid becomes zero.But that should result...