KT110535
KT110535
net/publication/283537363
CITATIONS READS
4 365
3 authors:
N. Ashoub
Egyptian Atomic Energy Authority
15 PUBLICATIONS 24 CITATIONS
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
All content following this page was uploaded by H. I. Saleh on 21 December 2017.
A. Arafa et al.: Design and implementation progress of multi-purpose simulator for nuclear research reactor using LabVIEW
A. Arafa et al.: Design and implementation progress of multi-purpose simulator for nuclear research reactor using LabVIEW
the NRR is the reactor core, which contains the fuel elements, Afterwards, the development process could be started and
coolant, moderator and control rods. The other part is the for this simulator it consists of three phases. Firstly, the mod-
coolant system, shown in Fig. 1, which contains the primary eling phase is the process of building the mathematical mod-
and secondary loops, heat exchanger and pumps. The core is els using design data of the NRR in question. Secondly, the
composed of channels that enclose the fuel elements. These building the simulation system phase develops the package it-
channels are connected to the lower plenum at the bottom of self. Finally, the verification and validation phase assures the
the core inlet, and to the upper plenum at the core exit. The constructed models as well as generates different accident
flow enters the channels at the bottom through an inlet ori- and transient scenarios. Figure 2 shows the different phases
fice. Under steady state (normal operation) conditions, the of the simulator developing process, and illustrates the impact
core inlet flow is sub-cooled in the cooling tower. As the flow of the simulator design on the first two phases.
in the channel passes by the fuel, it is heated and to preserve Many languages can be used to accomplish the develop-
the temperature at a certain predetermined value, pumps ment process. For this simulator, LabVIEW 11.0 installed on
were established in the recirculation loops. The simplest con- a personal computer with Windows OS was chosen as the de-
figuration was with external pumps suctioning the fluid from velopment platform for the mathematical models and GUI.
the down-comer region and injecting it into the lower plenum The \Modeling" phase utilizes the LabVIEW modeling lan-
of the core. guage for putting the mathematical models of different reac-
This work is a part of the nuclear research reactor simulator tor core processes, and instrumentation and thermo-hydraulic
project whose development environment and tools, GUI, ar- systems. Moreover, the \Building simulation system" phase
chitecture and design overview were introduced in [11]. The constructs the graphical user interface (GUI), which inter-
main goal of this paper is to illustrate the applied neutronic faces the user inputs/outputs of the modeling phase, through
and thermo-hydraulic models and the developing phases and the graphical building blocks (i. e. controls and indicators) of-
progress of the multipurpose NRR simulator. Section 2 intro- fered from LabVIEW. On the other hand, the \Verification"
duces the three development phases and the simulator struc- phase verifies the mathematical model according to the
ture. Section 3 provides the mathematical models and their NRR design and operation (experimental) data. Besides, it
LabVIEW implementations. In Section 4 the verification and validates the simulator software to assure that it satisfies the
validation of the simulator is discussed, and in Section 5 its fu- design requirements. The LabVIEW offers a \Real-Time Ex-
ture development is suggested. Conclusions are given in Sec- ecution Trace Toolkit", which helps in preparing error reports
tion 6. by capturing and displaying the timing and event data of vir-
tual instrument (VI) file and thread events for LabVIEW
modules [7]. In turns, these reported problems are used to
2 Simulator development perform source code analysis and software database analysis
to facilitate code maintenance and fixing those problems.
Before starting the development process, a full design vision The simulator structure consists of two main systems: the
should be founded. This design contains the NRR technical \Instrumentation and monitoring", and \Reactor Modeling"
data and actual specifications, and the simulation require- systems. The Reactor Modeling systems include the \Control
ments (i. e., the structured creation of software components & Reactivity", \Neutronic-Kinetic", and \Thermo-hydraulic"
to implement specific functionality and how a system will ac- systems. Figure 3 shows the simulator block diagram where
complish the desired functionality). the different building processes of each modeled systems are
2 80 (2015) 5
KT_KT110535 – 1.9.15/stm media köthen
A. Arafa et al.: Design and implementation progress of multi-purpose simulator for nuclear research reactor using LabVIEW
80 (2015) 5 3
KT_KT110535 – 1.9.15/stm media köthen
A. Arafa et al.: Design and implementation progress of multi-purpose simulator for nuclear research reactor using LabVIEW
described, and the relations between these different modeled 3.1 Neutronic-kinetic modeling
processes and systems are illustrated. Each process of mod-
eled systems has been built using LabVIEW as a subroutine The simulator uses the point kinetics equations to determine
program which is a separate virtual instrument (VI) file. Each the reactor power [13], [15] and allows up to 6 groups of de-
VI subroutine has its front panel which inputs or outputs the layed neutron data. The point kinetics model has been widely
related parameters and controls. used in the analysis of the transient behavior of a nuclear re-
actor. While, the point kinetics Eqs. (1-2) have been derived
from the general neutron balance equations where the time,
3 Reactor mathematical model based on LabVIEW space, energy and direction-dependent neutronics are taken
into account. These equations are the neutron dynamics and
A nuclear reactor is a high complex dynamic system, where the delayed neutron emitting nuclei dynamics as follow:
its components are physically interconnected and dynamically
dn q % b X6
interacting. During the simulator’s developing, it is important ¼ nþ ki Ci ð1Þ
that computer codes should be built to give support to other dt < i¼1
preliminary developments, like the system general concep-
tion, control and instrumentation possibilities, etc. These dCi b i
¼ n þ ki Ci ð2Þ
codes have to follow the dynamic models of the nuclear reac- dt <
tor system which were illustrated in the literatures [13, 14].
where t is the time, q is the total reactivity, b is the effective
One-dimensional thermal hydraulic equations coupled with
delayed neutron fraction, b i is the effective ith group delayed
point kinetics equations have been solved to provide an ap-
neutron fraction, < is the prompt neutron generation time,
propriate model for the NRR system. Simplified and faster
Ci is thtformate concentration of delayed neutron precursors
models are required to enable the utilization of more familiar
of group i, ki is the decay constant of delayed neutron precur-
languages such as MATLAB and LabVIEW.
sors of group i and n is the neutron density as a function of
For this work, the LabVIEW is selected as an appropriate
time.
alternative for the development of a dynamic computer code.
To solve the point kinetics equations, the reactor is consid-
LabVIEW provides advantages in the execution time, easy
ered to be a spatially homogeneous lumped parameter system
handling tools and objects, and flexible creation and modifica-
and its model is chosen to be a time-dependent, linear, single-
tion of models. As well as, LabVIEW allows scientific compu-
group neutron diffusion equation [13]. The point kinetics
tation through its powerful calculation tools. The following
equations implementation is shown in Fig. 4.
subsections discuss the point kinetics calculations with the
fully coupled thermal fluids and fuel temperature feedback Through the fission process, some of produced fission pro-
implemented into the LabVIEW framework. ducts have highly absorbed cross sections for thermal neu-
4 80 (2015) 5
KT_KT110535 – 1.9.15/stm media köthen
A. Arafa et al.: Design and implementation progress of multi-purpose simulator for nuclear research reactor using LabVIEW
> ; < ?
trons, so that it’s called poisons. There are two main types of 000 Le pz pL 1 pz
Dt ¼ Tout % Tin ¼ qo Ac sin þ sin þ cos . . . ð4Þ
poisons; the decay chains of Xenon-135 (135 Xe) and Samar- pmcp Le 2Le 2pRco h Le
ium-149 (149 Sm). 135 Xe and 149 Sm are very important and sig- 000
nificant neutron poisons as they have very large thermal ab- where qo is the Volumetric Heat Thermal source in kw/m3, Tin
sorption cross sections. The 135 Xe isotope is formed as a is the Inlet coolant temperature in 8C, Ac is the Channel flow
result of the decay of Iodine-135 (135 I) and also produced di- area in m2, L is the Core height, Le is the Extrapolated core
rectly from the fission. The 149 Sm isotope is of much less con- height, cp is the Specific heat of the coolant in KJ/kg · 8C, m
cern in a reactor than 135 Xe and often included separately in is the mass flow rate of the coolant (m3/h) and Tout is the out-
reactor calculations. 149 Sm is not formed directly in fission let coolant temperature. Figure 7 illustrates the LabVIEW im-
but appears as the result of the decay of the 149 Pm (pro- plementation of the general heat transfer between the fuel
methium-149). The forming of these fission products are mod- and the coolant which is given by Eq. (4) shows a graphical
eled in a LabVIEW VI file which was described in [11]. Fig- representation for the relations between the input and output
ure 5 shows LabVIEW block diagrams of the poisons variables of the two main cooling systems indicating different
equations that consist of the main two loops (i. e. Xenon and parameters used in calculating the temperature difference Dt.
Samarium loops) and the related subsystem for each loop.
3.3.2 Heat exchanger model
3.2 Control and reactivity modeling
Assuming a perfect thermal mixture between the fluid in the
Nuclear reactor dynamics concern with reactor power transi- primary loop and the fluid entering the secondary loop of
ents where the power change causes the reactivity change cooling; the heat exchanger equations are described as fol-
feedbacks. These reactivity feedbacks play an important role lows:
in overall reactor behavior. Tin ¼ Tout % gðTout % Tcold Þ ð5Þ
The reactivity dependence on the rod position is assumed
to be sinusoidal and it is assumed to be linear with the fuel wa
and moderator temperature. Then, the reactivity relates with Thot ¼ Tcold % ðTin % Tout Þ ð6Þ
wb
the temperature and the control rod position [16] as follows:
% & where
q ¼ (f Tf % Tfo þ (m ðTm % Tmo Þ þ
;; < < 1 % e%(
%x þ z 2:5 9 x þ z: g¼ wa %( ; ð7Þ
þb þ sin 2p þd ð3Þ 1% e
2 4p 2:5 wb
80 (2015) 5 5
KT_KT110535 – 1.9.15/stm media köthen
A. Arafa et al.: Design and implementation progress of multi-purpose simulator for nuclear research reactor using LabVIEW
Fig. 5. The poisons equations implementation with its subsystem block diagrams
6 80 (2015) 5
KT_KT110535 – 1.9.15/stm media köthen
A. Arafa et al.: Design and implementation progress of multi-purpose simulator for nuclear research reactor using LabVIEW
Fig. 6. The LabVIEW implementation of the reactivity showing its dependence on the rod position and the temperature of the fuel and moderator
Fig. 7. LabVIEW implementation of the general heat transfer between the fuel and the coolant
80 (2015) 5 7
KT_KT110535 – 1.9.15/stm media köthen
A. Arafa et al.: Design and implementation progress of multi-purpose simulator for nuclear research reactor using LabVIEW
The reactor mathematical models and simulator are verified In this exercise, the simulator is initialized to the full power at
and validated through the checking of the correctness of their the steady state normal operation. The axial variations of
outputs. The verification compares the conceptual model to coolant, cladding surface, fuel element center temperatures
the computer representation which implements that concep- at the steady state are illustrated in Fig. 10. Maximum tem-
tion. Besides, it ensures that the model is correctly implemen- peratures for coolant, clad and the fuel occur above the mid
ted in the computer and the input parameters and logical plane of the axial distance. The effect of the coolant is to shift
structures of the model are precisely represented. upward the radial temperatures along the fuel element. The
Various exercises are conducted to observe system behavior coolant temperature difference is achieving the designed val-
and its congruence to the literature standards. Other exercises ue. The thermal hydraulic parameters obtained from the si-
describe and explain the response of the system to certain op- mulators are coincidence with the theoretical model for the
eration modes and conditions. This section points to some ex- simulated reactor type. The operator keeps following these
ercises that conducted to operate all the inputs and find out screens illustrated in Fig. 10, which indicate the stability of
effects of such a system in action. the system and the nuclear safety margin of the core compo-
nents.
Fig. 10. Simulator’s screens displaying the distribution of Coolant, Clad-surface, and Fuel central-line Temperature with respect to axial distance in-
side the core
8 80 (2015) 5
KT_KT110535 – 1.9.15/stm media köthen
A. Arafa et al.: Design and implementation progress of multi-purpose simulator for nuclear research reactor using LabVIEW
Response to reactor trip and poison-out responding variations in the temperature of the coolant at
the core outlet.
This validation exercise shows the response of the simulator
to the reactor trip in which the simulator is set to \manually
reactor trip" and initialized to the full power. The response
of the overall unit is observed to be compared with the stan- 5 Future development
dard safety reports. After waiting five hours of operation the
reactor trip is shutting down. Therefore the rods statuses, The simulator was reevaluated and the possibilities for the fu-
shown in the simulator layout, are totally immersed. The reac- ture concepts were examined. Next steps of development of
tor could not be restarted unless the reactor poison-out time this simulator may be devoted to creation and implementa-
passed. It approximately takes 20 hours to re-set the reactor tion of additional interactive features for operator actions,
trip. Figure 11 shows the variations of Xenon and Iodine den- taking into account results of qualitative and quantitative as-
sities with respect to time during the operation of this scenar- sessments. Further accident scenarios \Mal functions" have
io, as well as the poison-out time is indicated. to be investigated and evaluated to achieve one of the aimed
purposes of this simulator. As well as, the simulation of other
Response to power maneuver reactor type designs will be developed.
Moreover, the simulator interfaces are planned to make
This validation exercise observes the behavior of the reactor use of the standard web browser, so that it opens up the possi-
for a power maneuver, and validates all parameters responses. bility of collaborative work on the computations side of the si-
Whereas, the simulator is initialized to the 100 % of the full mulator. Therefore, the calculations of the Neutronic and
power in which all observed parameters show their consis- thermal hydraulics models could be developed and run on dif-
tency. The power is reduced to 90 % of the full power at a rate ferent computers, allowing a pooling of programming and
of 1.0 % per sec. After the 10 seconds transient-time of power hardware resources. From the educational point of view, that
reduction, the simulator is running at the reduced power for feature will allow individual instructors to tailor the program
approximately one minute. The main change observed is the to specific learning objectives.
reducing of the coolant temperature output by one degree. The work in this simulator open many research fields such
The power is continuing raise to return to 100 % of the full as speeding up the simulation time in the calculation of math-
power at a rate of 1.0 %FP/sec. Therefore, the coolant tem- ematical models which need the parallel computing. The abil-
perature output is slowly increased to its previous saturated ity of using heterogeneous parallel computing will permit in-
value of 37.4. Figure 12 presents the power maneuver in terms vestigating more advanced mathematical model which in
of normalized excess power with respect to time and the cor- turn allows more accurate and convenient results.
80 (2015) 5 9
KT_KT110535 – 1.9.15/stm media köthen
A. Arafa et al.: Design and implementation progress of multi-purpose simulator for nuclear research reactor using LabVIEW
10 80 (2015) 5
KT_KT110535 – 1.9.15/stm media köthen
A. Arafa et al.: Design and implementation progress of multi-purpose simulator for nuclear research reactor using LabVIEW
80 (2015) 5 11