1. Introduction
Nuclear reactor power-level control is a nonlinear control problem that has been a hot topic in the fields of nuclear engineering for more than two decades, and some exciting results have been given. To design the autonomous control system of the space reactor TOPAZ II, the sliding mode control (SMC) method is adopted for the automatic control functions including reactor startup and shutdown as well as power-level maintaining, ramping, and load-following [
1]. In addition, some high-order SMC designs are presented for pressurized water reactors (PWRs) to mitigate the chattering effect during power maneuvering [
2,
3]. Feedback linearization (FL) is also an effective nonlinear control method, which is interconnected with a robust disturbance observer for the disturbance attenuation reactor power control in [
4]. Although the SMC and FL-based power control can provide strong closed-loop stability, these control laws usually have complicated expressions, inducing the difficulty in the deployment. Actually, proportional–integral–differential (PID)-like control laws still dominate practical control applications in nuclear power engineering. Hence, it is meaningful to develop simple reactor power-level control laws that are able to guarantee strong closed-loop stability performance.
Due to its simple expression and high performance, passivity-based control (PBC) methods have been intensively studied in the last two decades, which are mostly developed for mechanical, electrical, and magnetic systems and their interconnections. The basic idea of PBC is to reshape the storage function of open-loop dynamics to a convex positive-definite function with its setpoint as its minimal point by feedback control. The PBC called control by interconnection (CbI) is proposed in [
5], which gives simple static output control laws by reshaping the energy function via Casimir functions. However, given by its energy-balancing characteristics, the CbI cannot inject extra damping for enhancing the convergence of storage function, which is briefly called the dissipation obstacle [
6]. The interconnection and damping assignment PBC (IDA-PBC) proposed in [
7] avoids the dissipation obstacle based on state-feedback control design of solving a set of first-order partial differential equations, which is more complex in implementation and commissioning with comparison to the CbI. To remove the dissipation obstacle by a simple manner, the power-shaping CbI (PS-CbI) is given in [
8]. The key idea of PS-CbI is to extend the original dynamical system to a higher-order passive system, and the associated extend state vector is the cyclopassive output of original system defined as the function of the control input, the original state vector and its derivative. By adopting necessary adaptation mechanism and injecting extra damping terms, the PS-CbI can further provide asymptotically closed-loop stability [
9]. In addition, the power-shaping PBC given in [
10,
11] is another method being able to overcome the dissipation obstacle, which provides the design of a PBC by reshaping the power function instead of the energy function. Similar to the IDA-PBC, the power-shaping method still relies on solving partial differential equations. To sum up the PBC design methods, it can be seen that: (1) obtaining a storage function with its minimal point giving the expected balance of energy, mass, or momentum is the basis of applying PBC methods; (2) traditional CbI gives simple static output feedback control without extra damping; (3) IDA-PBC and power-shaping PBC can provide extra damping, whose design procedures rely on solving partial differential equations; and (4) PS-CbI provides extra damping by embedding the original dynamics to a higher-order passive system.
The PBC methods have been widely applied to industrial equipment and processes, where the classic CbI and IDA-PBC are the most two popular methods. The classical CbI has been applied to the frequency stabilization of multimachine power system [
12], the motion control of induction motors [
13], and the temperature control of heat exchanger networks [
14], whose popularity comes from its simplicity. The IDA-PBC method is applied for designing the state-feedback control laws for induction motors [
15], hydroturbine systems [
16], multimachine power systems considering hydroturbine with surge tank [
17], a fuel cell/super-capacitor hybrid system [
18], brushless DC motor [
19], and electric vehicles (EVs) [
20]. The popularity of the IDA-PBC method is given by its capability of injecting extra damping, and usually proper state observers are necessary for realizing practically implementable dynamic output feedback control strategies. Moreover, the power shaping PBC has been applied to the control of reaction systems such as the continuous stirred tan reactor (CSTR) [
21].
Since nuclear reactors are complex irreversible thermodynamic process with fission reaction, it is meaningful to develop PBC methods for simply controlling the power level of nuclear reactors, which might enhance the operation performance. Currently, there have been some interesting results on the PBC design of nuclear fission reactors and plants. For pressurized water reactors (PWRs), by giving the storage function from the shifted ectropies of neutron kinetics and reactor thermal hydraulics, a proportional–differential (PD) control is proposed by the use of classical CbI, being able to provide globally asymptotical closed-loop stability [
22]. The parameter uncertainties can be further compensated by coupling an adaptation law with this PD control [
23]. To further improve the performance of PD control, the port-Hamiltonian form of PWR dynamics is given, and an extra damping integration about the weighted sum of the errors in both neutron flux and primary coolant temperature is newly added to the PD law [
24]. Several PBC methods have been given for the modular high-temperature gas-cooled reactors (mHTGRs) [
25]. Based on the classical CbI, a simple static output feedback power control is designed for mHTGRs, which provides globally asymptotic closed-loop stability [
26]. By coupling the IDA-PBC with backstepping, a dynamic output feedback power-level control law is further given, which has stronger performance but is too complicated to be implemented practically [
27]. By considering the nonlinearity from dead zone and saturation as well as the coupled primary pressure dynamics, the simple control given in [
26] can be further enhanced [
28,
29]. Moreover, based on the basic idea of PBC, the coordinated control for a single mHTGR-based nuclear steam supply system (NSSS) module [
30,
31] and that for several NSSS modules coupled by a common steam turbine [
32] are also proposed for the control of the entire plant. The current PBC design for the power control of nuclear reactors are mainly based on the classic CbI method. The result based on the IDA-PBC should be implemented with an additional state observer, which largely strengthens the complexity and difficulty in deployment. To further enhance the closed-loop stability for better power ramping capability in the framework of PBC, extra damping terms related with specific performance indices should be injected.
In this paper, a novel passivity-based power-level control of nuclear reactors is proposed by extending the original power dynamics to a higher-order passive system, where the extra damping terms are given by the extended state variables. The storage function is given by not only the shifted ectropies of neutron kinetics and thermal hydraulics but also the squares of extended state variables. The sufficient condition for globally asymptotic closed-loop stability is given, and the PBC is then applied to the power-level control of a nuclear heating reactor (NHR). Simulation results in the case of power stepping show the feasibility and satisfactory performance.
2. State-Space Model for Control Design
The basic principle of nuclear reactor power level control by using control rods as the actuator is schematically shown in
Figure 1. The inputs of a power-level controller are the errors of neutron flux and average coolant temperature, where the errors are the deviations of the measurements from the setpoints. The output of a power-level controller is the control rod speed signal, which drives the control rods for proper reactivity. The neutron flux is measured by the instruments such as ion chamber detectors, and the measurement of average coolant temperature is the algebraic mean of the measurements of reactor inlet and outlet coolant temperatures. It is the basic requirement for a nuclear reactor power-level controller to suppress or eliminate the errors while guaranteeing satisfactory closed-loop stability. In this section, the nonlinear state-space model of point reactor dynamics is first given. Then, a port-Hamiltonian form of nuclear reactor dynamics is proposed, and the control problem to be solved is finally raised.
The nuclear reactor dynamics for control design is usually given by the point kinetics with an equivalent group of delayed neutrons and the reactivity feedback from the averaged fuel and coolant temperatures that can be given as
where
nr is the normalized neutron flux,
cr is the normalized concentration of delayed neutron precursors,
β is the fraction of delayed neutrons,
Λ is the effective prompt neutron lifetime,
λ is the radioactive decay constant of delayed neutron precursor,
αf and
αc are respectively the reactivity feedback coefficients of the fuel and coolant temperatures,
Tf and
Tcav are respectively the average temperatures of fuel assemblies and primary coolant,
Tcav,m and
Tf,m are respectively the initial steady values of
Tcav and
Tf,
Tcin is the primary coolant temperature at reactor core inlet,
Ω is the heat transfer coefficient between the fuel and primary coolant,
M is the primary heat capacity flow-rate, i.e., the multiplication of primary mass flowrate and specific heat capacity at constant pressure,
P0 is the rated reactor thermal power,
μf and
μc are respectively the total heat capacities of fuel assemblies and primary coolant, and
ρr,
Gr, and
zr are the induced reactivity, differential worth, and speed of control rods, respectively. It is not loss of generality to assume that scalar
αf is negative, which should be satisfied by all the current nuclear reactor designs to avoid severe accidents.
Let nr0, cr0, Tf0, Tcav0, Tcin0, and ρr0 be the setpoints of reactor process variables nr, cr, Tf, Tcav, Tcin, and ρr, respectively, where the setpoints are taken as the steady-state solution of (1) corresponding to the power setpoint. Then, the deviations of process variables from their setpoints can be given as δnr = nr − nr0, δcr = cr − cr0, δTf = Tf − Tf 0, δTcav = Tcav − Tcav0, δTcin = Tcin − Tcin0, δρr = ρr − ρr0.
We define
as the state variables while taking
as the control input. Here, it is assumed that δ
Tcin = 0, which means that the influence of secondary-loop is omitted. Based on (2)–(4), the nonlinear state-space model of nuclear reactor dynamics can be written as
where
y is the measurement output, and vector functions
f(
x),
g(
x) and
h(
x) are given by
From (5), the evolutional trend of state vector x is given by both the neutron kinetics and reactor thermal hydraulics, while that of the state variable ξ is determined by control rod dynamics. To guarantee desirable operation performance, it is meaningful to design a control law for nonlinear system (5) so that x→O as t→∞.
3. Passivity-Based Control Design
The PBC design for nuclear reactor power-level dynamics given by (5) is summarized as the following theorem.
Theorem 1. Consider nonlinear system (5) with control input u determined bywhereis the auxiliary state vector, ζi(0) = 0 (i = 1, 2), gains kD1, kD4, kP1 and kP4 as well as damping coefficients κ1 and κ2 are positive constants. The extended state vector z defined byconverges to the origin globally and asymptotically ifwhere ci (i = 1, 2) are arbitrarily given positive constants. Proof. From model (5), it is reasonable to regard
ξ as a virtual control, and design the stabilizer for subsystem
The shifted ectropies for the neutron kinetics and reactor thermal hydraulics are given by positive definite functions
and
respectively. By differentiating (15) and (16) along the trajectory of subsystem (14), it can be found that the time derivatives of
EN and
ET satisfy
and
respectively. □
Design virtual control
ξ as
where positive constants
kD1 and
kD4 are feedback gains, and both
ζ1 are
ζ2 are the auxiliary state variables to be designed.
Substitute control law (19) to Equation (17), and it can be seen that the time derivative of
EN is given by
Design auxiliary state variables
ζi (
i = 1, 2) as
where positive constants
κi (
i = 1, 2) are given damping coefficients and
ζi(0) = 0 (
i = 1, 2).
Define Lyapunov function candidate of extended system
as
Since both functions EN and ET are positive definite and feedback coefficient αf is strictly negative, it can be seen that function V(z) is positive definite. Differentiate V(z) defined by (24) along the trajectory of the closed-loop system formed by (23) and (19):
From (25), it can be seen that if both inequalities (12) and (13) are satisfied, then
which means that
z→
O as
t→∞, i.e., the extended system defined by (23) is globally asymptotically stable.
Furthermore, from the relationship
the control input
u can be given as
By interconnecting (27) with (21) and (22), passivity-based control (9) can then be obtained. This completes the proof of this theorem.
Remark 1. The system considered in the theorem is the nuclear reactor power dynamics, and the designed control law is for the regulation of reactor power-level. The control input u is just the control rod speed, i.e., the output of the control law. State variables x1 and x4 are just the error of normalized neutron flux and reactor coolant temperature. The control input u is essentially given by x1 and x4. By choosing storage function of nuclear reactor power dynamics asit can be derived thatwhereFrom (28) and (30), the reactor power dynamics is passive, which is the reason to call control law (27) as a PBC. Remark 2. From (21) and (22), it can be seen that auxiliary state variables ζ1 and ζ2 are the weighted integration of x1 and x4, respectively. Hence, virtual control law (19) can be seen as a proportional–weighted–integral control law.
Remark 3. From PBC (9), control input u can be written in S-domain aswhere It can be seen from (31)–(33) that the controller takes a simple form, which leads to easy implementation and commissioning in the practical engineering. Furthermore, it can be seen that PBC (9) degenerates to classical PD control if κ1 = κ2 = 0.
4. Simulation Results with Discussions
In this section, PBC (9) is applied to the power-level control of a 200MW
th nuclear heating reactor (NHR), which is a typical integral pressurized water reactor (iPWR) with a series of advanced design features such as integral arrangement, self-pressurization, full-power-range natural circulation, passively decay heat removal, and a built-in hydraulic control rod driven mechanism [
33,
34,
35]. The main reactor design parameters are given in
Table 1.
The simulation is performed on Matlab/Simulink. Feedback gains kP1, kD1, kP4, kD4, and damping coefficients κ1 and κ2 should be chosen so that both inequalities (12) and (13) are satisfied. According to (12) and (13), it is set that kP1 = 1, kD1 = 10, kP4 = 0.0001, kD4 = 0.001, and different values of damping coefficients κ1 and κ2 are adopted to show their influence on control performance.
In the simulation, the case of power stepping from 100% full power (FP) to 40% FP is considered. The transient responses of normalized neutron flux
nr, normalized precursor concentration
cr, averaged fuel temperature
Tf, average coolant temperature
Tcav, primary coolant temperatures at reactor core outlet
Tcout and inlet
Tcin, and primary coolant flowrate
Gp as well as control rod speed
vr with different
κ1 and
κ2=1 and those with different
κ2 and
κ1=1 are shown in
Figure 2 and
Figure 3, respectively.
From the simulation results, it can be seen that the closed-loop stability is well guaranteed even under different values of damping coefficients
κ1 and
κ2, which verifies the correctness of the theoretic result about closed-loop stability analysis. Furthermore, the dynamic responses given in
Figure 2 and
Figure 3 also show the influence of damping coefficients
κ1 and
κ2 quantitatively. From
Figure 2, it can be seen that the influence of
κ1 to the transient performance is limited, which is caused by the low inertia in neutron kinetics. It can be seen from
Figure 3 that the control performance is sensitive to the value of damping coefficient
κ2. Actually, the total heat capacity of a nuclear reactor is large, which gives a strong inertia in the reactor coolant temperature dynamics. Damping coefficient
κ2 is smaller, the integration effect of average coolant temperature error is stronger, i.e., auxiliary state variable
ζ2 is larger, which leads to a shorter transition period of average coolant temperature while causing a larger overshoot of neutron flux. Hence, in practical engineering, the value of
κ2 should be carefully tuned so as to obtain a satisfactory balance between the transition period of coolant temperature and the overshoot of neutron flux.
Finally, it can be seen from Equations (31)–(33) that the PBC has a simple expression that not only determines an easy implementation on the current digital control system platforms but also gives an easy tuning and commissioning procedure practically.