Book Description

This report documents work carried out during FY 2008 on further investigation of control strategies for supercritical carbon dioxide (S-CO2) Brayton cycle energy converters. The main focus of the present work has been on investigation of the S-CO2 cycle control and behavior under conditions not covered by previous work. An important scenario which has not been previously calculated involves cycle operation for a Sodium-Cooled Fast Reactor (SFR) following a reactor scram event and the transition to the primary coolant natural circulation and decay heat removal. The Argonne National Laboratory (ANL) Plant Dynamics Code has been applied to investigate the dynamic behavior of the 96 MWe (250 MWt) Advanced Burner Test Reactor (ABTR) S-CO2 Brayton cycle following scram. The timescale for the primary sodium flowrate to coast down and the transition to natural circulation to occur was calculated with the SAS4A/SASSYS-1 computer code and found to be about 400 seconds. It is assumed that after this time, decay heat is removed by the normal ABTR shutdown heat removal system incorporating a dedicated shutdown heat removal S-CO2 pump and cooler. The ANL Plant Dynamics Code configured for the Small Secure Transportable Autonomous Reactor (SSTAR) Lead-Cooled Fast Reactor (LFR) was utilized to model the S-CO2 Brayton cycle with a decaying liquid metal coolant flow to the Pb-to-CO2 heat exchangers and temperatures reflecting the decaying core power and heat removal by the cycle. The results obtained in this manner are approximate but indicative of the cycle transient performance. The ANL Plant Dynamics Code calculations show that the S-CO2 cycle can operate for about 400 seconds following the reactor scram driven by the thermal energy stored in the reactor structures and coolant such that heat removal from the reactor exceeds the decay heat generation. Based on the results, requirements for the shutdown heat removal system may be defined. In particular, the peak heat removal capacity of the shutdown heat removal loop may be specified to be 1.1 % of the nominal reactor power. An investigation of the oscillating cycle behavior calculated by the ANL Plant Dynamics Code under specific conditions has been carried out. It has been found that the calculation of unstable operation of the cycle during power reduction to 0 % may be attributed to the modeling of main compressor operation. The most probable reason for such instabilities is the limit of applicability of the currently used one-dimensional compressor performance subroutines which are based on empirical loss coefficients. A development of more detailed compressor design and performance models is required and is recommended for future work in order to better investigate and possibly eliminate the calculated instabilities. Also, as part of such model development, more reliable surge criteria should be developed for compressor operation close to the critical point. It is expected that more detailed compressor models will be developed as a part of validation of the Plant Dynamics Code through model comparison with the experiment data generated in the small S-CO2 loops being constructed at Barber-Nichols Inc. and Sandia National Laboratories (SNL). Although such a comparison activity had been planned to be initiated in FY 2008, data from the SNL compression loop currently in operation at Barber Nichols Inc. has not yet become available by the due date of this report. To enable the transient S-CO2 cycle investigations to be carried out, the ANL Plant Dynamics Code for the S-CO2 Brayton cycle was further developed and improved. The improvements include further optimization and tuning of the control mechanisms as well as an adaptation of the code for reactor systems other than the Lead-Cooled Fast Reactor (LFR). Since the focus of the ANL work on S-CO2 cycle development for the majority of the current year has been on the applicability of the cycle to SFRs, work has started on modification of the ANL Plant Dynamics Code to allow the dynamic simulation of the ABTR. The code modifications have reached the point where a transient simulation can be run in steady state mode; i.e., to determine the steady state initial conditions at full power without an initiating event. The results show that the steady state solution is maintained with minimal variations during at least 4,000 seconds of the transient. More SFR design specific modifications to the ANL Plant Dynamics Code are required to run the code in a full transient mode, including models for the sodium pumps and their control as well as models for reactivity feedback and control of the reactor power.