Integrated Robust and Resilient Control of Nuclear Power Plants for Operational Safety and High Performance
IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 57, NO. 2, APRIL 2010 807Integrated Robust and Resilient Control of NuclearPower Plants for Operational Safety andHigh PerformanceXin Jin, Student Member, IEEE, Asok Ray, Fellow, IEEE, and Robert M. EdwardsAbstractThis paper presents an integrated robust and resilientcontrol strategy to enhance the operational safety and performanceof nuclear power plants. The objective of robust control is to min-imize the sensitivity of plant operations to exogenous disturbancesand internal faults while achieving a guaranteed level of perfor-mance with a priori specified bounds of uncertainties. On the otherhand, the role of resilient control is to enhance plant recovery fromunanticipated adverse conditions and faults as well as from emer-gency situations by altering its operational envelope in real time. Inthis paper, the issues of real-time resilient control of nuclear powerplants are addressed for fast response during emergency opera-tions while the features of the existing robust control technologyare retained during normal operations under both steady-state andtransient conditions. The proposed control methodology has beenvalidated on the International Reactor Innovative & Secure (IRIS)simulator of nuclear power plants.Index TermsEmergency operation, nuclear power plant, oper-ational safety, resilient control, robust control.ACRONYMSBIBO Bounded-input bounded-output.FSM Finite state machine.IRIS International reactor innovative & secure.LFT Linear fractional transformation.LOFA Loss-of-flow accident.MIMO Multi-input multi-output.RCP Reactor coolant pump.SISO Single-input single-output.NOMENCLATUREThe variables that are used in the controller design procedure,described in Sections II and IV-B, are listed below.Manuscript received November 12, 2009; revised January 24, 2010. Cur-rent version published April 14, 2010. This work was supported in part by theU.S. Department of Energy under NERI-C Grant DE-FG07-07ID14895 and byNASA under Cooperative Agreement NNX07AK49A. Any opinions, findingsand conclusions or recommendations expressed in this publication are those ofthe authors and do not necessarily reflect the views of the sponsoring agencies.The authors are with the Department of Mechanical and Nuclear Engineering,The Pennsylvania State University, University Park, PA, 16802 USA (e-mail:firstname.lastname@example.org; email@example.com; firstname.lastname@example.org).Color versions of one or more of the figures in this paper are available onlineat http://ieeexplore.ieee.org.Digital Object Identifier 10.1109/TNS.2010.2042071Low pass filter.Field of complex and real numbers.Disturbances and uncertainties.Tracking error.Lower LFT of plant P and controller K.Strictly proper transfer function.Robust controller.Minimum-phase stable transfer functionmatrix.Nominal plant.Augmented plant.Solution of algebraic Lyapunov equation.Reference Signals.Temperatures of the nuclear power plant.Sampling time.Controller outputs.Outputs of robust and resilient controllers.Disturbances.Sensor noise.Uncertainty input of .Control action weighting function.Sensor noise weighting function.Tracking error weighting function.Uncertainty in plant modeling.Measurement of plant output and its estimate.Output estimation errorDesired plant output following a desiredmodel.Weighted control action.Weighted tracking error.Uncertainty output of .-norm of a transfer matrix operator.0018-9499/$26.00 2010 IEEE808 IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 57, NO. 2, APRIL 2010State transition function of the FSM.Set of all possible uncertainty.Block structure of plant uncertainties.Block structure of performance objectives.Structured singular value.Vector of adaptive parameters and itsestimate.Time constants.I. INTRODUCTIONN UCLEAR power plants are complex dynamical systemswith many variables that require dynamical adjustmentsto achieve safety and efficiency over the entire operational enve-lope because their stability and performance could be severelylimited by a wide variety of safety requirements, operating con-ditions, internal faults, and exogenous disturbances. To achievethe specified goals, multiple control variables are simultane-ously manipulated for generating the required power and en-abling the plant to exploit alternate decision and control strate-gies. These strategies are often dictated by economy and safetyof plant operations. For example, feedback regulation of non-linear dynamics (e.g., control of reactor power and thermal hy-draulics in the balance of plant) could lead to static bifurcation,which is linked to degeneracy in the systems zero dynamics .When faced with unanticipated situations, such as equipmentfailures or large exogenous disturbances to the plant controlsystem, the human operators are required to carry out diag-nostic and corrective actions. Even experienced operators couldbe overwhelmed by the sheer number of display devices andsensor outputs to be monitored. An intelligent decision and con-trol system with a large degree of autonomy could enhance theoperational safety and performance of nuclear plants by allevi-ating the burden of human operators and simultaneously miti-gating the adverse consequences at an incipient stage.Several researchers have reported intelligent decision andcontrol methods (e.g., optimal control , fuzzy logic ,neural networks , and model predictive control ) to en-hance operational safety and performance of nuclear powerplants. Along this line, robust control techniques have also beeninvestigated , , where the role of a robust controller isto achieve disturbance rejection by reducing sensitivity of theplant control system (i.e., plant plus controller) to exogenousdisturbances and internal faults. The task is to synthesize arobust decision and control law on an infinite-time horizon by: Mitigation of the detrimental effects of uncertainties andexogenous disturbances; Trade-offs between plant stability and performance withinspecified bounds of uncertainties .However, stability and performance robustness of such con-trol algorithms may not be assured beyond the a priori speci-fied bounds of uncertainties and disturbances; usually larger arethe bounds, lower is the plant performance and plant instabilityis less likely. Therefore, the bounds of structured and unstruc-tured uncertainties are usually specified design parameters fortrade-off between stability and performance, and are often basedon nominal and off-nominal plant operations that may include atmost a few anticipated abnormalities. In the event of a plant acci-dent, the deviations from the nominal plant operating conditionsmay significantly exceed these uncertainty bounds. Hence, im-mediate actions beyond the regime of robust control are neededfor operational safety and subsequent restoration of normalcy tothe original operational mode or to a gracefully degraded mode.Guo et al.  have investigated resilient propulsion controlof aircraft to determine on how engine control systems can im-prove safe-landing probabilities under adverse conditions. Thekey idea is as follows: In emergency situations, the conserva-tive procedure of engine control may not be suitable for aircraftsafety; it may be advantageous to compromise the engine healthto save the aircraft. The aim of their research is to develop adap-tive engine control methodologies to operate the engine beyondthe normal domain for emergency operations to enhance safelanding at the expense of possibly partial damage in the aircraft.Motivated by this idea, resilient controllers can be designed fornuclear power plants to ensure the operational safety by sacri-ficing performance under adverse conditions. Recently, Hov-akimyan and coworkers  have reported a control algorithmcalled -adaptive control to address the issues in Integrated Re-silient Aircraft Control (IRAC) that is an active area of researchin National Aeronautics and Space Administration (NASA). Itis noted that the notion of resilience, introduced in the presentcontext, is entirely different from being non-fragile or insensi-tive to some errors in the nominal state-space matrices of con-troller during implementation .The role of the proposed resilient control in a nuclear powerplant is to enhance recovery of the control system from unan-ticipated adverse conditions and faults as well as from emer-gency situations by altering its operational envelope in real time. Resilient decision and control laws are synthesized on afinite-time horizon as augmentation of robust decision and con-trol with the objectives of: Reliable and fast recovery from adverse conditions andemergency situations; Restoration of the control configuration upon returning tonormalcy or upon graceful degradation within design spec-ifications.This paper addresses the issues of real-time resilient con-trol of nuclear power plants for fast response during emergencyoperations while the features of the existing robust controllerare retained during normal operating conditions. The goal is toformulate and validate resilient and reconfigurable control al-gorithms toward deployment of real-time controllers for emer-gency recovery of nuclear plants from expected and unexpectedadverse conditions.The integrated robust and resilient control strategy, developedin this paper, has been tested on the International Reactor Inno-vative & Secure (IRIS) simulator ,  that is built uponone of the next generation nuclear reactor designs for a mod-ular pressurized water reactor with an integral configuration.Currently, IRIS is in the stage of pre-application licensing withNuclear Regulatory Commission (NRC); its safety testing forJIN et al.: INTEGRATED ROBUST AND RESILIENT CONTROL OF NUCLEAR POWER PLANTS 809Design Certification (DC) is expected to be completed by 2010with deployment in the 20152017 time frame.The major contributions of the paper are delineated below: Introduction of the innovative concept of resilient controlin nuclear plant control for fast recovery from unantici-pated adverse conditions and emergency situations. Extension of the concept of robustness by integration withresilience for control of nuclear power plants under bothnormal operations and emergency situations. Validation of the concept of integrated robust and resilientcontrol of nuclear power plants on the IRIS simulator.The paper is organized in five sections including the presentone. Section II presents underlying principles of robust con-trol and resilient control. Section III addresses integration of ro-bust control and resilient control strategies. Section IV presentstesting and validation of these strategies on the IRIS simulator.The paper is summarized and concluded in Section V.II. ROBUST AND RESILIENT CONTROLThis section introduces the underlying principles of multi-variable robust control and resilient control for nuclear powerplants, where robust controllers are designed by the -based-synthesis method  and resilient controllers are designedby -adaptive output feedback algorithm . A finite statemachine is then used to integrate the above two controllers fornormal operation and emergency recovery from both expectedand unexpected adverse conditions while bumpless transferbetween the control modes is assured by usage of smoothingfilters.A. Robust Multivariable Control Using -SynthesisIn this paper, a -synthesis robust control approach is used tosynthesize the feedback controller. The plant uncertainties, in-cluding the effects of unmodeled dynamics, linearization, andmodel reduction, are characterized and estimated. Based on thespecified uncertainties, robust multivariable controllers are de-signed using D-K iteration  based on the stability and perfor-mance specifications. The order of the synthesized controllersis then reduced to an acceptable level by Hankel norm approxi-mation .1) Uncertainty Modeling: Uncertainties due to the inabilityto model relevant dynamics in an actual plant and the simplifi-cation to achieve a mathematical representation, including lin-earization and model reduction, are taken into consideration forsynthesis of robust controllers. Consequently, robustness of thesynthesized controller is dependent on the type and size of theuncertainties.In the controller design for nuclear power plants, Shyu andEdwards  considered three types of uncertainties due to un-modeled dynamics, linearization, and model reduction, respec-tively. These uncertainties contribute to the deviation betweenthe real plant and its reduced-order linear model, based on whichthe robust controller is synthesized. The unmodeled uncertain-ties are obtained from the modeling process where the governingequations are derived to represent the plant dynamics includingthe errors induced by linearization and model reduction.Unstructured uncertainties, represented by the difference ofthe magnitude of input-output frequency response at each fre-Fig. 1. Closed-loop system interconnection diagram.quency point, are used in this paper. This bound is defined by thenorm of the uncertainty matrix. The overall uncertainty boundis chosen to cover the summation of all uncertainties consideredabove, i.e.,(1)where is bound of the overall uncertainty, is theoverall uncertainty which is the summation of the magnitudesof linearization uncertainty , model reduction uncertainty, and unmodeled uncertainty .The uncertainty bounds are defined as in (2) in a diagonalmatrix form to cover the overall uncertainty. For the purpose ofrobust controller design, it is advantageous to normalize uncer-tainty with a frequency dependent weighting function :(2)where , and.2) Performance Specifications: Selection of weighting func-tions and interconnection of the closed-loop system is an essen-tial step in the synthesis of a robust controller. The weightingfunctions that specify the plant uncertainty and performancespecifications invariably need to be adjusted iteratively. Fig. 1shows the block diagram of closed-loop control system with per-formance specifications, as explained below. The uncertainty weighting function represents the in-tegrated uncertainty bound. The tracking error weighting function specifies the per-formance requirements, which is chosen in such a way thatthe steady-state tracking errors in both channels should besmall (e.g., on the order of 0.01 or less). The control action weighting function is used to atten-uate the control action efforts. That is, if the control actionis excessive, is tuned offline to penalize the control ef-forts, and conversely if the plant response is sluggish. The sensor noise weighting function represents the fre-quency-dependent effects to filter the measurement noise.Upon selection of the above weighting functions, the nom-inal plant model can be interconnected with the weightingfunctions (i.e., , , , and ) to generate the aug-mented plant , as shown in Fig. 1, where the input andoutput of are:810 IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 57, NO. 2, APRIL 2010Fig. 2. LFT description of the robust control problem.is the control input from controller, is the measurement ofplant output, and is the tracking error.The robust control problem is formulated in a two-port frame-work using linear fractional transform (LFT) , as shown inFig. 2. Let and denote the numbers of outputs and inputsof the uncertainty weighting function, respectively, and anddenote the numbers of output and input of the performancespecification weighting functions, respectively. Then, the blockstructure of the uncertainty model is defined as:(3)The first block of the uncertainty matrix corresponds to theuncertainty block , used in modeling the plant uncertainty.The second block includes the performance objectives inthe framework of the -synthesis . The inputs to the secondblock are the weighted control action and weighted trackingerror , and the outputs are the reference , disturbanceand sensor noise signals .3) Controller Design Using -Synthesis: To meet the controlobjectives, a stabilizing controller is synthesized such that, ateach frequency , the structured singular value satis-fies the following condition:(4)where represents the linear fractional transformation(LFT)  of and . The fulfillment of this condition guaran-tees robust performance of the closed-loop control system. The-synthesis can be accomplished by using the D-K iteration toolin MATLAB . For faster computation, the controller ordercould be reduced by eliminating the insignificant states by bal-anced realization and Hankel norm approximation , .B. Resilient ControlThe role of resilient control is to enhance recovery of thecontrol system from unanticipated adverse conditions andfaults as well as from emergency situations. Therefore, theresilient controller should be re-configurable to accommodatewide-range operations and faulty conditions, and this requiresthe resilient controller to be both robust and adaptive. Recently,Hovakimyan et al.  have developed a novel control al-gorithm, called -adaptive control, to address the issues inIntegrated Resilient Aircraft Control (IRAC). The advantagesFig. 3. Closed-loop system with the -adaptive controller.of -adaptive controller are: (i) guaranteed fast adaptation,and (ii) simultaneous tracking of the input and output signalsthat are uniformly bounded. Since the -controller is bothrobust and adaptive, it is suitable for resilient control.1) -Adaptive Output Feedback Control: The -adaptivecontrol architecture was first presented by Cao and Hovakimyan using a state feedback approach in systems with constantunknown parameters. Later, the -adaptive control has beenextended to nonlinear time-varying systems in the presence ofmultiplicative and additive unmodeled dynamics . Its exten-sion to output feedback control has been presented for a class ofuncertain systems  that allows for tracking arbitrary refer-ence with guaranteed time-delay margin. In the work reported inthis paper, the -output feedback adaptive control architecturehas been employed to address the challenge of resilient controlof nuclear power plant, as shown in Fig. 3. While details onthe -adaptive controller are reported in recent literature ,, , the salient features are explained in terms of the fol-lowing single-input single-output (SISO) system model.(5)where is the control input, is the systemoutput, is a strictly proper unknown transfer function ofunknown relative degree for which only a known lowerbound is available, and is the time-de-pendent disturbances and uncertainties. With a slight abuse ofnotation, Laplace transforms of , , and are respec-tively denoted as , , and .Let be a given bounded continuous reference inputsignal. The control objective is to design an adaptive outputfeedback controller giving such that the system outputtracks the reference input following a desired model:(6)where is a minimum-phase stable transfer function of rel-ative degree . The system equations in terms of the de-sired model are rewritten as:(7)whereis the disturbance.JIN et al.: INTEGRATED ROBUST AND RESILIENT CONTROL OF NUCLEAR POWER PLANTS 811Closed-Loop Reference System: Let andbe the reference control input and reference disturbance, re-spectively, of the closed-loop reference system that defines anachievable control objective for the -adaptive controller suchthat:(8)where is a low pass filter with DC gain andis a (possibly) time-varying function ofthe reference output .The transfer matrices and are selected such that(9)is stable and that the -gain of the cascaded system is upperbounded as:(10)Then, the reference system in (8) is stable.Referring to Fig. 3, individual components of the -adaptivecontroller are described next.State Predictor (Passive Identifier): Let ( ,, ) be the minimal realization of the stable transfermatrix . Hence, is controllable and observ-able with being Hurwitz. Then, the system in (5) is rewrittenin the state-space setting as:(11)and the associated state predictor is given by:(12)where is the vector of adaptive parameters. No-tice that, in the state predictor equation may not belongto the space spanned by , while does belong to the spacespanned by in (11).Adaptation Law: Let be the solution of the following al-gebraic Lyapunov equation:(13)where . From the properties of it follows that therealways exists a nonsingular such that(14)Given the vector , let be the -dimen-sional nullspace of , i.e.,(15)and let(16)The update law for at sampling instant (with beingthe sampling time) as:(17)where the state transition matrix(18)and(19)where denotes a Cartesian basis vector in with its firstelement equal to 1 and other elements being 0.Control Law: The control law is defined via the output of thelow-pass filter:(20)The complete -adaptive controller consists of the state pre-dictor in (12), the adaptation law in (17), and the control lawin (20), subject to the -gain upper bound in (10). The perfor-mance bounds of the -adaptive output feedback controller aregiven by the following theorem:Theorem 2.1 (Theorem 1 and Lemma 3 in ):(21)where denotes the norm (i.e., essential supremumof the absolute value) of the function .III. INTEGRATION OF ROBUST AND RESILIENT CONTROLThis section describes how the robust controllers and resilientcontrollers are integrated. A finite-state machine (FSM)  hasbeen adopted for discrete decision-making, which monitors theconditions of the nuclear power plant via a fault detector andcontrols the bumpless transfer between the robust and resilient812 IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 57, NO. 2, APRIL 2010Fig. 4. Layout of the integrated robust and resilient control system.Fig. 5. Example of a finite state machine (FSM).controllers. Fig. 4 depicts the layout of an integrated robust andresilient control system that includes the robust controller, re-silient controller, fault detector, finite state machine (FSM), aset of reference points, and a filter bank for bumpless transfer.The fault detector is extrinsic to both robust and resilient con-trollers.A. Finite State MachineA finite state machine (FSM)  is a 5-tuple, where is a finite set called the states, isa finite set called the alphabet, is the statetransition function, is the start state, and is theset of accepted states or final states.Fig. 5 shows an example of a finite state machine (FSM)that consists of two states ,representing the normal operating condition and an anomalouscondition, respectively. The alphabet in this FSM is selected as, where the symbol 0 represents the normal condi-tion, and the symbol 1 represents the detection of anomaly (e.g.,abrupt temperature change in primary coolant flow). A brief ex-planation of the FSMs operation is presented below.When the nuclear power plant is in the normal operating con-dition (i.e., the symbol of the alphabet is 0), the FSM is in state. When an anomaly occurs in the plant, the symbol changesfrom 0 to 1, and the FSM makes a transition to the anomalousstate . Accordingly, the control system takes necessary actions(e.g., change of the feed water flow set point) are taken to keepthe nuclear power plant safe according to the above procedure.The FSM stays in state until the plant is restored to the normalcondition when the symbol changes back to 0, and thus the FSMmakes the transition back to the normal state . The role of re-silient control is to make the transition from state to stateas quickly and safely as possible to recover from unanticipatedadverse conditions/faults and emergency situations by alteringits operational envelope in real time. The control configurationis restored upon returning to normalcy if the plant is not dam-aged, or to a gracefully degraded condition if the plant is par-tially damaged but still operable within specified safety and per-formance criteria; this is achieved by returning to state fromstate .B. Bumpless TransferSwitching from one controller to another should entail as littleagitation (i.e., occurrence of undesirable transients) as possible.By parallel operation of a non-active controller one could try todrive its output signals towards the correct amplitude, so thatthe resulting transients of the closed loop due to a transfer ofauthority are as small as possible. This is known as the bump-less transfer problem for transition from one operating mode toanother.Many bumpless transfer techniques have been developedfor application in scenarios with different constraints, such asswitching between manual and automatic control, filter andcontroller tuning, scheduled and adaptive controller . Inthis paper, a piecewise linear filter is constructed to mitigate thetransients during switching between the robust controller andthe resilient controller. The key idea is explained below.While the robust controller is active to perform in both steadystate and transient conditions under normal operation, the re-silient controller becomes active to provide safety and recoveryto the plant during adverse conditions and emergency situations.In order to achieve bumpless transfer from robust control toresilient control and avoid abrupt changes in control actions,piecewise linear filters are added to the controller outputs. Theresulting control action is formulated as:(22)where and are outputs of the resilient controller androbust controller, respectively, and and are the filtertransfer functions for the respective control actions.C. Integrated Robust and Resilient Control StrategyFollowing Fig. 4, the next task is to incorporate the finitestate machine (FSM) and filter banks within the plant controlsystem for bumpless transfer between the robust and resilientcontrollers. For example, upon detection of a loss-of-flow acci-dent (LOFA), as the primary coolant temperature crosses a spec-ified threshold, this information activates the transition betweenstates of the FSM. The state transition function provides in-puts to both blocks, namely, set points and filter banks, as seenin Fig. 4. If a significant fault or emergency situation is detectedin the plant, the state transition in the FSM initiates a bump-JIN et al.: INTEGRATED ROBUST AND RESILIENT CONTROL OF NUCLEAR POWER PLANTS 813less transfer between the two controllers and may also alter thethresholds in the set points module. The altered thresholds arethen used as the new set points for the controller under the ab-normal condition. Upon return to normalcy, robust control withthe original set points are resumed. In practice, a fault detectionsystem is needed to identify an abnormal incident as early aspossible, which necessitates incorporation of a fault detectionscheme within the integrated robust and resilient system. Thisis a topic of future research.IV. TESTING AND VALIDATION ON THE IRIS SIMULATORThe proposed integrated robust and resilient control strategyhas been tested and validated on a simulator of nuclear plants.The results of simulation are presented and further experimentalvalidation is planned on research reactors (e.g., Penn StatesTRIGA research reactor ) in the future.A. The International Reactor Innovative & Secure SimulatorThe International Reactor Innovative & Secure (IRIS) sim-ulator of nuclear power plants is based on the design of anext-generation nuclear reactor. It is a modular pressurizedwater reactor (PWR) with an integral configuration of allprimary system components. Fig. 6 shows the layout of theprimary side of the IRIS system that is offered in configurationsof single or multiple modules, each having a power rating of1000 MWt (about 335 MWe) . The nominal reactor coreinlet and outlet temperatures are 557.6 (292 ) and 626(330 ), respectively. The pressurizer, eight steam generators,and the control rod mechanism are integrated into the pressurevessel with the reactor core. There is no huge pipe used toconnect these components. This design avoids the large loss ofcoolant accident (LOCA). The whole control rod mechanismis mounted inside the pressure vessel to avoid failures of thecontrol rod head penetration.As shown in Fig. 6, the integral design of the primary sidealso makes the containment vessel much smaller than the tra-ditional pressure vessel of a PWR. The IRIS reactor coolantpumps are of the spool type, and are located entirely within thereactor vessel; only small penetrations for the electrical powercables are required. The spool pump geometric configurationprovides high inertia/coastdown and high run-out flow capa-bility that contributes to minimize or even mitigate the negativeconsequences of loss-of-flow accidents (LOFAs).A simulation testbed for testing and validation of controlalgorithms has been developed under the project of NuclearEnergy Research Initiative (NERI). The testbed is built usingMATLAB/SIMULINK. This SIMULINK model includes a re-actor core model, a helical coil steam generator (HCSG) model.The turbine is not explicitly modeled in the testbed since thereactors operational safety is the major focus of this paper. Thesimulation testbed is implemented on a Quad Core 2.83 GHzCPU 8 GB RAM Workstation in the laboratory of Penn State.Another two workstations with the same configuration are alsoavailable for simulating a twin-unit plant operation scenario,in which two workstations host an individual IRIS module,and the third one hosts the controller to coordinate the plantmodules through a local network. This testbed is capable ofFig. 6. Layout of the primary side of the IRIS system .simulating normal operation conditions at different operationalmodes as well as various faulty scenarios including: Actuator failures: Feedwater pump trip, malfunctions ofreactor coolant pump and control rod mechanism; Sensor failures: Malfunctions of temperature, pressure,and flow-rate sensors; Internal faults: Uncertainties in fuel temperature coeffi-cient of reactivity, coolant heat capacity, and feedwaterheat capacity.B. Design of the Integrated Robust and Resilient ControlSystemThe details of the integrated control system design are intro-duced in this section.1) Robust Controller Design: The design of a robust con-troller requires judicious selection of the uncertainty weightingfunction and the performance weighting functions, ,and . These weighting functions essentially serve asuser-selected performance specifications.The uncertainty weighting function is selected to cap-ture the frequency-dependent uncertainties as:(23)In order that the plant output tracks the given referencesignal , the tracking error weighting functions with respect toeach plant output are chosen and put into a form of the diagonalmatrix as:(24)This weighting function indicates that the steady-state (i.e., low-frequency) tracking errors due to reference step-inputs in either814 IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 57, NO. 2, APRIL 2010TABLE IRESULTS OF THE -SYNTHESISchannel should be on the order of 0.01 or smaller. This perfor-mance requirement becomes less stringent at high frequencies.The control action weighting function is included totune the control action effort. For example, if the time responseof the controlled system displays excessive control action,the weighting function serves to penalize the energy of theoscillation; on the other hand, if the response is sluggish, theweighting function has a more benign role with a significantlyreduced penalty. The control action weighting function isselected as:(25)In the closed-loop interconnection, effects of frequency-de-pendent sensor noise are represented by the weighting functiondiagonal matrix as:(26)An increased sensor noise weight makes the controller morerobust to sensor noise possibly at the expense of relatively slowresponse.The -synthesis design of the robust controller is accom-plished by using the D-K iteration tool in MATLAB . Forthe IRIS plant, the results of D-K iterations are shown in Table I,where is the -norm of the transfer matrix operator that isa measure of the controllers robust stability, and is the struc-tured singular value that is a measure of the controllers robustperformance. For controllers with and , the robuststability and robust performance are guaranteed. The controlleris selected after five D-K iterations as it yields the smallestand values. The resulting robust controller of order with(sensor) inputs and (control action) outputs in the state spaceform, where , is shown below:(27)As seen in Table I, the synthesized controller yields theclosed-loop after five iterations and the robust per-formance of the control system is guaranteed for the prescribeduncertainty and performance. For faster computation, the con-troller order is reduced by eliminating the insignificant states,where balanced realization and Hankel norm approximationTABLE IIORDER REDUCTION OF CONTROLLER #5,  has been used; the results of controller order reductionare listed in Table II. It is seen in Table II that the -valuesdo not increase significantly after order-reduction until thecontroller order is reduced below 10. There is a clear increasein as the controller order is reduced from 7 to 6, and the valueof for controller order less than 6 is greater than 1. Therefore,the controller of reduced order 7 is selected.2) Resilient Controller Design: The resilient control algo-rithm in Section II-B is derived for a single-input single-output(SISO) system. However, the plant dynamics of the IRIS modelare coupled in two sensor channels, steam pressure and reactorpower, which require an extension of the SISO -adaptiveoutput feedback controller to multi-input multi-output (MIMO)systems. The components of the MIMO control system are de-scribed below.Desired System: The matrix transfer function of the desiredcontrol system is selected in the following form:(28)where and are the scalar transfer functions forcorresponding channels of secondary steam pressure and reactorpower, where zero non-diagonal elements of the matrix transferfunction imply decoupling of the channels in the desired transferfunction. In the current design, minimum-phase stable transferfunctions of relative degree are selected as:(29)where , . The parameter and arechosen in such a way that the frequency response of the desiredsystem is close to that of the linearized IRIS model.Low-Pass Filter: Following Fig. 4, the low-pass filter bank isselected as:(30)(31)where , . The structure and parameters of thelow-pass filter are chosen in such a way that system in(9) is stable and the -stability condition in (10) is satisfied.Consequently, the closed-loop reference system in (8) is stable.3) Smoothing Filter Design for Bumpless Transfer: Bump-less transfer between the robust and resilient controllers isJIN et al.: INTEGRATED ROBUST AND RESILIENT CONTROL OF NUCLEAR POWER PLANTS 815Fig. 7. Normalized plant outputs.achieved by incorporating smoothing filters. During the tran-sition from a normal condition (i.e., state ) to an abnormalcondition (i.e., state ), the filter transfer functions are chosenas:(32)and, during the transition from an abnormal condition (i.e., state) to an normal condition (i.e., state ), the filter transfer func-tions are chosen as:(33)such that the control actions during transition stages are actuallya linear combination of time-dependent weighted outputs fromthe controller actions, and of the resilient and robustcontrollers and vice versa. The respective time constants andare chosen according to performance specifications based onthe following rationale: While a large time constant yields slowand smooth transition, a small time constant ensures fast re-sponse. In this study, the time constants and are chosen as 2sec and 200 sec, respectively. The smaller value of representsrelatively rapid response that is needed for transition to resilientcontrol to deal with emergency situations. The larger value ofrepresents relatively slow response for transition back to robustcontrol as normalcy is restored.A -analysis test confirms that bounded-input bounded-output (BIBO) stability of the augmented closed-loop systemis retained after the addition of the filters, and .Synthesis of a more advanced filter is a topic of future research.C. Results of Testing and Validation on the IRIS SimulatorA scenario of loss-of-flow accident (LOFA) for the (closedloop) system has been simulated on the IRIS to evaluate theperformance of integrated robust and resilient control algo-rithm. Example of causes for LOFA are loss of off-site power,pump failure, heat exchanger blockage, pipe blockage, andfaulty valve closure. In this simulation, the LOFA is causedby reactor coolant pump (RCP) failure of 4 out of 8 reactorcoolant pumps, which is a Condition IV accident. Althoughthis accidental event can be mitigated by the Safe-by-Designfeature of IRIS  (i.e., natural circulation of coolant removesdecay heat from the core even when the coolant pump fails),the plant could still be damaged if appropriate control actionsare not taken.In the first 100 sec of the simulation exercise, the nuclearpower plant is in a normal operating condition with full outputpower load and steam pressure load. In the simulated scenario,when four of the eight primary flow pumps fail, the LOFA isdetected at and the resulting anomaly informationis generated as seen in Fig. 4. The fault detection system is ex-trinsic to both the robust and resilient controllers, and its incor-poration within the integrated control system is a topic of futureresearch.Upon detection of the LOFA incident, the control outputis switched quickly from the robust controller to the resilientcontroller that promptly reduces the set points of output powerload and feedwater flow by half. The plant is brought back tonormalcy at when all primary flow pumps becomefunctional. To make the plant return to normalcy, the outputpower load and feedwater flow set points are reset to theirrespective nominal values, and the control output is switchedbumplessly from the resilient controller back to the robustcontroller. Note that this bumpless transfer is ensured by thechoice of a relatively large filter time constantrelative to the tenure of the LOFA.Fig. 7 shows the normalized plant outputs, steam pressureand reactor output power . Under normal conditions, bothplant outputs are at their nominal values, which are normalizedto 1. When four of the eight primary reactor coolant pumpsfail, the plant is not able to generate 100% power, and thus theoutput power load is reduced to 50% by resilient control ac-tions. Upon returning to normalcy, the output power graduallyreturns to 100%. Note that there exist spikes in the reactor powerwhen the anomaly occurs and it is removed. Those spikesare generated due to temperature feedback from the reactor asa consequence of the abrupt change in primary flow when fourof the eight pumps fail. Note that, in practice, a change in thefeedwater flow may not occur as a step, but may have a rela-tively more gradual but sufficiently fast profile. Furthermore, ifthe output power in the control system is taken to be the turbineoutput power instead of the reactor core output power, the spikeswould be averaged out due to the mechanical inertia of the tur-bine. In the future work, a turbine model will be integrated intothe nuclear plant model.Fig. 8 shows the responses of the control actions, feedwaterflow and rod reactivity , after occurrence of the LOFA.When the plant is operated under normal conditions, the controlefforts are indicated as 0% implying zero deviations from thenominal values of and . Upon detection of the LOFA at, in order to save the plant, the feedwater flow isgradually increased by 10% of the nominal value and the rodreactivity is reduced by more than 80% of the nominal value.Fig. 9 shows the temperatures (i.e., primary coolant inlet, pri-mary coolant outlet, and secondary coolant steam temperatures)monitored by sensors inside the plant, and is the tempera-ture difference between the primary coolant inlet flow and outletflow, which is used to activate the FSM (see Figs. 5 and 4). Upondetection of the LOFA, the control system reduces the rod re-activity by more than 80%. Consequently, the primary coolantinlet temperature drops.As seen in Fig. 10, switching of the symbols (i.e., 0 and 1)of the alphabet are controlled by the FSM. When the plant isworking in normal conditions, the symbol is 0, the FSM staysin normal state , and the output power reference point is 1.816 IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 57, NO. 2, APRIL 2010Fig. 8. Control actions of integrated robust and resilient controller.Fig. 9. Temperatures of steam, primary inlet and outlet flows.Fig. 10. Finite state machine alphabets and output power set point changes.When the anomaly occurs, the symbol changes to 1 immedi-ately, which causes the FSM to transit from state to andthe output power reference point to change from 1 to 0.5.V. SUMMARY, CONCLUSIONS & FUTURE WORKThis paper investigates a new concept of integrated robust andresilient control for enhancement of the operational safety andperformance of nuclear power plants. The robust controller isdesigned by using the -synthesis tools with D-K iteration .The optimal Hankel approximation is used to reduce the orderof the robust controller . The concept of resilient controlin nuclear plants is derived from Integrated Resilient AircraftControl (IRAC)  that is an active area of research in Na-tional Aeronautics and Space Administration (NASA); this isa direct transfer of technology from aeronautics to nuclear en-gineering. The resilient control is implemented by -adaptiveoutput feedback controller, which guarantees fast adaptation,uniformly bounded transient and asymptotic tracking for bothinput and output signals of the control system simultaneously.The robust and resilient controllers are integrated by utilizinga finite state machine and smoothing filters to ensure bumplesstransfer from robust to resilient control modes and vice versa.In the simulation example of nuclear power plant control, thefinite state machine is designed to switch between normal andabnormal states and activate related control actions by moni-toring plant coolant temperatures. The filters ensure bumplesstransfer between robust controller and resilient controller andavoid impact to the nuclear reactor during transition stage.A simple scenario of loss-of-flow accident (LOFA) has beeninvestigated on the International Reactor Innovative & Secure(IRIS) simulator . Simulation results based on this LOFAscenario show that the proposed controller recovers from theemergency situation with a fast response, while the characteris-tics of the standard robust controller are retained during normaloperating conditions. Further analytical, simulation and exper-imental research is necessary before this novel concept of inte-grated robust and resilient control can be considered for appli-cation to commercial nuclear power plants. The following area few examples of future research in integrated robust and re-silient control of nuclear power plants: Incorporation of fault detection schemes within the inte-grated control system; Filter design for fast and stable switching between robustand resilient control modes; Construction of finite-state machines for incorporationwithin the integrated control system to represent variousfailure states and interstate switching; Demonstration of safe recovery from various real-lifeemergency situations by including balance of plant com-ponents in the IRIS simulator; Experimental validation of safe recovery under selectedaccident scenarios on research reactors (e.g., Penn StatesTRIGA research reactor ).ACKNOWLEDGMENTThe authors would like to thank Prof. Naira Hovakimyan ofthe University of Illinois at Urbana Champaign for her technicalsupport in developing the resilient controller. The authors arealso grateful to the thoughtful and meticulous comments of theanonymous reviewers.REFERENCES H. G. Kwatny, B.-C. Chang, and S.-P. Wang, Static bifurcation in me-chanical control systems, in Bifurcation Control. Berlin, Germany:Springer, 2004, vol. 293, pp. 6781. R. M. Edwards, K. Y. Lee, and A. Ray, Robust optimal control of nu-clear reactors and power plants, Nucl. Technol., vol. 98, pp. 137148,May 1992. Z.-Y. Huang, R. M. Edwards, and K. Y. Lee, Fuzzy-adapted recursive,sliding-mode controller design for a nuclear power plant control, IEEETrans. Nucl. Sci., vol. 51, no. 1, pp. 256266, Feb. 2004. K. Hadad, M. Mortazavi, M. Mastali, and A. A. Safavi, Enhancedneural network based fault detection of a VVER nuclear power plantwith the aid of principal component analysis, IEEE Trans. Nucl. Sci.,vol. 55, no. 6, pp. 36113619, Dec. 2008. M. G. Na and B. R. Upadhyaya, Application of model predictivecontrol strategy based on fuzzy identification to an SP-100 spacereactor, Ann. Nucl. Energy, vol. 33, no. 1718, pp. 14671478,Nov.Dec. 2006. R. N. Banavar and U. V. Deshpande, Robust controller design for anuclear power plant using optimization, IEEE Trans. Nucl. Sci.,vol. 45, no. 2, pp. 129140, Apr. 1998. S. S. Shyu and R. M. Edwards, A robust multivariable feedforward/feedback controller design for integrated power control of boiling waterreactor power plants, Nucl. Technol., vol. 140, no. 2, pp. 129140,Nov. 2002. K. Zhou, J. Doyle, and K. Glover, Robust and Optimal Control.Upper Saddle River, NJ: Prentice-Hall, 1996. T. H. Guo and J. S. Litt, Resilient propulsion control research for theNASA Integrated Resilient Aircraft Control (IRAC) project, in Proc.AIAA Conf. and Exhibit, Rohnert Park, CA, May 710, 2007, AIAA-20072802.JIN et al.: INTEGRATED ROBUST AND RESILIENT CONTROL OF NUCLEAR POWER PLANTS 817 C. Cao and N. Hovakimyan, Design and analysis of a novel adap-tive control architecture with guaranteed transient performance, IEEETrans. Autom. Control, vol. 53, no. 3, pp. 586591, 2008. M. S. Mahmoud, Resilient Control of Uncertain Dynamical Systems.Berlin, Germany: Springer, 2004. X. Jin, R. Edwards, and A. Ray, Integrated robust and resilient controlfor nuclear power plants, in Proc. 6th ANS Int. Topical Meeting NPIC& HMIT, Paper 3 Session I & C Grid Appropriate Reactor, Knoxville,TN, Apr. 59, 2009. M. Carelli, IRIS: A global approach to nuclear power renaissance,Nucl. News, vol. 46, no. 10, pp. 3242, Sep. 2003. M. D. Carelli, The design and safety features of the IRIS reactor,Nucl. Eng. Des., vol. 230, pp. 151167, 2004. G. J. Balas, J. C. Doyle, K. Glover, A. Packard, and R. Smith, -Anal-ysis and Synthesis Toolbox Users Guide. Natick, MA: MathWorks,2001. C. Cao and N. Hovakimyan, adaptive controller for systems in thepresence of unmodeled actuator dynamics, in Proc. 46th IEEE Conf.Decision and Control, New Orleans, LA, Dec. 2007, pp. 891896. C. Cao and N. Hovakimyan, adaptive output feedback controllerfor non strictly postive real multi-input multi-output systems in thepresence of unknown nonlinearities, in Proc. Amer. Control Conf., St.Louis, MO, Jun. 1012, 2009, pp. 51385143. M. Sipser, Introduction to the Theory of Computation. Boston, MA:PWS Publishing, 1997. W. S. Levine, Control Systems Fundamentals. Boca Raton, FL: CRC,2000. D. M. Fouquet, J. Razvi, and W. L. Whittemore, TRIGA researchreactors: A pathway to the peaceful applications of nuclear energy,Nucl. News, vol. 46, no. 12, pp. 4656, Nov. 2003.