Operational Resilience of Nuclear-Renewable Integrated-Energy Microgrids

The increasing prevalence and severity of wildfires, severe storms, and cyberattacks is driving the introduction of numerous microgrids to improve resilience locally. While distributed energy resources (DERs), such as small-scale wind and solar photovoltaics with storage, will be major components in future microgrids, today, the majority of microgrids are backed up with fossil-fuel-based generators. Small modular reactors (SMRs) can form synergistic mix with DERs due to their ability to provide baseload and flexible power. The heat produced by SMRs can also fulfill the heating needs of microgrid consumers. This paper discusses an operational scheme based on distributed control of flexible power assets to strengthen the operational resilience of SMR-DER integrated-energy microgrids. A framework is developed to assess the operational resilience of SMR-DER microgrids in terms of system adaptive real-power capacity quantified as a response area metric (RAM). Month-long simulation results are shown with a microgrid developed in a modified Institute of Electrical and Electronics Engineers (IEEE)-30 bus system. The RAM values calculated along the operational simulation reflect the system resilience in real time and can be used to supervise the microgrid operation and reactor’s autonomous control.


Introduction
Scientific evidence shows that the frequency and severity of extreme weather events are worsening day by day with global warming, which is harshly impacting the resiliency and integrity of the electrical grid structure [1]. Weather-related power outages account for more than 70% of the total power outages in the U.S. [2]. Between 2003-2012, there had been more than 679 widespread blackouts in the U.S. affecting at least 50,000 customers [3]. More recently, the winter storm Uri in February 2021 resulted in power outages for roughly 10 million people in Texas, causing cascading failure to other facilities relying upon electricity, including drinking water treatment and medical services [4]. It resulted in an estimated $130 billion in damage in Texas alone and reportedly the death of 210 people. Elsewhere outside the U.S., widespread flooding in Queensland, Australia in 2011-2012 resulted in severe damage to transmission facilities causing power disruptions to 150,000 customers [5]. In 2008, a severe ice storm in China resulted in the isolation of 2000 substations and the collapse of 8500 towers leading to power outages in 13 provinces and 170 cities [5].
With centralized electrical grids proving more vulnerable to extreme events, energy consumers are looking to develop microgrids with generation and control capacity to fulfill their demand. Furthermore, microgrids can add value and be compensated for providing power and services to the grid. By aggregating local resources, including energy storage, microgrids can decentralize larger grids into smaller entities, each capable of operating independently when the grid fails. While distributed renewable energy resources (DERs), such as wind and solar photovoltaics (PVs), are central to such microgrid architecture due frequent and amplified, it is apparent that this classical approach of system planning needs significant upgrades. Furthermore, the aggressive shift of the power systems towards renewable energy, which themselves are affected by weather conditions, necessitates a new approach capable of reflecting impacts of common-cause failures. Considering its emphasis on extreme disruptive events, the concept of resilience has received significant research attention in the power system domain in recent years.
The concept of resilience was first introduced by C.S. Holling in 1973 as a measure of a system's ability to absorb disturbances while maintaining unaltered relationships between its components and parameters [23]. Over the years, resilience has been adopted by different domains where it had numerous interpretations and evolutions. Conceptualizing the quantification of disaster resistance against extreme events, multidisciplinary research sponsored by the National Science Foundation utilized the resilience triangle to reflect the loss of functionality after an event and subsequent recovery and restoration [24]. The R4 framework describing robustness, redundancy, resourcefulness, and rapidity of the system against disturbances was introduced to describe system resilience. Cimellaro et al. adopted the resilience triangle to quantify the disaster resilience of hospital buildings [25]. Ouyang et al. extended the resilience triangle into a multi-stage framework by including one additional damage propagation stage between the inception of the event and the start of the recovery [26,27]. The framework was implemented to analyze the extreme event resilience of urban infrastructure systems in [26] and power systems in [27].
In the last two decades, several research works have attempted to define, qualify, quantify, and analyze power system resilience and develop metrics for comparison. References [22,28] define power system resilience as a system's capability to withstand abnormal high-impact disruptive events including extreme weather events and cyberthreats. References [29,30] further include the operational adaptation as a part of the theoretical definition of the power system resilience. Reference [31] distinguishes the resilience by design and resilience by real-time reaction focusing on disturbances related to cyberthreats. Panteli and Mancarella discussed the influence of extreme events on power system components to provide guidance to develop a modeling framework for power system resilience [32]. They also used the concept of a resilience trapezoid in the context of power systems, breaking the event into preventive, corrective, emergency, restorative, and adaptive stages [5,22]. To quantify the power system resilience with trapezoidal illustration, the "ΦΛEΠ" resilience metric, pronounced as "FLEP", was introduced in [22]. Panteli et al. discussed a unified resilience evaluation approach and utilized risk-based defensive islanding to mitigate cascading failures [28].
Chanda et al. discussed a generic code-based metric to quantify the temporal characteristics of power system resilience for an ongoing contingency in a power distribution system [33]. Rieger discussed the quantification of resilience in the context of control systems integrating its physical, cyber, and cognitive aspects [34]. The disturbance and impact resilience evaluation (DIRE) curve was introduced that helped classify a system as resilient or non-resilient based on the attributes reflected in the graph which include robustness, agility, adaptive capacity, adaptive insufficiency, resiliency, and brittleness. Eshghi et al. discussed an agent-based resilience model to develop metrics based on four power system attributes, namely small signal stability, transient stability, communication latency, and sensor physical degradation, recognized in terms of time and data of the measurements [35]. The metrics proposed were intended for the supervisory control and data acquisition system to improve the observability and controllability of the system against for three different classes of disturbances: physical, cyber, and cognitive.
In [29], McJunkin and Rieger used the concept of adaptive capacity to assess and quantify the operational resilience of a modern distribution system (MDS). Instead of focusing on system performance following a disruptive event, the adaptive capacity framework evaluates the system adaptive capacity to respond to such disruptive events. The adaptive capacity refers to the dynamic flexibility of the system expressed in terms of real and/or reactive power to respond to the anticipated disturbances. The adaptive capacity curve has both time and magnitude attributes influenced by the latency, ramp rates, and magnitude of power variation of system assets. Looking forward in time, the magnitude and temporal characteristics of energy-limited and energy-independent assets in a power system are modeled and aggregated to produce system-level adaptive power characteristics that reflect system resilience against anticipated disturbances. In addition to the high-impact disruptive events, the uncertainty in renewable generation (e.g., wind and PVs) and load demands are important issues to consider while analyzing the operational resilience of microgrids hosting large amount of DERs. The adaptive capacity framework analyzes the system reserve to respond to any form of generation-demand imbalance, whether it is from disruptive events or from DER intermittency and uncertainty. This feature of adaptive capacity framework that focuses on system capacity to respond, rather than the consequences of an event, makes it useful to analyze system resilience against uncertainties. In [29], the adaptive capacity curve generated for an MDS is mapped into the DIRE curve to understand the system capability at different stages following disturbances and the roles of different assets in those stages. The operational resilience framework in [29] revolves around the concept of adaptive capacity, which has been implemented for different generation assets. Phillips et al. [36] extended the work in [29] to develop a real-time asymmetric operational metrics by updating the adaptive capacity with control thresholds, contingencies, and environmental conditions and demonstrated the proposed model on a case system for both normal and abnormal scenarios. In [37], Phillips et al. analyzed assetlevel flexibility for three classes of hydropower plants, namely, run-off river hydropower, hydropower with reservoir, and pumped-storage hydropower, and developed a framework to integrate their resilience contribution to overall system adaptive capacity characteristics. In [38], the resilience contribution of a hybrid solar PV-battery system was analyzed in terms of adaptive capacity considering the generation uncertainty.
This paper aims to leverage the concept of adaptive capacity in the context of SMRs to develop operational resilience framework for SMR-DER integrated-energy microgrids. The balance between real-power flexibility needs and requirements of SMR-DER microgrids is one of the decisive factors for its operational feasibility. The reactive power, on the other hand, is a generator-centric capability and can also be compensated with localized compensation schemes, such as shunt capacitors and flexible alternating current transmission system (FACTS) devices [39]. Similarly, the power system inertia is based on generator mass and is independent of the energy source driving the prime mover. In other words, the reactive power and inertial model of the SMRs will be similar to other thermal power plants. For this paper, we neglect reactive power and inertia, making the assumption that sufficient inertia or synthetic inertia is present to maintain frequency stability, and sufficient voltage regulation capability is also available. The proposed framework is constructed to assess the operational resilience of an SMR-DER integrated-energy microgrid case system developed for electricity and district heating in a modified IEEE-30 bus network. The major contributions of this paper include:

1.
A resilience framework based on system adaptive real power of SMR-DER integratedenergy microgrids, taking into consideration the salient operational features of SMRs including flexibility added with cogeneration and the interaction between electrical and the heating system.

2.
A resilience metric termed as RAM is proposed for anticipated normal and compromised operational scenarios to provide resilience feedback to the microgrid's asset management system and the reactor's autonomous control.
The novelty of this research work lies in extending the adaptive capacity framework to model the characteristics of a nuclear system and the complex interaction between heat and power sides for a microgrid's operational resilience. Compared to earlier resilience works in power systems that mostly focused on modeling power system components, this research explores all three domains of SMR-DER integrated-energy microgrids (i.e., nuclear, heating, and electrical). Furthermore, this research also attempts to characterize the impact of interaction between the heating and electrical sub-systems on a system's operational resilience. The safety aspect of the nuclear system is added to the resilience model which includes fixed operational pattern and power maneuvering limits of SMR, and real-time decision-making to keep the reactor operating or shut it down based on the disturbance characteristics. The proposed resilience model aims to boost the operational awareness which is critical for SMR-DER microgrids where the flexible assets are limited, and the safety of the nuclear system cannot be compromised. This paper builds on the outcomes of papers that developed a simulation model of an SMR-DER integrated-energy system in [20] and analyzed its optimal design and operation in [21]. The system models and optimization framework from [21] are utilized in this paper to develop and test the operational resilience model of SMR-DER integrated-energy microgrids.

System Description
The proposed microgrid system, as depicted in Figure 1, consists of an SMR as the primary source of heat and power and hosts distributed small-scale wind and PV generation. The load in the microgrid can be divided into critical and non-critical loads. One of the assumptions here is that the critical loads cannot be shed, while non-critical loads can be curtailed as required to meet critical load demand. To have the capability to island, the microgrid's assets are sized to meet, at a minimum, the electricity demand of local critical loads. Additional capabilities can be considered if there is economic value in the potential export of power to the grid. An SMR may be a single-module or multimodule plant based on the availability of sizes to meet the system demand. The sizing of distributed wind and PVs and their pairing with battery energy storage (BES) or other storage mechanisms will be based on an assessment of expected critical load over time profile versus the expected availability and intermittency of the wind and sun. BES can support peak load shaving, support long-term storage, and provide synthetic inertia or primary frequency response. SMR-DER microgrids, as shown in Figure 1, can also meet local heating requirements along with electricity. The steam from an SMR is used to generate electricity and contribute to local heating as needed. The centralized thermal energy storage (TES) will act as a buffer for a heating network to absorb the fluctuations in heating demand and steam supply from a reactor.  Microgrids typically have a limited number of flexible assets to support system operation and respond appropriately to demand fluctuations or local disturbances. Therefore, the sizing of assets and operational strategy chosen will have significant impact on a microgrid's flexibility. Although SMRs are expected to support flexible operation, they still face technical limitations in terms of their capability to provide flexible power. Frequent operation of control rod for large power changes creates issues such as thermal fatigue and aging of the reactor and thermal components, erosion and corrosion of hydraulic components such as valves and pipes, the core power redistribution, and fission product poisoning [40]. Due to these technical challenges, the reactor control is restricted within the design limits typically defined in terms of the rate of change, the total change, and the total number of large power cycles. Therefore, an SMR's flexible operation should follow certain operational protocols based on its design limits. In general, an SMR's salient operational characteristics can be summarized as below [11,40,41]: Reactor power maneuvering restricted in terms of total power change, ramp rates, and the total number of power maneuver cycles.

2.
Non-zero limits on minimum power level in addition to the maximum limits typical to all power plants. 3.
Fixed standard operating schedule and daily power-cycle to reduce operational uncertainty to maximize system economics and reactor safety.

4.
Multiple means to provide electrical power change: steam extraction, steam bypass, and reactor control rod maneuvering.
These characteristics are common to SMRs with water/steam in secondary coolant circuit, which represents the majority of current SMR technologies. Note that the design configuration shown in Figure 1 is an example considered for the purposes of this paper. The resilience model developed in this paper for an SMR-DER microgrid with integrated heat application can be adapted to other system configurations, including power-only SMR-DER microgrids. Considering SMRs with the above characteristic features, an operational scheme was introduced in [20] decoupling SMR's power control into control rod maneuvering, steam extraction, and steam bypass. The simulation results discussed in [21] analyzed the adaptiveness of the operational scheme against uncertainty and intermittency of renewables and showed that an SMR can economically host as much as 50% renewables consisting of wind and PV generation. This paper adopts the operational scheme from [20] to develop the operational framework for microgrid control. The next three subsections discuss the operational framework, modeling of assets, and the resilience framework.

Operational Framework
The operational framework enables control of all the flexible generation assets to meet instantaneous generation-demand balance. Figure 2a illustrates the distributed control in which overall operation is divided into three parts: coarse-load shaping (CLS), load following (LF), and frequency control (FC). The CLS and LF are setpoint-based controls where operational setpoints are identified beforehand, whereas the FC is automatic control based on instantaneous power system frequency feedback. These three controls are decoupled and are executed in different timescales independently.   Figure 2b shows the real-time operational scheme for setpoint-based controls. The setpoint-based controls (i.e., CLS and LF) are repeated at fixed time intervals. During each interval, the setpoint prediction and operation blocks work in parallel. The prediction block identifies the setpoints for the next interval. The beginning of each interval is referred to as the decision point because it is the moment when the setpoint predicted over the earlier interval has to be executed. The operation block executes the setpoint evaluated by the prediction block at the decision point. The execution of CLS and LF setpoints will bring the generation close to the load demand. FC assets are then operated with frequency deviation feedback to meet the instantaneous generation-demand balance.

Coarse-Loading Shaping
The CLS is the outermost control to coarsely follow the shape of the daily electrical load demand. It is a semi-baseload operation with very few variations per day and with slower ramp rates. It is executed with the help of reactor maneuvering which has strict operational limitations in terms of ramp rate, total change, and the number of reactor maneuvers. In each CLS interval, denoted as t CLS , the reactor power level transitions to a new setpoint. The CLS interval or reactor maneuvering cycle consists of ramping and hold periods denoted as t R and t H , respectively. The reactor ramps to a new setpoint over the ramping period and stays at the new reactor power level throughout the hold duration. For the purposes of this paper, t CLS , t R , and t H are taken as 6 h, 2 h, and 4 h, respectively. The reactor power level setpoints are calculated based on the information analyzed by the decision point. The objective of CLS is to coarsely bring power generation close to demand such that LF and FC can cover short-term fluctuations to ensure instantaneous generation-demand balance. During each CLS interval, the prediction block calculates a new setpoint for the next CLS interval. The CLS prediction block is a mathematical function that captures information, such as predicted electricity and heating demand peaks for the next CLS interval and the loss of state of charge (SOC) of BES and TES, to calculate a new setpoint for the reactor power level. Equation (1) shows the calculation of reactor power setpoint for k + 1 th CLS interval as a function of different state variables and demand predictions: where P R,th k+1 is the reactor thermal power level setpoint at k + 1 th interval, P R,th k is the reactor thermal power level at the k th interval, P D k+1 and Q D k+1 are the predicted peak electrical and heating demand at the k + 1 th interval, respectively, and SOC TES k.t CLS and SOC BES k.t CLS are the SOC of the TES and BES at the decision point, respectively. The SOC of TES is represented by its temperature, and that of BES is represented by the stored energy in pu of the BES energy rating.
If the SMR plant has multiple reactor units, only one unit is maneuvered at a time in each CLS interval. This will help reduce the total number of reactor maneuvers per reactor unit [21]. However, if required, multiple units can operate simultaneously to provide larger power variations.

Load Following
The LF is another setpoint-based control with the interval of t LF . For the purposes of this paper, t LF is considered as 15 min. The generation-demand imbalance left after the CLS will be handled with the LF control. In the proposed microgrid configuration shown in Figure 1, the steam produced by the reactor is used for electricity and heat applications. The steam extraction to the heating side can be controlled to achieve power variation on the electrical side. This capability will be used for the LF. In each LF interval, the steam extraction level is changed to a new setpoint value. Meanwhile, the LF prediction block calculates steam extraction level setpoint for the next interval that will meet electrical demand predicted for the next interval. The steam extraction flow setpoint,ṁ EV j , is therefore a function of the electrical demand predicted over the next LF interval: where P D j is the electrical demand predicted for the j th LF interval. As mentioned earlier, the LF and CLS are decoupled from each other and controlled independently. Therefore, the sampling indices for CLS and LF (i.e., k and j, respectively) are independent of each other. Note that, unlike CLS interval, LF interval does not have a hold period.

Frequency Control
The setpoint-based controls do not consider the short-term variation within the sampling interval. Furthermore, there could be mismatches left due to the biases and uncertainties in predicted load demand and renewable generation. The short-term generationdemand imbalance left after setpoint-based controls is responded to with instantaneous FC with the help of BES. The BES control will be based on the power system frequency deviation feedback which represents the dynamic generation-demand imbalance. The instantaneous electrical power output of BES, denoted as P B t , is given as where K PB and K IB are the proportional and integral gain of the proportional-integral (PI) controller, respectively. ∆ f t is instantaneous frequency deviation from the rated frequency. If BES alone is unable to absorb short-term disturbance, the turbine bypass system, equipped in each reactor unit (see Figure 3), can support the BES by bypassing the process steam directly to the condenser prior to sending to turbine stages (i.e., high-pressure (HP) and low-pressure (LP) turbine stages). Note that the turbine bypass system has different setups than the steam extraction system. The steam extraction system extracts the steam from between turbine stages.
The operation of the turbine bypass system will reduce the reactor's electrical power and help meet instantaneous generation-demand balance on the power system side. The size and capability of the turbine bypass system may vary with the reactor designs. The NuScale SMR's turbine bypass system, for example, can bypass up to 100% steam at full reactor power directly to the condenser without shutting down [42]. For the purposes of this paper, the turbine bypass system of the NuScale SMR is assumed. Bypassing a large amount of high-temperature steam for long durations can cause wear and tear of condenser shell and tubes, steam bypass lines, and associated valves. Therefore, the turbine bypass is reserved for certain abnormal scenarios when BES is unable to absorb demand fluctuations. Depending on the disturbance and the anticipated recovery time to return to the normal operation, the reactor unit can remain operating or shutdown on bypass flow at full or reduced reactor power level [42]. The turbine bypass flow is controlled using a simple PI controller, as given by (4).
whereṁ BV t is instantaneous bypass flow; K Pbv and K Ibv are the proportional and integral gain of the bypass valve controller, respectively. ∆ f set is frequency deviation threshold or the dead band for bypass valve.

Small Modular Reactor
The steam produced by the SMR is used for electricity and heating applications. It flows through the bypass valve, HP turbine, extraction valve, and LP turbine before merging at the condenser. The condenser then supplies the feedwater back to the steam generator (SG) through the pump. Figure 3 shows the steam distribution of an SMR with water/steam in secondary coolant cycle. Note that there could be other components in the secondary coolant circuit in the actual SMR plant, including additional intermediate turbines, moisture separator and reheater between turbines, multiple pumps, and feedwater preheaters and deaerator in the feedwater line [42,43]. The resilience model of an SMR will vary based on the configuration of the reactor system. Specifically, the impact of steam extraction and turbine bypass systems on feedwater preheating and system efficiency needs to be analyzed. However, for the purposes of this paper, the reactor secondary coolant circuit is simplified without reheat and preheat stages assuming an ideal, basic Rankine cycle [44].
The SMR model discussed below is based on the dynamic models developed in [20,45] and is applicable to all SMRs with water/steam in secondary coolant cycle. The steam produced by the SG is a function of reactor thermal power level [45]. This statement is valid during reactor maneuverings due to slower ramp rates. The total steam for k th CLS interval, denoted asṁ SG k , is given asṁ and the instantaneous value of steam flow is given aṡ Similarly, based on the LF control, the instantaneous value of extraction steam is expressed asṁ Ifṁ BV t andṁ EV t are the steam flow rates to the bypass and extraction valves, respectively, the steam flow rates through HP turbine (ṁ HP t ) and LP turbine (ṁ LP t ) are given as The total electrical power output of the SMR is given as where P R,e t is the electrical power output of an SMR unit at time t; η g is the electrical generator efficiency; η HP and η LP are the efficiency of HP and LP turbines; and ∆h HP t and ∆h LP t are the enthalpy difference between the steam across the HP and LP turbines. Figure 4 shows the feasible operation region of an SMR unit. Reactor design bases specify the minimum power level that the reactor can operate at during normal operation. At points a and b, the reactor is operating at full power with 0% and 100% steam extraction, respectively. At points c and d, the reactor is at minimum power level with 0% and 100% steam extraction, respectively. The vectors along ac, bd, and corresponding interpolating vectors inside the envelope represent the reactor maneuvering control. The vectors along ab, cd, and the ones parallel to them inside the envelope represent the steam extraction control.
In LF control, a portion of steam generated by the SMR is extracted by using the extraction valve located between the HP and LP turbines to reduce its electrical power output. The turbine stages before the extraction valve receive the full steam while the turbine stages after the extraction valve receive only the portion of steam left after the steam extraction. Therefore, the LF control can only reduce the electrical power generated by LP turbine and not the HP turbine. Reducing the reactor power level reduces the total steam generated, consequently reducing the LF capability. The feasible operation region shown in Figure 4 provides the static limits of the reactor. Additional limits are imposed in terms of the maximum change allowed per power maneuvering cycle and ramp rates. and where ∆P R,th max is the maximum change allowed for reactor power level in a single maneuver under a normal operating scenario. R D th and R U th are the ramp rate limits for the reactor maneuvering. The steam flow rates through bypass and extraction valves are less than or equal to the total steam produced by the SG.
The operation of bypass valve for short-term response should abide by the valve operation ramp rates as given by (15): where R D BV and R U BV are the ramp rate limits for the bypass valve operation.

Heating System
The heating system model discussed here is based on [20]. Figure 4 shows the distribution steam produced by a single SMR unit. In a multi-module plant, steam extracted from multiple reactor units are collected and supplied to the heating side. Note that the heat exchanges shown in Figure 4 are part of the heating system and not from the reactor secondary coolant circuit. The heating system constitutes the components between the steam extraction valve and the condenser: heat exchangers (HX1 and HX2), TES, heating network, and the heat consumers in the microgrid. For simplicity, the heating demand is assumed to be lumped as a single heat load. Because the control of steam extraction is reserved for electrical side LF, the heating side receives the steam excess to the electrical side irrespective of the heating demand. The TES is responsible for meeting generation-demand balance on the heating side by acting as a buffer between the steam supply and heating demand. The change in the SOC level of the TES sustained while providing flexibility to the heating side is recovered by adjusting the reactor power level setpoint in the subsequent reactor maneuvering cycle. The heating demand will typically have slower and smaller short-term fluctuations which are absorbed by the TES. Therefore, the heating system is simulated using quasi-static models with LF interval as the simulation timestep. The total steam supplied to the heating sideṁ EX j in the j th LF interval will be the sum of the steam extracted from each reactor unit in the SMR plant: where N R is the total number of reactor units in the SMR plant. Similarly, the enthalpy of total extracted steam flow (∆h EX j ) can be expressed as the average of individual stream enthalpies.
where ∆h EV j,k is the useful enthalpy of extracted steam in the j th LF interval. It is assumed that the latent heat of the steam is utilized. Now, the total heat extracted to the heating side Q EX j in the j th LF interval is given as The extracted steam is typically either excess or insufficient for the heating system. It is further distributed among two heat exchangers-one supplying heat directly to the heating load and another storing heat in the TES. If the steam extracted is not sufficient, heat stored in the TES is used to fulfill the heating demand. If the extracted steam is excess to the heating demand, the surplus steam is forwarded to the TES. If the TES is at full capacity, it will not be able to absorb the heat from the excess steam and will be bypassed directly to the condenser: whereṁ HX1 j andṁ HX2 j are the steam flow rate to the HXs 1 and 2, respectively;ṁ EC j is the steam bypassed to the condenser from the extraction line; Q D j is the heating demand for the j th LF interval; T TES j is the temperature of TES representing its SOC during the j th LF interval; and T TES max is the maximum limit on TES temperature. These flow rates are further limited by the design flow limits of heat exchangers,ṁ HX1 max andṁ HX2 max . Assuming sensible heat TES with hot water storage tank, the temperature of TES in the j th LF interval is given as where Q TSi j is the heat received by the TES from HX2, and Q TSo j is the heat supplied by TES to the heating load during the j th LF interval. The heat input and output of TES can be expressed as If TES temperature reaches its lower limit, T TES min , it will not be able to supply heat. In such a case, the amount of steam supplied by SMRs to the heating side should at least meet the instantaneous heating demand, and the electrical side will partly lose the flexibility for LF for power-up capability. Similarly, if TES reaches its upper limit, it cannot absorb the excess heat, and, therefore, the steam extraction system can only extract steam equal to or less than what is needed to by the instantaneous heating demand. Mathematically, In the case of a multi-module plant, the coordination of steam extraction is necessary for reactor units to meet the required amount of steam on the heating side. The responsibility could be shared equally by multiple units, or, alternatively, some of the units could meet heating requirements while the rest offer the LF response.

Battery Energy Storage
The electrical power output of the BES is controlled using a PI controller with frequency deviation as feedback, as given by (3). The instantaneous energy state of BES or its SOC, E B t , is given as where E B 0 is the initial SOC, and η B is the battery efficiency, which is different for charging and discharging. BES can provide the FC response only if the SOC level is within minimum and maximum limits. If the BES SOC is below the lower operating limit, B < E B min , it will not be able to supply the power to the microgrid and can only either charge or remain at zero power (recovery state). Similarly, if the BES SOC is above the upper operating limit, B < E B max , it can either discharge or remain at zero power. Mathematically, where E B max and E B min are the upper and lower operational limits for BES SOC levels, respectively. The operation of BES for short-term response should obey the ramp rate limits as given by (28): where R D B and R U B are the ramp rate limits for BES power control.

Wind and PV Power Output Models
The power outputs of renewable assets are captured in full without any generation curtailment. In other words, the wind and PV plants will not add to the power flexibility but rather contribute to the disturbance due to their generation intermittency. The power output model for the wind plant is given as where P w t is the instantaneous electrical power output of the wind power plant; N WT is the number of wind turbines; P w r is the rated electrical power of a single wind turbine; ω t is the instantaneous wind speed; ω ci , ω co , and ω r are the cut-in, cut-out, and the rated speed of wind turbine, respectively; and A, B, and C are the turbine parameters calculated using cut-in, cut-out, and rated speed values [46]. The power output of a PV plant is given as [47] where P PV t is the instantaneous electrical power output of a PV plant; G PV t is instantaneous solar irradiance; K PV is an arbitrary coefficient that relates power output of a PV plant with the solar irradiance; η p is the panel efficiency, A p is the panel area; and T a is the ambient temperature in • C.
In normal operation, the objective would be to capture the maximum possible renewable generation without any curtailment. Therefore, in the case study, the power output of wind and PV plants will be used as given by (29) and (30), respectively. In abnormal circumstances, if necessary, the power generated by the renewables can be curtailed to meet generation-demand balance on the electrical side. The renewable curtailment, however, is not considered in this paper.

Resilience Framework
Any disturbance in a power system will require either an increase or decrease in power as a response from generation assets. For the purposes of this paper, the events requiring power increase are referred to as ∆P + , and the events requiring power decrease are referred to as ∆P − . In the proposed resilience framework, the real-power flexibility of microgrid assets is evaluated and aggregated to obtain overall adaptive capacity of a microgrid to respond to ∆P + and ∆P − . Note that the proposed framework does not evaluate the microgrid's response against the classes of disturbances or the performance of the microgrid following a disturbance. Instead, it identifies the adaptive capacity of the microgrid to respond to generation-demand imbalance created by any means, including large disruptive events and DER generation uncertainty.
The resilience framework is developed for the SMR-DER integrated-energy microgrid considering two operating scenarios: normal operation (NO) and compromised operation (CO). The NO includes operating conditions when the flexibility needs of a microgrid can be fulfilled by controlling assets normally based on the operational framework described in Section 3.1. The CO refers to abnormal operating conditions when NO of assets does not meet the flexibility needs of the microgrid, and additional flexibility schemes are activated.
During NO, the disturbances are typically smaller such that the assets operating normally will be able to handle them. BES alone can respond to immediate real-power imbalance without compromising the dynamic performance. Once the system reaches the start of the next LF cycle (i.e., LF decision point), steam extraction system provides additional flexibility. Once the system reaches the start of the next CLS cycle (i.e., CLS decision point), further flexibility will be added by the reactor's power maneuvering capability. The steam extraction system will be able to provide maximum LF response as long as the TES temperature is within minimum and maximum limits. The SOC levels of TES and BES are important factors to meet the flexibility requirement of a system against disturbances that last longer in the system. In NO, the reactor must comply with the ramp rate and total power change limits, as given by (11) and (12), respectively. Figure 5 illustrates the asset operation following the disturbance during NO along with flexibility limits for different assets in the system. Note that the figure is just an illustration and does not reflect the time and sequence of asset activation and the magnitude of variables in scale. The time of asset activation and the system adaptive capacity depends on the system's dynamic state at the onset of a disturbance.   Figure 6 shows the activation of additional flexibility during CO. During CO, the system would be prioritized to meet the demands of critical loads. For ∆P + , additional flexibility is added by curtailing non-critical electrical loads which can be activated im-mediately following a disturbance. Curtailing thermal loads would also be useful if the SOC level of the TES drops below the minimum limit. For ∆P − , the turbine bypass system can bypass excess steam to the condenser. The condenser can receive 100% of the process steam but only for a limited time. Once the system reaches the LF decision point, further flexibility could be added with curtailable thermal loads. Because the extraction steam has much less enthalpy compared to the directly bypassed steam, the full steam extraction capability can be used as long as needed. In other words, the maximum amount of steam that can be extracted to the heating side in response to ∆P − is not limited by the TES SOC. During CO, unlike during NO, the reactor power level of SMR can be maneuvered up to its maximum or minimum limits irrespective of the power change limit given by (12). In the case of ∆P − , if the generation-demand imbalance sustains even after continuous bypassing of steam for a certain allowed duration, the reactor will isolate from the microgrid and initiate emergency shutdown. This operating event is termed as emergency operation (EO) which arises due to large disturbances, such as sustained faults and total loss of electrical connection of SMR to microgrid. The EO is an event that originates out of CO rather than a separate operational mode. The ∆P + , on the other hand, may force the system to compromise the power supply to critical loads but not the reactor safety itself.
The event tree in Figure 7 shows the decision-making process following a disturbance. Based on the magnitude and form of the disturbance, the microgrid activates its flexible assets to respond to the generation-demand imbalance created. The microgrid switches between these operating modes in real time to respond to power fluctuation created by the disturbance. The seamless transition between the operational modes in real time is crucial. At the onset of disturbance, BES starts responding with the system assuming NO. If BES is unable to meet generation-demand balance, additional assets such as turbine bypass system or load curtailment activate switching system from NO to CO. If the disturbance dies out sooner, the system will first withdraw the assets that were activated as a part of CO followed by the assets of NO. In this way, the system will smoothly transition between NO and CO. Based on temporal characteristics of the disturbance, the system can switch between CO and NO at any time. To facilitate this seamless transition, the assets are decoupled to not disrupt the operational pattern of one another. For example, switching to CO does not allow the system to operate steam extraction or reactor control before the decision point. Rather, an equivalent change is achieved by other independent methods such as turbine bypass system or load curtailments. This allows operation at pre-optimized setpoints following the return of the system to NO after the disturbance dies out. Once the reactor shutdown sequence (EO mode) is activated, the system will not be able to recover. The microgrid should therefore continue operating in CO mode long enough to evaluate if the switching to EO mode is absolutely necessary.  In the proposed resilience model, the asset-level flexibility shown in Figures 5 and 6 are aggregated to obtain the adaptive real-power capacity of the microgrid in NO and CO mode. The adaptive capacity evaluated in NO and CO modes are analyzed to identify the disturbance threshold for microgrid control to facilitate quicker and efficient decisionmaking. Based on this threshold measurement, the microgrid selects between NO or CO mode, as shown in Figure 7. In actual operation, the resilience information obtained in real time can also be used to supervise the microgrid management system to reconfigure its assets to maximize operational resilience.

Case Study
A microgrid system is developed in a portion of the IEEE-30 bus network hosting generation in the form of an SMR plant and renewables. The district heating (DH) system is also integrated into the SMR plant. First, the case system is discussed along with the simulation results which include the 720 h quasi-static simulation and a 2.5 h dynamic simulation. The result is then shown in analyzing the microgrid's operational resilience in terms of adaptive real-power capacity for normal, compromised, and emergency operations. Note that the microgrid model and the operation simulation in this case study are replicated from [20]. Readers are highly encouraged to refer to [20,21] for background details on the case system and operational setup. In this paper, the system details and operational simulation results are summarized before moving into the resilience evaluation part.

Modified IEEE-30 Bus System
The IEEE-30 bus system consists of 33 kV and 230 kV transmission network [48]. The 33 kV portion can be isolated from the 230 kV system by taking four 230/33 kV transformers out of service. The resulting 33 kV network is used as a case system to develop the microgrid model for this paper. Figure 8 shows the proposed case system which has total electrical load of 102 MWe. The proposed case system does not have any generation source. Therefore, to meet the system demand, new power generation sources are needed. A microgrid with a cogenerating SMR plant, wind, PV, BES, and TES is suitable for the system. An SMR plant with two NuScale integral pressurized-water reactor (iPWR)-type SMRs, rated 50 MWe each, is hosted at node 10. A wind plant, sized 30 MWe, is hosted at node 15, and two PV plants, 10 MWe each, are hosted at node 17 and node 27, respectively. A centralized BES is hosted at node 10 to support SMR for short-term response. The system also has two synchronous condensers (SC) at nodes 11 and 13 to provide reactive power support. The DH system is sized based on the size of electrical system and lumped as a single load. The TES is collocated with the SMR plant along with HXs, as shown in Figure 3. Reference [21] provides the optimum sizing of generation assets based on operational and economic constraints including reactor maneuvering limits, heat waste in DH system, cost of BES, and power system dynamic performance. The sizing of assets is summarized in Table 1. The loads centers in Figure 8 are populated with residential, industrial, and commercial load profiles of 15 min resolution obtained from [49]. The heat demand profile is obtained by scaling practical heat load profile from [50] with 15 min resolution. The wind and PV generation profiles are obtained from [51][52][53], respectively.

Operational Results
In this paper, the operational simulation in [21] is replicated with the resilience framework functionality integrated in the microgrid model. In this paper, the resilience model is used for system observation rather than to execute any control actions. Therefore, the operational results will be identical to the ones in [21]. The description below will focus only on results essential for the resilience evaluation. The quasi-static simulation is performed for a month with a 15 min simulation timestep. The case system of Figure 8 developed in Power System Simulator Siemens Network Calculation (PSS/Sincal) [54] is cosimulated with the quasi-static model of SMR units and DH system developed in MATLAB. It captures the CLS and LF controls of the microgrid. The dynamic simulation is performed in Power System Simulator for Engineering (PSS/E) [55] for the 2.5 h period which contains the largest fluctuations in the system to analyze FC control and its interaction with CLS and LF controls. Figure 9 shows the results of quasi-static simulation. The electrical load variations along with the intermittency of wind and PV are the source of fluctuations on the electrical side. Figure 9a shows the demand at the SMR bus (Bus 10) which is taken as the slack bus for power flow simulations. On a daily basis, the PV generation goes through zero-power output periods during the nighttime. During the daytime, the PV generation faces large fluctuations due to the cloud movement covering the direct irradiance to the solar panels. Similarly, the wind power fluctuates based on the wind speed. The power outputs of wind and PV plants vary as high as 47% and 61%, respectively, of their rated output within a minute. This integration of wind and PV generation results in large power fluctuations in the microgrid which are managed by controlling the power outputs of flexible assets. As discussed earlier, the reactor power is controlled a fixed number of times per day based on the CLS algorithm. In the proposed scheme, usually only one unit is controlled in a CLS interval while another one remains on hold to reduce the total number of reactor power maneuvers. Out of the 120 CLS intervals, both units were required to provide power maneuvering on one occasion to meet large demand fluctuations. Similarly, for three CLS intervals, none of the units were required to change their power levels. Figure 9c shows the reactor power level for the two units over the 720 h duration that included 120 CLS intervals and total of 116 reactor power maneuvers.  Figure 9b shows DH steam extraction as a ratio of the total steam generated and the SMR plant's instantaneous electrical power in per unit (pu) of its rated electrical output. As seen in Figure 9a,b, when the power output from wind and PV plants increases, SMR's demand decreases. In response, the extraction flow to the DH system is increased to reduce the electrical output of SMR. The extraction steam is utilized to meet the heating demands with TES absorbing the surplus steam or supplying the deficient heat to the DH pipeline. Due to this, the TES temperature fluctuates during the simulation as shown in Figure 9d. For 2.4% of time, the TES was at upper limit and was not able to accept surplus steam, and the extraction steam excess to the DH system was directly bypassed to the condensor. The DH temperature never went below the lower limit and was always able to supply the deficient heat to the DH pipeline when the steam extracted was not sufficient to fulfill the DH demand.
In quasi-static simulation, BES and turbine bypass valves remain inactive, as they are operated to provide primary response to the frequency deviations. The quasi-static simulation only shows the CLS and LF controls. It does not account for the smaller fluctuations present within each LF intervals that occur due to short-term wind and PV intermittency.
The dynamic simulation is performed for a 2.5 h time segment of quasi-static simulation which includes largest fluctuations from quasi-static simulation (i.e., between 298.5 to 301 h). This 2.5 h period includes the reactor ramping period of the CLS control. Both reactor units were ramped up between 298.75 to 300.75 h as shown in Figure 10b. It also includes 10 consecutive LF controls which are repeated in the intervals of 15 min. Figure 10a plots the electrical power outputs of the two SMR units along with the load demand and BES power output. BES is operated to respond to the frequency deviation occurring due to generation-demand imbalance persisting after the CLS and LF control. When BES alone is not able to meet the FC requirements, the turbine bypass system operates as shown in Figure 10c.

Resilience Evaluation
The resilience model takes system dynamic states as inputs and analyzes the asset's capabilities at each simulation timestep. Looking forward in time, the magnitude and temporal constraints of microgrid assets are analyzed to obtain the flexibility characteristics that provide the magnitude shape of the real-power variation they can offer for a certain amount of time. These flexibility characteristics of assets are then aggregated to obtain the overall adaptive capacity characteristics of the microgrid.
The graph showing the temporal characteristics of asset-level flexibility and overall adaptive capacity is referred to as the resilience frame. Resilience frames are generated in real time at each simulation timestep. The time axis of a resilience frame begins at 0 s, corresponding to the simulation time at which resilience is being analyzed, and lasts until the end of subsequent CLS interval, assuming no substantial flexibility can be added after CLS is executed. Note that the resilience frame is an exploration of a future response to predict the adaptive capacity of microgrid in real time at each simulation timestep. The time axis of resilience frame is not a part of the actual operational simulation being executed.
Logarithmic scale is used for time axes of the resilience frame to illustrate the roles of different assets that differ significantly in terms of their activation and operation times. For example, BES can ramp up to maximum power in either direction within 20 s, whereas steam extraction system takes 15 min and reactor ramping takes 2 h. Logarithmic scale adequately represents operation of assets with significantly different temporal characteristics in a single graph.
The resilience model is integrated into both quasi-static and dynamic simulation models of the microgrid. Sample calculations are shown for three different cases: two for quasi-static simulation and one for dynamic simulation. The system resilience is then quantified in terms of the normalized area under the adaptive capacity curve, termed as RAM, for both NO and CO modes. The ∆P + and ∆P − flexibilities are analyzed separately for each operational mode. Note that the proposed resilience model calculates the maximum adaptive capacity instead of analyzing the response to a disturbance. Therefore, the asset will be operated up to their limits allowed in each operational mode. Table 2 provides the description of variables shown in the resilience frames in Figures 11-13.

Variable Description
Rx_t_pu, x = 1,2 Reactor thermal power level of unit x in pu Rx_e_pu, x = 1,2 Reactor electrical power level of unit x in pu Rx_ex_pu, x = 1,2 Extraction steam flow rate of unit x in pu Rx_bp_pu, x = 1,2 Bypass steam flow rate of unit x in pu PB_e_MWe Battery electrical power output in MWe PD_crt_MWe Electrical load curtailed in MWe QD_crt_MWt Thermal load curtailed in MWt P_ad_MWe Microgrid adaptive real-power capacity in MWe The first case selected is 0 h of the simulation. At 0 h, TES is at 90 • C, and BES SOC is at 0.65 pu. Figure 11a-d shows the resilience frames for four cases: ∆P + and ∆P − for NO and CO. The flexibility offered by different assets and their aggregated, in terms of the adaptive real-power capacity changing with time, are shown in each resilience frame. The BES starts ramping up or down immediately and lasts until it reaches its upper or lower SOC limits (given in Table 1). For ∆P + , after BES reaches its lower SOC limit, it can no longer discharge, and BES power returns to zero at 1306 s, shown in Figure 11a,c. Similarly, for ∆P − , BES power returns to zero after reaching the upper SOC limit at 1144 s, shown in Figure 11b,d. The LF decision point arrives at 900 s, and the CLS decision point arrives at 2700 s. At the LF decision point, the steam extraction valves can be controlled to their limits in either direction (i.e., completely closed extraction valve or 0 kg/s extraction for ∆P − ) and completely open extraction valve or fullṁ SG t extraction for ∆P − . Similarly, both reactor units are allowed to maneuver at once but up to the reactor power change limits (see Table 1). The reactor power ramping completes at 9900 s, and the reactor stays on hold until the next CLS decision point at 24,300 s. The TES did not fully discharge over the duration considered, and, therefore, the flexibility offered by the steam extraction for ∆P + remained over the resilience frame.
For ∆P + in CO, the non-critical electrical loads were curtailed immediately to meet the short-term response. For ∆P − in CO, the turbine bypass system operates which brings electrical power output of SMRs to zero, as shown in Figure 11d. After two LF cycles, the turbine bypass system is withdrawn. This significantly reduces the adaptive capacity in CO mode. Although not analyzed in this paper, if a power imbalance persists even after bypassing of steam for two LF cycles, the reactor will initiate the emergency shutdown. For ∆P − in CO, the steam extraction can be fully operated irrespective of the TES temperature limits as the excess low-temperature extracted steam can be continuously bypassed to the condenser.
The second case is for 216.75 h of the simulation, as shown in Figure 12a-d. Similar to the first case, the BES is activated immediately for NO, and electrical load curtailment or turbine bypass system are activated for CO mode. Steam extraction system and reactor maneuvering both started at 900 s after the disturbance. Initially, BES SOC was 0.65 pu, and the TES temperature was 74.15 • C. This case with low starting TES temperature is selected to demonstrate how the quick depletion of TES SOC will impact the operational resilience on the power side. In this case, the TES temperature reaches its lower limit of 70 • C at 9900 s. After TES reached its lower temperature, the steam extraction system could not provide LF control for ∆P + in NO mode. This results in sharp reduction of adaptive real-power capacity at 9900 s, as shown in Figure 12a. For ∆P + in CO mode, the non-critical thermal load is curtailed to provide additional flexibility on the electrical side, as shown in Figure 12b. All other results are similar to the first case.   The third case shown in Figure 13a-d is for 300.97769 h of the simulation where the resilience model uses system state values directly from dynamic simulation. Initially, BES SOC was at 0.463 pu and TES at 92.72 • C. Compared to the second case, BES discharged quicker for ∆P + . Similarly, the bypass valve flows were non-zero at the start: 0.182 pu/s (i.e., 12.11 kg/s) for unit 1 and 0.129 pu/s (i.e., 8.59 kg/s) for unit 2. For ∆P + , the bypass flows came back to zero, as shown in Figure 13a,c, whereas for ∆P − , the bypass flows stayed at the initial non-zero value for NO, as shown in Figure 13b. For ∆P − in CO, the bypass flows ramp up to maximum, as shown in Figure 13d. All other results are similar to the earlier two cases. Figures 11-13 show the use of proposed resilience frame to evaluate system adaptive real-power capacity. These three cases show that the resilience model can be used with both quasi-static and dynamic simulations. The adaptive real-power capacity changes over time due to ramping delays of the assets, activation delays of new assets, and exhaustion of energy-limited assets. The adaptive capacity curves produced in Figures 11-13 give an idea of the disturbance characteristics the microgrid can respond to in different modes of operation.  Figure 14 maps the adaptive capacity curve to the DIRE curve for the first case of 0 h to analyze the Rs of resilience (i.e., resist, response, and recover). Although all assets are represented in terms of their capacity to respond to disturbances, they have different principal functions. BES, for example, is operated immediately and helps to arrest disturbance and resist the dynamic performance degradation. Therefore, a part of the BES response will fall under the resist epoch which is considered in the first 10 s of the resilience frame for the purpose of this paper. It is the time when a disturbance typically grows. The steam extraction system operates to maintain generation-demand balance for steady-state operation. It falls mainly during the response epoch between 10-1000 s during which disturbance reaches its peak and persists. The reactor power control with CLS falls under recover epoch (>1000 s), and it helps recover steam extraction setpoint and energy lost of energy-limited assets while responding to the disturbance. As discussed earlier, the adaptive capacity is a temporal characteristic rather than a single point value. Therefore, it is not computationally viable to analyze the whole characteristics in real-time decision-making. To compare the system resilience over time, a metric is introduced quantifying the temporal change of adaptive real-power capacity with a single variable. The metric termed as RAM is given as

RAM =
Area under adaptive capacity curve Total duration = where P − ad − MWe is time-varying system adaptive capacity as obtained from the resilience frame, and t FF is the total duration of the resilience frame. RAM represents the average adaptive capacity over the resilience time frame. RAM evaluated for simulation period 0 h to 720 h is plotted in Figure 15. It includes RAM values for each quasi-static simulation time steps and 2.5 h dynamic simulation period. The RAM values for the 2.5 h dynamic simulation are plotted in the zoomed section in Figure 15b. The RAM values provide information to the system regarding its capacity to respond to any disturbance that may arise. Figure 15a,b shows the calculation of RAM values in real time, parallel to the simulation. This information can be provided as a feedback to the microgrid management system in real time to execute precautionary actions such as reconfiguration of assets to improve operational resilience. Figure 15 clearly shows that the resilience values are inconsistent over simulation duration. For example, the RAM values in NO for ∆P + reached its lowest of 6.05 MWe at 258.75 h, whereas for ∆P − , it was −6.56 MWe at 300.9656 h. If the system were to face disturbances larger than these values, it may have to switch to CO. In CO, RAM was a minimum of 28.92 MWe at 258.75 h for ∆P + and −13.99 MWe for ∆P + at 298.7347 h.
Although RAM provides high-level information on system resilience as discussed, it has some limitations. Firstly, it does not properly credit the contributions of short-term assets. The short-term assets such as BES will not have a large area contribution to the RAM due to energy limitation, but they are crucial to arrest the disturbance immediately after its inception. Secondly, RAM also loses the temporal feature of adaptive capacity curve which is crucial to analyze dynamic characteristics of disturbances. Therefore, in this paper, we intend to present RAM as a comparison metric for studies and an initial guidance for metric development rather than a final product. The next steps are to resolve these limitations of the RAM model, identify standard limits for RAM for comparison in real-time operation, and develop control algorithms for microgrid to take precautionary actions based on RAM feedbacks.

Conclusions and Future Work
Amid increasing interest in SMR use in microgrids, this paper analyzed SMR salient features and modeled the interaction between heat and power systems to develop a framework that evaluates adaptive capacity for the resilient operation of SMR-DER integratedenergy microgrids. An operational scheme is designed for microgrid control that decoupled overall microgrid operation into three separate controls of CLS, LF, and FC. In addition to normal operation, operational algorithms were discussed for CO. The proposed resilience model aggregates asset-level real-power flexibility to evaluate the operational resilience indicator in terms of adaptive real-power capacity characteristics. The adaptive capacity curve produced has both time and magnitude characteristics which reflect the operational margin of the system to respond to generation-demand imbalance created by any form of disturbances, including large disruptive events and generation uncertainties. A single point operational metric known as RAM was also obtained from the adaptive capacity characteristic to compare the operational resilience over time and provide high-level guidance to microgrid operators. The quasi-static simulation results for 720 h and the dynamic simulation results for 2.5 h show that the system with the proposed operational framework is resilient enough to withstand and operate normally with DER penetration level as high as 50% (50 MWe, wind, and PV collectively). The RAM results shown reflect the system resilience in real time while operating with the proposed operational framework. The case study clearly shows that the proposed resilience model can improve the operational awareness of the microgrid and help optimize the system control.
The adaptive capacity framework, previously limited to resilience modeling of power system components, is extended in this paper to characterize the impact of nuclear and heating systems on overall system resilience. The proposed model improves the awareness of power side to the nuclear and heating sides and vice versa for resilience planning, operational decision-making, and real-time control. Future work planned with this resilience framework includes incorporating the inertial response as a part of the resist epoch, including more detailed operational procedure of SMRs based on data from manufacturers and operators, including curtailment models of wind and PVs to support resilience, synthetic inertia by inverter-based resources in microgrids, and analyzing the reactive power capabilities of system assets. Further work is necessary to resolve the limitations with the RAM model and develop its standards. Future work will also demonstrate the application of the proposed metric for real-time system control with hardware-in-loop or software-in-loop simulations.

Data Availability Statement:
The data presented in this study is available on request from the corresponding author. It is not publicly available due to its large size. The description of the case system and parameters necessary for obtaining the presented results using commercially available simulation softwares are provided in this paper.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: