# Dynamic Modeling and Control Characteristics of the Two-Modular HTR-PM Nuclear Plant

Authors: Zhe Dong; Yifei Pan; Maoxuan Song; Xiaojing Huang; Yujie Dong; Zuoyi Zhang
Journal: Science and Technology of Nuclear Installations
 (2017)

---

## Abstract

The modular high temperature gas-cooled reactor (MHTGR) is a typical small modular reactor (SMR) with inherent safety feature. Due to its high reactor outlet coolant temperature, the MHTGR can be applied not only for electricity production but also as a heat source for industrial complexes. Through multimodular scheme, that is, the superheated steam flows produced by multiple MHTGR-based nuclear supplying system (NSSS) modules combined together to drive a common thermal load, the inherent safety feature of MHTGR is applicable to large-scale nuclear plants at any desired power ratings. Since the plant power control technique of traditional single-modular nuclear plants cannot be directly applied to the multimodular plants, it is necessary to develop the power control method of multimodular plants, where dynamical modeling, control design, and performance verification are three main aspects of developing plant control method. In this paper, the study in the power control for two-modular HTR-PM plant is summarized, and the verification results based on numerical simulation are given. The simulation results in the cases of plant power step and ramp show that the plant control characteristics are satisfactory.


---

## Body

## 1. Introduction

With comparison to burning fossil fuels, there is nearly no greenhouse gas emission in nuclear fission reaction. Thus, nuclear fission energy is a crucial clean energy for giving basis products such as electricity, fresh water, and process heat and for addressing the challenges associated with global climate and environmental impact [1, 2]. Actually, nuclear energy may substitute the fossil in a centralized way and in a great amount with commercial availability and economic competitiveness [3]. After the successful development of small (tens of megawatts) light water reactors (LWRs) for propulsion by the US Navy, the commercial nuclear reactors began to be commissioned since the late 1950s. These early-stage nuclear reactors were the simply scaled up versions of the naval reactor with high power density [4]. Although the increase of the reactor size induced the economic competitiveness with those fossil plants, the complicated operational issues began to moderate the industry’s confidence in the plant safety. More stringent safety requirements were then imposed, which induced a complex layering of redundant safety and auxiliary systems to the original simple LWRs. The escalation of plant complexity contributed to the rapidly increasing costs, licensing periods, and construction delays. However, the plant safety still cannot be guaranteed satisfactorily. After the severe nuclear accidents, that is, Three Mile Island, Chernobyl, and Fukushima, the safety issues of nuclear reactors have become much more significant than before. Thus, there should be a novel reactor technology providing a satisfactory tradeoff between safety and economics.Small modular reactors (SMRs) are those nuclear fission reactors whose electrical output power is less than 300MWe [4–6]. Due to the low power density and large heat capacity, some SMRs even have the inherent safety feature which prevents SMRs from the hazards of core-melting, radiological release, and LOCA (Loss of Coolant Accident). With load-following function, SMRs can be incorporated with renewable energy sources to build hybrid energy systems (HESs) having the virtues such as persistent power supply and free refueling [1, 2]. By adopting multimodular operation strategy [7], that is, multiple nuclear steam supplying system (NSSS) modules feeding steam to common thermal load equipment such as the turbine/generator set, the inherent safety feature of a SMR can be applicable to large-scale nuclear plants at any desired power ratings, which could be beneficial in providing electricity to remote 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 the industrial complexes. Furthermore, the SMRs can offer simpler, safer, and standardized modular design by being factory-built, requiring smaller initial capital investment, and having shorter construction period, and have been regarded by International Atomic Energy Agency (IAEA) as one of the active developing trends of nuclear energy.Due to its inherent safety feature, the modular high temperature gas-cooled reactor (MHTGR, such as HTR-Module designed in Germany and MHTGR designed in US) has been seen as one of the best candidates for the next-generation nuclear plants. The MHTGR is a typical SMR which uses helium as coolant and graphite as moderator and structural material, and its inherent safety is determined by the low power density, strong negative temperature feedback effect, and slim reactor shape [8–10]. China began to study the MHTGR at the end of 1970s, and a 10 MWth pebble-bed high temperature gas-cooled reactor HTR-10, which was designed by the Institute of Nuclear and New Energy Technology (INET) of Tsinghua University, achieved its criticality in December 2000 and full power in January 2003 [11]. Then, six safety demonstration tests were done on the HTR-10, which manifested its inherent safety and self-stabilizing features [12]. Based on the experience of the HTR-10 project, a high temperature gas-cooled reactor pebble-bed module (HTR-PM) plant was then proposed, which consists of two pebble-bed one-zone MHTGRs with combined 2 × 250 MWth power and adopts the operation scheme of two NSSS modules driving one steam turbine [13, 14]. Figure 1 shows that a NSSS module of the HTR-PM plant is composed of an MHTGR, a helically coiled once-through steam generator (OTSG), a primary helium blower, and some connecting pipes. Figure 2 illustrates a simplified flow diagram of the entire two-modular HTR-PM plant. The HTR-PM plant is now in the stage of plant construction and commissioning. By incorporating more MHTGR-based NSSS modules, more large-scale multimodular high temperature gas-cooled nuclear plant can be realized.Figure 1
Composition of the MHTGR-based NSSS module.Figure 2
Simplified flow diagram for the HTR-PM plant.Satisfactory plant dynamic behavior is important for realizing load-following function and for further forming HESs with fossil or renewable energy sources, which is given by not only the open-loop dynamics but also the control strategy. There have been some promising results in dynamical modeling and control design for entire nuclear plants. In the field of plant modeling, Han gave a mathematical model of pressurized water reactors for thermal-hydraulic analysis [15]. Fazekas et al. proposed a simple dynamic model for the primary circuit in VVER plants for controller design purposes [16]. Dong et al. presented a lumped parameter dynamic model for control design and verification for the NSSS of the nuclear heating reactors (NHRs) [17]. With comparison to the abundant results in reactor power-level regulation such as the linear control, nonlinear control, and intelligent method [18–20], there are still less results in the field of entire nuclear plant control. Shtessel proposed a sliding-mode plant control and state-observation strategy for space nuclear energy system TOPAZ II [21]. Huang et al. gave a multi-input-multioutput (MIMO) fuzzy-adapted recursive sliding-mode plant controller for an advanced boiling water reactor (ABWR) nuclear power plant [22]. In order to achieve safe, stable, and efficient operation of the HTR-PM plant, it is necessary to give the lumped parameter model of the main plant dynamics and to propose proper plant control strategy. It is worth noting that, due to the two-modular scheme, there is much difference between the modeling and control of the HTR-PM plant and that of those single-modular nuclear plant.In this paper, both the dynamical modeling and control design method for HTR-PM plant are first summarized or reviewed, the realization of feedback loops and control laws are proposed, and then the verification results based upon numerical simulation are given, which shows the main transient performance of the HTR-PM plant is satisfactory for building HESs with the renewables.


## 2. Dynamical Modeling

In this section, the main parts of dynamic model of the HTR-PM plant, that is, the pebble-bed one-zone MHTGR, OTSG, fluid flow network (FFN) of the common secondary-loop system, and turbine/generator unit, are summarized from [23–26]. This model is proposed for control system design and verification, whose parameters are determined by the data collected from physical and thermal-hydraulic design.


### 2.1. MHTGR

From Figure1, the cold helium enters the primary blower where it is pressurized before flowing into the cold gas duct. Then, the cold helium gets into the channels (called riser) inside the side reflector from bottom to top and passes through the pebble-bed (called downcomer) from top to bottom where it is heated to a high temperature about 750°C. The hot helium leaves the hot gas chamber inside the bottom reflector and flows into the OTSG primary side where it transfers the heat to the secondary feedwater. To guarantee the inherent safety feature, the MHTGR of the HTR-PM plant has a slim core shape. Relative to using a single node [23], it is more proper to divide the reactor core intoN nodes along the axial direction [24]. The nodalization is shown in Figure 3, where the riser is some tubes in the side reflector through which helium can flow upward to the upper entry of the pebble-bed, and the downcomer is the volume gap among the fuel elements inside the pebble-bed through which helium can flow down. Moreover, the pressure loss of the primary helium flow is very small with comparison to the primary pressure and is omitted here.Figure 3
Nodalization of the MHTGR.


#### 2.1.1. Neutron Kinetics

The nodalized neutron kinetics can be written as [24](1)dnr,1dt=ρ1-β-α1,1Λ1nr,1+1Λ1α1,2nr,2+∑k=16βkΛ1Cr,1,k,dnr,idt=ρi-β-αi,iΛinr,i+1Λiαi,i-1nr,i-1+αi,i+1nr,i+1+∑k=16βkΛiCr,i,k,i=2,…,N-1,dnr,Ndt=ρN-β-αN,NΛNnr,N+1ΛNαN,N-1nr,N-1+∑k=16βkΛNCr,N,k,dCr,i,kdt=λknr,i-Cr,i,k,i=1,…,N,k=1,2,…,6,where i=1,…,N, nr,i is the relative neutron flux of the ith node, Cr,i,k is the relative concentration of the kth delayed neutron group of the ith node, αi,i and αi,i (j=1,…,Ni) are the coupling coefficients corresponding to the ith node, Ni is the number of adjacent nodes of the ith node, Λi is the prompt neutron lifetime in the ith node, ρi is the reactivity of the ith node, β is the fraction of all the delayed neutrons in all the nodes, and βk is the fraction of thekth delayed neutron group in every node.The nodal reactivity feedback is expressed as function of the nodal fuel temperatureTc,i (i=1,…,N) and temperature of the reflector Tr. The total nodal reactivity is given as [24] (2)ρi=ρr,i+αf,i+αm,iTc,i-Tc0,i+αrTr-Tr0,where ρr,i is the nodal reactivity provided by the control rods, αf,i and αm,i are, respectively, the reactivity coefficients of the fuel and moderator inside nodei, αr is the reactivity coefficient of the average reflector temperature, and Tc0,i and Tr0 are, respectively, the initial temperature of the fuel inside node i and that of the reflector. Since the reflector is very thick, and since the metal core support structure is cooled by the cold helium with its temperature no more than 250°C, the reactivity feedback induced by the support structure is omitted here, which has been verified by the MHTGR physical and thermal-hydraulic design.Moreover, it can be derived that [24] (3)αi,j=nj0ni0DAυΣf,iVd,j=2,i=1i-1,i+1,i=2,…,N-1N-1,i=N,αi,i=DAυΣf,2Vd,i=1,DAυVd1Σf,i-1+1Σf,i+1,i=2,…,N-1,DAυΣf,N-1Vd,i=N,whereA is the cross-sectional area of the active core region, V is the volume of every node, d is the height of every node, υ is the fission number, Σa,i and Σf,i (i=1,…,N) are, respectively, the macroabsorption and fission cross sections in the ith node.


#### 2.1.2. Thermal Hydraulics

(
1) Pebble-Bed. The thermal dynamics of ith node of the pebble-bed can be written as [24] (4)1-ερc,iVCc,idTc,idt=P0,inr,i-Kd,iAd,iTc,i-Td,i-Kcr,iAcr,iTc,i-Tr+-Kc,1ATc,1-Tc,2,i=1,-Kc,iATc,i-Tc,i+1+Kc,i-1ATc,i-1-Tc,i,i=2,…,N-1,Kc,N-1ATc,N-1-Tc,N,i=N,where ρc,i, Cc,i, Tc,i, and  P0,i are, respectively, the density, specific heat, temperature, and rated power of fuel pile inside node i, Kd,i and Ad,i are, respectively, the heat transfer coefficient and heat transfer area between the fuel pile and helium inside node i, Td,i is the temperature of the helium inside nodei, Kcr,i and Acr,i are, respectively, the heat transfer coefficient and heat transfer area between the fuel pile in node i and the reflector, Tr is the temperature of the reflector, Kc,i (i=1,…,N-1) is the heat transfer coefficient between the fuel pile inside node i and that inside node i+1, and ε is the porosity of the pebble-bed.(
2) Reflector. To improve the heat capacity and to weaken the neutron leakage, the reflector of the MHTGR is very thick, which means that the thermal inertial of the reflector is so large that the dynamic response of the reflector is in different scale with that of the fuel elements and helium flow. Moreover, the side reflector is cooled by the cold helium flow injected into the reactor with a high speed. Therefore, for the simplicity of dynamic model, it is reasonable to regard the reflector as one nodal.From Figure3, there are two heat transfer processes with the reflector. This first one is the fission heat conducted from the fuel pile inside every node, which heats up the reflector. The second one is that the cold helium inside the riser cools down the reflector. The heat transfer about the reflector can be written as [24](5)ρrVrCrdTrdt=∑i=1NKcr,iAcr,iTc,i-Tr-KurAurTr-Tu,where ρr, Vr, and Cr are, respectively, the density, volume, and specific heat of the reflector, Kur and Aur are the heat transfer coefficient and heat transfer area between the reflector and the coolant in the riser, and Tu is the average coolant temperature inside the riser.(
3) Downcomer. The cold helium enters the downcomer from its upper entry and flows downward to cool the pebble-bed. The heat is mainly transferred to the helium inside the downcomer, which can be described as [24](6)Vidρd,idt=Wu-Wd,1,i=1,Wd,i-1-Wd,i,i=2,…,N,εViCpdρd,iTd,idt=Kd,iAd,iTc,i-Td,i+CpWuTu-Wd,1Td,1,i=1,CpWd,i-1Td,i-1-Wd,iTd,i,i=2,…,N,where ρd,i and Wd,i are, respectively, the density and flowrate of the helium inside node i of the downcomer and Cp is the helium specific heat at constant pressure. Then, substitute (6), the thermodynamics of the helium in the downcomer can be written as [24](7)εViρd,iCpdTd,idt=Kd,iAd,iTc,i-Td,i+CpWuTu-Td,1,i=1,CpWd,i-1Td,i-1-Td,i,i=2,…,N.(
4) Riser. In the channels of the riser, the cold helium flows upward to cool the reflector, and the heat transfer process can be described by [23, 24](8)Vudρudt=1-κWlh-Wu,VuCpdρuTudt=Cp1-κWlhTlh-WuTu+KurAurTr-Tu,where ρu, Wu, and Tu are, respectively, the density, flowrate, and temperature of the helium inside the riser, Wlh and Tlh are, respectively, the flowrate and temperature of the helium inside the lower header, and κ is the leakage ration of the lower header. Substitute (8), the thermodynamics of the helium inside the riser can be written as(9)VuρuCpdTudt=Cp1-κWlhTlh-Tu+KurAurTr-Tu.(
5) Lower Header. The lower header is assumed to be a well-mixed chamber which has no heat exchange with the environment, and the thermodynamics corresponding to the lower header can be expressed as [23, 24](10)VlhρlhdTlhdt=WlhTin-Tlh,where Tin is the temperature of the inlet cold helium, ρlh, Wlh, and Tlh are, respectively, the density, flowrate, and temperature of the helium inside the lower header, and Vlh is the volume of the low header.(
6) Outlet Header. The output header is assumed to be a well-mixed chamber, and its dynamics can be given as [23, 24](11)VohρohdTohdt=κWlhTlh-Toh+Wd,NTd,10-Toh,where ρoh, Wlh, and Toh are, respectively, the density, flowrate, and temperature of the helium inside the outlet header and Voh is the volume of the outlet header.


#### 2.1.3. MHTGR Dynamics

The dynamic model of MHTGR is composed of neutron kinetics determined by (1)–(3) as well as thermal hydraulics given by (4), (5), (7), and (9)–(11). The dynamic behavior of this MHTGR model can be influenced by nodal numberN, and an improper nodal number leads to the fluctuation of the responses. Actually, this model is only utilized for control design and verification and is not for physical and thermal-hydraulic design. The nodal numbers are determined by comparing the response of this model with that of the model for physical and thermal-hydraulic design. When the error between responses is satisfactorily small, the corresponding nodal number is suitable. Here, the nodal numberN is chosen to be 10. Based on the simulation results to be given in Section 4, it will be seen that there is no numerical fluctuation in the responses, and all the transience is determined by power-level maneuver which show that the nodal number is properly selected.Since the dynamic model in this paper is tuned based upon the comparison with the model for reactor physical and thermal-hydraulic design, the stability of this model is also guaranteed by the model for reactor design. Moreover, it was shown in [27] that the MHTGR is open-loop globally asymptotically stable, and by the use of the shifted-ectropy-based self-stability analysis method given in [27], it can be shown that the model coupled by neutron kinetics and thermal hydraulics are open-loop globally asymptotically stable.


### 2.2. OTSG

Since the OTSG can provide superheated steam which does not satisfy the one-to-one map from pressure to temperature, the OTSG is more proper for building multi-SMR-based nuclear plants. Therefore, the dynamic model of the OTSG is very crucial for studying behavior of the SMR-based plant and also for verifying the operation and control strategies. In this section, a moving-boundary model [25] is introduced for describing OTSG dynamics in the cases of generating saturated and superheated steam [26]. The nodalization of the OTSG is given in Figure 4.Figure 4
Nodalization scheme of the OTSG.


#### 2.2.1. Secondary Side

The dynamics of the OTSG secondary side is given by both the thermodynamics and pressure drops of the secondary fluid flow.(
1) Thermodynamics. The thermodynamics of the secondary side can be given by the following vector-valued nonlinear differential equation [26]:(12)Esxsdxsdt=fsxs,us,where (13)xs=l13hs5hs1Gs3T,us=Gs1hfwQs2Qs4T,Esxs=Es11xs012ρs2l13Es14xsEs21xs12ρs4l350Es24xs00τH0000τG,Es11xs=12ρs3hs1-hs3+Ps3-Ps1+l131-ρs2∂hs3∂Ps3Ps3ρs2gsin⁡θ+Fs2Gs22As22,Es14xs=Fs2Gs22l1322As221-ρs2∂hs3∂Ps3Ps3,Es21xs=12ρs3hs3-hs5+Ps5-Ps3+l352-ρs4∂hs3∂Ps3Ps3ρs2gsin⁡θ+Fs2Gs22As22-ρs4gsin⁡θ+Fs4Gs42As42,Es24xs=l351-12ρs4∂hs3∂Ps3Ps3Fs2Gs2l13As22+2Fs4Gs4l13As22,fsxs,us=Qs2As2+12Gs1As1+Gs3As3hs1-hs3Qs4+Gs4hs3-hs5As4hfw-hs1τHGs1-Gs3τGT,g is the gravitational acceleration (N·m/s2), both τG and τH are given positive constants denoting inertial effect, hfw is the enthalpy per unit mass (J/kg) of the feedwater flow, lij (i,j=1,…,5 and i<j) is the length between points i and j in Figure 4, and ρsi, Gsi, hsiQsi, and Asi are, respectively, the flow density (kg/m3), flowrate (kg/s), enthalpy per unit mass (J/kg), heat flux from the primary side, and cross section (m2) of node i of the secondary side. Here, l13 is the length of the subcooled section, l35 is that of the boiling section, and h1s is specific enthalpy of the feedwater. Due to the high speed of the steam flow, it is assumed that Gs3=Gs4=Gs5.(
2) Pressure Drops. The pressure drop of the fluid flow in the OTSG secondary side cannot be omitted, and the pressure of each node is given by [26](14)Ps3=Ps1-ρs2gl13sin⁡θ-Fs2l13Gs22,Ps5=Ps3-ρs4gl35sin⁡θ-Fs4l35Gs42,Ps4=Ps3+Ps52,Ps2=Ps1+Ps32,where θ is the helically coiled angle of the tubes, Fsi (i=2,4) is the friction factor of the two sections, and Ps1 is the feedwater pressure.


#### 2.2.2. Tube Wall

There is no flow inside the metal tube wall between the primary and secondary sides, and the dynamics of the tube wall are determined by its thermodynamics, which can be expressed as [26](15)dTm2dt=Qp2-Qs2CmρmAm2l13+Tm4-Tm2l15dl13dt,dTm4dt=Qp4-Qs4CmρmAm4l35+Tm4-Tm2l15dl13dt,where Qp2 and Qp4 are, respectively, the heat flux from the primary helium flow to tube wall in the subcooled and boiling sections.


#### 2.2.3. Primary Side

Since velocity of the primary helium flow is much higher than that of the two-phase flow inside the secondary side, the energy and temperature relationship of the primary side can be given by the following algebraic equation [26]:(16)Apxsxp=bpxm,xp,where(17)xp=Tp1Tp2Tp3Tp4T,Apxm=1-2100l35-l15l13-CpGp-Kp2l13CpGp000-CpGp-Kp4l35,bpxm,xp=00-Kp2l13Tm2-Kp4l35Tm4-CpGpTp5T,where Gp is the primary helium flowrate, Cp is the helium specific heat, Tpi is the helium temperature at point i (i=1,…,5), Kp2 and Kp4 are the heat transfer coefficients between the primary side and tube wall of the subcooled and boiling sections, respectively, and Ks2 and Ks4 are, respectively, the heat transfer coefficients between the tube wall and secondary side of the subcooled and boiling sections.


### 2.3. Secondary-Loop FFN

After the hot helium leaves the reactor, it flows into the OTSG primary side where it transfers the heat to the secondary feedwater. By absorbing the heat of the primary loop, the feedwater is converted to the superheated steam in 571°C and about 13.9 MPa or so. The main steam flow at 13.24 MPa combined from the two NSSS modules then enters the steam turbine for producing electricity. The condensed flow is fed into the OTSG secondary side again after it goes through the grand heater, low pressure heaters (LPHs), deaerator, and high pressure heaters (HPHs). Thus, the hydraulics of the secondary-loop system can be modeled as the FFN shown in Figure5, where branches f1 and f2 are the fan branches of which everyone has a feedwater pump, branches 1 and 2 denote the secondary side of the OTSGs corresponding to the first and second NSSS module respectively, and branch 3 contains the equipment of the plant secondary loop such as the steam turbine, condenser, low pressure heaters, and deaerator.Figure 5
Fluid flow network of the secondary loop of the HTR-PM plant.The flowrate dynamics of branches 1, 2, and 3 can be given as [28, 29] (18)dQdt+KQD2R=KH,where Q=Q1,Q2,Q3T, QD2=diag⁡Q2, K=diag⁡K1,K2,K3, R=R1,R2,R3T, H=H1,H2,H3T, and Qj, Rj, Hj, and Kj are the flowrate, resistance, pressure drop, and inertia coefficient of branch j (j=1,2,3), respectively.Every FFN can be divided into a tree containing the fan branches and its complement, that is, the cotree whose branches are referred to as the links. From Figure5, the branches with solid line, that is, branches 3, f1, and f2, constitute the tree of this FFN, and the branches with dashed line, that is, branches 1 and 2, are the links which constitute the cotree. Based on the network property, the pressure drops and the flowrates of the tree branches are not independent and obey the following algebraic equations given by Kirchhoff’s voltage law (KVL) and Kirchhoff’s current law (KCL) [28]: (19)H1H2T=-H3+Hf1H3+Hf2T,Q3Qf1Qf2T=Q1+Q2Q1Q2T, where Hfj is the pressure drop of fan branch j (j=1,2) and can be further expressed as(20)Hfj=-Hdj+RfjQfj,j=1,2.Then, based on (18) and (19), the hydraulics of the FFN shown in Figure 4 can be given by differential-algebraic system [29]:(21)dQcdt=-KcRfQc+QcD2Rc+KcHd-H3η,K-H3=kcTHd-kfcTQc-kcTQcD2Rc+K3Q32R3,where Qc=Q1,Q2T, QcD2=diag⁡Q1,Q22, kc=K1,K2T, kfc=K1Rf1,K2Rf2T, K-=K1+K2+K3, Kc=diag⁡K1,K2, Rf=diag⁡Rf1,Rf2, and η=11T.


### 2.4. Model Integration

The above dynamic models of the MHTGR, OTSG, and secondary-loop FFN can be integrated together with other commonly utilized models such as the turbine and the synchronous generator to form the dynamic model for the whole plant. Based upon this model, a simulation code for the dynamic behavior and control characteristics of the HTR-PM plant is developed under MATLAB/Simulink platform, and the composition of this Simulink model is shown in Figure6. This model can incorporate the plant control laws proposed in the next section for giving the simulations of the operational and control characteristics of the HTR-PM plant, where the program modules named as “1# NSSS” and “2# NSSS” are used to simulate the dynamic behavior of the two NSSS modules, the program module named as “secondary-loop FFN” is used to simulate the thermal hydraulics for the secondary-loop FFN, the module “Secondary Thermal Feature” is used to simulate the thermodynamics of secondary-loop systems, and module “Generator” is used to simulate the electric generator dynamics. Furthermore, some important model parameters are given in Table 1.Table 1
Main parameters of the dynamic model.
Parameter
Name
Unit
Value
Λ
i (i=1,…,10)
Prompt neutron lifetime
s
7.66
e
-
4
β
Fraction of all the delayed neutrons
%
0.67
β1
Fraction of 1st group delayed neutrons
%
0.0256
β2
Fraction of 2nd group delayed neutrons
%
0.14
β3
Fraction of 3rd group delayed neutrons
%
0.13
β4
Fraction of 4th group delayed neutrons
%
0.27
β5
Fraction of 5th group delayed neutrons
%
0.086
β6
Fraction of 6th group delayed neutrons
%
0.017
λ1
Decay constant of 1st group precursor
s−1
0.0256
λ2
Decay constant of 2nd group precursor
s−1
0.14
λ3
Decay constant of 3rd group precursor
s−1
0.13
λ4
Decay constant of 4th group precursor
s−1
0.27
λ5
Decay constant of 5th group precursor
s−1
0.086
λ6
Decay constant of 6th group precursor
s−1
0.017
A=αij10×10
Coupling coefficients
3.5
7.0
1.8
7.1
5.5
2.2
6.9
3.9
3.1
7.0
3.3
3.7
7.0
2.9
4.2
7.0
2.8
4.5
7.0
2.6
4.7
7.0
2.4
5.0
7.0
2.2
5.6
3.5
×
1
0
-
3
α
f
,
i (i=1,…,10)
Reactivity coefficient of fuel temperature
Δ
K
/
K/°C
-
4.36
e
-
5
α
m
,
i (i=1,…,10)
Reactivity coefficient of moderator temperature
Δ
K
/
K/°C
-
0.94
e
-
5
αr
Reactivity coefficient of reflector temperature
Δ
K
/
K/°C
1.49
e
-
5
ε
Porosity of the pebble-bed
0.39Figure 6
Composition of the SIMULINK code.


## 2.1. MHTGR

From Figure1, the cold helium enters the primary blower where it is pressurized before flowing into the cold gas duct. Then, the cold helium gets into the channels (called riser) inside the side reflector from bottom to top and passes through the pebble-bed (called downcomer) from top to bottom where it is heated to a high temperature about 750°C. The hot helium leaves the hot gas chamber inside the bottom reflector and flows into the OTSG primary side where it transfers the heat to the secondary feedwater. To guarantee the inherent safety feature, the MHTGR of the HTR-PM plant has a slim core shape. Relative to using a single node [23], it is more proper to divide the reactor core intoN nodes along the axial direction [24]. The nodalization is shown in Figure 3, where the riser is some tubes in the side reflector through which helium can flow upward to the upper entry of the pebble-bed, and the downcomer is the volume gap among the fuel elements inside the pebble-bed through which helium can flow down. Moreover, the pressure loss of the primary helium flow is very small with comparison to the primary pressure and is omitted here.Figure 3
Nodalization of the MHTGR.


### 2.1.1. Neutron Kinetics

The nodalized neutron kinetics can be written as [24](1)dnr,1dt=ρ1-β-α1,1Λ1nr,1+1Λ1α1,2nr,2+∑k=16βkΛ1Cr,1,k,dnr,idt=ρi-β-αi,iΛinr,i+1Λiαi,i-1nr,i-1+αi,i+1nr,i+1+∑k=16βkΛiCr,i,k,i=2,…,N-1,dnr,Ndt=ρN-β-αN,NΛNnr,N+1ΛNαN,N-1nr,N-1+∑k=16βkΛNCr,N,k,dCr,i,kdt=λknr,i-Cr,i,k,i=1,…,N,k=1,2,…,6,where i=1,…,N, nr,i is the relative neutron flux of the ith node, Cr,i,k is the relative concentration of the kth delayed neutron group of the ith node, αi,i and αi,i (j=1,…,Ni) are the coupling coefficients corresponding to the ith node, Ni is the number of adjacent nodes of the ith node, Λi is the prompt neutron lifetime in the ith node, ρi is the reactivity of the ith node, β is the fraction of all the delayed neutrons in all the nodes, and βk is the fraction of thekth delayed neutron group in every node.The nodal reactivity feedback is expressed as function of the nodal fuel temperatureTc,i (i=1,…,N) and temperature of the reflector Tr. The total nodal reactivity is given as [24] (2)ρi=ρr,i+αf,i+αm,iTc,i-Tc0,i+αrTr-Tr0,where ρr,i is the nodal reactivity provided by the control rods, αf,i and αm,i are, respectively, the reactivity coefficients of the fuel and moderator inside nodei, αr is the reactivity coefficient of the average reflector temperature, and Tc0,i and Tr0 are, respectively, the initial temperature of the fuel inside node i and that of the reflector. Since the reflector is very thick, and since the metal core support structure is cooled by the cold helium with its temperature no more than 250°C, the reactivity feedback induced by the support structure is omitted here, which has been verified by the MHTGR physical and thermal-hydraulic design.Moreover, it can be derived that [24] (3)αi,j=nj0ni0DAυΣf,iVd,j=2,i=1i-1,i+1,i=2,…,N-1N-1,i=N,αi,i=DAυΣf,2Vd,i=1,DAυVd1Σf,i-1+1Σf,i+1,i=2,…,N-1,DAυΣf,N-1Vd,i=N,whereA is the cross-sectional area of the active core region, V is the volume of every node, d is the height of every node, υ is the fission number, Σa,i and Σf,i (i=1,…,N) are, respectively, the macroabsorption and fission cross sections in the ith node.


### 2.1.2. Thermal Hydraulics

(
1) Pebble-Bed. The thermal dynamics of ith node of the pebble-bed can be written as [24] (4)1-ερc,iVCc,idTc,idt=P0,inr,i-Kd,iAd,iTc,i-Td,i-Kcr,iAcr,iTc,i-Tr+-Kc,1ATc,1-Tc,2,i=1,-Kc,iATc,i-Tc,i+1+Kc,i-1ATc,i-1-Tc,i,i=2,…,N-1,Kc,N-1ATc,N-1-Tc,N,i=N,where ρc,i, Cc,i, Tc,i, and  P0,i are, respectively, the density, specific heat, temperature, and rated power of fuel pile inside node i, Kd,i and Ad,i are, respectively, the heat transfer coefficient and heat transfer area between the fuel pile and helium inside node i, Td,i is the temperature of the helium inside nodei, Kcr,i and Acr,i are, respectively, the heat transfer coefficient and heat transfer area between the fuel pile in node i and the reflector, Tr is the temperature of the reflector, Kc,i (i=1,…,N-1) is the heat transfer coefficient between the fuel pile inside node i and that inside node i+1, and ε is the porosity of the pebble-bed.(
2) Reflector. To improve the heat capacity and to weaken the neutron leakage, the reflector of the MHTGR is very thick, which means that the thermal inertial of the reflector is so large that the dynamic response of the reflector is in different scale with that of the fuel elements and helium flow. Moreover, the side reflector is cooled by the cold helium flow injected into the reactor with a high speed. Therefore, for the simplicity of dynamic model, it is reasonable to regard the reflector as one nodal.From Figure3, there are two heat transfer processes with the reflector. This first one is the fission heat conducted from the fuel pile inside every node, which heats up the reflector. The second one is that the cold helium inside the riser cools down the reflector. The heat transfer about the reflector can be written as [24](5)ρrVrCrdTrdt=∑i=1NKcr,iAcr,iTc,i-Tr-KurAurTr-Tu,where ρr, Vr, and Cr are, respectively, the density, volume, and specific heat of the reflector, Kur and Aur are the heat transfer coefficient and heat transfer area between the reflector and the coolant in the riser, and Tu is the average coolant temperature inside the riser.(
3) Downcomer. The cold helium enters the downcomer from its upper entry and flows downward to cool the pebble-bed. The heat is mainly transferred to the helium inside the downcomer, which can be described as [24](6)Vidρd,idt=Wu-Wd,1,i=1,Wd,i-1-Wd,i,i=2,…,N,εViCpdρd,iTd,idt=Kd,iAd,iTc,i-Td,i+CpWuTu-Wd,1Td,1,i=1,CpWd,i-1Td,i-1-Wd,iTd,i,i=2,…,N,where ρd,i and Wd,i are, respectively, the density and flowrate of the helium inside node i of the downcomer and Cp is the helium specific heat at constant pressure. Then, substitute (6), the thermodynamics of the helium in the downcomer can be written as [24](7)εViρd,iCpdTd,idt=Kd,iAd,iTc,i-Td,i+CpWuTu-Td,1,i=1,CpWd,i-1Td,i-1-Td,i,i=2,…,N.(
4) Riser. In the channels of the riser, the cold helium flows upward to cool the reflector, and the heat transfer process can be described by [23, 24](8)Vudρudt=1-κWlh-Wu,VuCpdρuTudt=Cp1-κWlhTlh-WuTu+KurAurTr-Tu,where ρu, Wu, and Tu are, respectively, the density, flowrate, and temperature of the helium inside the riser, Wlh and Tlh are, respectively, the flowrate and temperature of the helium inside the lower header, and κ is the leakage ration of the lower header. Substitute (8), the thermodynamics of the helium inside the riser can be written as(9)VuρuCpdTudt=Cp1-κWlhTlh-Tu+KurAurTr-Tu.(
5) Lower Header. The lower header is assumed to be a well-mixed chamber which has no heat exchange with the environment, and the thermodynamics corresponding to the lower header can be expressed as [23, 24](10)VlhρlhdTlhdt=WlhTin-Tlh,where Tin is the temperature of the inlet cold helium, ρlh, Wlh, and Tlh are, respectively, the density, flowrate, and temperature of the helium inside the lower header, and Vlh is the volume of the low header.(
6) Outlet Header. The output header is assumed to be a well-mixed chamber, and its dynamics can be given as [23, 24](11)VohρohdTohdt=κWlhTlh-Toh+Wd,NTd,10-Toh,where ρoh, Wlh, and Toh are, respectively, the density, flowrate, and temperature of the helium inside the outlet header and Voh is the volume of the outlet header.


### 2.1.3. MHTGR Dynamics

The dynamic model of MHTGR is composed of neutron kinetics determined by (1)–(3) as well as thermal hydraulics given by (4), (5), (7), and (9)–(11). The dynamic behavior of this MHTGR model can be influenced by nodal numberN, and an improper nodal number leads to the fluctuation of the responses. Actually, this model is only utilized for control design and verification and is not for physical and thermal-hydraulic design. The nodal numbers are determined by comparing the response of this model with that of the model for physical and thermal-hydraulic design. When the error between responses is satisfactorily small, the corresponding nodal number is suitable. Here, the nodal numberN is chosen to be 10. Based on the simulation results to be given in Section 4, it will be seen that there is no numerical fluctuation in the responses, and all the transience is determined by power-level maneuver which show that the nodal number is properly selected.Since the dynamic model in this paper is tuned based upon the comparison with the model for reactor physical and thermal-hydraulic design, the stability of this model is also guaranteed by the model for reactor design. Moreover, it was shown in [27] that the MHTGR is open-loop globally asymptotically stable, and by the use of the shifted-ectropy-based self-stability analysis method given in [27], it can be shown that the model coupled by neutron kinetics and thermal hydraulics are open-loop globally asymptotically stable.


## 2.1.1. Neutron Kinetics

The nodalized neutron kinetics can be written as [24](1)dnr,1dt=ρ1-β-α1,1Λ1nr,1+1Λ1α1,2nr,2+∑k=16βkΛ1Cr,1,k,dnr,idt=ρi-β-αi,iΛinr,i+1Λiαi,i-1nr,i-1+αi,i+1nr,i+1+∑k=16βkΛiCr,i,k,i=2,…,N-1,dnr,Ndt=ρN-β-αN,NΛNnr,N+1ΛNαN,N-1nr,N-1+∑k=16βkΛNCr,N,k,dCr,i,kdt=λknr,i-Cr,i,k,i=1,…,N,k=1,2,…,6,where i=1,…,N, nr,i is the relative neutron flux of the ith node, Cr,i,k is the relative concentration of the kth delayed neutron group of the ith node, αi,i and αi,i (j=1,…,Ni) are the coupling coefficients corresponding to the ith node, Ni is the number of adjacent nodes of the ith node, Λi is the prompt neutron lifetime in the ith node, ρi is the reactivity of the ith node, β is the fraction of all the delayed neutrons in all the nodes, and βk is the fraction of thekth delayed neutron group in every node.The nodal reactivity feedback is expressed as function of the nodal fuel temperatureTc,i (i=1,…,N) and temperature of the reflector Tr. The total nodal reactivity is given as [24] (2)ρi=ρr,i+αf,i+αm,iTc,i-Tc0,i+αrTr-Tr0,where ρr,i is the nodal reactivity provided by the control rods, αf,i and αm,i are, respectively, the reactivity coefficients of the fuel and moderator inside nodei, αr is the reactivity coefficient of the average reflector temperature, and Tc0,i and Tr0 are, respectively, the initial temperature of the fuel inside node i and that of the reflector. Since the reflector is very thick, and since the metal core support structure is cooled by the cold helium with its temperature no more than 250°C, the reactivity feedback induced by the support structure is omitted here, which has been verified by the MHTGR physical and thermal-hydraulic design.Moreover, it can be derived that [24] (3)αi,j=nj0ni0DAυΣf,iVd,j=2,i=1i-1,i+1,i=2,…,N-1N-1,i=N,αi,i=DAυΣf,2Vd,i=1,DAυVd1Σf,i-1+1Σf,i+1,i=2,…,N-1,DAυΣf,N-1Vd,i=N,whereA is the cross-sectional area of the active core region, V is the volume of every node, d is the height of every node, υ is the fission number, Σa,i and Σf,i (i=1,…,N) are, respectively, the macroabsorption and fission cross sections in the ith node.


## 2.1.2. Thermal Hydraulics

(
1) Pebble-Bed. The thermal dynamics of ith node of the pebble-bed can be written as [24] (4)1-ερc,iVCc,idTc,idt=P0,inr,i-Kd,iAd,iTc,i-Td,i-Kcr,iAcr,iTc,i-Tr+-Kc,1ATc,1-Tc,2,i=1,-Kc,iATc,i-Tc,i+1+Kc,i-1ATc,i-1-Tc,i,i=2,…,N-1,Kc,N-1ATc,N-1-Tc,N,i=N,where ρc,i, Cc,i, Tc,i, and  P0,i are, respectively, the density, specific heat, temperature, and rated power of fuel pile inside node i, Kd,i and Ad,i are, respectively, the heat transfer coefficient and heat transfer area between the fuel pile and helium inside node i, Td,i is the temperature of the helium inside nodei, Kcr,i and Acr,i are, respectively, the heat transfer coefficient and heat transfer area between the fuel pile in node i and the reflector, Tr is the temperature of the reflector, Kc,i (i=1,…,N-1) is the heat transfer coefficient between the fuel pile inside node i and that inside node i+1, and ε is the porosity of the pebble-bed.(
2) Reflector. To improve the heat capacity and to weaken the neutron leakage, the reflector of the MHTGR is very thick, which means that the thermal inertial of the reflector is so large that the dynamic response of the reflector is in different scale with that of the fuel elements and helium flow. Moreover, the side reflector is cooled by the cold helium flow injected into the reactor with a high speed. Therefore, for the simplicity of dynamic model, it is reasonable to regard the reflector as one nodal.From Figure3, there are two heat transfer processes with the reflector. This first one is the fission heat conducted from the fuel pile inside every node, which heats up the reflector. The second one is that the cold helium inside the riser cools down the reflector. The heat transfer about the reflector can be written as [24](5)ρrVrCrdTrdt=∑i=1NKcr,iAcr,iTc,i-Tr-KurAurTr-Tu,where ρr, Vr, and Cr are, respectively, the density, volume, and specific heat of the reflector, Kur and Aur are the heat transfer coefficient and heat transfer area between the reflector and the coolant in the riser, and Tu is the average coolant temperature inside the riser.(
3) Downcomer. The cold helium enters the downcomer from its upper entry and flows downward to cool the pebble-bed. The heat is mainly transferred to the helium inside the downcomer, which can be described as [24](6)Vidρd,idt=Wu-Wd,1,i=1,Wd,i-1-Wd,i,i=2,…,N,εViCpdρd,iTd,idt=Kd,iAd,iTc,i-Td,i+CpWuTu-Wd,1Td,1,i=1,CpWd,i-1Td,i-1-Wd,iTd,i,i=2,…,N,where ρd,i and Wd,i are, respectively, the density and flowrate of the helium inside node i of the downcomer and Cp is the helium specific heat at constant pressure. Then, substitute (6), the thermodynamics of the helium in the downcomer can be written as [24](7)εViρd,iCpdTd,idt=Kd,iAd,iTc,i-Td,i+CpWuTu-Td,1,i=1,CpWd,i-1Td,i-1-Td,i,i=2,…,N.(
4) Riser. In the channels of the riser, the cold helium flows upward to cool the reflector, and the heat transfer process can be described by [23, 24](8)Vudρudt=1-κWlh-Wu,VuCpdρuTudt=Cp1-κWlhTlh-WuTu+KurAurTr-Tu,where ρu, Wu, and Tu are, respectively, the density, flowrate, and temperature of the helium inside the riser, Wlh and Tlh are, respectively, the flowrate and temperature of the helium inside the lower header, and κ is the leakage ration of the lower header. Substitute (8), the thermodynamics of the helium inside the riser can be written as(9)VuρuCpdTudt=Cp1-κWlhTlh-Tu+KurAurTr-Tu.(
5) Lower Header. The lower header is assumed to be a well-mixed chamber which has no heat exchange with the environment, and the thermodynamics corresponding to the lower header can be expressed as [23, 24](10)VlhρlhdTlhdt=WlhTin-Tlh,where Tin is the temperature of the inlet cold helium, ρlh, Wlh, and Tlh are, respectively, the density, flowrate, and temperature of the helium inside the lower header, and Vlh is the volume of the low header.(
6) Outlet Header. The output header is assumed to be a well-mixed chamber, and its dynamics can be given as [23, 24](11)VohρohdTohdt=κWlhTlh-Toh+Wd,NTd,10-Toh,where ρoh, Wlh, and Toh are, respectively, the density, flowrate, and temperature of the helium inside the outlet header and Voh is the volume of the outlet header.


## 2.1.3. MHTGR Dynamics

The dynamic model of MHTGR is composed of neutron kinetics determined by (1)–(3) as well as thermal hydraulics given by (4), (5), (7), and (9)–(11). The dynamic behavior of this MHTGR model can be influenced by nodal numberN, and an improper nodal number leads to the fluctuation of the responses. Actually, this model is only utilized for control design and verification and is not for physical and thermal-hydraulic design. The nodal numbers are determined by comparing the response of this model with that of the model for physical and thermal-hydraulic design. When the error between responses is satisfactorily small, the corresponding nodal number is suitable. Here, the nodal numberN is chosen to be 10. Based on the simulation results to be given in Section 4, it will be seen that there is no numerical fluctuation in the responses, and all the transience is determined by power-level maneuver which show that the nodal number is properly selected.Since the dynamic model in this paper is tuned based upon the comparison with the model for reactor physical and thermal-hydraulic design, the stability of this model is also guaranteed by the model for reactor design. Moreover, it was shown in [27] that the MHTGR is open-loop globally asymptotically stable, and by the use of the shifted-ectropy-based self-stability analysis method given in [27], it can be shown that the model coupled by neutron kinetics and thermal hydraulics are open-loop globally asymptotically stable.


## 2.2. OTSG

Since the OTSG can provide superheated steam which does not satisfy the one-to-one map from pressure to temperature, the OTSG is more proper for building multi-SMR-based nuclear plants. Therefore, the dynamic model of the OTSG is very crucial for studying behavior of the SMR-based plant and also for verifying the operation and control strategies. In this section, a moving-boundary model [25] is introduced for describing OTSG dynamics in the cases of generating saturated and superheated steam [26]. The nodalization of the OTSG is given in Figure 4.Figure 4
Nodalization scheme of the OTSG.


### 2.2.1. Secondary Side

The dynamics of the OTSG secondary side is given by both the thermodynamics and pressure drops of the secondary fluid flow.(
1) Thermodynamics. The thermodynamics of the secondary side can be given by the following vector-valued nonlinear differential equation [26]:(12)Esxsdxsdt=fsxs,us,where (13)xs=l13hs5hs1Gs3T,us=Gs1hfwQs2Qs4T,Esxs=Es11xs012ρs2l13Es14xsEs21xs12ρs4l350Es24xs00τH0000τG,Es11xs=12ρs3hs1-hs3+Ps3-Ps1+l131-ρs2∂hs3∂Ps3Ps3ρs2gsin⁡θ+Fs2Gs22As22,Es14xs=Fs2Gs22l1322As221-ρs2∂hs3∂Ps3Ps3,Es21xs=12ρs3hs3-hs5+Ps5-Ps3+l352-ρs4∂hs3∂Ps3Ps3ρs2gsin⁡θ+Fs2Gs22As22-ρs4gsin⁡θ+Fs4Gs42As42,Es24xs=l351-12ρs4∂hs3∂Ps3Ps3Fs2Gs2l13As22+2Fs4Gs4l13As22,fsxs,us=Qs2As2+12Gs1As1+Gs3As3hs1-hs3Qs4+Gs4hs3-hs5As4hfw-hs1τHGs1-Gs3τGT,g is the gravitational acceleration (N·m/s2), both τG and τH are given positive constants denoting inertial effect, hfw is the enthalpy per unit mass (J/kg) of the feedwater flow, lij (i,j=1,…,5 and i<j) is the length between points i and j in Figure 4, and ρsi, Gsi, hsiQsi, and Asi are, respectively, the flow density (kg/m3), flowrate (kg/s), enthalpy per unit mass (J/kg), heat flux from the primary side, and cross section (m2) of node i of the secondary side. Here, l13 is the length of the subcooled section, l35 is that of the boiling section, and h1s is specific enthalpy of the feedwater. Due to the high speed of the steam flow, it is assumed that Gs3=Gs4=Gs5.(
2) Pressure Drops. The pressure drop of the fluid flow in the OTSG secondary side cannot be omitted, and the pressure of each node is given by [26](14)Ps3=Ps1-ρs2gl13sin⁡θ-Fs2l13Gs22,Ps5=Ps3-ρs4gl35sin⁡θ-Fs4l35Gs42,Ps4=Ps3+Ps52,Ps2=Ps1+Ps32,where θ is the helically coiled angle of the tubes, Fsi (i=2,4) is the friction factor of the two sections, and Ps1 is the feedwater pressure.


### 2.2.2. Tube Wall

There is no flow inside the metal tube wall between the primary and secondary sides, and the dynamics of the tube wall are determined by its thermodynamics, which can be expressed as [26](15)dTm2dt=Qp2-Qs2CmρmAm2l13+Tm4-Tm2l15dl13dt,dTm4dt=Qp4-Qs4CmρmAm4l35+Tm4-Tm2l15dl13dt,where Qp2 and Qp4 are, respectively, the heat flux from the primary helium flow to tube wall in the subcooled and boiling sections.


### 2.2.3. Primary Side

Since velocity of the primary helium flow is much higher than that of the two-phase flow inside the secondary side, the energy and temperature relationship of the primary side can be given by the following algebraic equation [26]:(16)Apxsxp=bpxm,xp,where(17)xp=Tp1Tp2Tp3Tp4T,Apxm=1-2100l35-l15l13-CpGp-Kp2l13CpGp000-CpGp-Kp4l35,bpxm,xp=00-Kp2l13Tm2-Kp4l35Tm4-CpGpTp5T,where Gp is the primary helium flowrate, Cp is the helium specific heat, Tpi is the helium temperature at point i (i=1,…,5), Kp2 and Kp4 are the heat transfer coefficients between the primary side and tube wall of the subcooled and boiling sections, respectively, and Ks2 and Ks4 are, respectively, the heat transfer coefficients between the tube wall and secondary side of the subcooled and boiling sections.


## 2.2.1. Secondary Side

The dynamics of the OTSG secondary side is given by both the thermodynamics and pressure drops of the secondary fluid flow.(
1) Thermodynamics. The thermodynamics of the secondary side can be given by the following vector-valued nonlinear differential equation [26]:(12)Esxsdxsdt=fsxs,us,where (13)xs=l13hs5hs1Gs3T,us=Gs1hfwQs2Qs4T,Esxs=Es11xs012ρs2l13Es14xsEs21xs12ρs4l350Es24xs00τH0000τG,Es11xs=12ρs3hs1-hs3+Ps3-Ps1+l131-ρs2∂hs3∂Ps3Ps3ρs2gsin⁡θ+Fs2Gs22As22,Es14xs=Fs2Gs22l1322As221-ρs2∂hs3∂Ps3Ps3,Es21xs=12ρs3hs3-hs5+Ps5-Ps3+l352-ρs4∂hs3∂Ps3Ps3ρs2gsin⁡θ+Fs2Gs22As22-ρs4gsin⁡θ+Fs4Gs42As42,Es24xs=l351-12ρs4∂hs3∂Ps3Ps3Fs2Gs2l13As22+2Fs4Gs4l13As22,fsxs,us=Qs2As2+12Gs1As1+Gs3As3hs1-hs3Qs4+Gs4hs3-hs5As4hfw-hs1τHGs1-Gs3τGT,g is the gravitational acceleration (N·m/s2), both τG and τH are given positive constants denoting inertial effect, hfw is the enthalpy per unit mass (J/kg) of the feedwater flow, lij (i,j=1,…,5 and i<j) is the length between points i and j in Figure 4, and ρsi, Gsi, hsiQsi, and Asi are, respectively, the flow density (kg/m3), flowrate (kg/s), enthalpy per unit mass (J/kg), heat flux from the primary side, and cross section (m2) of node i of the secondary side. Here, l13 is the length of the subcooled section, l35 is that of the boiling section, and h1s is specific enthalpy of the feedwater. Due to the high speed of the steam flow, it is assumed that Gs3=Gs4=Gs5.(
2) Pressure Drops. The pressure drop of the fluid flow in the OTSG secondary side cannot be omitted, and the pressure of each node is given by [26](14)Ps3=Ps1-ρs2gl13sin⁡θ-Fs2l13Gs22,Ps5=Ps3-ρs4gl35sin⁡θ-Fs4l35Gs42,Ps4=Ps3+Ps52,Ps2=Ps1+Ps32,where θ is the helically coiled angle of the tubes, Fsi (i=2,4) is the friction factor of the two sections, and Ps1 is the feedwater pressure.


## 2.2.2. Tube Wall

There is no flow inside the metal tube wall between the primary and secondary sides, and the dynamics of the tube wall are determined by its thermodynamics, which can be expressed as [26](15)dTm2dt=Qp2-Qs2CmρmAm2l13+Tm4-Tm2l15dl13dt,dTm4dt=Qp4-Qs4CmρmAm4l35+Tm4-Tm2l15dl13dt,where Qp2 and Qp4 are, respectively, the heat flux from the primary helium flow to tube wall in the subcooled and boiling sections.


## 2.2.3. Primary Side

Since velocity of the primary helium flow is much higher than that of the two-phase flow inside the secondary side, the energy and temperature relationship of the primary side can be given by the following algebraic equation [26]:(16)Apxsxp=bpxm,xp,where(17)xp=Tp1Tp2Tp3Tp4T,Apxm=1-2100l35-l15l13-CpGp-Kp2l13CpGp000-CpGp-Kp4l35,bpxm,xp=00-Kp2l13Tm2-Kp4l35Tm4-CpGpTp5T,where Gp is the primary helium flowrate, Cp is the helium specific heat, Tpi is the helium temperature at point i (i=1,…,5), Kp2 and Kp4 are the heat transfer coefficients between the primary side and tube wall of the subcooled and boiling sections, respectively, and Ks2 and Ks4 are, respectively, the heat transfer coefficients between the tube wall and secondary side of the subcooled and boiling sections.


## 2.3. Secondary-Loop FFN

After the hot helium leaves the reactor, it flows into the OTSG primary side where it transfers the heat to the secondary feedwater. By absorbing the heat of the primary loop, the feedwater is converted to the superheated steam in 571°C and about 13.9 MPa or so. The main steam flow at 13.24 MPa combined from the two NSSS modules then enters the steam turbine for producing electricity. The condensed flow is fed into the OTSG secondary side again after it goes through the grand heater, low pressure heaters (LPHs), deaerator, and high pressure heaters (HPHs). Thus, the hydraulics of the secondary-loop system can be modeled as the FFN shown in Figure5, where branches f1 and f2 are the fan branches of which everyone has a feedwater pump, branches 1 and 2 denote the secondary side of the OTSGs corresponding to the first and second NSSS module respectively, and branch 3 contains the equipment of the plant secondary loop such as the steam turbine, condenser, low pressure heaters, and deaerator.Figure 5
Fluid flow network of the secondary loop of the HTR-PM plant.The flowrate dynamics of branches 1, 2, and 3 can be given as [28, 29] (18)dQdt+KQD2R=KH,where Q=Q1,Q2,Q3T, QD2=diag⁡Q2, K=diag⁡K1,K2,K3, R=R1,R2,R3T, H=H1,H2,H3T, and Qj, Rj, Hj, and Kj are the flowrate, resistance, pressure drop, and inertia coefficient of branch j (j=1,2,3), respectively.Every FFN can be divided into a tree containing the fan branches and its complement, that is, the cotree whose branches are referred to as the links. From Figure5, the branches with solid line, that is, branches 3, f1, and f2, constitute the tree of this FFN, and the branches with dashed line, that is, branches 1 and 2, are the links which constitute the cotree. Based on the network property, the pressure drops and the flowrates of the tree branches are not independent and obey the following algebraic equations given by Kirchhoff’s voltage law (KVL) and Kirchhoff’s current law (KCL) [28]: (19)H1H2T=-H3+Hf1H3+Hf2T,Q3Qf1Qf2T=Q1+Q2Q1Q2T, where Hfj is the pressure drop of fan branch j (j=1,2) and can be further expressed as(20)Hfj=-Hdj+RfjQfj,j=1,2.Then, based on (18) and (19), the hydraulics of the FFN shown in Figure 4 can be given by differential-algebraic system [29]:(21)dQcdt=-KcRfQc+QcD2Rc+KcHd-H3η,K-H3=kcTHd-kfcTQc-kcTQcD2Rc+K3Q32R3,where Qc=Q1,Q2T, QcD2=diag⁡Q1,Q22, kc=K1,K2T, kfc=K1Rf1,K2Rf2T, K-=K1+K2+K3, Kc=diag⁡K1,K2, Rf=diag⁡Rf1,Rf2, and η=11T.


## 2.4. Model Integration

The above dynamic models of the MHTGR, OTSG, and secondary-loop FFN can be integrated together with other commonly utilized models such as the turbine and the synchronous generator to form the dynamic model for the whole plant. Based upon this model, a simulation code for the dynamic behavior and control characteristics of the HTR-PM plant is developed under MATLAB/Simulink platform, and the composition of this Simulink model is shown in Figure6. This model can incorporate the plant control laws proposed in the next section for giving the simulations of the operational and control characteristics of the HTR-PM plant, where the program modules named as “1# NSSS” and “2# NSSS” are used to simulate the dynamic behavior of the two NSSS modules, the program module named as “secondary-loop FFN” is used to simulate the thermal hydraulics for the secondary-loop FFN, the module “Secondary Thermal Feature” is used to simulate the thermodynamics of secondary-loop systems, and module “Generator” is used to simulate the electric generator dynamics. Furthermore, some important model parameters are given in Table 1.Table 1
Main parameters of the dynamic model.
Parameter
Name
Unit
Value
Λ
i (i=1,…,10)
Prompt neutron lifetime
s
7.66
e
-
4
β
Fraction of all the delayed neutrons
%
0.67
β1
Fraction of 1st group delayed neutrons
%
0.0256
β2
Fraction of 2nd group delayed neutrons
%
0.14
β3
Fraction of 3rd group delayed neutrons
%
0.13
β4
Fraction of 4th group delayed neutrons
%
0.27
β5
Fraction of 5th group delayed neutrons
%
0.086
β6
Fraction of 6th group delayed neutrons
%
0.017
λ1
Decay constant of 1st group precursor
s−1
0.0256
λ2
Decay constant of 2nd group precursor
s−1
0.14
λ3
Decay constant of 3rd group precursor
s−1
0.13
λ4
Decay constant of 4th group precursor
s−1
0.27
λ5
Decay constant of 5th group precursor
s−1
0.086
λ6
Decay constant of 6th group precursor
s−1
0.017
A=αij10×10
Coupling coefficients
3.5
7.0
1.8
7.1
5.5
2.2
6.9
3.9
3.1
7.0
3.3
3.7
7.0
2.9
4.2
7.0
2.8
4.5
7.0
2.6
4.7
7.0
2.4
5.0
7.0
2.2
5.6
3.5
×
1
0
-
3
α
f
,
i (i=1,…,10)
Reactivity coefficient of fuel temperature
Δ
K
/
K/°C
-
4.36
e
-
5
α
m
,
i (i=1,…,10)
Reactivity coefficient of moderator temperature
Δ
K
/
K/°C
-
0.94
e
-
5
αr
Reactivity coefficient of reflector temperature
Δ
K
/
K/°C
1.49
e
-
5
ε
Porosity of the pebble-bed
0.39Figure 6
Composition of the SIMULINK code.


## 3. Plant Power Control

Plant control is a key technique to provide safe, stable, and efficient operation for every nuclear plant and to balance the power supply and demand. Since the two NSSS modules of the HTR-PM plant are coupled together by the common secondary loop including the turbine/generator set, and since the side-by-side arranged MHTGR and OTSG of a NSSS module are tightly coupled with each other through the connecting pipes, it is more difficult to design a power control strategy for the two-modular HTR-PM plant than for those traditional single-modular traditional nuclear plants. Furthermore, the HTR-PM plant is essentially a large-scale and multi-input-multioutput (MIMO) nonlinear system whose complexity in its dynamics certainly leads to the complexity in its plant control strategy, while, in the practical engineering, a single control loop is usually utilized to regulate a single process parameter such as the nuclear power, primary flow, feedwater flow, helium temperature, and steam temperature. Since the combination of all these single control loops deeply affects the plant dynamic behavior, it is necessary to coordinate these control loops for satisfactory closed-loop stability as well as transient performance. In this section, the control design method is first reviewed, and the control realization including feedback loops and control laws is then given.


### 3.1. Control Design Method

The control design mainly lies in two aspects, that is, feedback loop design and control law design, which are tightly related to each other. There have been many process control system design approaches such as the relative gain analysis- (RGA-) based method. In this paper, both the feedback loops and control laws of HTR-PM control system are designed by the use of the physics-based NSSS control method proposed in [30–36] and module coordination control method proposed in [28, 29], which provide the globally asymptotic closed-loop stability through guaranteeing the convergence of Lyapunov functions determined by the shifted-ectropies of neutron kinetics and thermodynamics as well as the kinetic energy stored in secondary-loop FFN.


#### 3.1.1. NSSS Module Control

Both the MHTGR and OTSG are complex nonlinear dynamic systems, and their dynamics are tightly coupled with each other due to their side-by-side arrangement. To improve both the steady and transient performance, it is necessary to give proper control laws for the NSSS module composed of the MHTGR and OTSG by considering the nonlinear dynamics. The physics-based control (PBC) method is 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 PWR [30, 31], MHTGR [32–34], and OTSG [35]. Very recently, based upon the PBC method, a novel model-free MHTGR-based NSSS module control strategy was proposed in [36], which provides globally asymptotic stability (GAS) for the NSSS module if the feedwater temperature is constant and both the helium and feedwater flowrate are well regulated. This NSSS module control strategy can be realized by feedback loops given in Section 3.2.1 and simple control laws given in Section 3.3.1.


#### 3.1.2. Coordination Control of Multiple NSSS Modules

Through the multimodular scheme, the inherent safety of the MHTGR is applicable to the large-scale plants with any desired power ratings. Module coordination control is one of the most important characteristics of the multimodular nuclear plants such as the HTR-PM with comparison to the control of those traditional single module nuclear plants. It has been revealed in [29] that the module coordination problem is essentially a flowrate-pressure regulation problem of the secondary-loop FFN. Moreover, based on the dissipation structure of the FFNs, it has also been proved in [29] that feedback loops given in Section 3.2.2 and control laws given in Section 3.3.2 can provide globally input-to-state stability (ISS) and satisfactory transient performance for the cotree flowrates and tree branch pressure drops for the secondary-loop FFN.


### 3.2. Feedback Loops and Controller Functions

All the feedback loops of the plant control system of the HTR-PM plant are shown in Figure7. Since the feedback loops of the two NSSS modules are the same with each other, only those of the first NSSS module are illustrated. From Figure 2, we can see that the plant control system is constituted by the controllers for regulating the reactor nuclear power, primary helium flowrate, secondary feedwater flowrate, outlet helium temperature, steam temperature, reactor thermal power, and main steam pressure as well as the controller for discharge and bypass. Based on their functions, these controllers can be classified to the layers of module control, module coordination, and supervisory as well as bypass and discharge.Figure 7
Composition of HTR-PM plant control system.


#### 3.2.1. Module Control

The layer of module control is composed of the controllers for regulating the nuclear power, helium flowrate, reactor outlet helium temperature, steam temperature, and NSSS output thermal power. The functions of these five controllers are given as follows.(
1) Nuclear Power Controller. It regulates the neutron flux by generating both the given direction and velocity signals that determine the control rod movement. These signals are utilized to drive the stepping motors belonging to the reactor control rod mechanism for realizing the given movement. The movement of the control rods induces satisfactory reactivity to the reactor, which guarantees the stability and transient performance of the reactor neutron flux.(
2) Helium Flow Controller. Satisfactory primary helium circulation is the precondition of transferring the heat produced by the fission reaction to the secondary loop for steam generation. In order to provide enough pressure header for the primary helium circulation, the helium flow controller gives the referenced frequency for the variable frequency driver (VFD) of the induction motor that drives the primary helium blower.(
3) Feedwater Flow Controller.It regulates the secondary feedwater flowrates of the OTSGs through adjusting the rotational speed of the feedwater pumps by their hydraulic couplings. It is worth noting that the feedwater flowrate controller is responsible for the stability and transient performance of not only the NSSS modules but also the secondary-loop FFN of the HTR-PM plant.(
4) Helium Temperature Controller. It regulates the reactor outlet helium temperature by properly adjusting the setpoints of the nuclear power controller and (or) helium flow controller. Specifically, if there is only modification of the referenced relative nuclear power, then the helium temperature controller and the nuclear power controller constitute a cascaded feedback loop, where the former one is in the outer loop as the major controller, and the latter one is in the inner loop as minor controller.(
5) Steam Temperature Controller. It regulates the outlet steam temperature of the OTSG by properly adjusting the setpoints of the helium flowrate controller and (or) feedwater flowrate controller. Specifically, if its output is only the revision of the setpoints of the primary helium flowrate, then the steam temperature controller serves as the major regulator of the cascaded loop formed with the helium flow controller. Moreover, if its output is only the revision for the referenced feedwater flowrate, then the steam temperature controller and feedwater flow controller form a cascaded loop where the former one is the major regulator. In the power operation stage of the HTR-PM plant, the OTSG outlet steam temperature should be kept to be 571°C constantly.(
6) Thermal Power Controller. It is essentially a feedforward controller that gives the referenced feedwater flowrate based on both the setpoint of NSSS thermal power and the measurement of the temperatures and pressures of the feedwater and outlet superheated steam of each OTSG.


#### 3.2.2. Module Coordination

The two NSSS modules of the HTR-PM plant are tightly coupled with each other through the common secondary-loop FFN. To regulate the feedwater flowrate of 1# NSSS module, it is necessary to adjustHd1 by revising the rotating speed of the pump. However, from model (21), the variation of Hd1 influences not only the feedwater flowrate of 1# NSSS module but also pressure dropH3 which can be utilized to measure the main steam pressure. Then, from (21), the variation inH3 influences the dynamic behavior of the feedwater flowrate of 2# NSSS module, which induces the need of adjusting Hd2 for keeping this feedwater flowrate. Therefore, regulating the feedwater flowrate of one NSSS module needs to adjust the rotating speeds of the two pumps, and each NSSS feedwater flowrate can influence the main steam pressure. The module coordination is just the flowrate-pressure control of the secondary-loop FFN.It has been proved theoretically in [29] that the FFN flowrate regulation can be realized by distributed controllers of which each one only needs the measurement of the flowrate to be regulated. Thus, the NSSS feedwater flowrate controllers can well guarantee the stability and transient performance of the flowrates of the HTR-PM secondary-loop FFN shown in Figure 7. Moreover, in the power operation stage of the HTR-PM plant, the main steam pressure should be kept to be 13.24 MPa constantly. The stability of the main steam pressure is a basic operation requirement of the OTSG and also indicates the power balance between NSSS modules and the turbine/generator. Here, the main steam pressure of the HTR-PM plant is regulated by the turbine controller which provides the main steam pressure stability through adjusting the opening of the main steam valve. Thus, the collaboration of the NSSS feedwater flowrate controllers and the turbine controller realizes the flowrate-pressure regulation of the secondary-loop FFN, which achieves the function of NSSS module coordination and keeps the power balance between the NSSS modules and turbine.


#### 3.2.3. Supervisor

The supervisor gives the setpoint of the thermal power for every NSSS module according to the requirement of the power grid and has no feedback regulation function. Moreover, it also contains the map from the referenced NSSS thermal power to the referenced values of the NSSS process variables including the reactor nuclear power, primary helium flowrate, secondary feedwater flowrate, reactor outlet helium temperature, OTSG outlet steam temperature, and main steam pressure, which gives the setpoints of all the controllers. The supervisor can be viewed as an input-output (IO) interface between the power plant and the power grid so as to improve the power balance quality between the grid and plant.


#### 3.2.4. Bypass and Discharge

The function bypass and discharge guarantee the operationability of the HTR-PM plant in some tough cases. The bypass controller opens the bypass valve so that the OTSG outlet steam flows into the condenser in the cases of load rejection, reactor startup, turbine bypass, and so forth. The function of the bypass controller is essentially providing a virtual load to compensate the large difference between the thermal inertias of the NSSS and turbine. The discharge controller opens the discharge valve if the pressure of the outlet steam is larger than its maximal value and closes the valve when the steam pressure is lower than a given value.


### 3.3. Control Laws

The control laws for the nuclear power, primary helium flowrate, feedwater flowrate, MHTGR outlet helium temperature, OTSG outlet steam temperature, and main steam pressure, which have the feedback regulation function, are all given in this subsection. These control laws can be classified into control laws for NSSS modules and those for module coordination.


#### 3.3.1. Module Control Laws

The module control laws include the feedback regulation algorithms of the nuclear power, primary helium flowrate, feedwater flowrate, MHTGR outlet helium temperature, and OTSG outlet steam temperature, which are designed based on the physics-based NSSS control method proposed in [30–36].(
1) Nuclear Power Control Law. The nuclear power control laws is designed to be proportional-differential (PD) feedback regulation algorithm with the transfer function given by(22)GNs=uNsnr,es=kn,p+kn,dsτn,ds+1,where positive scalars kn,p and τn,d are, respectively, the proportional and differential feedback gains, τn,d is a small inertia time constant, uN is the setpoint of the control rod speed signal used to drive stepping motors of the control rod driving mechanism, and nr,e is the mismatch of the referenced relative nuclear power nr,r from its actual value nr; that is, nr,e=nr,r-nr. It is necessary to note that relative nuclear power is the normalized reactor neutron concentration, and the total thermal power is the outlet thermal power of the NSSS module. The steady values of the relative nuclear power is calibrated by that of total thermal power through experiment.(
2) Helium Flow Control Law. The feedback control of the primary helium flowrate adopts the proportional-integral (PI) law with its transfer function given by(23)GHFs=uHFsWH,es=khf,p+khf,is,where positive constants khf,p and khf,i are, respectively, the proportional and integral feedback gains, uHF is the referenced rotating speed of the primary helium blower, and WH,e is the error between the referenced and actual values of the primary helium flowrate.(
3) Feedwater Flow Control Law. The NSSS module feedwater flowrate is also regulated by PI control law, which is given by (24)GFWs=uFWsWfw,es=kfw,p+kfw,is,where positive scalars kfw,p and kfw,i are, respectively, the corresponding proportional and integral gains, uFW is the driving signal of the hydraulic coupling, and Wfw,e is the error between the referenced and actual values of the feedwater flowrate.(
4) MHTGR Outlet Helium Temperature Control Law. The reactor outlet helium temperature is regulated by the lead-lag feedback law that is given by transfer functions:(25)GHT1s=uHT1sTHO,es=kht1,p+kht,dskn,p+kn,ds,GHT2s=uHT2sTHO,es=-kht2,p+kht,is,where kn,p and kn,d are, respectively, the proportional and differential gains of the nuclear power control law, gains khtj,p (j=1,2), kht,d, and kht,i are positive constants determining the lead-lag performance, uHT1 is the revision to the setpoint of the relative nuclear power, uHT2 is the revision to the referenced primary helium flowrate, THO,e is the error between the setpoint and measurement of the reactor outlet helium temperature.(
5) OTSG Outlet Steam Temperature Control Law. The OTSG outlet steam temperature of every NSSS module is regulated by PI feedback laws:(26)GS1s=uS1sTS,es=ks1,p+ks1,is,GS2s=uS2sTS,es=-ks2,p+ks2,is,where positive constants ksj,p and ksj,i (j=1,2) are, respectively, the proportional and integral feedback gains, uS1 is the revision of the setpoint of the primary helium flowrate, uS2 is the revision of the setpoint of the secondary feedwater flowrate, and Ts,e is the error between the setpoint and measurement of the OTSG outlet steam temperature.


#### 3.3.2. Module Coordination Laws

The NSSS module coordination laws are just the flowrate-pressure feedback control laws for the secondary-loop FFN flowrates and pressure drop, which are designed based on the module coordination control method proposed in [28, 29].(
1) FFN Flowrate Control Law. As discussed in Section 3.1.2 and [28, 29], the flowrates of a FFN can be well regulated by controlling the flowrates of its cotree branches in a distributed way. From the structure of HTR-PM secondary-loop FFN shown in Figure 5, the OTSG secondary sides of the first and the second NSSS modules form branches 1 and 2 of this FFN which are just the cotree branches. Thus, from the graph properties of the FFN, all the branch flowrates can be well controlled or stabilized by regulating the feedwater flowrates of the OTSGs. Then, it is clear that OTSG feedwater flowrate control law (24) can realize the FFN flowrate control function, which means that flowrate control law (24) belongs to both the layers of module control and module coordination.(
2) FFN Pressure Drop Control Law. The pressure drop of the FFN tree branches (excluding the branches with pump) is also a crucial issue to guarantee the stability and operationability. The pressure drop H3 of the HTR-PM secondary-loop FFN shown in Figure 5 can be measured by the main steam pressure which is regulated by the turbine control law with transfer function given by(27)GSPs=uSPsPSs=kp,p+kp,is.Here, in (27), positive scalars kp,p and kp,i are, respectively, the proportional and integral gains, uSP is the setpoint of the steam valve opening, and PS,e is the mismatch of the referenced main steam pressure from its actual value.


#### 3.3.3. Parameters of Control Laws

The plant control performance relies on the two aspects, that is, the performances of both the module control and module coordination. The former one is guaranteed by the physics-based NSSS control approach, and the latter one is provided by the FFN flowrate-pressure control method. All the control laws take the form as PI or PD feedback laws. For PI control laws, the proportional gain is larger, and the response is faster, but if it is too large, the control action will hit the saturation margin, which leads to the deterioration of closed-loop stability. Moreover, the integral gain is larger, and the transition period to the steady state is shorter, but the overshoot is larger. Usually, proportional gain is chosen not larger than 1.0, and the integral gain is chosen not larger than 0.1. For the PD law for nuclear power regulation, the proportional gain is chosen not larger than 0.5, and the differential gain is chosen not larger than 3.0.


## 3.1. Control Design Method

The control design mainly lies in two aspects, that is, feedback loop design and control law design, which are tightly related to each other. There have been many process control system design approaches such as the relative gain analysis- (RGA-) based method. In this paper, both the feedback loops and control laws of HTR-PM control system are designed by the use of the physics-based NSSS control method proposed in [30–36] and module coordination control method proposed in [28, 29], which provide the globally asymptotic closed-loop stability through guaranteeing the convergence of Lyapunov functions determined by the shifted-ectropies of neutron kinetics and thermodynamics as well as the kinetic energy stored in secondary-loop FFN.


### 3.1.1. NSSS Module Control

Both the MHTGR and OTSG are complex nonlinear dynamic systems, and their dynamics are tightly coupled with each other due to their side-by-side arrangement. To improve both the steady and transient performance, it is necessary to give proper control laws for the NSSS module composed of the MHTGR and OTSG by considering the nonlinear dynamics. The physics-based control (PBC) method is 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 PWR [30, 31], MHTGR [32–34], and OTSG [35]. Very recently, based upon the PBC method, a novel model-free MHTGR-based NSSS module control strategy was proposed in [36], which provides globally asymptotic stability (GAS) for the NSSS module if the feedwater temperature is constant and both the helium and feedwater flowrate are well regulated. This NSSS module control strategy can be realized by feedback loops given in Section 3.2.1 and simple control laws given in Section 3.3.1.


### 3.1.2. Coordination Control of Multiple NSSS Modules

Through the multimodular scheme, the inherent safety of the MHTGR is applicable to the large-scale plants with any desired power ratings. Module coordination control is one of the most important characteristics of the multimodular nuclear plants such as the HTR-PM with comparison to the control of those traditional single module nuclear plants. It has been revealed in [29] that the module coordination problem is essentially a flowrate-pressure regulation problem of the secondary-loop FFN. Moreover, based on the dissipation structure of the FFNs, it has also been proved in [29] that feedback loops given in Section 3.2.2 and control laws given in Section 3.3.2 can provide globally input-to-state stability (ISS) and satisfactory transient performance for the cotree flowrates and tree branch pressure drops for the secondary-loop FFN.


## 3.1.1. NSSS Module Control

Both the MHTGR and OTSG are complex nonlinear dynamic systems, and their dynamics are tightly coupled with each other due to their side-by-side arrangement. To improve both the steady and transient performance, it is necessary to give proper control laws for the NSSS module composed of the MHTGR and OTSG by considering the nonlinear dynamics. The physics-based control (PBC) method is 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 PWR [30, 31], MHTGR [32–34], and OTSG [35]. Very recently, based upon the PBC method, a novel model-free MHTGR-based NSSS module control strategy was proposed in [36], which provides globally asymptotic stability (GAS) for the NSSS module if the feedwater temperature is constant and both the helium and feedwater flowrate are well regulated. This NSSS module control strategy can be realized by feedback loops given in Section 3.2.1 and simple control laws given in Section 3.3.1.


## 3.1.2. Coordination Control of Multiple NSSS Modules

Through the multimodular scheme, the inherent safety of the MHTGR is applicable to the large-scale plants with any desired power ratings. Module coordination control is one of the most important characteristics of the multimodular nuclear plants such as the HTR-PM with comparison to the control of those traditional single module nuclear plants. It has been revealed in [29] that the module coordination problem is essentially a flowrate-pressure regulation problem of the secondary-loop FFN. Moreover, based on the dissipation structure of the FFNs, it has also been proved in [29] that feedback loops given in Section 3.2.2 and control laws given in Section 3.3.2 can provide globally input-to-state stability (ISS) and satisfactory transient performance for the cotree flowrates and tree branch pressure drops for the secondary-loop FFN.


## 3.2. Feedback Loops and Controller Functions

All the feedback loops of the plant control system of the HTR-PM plant are shown in Figure7. Since the feedback loops of the two NSSS modules are the same with each other, only those of the first NSSS module are illustrated. From Figure 2, we can see that the plant control system is constituted by the controllers for regulating the reactor nuclear power, primary helium flowrate, secondary feedwater flowrate, outlet helium temperature, steam temperature, reactor thermal power, and main steam pressure as well as the controller for discharge and bypass. Based on their functions, these controllers can be classified to the layers of module control, module coordination, and supervisory as well as bypass and discharge.Figure 7
Composition of HTR-PM plant control system.


### 3.2.1. Module Control

The layer of module control is composed of the controllers for regulating the nuclear power, helium flowrate, reactor outlet helium temperature, steam temperature, and NSSS output thermal power. The functions of these five controllers are given as follows.(
1) Nuclear Power Controller. It regulates the neutron flux by generating both the given direction and velocity signals that determine the control rod movement. These signals are utilized to drive the stepping motors belonging to the reactor control rod mechanism for realizing the given movement. The movement of the control rods induces satisfactory reactivity to the reactor, which guarantees the stability and transient performance of the reactor neutron flux.(
2) Helium Flow Controller. Satisfactory primary helium circulation is the precondition of transferring the heat produced by the fission reaction to the secondary loop for steam generation. In order to provide enough pressure header for the primary helium circulation, the helium flow controller gives the referenced frequency for the variable frequency driver (VFD) of the induction motor that drives the primary helium blower.(
3) Feedwater Flow Controller.It regulates the secondary feedwater flowrates of the OTSGs through adjusting the rotational speed of the feedwater pumps by their hydraulic couplings. It is worth noting that the feedwater flowrate controller is responsible for the stability and transient performance of not only the NSSS modules but also the secondary-loop FFN of the HTR-PM plant.(
4) Helium Temperature Controller. It regulates the reactor outlet helium temperature by properly adjusting the setpoints of the nuclear power controller and (or) helium flow controller. Specifically, if there is only modification of the referenced relative nuclear power, then the helium temperature controller and the nuclear power controller constitute a cascaded feedback loop, where the former one is in the outer loop as the major controller, and the latter one is in the inner loop as minor controller.(
5) Steam Temperature Controller. It regulates the outlet steam temperature of the OTSG by properly adjusting the setpoints of the helium flowrate controller and (or) feedwater flowrate controller. Specifically, if its output is only the revision of the setpoints of the primary helium flowrate, then the steam temperature controller serves as the major regulator of the cascaded loop formed with the helium flow controller. Moreover, if its output is only the revision for the referenced feedwater flowrate, then the steam temperature controller and feedwater flow controller form a cascaded loop where the former one is the major regulator. In the power operation stage of the HTR-PM plant, the OTSG outlet steam temperature should be kept to be 571°C constantly.(
6) Thermal Power Controller. It is essentially a feedforward controller that gives the referenced feedwater flowrate based on both the setpoint of NSSS thermal power and the measurement of the temperatures and pressures of the feedwater and outlet superheated steam of each OTSG.


### 3.2.2. Module Coordination

The two NSSS modules of the HTR-PM plant are tightly coupled with each other through the common secondary-loop FFN. To regulate the feedwater flowrate of 1# NSSS module, it is necessary to adjustHd1 by revising the rotating speed of the pump. However, from model (21), the variation of Hd1 influences not only the feedwater flowrate of 1# NSSS module but also pressure dropH3 which can be utilized to measure the main steam pressure. Then, from (21), the variation inH3 influences the dynamic behavior of the feedwater flowrate of 2# NSSS module, which induces the need of adjusting Hd2 for keeping this feedwater flowrate. Therefore, regulating the feedwater flowrate of one NSSS module needs to adjust the rotating speeds of the two pumps, and each NSSS feedwater flowrate can influence the main steam pressure. The module coordination is just the flowrate-pressure control of the secondary-loop FFN.It has been proved theoretically in [29] that the FFN flowrate regulation can be realized by distributed controllers of which each one only needs the measurement of the flowrate to be regulated. Thus, the NSSS feedwater flowrate controllers can well guarantee the stability and transient performance of the flowrates of the HTR-PM secondary-loop FFN shown in Figure 7. Moreover, in the power operation stage of the HTR-PM plant, the main steam pressure should be kept to be 13.24 MPa constantly. The stability of the main steam pressure is a basic operation requirement of the OTSG and also indicates the power balance between NSSS modules and the turbine/generator. Here, the main steam pressure of the HTR-PM plant is regulated by the turbine controller which provides the main steam pressure stability through adjusting the opening of the main steam valve. Thus, the collaboration of the NSSS feedwater flowrate controllers and the turbine controller realizes the flowrate-pressure regulation of the secondary-loop FFN, which achieves the function of NSSS module coordination and keeps the power balance between the NSSS modules and turbine.


### 3.2.3. Supervisor

The supervisor gives the setpoint of the thermal power for every NSSS module according to the requirement of the power grid and has no feedback regulation function. Moreover, it also contains the map from the referenced NSSS thermal power to the referenced values of the NSSS process variables including the reactor nuclear power, primary helium flowrate, secondary feedwater flowrate, reactor outlet helium temperature, OTSG outlet steam temperature, and main steam pressure, which gives the setpoints of all the controllers. The supervisor can be viewed as an input-output (IO) interface between the power plant and the power grid so as to improve the power balance quality between the grid and plant.


### 3.2.4. Bypass and Discharge

The function bypass and discharge guarantee the operationability of the HTR-PM plant in some tough cases. The bypass controller opens the bypass valve so that the OTSG outlet steam flows into the condenser in the cases of load rejection, reactor startup, turbine bypass, and so forth. The function of the bypass controller is essentially providing a virtual load to compensate the large difference between the thermal inertias of the NSSS and turbine. The discharge controller opens the discharge valve if the pressure of the outlet steam is larger than its maximal value and closes the valve when the steam pressure is lower than a given value.


## 3.2.1. Module Control

The layer of module control is composed of the controllers for regulating the nuclear power, helium flowrate, reactor outlet helium temperature, steam temperature, and NSSS output thermal power. The functions of these five controllers are given as follows.(
1) Nuclear Power Controller. It regulates the neutron flux by generating both the given direction and velocity signals that determine the control rod movement. These signals are utilized to drive the stepping motors belonging to the reactor control rod mechanism for realizing the given movement. The movement of the control rods induces satisfactory reactivity to the reactor, which guarantees the stability and transient performance of the reactor neutron flux.(
2) Helium Flow Controller. Satisfactory primary helium circulation is the precondition of transferring the heat produced by the fission reaction to the secondary loop for steam generation. In order to provide enough pressure header for the primary helium circulation, the helium flow controller gives the referenced frequency for the variable frequency driver (VFD) of the induction motor that drives the primary helium blower.(
3) Feedwater Flow Controller.It regulates the secondary feedwater flowrates of the OTSGs through adjusting the rotational speed of the feedwater pumps by their hydraulic couplings. It is worth noting that the feedwater flowrate controller is responsible for the stability and transient performance of not only the NSSS modules but also the secondary-loop FFN of the HTR-PM plant.(
4) Helium Temperature Controller. It regulates the reactor outlet helium temperature by properly adjusting the setpoints of the nuclear power controller and (or) helium flow controller. Specifically, if there is only modification of the referenced relative nuclear power, then the helium temperature controller and the nuclear power controller constitute a cascaded feedback loop, where the former one is in the outer loop as the major controller, and the latter one is in the inner loop as minor controller.(
5) Steam Temperature Controller. It regulates the outlet steam temperature of the OTSG by properly adjusting the setpoints of the helium flowrate controller and (or) feedwater flowrate controller. Specifically, if its output is only the revision of the setpoints of the primary helium flowrate, then the steam temperature controller serves as the major regulator of the cascaded loop formed with the helium flow controller. Moreover, if its output is only the revision for the referenced feedwater flowrate, then the steam temperature controller and feedwater flow controller form a cascaded loop where the former one is the major regulator. In the power operation stage of the HTR-PM plant, the OTSG outlet steam temperature should be kept to be 571°C constantly.(
6) Thermal Power Controller. It is essentially a feedforward controller that gives the referenced feedwater flowrate based on both the setpoint of NSSS thermal power and the measurement of the temperatures and pressures of the feedwater and outlet superheated steam of each OTSG.


## 3.2.2. Module Coordination

The two NSSS modules of the HTR-PM plant are tightly coupled with each other through the common secondary-loop FFN. To regulate the feedwater flowrate of 1# NSSS module, it is necessary to adjustHd1 by revising the rotating speed of the pump. However, from model (21), the variation of Hd1 influences not only the feedwater flowrate of 1# NSSS module but also pressure dropH3 which can be utilized to measure the main steam pressure. Then, from (21), the variation inH3 influences the dynamic behavior of the feedwater flowrate of 2# NSSS module, which induces the need of adjusting Hd2 for keeping this feedwater flowrate. Therefore, regulating the feedwater flowrate of one NSSS module needs to adjust the rotating speeds of the two pumps, and each NSSS feedwater flowrate can influence the main steam pressure. The module coordination is just the flowrate-pressure control of the secondary-loop FFN.It has been proved theoretically in [29] that the FFN flowrate regulation can be realized by distributed controllers of which each one only needs the measurement of the flowrate to be regulated. Thus, the NSSS feedwater flowrate controllers can well guarantee the stability and transient performance of the flowrates of the HTR-PM secondary-loop FFN shown in Figure 7. Moreover, in the power operation stage of the HTR-PM plant, the main steam pressure should be kept to be 13.24 MPa constantly. The stability of the main steam pressure is a basic operation requirement of the OTSG and also indicates the power balance between NSSS modules and the turbine/generator. Here, the main steam pressure of the HTR-PM plant is regulated by the turbine controller which provides the main steam pressure stability through adjusting the opening of the main steam valve. Thus, the collaboration of the NSSS feedwater flowrate controllers and the turbine controller realizes the flowrate-pressure regulation of the secondary-loop FFN, which achieves the function of NSSS module coordination and keeps the power balance between the NSSS modules and turbine.


## 3.2.3. Supervisor

The supervisor gives the setpoint of the thermal power for every NSSS module according to the requirement of the power grid and has no feedback regulation function. Moreover, it also contains the map from the referenced NSSS thermal power to the referenced values of the NSSS process variables including the reactor nuclear power, primary helium flowrate, secondary feedwater flowrate, reactor outlet helium temperature, OTSG outlet steam temperature, and main steam pressure, which gives the setpoints of all the controllers. The supervisor can be viewed as an input-output (IO) interface between the power plant and the power grid so as to improve the power balance quality between the grid and plant.


## 3.2.4. Bypass and Discharge

The function bypass and discharge guarantee the operationability of the HTR-PM plant in some tough cases. The bypass controller opens the bypass valve so that the OTSG outlet steam flows into the condenser in the cases of load rejection, reactor startup, turbine bypass, and so forth. The function of the bypass controller is essentially providing a virtual load to compensate the large difference between the thermal inertias of the NSSS and turbine. The discharge controller opens the discharge valve if the pressure of the outlet steam is larger than its maximal value and closes the valve when the steam pressure is lower than a given value.


## 3.3. Control Laws

The control laws for the nuclear power, primary helium flowrate, feedwater flowrate, MHTGR outlet helium temperature, OTSG outlet steam temperature, and main steam pressure, which have the feedback regulation function, are all given in this subsection. These control laws can be classified into control laws for NSSS modules and those for module coordination.


### 3.3.1. Module Control Laws

The module control laws include the feedback regulation algorithms of the nuclear power, primary helium flowrate, feedwater flowrate, MHTGR outlet helium temperature, and OTSG outlet steam temperature, which are designed based on the physics-based NSSS control method proposed in [30–36].(
1) Nuclear Power Control Law. The nuclear power control laws is designed to be proportional-differential (PD) feedback regulation algorithm with the transfer function given by(22)GNs=uNsnr,es=kn,p+kn,dsτn,ds+1,where positive scalars kn,p and τn,d are, respectively, the proportional and differential feedback gains, τn,d is a small inertia time constant, uN is the setpoint of the control rod speed signal used to drive stepping motors of the control rod driving mechanism, and nr,e is the mismatch of the referenced relative nuclear power nr,r from its actual value nr; that is, nr,e=nr,r-nr. It is necessary to note that relative nuclear power is the normalized reactor neutron concentration, and the total thermal power is the outlet thermal power of the NSSS module. The steady values of the relative nuclear power is calibrated by that of total thermal power through experiment.(
2) Helium Flow Control Law. The feedback control of the primary helium flowrate adopts the proportional-integral (PI) law with its transfer function given by(23)GHFs=uHFsWH,es=khf,p+khf,is,where positive constants khf,p and khf,i are, respectively, the proportional and integral feedback gains, uHF is the referenced rotating speed of the primary helium blower, and WH,e is the error between the referenced and actual values of the primary helium flowrate.(
3) Feedwater Flow Control Law. The NSSS module feedwater flowrate is also regulated by PI control law, which is given by (24)GFWs=uFWsWfw,es=kfw,p+kfw,is,where positive scalars kfw,p and kfw,i are, respectively, the corresponding proportional and integral gains, uFW is the driving signal of the hydraulic coupling, and Wfw,e is the error between the referenced and actual values of the feedwater flowrate.(
4) MHTGR Outlet Helium Temperature Control Law. The reactor outlet helium temperature is regulated by the lead-lag feedback law that is given by transfer functions:(25)GHT1s=uHT1sTHO,es=kht1,p+kht,dskn,p+kn,ds,GHT2s=uHT2sTHO,es=-kht2,p+kht,is,where kn,p and kn,d are, respectively, the proportional and differential gains of the nuclear power control law, gains khtj,p (j=1,2), kht,d, and kht,i are positive constants determining the lead-lag performance, uHT1 is the revision to the setpoint of the relative nuclear power, uHT2 is the revision to the referenced primary helium flowrate, THO,e is the error between the setpoint and measurement of the reactor outlet helium temperature.(
5) OTSG Outlet Steam Temperature Control Law. The OTSG outlet steam temperature of every NSSS module is regulated by PI feedback laws:(26)GS1s=uS1sTS,es=ks1,p+ks1,is,GS2s=uS2sTS,es=-ks2,p+ks2,is,where positive constants ksj,p and ksj,i (j=1,2) are, respectively, the proportional and integral feedback gains, uS1 is the revision of the setpoint of the primary helium flowrate, uS2 is the revision of the setpoint of the secondary feedwater flowrate, and Ts,e is the error between the setpoint and measurement of the OTSG outlet steam temperature.


### 3.3.2. Module Coordination Laws

The NSSS module coordination laws are just the flowrate-pressure feedback control laws for the secondary-loop FFN flowrates and pressure drop, which are designed based on the module coordination control method proposed in [28, 29].(
1) FFN Flowrate Control Law. As discussed in Section 3.1.2 and [28, 29], the flowrates of a FFN can be well regulated by controlling the flowrates of its cotree branches in a distributed way. From the structure of HTR-PM secondary-loop FFN shown in Figure 5, the OTSG secondary sides of the first and the second NSSS modules form branches 1 and 2 of this FFN which are just the cotree branches. Thus, from the graph properties of the FFN, all the branch flowrates can be well controlled or stabilized by regulating the feedwater flowrates of the OTSGs. Then, it is clear that OTSG feedwater flowrate control law (24) can realize the FFN flowrate control function, which means that flowrate control law (24) belongs to both the layers of module control and module coordination.(
2) FFN Pressure Drop Control Law. The pressure drop of the FFN tree branches (excluding the branches with pump) is also a crucial issue to guarantee the stability and operationability. The pressure drop H3 of the HTR-PM secondary-loop FFN shown in Figure 5 can be measured by the main steam pressure which is regulated by the turbine control law with transfer function given by(27)GSPs=uSPsPSs=kp,p+kp,is.Here, in (27), positive scalars kp,p and kp,i are, respectively, the proportional and integral gains, uSP is the setpoint of the steam valve opening, and PS,e is the mismatch of the referenced main steam pressure from its actual value.


### 3.3.3. Parameters of Control Laws

The plant control performance relies on the two aspects, that is, the performances of both the module control and module coordination. The former one is guaranteed by the physics-based NSSS control approach, and the latter one is provided by the FFN flowrate-pressure control method. All the control laws take the form as PI or PD feedback laws. For PI control laws, the proportional gain is larger, and the response is faster, but if it is too large, the control action will hit the saturation margin, which leads to the deterioration of closed-loop stability. Moreover, the integral gain is larger, and the transition period to the steady state is shorter, but the overshoot is larger. Usually, proportional gain is chosen not larger than 1.0, and the integral gain is chosen not larger than 0.1. For the PD law for nuclear power regulation, the proportional gain is chosen not larger than 0.5, and the differential gain is chosen not larger than 3.0.


## 3.3.1. Module Control Laws

The module control laws include the feedback regulation algorithms of the nuclear power, primary helium flowrate, feedwater flowrate, MHTGR outlet helium temperature, and OTSG outlet steam temperature, which are designed based on the physics-based NSSS control method proposed in [30–36].(
1) Nuclear Power Control Law. The nuclear power control laws is designed to be proportional-differential (PD) feedback regulation algorithm with the transfer function given by(22)GNs=uNsnr,es=kn,p+kn,dsτn,ds+1,where positive scalars kn,p and τn,d are, respectively, the proportional and differential feedback gains, τn,d is a small inertia time constant, uN is the setpoint of the control rod speed signal used to drive stepping motors of the control rod driving mechanism, and nr,e is the mismatch of the referenced relative nuclear power nr,r from its actual value nr; that is, nr,e=nr,r-nr. It is necessary to note that relative nuclear power is the normalized reactor neutron concentration, and the total thermal power is the outlet thermal power of the NSSS module. The steady values of the relative nuclear power is calibrated by that of total thermal power through experiment.(
2) Helium Flow Control Law. The feedback control of the primary helium flowrate adopts the proportional-integral (PI) law with its transfer function given by(23)GHFs=uHFsWH,es=khf,p+khf,is,where positive constants khf,p and khf,i are, respectively, the proportional and integral feedback gains, uHF is the referenced rotating speed of the primary helium blower, and WH,e is the error between the referenced and actual values of the primary helium flowrate.(
3) Feedwater Flow Control Law. The NSSS module feedwater flowrate is also regulated by PI control law, which is given by (24)GFWs=uFWsWfw,es=kfw,p+kfw,is,where positive scalars kfw,p and kfw,i are, respectively, the corresponding proportional and integral gains, uFW is the driving signal of the hydraulic coupling, and Wfw,e is the error between the referenced and actual values of the feedwater flowrate.(
4) MHTGR Outlet Helium Temperature Control Law. The reactor outlet helium temperature is regulated by the lead-lag feedback law that is given by transfer functions:(25)GHT1s=uHT1sTHO,es=kht1,p+kht,dskn,p+kn,ds,GHT2s=uHT2sTHO,es=-kht2,p+kht,is,where kn,p and kn,d are, respectively, the proportional and differential gains of the nuclear power control law, gains khtj,p (j=1,2), kht,d, and kht,i are positive constants determining the lead-lag performance, uHT1 is the revision to the setpoint of the relative nuclear power, uHT2 is the revision to the referenced primary helium flowrate, THO,e is the error between the setpoint and measurement of the reactor outlet helium temperature.(
5) OTSG Outlet Steam Temperature Control Law. The OTSG outlet steam temperature of every NSSS module is regulated by PI feedback laws:(26)GS1s=uS1sTS,es=ks1,p+ks1,is,GS2s=uS2sTS,es=-ks2,p+ks2,is,where positive constants ksj,p and ksj,i (j=1,2) are, respectively, the proportional and integral feedback gains, uS1 is the revision of the setpoint of the primary helium flowrate, uS2 is the revision of the setpoint of the secondary feedwater flowrate, and Ts,e is the error between the setpoint and measurement of the OTSG outlet steam temperature.


## 3.3.2. Module Coordination Laws

The NSSS module coordination laws are just the flowrate-pressure feedback control laws for the secondary-loop FFN flowrates and pressure drop, which are designed based on the module coordination control method proposed in [28, 29].(
1) FFN Flowrate Control Law. As discussed in Section 3.1.2 and [28, 29], the flowrates of a FFN can be well regulated by controlling the flowrates of its cotree branches in a distributed way. From the structure of HTR-PM secondary-loop FFN shown in Figure 5, the OTSG secondary sides of the first and the second NSSS modules form branches 1 and 2 of this FFN which are just the cotree branches. Thus, from the graph properties of the FFN, all the branch flowrates can be well controlled or stabilized by regulating the feedwater flowrates of the OTSGs. Then, it is clear that OTSG feedwater flowrate control law (24) can realize the FFN flowrate control function, which means that flowrate control law (24) belongs to both the layers of module control and module coordination.(
2) FFN Pressure Drop Control Law. The pressure drop of the FFN tree branches (excluding the branches with pump) is also a crucial issue to guarantee the stability and operationability. The pressure drop H3 of the HTR-PM secondary-loop FFN shown in Figure 5 can be measured by the main steam pressure which is regulated by the turbine control law with transfer function given by(27)GSPs=uSPsPSs=kp,p+kp,is.Here, in (27), positive scalars kp,p and kp,i are, respectively, the proportional and integral gains, uSP is the setpoint of the steam valve opening, and PS,e is the mismatch of the referenced main steam pressure from its actual value.


## 3.3.3. Parameters of Control Laws

The plant control performance relies on the two aspects, that is, the performances of both the module control and module coordination. The former one is guaranteed by the physics-based NSSS control approach, and the latter one is provided by the FFN flowrate-pressure control method. All the control laws take the form as PI or PD feedback laws. For PI control laws, the proportional gain is larger, and the response is faster, but if it is too large, the control action will hit the saturation margin, which leads to the deterioration of closed-loop stability. Moreover, the integral gain is larger, and the transition period to the steady state is shorter, but the overshoot is larger. Usually, proportional gain is chosen not larger than 1.0, and the integral gain is chosen not larger than 0.1. For the PD law for nuclear power regulation, the proportional gain is chosen not larger than 0.5, and the differential gain is chosen not larger than 3.0.


## 4. Simulation Results with Discussions

In this section, the dynamic behavior of the HTR-PM plant in both the cases of power step and ramp are studied through the numerical simulation. This study is performed based on the MATLAB/Simulink simulation code shown in Figure6 that integrates both the plant model introduced in Section 2 and control strategy given in Section 3. Moreover, necessary discussions are also given in this section.


### 4.1. Simulation Results

#### 4.1.1. Power Step

The dynamic responses of the power step between 100% and 95% plant full power (PFP) are obtained through numerical simulation.(
1) Power Step-Down. Initially, the HTR-PM plant remains at the steady state of 100% PFP which means that the two NSSS modules are running at the 100% reactor full power (RFP). At 5000 s, the thermal power setpoint of the first (i.e., 1#) NSSS module steps down from 100% to 95% RFP, and then at 10000 s, the referenced thermal power of the second (i.e., 2#) NSSS module steps down from 100% to 95% RFP. The dynamic responses of the relative nuclear powers, reactor outlet helium temperatures, OTSG outlet steam temperatures, and thermal powers of 1# and 2# NSSS modules are given in Figure 8. The responses of flowrates Qi (i=1,2,3), feedwater pressures -Hfk (k=1,2), and main steam pressure H3 are illustrated in Figure 9.Figure 8
Responses of 1# and 2# NSSS in the case of power step-down:nri, Tcouti, Tsi, and Pthi are, respectively, the relative nuclear power, MHTGR outlet helium temperature, OTSG outlet steam temperature, and total thermal power of the ith NSSS, i=1,2.Figure 9
Responses of secondary-loop FFN flowrates and pressure drops in the case of power step-down.(
2) Power Step-Up. Initially, the HTR-PM plant remains at the steady state of 95% PFP; that is, the two NSSS modules operate at 95% RFP. The thermal power setpoints of 1# and 2# NSSS modules step up from 95% to 100% RFP at 5000 s and 10000 s, respectively. The dynamic responses of the relative nuclear powers, reactor outlet helium temperatures, OTSG outlet steam temperatures, and thermal powers of 1# and 2# NSSS modules are given in Figure 10. The responses of Qi (i=1,2,3), -Hfk (k=1,2), andH3 are shown in Figure 11.Figure 10
Responses of 1# and 2# NSSS in the case of power step-up:nri, Tcouti, Tsi, and Pthi are, respectively, the relative nuclear power, MHTGR outlet helium temperature, OTSG outlet steam temperature, and total thermal power of the ith NSSS, i=1,2.Figure 11
Responses of secondary-loop FFN flowrates and pressure drops in the case of power step-up.


#### 4.1.2. Power Ramp

The responses of the steady power ramp between 100% and 50% plant full power (PFP) are given by the simulation study.(
1) Power Ramp-Down. Initially, the HTR-PM plant remains at the steady state of 100% PFP; that is, the two NSSS modules operate at 100% RFP. At 5000 s, 1# NSSS module maneuvers from 100% to 50% RFP steadily in 10 minutes, and then at 10000 s, the referenced thermal power of 2# NSSS module changes steadily from 100% to 50% RFP in 10 minutes, which is just the power decrease process from 100% to 50% PFP of the HTR-PM plant. The responses of the relative nuclear powers, reactor outlet helium temperatures, OTSG outlet steam temperatures, and thermal powers of 1# and 2# NSSS modules are given in Figure 12. The responses of flowrates Qi (i=1,2,3), feedwater headers -Hfk (k=1,2), and main steam pressureH3 are illustrated in Figure 13.Figure 12
Responses of 1# and 2# NSSS in the case of power ramp-down:nri, Tcouti, Tsi, and Pthi are, respectively, the relative nuclear power, MHTGR outlet helium temperature, OTSG outlet steam temperature, and total thermal power of the ith NSSS, i=1,2.Figure 13
Responses of secondary-loop FFN flowrates and pressure drops in the case of power ramp-down.(
2) Power Ramp-Up. Initially, the HTR-PM plant remains at the steady state of 95% PFP; that is, the two NSSS modules operate at 95% RFP. At 5000 s and 10000 s, 1# and 2# NSSS modules, respectively, maneuver steadily from 50% to 100% RFP in 10 minutes, which is the process of HTR-PM power increase from 50% to 100% PFP. The responses of the relative nuclear powers, reactor outlet helium temperatures, OTSG outlet steam temperatures, and thermal powers of 1# and 2# NSSS modules are given in Figure 14. The responses of flowrates Qi (i=1,2,3), -Hfk (k=1,2), andH3 are illustrated in Figure 15.Figure 14
Responses of 1# and 2# NSSS in the case of power ramp-up:nri, Tcouti, Tsi, and Pthi are, respectively, the relative nuclear power, MHTGR outlet helium temperature, OTSG outlet steam temperature, and total thermal power of the ith NSSS, i=1,2.Figure 15
Responses of secondary-loop FFN flowrates and pressure drops in the case of power ramp-up.


### 4.2. Discussions

The variation of the plant power is realized by the successive power maneuver of 1# and 2# NSSS modules, and at most one NSSS module can be in the state of power maneuver at any time. Generally, from Figures8–15, the control laws for the modules and those for the module coordination guarantee satisfactory closed-loop stability and transient performance of the process variables of the NSSS modules as well as the secondary-loop FFN.The variation of a NSSS power-level setpoint directly induces the variations of the setpoints of nuclear power, reactor outlet helium temperature, primary helium flowrate, and secondary feedwater flowrate of this NSSS, which enlarge the errors between the practical values and setpoints of these variables. These errors drive the nuclear power controller and helium temperature controller to insert or withdraw the control rods and drive the flowrate controllers to slow down or speed up the helium blower and feedwater pump. Due to the negative feedback laws, the movement of control rods, helium blower, and feedwater pump can suppress the errors. Moreover, the variations of the flowrates and helium temperature induce the variations of steam temperature and main steam pressure, which drive the steam temperature controller to revise the setpoints of helium or feedwater flowrates for maintaining the steam temperature and drive the pressure controller to change the steam valve opening for keeping the main steam pressure. The variation of the steam valve opening can also reduce the feedwater flowrate of the other NSSS, which leads to the successive variations of steam temperature, helium temperature, fuel temperature, and nuclear power of the other NSSS. Then, the controllers of the other NSSS are active to suppress the errors of process variables caused by the variation of steam valve opening. Furthermore, since the steam flow extracted from the turbine is equally distributed to the two HPHs, the power maneuver of one module changes feedwater temperatures of both two modules. To maintain the thermal powers, the setpoints of feedwater flowrates of the two modules should be adjusted by the thermal power controllers.Specifically to the NSSS module control performance, we can clearly see from the simulations results illustrated in Figures8–15 that simple module control laws (22)–(26) can guarantee satisfactory stability and transient performance for the NSSS modules in the cases of power-level maneuver and maintenance. Here, the expressions of these module control laws are given by guaranteeing proper dissipation feature for the NSSS module with its shifted-ectropy as the storage function, which is the key idea of physics-based nuclear plant control method with the ability of giving simple control laws. These simulation results effectively show the feasibility of these module control laws which can be implemented and commissioned easily in the practical engineering.Furthermore, the key difference between the HTR-PM plant and those single-modular nuclear plants is the module coordination control. There is a tight hydraulic coupling effect of the two NSSS modules that is induced by the common secondary-loop FFN. There is also coupling effect between the feedwater temperatures of the two NSSS modules, which is induced by the equal distribution of the steam flow extracted from the turbine to the two HPHs belonging to the two modules, respectively. Because of the existence of the hydraulic and thermal module coupling effects, when one NSSS module is in the state of power maneuvering, both the feedwater flowrate and output thermal power of the other NSSS module must vary simultaneously if there is no control input. Due to the proper referenced feedwater flowrates given by thermal power controllers, and due to the effectiveness of the feedwater flowrate control law, both the secondary flowrates and thermal power of the two NSSS modules are statically decoupled with each other, which means that the steady values of the feedwater flowrates (outlet thermal power) of the two NSSS modules are decoupled with each other. For example, from columns 1 and 3 of Figures8 and 10, the power step of 1# NSSS module does not influence the steady values of the feedwater flowrate and outlet thermal power of 2# NSSS module. Thanks to the plant control strategy proposed in this paper, we can see from the simulation results that the power step and ramp of one NSSS module only induce small transient disturbance to the other module and do not induce steady-point perturbations of the secondary flowrate and thermal power. The feedback laws of the module coordinated control given by (24) and (27) are also simple for engineering implementation and conditioning.In summary, the power maneuver of the HTR-PM plant can be realized by the plant control strategy given in Section3 which has already been applied to the HTR-PM power control project. The numerical simulation results show that the HTR-PM plant is flexible enough for constituting HESs with those renewable energy sources such as the wind and solar. Moreover, the technical issues of integrating more NSSS modules are nearly the same with integrating two modules by allocating a pump to each NSSS module, implementing module control laws in Section 3.2.2 to each NSSS module and applying module coordination laws in Section 3.2.3 to the secondary FFN.


## 4.1. Simulation Results

### 4.1.1. Power Step

The dynamic responses of the power step between 100% and 95% plant full power (PFP) are obtained through numerical simulation.(
1) Power Step-Down. Initially, the HTR-PM plant remains at the steady state of 100% PFP which means that the two NSSS modules are running at the 100% reactor full power (RFP). At 5000 s, the thermal power setpoint of the first (i.e., 1#) NSSS module steps down from 100% to 95% RFP, and then at 10000 s, the referenced thermal power of the second (i.e., 2#) NSSS module steps down from 100% to 95% RFP. The dynamic responses of the relative nuclear powers, reactor outlet helium temperatures, OTSG outlet steam temperatures, and thermal powers of 1# and 2# NSSS modules are given in Figure 8. The responses of flowrates Qi (i=1,2,3), feedwater pressures -Hfk (k=1,2), and main steam pressure H3 are illustrated in Figure 9.Figure 8
Responses of 1# and 2# NSSS in the case of power step-down:nri, Tcouti, Tsi, and Pthi are, respectively, the relative nuclear power, MHTGR outlet helium temperature, OTSG outlet steam temperature, and total thermal power of the ith NSSS, i=1,2.Figure 9
Responses of secondary-loop FFN flowrates and pressure drops in the case of power step-down.(
2) Power Step-Up. Initially, the HTR-PM plant remains at the steady state of 95% PFP; that is, the two NSSS modules operate at 95% RFP. The thermal power setpoints of 1# and 2# NSSS modules step up from 95% to 100% RFP at 5000 s and 10000 s, respectively. The dynamic responses of the relative nuclear powers, reactor outlet helium temperatures, OTSG outlet steam temperatures, and thermal powers of 1# and 2# NSSS modules are given in Figure 10. The responses of Qi (i=1,2,3), -Hfk (k=1,2), andH3 are shown in Figure 11.Figure 10
Responses of 1# and 2# NSSS in the case of power step-up:nri, Tcouti, Tsi, and Pthi are, respectively, the relative nuclear power, MHTGR outlet helium temperature, OTSG outlet steam temperature, and total thermal power of the ith NSSS, i=1,2.Figure 11
Responses of secondary-loop FFN flowrates and pressure drops in the case of power step-up.


### 4.1.2. Power Ramp

The responses of the steady power ramp between 100% and 50% plant full power (PFP) are given by the simulation study.(
1) Power Ramp-Down. Initially, the HTR-PM plant remains at the steady state of 100% PFP; that is, the two NSSS modules operate at 100% RFP. At 5000 s, 1# NSSS module maneuvers from 100% to 50% RFP steadily in 10 minutes, and then at 10000 s, the referenced thermal power of 2# NSSS module changes steadily from 100% to 50% RFP in 10 minutes, which is just the power decrease process from 100% to 50% PFP of the HTR-PM plant. The responses of the relative nuclear powers, reactor outlet helium temperatures, OTSG outlet steam temperatures, and thermal powers of 1# and 2# NSSS modules are given in Figure 12. The responses of flowrates Qi (i=1,2,3), feedwater headers -Hfk (k=1,2), and main steam pressureH3 are illustrated in Figure 13.Figure 12
Responses of 1# and 2# NSSS in the case of power ramp-down:nri, Tcouti, Tsi, and Pthi are, respectively, the relative nuclear power, MHTGR outlet helium temperature, OTSG outlet steam temperature, and total thermal power of the ith NSSS, i=1,2.Figure 13
Responses of secondary-loop FFN flowrates and pressure drops in the case of power ramp-down.(
2) Power Ramp-Up. Initially, the HTR-PM plant remains at the steady state of 95% PFP; that is, the two NSSS modules operate at 95% RFP. At 5000 s and 10000 s, 1# and 2# NSSS modules, respectively, maneuver steadily from 50% to 100% RFP in 10 minutes, which is the process of HTR-PM power increase from 50% to 100% PFP. The responses of the relative nuclear powers, reactor outlet helium temperatures, OTSG outlet steam temperatures, and thermal powers of 1# and 2# NSSS modules are given in Figure 14. The responses of flowrates Qi (i=1,2,3), -Hfk (k=1,2), andH3 are illustrated in Figure 15.Figure 14
Responses of 1# and 2# NSSS in the case of power ramp-up:nri, Tcouti, Tsi, and Pthi are, respectively, the relative nuclear power, MHTGR outlet helium temperature, OTSG outlet steam temperature, and total thermal power of the ith NSSS, i=1,2.Figure 15
Responses of secondary-loop FFN flowrates and pressure drops in the case of power ramp-up.


## 4.1.1. Power Step

The dynamic responses of the power step between 100% and 95% plant full power (PFP) are obtained through numerical simulation.(
1) Power Step-Down. Initially, the HTR-PM plant remains at the steady state of 100% PFP which means that the two NSSS modules are running at the 100% reactor full power (RFP). At 5000 s, the thermal power setpoint of the first (i.e., 1#) NSSS module steps down from 100% to 95% RFP, and then at 10000 s, the referenced thermal power of the second (i.e., 2#) NSSS module steps down from 100% to 95% RFP. The dynamic responses of the relative nuclear powers, reactor outlet helium temperatures, OTSG outlet steam temperatures, and thermal powers of 1# and 2# NSSS modules are given in Figure 8. The responses of flowrates Qi (i=1,2,3), feedwater pressures -Hfk (k=1,2), and main steam pressure H3 are illustrated in Figure 9.Figure 8
Responses of 1# and 2# NSSS in the case of power step-down:nri, Tcouti, Tsi, and Pthi are, respectively, the relative nuclear power, MHTGR outlet helium temperature, OTSG outlet steam temperature, and total thermal power of the ith NSSS, i=1,2.Figure 9
Responses of secondary-loop FFN flowrates and pressure drops in the case of power step-down.(
2) Power Step-Up. Initially, the HTR-PM plant remains at the steady state of 95% PFP; that is, the two NSSS modules operate at 95% RFP. The thermal power setpoints of 1# and 2# NSSS modules step up from 95% to 100% RFP at 5000 s and 10000 s, respectively. The dynamic responses of the relative nuclear powers, reactor outlet helium temperatures, OTSG outlet steam temperatures, and thermal powers of 1# and 2# NSSS modules are given in Figure 10. The responses of Qi (i=1,2,3), -Hfk (k=1,2), andH3 are shown in Figure 11.Figure 10
Responses of 1# and 2# NSSS in the case of power step-up:nri, Tcouti, Tsi, and Pthi are, respectively, the relative nuclear power, MHTGR outlet helium temperature, OTSG outlet steam temperature, and total thermal power of the ith NSSS, i=1,2.Figure 11
Responses of secondary-loop FFN flowrates and pressure drops in the case of power step-up.


## 4.1.2. Power Ramp

The responses of the steady power ramp between 100% and 50% plant full power (PFP) are given by the simulation study.(
1) Power Ramp-Down. Initially, the HTR-PM plant remains at the steady state of 100% PFP; that is, the two NSSS modules operate at 100% RFP. At 5000 s, 1# NSSS module maneuvers from 100% to 50% RFP steadily in 10 minutes, and then at 10000 s, the referenced thermal power of 2# NSSS module changes steadily from 100% to 50% RFP in 10 minutes, which is just the power decrease process from 100% to 50% PFP of the HTR-PM plant. The responses of the relative nuclear powers, reactor outlet helium temperatures, OTSG outlet steam temperatures, and thermal powers of 1# and 2# NSSS modules are given in Figure 12. The responses of flowrates Qi (i=1,2,3), feedwater headers -Hfk (k=1,2), and main steam pressureH3 are illustrated in Figure 13.Figure 12
Responses of 1# and 2# NSSS in the case of power ramp-down:nri, Tcouti, Tsi, and Pthi are, respectively, the relative nuclear power, MHTGR outlet helium temperature, OTSG outlet steam temperature, and total thermal power of the ith NSSS, i=1,2.Figure 13
Responses of secondary-loop FFN flowrates and pressure drops in the case of power ramp-down.(
2) Power Ramp-Up. Initially, the HTR-PM plant remains at the steady state of 95% PFP; that is, the two NSSS modules operate at 95% RFP. At 5000 s and 10000 s, 1# and 2# NSSS modules, respectively, maneuver steadily from 50% to 100% RFP in 10 minutes, which is the process of HTR-PM power increase from 50% to 100% PFP. The responses of the relative nuclear powers, reactor outlet helium temperatures, OTSG outlet steam temperatures, and thermal powers of 1# and 2# NSSS modules are given in Figure 14. The responses of flowrates Qi (i=1,2,3), -Hfk (k=1,2), andH3 are illustrated in Figure 15.Figure 14
Responses of 1# and 2# NSSS in the case of power ramp-up:nri, Tcouti, Tsi, and Pthi are, respectively, the relative nuclear power, MHTGR outlet helium temperature, OTSG outlet steam temperature, and total thermal power of the ith NSSS, i=1,2.Figure 15
Responses of secondary-loop FFN flowrates and pressure drops in the case of power ramp-up.


## 4.2. Discussions

The variation of the plant power is realized by the successive power maneuver of 1# and 2# NSSS modules, and at most one NSSS module can be in the state of power maneuver at any time. Generally, from Figures8–15, the control laws for the modules and those for the module coordination guarantee satisfactory closed-loop stability and transient performance of the process variables of the NSSS modules as well as the secondary-loop FFN.The variation of a NSSS power-level setpoint directly induces the variations of the setpoints of nuclear power, reactor outlet helium temperature, primary helium flowrate, and secondary feedwater flowrate of this NSSS, which enlarge the errors between the practical values and setpoints of these variables. These errors drive the nuclear power controller and helium temperature controller to insert or withdraw the control rods and drive the flowrate controllers to slow down or speed up the helium blower and feedwater pump. Due to the negative feedback laws, the movement of control rods, helium blower, and feedwater pump can suppress the errors. Moreover, the variations of the flowrates and helium temperature induce the variations of steam temperature and main steam pressure, which drive the steam temperature controller to revise the setpoints of helium or feedwater flowrates for maintaining the steam temperature and drive the pressure controller to change the steam valve opening for keeping the main steam pressure. The variation of the steam valve opening can also reduce the feedwater flowrate of the other NSSS, which leads to the successive variations of steam temperature, helium temperature, fuel temperature, and nuclear power of the other NSSS. Then, the controllers of the other NSSS are active to suppress the errors of process variables caused by the variation of steam valve opening. Furthermore, since the steam flow extracted from the turbine is equally distributed to the two HPHs, the power maneuver of one module changes feedwater temperatures of both two modules. To maintain the thermal powers, the setpoints of feedwater flowrates of the two modules should be adjusted by the thermal power controllers.Specifically to the NSSS module control performance, we can clearly see from the simulations results illustrated in Figures8–15 that simple module control laws (22)–(26) can guarantee satisfactory stability and transient performance for the NSSS modules in the cases of power-level maneuver and maintenance. Here, the expressions of these module control laws are given by guaranteeing proper dissipation feature for the NSSS module with its shifted-ectropy as the storage function, which is the key idea of physics-based nuclear plant control method with the ability of giving simple control laws. These simulation results effectively show the feasibility of these module control laws which can be implemented and commissioned easily in the practical engineering.Furthermore, the key difference between the HTR-PM plant and those single-modular nuclear plants is the module coordination control. There is a tight hydraulic coupling effect of the two NSSS modules that is induced by the common secondary-loop FFN. There is also coupling effect between the feedwater temperatures of the two NSSS modules, which is induced by the equal distribution of the steam flow extracted from the turbine to the two HPHs belonging to the two modules, respectively. Because of the existence of the hydraulic and thermal module coupling effects, when one NSSS module is in the state of power maneuvering, both the feedwater flowrate and output thermal power of the other NSSS module must vary simultaneously if there is no control input. Due to the proper referenced feedwater flowrates given by thermal power controllers, and due to the effectiveness of the feedwater flowrate control law, both the secondary flowrates and thermal power of the two NSSS modules are statically decoupled with each other, which means that the steady values of the feedwater flowrates (outlet thermal power) of the two NSSS modules are decoupled with each other. For example, from columns 1 and 3 of Figures8 and 10, the power step of 1# NSSS module does not influence the steady values of the feedwater flowrate and outlet thermal power of 2# NSSS module. Thanks to the plant control strategy proposed in this paper, we can see from the simulation results that the power step and ramp of one NSSS module only induce small transient disturbance to the other module and do not induce steady-point perturbations of the secondary flowrate and thermal power. The feedback laws of the module coordinated control given by (24) and (27) are also simple for engineering implementation and conditioning.In summary, the power maneuver of the HTR-PM plant can be realized by the plant control strategy given in Section3 which has already been applied to the HTR-PM power control project. The numerical simulation results show that the HTR-PM plant is flexible enough for constituting HESs with those renewable energy sources such as the wind and solar. Moreover, the technical issues of integrating more NSSS modules are nearly the same with integrating two modules by allocating a pump to each NSSS module, implementing module control laws in Section 3.2.2 to each NSSS module and applying module coordination laws in Section 3.2.3 to the secondary FFN.


## 5. Conclusions

Due to its inherent safety feature and potential economic competitiveness, the modular high temperature gas-cooled reactor (MHTGR) has already been seen as one of the best candidates in building the next generation nuclear plants. Through the multimodular scheme, the inherent safety feature of the MHTGR can be applicable to the large-scale nuclear plants at any desired power ratings. Plant power control is the one of the crucial techniques for guaranteeing the safe, stable, and efficient operation of every nuclear plant, and the plant power control strategy of the traditional single-modular nuclear plants cannot be directly applied to multimodular plants. Thus, it is very necessary to study the plant power control of multimodular nuclear plants. Here, dynamic modeling, control strategy design, and performance verification are the three main aspects of plant power control studying, where the former one is the basis for the latter two. In this paper, the dynamical modeling of the main parts of the HTR-PM plant is first given, and the corresponding plant control strategy including the feedback loops and control laws is then proposed. By integrating the dynamic model and power control strategy for a simulation code in MATLAB/Simulink environment, numerical verification is performed. Simulation results in the cases of plant power step and ramp show that the plant control strategy provides satisfactory closed-loop stability and transient performance for the HTR-PM plant. The future work mainly lies in commissioning the power control system of the HTR-PM plant, which is implemented based on the control strategy given in this paper, practically for a practical verification.



---

*Source: 6298037-2017-05-22.xml*