Harmonic oscillator (classical): Difference between revisions

From Citizendium
Jump to navigation Jump to search
imported>Paul Wormer
mNo edit summary
 
(9 intermediate revisions by 3 users not shown)
Line 1: Line 1:
{{subpages}}
{{subpages}}
[[Image:Springs.png|right|thumb|350px|Harmonic oscillator: mass ''m'' oscillates between −''x'' and ''x''. Equilibrium at ''x'' = 0.]]
{{Image|Springs.png|right|350px|Harmonic oscillator: mass ''m'' oscillates between &minus;''x'' and ''x''. Equilibrium at ''x'' <nowiki>= 0 </nowiki>.}}
In [[physics]], a '''harmonic oscillator''' appears frequently as a simple model for  many different types of phenomena. The simplest physical realization of a harmonic oscillator consists of a mass ''m'' on which a force acts that is linear in a displacement from equilibrium. By [[Hooke's law]]  a spring gives  to very good approximation a force that is linear  for small displacements and so figure 1 shows a simple realization of a harmonic oscillator.
In [[physics]], a '''harmonic oscillator''' appears frequently as a simple model for  many different types of phenomena. The simplest physical realization of a harmonic oscillator consists of a mass ''m'' on which a force acts that is linear in a displacement from equilibrium. By [[Hooke's law]]  a spring gives  a force that is linear  for small displacements and hence figure 1 shows a simple realization of a harmonic oscillator.
The uppermost mass ''m'' feels a force acting to the right equal to ''k x'', where ''k'' is Hooke's spring constant (a positive number). The mass in the middle is in equilibrium, and the mass at the bottom feels a force to the left equal to &minus;''k x''.  The force executed by the springs always acts so as to restore the mass back toward its equilibrium position. If the mass is pulled out of equilibrium to ''x'' = ''A'' and then let go, the mass will (in an ideal setup) oscillate forever between &plusmn;''A''. In practice, energy will be lost due to friction of the atoms in the spring, so that the oscillation will die out after some time.
The uppermost mass ''m'' feels a force acting to the right equal to ''k x'', where ''k'' is Hooke's spring constant (a positive number). The mass in the middle is in equilibrium, and the mass at the bottom feels a force to the left equal to &minus;''k x''.  The force executed by the springs always acts so as to restore the mass back toward its equilibrium position. If the mass is pulled out of equilibrium to ''x'' = ''A'' and then let go, the mass will (in an ideal setup) oscillate forever between &plusmn;''A''. In practice, energy will be lost due to air drag, so that the oscillation will die out after some time.


==Mathematical description==
==Mathematical description==
Line 19: Line 19:
-Am\omega^2 \cos(\omega t +\phi) + k A \cos(\omega t +\phi) = 0\;,
-Am\omega^2 \cos(\omega t +\phi) + k A \cos(\omega t +\phi) = 0\;,
</math>
</math>
which fixes <math>\omega = \sqrt{k/m}</math>. The quantity &omega;,  called ''angular frequency'', has dimension of frequency (1/time)  and is often written as  &omega; = 2&pi; &nu;, where &nu; is the ''frequency'' (number of cycles per second). The time ''T'' &equiv; 1/&nu; is called the ''period'' of the cycle (duration of one cycle),  
which fixes <math>\omega = \sqrt{k/m}</math>. The quantity &omega;,  called ''angular frequency'', has dimension of frequency (1/time)  and is often written as  &omega; = 2''&pi;'' ''&nu;'', where ''&nu''; is the ''frequency'' (number of oscillation cycles per second). The time ''T'' &equiv; 1/''&nu;'' is called the ''period'' of the cycle (duration of one cycle),  
<table border="2" align="center" rules="all" cellpadding="10%">
<table border="2" align="center" rules="all" cellpadding="10%">
<tr><td><math>T =  \frac{1}{\nu} = \frac{2\pi}{\omega} = 2\pi \sqrt{\frac{m}{k}}</math>
<tr><td><math>T =  \frac{1}{\nu} = \frac{2\pi}{\omega} = 2\pi \sqrt{\frac{m}{k}}</math>
Line 67: Line 67:
  V(x)-V(x_0) = \int_{x_0}^{x} \; ky\;\mathrm{d}y = \frac{1}{2}k x^2 - \frac{1}{2}k x_0^2,  
  V(x)-V(x_0) = \int_{x_0}^{x} \; ky\;\mathrm{d}y = \frac{1}{2}k x^2 - \frac{1}{2}k x_0^2,  
</math>
</math>
where the zero of the potential  is at the arbitrary point ''x''<sub>0</sub> (the lower limit of the integral). When we choose ''x''<sub>0</sub> = 0,  we find  ''V''(''x''<sub>0</sub>) = 0, so that ''V''(''x'') = (''k''/2)''x''<sup>2</sup>.  
where the zero of the potential  is at the arbitrary point ''x''<sub>0</sub> (the lower limit of the integral). Upon choosing ''x''<sub>0</sub> = 0,  we find  ''V''(''x''<sub>0</sub>) = 0, so that ''V''(''x'') = (''k''/2)''x''<sup>2</sup>.  


Introducing a minor generalization, we shift the origin  on the ''x''-axis from 0 to the arbitrary point ''x''<sub>0</sub>, i.e. we make the substitution ''x'' &rarr; ''x''&minus;''x''<sub>0</sub>.  Then the ''harmonic potential'' becomes the following quadratic function in the displacement ''x'' &minus; ''x''<sub>0</sub>,
Introducing a minor generalization, we shift the origin  on the ''x''-axis from 0 to the arbitrary point ''x''<sub>0</sub>, i.e. we make the substitution ''x'' &rarr; ''x''&minus;''x''<sub>0</sub>.  Then the ''harmonic potential'' becomes the following quadratic function in the displacement ''x'' &minus; ''x''<sub>0</sub>,
Line 78: Line 78:
V(x-x_0) =  V(x_0)\; +\; (x-x_0)\;\left[\frac{d V(x)}{dx}\right]_{x_0}\; +\; \frac{1}{2}\;(x-x_0)^2\; \left[\frac{d^2 V(x)}{dx^2}\right]_{x_0}\; + \cdots
V(x-x_0) =  V(x_0)\; +\; (x-x_0)\;\left[\frac{d V(x)}{dx}\right]_{x_0}\; +\; \frac{1}{2}\;(x-x_0)^2\; \left[\frac{d^2 V(x)}{dx^2}\right]_{x_0}\; + \cdots
</math>
</math>
We truncate the expansion after these three terms. We take ''V''(''x''<sub>0</sub>) as zero of potential  and  assume that ''x''<sub>0</sub> is a [[stationary point]] of the potential, i.e., that the potential at ''x''<sub>0</sub> has a minimum, a maximum, or a point of inflection (saddle point). This means that we assume that
The expansion is truncated after these three terms. Take ''V''(''x''<sub>0</sub>) as zero of potential  and  assume that ''x''<sub>0</sub> is a [[stationary point]] of the potential, i.e., that the potential at ''x''<sub>0</sub> has a minimum, a maximum, or a point of inflection (saddle point). This means that we assume that
:<math>
:<math>
\left[\frac{d V(x)}{dx}\right]_{x_0} = 0.
\left[\frac{d V(x)}{dx}\right]_{x_0} = 0.
Line 86: Line 86:
V(x-x_0) = \frac{1}{2}\; (x-x_0)^2\;  \left[\frac{d^2 V(x)}{dx^2}\right]_{x_0}.
V(x-x_0) = \frac{1}{2}\; (x-x_0)^2\;  \left[\frac{d^2 V(x)}{dx^2}\right]_{x_0}.
</math>
</math>
If we identify Hooke's (spring) constant thus:  
If Hooke's (spring) constant is identified thus:  
:<math>
:<math>
k = \left[\frac{d^2 V(x)}{dx^2}\right]_{x_0},
k = \left[\frac{d^2 V(x)}{dx^2}\right]_{x_0},
</math>
</math>
we have found the harmonic potential from the truncated Taylor series. Note that since ''k'' > 0, the potential has a minimum at ''x''<sub>0</sub> (the stationary point is a minimum).
we have found the harmonic potential from the truncated Taylor series. Note that since ''k'' > 0, the potential has a minimum at ''x''<sub>0</sub> (the stationary point is a minimum). It is exactly by his truncated Taylor expansion of the potential near a minimum that the harmonic oscillator enters many physical problems.


We went from the force to the potential by integrating over distance. Conversely, given a potential ''V'', it is known from [[classical mechanics]] that the corresponding force ''F'' is minus the derivative <i>V</i>',
The derivation went from the force to the potential by integrating over distance. Conversely, given a potential ''V'', it is known from [[classical mechanics]] that the corresponding force ''F'' is minus the derivative <i>V</i>',
:<math>
:<math>
F = -V' = - (x-x_0)\;\left[\frac{dV(x)}{dx}\right]_{x_0}  = - (x-x_0)\; k.
F = -V' = - (x-x_0)\;\left[\frac{dV(x)}{dx}\right]_{x_0}  = - (x-x_0)\; k.
</math>   
</math>   
So, one could say that we proved [[Hooke's law]] (force linear in displacement) by use of a truncated Taylor expansion.
So, one could say that [[Hooke's law]] (force linear in displacement) is proved by use of a truncated Taylor expansion.
 
==Energy==
==Energy==
[[Image:Harmonic field.png|right|thumb|350px|''E''-''x'' diagram. Vertical axis energy; horizontal axis displacement ''x''. The potential energy ''V''(''x'') = &frac12; ''kx''<sup>2</sup> is shown in red.]]
{{Image|Harmonic field.png|right|350px|''E''-''x'' diagram. Vertical axis energy; horizontal axis displacement ''x''. The potential energy ''V''(''x'')<nowiki> =</nowiki> &frac12; ''kx''<sup>2</sup> is shown in red.}}
The oscillating mass has kinetic plus potential [[energy]]:
The oscillating mass has kinetic plus potential [[energy]]:
:<math>
:<math>
E = \frac{1}{2} m \dot{x}(t)^2 + \frac{1}{2} k x(t)^2
E = \frac{1}{2} m \dot{x}^2 + \frac{1}{2} k x^2, \quad\hbox{with}\quad x=x(t).
</math>
</math>
Above we saw that
Above it was shown that
:<math>
:<math>
x(t) = A\cos\omega t\quad \hbox{and}\quad \dot{x}(t) = -A\omega \sin\omega t
x(t) = A\cos\omega t\quad \hbox{and}\quad \dot{x}(t) = -A\omega \sin\omega t
Line 115: Line 116:
\end{align}
\end{align}
</math>
</math>
where we used &omega;<sup>2</sup> = ''k''/''m'' and cos<sup>2</sup> + sin<sup>2</sup> = 1.  Hence, the total energy of the harmonic oscillator is constant (independent of time). It is quadratic in the amplitude ''A'' and proportional to the spring constant ''k''. The stiffer the string the more energy.
where it was used that &omega;<sup>2</sup> = ''k''/''m'' and cos<sup>2</sup> + sin<sup>2</sup> = 1.  Hence, the total energy of the harmonic oscillator is constant (independent of time). It is quadratic in the amplitude ''A'' and proportional to the spring constant ''k''. The stiffer the string the more energy.


During the oscillation the energy is converted from potential energy to kinetic energy and back.
During the oscillation the energy is converted from potential energy to kinetic energy and back.
Consider the turning point ''x'' = ''A''. We saw above that this point occurs for time ''t''<sub>max</sub> = 2&pi;n / &omega;, and since sin(2&pi;n) = 0, the velocity at this turning point is zero. Hence all energy &frac12;''k A''<sup>2</sup> is potential energy at this point. The mass travels to the left and picks up speed and therewith kinetic energy. At the equilibrium point ''x'' = 0 we have ''V''(''x'') = 0 and all energy &frac12;''k A''<sup>2</sup> is turned into kinetic energy. The mass loses speed and comes to a momentary standstill at the left turning point ''x'' = &minus;''A'', where again all energy is converted into potential energy.
Consider the turning point ''x'' = ''A''. Above it was shown that the velocity at the turning points is zero. Hence all energy &frac12;''k A''<sup>2</sup> is potential energy at these points. The mass travels to the left and picks up speed and therewith kinetic energy. At the equilibrium point ''x'' = 0 the potential energy ''V''(''x'') = 0, all energy &frac12;''k A''<sup>2</sup> is turned into kinetic energy. Moving beyond ''x'' = 0 the mass loses speed and comes to a momentary standstill at the left turning point ''x'' = &minus;''A'', where again all energy is converted into potential energy.[[Category:Suggestion Bot Tag]]

Latest revision as of 17:00, 25 August 2024

This article is developing and not approved.
Main Article
Discussion
Related Articles  [?]
Bibliography  [?]
External Links  [?]
Citable Version  [?]
 
This editable Main Article is under development and subject to a disclaimer.
CC Image
Harmonic oscillator: mass m oscillates between −x and x. Equilibrium at x = 0 .

In physics, a harmonic oscillator appears frequently as a simple model for many different types of phenomena. The simplest physical realization of a harmonic oscillator consists of a mass m on which a force acts that is linear in a displacement from equilibrium. By Hooke's law a spring gives a force that is linear for small displacements and hence figure 1 shows a simple realization of a harmonic oscillator. The uppermost mass m feels a force acting to the right equal to k x, where k is Hooke's spring constant (a positive number). The mass in the middle is in equilibrium, and the mass at the bottom feels a force to the left equal to −k x. The force executed by the springs always acts so as to restore the mass back toward its equilibrium position. If the mass is pulled out of equilibrium to x = A and then let go, the mass will (in an ideal setup) oscillate forever between ±A. In practice, energy will be lost due to air drag, so that the oscillation will die out after some time.

Mathematical description

The motion of the mass (the blue ball in figure 1) as a function of time t may be obtained from Newton's second law (F = ma). Bringing the force F and the acceleration a (times mass m) to the same side of Newton's equation, the harmonic oscillator equation becomes the following equation on the left:

We see that the second derivative of x is proportional to −x. There are several functions known to have this property (sine, cosine, exponent with imaginary argument), we choose fairly arbitrarily the cosine as a trial function. The oscillator equation is a second order ordinary differential equation and its general solution contains two integration constants. Further the unit of time is undetermined and hence we scale time by a real number ω. In total the trial function is

where the integration constants A and φ are to be determined from the initial conditions. Substitution into the harmonic oscillator equation gives

which fixes . The quantity ω, called angular frequency, has dimension of frequency (1/time) and is often written as ω = 2π ν, where &nu; is the frequency (number of oscillation cycles per second). The time T ≡ 1/ν is called the period of the cycle (duration of one cycle),

As stated, the amplitude A and the phase angle φ are integration constants. Often it is possible to choose the zero of time such that ωt + φ = 0, implying that φ is zero:

The amplitude A is determined by how far we move initially (at t = 0) the mass (in figure 1) away from equilibrium

.

The value of x is maximum (mass in point x = A in figure 1) for multiples of full periods,

It takes time T for a full cycle, i.e., for a round trip of the mass from A back to A.

From the form of the cosine function we know that x goes through zero (mass in equilibrium x = 0 in figure 1) at the following points in time:

The displacement is minimum (mass in point x = −A in figure 1, i.e., furthest to the left) for multiples of half period T/2,

The speed of the mass is the first derivative of position with respect to time,

,

which for the specific values of time becomes

Physically, after the mass is displaced from equilibrium over a distance A to the right, the restoring force −kx pushes the mass back toward its equilibrium position, causing it to accelerate to the left. When it reaches equilibrium (x = 0), there is no force acting on it at that instant, but it is moving at maximum speed −Aω (to the left), and by Newton's first (inertia) law it persists in its motion moving beyond equilibrium position. Before it is stopped it reaches position −A, and by this time there is a force acting on it again, pulling it back toward equilibrium (to the right in figure 1). The direction of the motion is turned around in the points x = ±A, which is why these are called turning points.

The whole process, known as simple harmonic motion, repeats itself endlessly with a period given by T = 2π √m/k. This relation means that the stiffer the springs (i.e., the larger k), the shorter the period (the faster the oscillations). Making the mass m greater has exactly the opposite effect, slowing the motion down.

One of the most important features of harmonic motion is the fact that the frequency (ω, ν) and duration (T) of a cycle depend only on the mass m and the stiffness k of the spring. That is, m and k only affect the argument of the cosine, but not the amplitude A. If the amplitude is increased, the mass moves faster, but the argument of the cosine remains the same. This invariance of frequency is used in timekeeping and in musical string instruments (the pitch does not depend on the amplitude of the string).

Harmonic potential

From classical mechanics it is known that the potential energy is minus the force integrated over distance:

where the zero of the potential is at the arbitrary point x0 (the lower limit of the integral). Upon choosing x0 = 0, we find V(x0) = 0, so that V(x) = (k/2)x2.

Introducing a minor generalization, we shift the origin on the x-axis from 0 to the arbitrary point x0, i.e. we make the substitution xxx0. Then the harmonic potential becomes the following quadratic function in the displacement xx0,

It is of interest to show the connection with the Taylor expansion of an arbitrary potential around x0

The expansion is truncated after these three terms. Take V(x0) as zero of potential and assume that x0 is a stationary point of the potential, i.e., that the potential at x0 has a minimum, a maximum, or a point of inflection (saddle point). This means that we assume that

Then the truncated Taylor series shrinks to

If Hooke's (spring) constant is identified thus:

we have found the harmonic potential from the truncated Taylor series. Note that since k > 0, the potential has a minimum at x0 (the stationary point is a minimum). It is exactly by his truncated Taylor expansion of the potential near a minimum that the harmonic oscillator enters many physical problems.

The derivation went from the force to the potential by integrating over distance. Conversely, given a potential V, it is known from classical mechanics that the corresponding force F is minus the derivative V',

So, one could say that Hooke's law (force linear in displacement) is proved by use of a truncated Taylor expansion.

Energy

CC Image
E-x diagram. Vertical axis energy; horizontal axis displacement x. The potential energy V(x) = ½ kx2 is shown in red.

The oscillating mass has kinetic plus potential energy:

Above it was shown that

so that

where it was used that ω2 = k/m and cos2 + sin2 = 1. Hence, the total energy of the harmonic oscillator is constant (independent of time). It is quadratic in the amplitude A and proportional to the spring constant k. The stiffer the string the more energy.

During the oscillation the energy is converted from potential energy to kinetic energy and back. Consider the turning point x = A. Above it was shown that the velocity at the turning points is zero. Hence all energy ½k A2 is potential energy at these points. The mass travels to the left and picks up speed and therewith kinetic energy. At the equilibrium point x = 0 the potential energy V(x) = 0, all energy ½k A2 is turned into kinetic energy. Moving beyond x = 0 the mass loses speed and comes to a momentary standstill at the left turning point x = −A, where again all energy is converted into potential energy.