22.05 Reactor Physics - Part Twenty-Nine: Reactor Operation With Feedback Effects
22.05 Reactor Physics - Part Twenty-Nine: Reactor Operation With Feedback Effects
22.05 Reactor Physics - Part Twenty-Nine: Reactor Operation With Feedback Effects
1. Reference Material: See pp. 368 – 372 in “Light Water Reactor Control
Systems,” in Wiley Encyclopedia of Electrical and Electronics Engineering.
3. C
oolant Temperature
4. Void Coefficient
The primary coolant in both PWRs and BWRs performs several functions. These
are:
— Removal of heat.
— Moderation of neutrons so as to continue the neutron chain reaction.
— Shielding during maintenance.
If the coolant becomes less dense, it becomes less efficient as a moderator. One
way to decrease coolant density is to decrease its temperature. Another is to
create voids in it. This occurs in BWRs when steam bubbles form. BWRs can be
controlled by adjusting the recirculation flow which in turn controls the rate at
which voids are swept out of the core.
5. R
eactor Regulation
Negative coefficients of reactivity promote self-regulation of a reactor. Consider
a PWR and suppose that the demand on the turbine increases. The following
sequence then occurs:
— Turbine first stage steam pressure decreases.
— Steam flow from the steam generator increases. This causes steam
generator pressure to drop.
— The steam generator is a saturated system. So, its temperature also drops.
— The decrease in steam generator temperature causes a decrease in the cold
leg temperature of the primary coolant.
— Cooler primary coolant enters the reactor core. This denser coolant
increases neutron moderation.
— Reactor power increases and so does the temperature of the hot leg.
— Hotter primary coolant reaches the steam generator. Steam generator
temperature and pressure rise and the steam supply equals the demand.
The final result is that the reactor power has increased to equal the demand. Also
the difference between the hot and cold leg temperatures has increased but the
average of these two temperatures is unchanged.
1 0 7 0 p s ia , 5 5 3 ° F
S tea m S te a m
R e lie f
S p ra y S e p a r a t o rs HP
Turb ine
F e e d w a te r
S te a m 4 50 °F
P r e s s u r iz e r T ub e
B u n d le
LP
L iq u id Tu rb ine
HTRS P um p
H e a te r
F lo w G u id e
Hot S te a m C o ld L e g C on den sor
L eg G e n e ra to r
C o re (6 22 ° F )
C o ld L e g ( 5 6 5 ° F ) S a m p le S ec o n d a ry
R e a c to r Ve s s e l
Pum ps
2 2 5 0 p s ia L in e a n d
H ig h P re s s u r e C le a n u p
E m e rg e n c y C o o lin g S y ste m S e a W a te r
L o w P r e s su r e
E m e rg e n c y C o o lin g
8
P r im a r y F l o w 1. 6 ⋅1 0 lb m / h r
6
S te a m F lo w 9 ⋅10 lb m / h r
P r e s s u r i z e d W a te r R e a c to r ( P W R )
6. D
oppler Effect
As the temperature of the fuel increases, the U-238 resonances broaden
and capture more neutrons in reactions that do not lead to fission. Hence,
negative reactivity is generated. Time scale is seconds or less.
The Doppler Effect is an inherent safety feature that may prevent fuel
damage during an accident. Some research reactors (TRIGAs) are
designated to eject control rods thereby causing very rapid power
increases to hundreds of MWs. The Doppler Effect shuts these reactors
down so quickly that the total energy produced is very small. (Note:
TRIGAs use a special fuel that accentuates the U-238 resonance
absorption.)
neutron, the relative speed between the two goes up. Hence, that neutron will be
absorbed.
But for every neutron that is now newly absorbed, one that was previously at the
resonance energy is now too high and it only scatters off the nucleus. So, why is
there a net increase in absorption? The reason is that scattered neutron loses only
a slight amount of energy (small object bouncing off a large one) and on its next
collision, which will likely be with the fuel again, it will be absorbed. Thus, U
238 resonances broaden because (1) the U-238 nuclei vibrate more rapidly on heat
up and (2) the fuel is separate from the coolant so that successive interactions
occur in the fuel.
8. Xenon
Xenon is a fission product that absorbs neutrons. It is produced whenever
a reactor is at power and, after about 40 hours, reaches an equilibrium
value. Xenon peaks 11 hours after shutdown and then decays away over
several days.
Reactors must be designed with enough fuel to offset the effect of Xenon.
This increases the cost of the reactor and the complexity of the control
system.
At some points in a refueling cycle, there may not be enough excess fuel
to restart during peak Xenon. Such reactors are referred to as 'Xenon
precluded.'
9. Analysis of Xenon-135
Background:
— Xe-135 (and I-135 which decays to Xe-135) are at the peak of the
right side of the fission product yield curve.
— Xe-135 has a large resonance at 1.082 eV
— Xe-135’s cross-section is 2.7 x 106 barns for thermal neutrons. Its
cross-section for fast neutrons is not significant.
The following figure below shows the production/removal sequence.
Fission
γ Xe = 0.003
γ Te = 0.064 β β-
Fission 135 135 β 135 β 135 Ba135
Te I Xe Cs
λTe = 0.023 λI = 2.87 x10 −5 λXe = 2.09 x10 −5 λCs = 8.8 x10 −14
Xe136
There are two sources of xenon. The first is from a fission product decay
chain that begins with Te-135 which has a 6.4% yield. The second is
direct from fission with a 0.3% yield. So, the isotope decay chain is the
major source. There are also two sinks. The first, which is minor, is
decay to Cs-135. The second is a neutron absorption reaction that
produces Xe-136. The above diagram can be simplified:
— The decay of Te-135 is so fast that it can be assumed that I-135 is
formed directly from fission with a yield of 6.4%
— The decay of Cs-135 is so slow (2.6 million year half-life) that it
can be taken as stable.
The equations that describes iodine and xenon behavior are therefore:
∂I( t)
= γ I Σ f1Φ1 (t) − λI I(t),
∂t
∂X(t)
∂t
= γ Xe Σ f1Φ1 (t) + λI I(t) − σ a1
Xe
[
Φ1 (t) + λXe X(t) ]
Xe
Where Σ f1, Φ1 (t), and σa1 are the one-group macroscopic fission cross-
section, flux, and microscopic absorption cross-section for Xe135,
respectively.
∂I(t) ∂X(t)
If we set and to zero, then we can obtain the equilibrium
∂t ∂t
concentrations of iodine and xenon. They are:
γ I Σ f1Φ 1
I(∞) =
λ1
γΣ f1Φ 1
X(∞) =
σ a1 Φ 1 + λXe
Xe
Where γ = λ1 + γ Xe
Note that the iodine concentration is directly proportional to the flux level.
The xenon concentration is a more complex function of flux:
Xe
— For high power levels, σa1 Φ1 >> λXe and the equilibrium xenon
Xe
level approaches a constant which is given by γΣ f1 / σ a1 .
The time dependent behavior of both iodine and xenon is, for a constant
flux level given by:
I(t) =
γ I Σ f1Φ 1
λ I
[1 − exp(−λ t)]+ I(0) exp(−λ t),
1 1
Xe
X(t) = X(0) exp(−(σ a1 Φ 1 + λXe )t) +
γΣ f1Φ 1
σ a1 Φ 1 + λXe
Xe
[1 − exp(−(σ Xe
a1 Φ 1 )]
+ λXe t
−
γΣ f1Φ 1 − λI I(0)
I
λ −λ Xe Xe
− σ a1 Φ1
[exp(−(σ Xe
a1 Φ 1 + λXe )t) − exp(−λI t) ]
Where γ = γ I + γ Xe
Xenon Behavior:
Figures showing xenon behavior for various power profiles are shown
later in these notes.
C
P P
P P
C
the central rod will be worth very little and the peripheral ones will be
worth a lot. This situation has occurred in actual reactors and operators
have inadvertently withdrawn peripheral rods that placed the reactor on
too short a period thereby causing an automatic shutdown.
Xenon can cause problems with detector readings because changes in the
flux shape affect the number of neutrons that reach a given detector.
Xenon can also cause spatial flux oscillations if a core is large. This
happens when the xenon in one quadrant is out-of-balance with that in the
other. A localized xenon peak slowly moves azimuthally around the core.
— There are two sources of Xenon: direct from fission and indirect
from the decay of iodine. The latter is the major source.
— There are two sinks for Xenon: burnup and decay to Cesium
which does not act as a neutron poison. Burnup is the major sink.
Xenon) is below its equilibrium value because the reactor has been
off-line. So, Xenon decreases.
10. S
amarium
On reactor shutdown, the only sink is removed. The source remains. So,
Samarium peaks and remains at peak value until operation resumes.
Background
γ Nd = .0113 β- β
Fission Nd 149
Pm149 Sm149 (Stable)
2 hr 54 hr
Thus, there is only one means of production and one means of removal. It
is common practice to ignore the Nd-149 because its decay time is fast
relative to that of Pm-149. Thus,
∂P(r, t) ∞
= γ Nd ∫ Σ f (r, E, t)Φ(r, E, t)dE − λPm P(r, t),
∂t 0
∂S(r, t) ∞
= λPm P(r, t) − S(r, t) ∫ σSm
a (E)Φ(r, E, t)dE
∂t 0
At equilibrium, we have:
γ Nd Σ f1Φ1
P(∞) =
λPm
γ Nd Σ f1
S(∞) =
σ Sm
a1
10
γ Nd Σ f1Φ1
P(t) = Pm
(1 − exp(−λPm t)) + P(0) exp(−λPm t)
λ
and
−
γ Nd Σ f1Φ 1 − λPm P(0)
λPm
− σ Sm
[exp(− σ Sm
a1 Φ 1 t )− exp(−λ t)]
Pm
a1 Φ 1
Where P(0) and S(0) are the concentrations of Pm149 and Sm149 at time
zero.
On startup of a new reactor, Sm-149 builds into its equilibrium level over
a period of several weeks. Its reactivity worth is about -1 Beta. On
shutdown, Sm-149 peaks (Pm inventory remains a source; burnup is gone
as a sink) and remains peaked until operation resumes.
Fuel vendors often add equilibrium Sm to the fuel so that its effect is not
seen at the initial startup of a reactor.
11