Parametric oscillator

From Wikipedia, the free encyclopedia

A parametric oscillator is a simple harmonic oscillator whose parameters (its resonant frequency ω and damping β) vary in time in a defined way

\frac{d^{2}x}{dt^{2}} + \beta(t) \frac{dx}{dt} + \omega^{2}(t) x = 0

This equation is linear in x(t). By assumption, the parameters ω2 and β depend only on time and do not depend on the state of the oscillator. In general, β(t) and/or ω2(t) are assumed to vary periodically with the same period T.

Remarkably, if the parameters vary at roughly twice the natural frequency of the oscillator (defined below), the oscillator phase-locks to the parametric variation and absorbs energy at a rate proportional to the energy it already has. Without a compensating energy-loss mechanism, the oscillation amplitude grows exponentially. (This phenomenon is called parametric excitation, parametric resonance or parametric pumping.) However, if the initial amplitude is zero, it will remain so; this distinguishes it from the non-parametric resonance of driven simple harmonic oscillators, in which the amplitude grows linearly in time regardless of the initial state.

A familiar experience of parametric oscillation is playing on a swing. By alternately raising and lowering their center of mass (changing their moment of inertia and, thus, the resonant frequency) at key points in the swing, children can quickly reach large amplitudes provided that they have some amplitude to start with (e.g., get a push). Doing so at rest, however, goes nowhere.

Contents

[edit] Transformation of the equation

We begin by making a change of variables

q(t) \ \stackrel{\mathrm{def}}{=}\   e^{D(t)} x(t)

where D(t) is a time integral of the damping

D(t) \ \stackrel{\mathrm{def}}{=}\   \frac{1}{2} \int^{t} d\tau \ \beta(\tau).

This change of variables eliminates the damping term

\frac{d^{2}q}{dt^{2}} + \Omega^{2}(t) q = 0

where the transformed frequency is defined

\Omega^{2}(t) = \omega^{2}(t) -  \frac{1}{2} \left( \frac{d\beta}{dt} \right) - \frac{1}{4} \beta^{2}.

In general, the variations in damping and frequency are relatively small perturbations

\beta(t) = \omega_{0} \left[b + g(t) \right]
\omega^{2}(t) = \omega_{0}^{2} \left[1 + h(t) \right]

where ω0 and bω0 are constants, namely, the time-averaged oscillator frequency and damping, respectively. The transformed frequency can be written in a similar way:

\Omega^{2}(t) = \omega_{n}^{2} \left[1 + f(t) \right],

where ωn is the natural frequency of the damped harmonic oscillator

\omega_{n}^{2} \ \stackrel{\mathrm{def}}{=}\   \omega_{0}^{2} \left( 1 - \frac{b^{2}}{4} \right)

and

\omega_{n}^{2} f(t) \ \stackrel{\mathrm{def}}{=}\   \omega_{0}^{2} h(t) -  \frac{1}{2\omega_{0}} \left( \frac{dg}{dt} \right) - \frac{b}{2} g(t) - \frac{1}{4} g^{2}(t).

Thus, our transformed equation can be written

\frac{d^{2}q}{dt^{2}} + \omega_{n}^{2} \left[1 + f(t) \right] q = 0.

Remarkably, the independent variations g(t) and h(t) in the oscillator damping and resonant frequency, respectively, can be combined into a single pumping function f(t). The converse conclusion is that any form of parametric excitation can be accomplished by varying either the resonant frequency or the damping, or both.

[edit] Solution of the transformed equation

Let us assume that f(t) is sinusoidal, specifically

f(t) = f0sin2ωpt

where the pumping frequency 2\omega_{p} \approx 2\omega_{n} but need not equal n exactly. The solution q(t) of our transformed equation may be written

q(t) = A(t)cosωpt + B(t)sinωpt

where we have factored out the rapidly varying components (cosωpt and sinωpt) to isolate the slowly varying amplitudes A(t) and B(t). This corresponds to Laplace's variation of parameters method.

Substituting this solution into the transformed equation and retaining only the terms first-order in f_{0} \ll 1 yields two coupled equations

2\omega_{p} \frac{dA}{dt} =  \left( \frac{f_{0}}{2} \right) \omega_{n}^{2} A -  \left( \omega_{p}^{2} - \omega_{n}^{2} \right) B
2\omega_{p} \frac{dB}{dt} =  -\left( \frac{f_{0}}{2} \right) \omega_{n}^{2} B +  \left( \omega_{p}^{2} - \omega_{n}^{2} \right) A

We may decouple and solve these equations by making another change of variables

A(t) \ \stackrel{\mathrm{def}}{=}\   r(t) \cos \theta(t)
B(t) \ \stackrel{\mathrm{def}}{=}\   r(t) \sin \theta(t)

which yields the equations

\frac{dr}{dt} = \left( \alpha_{\mathrm{max}} \cos 2\theta \right) r
\frac{d\theta}{dt} = -\alpha_{\mathrm{max}}  \left[\sin 2\theta - \sin 2\theta_{\mathrm{eq}} \right]

where we have defined for brevity

\alpha_{\mathrm{max}} \ \stackrel{\mathrm{def}}{=}\   \frac{f_{0} \omega_{n}^{2}}{4\omega_{p}}
\sin 2\theta_{\mathrm{eq}} \ \stackrel{\mathrm{def}}{=}\   \left( \frac{2}{f_{0}} \right) \epsilon

and the detuning

\epsilon \ \stackrel{\mathrm{def}}{=}\   \frac{\omega_{p}^{2} - \omega_{n}^{2}}{\omega_{n}^{2}}

The θ equation does not depend on r, and linearization near its equilibrium position θeq shows that θ decays exponentially to its equilibrium

\theta(t) = \theta_{\mathrm{eq}} +  \left( \theta_{0} - \theta_{\mathrm{eq}} \right) e^{-2\alpha t}

where the decay constant

\alpha \ \stackrel{\mathrm{def}}{=}\   \alpha_{\mathrm{max}} \cos 2\theta_{\mathrm{eq}}.

In other words, the parametric oscillator phase-locks to the pumping signal f(t).

Taking θ(t) = θeq (i.e., assuming that the phase has locked), the r equation becomes

\frac{dr}{dt} = \alpha r

whose solution is r(t) = r0eαt; the amplitude of the q(t) oscillation diverges exponentially. However, the corresponding amplitude R(t) of the untransformed variable x \ \stackrel{\mathrm{def}}{=}\   q e^{-D} need not diverge

R(t) = r(t)e D = r0eαtD

The amplitude R(t) diverges, decays or stays constant, depending on whether αt is greater than, less than, or equal to D, respectively.

The maximum growth rate of the amplitude occurs when ωp = ωn. At that frequency, the equilibrium phase θeq is zero, implying that cos2θeq = 1 and α = αmax. As ωp is varied from ωn, θeq moves away from zero and α < αmax, i.e., the amplitude grows more slowly. For sufficiently large deviations of ωp, the decay constant α can become purely imaginary since

\alpha = \alpha_{\mathrm{max}}  \sqrt{1- \left( \frac{2}{f_{0}} \right)^{2} \epsilon^{2}}

If the detuning ε exceeds f0 / 2, α becomes purely imaginary and q(t) varies sinusoidally. Using the definition of the detuning ε, the pumping frequency p must lie between 2\omega_{n} \sqrt{1 - \frac{f_{0}}{2}} and 2\omega_{n} \sqrt{1 + \frac{f_{0}}{2}}. Expanding the square roots in a binomial series shows that the spread in pumping frequencies that result in exponentially growing q is approximately ωnf0.

[edit] Intuitive understanding of parametric excitation

The above derivation may seem like a mathematical sleight-of-hand, so it may be helpful to give an intuitive derivation. The q equation may be written in the form

\frac{d^{2}q}{dt^{2}} + \omega_{n}^{2} q = -\omega_{n}^{2} f(t) q

which represents a simple harmonic oscillator (or, alternatively, a bandpass filter) being driven by a signal -\omega_{n}^{2} f(t) q that is proportional to its response q.

Assume that q(t) = Acosωpt already has an oscillation at frequency ωp and that the pumping f(t) = f0sin2ωpt has double the frequency and a small amplitude f_{0} \ll 1. Their product q(t)f(t) produces two driving signals, one at frequency ωp and the other at frequency p

f(t)q(t) = \frac{f_{0}}{2} A  \left( \sin \omega_{p} t + \sin 3\omega_{p} t \right)

Being off-resonance, the p signal is attentuated and can be neglected initially. By contrast, the ωp signal is on resonance, serves to amplify q and is proportional to the amplitude A. Hence, the amplitude of q grows exponentially unless it is initially zero.

Expressed in Fourier space, the multiplication f(t)q(t) is a convolution of their Fourier transforms \tilde{F}(\omega) and \tilde{Q}(\omega). The positive feedback arises because the + 2ωp component of f(t) converts the − ωp component of q(t) into a driving signal at + ωp, and vice versa (reverse the signs). This explains why the pumping frequency must be near n, twice the natural frequency of the oscillator. Pumping at a grossly different frequency would not couple (i.e., provide mutual positive feedback) between the − ωp and + ωp components of q(t).

[edit] Parametric amplifiers

The parametric oscillator equation can be extended by adding an external driving force E(t):

\frac{d^{2}x}{dt^{2}} + \beta(t) \frac{dx}{dt} + \omega^{2}(t) x = E(t).

We assume that the damping D is sufficiently strong that, in the absence of the driving force E, the amplitude of the parametric oscillations does not diverge, i.e., that αt < D. In this situation, the parametric pumping acts to lower the effective damping in the system. For illustration, let the damping be constant β(t) = ω0b and assume that the external driving force is at the mean resonant frequency ω0, i.e., E(t) = E0sinω0t. The equation becomes

\frac{d^{2}x}{dt^{2}} + b \omega_{0} \frac{dx}{dt} +  \omega_{0}^{2} \left[1 + h_{0} \sin 2\omega_{0} t \right] x =  E_{0} \sin \omega_{0} t

whose solution is roughly

x(t) = \frac{2E_{0}}{\omega_{0}^{2} \left( 2b - h_{0} \right)} \cos \omega_{0} t.

As h0 approaches the threshold 2b, the amplitude diverges. When h \geq 2b, the system enters parametric resonance and the amplitude begins to grow exponentially, even in the absence of a driving force E(t).

[edit] Other relevant mathematical results

If the parameters of any second-order linear differential equation are varied periodically, Floquet analysis shows that the solutions must vary either sinusoidally or exponentially.

The q equation above with periodically varying f(t) is an example of a Hill equation. If f(t) is a simple sinusoid, the equation is called a Mathieu equation.

[edit] History

Faraday (1831) was the first to notice oscillations of one frequency being excited by forces of double the frequency, in the crispations (ruffled surface waves) observed in a wine glass excited to "sing". Melde (1859) generated parametric oscillations in a string by employing a tuning fork to periodically vary the tension at twice the resonant frequency of the string. Parametric oscillation was first treated as a general phenomenon by Rayleigh (1883,1887), whose papers are still worth reading today.

Parametric oscillators have been developed as low-noise amplifiers, especially in the radio and microwave frequency range. Thermal noise is minimal, since a reactance (not a resistance) is varied. Another common use is frequency conversion, e.g., conversion from audio to radio frequencies. Parametric amplifiers (paramps) were first used in 1913-1915 for radio telephony from Berlin to Vienna and Moscow, and were predicted to have a useful future (Alexanderson, 1916). The early paramps varied inductances, but other methods have been developed since, e.g., the varactor diodes, klystron tubes, Josephson junctions and optical methods. A practical parametric oscillator needs two three connections. One for common, one to feed the pump, and one to retrieve the oscillation and maybe a fourth one for biasing. A parametric amplifier needs a fifth port for the seed. Since a varactor diode has only two connections it can only be a part of a LC network with four eigenvectors with nodes at the connections. This can be implemented as a transimpedance amplifier, a traveling wave amplifier or by means of a circulator.

[edit] References

  • Faraday M. (1831) "On a peculiar class of acoustical figures; and on certain forms assumed by a group of particles upon vibrating elastic surfaces", Phil. Trans. Roy. Soc. (London), 121, 299-318.
  • Melde F. (1859) "Über Erregung stehender Wellen eines fadenförmigen Körpers", Ann. Phys. Chem. (ser. 2), 109, 193-215.
  • Strutt JW (Lord Rayleigh). (1883) "On Maintained Vibrations", Phil. Mag., 15, 229-235.
  • Strutt JW (Lord Rayleigh). (1887) "On the maintenance of vibrations by forces of double frequency, and on the propagation of waves through a medium endowed with periodic structure", Phil. Mag., 24, 145-159.
  • Strutt JW (Lord Rayleigh). (1945) The Theory of Sound, 2nd. ed., Dover.
  • Kühn L. (1914) Elektrotech. Z., 35, 816-819.
  • Pungs L. DRGM Nr. 588 822 (24 October 1913); DRP Nr. 281440 (1913); Elektrotech. Z., 44, 78-81 (1923?); Proc. IRE, 49, 378 (1961).
  • Mumford WW. (1961) "Some Notes on the History of Parametric Transducers", Proc. IRE, 48, 848-853.

[edit] See also