<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JPEE</journal-id><journal-title-group><journal-title>Journal of Power and Energy Engineering</journal-title></journal-title-group><issn pub-type="epub">2327-588X</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jpee.2016.47003</article-id><article-id pub-id-type="publisher-id">JPEE-69355</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Engineering</subject></subj-group></article-categories><title-group><article-title>
 
 
  Dynamic Characteristics Analysis on MHTGR Plant’s Secondary Side Fluid Flow Network
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Maoxuan</surname><given-names>Song</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Zhe</surname><given-names>Dong</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Institute of Nuclear and New Energy Technology, Collaborative Innovation Center of Advanced Nuclear Energy Technology of China, Key Laboratory of Advanced Reactor Engineering and Safety of Ministry of Education, Tsinghua University, Beijing, China</addr-line></aff><pub-date pub-type="epub"><day>25</day><month>07</month><year>2016</year></pub-date><volume>04</volume><issue>07</issue><fpage>15</fpage><lpage>22</lpage><history><date date-type="received"><day>30</day>	<month>March</month>	<year>2016</year></date><date date-type="rev-recd"><day>accepted</day>	<month>26</month>	<year>July</year>	</date><date date-type="accepted"><day>29</day>	<month>July</month>	<year>2016</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  Multipe NSSS (Nuclear Steam Supply System) modules use the common feeding-water system to drive the common turbine power generation set. The SSFFN (secondary side fluid flow network) of MHTGR plant has features i.e. strong-coupling and nonlinearity. A wide range of power switching operation will cause unsteady flow, which may destroy the working elements and will be a threat for normal operation. To overcome those problems, a differential-algebraic model and PI controllers are designed for the SSFFN. In MATLAB\SIMULINK environment, a simulation platform is established and used to make a simulation of SSFFN of a MHTGR plant with two NSSS modules, which uses feedwater valves to control the mass flow rate in each module instead of feedwater pump. Results reflect good robustness of controllers.
 
</p></abstract><kwd-group><kwd>MHTGR Plant</kwd><kwd> Secondary Side Fluid Flow Network</kwd><kwd> a Differential-Algebraic Model</kwd><kwd> PI Controllers</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>As the fourth generation of nuclear power system, the modular high-temperature gas-cooled reactor (MHTGR) plant receives extensive attention from the world’s nuclear power industry due to its inherent safety and economic efficiency. As characterized by modularization and miniaturization, MHTGR uses helium as coolant and uses graphite as moderator and structural materials. Its fuel elements contain thousands of ceramic mould package particles. Its core outlet helium temperature can reach 700˚C - 950˚C [<xref ref-type="bibr" rid="scirp.69355-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.69355-ref2">2</xref>]. In order to own inherent safety, the core density and geometry size is designed smaller, thus the highest temperature of fuel elements will not surpass the limit temperature 1620˚C [<xref ref-type="bibr" rid="scirp.69355-ref3">3</xref>]. In order to improve economic competitiveness, MHTGR plant adopts modularization construction, using multiple NSSS (Nuclear Steam Supply System) modules to feed one steam turbine generator set, e.g. HTR-PM plant. HTR-PM plant uses two NSSS modules to feed one steam turbine set. Each NSSS module contains a pebble-bed reactor, whose heat power is 250 MW<sub>th</sub>, a helium fan, a once-through steam generator and corresponding pipes. Steam generated by each NSSS flows into the same steam header and then flows into the steam turbine generator set after uniform mixed in the steam header. Because the two modules use the same feed-water system to feed the same steam turbine generator set, the secondary side fluid flow network (SSFFN) between the two modules has features e.g. strong-coupling, nonlinearity, etc. A wide range of power switching operation may cause unsteady fluid flow, which may destroy working components. In order to control modules coordinately and ensure MHTGR plant’s safety, it is meaningful to make a research on the dynamic characteristics of the SSFFN.</p><p>Fluid Flow networks (FFNs) widely exist in industrial engineering, e.g. mine ventilation, water supply and drainage system and oil gas industry, etc. [<xref ref-type="bibr" rid="scirp.69355-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.69355-ref5">5</xref>]. Research on FFNs began in study on controlling mine ventilation system in the 1970s-1980s. It is based on the fluid dynamics and can be divided into two categories, active control and passive control methods. The difference between the two is that active control method uses the fluid resistance characteristics and relies on auxiliary control operations [<xref ref-type="bibr" rid="scirp.69355-ref6">6</xref>] [<xref ref-type="bibr" rid="scirp.69355-ref7">7</xref>]. Kocic [<xref ref-type="bibr" rid="scirp.69355-ref8">8</xref>] pioneered the research on FFNs. He approximately regarded the mine ventilation system as a lumped parameter dynamic model and designed a linear feedback flow controller by using branch fluid dynamics. However, lacking of graph properties, he ignored nonlinear features of FFNs. By incorporating graph properties and branch fluid dynamics, Krstić [<xref ref-type="bibr" rid="scirp.69355-ref9">9</xref>] put forward the nonlinear non-minimal model and minimal model. Both models can guarantee the closed-loop global exponential stability of the FFNs. Dong Zhe [<xref ref-type="bibr" rid="scirp.69355-ref10">10</xref>] proposed distributed control method based on the minimal model [<xref ref-type="bibr" rid="scirp.69355-ref9">9</xref>]. The method is formed by the flow controller on the network of independent branch. He realized distributed control of FFNs and guaranteed the global asymptotic stability of the closed-loop system.</p><p>In this paper, a differential-algebraic model of MHGTR plant’s SSFFN, based on Kirchhoff’s law and branch dynamics, and PI controller sensuring globally asymptotic stability are given. What’s more, a simulation platform was established, which can be used for simulation experiments on SSFFN, which has more NSSS modules.</p></sec><sec id="s2"><title>2. Model of SSFFN</title><sec id="s2_1"><title>2.1. Structure of the SSFFN</title><p>From <xref ref-type="fig" rid="fig1">Figure 1</xref>, We can see two type plants with similar structure. NSSS modules in left are fed by separate feedwater pumps. The SSFFN comprises a steam heater, a steam turbine generator set, condenser, LP heater, HP heater, deaerator, pipes and valves etc. The mass flow rate in each NSSS module is controlled by feedwater pumps and flow changes in one NSSS module doesn’t influence other modules’. The plant in right uses the same feedwater pump to feed water for all NSSS modules and uses valves instead of pumps to control the mass flow rate in NSSS modules. Like SSFFN in left, the flow’s change in one NSSS module doesn’t influence other modules’, but it does impact on other pipes’ flow in the SSFFN. Thus, control on SSFFN is a multivariable control problem as the SSFFN has strong-coupling and nonlinear features. This paper will make analysis on SSFFN in right.</p></sec><sec id="s2_2"><title>2.2. Branch Dynamic Equation</title><p>Assuming that all flow dynamic parameters in pipes are evenly distributed on the radial cross section. Under the assumption, the flow in pipes can be seen as a one-dimensional flow. We can use the average values of dynamic parameters of the radial cross section and branch dynamic equation is deduced as Equation (1).</p><disp-formula id="scirp.69355-formula404"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/69355x4.png"  xlink:type="simple"/></disp-formula><p>where Q<sub>j</sub>, R<sub>j</sub>, H<sub>j</sub> and K<sub>j</sub> is the mass flow rate, the resistance, the pressure drop, the inertia coefficient of the branch j, j = 1, ∙∙∙, N and N is the branch number of network.</p></sec><sec id="s2_3"><title>2.3. Graph Features</title><p>Making a topology analysis on the structure of the SSFFN of MHTGR plant shown in right in <xref ref-type="fig" rid="fig1">Figure 1</xref>.</p><p>There are n + 3 nodes and 2 &#215; n + 2 branches (including 1 pump branch) in this network shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>. The SSFFN can be divided into a tree containing the feedwater pump branch and its complement, the co-tree whose branches are referred to as links. From <xref ref-type="fig" rid="fig2">Figure 2</xref>, it is clear that the number of the link is n. In <xref ref-type="fig" rid="fig2">Figure 2</xref>, branches Q<sub>1</sub>-Q<sub>n</sub> representing 1#-n# NSSS module are links and the rest are trees. Branch Q<sub>n</sub><sub>+1</sub> represents steam turbine generator set. Branch Q<sub>f</sub> represents feedwater pump branch which is connected to the node 1. BranchesQ<sub>v</sub><sub>1</sub>-Q<sub>v</sub><sub>n</sub> represents feeding valve branches. Like an electrical network, the SSFFN must satisfy Kirchhoff’s current law (KCL) and voltage law (KVL) [<xref ref-type="bibr" rid="scirp.69355-ref9">9</xref>]. Based on Kirchhoff’s current law, we have</p><fig-group id="fig1"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> MHTGR plant.</title></caption><fig id ="fig1_1"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/69355x5.png"/></fig><fig id ="fig1_2"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/69355x6.png"/></fig></fig-group><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Topology of the SSFFN</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/69355x7.png"/></fig><disp-formula id="scirp.69355-formula405"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/69355x8.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.69355-formula406"><graphic  xlink:href="http://html.scirp.org/file/69355x9.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.69355-formula407"><graphic  xlink:href="http://html.scirp.org/file/69355x10.png"  xlink:type="simple"/></disp-formula><p>Based on Kirchhoff’s voltage law, we have</p><disp-formula id="scirp.69355-formula408"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/69355x11.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.69355-formula409"><graphic  xlink:href="http://html.scirp.org/file/69355x12.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.69355-formula410"><graphic  xlink:href="http://html.scirp.org/file/69355x13.png"  xlink:type="simple"/></disp-formula><p>H<sub>f</sub> is the pressure drop of the feedwater pump branch. H<sub>f</sub> can be expressed as<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/69355x14.png" xlink:type="simple"/></inline-formula>. So we can obtain</p><disp-formula id="scirp.69355-formula411"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/69355x15.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.69355-formula412"><graphic  xlink:href="http://html.scirp.org/file/69355x16.png"  xlink:type="simple"/></disp-formula><p>From Equation (4), we can obtain the mass flow rates and pressure drops of the branch in the SSFFN, it is the differential-algebraic model of the SSFFN. If modules’ number n is given, it is clear that the topology of the SSFFN can be determined. So modules’ number n is important to determine the SSFFN dynamics.</p></sec></sec><sec id="s3"><title>3. Control Scheme</title><p>In actual engineering, the mass flow rates and pressure drops should reach the expectation values. From Equation (4), a differential-algebraic model of the SSFFN has been established. R, K and H<sub>d</sub> is the input of the model. We can adjust the opening of feedwater valves to change mass flow rates in the NSSS modules and control the main steam pressure drop of the steam turbine generator set by adjusting the pressuer head of feedwater pump. A PI control scheme based on above has been established and its structure can be seen in <xref ref-type="fig" rid="fig3">Figure 3</xref>. We can use controller H<sub>d</sub> to control main pressure drop H<sub>n</sub><sub>+1</sub> of the steam turbine generator set and controllers R<sub>vi</sub> (i = 1,∙∙∙, n) to control mass flow rates Q<sub>i</sub> (i = 1,∙∙∙, n) in NSSS modules.</p><p>To overcome the pipe sensor noise, a section of inertia link is added in controllers, which plays an role of low pass filter. Controllers shown in <xref ref-type="fig" rid="fig3">Figure 3</xref> can guarantee the global asymptotic stability of the system, which can can be proved by Lyapunov energy function. In addition, the simulation results in this paper can also verify the above control scheme.</p></sec><sec id="s4"><title>4. Implementation and Simulation</title><sec id="s4_1"><title>4.1. Implementation</title><p>In MATLAB\SIMULINK environment, the model has been implemented by using the programming language M code based on the principles of Level-2 S-Function. The model provides parameters interface in order to let users to define the number of NSSS modules and structure of the SSFFN. The users must input the parameters, such as the nodes’ number, the feeding-pump branches’ number, the rest branches’ number and the matrices of KCL, which depends on the number of NSSS modules, to initialize the model. The users can add other simulation modules to make real-time simulation.</p></sec><sec id="s4_2"><title>4.2. Simulation</title><p>In MATLAB\SIMULINK environment, the model has been implemented by using the programming language M code based on the principles of Level-2 S-Function (<xref ref-type="fig" rid="fig4">Figure 4</xref>). The model provides parameters interface in order to let users to define the number of NSSS modules and structure of the SSFFN. The users must input the parameters, such as the nodes’ number, the feeding-pump branches’ number, the rest branches’ number and the matrices of KCL, which depends on the number of NSSS modules, to initialize the model. The users can add other simulation modules to make real-time simulation.</p><p>This paper has made simulation of SSFFN of MHTGR plant, which has two NSSS modules and uses the common feedwater pump to feed for the two. In this network, the mass flow rate is controlled by the feeding valve, making topology analysis on the structure of the SSFFN and setting inertia coefficients of the branch excluding feeding-pump branch K as diag (0.1; 0.1; 0.1; 0.1; 0.1). To simulate the real operation ,conditions can be set as follow, the time of simulation t is 3000s,the initial value of pressure head of feedwater pump H<sub>d</sub> is 25 Mpa, the resistance of branch R<sub>1</sub>, R<sub>2</sub>, R<sub>3</sub>, R<sub>f</sub> is 0.00118, 0.00118, 0.0038, 0.00044, the initial value of R<sub>v</sub><sub>1</sub>, R<sub>v</sub><sub>2</sub> is</p><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> PI controllers</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/69355x17.png"/></fig><fig-group id="fig4"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> Encapsulation model and its parameters interface.</title></caption><fig id ="fig4_1"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/69355x18.png"/></fig><fig id ="fig4_2"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/69355x19.png"/></fig></fig-group><p>0.00118, 0.00118. The main steam pressure drop of the steam turbine generator set remains 13.9 Mpa. A case study is given to show the performance of the model. In this case, the mass flow rate of 1#NSSS remains 96 kg/s in 1000 s, it will drop to 90 kg/s at 1000 s. Controller H<sub>d</sub>: k<sub>p</sub> = 1, k<sub>i</sub> = 1; Controller R<sub>v</sub><sub>1</sub>: k<sub>p</sub> = 0.0001, k<sub>i</sub> = 0.0001; The mass flow rates, pressure drops and control inputs are respectively shown in <xref ref-type="fig" rid="fig5">Figure 5</xref>.</p><p>The mass flow rate of 1#NSSS module, Q<sub>1</sub>, is controlled by the feeding valves. The drop of mass flow rate is caused by the decrease of the valve opening, thus resulting in increase in the branches’ resistance and pressure</p><p>head of feedwater pump. Because of the controller, the main steam pressure drop, H<sub>3</sub> remians the initial value. The results show good robustness of controllers.</p></sec></sec><sec id="s5"><title>5. Conclusion</title><p>The SSFFN of MHTGR plant has features, i.e. strong-coupling and nonlinear, because multiple NSSS modules use the common feedwater system to drive the same steam turbine generator set. A wide range of power switching operation may cause unsteady fluid flow in the network, which may destroy working components. This problem must be avoided in order to ensure the system’s safety. It is meaningful to make analysis on dynamics of the SSFFN. Branch dynamics and network’s structure features plays an important role in modeling SSFFN.</p><fig-group id="fig5"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Case results.</title></caption><fig id ="fig5_1"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/69355x20.png"/></fig><fig id ="fig5_2"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/69355x21.png"/></fig><fig id ="fig5_3"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/69355x22.png"/></fig><fig id ="fig5_4"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/69355x23.png"/></fig><fig id ="fig5_5"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/69355x24.png"/></fig></fig-group><p>Base on KCL and KVL, this paper designed a differential-algebraic model for SSFFN. Based on the model, this paper designed PI controllers to promise the closed system globally asymptotically stable. In MATLAB\ SIMULINK environment, the paper made implementation of the model and designed a good human-machine interface for users to define the structure of the SSFFN at ease. To proof the model and the controller, this paper made simulation of SSFFN of a MHTGR plant which has two modules and uses one common feedwater pump feeding for the two. The results efficiently reflect the controllers’ good performance. The simulation platform established in this paper may help engineers to make integral coefficient test and choose the best one, which is actually based on engineers’ experiences. Also, the simulation platform can be used to control SSFFN of MHTGR plant which has more modules.</p></sec><sec id="s6"><title>Acknowledgements</title><p>This work was supported in part by National S&amp;T Major Project under Grant ZX06901, and in part by Natural Science Foundation of China (NSFC) under Grant 61374045.</p></sec><sec id="s7"><title>Cite this paper</title><p>Maoxuan Song,Zhe Dong, (2016) Dynamic Characteristics Analysis on MHTGR Plant’s Secondary Side Fluid Flow Network. Journal of Power and Energy Engineering,04,15-22. doi: 10.4236/jpee.2016.47003</p></sec></body><back><ref-list><title>References</title><ref id="scirp.69355-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Zhang, Z., Wu, Z., Wang, D., et al. (2009) Current Status and Technical Description of Chinese 2 × 250 MWth HTR- PM Demonstration Plant. Nuclear Engineering and Design, 239, 1212-1219.  
http://dx.doi.org/10.1016/j.nucengdes.2009.02.023</mixed-citation></ref><ref id="scirp.69355-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Lohnert, G.H. (1990) Technical Design Features and Essential Safety-Related Properties of the HTR-Module. Nuclear Engineering and Design, 121, 259-275. http://dx.doi.org/10.1016/0029-5493(90)90111-A</mixed-citation></ref><ref id="scirp.69355-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Economic Potential of Modular Reactor Nuclear Power Plants Based on the Chinese HTR-PM Project</mixed-citation></ref><ref id="scirp.69355-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Hu, Y., Koroleva, O.I. and Krsti?, M. (2003) Nonlinear Control of Mine Ventilation Networks. Systems &amp; control letters, 49, 239-254. http://dx.doi.org/10.1016/S0167-6911(02)00336-5 </mixed-citation></ref><ref id="scirp.69355-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Izadi, G. and Elsworth, D. (2014) Reservoir Stimulation and Induced Seismicity: Roles of Fluid Pressure and Thermal Transients on Reactivated Fractured Networks. Geothermics, 51, 368-379.  
http://dx.doi.org/10.1016/j.geothermics.2014.01.014</mixed-citation></ref><ref id="scirp.69355-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Aldridge, M.D., Swartwout, R.E., Smith Jr., N.S., et al. (1976) Electronic Monitoring and Control of Mine Ventilation. Proceedings of the 3rd WVU Conference Coal Mine Electrotechnology, Morgantown.</mixed-citation></ref><ref id="scirp.69355-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Hu, Y., Koroleva, O.I. and Krstic, M. (2002) Control Design for Mine Ventilation Network Systems. Proceedings of the 41st IEEE Conference on Decision and Control, 1, 543-548.</mixed-citation></ref><ref id="scirp.69355-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Kocic, D.D. (1979) On the Autonomy of Local Systems in Mine Ventilation Control. 2nd Mine Ventilation Congress, Reno.</mixed-citation></ref><ref id="scirp.69355-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Hu, Y., Koroleva, O.I. and Krsti?, M. (2003) Nonlinear Control of Mine Ventilation Networks. Systems &amp; control letters, 49, 239-254. http://dx.doi.org/10.1016/S0167-6911(02)00336-5 </mixed-citation></ref><ref id="scirp.69355-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Dong, Z. (2015) Dissipation Analysis and Adaptive Control of Fluid Networks. The 2015 American Control Conference, Chicago.</mixed-citation></ref></ref-list></back></article>