# Neutron Noise Analysis with Flash-Fourier Algorithm at the IBR-2M Reactor

Authors: Mihai O. Dima; Yuri N. Pepelyshev; Lachin Tayibov
Journal: Advances in High Energy Physics
 (2013)

---

## Abstract

Neutron noise spectra in nuclear reactors are a convolution of multiple effects. For the IBR-2M pulsed reactor (JINR, Dubna), one part of these is represented by the reactivities induced by the two moving auxiliary reflectors and another part of these by other sources that are moderately stable. The study of neutron noise involves, foremostly, knowing its frequency spectral distribution, hence Fourier transforms of the noise. Traditional methods compute the Fourier transform of the autocorrelation function. We show in the present study that this is neither natural nor real-time adapted, for both the autocorrelation function and the Fourier transform are highly CPU intensive. We present flash algorithms for processing the Fourier-like transforms of the noise spectra.


---

## Body

## 1. Introduction

Neutronic processes in nuclear reactors have a probabilistic character due to the quantum mechanics of scattering and the stochastic of propagation in materials. Design is mostly performed on the equations of neutron flux transport (in energy and space) and associated effects (fission, thermal fluxes, etc.). Statistical deviations from average quantities give however the complete image of neutron physics in the reactor—the so-termed neutron noise, described by large using Markov-chain theories. The theories associated with the underlying stochasticity that produces the said fluctuations are actually century old [1, 2], stemming from population studies. It was thus shown that (Alphonse) de Candolle’s conjecture on the extinction of family names [3–5] leads to nonergodic behavior in the case of neutron chains. Such theories were picked up in nuclear physics by Feynman, de Hoffmann and others to describe neutron processes in fission [6, 7], leading to the Feynman α-formula:
(1)σZ2(t)〈Z(t)〉=1+ϵ(1-1-e-αtαt)
which shows the fluctuations (1) being over Poissonian—due to correlations of neutrons in the same chain—and (2) the neutrons produced in the same group statistically disappearing (exponentially) all at the same time. The modern theoretical approach is given by the Pal-Bell equation [8–14], as an applied case of the Chapman-Kolmogorov master equation to Markov-chain neutron processes.On top of said, neutron stochastic behavior is the modulation of the neutron flux by various reactivities: some due to 2-phase liquid flow (bubbling), or fuel embrittlement, mechanically induced reactivities, and so forth. Any addition to the spectrum may be detected and classified, issuing a specific warning. In this respect, neutron noise spectrum analysis is a very far reaching tool in nuclear safety.To analyze the noise, usually a Fourier transform of the autocorrelation function is performed. This is basically the same as the square of the Fourier transformed signal. In practice, one also needs a time window, given by ad hoc apodisation functions. This is rather arbitrary and nonnatural. We show a more natural approach passing the noise through a resonant analyzing instrument together with flash algorithms for processing the Fourier-like transforms of the noise spectra.


## 2. The IBR-2M Reactor

The IBR-2M [15–17] is a modification of the IBR-2 (2 MW nominal power) pulsed research fast reactor with PuO2 fuel elements. The reactor coolant is liquid sodium. The pulsed mode operation is enabled by a reactivity modulator consisting of two rotating parts—a main movable reflector (OPO, at 1500 rpm) and an auxiliary movable reflector (DPO, at 300 rpm)—shown in Figure 1, that create reactivity pulses. For nominal DPO rotation speed, the reactivity of every fifth pulse is positive; that is, the reactor becomes prompt neutron supercritical for ca. 0.400 ms (0.215 ms half-width), with a repetition frequency of 5 Hz (Figure 2). As a result, powerful power pulses of 5 Hz repetition frequency take place in the reactor. The energy range of the fast neutrons is up to 10 MeV, with an almost linear logarithmic distribution versus energy.Figure 1
Details of the IBR-2M reactor showing the active core and the two movable reflectors.Figure 2
Power versus time between two subsequent pulses. Data is normalized to the maximum pulse.The pulsed operation mode of the reactor is established when the prompt neutron supercriticality(δk-β) reaches the “equilibrium” value εm=εm0~1·10-3 (at 5 Hz) at which the reactor can be periodically pulsed. For supercriticality smaller than the “equilibrium” value, the amplitude (and consequently the energy) of each subsequent pulse is smaller than that of the previous one, which means that the reactor is attenuating. There are two main causes for pulse energy fluctuation in the reactor: the stochastic character of fission and neutron multiplication processes, and the fluctuation of external reactivity. Stochastic noise dominates power fluctuations at low neutron intensity for powers below 1 W. Pulse energy fluctuations at high power have adverse effects on the operation of the reactor: in the dynamics, startup and adjustment process, performance of the experimental equipment, and so forth. However, power fluctuations also have a positive aspect, being as a tool for reactor diagnosis. The most important characteristic of a pulsed reactor is the relative dispersion of pulse energy fluctuations:
(2)(σQQ-)2=Δst2(1+δ02)+δ02,
where Δst2=νΓ/2Sτ is the relative dispersion of stochastic fluctuations; δ02 is the relative dispersion of noise caused by external reactivity fluctuations; Γ and ν are the dispersion of the number of prompt neutrons, and the average number of neutrons in a single fission event, respectively; τ is the lifetime of prompt neutrons in the core; S=Ssp+Sdel is the neutron source intensity (that operates continuously during the pulse); Ssp is the intensity of spontaneous neutrons and those from the (α,n) reaction on oxygen 18O (part of the oxide fuel); Sdel=βeffνFandβeff are the intensity, and effective fraction of delayed neutrons respectively; F=L·W is the fission rate; L is the number of fission events in a second per watt; and W is the absolute reactor power. At high power (for IBR-2M above 100 W) the fluctuations of the external reactivity (component δ02) dominates. Basically, the mechanical vibrations of the two moving auxiliary reflectors cause this. All noise diagnostics of the reactor is based on research of this noise component, which naturally degrades in time. The power is dependent on the neutron flux; hence, monitoring reactor power indicates the size of the neutron flux, together with its fluctuations.


## 3. Power Noise Baseline and Fourier Transform

For a power noise analysis [18–21], we must first extract any reactor dynamics influences from the recorded waveforms (overcompensation by the AR power control rod, etc.). In this respect, we want to determine a baseline with respect to which reference the instantaneous power by subtraction obtains the power noise. Such method, generally for any power-generating systems, has been presented elsewhere [22]. The method applied here is a spline of third order over ±30 bins around each current point. Basically, the order of the spline and the number of bins were so determined such that the effective cut-off frequency is below the natural cut-off of the mechanical vibrations of the moving reflectors in our case 0.2 Hz. The result is presented in Figure 3.Figure 3
Reactor power during power up (red line) with fitted baseline (blue curve of spline orderk=3 over a time interval [-30,+30] bins).The concept of Fourier transform in engineering can be quite different from its mathematical analogue. In mathematics, a given signal (function of time) is transformed to another function of frequency, basically the coefficients of the signal in the canonically conjugate image.In engineering, a signal is long, and we need an “instantaneous” Fourier Transform at each given time, thus a function that provides a time—window-termed apodisation function. The question is, for discrete Fourier transforms (DFFT-[23]), how to perform said clipping, as the apodisation function introduces sidelobes to the central lobe describins the signal-best seen for the Parzen apodisation function, in Figure 4.Fourier spectrum for various apodising functions. The Parzen function shows best the sidelobes that occur, while the SRS function provides a relatively flattop spectrum with good side-rejection of the lobes.
(a)
(b)
(c)
(d)Traditionally, the Fourier transform of the power noise is computed from the autocorrelation function:(3)Pw~Vw2=12π∫-∞+∞eiwτrV(τ)dτ,
where rV(τ) is the autocorrelation function:
(4)rV(τ)=∫-∞+∞V*(t)V(t+τ)dt.
It is evident that although the Fourier transform is a single integral, the problem of calculating a double integral (high CPU cost), has not been eliminated, neither that of calculating the autocorrelation function.Beyond the above considerations, special tools are needed to compute—what is mathematically improper termed—Fourier Transform of noise signals. The basic image is that of sound noise, perceived as higher or lower pitch. In categorizing noise as such, the ear for instance (as one example of analysis instrument) does not work with orthogonal wave-vectors, or mathematically complete spaces, nor does it obtain coefficients for said wave-vectors. What such an analysis instrument does is to use a set of resonance phenomena, tuned on various frequencies. The instantaneous mechanical quantity producing measurable effect is (in the case of sound) the air velocity, or in electrical circuits, the current:(5)i=2w0R4Q2-1∫0tϵ(τ)ew0((τ-t)/2Q)=2w0R4Q2-1×sin(w0(τ-t)1-14Q2-ϕ)dτ,sinϕ=1-14Q2,cosϕ=-12Q,
where Q is the quality factor, ω0 the resonant pulsation, R the electrical circuit’s resistance, and ε(t) is the driving voltage. It can be seen that the “analysis instrument” performs “something like” a Fourier Transform with an asymmetrical Poisson apodising function, looking into the “past.”As shown above, this is a natural choice—with a definite model in sight—for an apodising function, containing what is natural asanalysis instrument.The remaining problem is that such an instrument is highly computationally the intensive.It can be recognized that the sine function can be written as sum of exponentials and that the mathematical construct to be computed is principally(6)∫0Tf(t)eiwtdt=∑f(nΔ)eiw(n-n0)Δ·Δ,
where Δ is the bin size and n the bin number over which summation occurs.There are two cases for how the summation padding is advantageous, depending on how{2π/ωΔ} compares to 1/2, where {⋯} and [⋯] are the fractional, respectively integer part:(i)
{2π/ωΔ}≥1/2: the sum at a given bin n0 is
(7)Sn0=∑n=n0n0+[2π/wΔ]f(nΔ)eiwΔ(n-n0+(1/2){2π/wΔ}).
The recurrence formula for obtaining the sum at the next point is obtained by shifting the exponentials with 1 bin, subtracting the lower end and adding the upper end:
(8)Sn0+1=e-iwΔSn0+f(n0+1+[2πwΔ])eiϕ-f(n0)e-iϕ,
where ϕ=ωΔ/2(2-{2π/ωΔ}). In doing so, most of the terms are known and can be computed in advance; the previous bin’s sum is also known, so very little needs to be computed for the current bin.(ii)
{2π/ωΔ}<1/2: the sum at a given bin n0 is
(9)Sn0=∑n=n0n0-1+[2π/wΔ]f(nΔ)eiwΔ(n-n0+1/2+(1/2){2π/wΔ}).
The recurrence formula for obtaining the sum at the next point is obtained by shifting the exponentials with 1 bin, subtracting the lower end and adding the upper end:
(10)Sn0+1=e-iwΔSn0+f(n0+[2πwΔ])eiϕ-f(n0)e-iϕ,
where ϕ=ωΔ/2(1-{2π/ωΔ}). In doing so, most of the terms are known and can be computed in advance; the previous bin’s sum is also known, so very little needs to be computed for the current bin.As mentioned in the introduction, there are beneficial aspects in detecting and subtracting the baseline. Apart from extracting from the low frequency end reactor dynamics features, also technical aspects are achieved, such as rejecting aliasing. In general, such unwanted phenomenon is due to the fact that some frequency is much higher than the Nyquist sampling frequency; in our case,fNy=1/2Δt=2.5Hz. Imagine the effect of 50 Hz ripple for instance, or that of certain reactor chirps. Active analog filters would come to mind as first resource since they have a somewhat higher Q-factor than passive filters. (Passive filters in the range [0, 2.5 Hz] are too bulky in our view of a DAQ setup.) Still, even active filters do not have that high Q-factors, and if they do, 50 Hz sources do have a certain bandwidth, rendering very high-Q filters useless. As for reactor chirps, they are of broad bandwidth and impossible to filter. All these phenomena create unwanted aliasing that is impossible to remove in any traditional way. Fortunately, the effect is present for both raw data and for the baseline fit. Subtracting one from the other yields a clean noise signal reliable for physics studies of the reactor. Figure 5 shows the Fourier transform calculated with the above method for the reactor power and baseline and noise. It can be seen that the raw power and baseline both are affected by aliasing in the same amount. The strong peaks at 1 Hz and 2 Hz disappear tracelessly in the difference (the power noise spectrum), proving the efficiency of the method.Fourier transform of reactor power, baseline, and noise. It can be seen that the raw power and baseline both are affected by aliasing in the same amount. The strong peaks at 1 Hz and 2 Hz disappear tracelessly in the difference (power noise spectrum), proving the efficiency of the method.
(a)
(b)A toolkit with spectrum processing various utility codes has been created to help do all the steps required in a modular fashion.


## 4. Conclusions

For low frequency power noise monitoring, it is in certain contexts impossible to remove aliasing frequency components, presampling. Although aliasing collision is not present in the high end of the spectrum, close tofNy—this is being a naturally free region—in the low region, the post-DAQ processing method presented here of raw-baseline subtraction is the only effective means of aliasing collision suppression. Such context appears for instance when the sensors delivering the data have saturation and blind-times after high pulses (for instance of radiation), acting as an effective sampling that is not avoidable through filters. Also, very low-pass filters are difficult to design for 2.5 Hz in the case presented. Digital filters need high sampling rates to avoid the exact condition that is being presented. Analog active filters have a high Q; thus, broadband signals have components that bypass them. Passive very low-pass filters are extremely bulky, as the inductors need to reach into thousands of mH range.The method presented is feasible in real time, via a flash updating procedure of the data on the fly, allowing a physically credible estimate of noise spectral density.It is interesting to note that FFT of noise is a very useful tool also in monitoring the electrical grid, for instance in observing hysteresis noise in transformer cores (leading to magnetic embrittlement and anomalous heating of the core, resp., transformer power losses). Other types of losses (dielectric polarization losses in areas with high humidity) can also be well monitored with this method.Other nuclear reactor types could also benefit from the method, in principle power reactors which have significant water cooling systems, as these reactors produce vacuum reactivity noise if boiling occurs in the pipes.



---

*Source: 821709-2013-12-24.xml*