Next Article in Journal
Analysis of the Channel Influence to Power Line Communications Based on ITU-T G.9904 (PRIME)
Previous Article in Journal
Developing a New HSR Switching Node (SwitchBox) for Improving Traffic Performance in HSR Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Model-Free Coordinated Control for MHTGR-Based Nuclear Steam Supply Systems

1
Institute of Nuclear and New Energy Technology (INET), Tsinghua University, Beijing 100084, China
2
Collaborative Innovation Center of Advanced Nuclear Energy Technology of China, Beijing 100084, China
Submission received: 1 December 2015 / Revised: 5 January 2016 / Accepted: 6 January 2016 / Published: 11 January 2016

Abstract

:
The modular high temperature gas-cooled reactor (MHTGR) is a typical small modular reactor (SMR) that offers simpler, standardized and safer modular design by being factory built, requiring smaller initial capital investment, and having a shorter construction period. Thanks to its small size, the MHTGRs could be beneficial in providing electric power to remote areas that are deficient in transmission or distribution and in generating local power for large population centers. Based on the multi-modular operation scheme, the inherent safety feature of the MHTGRs can be applicable to large nuclear plants of any desired power rating. The MHTGR-based nuclear steam supplying system (NSSS) is constituted by an MHTGR, a side-by-side arranged helical-coil once-through steam generator (OTSG) and some connecting pipes. Due to the side-by-side arrangement, there is a tight coupling effect between the MHTGR and OTSG. Moreover, there always exists the parameter perturbation of the NSSSs. Thus, it is meaningful to study the model-free coordinated control of MHTGR-based NSSSs for safe, stable, robust and efficient operation. In this paper, a new model-free coordinated control strategy that regulates the nuclear power, MHTGR outlet helium temperature and OTSG outlet overheated steam temperature by properly adjusting the control rod position, helium flowrate and feed-water flowrate is established for the MHTGR-based NSSSs. Sufficient conditions for the globally asymptotic closed-loop stability is given. Finally, numerical simulation results in the cases of large range power decrease and increase illustrate the satisfactory performance of this newly-developed model-free coordinated NSSS control law.

1. Introduction

Small modular reactors (SMRs) are those nuclear fission reactors with electrical output power less than 300 MWe, and could be beneficial in providing electricity to the areas without transmission or distribution infrastructure in generating local power for a large population center and in being viable for specific applications such as heat sources for industrial complexes [1,2]. Due to the low power density and large heat capacity, some SMRs even have the inherent safety feature which is the most advanced feature of these SMR-based plants relative to those conventional nuclear plants and also prevents these SMRs from the hazards of core-melting, radiological release and loss of coolant accident (LOCA) [3,4]. Due to the ability of offering simpler, safer and standardized modular design by factory-building, less initial capital investment, and shorter construction period, SMRs have been viewed by International Atomic Energy Agency (IAEA) as a developing trend of nuclear energy. Through adopting multi-modular scheme, i.e., multi-SMRs providing steam for one turbine/generator set, the inherent safety feature can be applicable to the large power plants of any desired power ratings. It is well known that the modular high temperature gas-cooled reactor (MHTGR, such as the HTR (High Temperature Reactor)-Module designed in Germany [5]), which uses helium as coolant and graphite as both moderator and structural materials, is a typical SMR with a strong inherent safety feature [6,7]. An MHTGR-based nuclear steam supplying system (NSSS) module composed of a one-zone pebble-bed MHTGR, a helical-coil once-through steam generator (OTSG) arranged side-by-side with the MHTGR and some connecting pipes is shown in Figure 1, which is one module of the under-constructed two-modular plant HTR-PM (Pebble-Bed Module) designed by Institute of Nuclear and New Energy Technology in Tsinghua University [8].
Figure 1. Schematic view of a modular high temperature gas-cooled reactor (MHTGR)-based nuclear steam supplying system (NSSS) module.
Figure 1. Schematic view of a modular high temperature gas-cooled reactor (MHTGR)-based nuclear steam supplying system (NSSS) module.
Energies 09 00037 g001
It is very clear that both the MHTGR and OTSG are complex nonlinear dynamic systems, and their dynamics are tightly coupled with each other. To improve both the steady and transient performance, it is necessary to develop coordinated control methods for the NSSS composed of the MHTGR and OTSG by considering the nonlinear dynamics. There are still very few results in the coordinated control of the reactor and steam generator. However, there have been some promising results in nonlinear control design methods for the reactors or the steam generators (SGs), which are helpful for developing coordinated control strategies. The sliding mode control (SMC) method is an effective approach to realize nonlinear control design for nuclear reactors. Shtessel gave a nonlinear power-level controller for space reactor TOPAZ II through the SMC method [9]. Due to the intrinsic chattering effect of the SMC, the recursive sliding mode control (RSMC) method was introduced by Huang, Edwards and Lee to solve the chattering problem and to improve the robustness. Based on the RSMC method, a fuzzy-adapted nonlinear power-level control law was designed for the advanced boiling water reactors (ABWR) [10]. The physics-based control (PBC) method is also an effective way to design nonlinear reactor control laws by retaining or strengthening stable subdynamics and by cancelling or suppressing unstable subdynamics, which has been applied to the load-following control design for the PWRs (Pressurized Water Reactors) [11,12,13] and MHTGRs [14,15] recently. From the aspect of SG controller design, some advanced control approaches such as the model predictive control (MPC) [16], neurofuzzy [17] and feedback-dissipation [18] methods have all been applied to improve the regulation performance of SGs.
Although the above control strategies can provide large-scale closed-loop stability for the nuclear reactors or steam generators. However, these control laws have the drawback of heavy model-reliance. Some control laws needs the state prediction or observation directly based upon the reactor or SG models [9,10,11,12,14,15,16,17], and the control parameters also vary with the reactor or SG physical and thermal-hydraulic parameters [9,10,11,12,13,14,15,16,17,18], which leads to the difficulty in practically tuning and implementation. Therefore, it is necessary to develop model-free control methods to relieve the heavy model reliance of the designed control strategies. In the early 2000s, by applying the stable adaptive fuzzy control theory presented in [19], Marseguerra and Zio gave a model-free reactor control [20]. However, the fuzzy system of this model-free controller should be either suggested by experienced operators or trained online by measurement data, which is expensive especially for the newly-built reactors. Moreover, the fuzzy adaptive nonlinear control approach given in [19] is only suitable to those nonlinear systems in Brunovsky normalized form, which lose the prior knowledge about the structure of the plant dynamics given by both neutron kinetics and thermal-hydraulics. Thus, it is necessary and meaningful to study the model-free NSSS coordinated control design approach based on the inherent physical structure of the reactor and SG dynamics, which is not only for the freedom from the expression of dynamic model but also for insensitivity to the physical and thermal-hydraulic parameters. Actually, using inherent dynamical structure given by the energy or ectropy is a crucial way to design model-free control laws. For example, based upon the energy method, Rigatos gave a model-free controller of flexible-link robots [21], which is free from the dynamic model, and is also insensitive to the robot parameters. Thus, a physics-based control design method may also be an effective way to realize model-free coordinated control for MHTGR-based NSSSs.
Motivated by the necessity of developing model-free coordinated control method for NSSSs, a model-free coordinated controller of the NSSS composed of one MHTGR and one OTSG are proposed through the physics-based approach in this paper. This NSSS coordinated control law has three key features: (a) there is no representation of any dynamic models for generating control input, state-observation or prediction; (b) the controller parameters are insensitive to the NSSS physical or thermal-hydraulic parameters; and (c) the closed-loop is globally asymptotically stable. Numerical simulation results in the cases of large range power increase and decrease illustrate the satisfactory regulating performance of this NSSS coordinated controller.

2. Modeling and Problem Formulation

In this section, the state-space model of the MHTGR-based NSSS for coordinate control design is firstly given, and then the theoretic problem to be solved in the following sections is formed.
From Figure 1, the MHTGR and OTSG in one NSSS are arranged side by side, and are connected with each other by a horizontal coaxial hot/cold gas duct. The cold helium enters the blower that mounted on the upper part of the OTSG, and is pressurized before flowing into the cold gas duct. It enters the channels inside the reflector from bottom to top for cooling the reflector, and then passes through the pebble-bed from top to bottom where it is heated to a very high temperature. The hot helium leaves the hot chamber inside the bottom reflector, and then flows into the OTSG primary side where it is cooled to be the cold helium by transferring its heat to the secondary water/steam flow. By adopting the point kinetics with one equivalent delayed neutron group and by considering the temperature feedback effect of the reactor core composed of the pebble-bed and reflector, the dynamic model for coordinated control design is:
{ Λ n ˙ r = ( ρ r β ) n r + β c r + α R n r ( T R T R , m ) , c ˙ r = λ ( n r c r ) , μ R T ˙ R = Ω P ( T R T H ) + P 0 n r , μ H T ˙ H = Ω P ( T R T H ) Ω S ( T H T S ) , μ S T ˙ S = Ω S ( T H T S ) 2 M S ( T S T Sin ) , ρ ˙ r = G r v r ,
where the meaning of the variables are given in Nomenclature. Here, it is worthy to be noted that αR is guaranteed to be negative by the physical design.
Define deviations of nr, cr, TR, TH, TS, TSin and ρr from their equilibrium values, i.e., nr0, cr0, TR0, TH0, TS0, TSin0 and ρr0, respectively, as δnr = nrnr0, δcr = crcr0, δTR = TRTR0, δTH = THTH0, δTS = TSTS0, δTSin = TSinTSin0, δρr = ρr − ρr0. Since TSin is usually given by the environment, assume that:
δ T Sin 0
Because the helium flowrate can be adjusted by changing the main blower speed, it is not a loss of generality to assume that:
δ Ω S = Ω S Ω S 0 = κ S δ W P
where ΩS0 is the steady value of ΩS, and δWP is the variation of the primary helium flowrate relative to its steady value, and κS is a positive scalar.
Moreover, the feedwater flowrate can also be adjusted, and thus we assume that:
δ M S = M S M S 0 = C S δ W S
where MS0 is the steady value of MS, δWS is the error of the feedwater flowrate relative to its steady value, and CS is a positive scalar denoting the heat capacity of secondary coolant.
Define:
x = [ x 1 x 2 x 3 x 4 x 5 ] T = [ δ n r δ c r δ T R δ T H δ T S ] T
ξ = δ z r = G r 1 δ ρ r
u = [ u 1 u 2 u 3 ] T = [ v r δ W P δ W S ] T
and
y = [ x 1 x 4 x 5 ] T
where x is called the NSSS state-vector, u is called the control input, y is the output that can be obtained by measurement, and δzr is the total displacement of the control rods.
Then, nonlinear state-space model for coordinated control design can be written as:
{ x ˙ = f ( x ) + g 1 ( x ) ξ + g 2 ( x ) u 2 + g 3 ( x ) u 3 , ξ ˙ = u 1 , y = h ( x ) ,
where
f ( x ) = [ Λ 1 β ( x 1 x 2 ) + Λ 1 α R ( n r 0 + x 1 ) x 3 λ ( x 1 x 2 ) μ R 1 P 0 x 1 μ R 1 Ω P ( x 3 x 4 ) μ H 1 Ω P ( x 3 x 4 ) μ H 1 Ω S 0 ( x 4 x 5 ) μ S 1 Ω S 0 ( x 4 x 5 ) 2 μ S 1 M S 0 x 5 ]
g 1 ( x ) = [ Λ 1 G r ( n r 0 + x 1 ) O 1 × 4 ] T
g 2 ( x ) = [ O 1 × 3 μ H 1 κ S ( T H T S ) μ S 1 κ S ( T H T S ) ] T
g 3 ( x ) = [ O 1 × 4 2 μ S 1 C S T S ] T
and
h ( x ) = [ x 1 x 4 x 5 ] T
From the MHTGR-based NSSS state-space model Equation (9), the theoretic problem to be solved in the next sections is summarized as:
Problem 1. How can a model-free control input u defined by:
u = [ u 1 u 2 u 3 ] T
for nonlinear system Equation (9) be designed so that the closed-loop is globally asymptotically stable?

3. Model-Free Coordinated Control Design

The following Theorem 1, which is the main result of this paper, summarizes the design of model-free coordinated control input u for the MHTGR-based NSSS. The corresponding design process can be seen in the proof Theorem 1.
Theorem 1. Consider nonlinear system Equation (9). The coordinated control given by:
u 1 = ( k p 11 x 1 + k d 11 x ˙ 1 + k p 14 x 4 + k d 14 x ˙ 4 + k p 15 x 5 + k d 15 x ˙ 5 )
u 2 = Ω S 0 κ S ( T H T S ) ( k p 24 x 4 k p 25 x 5 )
and
u 3 = M S 0 C S T S k p 35 x 5
where feedback gains kp11, kd11, kp14, kd14, kp15, kd15, kp24, kp25 and kp35 are all positive constants satisfying:
k d 14 = ( 1 η 1 ) ( γ 1 + γ 2 ) | α R | μ H G r μ R
k d 15 = ( 1 η 2 ) [ γ 1 + Ω S 0 ( 1 + k p 25 ) Ω S 0 ( 1 + k p 25 ) + 2 M S 0 ( 1 + k p 35 ) γ 2 ] | α R | μ S G r μ R
k p 14 = γ 2 ( 1 + k p 24 ) 2 M S 0 ( 1 + k p 35 ) Ω S 0 ( 1 + k p 25 ) + 2 M S 0 ( 1 + k p 35 ) | α R | Ω S 0 G r μ R
and
k p 15 = γ 1 ( 1 + k p 35 ) 2 | α R | M S 0 G r μ R
respectively, and γi and ηi (i = 1, 2) are all given positive constants with 0 < ηi < 1. Then, the closed-loop system formed by Equations (9) and (16)–(18) is globally asymptotically stable if:
k d 11 > | α R | P 0 G r Ω P ( γ 1 + γ 2 ) 2 [ 1 + η 1 μ H 4 μ R + η 2 μ S ( 1 + k p 24 ) 4 μ R ( 1 + k p 25 ) ]
k p 24 > 2 Ω P Ω S 0 ( 1 + η 1 μ H μ R ) 1
and
k p 35 > Ω S 0 2 M S 0 ( 1 + k p 25 + η 2 Ω P μ S Ω S 0 μ R ) 1
Proof: See Appendix.
Remark 1. Define
κ p 24 = Ω S 0 κ S ( T H T S ) k p 24
κ p 25 = Ω S 0 κ S ( T H T S ) k p 25
and
κ p 35 = M S 0 C S T S k p 35
From Equations (26)–(28), as well as inequalities Equations (24) and (25), we can easily obtain that both
κ p 24 > [ 2 Ω P Ω S 0 ( 1 + η 1 μ H μ R ) 1 ] Ω S 0 κ S ( T H T S )
and
κ p 35 > { Ω S 0 2 M S 0 [ 1 + κ S Ω S 0 κ p 25 ( T H T S ) + η 2 Ω P μ S Ω S 0 μ R ] 1 } Ω S 0 κ S ( T H T S )
should be satisfied for guaranteeing globally asymptotic closed-loop stability. Moreover, control input u2 and u3 can be certainly rewritten as:
u 2 = κ p 24 x 4 κ p 25 x 5
and
u 3 = κ p 35 x 5
respectively.
Remark 2. Since γi (i = 1, 2) are arbitrarily given positive constants, from Equations (21) and (22), it is clear that feedback gains kp14 and kp15 can be arbitrary positive values. Moreover, since ηi (i = 1, 2) are arbitrary given constants satisfying 0 < ηi < 1, from Equations (19) and (20), feedback gains kd14 and kd15 can be arbitrarily in the sets
Π d 14 = { θ R | 0 < θ < ( γ 1 + γ 2 ) | α R | μ H G r μ R }
and
Π d 15 = { θ R | 0 < θ < [ γ 1 + Ω S 0 ( 1 + k p 25 ) Ω S 0 ( 1 + k p 25 ) + 2 M S 0 ( 1 + k p 35 ) γ 2 ] | α R | μ S G r μ R }
respectively. If constants γi and ηi (i = 1, 2) are given, then the values of feedback gains kp14, kp15, kd14 and kd15 are given. From Theorem 1, we can see that the globally asymptotic closed-loop stability can be guaranteed if inequalities Equations (23)–(25) are all well satisfied. Therefore, from Equations (19)–(22) as well as inequalities Equations (23), (29) and (30), the values of feedback gains kd11, kd14, kd15, kp11, kp14, kp15, κp24, κp25, κp35 are all not directly determined by the NSSS physical and thermal parameters. Actually, the values of feedback gains kp11, kp14, kp15 and κp35 can be arbitrarily given, and the values of feedback gains kd11, kd14, kd15, κp24 and κp25 can be chosen in the scopes determined by Equations (33), (34), (23), (29) and (30). Moreover, from the expressions of this newly-built coordinated control law composed of Equations (16), (31) and (32), there is no expressions about the dynamic model in this controller utilized for either state-prediction or control input generation. From the above discussion, MHTGR-based NSSS coordinated control given by Equations (16), (31) and (32) is model-free, which leads to easy implementation and tuning.

4. Simulation Results with Discussion

To show the feasibility and performance of the newly-built nonlinear adaptive coordinated control strategy given in Theorem 1, it is applied to the regulation of an NSSS of HTR-PM plant in this section.
The dynamic model of the MHTGR that is used for numerical simulation is the nodal model given in [22]. The OTSG model adopted the moving boundary model given in [23]. The code was developed based on VC++ [24], and this simulation was done by using this software. Actually, both the space-dependent reactor models given in [22,23] and point model Equation (1) can describe the NSSS dynamics. The former one gives a more detailed description and is more suitable for verification, and the later model is much simpler and is more suitable for control law design. The dynamic model of the second-loop system gives the feedwater temperature and pressure for the dynamic model of the NSSS. In the numerical verification, the simulation software adopts the software given in [24] which contains the dynamics models of the NSSS as well as the secondary-loop. Moreover, from the above Theorem 1 and Comment 1, the NSSS coordinated control law can be written as Equations (16)–(18), and can also take the form as Equations (16), (31) and (32). If the first form is adopted, then control gains kp24 and kp35 should satisfy inequalities Equations (24) and (25). Moreover, if the second form is adopted, then control gains κp24 and κp35 should satisfy inequalities Equations (29) and (30). Here, the second form is adopted, and control gains are chosen as kd11 = 2.0, kp11 = 0.2, kd14 = 0.5, kp14 = 0.05, kd15 = 0.02, kp15 = 0.001, κp24 = 0.5, κp25 = 0.05 and κp35 = 0.3 so that Equations (23), (29), (30), (33) and (34) are satisfied.
To illustrate the coordinate control performance, the following two case studies are done:
  • Case A: Power-demand decreases from 100% to 50% RFP (Reactor Full Power) linearly in 10 min.
  • Case B: Power-demand increases from 50% to 100% RFP linearly in 10 min.
Both cases A and B represent hard operations of the NHR. Dynamic responses of the relative nuclear power nr, reactor outlet helium temperature Tcout and OTSG outlet steam temperature Ts as well as the control rod speed vr, helium flowrate Gp and feedwater flowrate Gs generated by the newly-built coordinated controller in the above two cases are shown in Figure 2 and Figure 3, respectively. Here, the power-demand decrease starts at the 1000 s, and the power-demand increase starts at the 2000 s. The difference in the start-time of power-demand variation does not affect the verification result. Since the OTSG outlet steam temperature is constant (571 °C), and since the feed-water temperature of the OTSG is different at different power-levels, it is very clear that the feedwater flowrate in a 50% reactor full power-level (RFP) is not the 50% flowrate at 100% RFP.
Since the power-demand changes, the errors between the actual and referenced values of the relative nuclear power, reactor outlet helium temperature and OTSG outlet steam temperature become larger. These three errors make coordinated control law to drive the control rods, primary helium blower and feedwater pump to suppress these error signals. The system enters into a steady state when the supplied and demanded power balances with each other. From Figure 2 and Figure 3, it is clear that closed-loop stability can be well guaranteed by properly choosing the feedback gains kd11, kd14, kd15, kp11, kp14, kp15, κp24, κp25 and κp35 so that feedback gains kd14 and kd15 are in the sets given by Equations (33) and (34), and inequalities Equations (23), (29) and (30) are satisfied. Thus, it is very clear that the numerical phenomena in this simulation fit well with the theoretic results. Finally, this newly-built model-free coordinated control for MHTGR-based NSSSs proposed in this paper has a simple form and can be easily implemented in the practical engineering.
Due to the HTR-PM plant still being under construction, we cannot verify the result in this paper by some real operating data in this stage. However, the obtained control parameters, i.e., controller structures given by Equations (16), (31) and (32) as well as control gain scopes given by Equations (23), (29), (30), (33) and (34) have deep implications to the future realization of controlling the real HTR-PM nuclear system. Here, the implication lies in the following two aspects:
(1)
Since the coordinated control law given by Equations (16), (31) and (32) has been proven to have the globally asymptotically stabilization performance, the coordinated control law of the real HTR-PM NSSS module should take the same form as Equations (16), (31) and (32) in the engineering.
(2)
The control gains in the practical engineering should be chosen to satisfy conditions Equations (23), (29), (30), (33) and (34) generally. However, due to the mismatch between the real and designed NSSS parameters as well as that between the dynamic features of the real and designed actuators (such as the control rods, helium bowers and feedwater pumps), some experiments or tests should be given practically to verify that the overshoots and transition periods of the real response of key process variables meets the engineering requirements.
Figure 2. Simulation results in Case A: (a) Relative nuclear power; (b) Outlet helium temperature; (c) steam temperature; (d) Control rod speed signal; (e) Primary helium flowrate; (f) Feedwater flowrate.
Figure 2. Simulation results in Case A: (a) Relative nuclear power; (b) Outlet helium temperature; (c) steam temperature; (d) Control rod speed signal; (e) Primary helium flowrate; (f) Feedwater flowrate.
Energies 09 00037 g002
Furthermore, the parameters of the NSSS models for control design and simulation are given by physical and thermal-hydraulic design. It was verified that the model and software given in [22,23,24] can provide credible dynamical simulation for the HTR-PM plant. The theoretical and simulation results in this paper show at least that the coordinated control law determined by Equations (16), (31) and (32) can be implemented in the practical engineering. Verifying the results in this paper by the operating data of the HTR-PM plant will be very meaningful work in the future.
Figure 3. Simulation results in Case B: (a) Relative nuclear power; (b) Outlet helium temperature; (c) steam temperature; (d) Control rod speed signal; (e) Primary helium flowrate; and (f) Feedwater flowrate.
Figure 3. Simulation results in Case B: (a) Relative nuclear power; (b) Outlet helium temperature; (c) steam temperature; (d) Control rod speed signal; (e) Primary helium flowrate; and (f) Feedwater flowrate.
Energies 09 00037 g003

5. Conclusions

Due to the simpler, standardized, and safer modular design, SMRs have the virtues of being factory built, requiring smaller initial capital investment, and having shorter construction time. Because of its persistency of power supply, SMRs can be incorporated with new energy resources to build micro-grids and can also be used to build large-scale nuclear plants by the multi-modular scheme. Coordinated control of SMR-based NSSSs is meaningful to realize the load-following function of SMR-based nuclear plants. However, until now, there has still been very little work in the field of NSSS coordinated control. MHTGR is a typical SMR having the inherent safety features and economic competitive power, and has also been seen as a candidate for a next generation nuclear reactor. Coordinated control of MHTGR-based NSSSs is quite important in providing the stable and efficient operation for the MHTGR-based nuclear energy systems. In this paper, by considering the exterior disturbances, a model-free coordinated control law for the MHTGR-based NSSS is proposed, which does not contain any expressions of the NSSS dynamic model and is insensitive to the physical or thermal parameters. The sufficient conditions for guaranteeing globally asymptotic closed-loop stability are given. Numerical simulation results show the satisfactory regulation performance of this newly-built model-free coordinated controller. Due to its simplicity, this control law can be very easily implemented in those digital control system platforms.
Moreover, the MHTGR-based NSSS is quite different from the PWR-based NSSS. As stated in Section 1 of this paper, the MHTGR-based NSSS consists of a modular high temperature gas-cooled reactor (MHTGR), a helical-coil once-through steam-generator (OTSG), primary helium blower and some connecting pipes. However, the PWR-based NSSS usually consists of a pressurized water reactor (PWR), a U-tube steam-generator (UTSG), primary feed-water pump and some pipes. The main difference of the dynamics of the MHTGR and that of the PWR is determined by the reactivity effect of the coolant. For the MHTGR, the coolant is helium, which is transparent to the fission reaction and has NO reactivity feedback effect on the neutron kinetics. For the PWR, the coolant is the pressurized water, which is also utilized as the moderator and does have reactivity feedback effect on the neutron kinetics. The difference in the reactor dynamics is only the first difference between the MHTGR and PWR based NSSSs. The second difference is that the MHTGR-based NSSS adopts the helical-coil OTSG, and the PWR-based NSSS adopts the UTSG. Since natural circulation in the UTSG exists, its dynamic behavior is quite different from that of the OTSG. Therefore, the dynamic behavior of the PWR-based NSSS cannot be described by model Equation (1) in this paper, and the coordinated control proposed in this paper for the MHTGR-based NSSS cannot be directly applied to the other types of NSSSs such as the PWR-based NSSS. However, the key idea of control design approach in this paper, i.e., the physics-based control design approach is helpful to realize the coordinated control laws for other types of NSSSs.

Acknowledgments

The author would like to thank the anonymous reviewers for the constructive comments. Moreover, the author would like to gratefully acknowledge the contributions of Xiaojin Huang for the valuable discussions with him and the positive suggestions from him.

Conflicts of Interest

The author declares no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

Nomenclature

crrelative concentration of the delayed neutron precursor
nrrelative neutron power
vrdesigned control rod speed signal
Grdifferential worth of a control rod
MSmass flowrate times the heat capacity of the secondary water/steam flow
P0rated thermal power
THaverage helium temperature of the primary loop
TRaverage temperature of the reactor core constituted by the fuel elements and reflector
TR,minitial equilibrium value of TR
TSaverage coolant temperature of the OTSG secondary side
TSininlet coolant temperature of the OTSG secondary side
αRreactivity feedback coefficient of the reactor core temperature
βfraction of the delayed neutrons
λeffective radioactive decay constant of the delayed neutron precursor
ρrreactivity given by the control rods
μRtotal heat capacity of the reactor core
μHtotal heat capacity of the helium inside the primary loop
μStotal heat capacity of the water/steam inside the OTSG secondary side
Λeffective prompt neutron life time
ΩPheat transfer coefficient between the helium and the reactor core
ΩSheat transfer coefficient between the two sides of the OTSG

Appendix: Proof of Theorem 1

The proof the Theorem 1 presented in Section 3 is given in the following two steps.
First, according to the backstepping technique, control inputs u2, u3, and virtual control input ξr to stabilize subsystem:
{ x ˙ = f ( x ) + g 1 ( x ) ξ r + g 2 ( x ) u 2 + g 3 ( x ) u 3 , y = h ( x ) ,
where ξr is called the reference signal of ξ.
By substituting u2 given by Equation (17) and u3 given by Equation (18) to subsystem, we have
{ Λ x ˙ 1 = β ( x 1 x 2 ) + ( n r 0 + x 1 ) ( α R x 3 + G r ξ ) , x ˙ 2 = λ ( x 1 x 2 ) , μ R x ˙ 3 = Ω P ( x 3 x 4 ) + P 0 x 1 , μ H x ˙ 4 = Ω P ( x 3 x 4 ) Ω S 0 [ ( 1 + k p 24 ) x 4 ( 1 + k p 25 ) x 5 ] , μ S x ˙ 5 = Ω S 0 [ ( 1 + k p 24 ) x 4 ( 1 + k p 25 ) x 5 ] 2 M S 0 ( 1 + k p 35 ) x 5 ,
From Equation (A2), it is clear that u2 is utilized to reinforce the heat transfer between the primary helium flow and secondary water/steam flow, and u3 is used to strengthen the dissipation term in the OTSG dynamics so that the secondary loop of the OTSG can be more stable.
Construct the Lyapunov function of subsystem:
V 1 ( x ) = n r 0 ζ N ( x 1 , x 2 ) + | α R | P 0 [ ζ T ( x 3 , x 4 , x 5 ) + γ 1 ς T 1 ( x 3 , x 4 , x 5 ) + γ 2 ς T 2 ( x 3 , x 4 , x 5 ) ]
where
ζ N ( x 1 , x 2 ) = Λ [ ( 1 + x 1 n r 0 ) ln ( 1 + x 1 n r 0 ) ] + β λ [ ( 1 + x 2 n r 0 ) ln ( 1 + x 2 n r 0 ) ]
and
ζ T ( x 3 , x 4 , x 5 ) = 1 2 ( μ R x 3 2 + μ H x 4 2 + 1 + k p 25 1 + k p 24 μ S x 5 2 )
are the shifted-ectropies of the neutron kinetics and thermal-hydraulics respectively, and both
ς T 1 ( x 3 , x 4 , x 5 ) = 1 2 μ R [ μ R x 3 + μ H x 4 + μ S x 5 + 2 M S 0 ( 1 + k p 35 ) 0 t x 5 ( s ) d s ] 2
and
ς T 2 ( x 3 , x 4 , x 5 ) = 1 2 μ R [ μ R x 3 + μ H x 4 + Ω S 0 ( 1 + k p 25 ) μ S x 5 + 2 M S 0 Ω S 0 ( 1 + k p 24 ) ( 1 + k p 35 ) 0 t x 4 ( s ) d s Ω S 0 ( 1 + k p 25 ) + 2 M S 0 ( 1 + k p 35 ) ] 2
are performance index functions. From Eqautions (A4)–(A7), it is clear that function Eqaution (A3) is positive definite.
Differentiate function V1 defined by Eqaution (A3) along the trajectory determined by model Equation (A1),
V ˙ 1 ( x ) = n r 0 β ( x 1 x 2 ) 2 ( n r 0 + x 1 ) ( n r 0 + x 2 ) + x 1 ( G r ξ + α R x 3 ) + | α R | [ x 1 x 3 Ω P P 0 ( x 3 x 4 ) 2 Ω S 0 P 0 ( 1 + k p 24 ) ( x 4 1 + k p 25 1 + k p 24 x 5 ) 2 2 M S 0 P 0 ( 1 + k p 35 ) 1 + k p 25 1 + k p 24 x 5 2 ] + γ 1 | α R | x 1 [ x 3 + μ H μ R x 4 + μ S μ R x 5 + 2 M S 0 μ R ( 1 + k p 35 ) 0 t x 5 ( s ) d s ] + γ 2 | α R | x 1 [ x 3 + μ H μ R x 4 + Ω S 0 μ R μ S ( 1 + k p 25 ) x 5 + 2 M S 0 ( 1 + k p 35 ) ( 1 + k p 24 ) 0 t x 4 ( s ) d s Ω S 0 ( 1 + k p 25 ) + 2 M S 0 ( 1 + k p 35 ) ] n r 0 β ( x 1 x 2 ) 2 ( n r 0 + x 1 ) ( n r 0 + x 2 ) + G r x 1 ξ | α R | P 0 { Ω P 2 x 3 2 + [ Ω S 0 2 ( 1 + k p 24 ) Ω P ] x 4 2 + 1 + k p 25 1 + k p 24 [ 2 M S 0 ( 1 + k p 35 ) Ω S 0 ( 1 + k p 25 ) ] x 5 2 } + | α R | ( γ 1 + γ 2 ) ( x 1 x 3 + μ H μ R x 1 x 4 ) + | α R | [ γ 1 + Ω S 0 ( 1 + k p 25 ) Ω S 0 ( 1 + k p 25 ) + 2 M S 0 ( 1 + k p 35 ) γ 2 ] μ S μ R x 1 x 5 + | α R | γ 1 ( 1 + k p 35 ) 2 M S 0 μ R x 1 0 t x 5 ( s ) d s + | α R | γ 2 ( 1 + k p 24 ) 2 M S 0 ( 1 + k p 35 ) Ω S 0 ( 1 + k p 25 ) + 2 M S 0 ( 1 + k p 35 ) Ω S 0 μ R x 1 0 t x 4 ( s ) d s n r 0 β ( x 1 x 2 ) 2 ( n r 0 + x 1 ) ( n r 0 + x 2 ) | α R | P 0 { Ω P 4 x 3 2 + Ω S 0 2 [ k p 24 2 Ω P Ω S 0 ( 1 + η 1 μ H μ R ) + 1 ] x 4 2 + 2 M S 0 1 + k p 25 1 + k p 24 [ k p 35 Ω S 0 2 M S 0 ( 1 + k p 25 + η 2 Ω P μ S Ω S 0 μ R ) + 1 ] x 5 2 } + | α R | P 0 Ω P ( γ 1 + γ 2 ) 2 [ 1 + η 1 μ H 4 μ R + η 2 μ S ( 1 + k p 24 ) 4 μ R ( 1 + k p 25 ) ] x 1 2 + G r x 1 { ξ + ( 1 η 1 ) ( γ 1 + γ 2 ) | α R | μ H G r μ R x 4 + + ( 1 η 2 ) [ γ 1 + Ω S 0 ( 1 + k p 25 ) Ω S 0 ( 1 + k p 25 ) + 2 M S 0 ( 1 + k p 35 ) γ 2 ] | α R | μ S G r μ R x 5 + γ 1 ( 1 + k p 35 ) 2 | α R | M S 0 G r μ R 0 t x 5 ( s ) d s + γ 2 ( 1 + k p 24 ) 2 M S 0 ( 1 + k p 35 ) Ω S 0 ( 1 + k p 25 ) + 2 M S 0 ( 1 + k p 35 ) | α R | Ω S 0 G r μ R 0 t x 4 ( s ) d s } = n r 0 β ( x 1 x 2 ) 2 ( n r 0 + x 1 ) ( n r 0 + x 2 ) σ d 11 x 1 2 | α R | P 0 ( Ω P 4 x 3 2 + σ p 24 x 4 2 + σ p 35 x 5 2 ) + G r x 1 [ ξ + k d 11 x 1 + k d 14 x 4 + k p 14 0 t x 4 ( s ) d s + k d 15 x 5 + k p 15 0 t x 5 ( s ) d s ] ,
where
σ d 11 = G r k d 11 | α R | P 0 Ω P ( γ 1 + γ 2 ) 2 [ 1 + η 1 μ H 4 μ R + η 2 μ S ( 1 + k p 24 ) 4 μ R ( 1 + k p 25 ) ]
σ p 24 = Ω S 0 2 [ k p 24 2 Ω P Ω S 0 ( 1 + η 1 μ H μ R ) + 1 ]
and
σ p 35 = 2 M S 0 1 + k p 25 1 + k p 24 [ k p 35 Ω S 0 2 M S 0 ( 1 + k p 25 + η 2 Ω P μ S Ω S 0 μ R ) + 1 ]
Choose virtual control ξ as:
ξ = [ k d 11 x 1 + k d 14 x 4 + k p 14 0 t x 4 ( s ) d s + k d 15 x 5 + k p 15 0 t x 5 ( s ) d s ]
where feedback gains kd11, kd14, kd15, kp14 and kp15 are all positive scalars satisfying Equations (19)–(22) respectively.
By substituting virtual control law Equation (A12) for inequality Equation (A8), and by considering inequalities Equations (23)–(25), we have
V ˙ 1 ( x ) n r 0 β ( x 1 x 2 ) 2 ( n r 0 + x 1 ) ( n r 0 + x 2 ) σ d 11 x 1 2 | α R | P 0 ( Ω P 4 x 3 2 + σ p 24 x 4 2 + σ p 35 x 5 2 )
From inequalities Equations (23)–(25), it is clear that σd11, σp24 and σp35 are positive scalars. Thus, the closed-loop constituted by subsystem Equation (A1), virtual control Equation (A12) and control laws Equations (17) and (18) is certainly globally asymptotically stable. Moreover, from inequalities Equation (A13) and Equations (A9)–(A11), feedback gains kd11, kp24 and kp35 are higher, the closed-loop system is more stable.
Second, the control input u1 will be designed. Choose the Lyapunov function for the entire system Equation (9) as:
V 2 ( x , e ξ ) = V 1 ( x ) + G r 2 k p 11 e ξ 2
where
e ξ = ξ ξ r
Differentiate V2 along the trajectory given by entire system dynamics Equation (9),
V ˙ 2 ( x , e ξ ) = V ˙ 1 | e ξ 0 + G r k p 11 e ξ ( u 1 + k p 11 x 1 ξ ˙ r )
from which we can see that if control u1 is designed as:
u 1 = k p 11 x 1 + ξ ˙ r = ( k p 11 x 1 + k d 11 x ˙ 1 + k p 14 x 4 + k d 14 x ˙ 4 + k p 15 x 5 + k d 15 x ˙ 5 )
then
V ˙ 2 ( x , e ξ ) = V ˙ 1 | ξ ξ r n r 0 β ( x 1 x 2 ) 2 ( n r 0 + x 1 ) ( n r 0 + x 2 ) σ d 11 x 1 2 | α R | P 0 ( Ω P 4 x 3 2 + σ p 24 x 4 2 + σ p 35 x 5 2 )
From inequality Equation (A18), the entire system state z defined by
z = [ x T ξ ] T
and enters the set
Ξ = { z R 5 | x O }
Moreover, from the first differential equation of model Equation (A2), xO induces that ξ ≡ 0. Thus, if state z enters set Ξ, then zO finally, which means that the closed-loop system formed by system Equation (9) and control laws Equation (16)—(18) is globally asymptotic stability. This completes the proof of Theorem 1.

References

  1. Vujić, J.; Bergmann, R.M.; Škoda, R.; Miletić, M. Small modular reactors: Simpler, safer, cheaper? Energy 2012, 45, 288–295. [Google Scholar] [CrossRef]
  2. Rowinski, M.K.; White, T.J.; Zhao, J. Small and medium sized reactors (SMR): A review of technology. Renew. Sustain. Energy Rev. 2015, 44, 643–656. [Google Scholar] [CrossRef]
  3. Reutler, H.; Lohnert, G.H. The modular high-temperature reactor. Nucl. Technol. 1983, 62, 22–30. [Google Scholar]
  4. Reutler, H.; Lohnert, G.H. Advantages of going modular in HTRs. Nucl. Eng. Des. 1984, 78, 129–136. [Google Scholar] [CrossRef]
  5. Lohnert, G.H. Technical design features and essential safety-related properties of the HTR-MODULE. Nucl. Eng. Des. 1990, 121, 259–275. [Google Scholar] [CrossRef]
  6. Lanning, D.D. Modularized high-temperature gas-cooled reactor systems. Nucl. Technol. 1989, 88, 139–156. [Google Scholar]
  7. Wu, Z.; Lin, D.; Zhong, D. The design features of the HTR-10. Nucl. Eng. Des. 2002, 218, 25–32. [Google Scholar] [CrossRef]
  8. Zhang, Z.; Wu, Z.; Wang, D.; Xu, Y.; Sun, Y.; Li, F.; Dong, Y. Current status and technical description of Chinese 2 × 250 MWth HTR-PM demonstration plant. Nucl. Eng. Des. 2009, 239, 2265–2274. [Google Scholar] [CrossRef]
  9. Shtessel, Y.B. Sliding mode control of the space nuclear reactor system. IEEE Trans. Aerosp. Electron. Syst. 1998, 34, 579–589. [Google Scholar] [CrossRef]
  10. Huang, Z.; Edwards, R.M.; Lee, K.Y. Fuzzy-adapted recursive sliding-mode controller design for a nuclear power plant control. IEEE Trans. Nucl. Sci. 2004, 51, 256–266. [Google Scholar] [CrossRef]
  11. Dong, Z. Nonlinear state-feedback dissipation power level control for nuclear reactors. IEEE Trans. Nucl. Sci. 2011, 58, 241–257. [Google Scholar] [CrossRef]
  12. Dong, Z. Nonlinear dynamic output-feedback power-level control for PWRs: A shifted-ectropy based design approach. Prog. Nucl. Sci. 2013, 68, 223–234. [Google Scholar] [CrossRef]
  13. Dong, Z. PD power-level control design for PWRs: A physically-based approach. IEEE Trans. Nucl. Sci. 2013, 60, 3889–3898. [Google Scholar] [CrossRef]
  14. Dong, Z. Dynamic output feedback power-level control for the MHTGR based on iterative damping assignment. Energies 2012, 5, 1782–1815. [Google Scholar] [CrossRef]
  15. Dong, Z. Physically-based power-level control for modular high temperature gas-cooled reactors. IEEE Trans. Nucl. Sci. 2012, 59, 2531–2549. [Google Scholar] [CrossRef]
  16. Na, M.G.; Shim, Y.R.; Lee, Y.J. Design of an adaptive predictive controller for steam generator. IEEE Trans. Nucl. Sci. 2003, 50, 186–193. [Google Scholar]
  17. Munasinghe, S.R.; Kim, M.-S.; Lee, J.-J. Adaptive neurofuzzy controller to regulate UTSG water level in nuclear plants. IEEE Trans. Nucl. Sci. 2005, 52, 421–429. [Google Scholar] [CrossRef]
  18. Dong, Z.; Huang, X.; Zhang, L. Saturated output feedback dissipation steam temperature control for the OTSG of MHTGRs. IEEE Trans. Nucl. Sci. 2011, 58, 1277–1289. [Google Scholar] [CrossRef]
  19. Wang, L.-X. Stable adaptive fuzzy control of nonlinear systems. IEEE Trans. Fuzzy Syst. 1993, 1, 146–155. [Google Scholar] [CrossRef]
  20. Marseguerra, M.; Zio, E. Model-free fuzzy tracking control of a nuclear reactor. Ann. Nucl. Energy 2003, 30, 953–981. [Google Scholar] [CrossRef]
  21. Rigatos, G.G. Model-based and model-free control of flexible-link robots: A comparison between representative methods. Appl. Math. Model. 2009, 33, 3906–3925. [Google Scholar] [CrossRef]
  22. Dong, Z.; Huang, X.; Zhang, L. A nodal dynamic model for control system design and simulation of an MHTGR core. Nucl. Eng. Des. 2010, 240, 1251–1261. [Google Scholar] [CrossRef]
  23. Dong, Z. A Differential-algebraic model for the once-through steam generator of MHTGR-based multi- modular nuclear plants. Math. Probl. Eng. 2015, 2015, 370101. [Google Scholar] [CrossRef]
  24. Dong, Z.; Huang, X. Real-time simulation platform for the design and verification of the operation strategy of the HTR-PM. In Proceedings of the 21th International Conference on Nuclear Engineering, Chengdu, China, 29 July–2 August 2013.

Share and Cite

MDPI and ACS Style

Dong, Z. Model-Free Coordinated Control for MHTGR-Based Nuclear Steam Supply Systems. Energies 2016, 9, 37. https://0-doi-org.brum.beds.ac.uk/10.3390/en9010037

AMA Style

Dong Z. Model-Free Coordinated Control for MHTGR-Based Nuclear Steam Supply Systems. Energies. 2016; 9(1):37. https://0-doi-org.brum.beds.ac.uk/10.3390/en9010037

Chicago/Turabian Style

Dong, Zhe. 2016. "Model-Free Coordinated Control for MHTGR-Based Nuclear Steam Supply Systems" Energies 9, no. 1: 37. https://0-doi-org.brum.beds.ac.uk/10.3390/en9010037

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop