1937-Analysis of Series Capacitor Application Problems
1937-Analysis of Series Capacitor Application Problems
1937-Analysis of Series Capacitor Application Problems
I
N T H E past several years series capacitors for the com tive ratio 6Ï line resistance to line reactance increases, and
pensation of line drop in power circuits have found ^putting in a series capacitor to decrease the reactance ob
increasing use, 1 - 9 because improved and automatic' viously increases this ratio. Consequently, if the machines
voltage regulation can, in many cases, be obtained more are susceptible to hunting, the condition is further ag
economically by this method than by any other means. gravated by the addition of capacitors. The case of hunt
Most of the applications have been entirely successful ing of a single synchronous machine operating from an
in improving the system performance. However in a few infinite bus was investigated several years ago.11,12 In
instances unforseen difficulties of an unusual character section II of the present paper, further analytical methods
have been encountered. These difficulties were the oc are developed involving 2 or more machines which enable
casion of a systematic and detailed analytical study of the size and location of a capacitor to be so chosen as to
series capacitor performance. This study has been paral product the least tendency to hunt consistent with proper
leled, as far as possible with laboratory tests and correlated voltage regulation.
with field experience. It is believed that a satisfactory In case 3, under certain circuit conditions an induction
explanation and understanding of the abnormal system motor operating with series capacitors will be self-excited
behavior have now been obtained and that adequate and will generate other voltages of lower than normal fre
methods for predetermining the system performance have quency and thus produce objectionable voltage fluctu
been developed, although it is recognized that the investi ations. Analytical methods are presented in section III,
gation cannot be regarded as absolutely final and complete. which make it possible to determine the effectiveness of
The objects of this paper are to present a survey of some methods for eliminating this trouble.
of the troubles which may be encountered in the applica It has been found in this investigation that the addition
tion of series capacitors to power circuits, to present cri of a resistor in shunt with the series capacitor will eliminate
teria for determining the range and probability of abnormal the difficulties of cases 1 and 3, while the synchronous ma
operation, and to describe methods for minimizing the chine hunting of case 2 is minimized or not obtained by
difficulties. proper choice and location of the capacitor. The shunt
Three types of abnormality have been investigated : resistor is of such a high value as not to interfere with the
capacitor's effective operation as a line drop compensator
1. Distorted and excessively large transformer exciting currents, and the value of the power loss in the resistor necessary to
due to saturation. ensure complete freedom from these circuit troubles will
2. Hunting of synchronous machines. in general be less than 10 per cent of the resistance loss
3. Self-excitation of induction motors. in the line tô be compensated. The proper value of re
sistance to be used in any given case may be determined
In case 1, the application of voltage to an unloaded or from the analysis presented here. It is concluded, therefore,
very lightly loaded transformer through a series capacitor that the material offered in this paper, will enable prac
may result in abnormally large and distorted exciting cur tical and economically sound series capacitor applications
rents. These large currents are not simply transients but to be made with confidence that the operating results will
persist as a steady-state condition. This type of circuit, be in accordance with predictions.
that is, a saturable inductance in series with a capacitor,
has been found to have at least 2, and probably more,
different steady states, depending on how the voltage is Section I—Distorted and Excessively Large
applied. Transformer Exciting Currents
Solutions of the equations for this type of circuit have
been found by the differential analyzer10 at the Moore A. General
School of Electrical Engineering and are discussed in sec
tion I. In the usual case of application of voltage to an unloaded
In case 2, interconnected synchronous machines will transformer it is known that if the voltage is applied at or
near the zero point of the voltage wave a high inrush cur
A paper recommended for publication by the AIEE committee on power trans
mission and distribution. Manuscript submitted June 1, 1937; released for rent may result. This is because the transformer flux
publication June 28, 1937. tends to rise, in the first half cycle, to double its normal
J. W. BUTLER and C. CONCORDIA are in the engineering division, central station
department, General Electric Company, Schenectady, N. Y. value, with consequent excessive magnetizing current due
1. For all numbered references, see list at end of paper. to the very high degree of saturation. Stated in another
W M WM
at the zero point of the wave. During the first part of the voltages obtained
cycle the flux in the core builds up just as though the when transformer
capacitor were not present, since the current drawn is so is excited through
small that there is no appreciable voltage drop across the series capacitor
n
ISi
V.Q.Q.Q.QQ&-
-TJUirLT-
D. Conclusions
An unloaded transformer, operating at normal densities,
if excited through a series capacitor of reasonable reactance,
is susceptible to conditions of large distorted magnetizing
currents. These abnormal conditions can be restored to
normal by the addition of a suitable shunting resistor to
the capacitor, as indicated by figure 7, or by the locating
of a certain amount of load on the load side of the capaci
tor.
The resistor shunting the capacitor, introduces addi
tional losses. However, in general, this loss is less than
10 per cent of the normal line loss and this can usually be
justified by the improvement in voltage regulation. More
over, since the presence of the abnormal condition depends
^V C =0.Q33
on the initial transient, it may be eliminated by shorting
out the capacitor only while the starting switch is being
closed, or by having the resistor connected only during the
Figure 4. Effect of applied voltage on the apparent perio
energizing period. The shorting switch or temporary re
dicity of exciting current
sistor may readily be arranged so as to function automati
h = 0.049 1 / C = 0.214 (a) E - 1.0 cally whenever the abnormal condition tends to develop.
Rc = CO Magnetization curve 1 (b) E = 0.75
L = 0.197 a = 0 (c) E - 0.6
Section II—Hunting of Synchronous Machines
between the transformer and the capacitor, was also in
vestigated on the differential analyzer. The critical points A· General
obtained with this circuit agreed very well with the points
for the same resistance when it was shunting the capacitor, An electric power system consists in general of a group
indicating that it made little difference where the resistor of interconnected synchronous machines together with
was located. impedance and induction motor loads. Any such system
From figure 2 it may be seen that the 2 locations of the is subject at practically every instant of its existence to
resistor are equivalent to 2 locations of the applied voltage, disturbances caused by changes in load, switching, pul
the circuit remaining unchanged. Further, in the canon sating loads, pulsating driving torques, changes in excita-
*ψ- CURVE 1
, φ - CURVE 2
Figure 5. Effect of voltage on
excitation
Curve 1 E = 1.0. Normal excitation,
capacitor voltage decreases
Curve 2 E = 0.75. Abnormal excita V c -CURVET
tion, capacitor voltase builds up V c - CURVE 2
1 / C = 0.35 L = 0.394
Re = 6.3 Γι = 0.197
a =0 Masnetization curve 2
M
8 0.394 0.197 1 tion as a function
3. | CC=0 A. General
CURV E 2 | of Rc and 1/C.
o
y 1 /r Area below each It is well known that an induction machine will operate
1 / / curve is normal.
6 1 as a generator and that such operation depends on the
\ 1 1 Area above each
1
supplying of excitation from either a synchronous machine
\ \ // 1/
ί
curve is abnormal or a capacitor even though in case of the capacitor the
\ \ 1 f
operation has not previously been very well understood.
MAG. Now consider an induction motor supplied through a
/ ClJRVE
line containing a series capacitor. Then in addition to the
1
^ currents flowing due to normal operation of the motor, un
Figure 8 (below).
2 / / der certain conditions the motor may act as an induction
Magnetization
\ (y curves. Curve 2 is generator of current of lower than normal frequency. This
curve 1 with cur
low-frequency current is limited only by the impedance of
rent scale multi
the supply circuit at the low frequency, and may reach
0
0.2 0.4 0.6 plied by 3
relatively large values. These large low-frequency cur
CAPACITIVE REACTANCE j / C rents manifest themselves not only as current surges and
voltage swings but also as strong oscillations of the motor
I.6 M\\\ 1 ΠΤΤ7ΠΤΤ 7ΡΓ TTI illfff rotor caused by the large pulsating torques produced.
/ lifitn
W"ΊΊΊΙΙΙΙ 11 This phenomenon of self-excitation will not always take
'U M1
I.4
place and may always be entirely eliminated by sufficient
r* uL IP--· Tl ffl
4
jjlll
Ίmn
11 illll 11 line or shunt resistance, or by judiciously locating the
y
ma\
mmvC* capacitor. It may, moreover, be calculated with reason
§0.8 ^ i™ able certainty and accuracy by equations given in ap
Ί
I 0.6 ΊI lU-rrU
Iif / 1 if pendix C, if the circuit constants and operating conditions
ml /
IJ Jim
£o.4 are known.
K 11 Hill j11
0.2 Ml
A l fl
Um
Figures 10 to 13 show the regions in which self-excita
0 y
\M\\ ( 1 Illll 11 tion is likely to occur for the case of a purely induction
0.OOI 0.01 0.1 i.O K
P£« UNIT CURRENT-PEAK motor load. Figures 10 and 11 show the effect of series
ift(i
0.4 \ excitation.
7 0.6 To-
/> f
V the
S. B.
line charging
Crary, 21
M.
capacity of synchronous generators.
Takahashi, 22
and Y. H. Ku 23 have made
I7.I4 fl
■ \
cc
a study of this subject and have shown the tendency of
/ UN 5TA 3LE
II
8 0.6
/ ' /// \ line resistance to reduce the self-excitation. By means of
\ the theory given in reference 20, it may further be shown
0.4 1 'h
.D=I7 .4
\
that a resistance in shunt with the capacitor may also be
used to eliminate the self-excitation (see appendix D).
/ 1 \
Figure 18 shows, for the same generator as in figure 1 of
0.2 /
L1
reference 21, the effect of shunt resistance.
UN ST/û BLE
k
Since the effect of rotor circuits in addition to the main
1
1
for stable operation XURVEi \
02 \ / \ as a function of the
A k \ line capacitance for > | x'=0.4,To='7l
LÛLL
0.1
_L 0.2 0.3
Γ
0.4 0.5 0.6 0.7 0.8 0.9 1.0
various systems
x = 3 r = 0
X'=0.I6JO=I7I
SERIES LINE CAPACITIVE REACTANCE-Xc χ'=0.4,Το=Ι7.Ι
ω = 0.975
Figure 13. Self-excitation of induction motors with series ! I x'=0.4,To=m
0.8|
capacitors H CURVES X'=0.4,TO=I7.I
0.6
Showing effect of varying motor speed
-ω = 0.25 Β—ω = 0.5 C—ω = 0.707 D—ω = 1.0 0.4
X'=0.I6,TO=I7I
0.4 3.24 To = 1 7 1 . 4 r = 0
0.2
0.2 0.4 0.6 0.8 I.O
RATIO χ ς / χ '
would be to set up the equations of motion of the circuit
including both electrical and mechanical motions, and
then to find out whether the system as a whole was stable R — resistance in shunt with the capacitor
or unstable. This method brings out very clearly the basic k
identity of the phenomena of hunting and self-excitation
m, n = subscripts given to currents, voltages, impedances, and
but is, from an engineering point of view, at present im
angles in relating them to certain machines and circuits
practicable because of the excessive labor involved in the as used in section I I
calculations. It is believed that the results obtained by
the approximate studies are practically the same as would
be found by the unified method. Appendix A—Derivation of Equations
Set Up on the Differential Analyzer
Nomenclature
The single-phase circuit analyzed is represented in figure 2. T h e
Section I magnetization curves shown in figure 8 are the 2 iron characteristics
t h a t were investigated and discussed in the body of the paper.
E — applied voltage, maximum value The differential equations applying to this circuit are:
/ = time in radians
a — phase angle of applied voltage dii 1 Ç,.
E sin (t + a) = nil + L i2)dt + (la)
fi = line resistance dt
i\ = line current
ii = current through Rc HR> i<ù dt (2a)
Rc = resistance shunting the capacitor
L = line inductance Rearranging and integrating ( l a ) into a form applicable to the
C = capacitance differential analyzer, there is:
=S[Esi
ψ — flux linkage
Vc = voltage across capacitor Lix + ψ sin (t + ex.) — riii I (ii — it)dt «lit (3a)
C
Sections II and III The schematic hook-up of the analyzer, giving the solution t o
E = internal voltage equations 2a and 3a is shown in figure 3. The input voltage was
ea, eq — direct and quadrature axis terminal voltages actually not cranked in by an operator. Instead an auxiliary set-up
Ψα> Ψα — direct and quadrature axis flux linkages was made t h a t gave the solution to the differential equation
id, iq = direct and quadrature axis currents <Px
xd, xq = direct and quadrature axis synchronous reactances (4a)
dt2
(section II)
r — stator line resistance per phase which is a sinusoid, and it was impressed automatically into the
Ô = angle of voltage E from some reference point hook-up. An operator was needed to crank t h e current into t h e
0 = angular displacement of rotor system by keeping a pointer representing ψ, t h e linkage, on t h e
ω — ρθ — rotor speed linkage-current characteristic. Values of applied voltage, linkage,
xc — line series capacitive reactance capacitor voltage, and total current were plotted by the output
T0 — rotor time constant with open-circuited stator tables.
x = direct or quadrature axis synchronous reactance of induc The procedure used to determine the points on the résistance-
tion motor including line reactance capacitive reactance curves of figure 7 was as follows. The differen
x' — direct or quadrature axis transient reactance of induction tial analyzer was started at a point corresponding to application of
motor including line reactance the voltage at the zero point of the wave and t h e several dependent
a = xc/R variables of the circuit plotted automatically on an output table.
2-0.15
V V nid2
Xiidi
—nidZ
+#4%4
~X2iq2+Xzig3
-\~niq2
= E3 sin δι
+ 7-4% = E4COS(5i — δ4)
nidi +nidi —xdqi — £ 4 sin(5i—δ 4 )
0.10 (lb)
/
0.05 If the effect of the amortisseur windings is neglected it has been
shown 12 t h a t the load angle below which negative damping of
i J
0.2 0.4 0.6 0.fi I.O I.2 I.4 I.6 I.8 2.0
\ j i
2.2 2.4
generator number 1 may occur is accurately given by the condition
PER UNIT x c (when the frequency of the mechanical oscillation is small compared
to normal-voltage frequency),
As indicated by the equation above there are several variables which (2b)
can be used in the interpretation of the analyzer results. The
quantities which would ordinarily be measured in test of an actual Solving equations lb for idl, we find
circuit are t h e transformer voltage and current. However, it was
0 0 0 0 -1
found by experience t h a t the most significant quantities, as far as
the efficient determination of a criterion for proper circuit operation Γ 1 #2 0 0 n n 0 0
Es cos δι x2 -Xz 0 0 n —n 0
was concerned, were the capacitor voltage and transformer flux.
0 1 1 -1 0 0 0 0
Since, as shown in the oscillograms, the current only appears as an àidi (3b)
0 -n 0 0 xq x2 0 0
appreciable magnitude at widely separated intervals of time (e.g.,
Es sin δι n —n 0 0 — %2 xz 0
figures la and 4c) and, when it does appear, it is likely to bear no
obvious relation to its previous value, thus it is very difficult to use
Ei cos ( δι - δ 4 ) X2 0 XA 0 n 0 n
E4sin(5i - δ 4 ) n 0 n 0 — %2 0 — Xi
current as a continuous indication of the behavior of the circuit.
On the other hand, the transformer flux, which bears a direct func
where Δ is t h e determinant of the coefficients of the currents in
tional relation to the current as given by the saturation curve, is equations lb.
always of a reasonably large magnitude and of a nearly sinusoidal
Differentiating with respect to δι, there is
wave shape. I t appears, as shown in figures 4 to 6, as a sinusoidal
wave having a slowly varying displacement, and the rate of change 0 0 0 0 -1 1 I -il
and magnitude of this displacement were found to be useful factors 0 x2 0 0 n n 0 0
in determining as quickly as possible from a given analyzer run, —Ez sin δι
whether or not the circuit would approach its proper steady state.
Χ2 — Xi 0 0 n -n 0
d 0 1 1 -1 0 0 0 0
Δ
The capacitor voltage was found to be still more useful in this τ ~ idi = 0
respect; t h a t is, if one had to judge from the behavior of only one
ddi -n 0 0 x< x2 0 0
Ez cos δι n -n 0 0 ~~x2 Xz 0
quantity whether the circuit was to approach a normal or abnormal - £ 4 ( l - £ ) s i n ( o i - -δ 4 ) X2 0 XA 0 n 0 n\
condition, one would choose t h e capacitor voltage. This voltage Eiil-k)cos(ôi- -δ 4 ) n 0 TA 0 -Xi 0 —Xi\
determined by its magnitude the rate of change of flux wave dis (4b)
placement, b y its rate of discharge the effectiveness of the shunt
resistor, by its duration the continuance of the shifting of displace 5δ 4
ment of the flux wave. I n general, it was discovered t h a t if ever where k = — , i.e., k is the ratio of the motion of machine number
Oui
the capacitor voltage increased in magnitude from cycle to cycle
4 to t h a t of machine number 1. For sinusoidal oscillations k may
(it will be understood here t h a t by "cycle" is meant the apparent
be taken as the ratio of the amplitude of the machine oscillations.
cycle of the low frequency current pulsations and not that of the Equation 4δ, together with equation 2b, may be written
applied voltage), for example, if it ever went higher the second time
than it did the first, the system behavior was abnormal. On the
0= Δ -Ez(Mzi sin δι + M6i cos δι) —
other hand, if it continually decreased in magnitude from cycle to ;*Γ-
cycle the circuit operation was sure to be normal. This is illus E4[Mn sin (0! - δ4) + Mm cos (δχ - δ4)](1 - k)
trated by figure 6 which shows curves of capacitor voltage taken a t (5b.l)
" ■ ^
J
CIRCUIT OPENED OR CLOS ED
Ez
t a n δι = — (5b.2) ^"
E4
Mzi + ( 1 — k) — (Mn c o s h + Λίβι s i n δ4) 3 0.6
■Ë3
*2 0 f\ — r2 0
-r2 0 Xq — x2 0 0 0.04 0.08 0.I2 0.I6 0.20 0.24 0.28 0.32 0.36
Mzi = ^2 + ^3 —r% xz x2 + x 3 #8 (6b. 1) œ=WH
x2 xA 0 — r2 7-4
r Figure 18. Regions of self-excitation of synchronous
* TA 0 x2 — *4
machine with armature capacitance as affected by shunt
x2 0 n r2 0 resistance R
x2 4 - xz - x 3 —n r2 + r 3 —rs
M,i = —r2 0 *« X2 0 (6b.2) xd = 1.0 xQ = 0.6 xd' = 0.30 T0 = 1,000 r = 0
X2 XA 0 r2 r4
r
2 TA 0 — X2 — XA
X2 0 *Ί — ?-2 0
x2 #3 0 —r2 — r3 Z2\ri + rz) + xq(nx2 — Xzr2) + ^1(^2^3 + r2rz)
tan 6i = (8b.2)
M 7i —r2 0 *« —x2 0 (6b.3) Z22(xg + Xz) + xq(x2Xz -f r2rz) + ri(x3r2 — x2rz)
r2 n 0 x2 Xz
where
r2 - f U —rA x4 x2 - f x 4 — XA
plification. Xz ri + n rz
If machine number 4 is eliminated by letting xA and/or YA be -Mn — XA r\ -r4 (9b.2)
infinite, Mn, Mu become negligible compared to M3U M*i, and x
YA a — XA
equations 5b reduce t o :
-Xz fi rz
x2 — r\ r2 -Mn = rz Xq Xz (9b.3)
X2 + Xz rz r2 + r z —rA Xq + Xi — XA
—r2 ~xQ x2
tan δι = (8b.l) xz n rz
x2 r\ —r2 — M&1 -rz Xq Xz (9b.4)
— r2 xg —x2 —XA f\ + ^A — rA
r2 + r 3 Xz x2 + #3 ons 5b EA = 0, or if x 00 t h e r e is,
_, Mn
δι = —tan (10b)
Mzi
which is equation 7b if EA = 0, and is equation 8b if XA or r4 = œ.
If EA = 0 b u t XA and r4 are finite, equation 10b m a y of course be put
in the form Sb by paralleling Z 2 and Z 4 , i.e., in equation 8b replace
r2 and x2 by the corresponding real and imaginary parts of the
parallel values of Z 2 and Z 4 . This reasoning shows that since tan δ
in (86.2) equals the ratio of t h e real and imaginary components of
the transfer impedance from machine number 1 to the infinite bus
(number 3) therefore —Mu/Mzi, where Mu and Mzi are given by
(66.1) and (66.2), must also be the ratio of the real and imaginary
parts of the transfer impedance of machine number 1 to the infinite
bus with ZA finite. In fact, directly from equation 36 it may be
seen that since :
Figure 17. Starting and self-excited current of induction Δίάι = —Μ2ϊΕι -f MziEz c o s δι — Μ6ιΕζ s i n δι +
motor with series capacitance in primary circuit M7IEA COS (δι - 6A) - MSIEA s i n (fc - δ4) (lib)
— Μ21/Δ is the imaginary part of the driving point admittance t e + a)[ed — ρψα + rid + ωψϋ] + xcid = u[eQ — ρψϋ + riq - ωψά]
at machine number 1 (p + a)[eq — ρφα + riq - ωψα] + Xch = —<a[ed — ρψα + ηα+ωψ9]
Mn/ Δ is the imaginary part of the transfer admittance from 1 to 3 (le)
Mn/ Δ is the corresponding real part where, for an induction motor with completely symmetric rotor and
Mn/ A is the imaginary part of the transfer admittance from 1 to 4 stator,
Msi/ Δ is the corresponding real part Ψα = -x(p)id (2c)
cates that instability lies wholly in the range, V(x + 3 x ' ) V 2 + (4xx' - r'2)[(x x') 2 - 4r' 2 ] (12c)
where r' — r/ω.
A ' < Xc < ω2Χ (9c)
If r = 0, there is
so that if xc < ω2χ' the system is stable when in fact there is required
a certain finite series or shunt resistance to assure stability. Equa M = (*-_*2 (13c.l)
tion 8c is therefore not a good criterion when xc is small as in the \ω/ιηα Χ 2Λ/ΧΧ'
case of series capacitors for line reactance compensation except
when the motor is running at low speeds. which occurs a t
Xç X ~j* X
(13c.2)
EXACT CRITERION
Z2 " ~~2
To find the true boundary in the region xc < ω2χ' the criterion of Similarly if a = 0, there is
Routh 24 is used. This is carried out numerically b y assigning values
to all of the system constants and testing equations 6c or 7c for fr\ = (x - x')
(13c.3)
stability with various values of line or shunt resistance. The boun \w/max 2
daries found in this manner are shown in figures 10 and 12 (solid
lines). These figures are calculated for a motor of 0.16 per unit which occurs a t the same xc as in equation 13c.2. These last rela
transient reactance and an external reactance of 0.24. In figures tions determine equations,
10 and 12 there are also plotted boundaries corresponding to opera
tion at the same speed but with some external rotor resistance, the (x — x')2xc
for < (14c.l)
total rotor resistance being 10 times the rotor winding resistance.* 2(x + x')
The effects of series and shunt resistors have been determined sepa
and
rately in order to show as simply as possible the orders of magnitude
of the various quantities involved. However, since there is always (x — x')2xc
some line resistance present, the value of a required need not be as f or r = 0 (14C.2)
2(x + x')
large as that given by the curves of figure 12. Figure 15 shows an
example of how the a necessary for stability can be decreased with which are the equations of the loci of the maxima of the r versus xc
increasing line resistance, and has been computed from equation Qc and a versus xc curves as the speed ω is varied. Then, if the capaci
by Routh's 24 method. tive reactance xc in a given circuit is known equations 14c determine
the value of r or a needed to eliminate self-excitât ion at all speeds.
EFFECT OF ROTOR S P E E D However, as mentioned in the body of the paper, it will not in general
be necessary to do this.
Equations 6c or 7c may be used to determine exactly the regions
of unstable operation at any rotor speed. However, the effect of FREQUENCY AND T I M E CONSTANT
speed may be more clearly and easily seen from the approximate OF THE AMPLIFYING OSCILLATION
expression (8c). Since the rotor open-circuit time constant is
usually rather large, equation 8c may be further simplified by neg The frequencies and time constant of all these components of the
lecting terms containing T0, whence transient currents may of course be found directly by solving equa
tions 6c or 7c. However, it has been found that within the normal
ï[-+-fe)Q]
range of circuit constants only one of these currents is amplifying
and so likely to be observed as a large current pulsation. Moreover,
Equation 10c shows that as the speed ω is varied the scales of the (xc — ω2χ) -{-jœ(ax + r)
r versus xc or a versus xc diagram are simply changed so that Pi 2
(15c)
[(xc — ω χ' + ar)T0 -\- ax -\- r] -\-jœ [(αχ' -f- r)T0 + 2#]
xc is proportional toco 2
(lie) A further obvious simplification is obtained by neglecting the
r and a are proportional to ω
terms in the denominator which do not contain T0. If pi — pr + jq
* If rotor resistance is zero (To — °») there is no instability. Also if T0 then pr is the decrement factor and (ω + ç) the frequency of the
or ω = 0 there is obviously no instability. possibly amplifying oscillation. The frequency is ω + 2 rather than