Power absorption in a lightly damped oscillator approaching Lorentzian form

Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
2 replies · 3K views
mbigras
Messages
61
Reaction score
2

Homework Statement


For a lightly damped harmonic oscillator and driving frequencies close to the natural frequency [itex]\omega \approx \omega_{0}[/itex], show that the power absorbed is approximately proportional to
[tex] \frac{\gamma^{2}/4}{\left(\omega_{0}-\omega\right)^{2}+\gamma^{2}/4}[/tex]
where [itex]\gamma[/itex] is the damping constant. This is the so called Lorentzian function.



Homework Equations


[tex] \text{Average power absorbed} = P_{avg} = \frac{F_{0}^{2} \omega_{0}}{2k Q} \frac{1}{\left(\frac{\omega_{0}}{\omega}-\frac{\omega}{\omega_{0}}\right)^{2}+\frac{1}{Q^{2}}} \\<br /> <br /> \omega_{0} = \sqrt{\frac{k}{m}}\\<br /> <br /> m = \frac{b}{\gamma}\\<br /> \text{where $b$ is the damping constant and $m$ is the mass}\\<br /> <br /> \Delta \omega = \frac{\gamma}{2}[/tex]





The Attempt at a Solution


The course of action that I took goes like:
1.Find [itex]k[/itex] and [itex]Q[/itex] in terms of [itex]\omega_{0}[/itex] and [itex]\gamma[/itex].
2. Chug through and do some algebra (and it is here that its very possible that a mistake was made, but I'll put my result not all the steps).
3. Expand a function about [itex]w_{0}[/itex] and make approximations so that [itex]\Delta \omega[/itex] is small.
(4) See the above equation fall out. This is the stage that I'm stuck at.

[tex] k = b \frac{\omega_{0}^{2}}{\gamma}\\<br /> Q = \frac{\omega_{0}}{\gamma}\\<br /> 2 \Delta \omega = \gamma<br /> \\<br /> P_{avg} = \text{plug in and do lots of algebra...}\\<br /> P_{avg} = \frac{\frac{\omega^{2}\gamma^{2}}{(\omega+\omega_{0})^{2}}}{(\omega_{0}-\omega)^{2}+\frac{\omega^{2}\gamma^{2}}{(\omega+\omega_{0})^{2}}}[/tex]
Then taylor expanding [itex]f(\omega) = \frac{\omega^{2}}{(\omega+\omega_{0})^{2}}[/itex] about [itex]\omega_{0}[/itex]...

Am I on the right try here? I'd like that taylor expansion to equal [itex]\frac{1}{4}[/itex] because then the equation would match the one described in the question but I'm trying it by hand and with mathematica and I'm not seeing them match.
 
Physics news on Phys.org
mbigras said:

Homework Statement


For a lightly damped harmonic oscillator and driving frequencies close to the natural frequency [itex]\omega \approx \omega_{0}[/itex], show that the power absorbed is approximately proportional to
[tex] \frac{\gamma^{2}/4}{\left(\omega_{0}-\omega\right)^{2}+\gamma^{2}/4}[/tex]
where [itex]\gamma[/itex] is the damping constant. This is the so called Lorentzian function.

Homework Equations


[tex] \text{Average power absorbed} = P_{avg} = \frac{F_{0}^{2} \omega_{0}}{2k Q} \frac{1}{\left(\frac{\omega_{0}}{\omega}-\frac{\omega}{\omega_{0}}\right)^{2}+\frac{1}{Q^{2}}} \\<br /> <br /> \omega_{0} = \sqrt{\frac{k}{m}}\\<br /> <br /> m = \frac{b}{\gamma}\\<br /> \text{where $b$ is the damping constant and $m$ is the mass}\\<br /> <br /> \Delta \omega = \frac{\gamma}{2}[/tex]

The Attempt at a Solution


The course of action that I took goes like:
1.Find [itex]k[/itex] and [itex]Q[/itex] in terms of [itex]\omega_{0}[/itex] and [itex]\gamma[/itex].
2. Chug through and do some algebra (and it is here that its very possible that a mistake was made, but I'll put my result not all the steps).
3. Expand a function about [itex]w_{0}[/itex] and make approximations so that [itex]\Delta \omega[/itex] is small.
(4) See the above equation fall out. This is the stage that I'm stuck at.

[tex] k = b \frac{\omega_{0}^{2}}{\gamma}\\<br /> Q = \frac{\omega_{0}}{\gamma}\\<br /> 2 \Delta \omega = \gamma<br /> \\<br /> P_{avg} = \text{plug in and do lots of algebra...}\\<br /> P_{avg} = \frac{\frac{\omega^{2}\gamma^{2}}{(\omega+\omega_{0})^{2}}}{(\omega_{0}-\omega)^{2}+\frac{\omega^{2}\gamma^{2}}{(\omega+\omega_{0})^{2}}}[/tex]
Then taylor expanding [itex]f(\omega) = \frac{\omega^{2}}{(\omega+\omega_{0})^{2}}[/itex] about [itex]\omega_{0}[/itex]...

Am I on the right try here? I'd like that taylor expansion to equal [itex]\frac{1}{4}[/itex] because then the equation would match the one described in the question but I'm trying it by hand and with mathematica and I'm not seeing them match.
That was very good and basically done. The key point is that we Taylor expand around $\omega \approx \omega_0$, in which case
##\frac{\omega^{2}}{(\omega+\omega_{0})^{2}} \approx 1/4 ## and the answer is obtained.
 
mbigras said:
[tex]P_{avg} = \frac{\frac{\omega^{2}\gamma^{2}}{(\omega+\omega_{0})^{2}}}{(\omega_{0}-\omega)^{2}+\frac{\omega^{2}\gamma^{2}}{(\omega+\omega_{0})^{2}}}[/tex]
Multiply both the numerator and the denominator by (ω+ω0)2.
[tex]P_{avg} = \frac{\omega^{2}\gamma^{2}}{(\omega_{0}-\omega)^{2}(\omega+\omega_{0})^{2}+\omega^{2}\gamma^{2}}[/tex]
When ω is near to ω0 you can replace ω+ω0 by 2ω, and the formula for Pavg can be simplified by ω2.
[tex]P_{avg} = \frac{\gamma^{2}}{4(\omega_{0}-\omega)^2+\gamma^{2}}[/tex]