Search NASA⌕ Search

SEARCH · Search NASA

Results for “grid simulator testing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,315 records · Page 73

Real-Time Digital Simulation of Inertial Response with Hardware-in-the-Loop Implementation on the CART3 Wind Turbine at the National Wind Technology Center

With increasing penetrations of wind power on electric grids, the stability and reliability of interconnected power systems may be impacted. In some countries that have developed renewable energy sources and systems, grid codes have been revised to require wind power plants (WPPs) to provide ancillary services to support the power system frequency in case of severe grid events. To do this, wind turbine generators (WTGs) should be deloaded to reserve a certain amount of active power for primary frequency response; however, deloading curtails annual energy production, and the market for this type of service needs to be further developed. In this report, we focus on the temporary frequency support provided by WTGs through inertial response. WTGs have potential to provide inertial response, but appropriate control methods should be implemented. With the implemented inertial control methods, wind turbines are capable of increasing their active power output by releasing some of their stored kinetic energy when a frequency excursion occurs. Active power can be temporarily boosted above the maximum power points, after which the rotor speed decelerates, and subsequently an active power output reduction restores the kinetic energy. In this report, we develop two types of models for wind power systems: the first is common, based on the wind power aerodynamic equation, and the power coefficient can be regressed using nonlinear functions; the second is much more complicated, wherein the wind turbine system is modeled using the Fatigue, Aerodynamics, Structures, and Turbulence Modeling (FAST) tool with several degrees of freedoms. A nine-bus test power system is built in Simulink and the Real-Time Digital Simulator, respectively, which are used to evaluate the frequency support performance of the WPPs. We implement two distinct types of inertial control methods in the modeled wind turbines: frequency-based inertial control (FBIC) and stepwise inertial control (SIC). We compare the performances of the two methods in terms of their frequency nadirs, rates of change of frequency, and recovery times. We conclude the results under various wind speeds and penetration cases, which provide insight into designing the inertial response of WTGs. Further, we discuss the impact of the parameters on the performance of the inertial control methods. We evaluate both the scaling factors for the FBIC method and the slope values for the TLIC methods. The simulation work shows the characteristics of different inertial responses compared to conventional synchronous generators. Based on the simulation results, we modify, improve, and test the inertial control methods under a more realistic wind turbine model based on FAST. We then validate the inertial responses under highly turbulent wind conditions generated by TurbSim, and we examine their influences on the turbine mechanical components. The extensive simulation proves the effectiveness of the proposed inertial control methods as well as the nine-bus test power system. We then reconsider the parameters. We rebuild the same test power system using Real time Simulator Computer Aided Design (RSCAD), and we implement the inertial control methods in the real Controls Advanced Research Turbine (CART3), which is prepared for the hardware-in-the-loop field-test simulation. After the setups for the hardware and software hybrid simulation platform are complete, the inertial response is further tested on a real wind turbine for the first time, in which CART3 release the controlled inertial response against the emulated frequency excursion, provided by the real-time simulated power system test bed in RTDS.

17 WIND ENERGY↗

A purely Lagrangian method for simulating the shallow water equations on a sphere using smooth particle hydrodynamics

It has long been suggested that a purely Lagrangian solution to global-scale atmospheric/oceanic flows can potentially outperform tradition Eulerian schemes. Meanwhile, a demonstration of a scalable and practical framework remains elusive. Motivated by recent progress in particle-based methods when applied to convection dominated flows, this work presents a fully Lagrangian method for solving the inviscid shallow water equations on a rotating sphere in a smooth particle hydrodynamics framework. To avoid singularities at the poles, the governing equations are solved in Cartesian coordinates, augmented with a Lagrange multiplier to ensure that fluid particles are constrained to the surface of the sphere. An underlying grid in spherical coordinates is used to facilitate efficient neighbor detection and parallelization. The method is applied to a suite of canonical test cases, and conservation, accuracy, and parallel performance are assessed..

Computer Science↗

High Order Anchoring and Reinitialization of Level Set Function for Simulating Interface Motion

A second order interface anchoring method has been developed and used with fast sweeping algorithm for reinitialization of a level set function. The algebraic anchoring formulation ensures that the location of the actual interface is preserved, leading to better mass conservation property. It also provides high order accurate algebraic constraint for solving the Eikonal equation on a finite difference grid. Geometric properties of the interface such as surface normal and curvature are subsequently computed from the reinitialized distance function. Various analytical functions for modeling distortion in level set field are considered and accuracy of reinitialization is evaluated using first and second order anchoring schemes. It is also shown that accurate computation of interface curvature requires a high order anchor in addition to a high order fast sweeping method. Mass conservation property of reinitialization is also analyzed by considering test problems from literature including the classic Rider–Kothe single vortex problem. The formulation is suitable for efficient parallelization for both distributed memory and on-node shared memory parallel systems. In conclusion, scalability and performance of the reinitialization scheme on multiple architectures are demonstrated.

97 MATHEMATICS AND COMPUTING↗

A weighted shifted boundary method for immersed moving boundary simulations of Stokes' flow

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. The surrogate domain is constructed so as to avoid cut cells and the associated problematic implementation and numerical integration issues. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions: hence the name of the method, that shifts the location and values of the boundary conditions. Here, in this article, we extend the SBM to the simulation of incompressible Stokes flow, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach allows to drastically reduce spurious pressure oscillations in time, which are produced if the total volume of active fluid were to change abruptly over a time step. The proposed Weighted SBM (W-SBM) exactly preserves states of hydrostatic equilibrium, and induces small mass and momentum conservation errors, which converge as the grid is refined. This is in analogy to cutFEMs and related unfitted approaches, which rely on an affine representation of cut boundaries. We demonstrate the robustness and accuracy of the proposed method with an extensive suite of two-dimensional tests.

97 MATHEMATICS AND COMPUTING↗

Flexible FlueCO2

Carbon dioxide (CO2) emission reductions remain a significant challenge on the path to clean energy. There are increasing legislative, social, and environmental factors motivating CO2 emissions reduction from power plants with carbon capture and storage (CCS). CCS in natural gas combined cycle (NGCC) power plants is critical to achieve a net-zero carbon electricity grid. Enhanced 45Q tax credits provide new incentives, but currently available technologies are unable to profitably operate in grids with deep variable renewable penetration which require flexible NGCC operation. Luna Labs has developed the FlueCO2 membrane to enable a profitable NGCC-CCS process. The FlueCO2 membrane couples steam transport across the membrane to CO2 transport in the opposite direction, enabling high capture efficiencies and low energy costs even at low CO2 concentrations. The dual-phase membrane can operate in the range of typical flue gas temperatures and pressures and does not require temperature or pressure cycling. Luna Labs’ FlueCO2 technology enables flexible and profitable operation of NGCC plants with lower capital investment and impact on electricity prices. In this Phase 1 project, Luna Labs utilized experimental testing, modeling, process simulation, and standardized costing methodologies to evaluate the techno-economic value of a 650 MW greenfield NGCC plant with FlueCO2 (NGCC-FlueCO2). Key design requirements for operation were established and plant performance under load-leveling conditions was validated through computational fluid dynamics and process modeling. Luna Labs developed a dynamic modeling tool which modeled plant operational modes across a variety of tax structures and electricity pricing scenarios to project the overall Net Present Value (NPV) of the NGCC-FlueCO2. FlueCO2 minimizes the impact of CCS integration on plant operation by integrating directly into the NGCC heat recovery steam generator (HRSG). By tapping into the plant’s low-pressure (LP) steam, operators can divert LP steam to the greenfield NGCC and/or CCS process in response to dynamic markets. Since FlueCO2 will not significantly affect HRSG (or NGCC) operation, CCS only turns off during peak power demand (>$250/MWh). Under baseload conditions, FlueCO2 lowers the capital (37%), energy (36%) and carbon capture (<$40/tonne) costs and can increase the overall plant lifetime NPV by approximately ~$1B in comparison with NGCC solvent-based capture reference cases (NETL Case 31B). Luna Labs has shared its costing tools with several interested partners and customers, which follows a generalizable approach to costing analysis.

Kelly, Jesse↗

A framework to analyze the stochastic harmonics and resonance of wind energy grid interconnection

This study addresses a modeling and analysis methodology for investigating the stochastic harmonics and resonance concerns of wind power plants (WPPs). Wideband harmonics from modern wind turbines are observed to be stochastic, associated with real power production, and they may adversely interact with the grid impedance and cause unexpected harmonic resonance if not comprehensively addressed in the planning and commissioning of the WPPs. These issues should become more critical as wind penetration levels increase. We thus propose a planning study framework comprising the following functional steps: First, the best-fitted probability density functions (PDFs) of the harmonic components of interest in the frequency domain are determined. In operations planning, maximum likelihood estimations followed by a chi-square test are used once field measurements or manufacturers' data are available. Second, harmonic currents from the WPP are represented by randomly-generating harmonic components based on their PDFs (frequency spectrum) and then synthesized for time-domain simulations via inverse Fourier transform. Finally, we conduct a comprehensive assessment by including the impacts of feeder configurations, harmonic filters, and the variability of parameters. We demonstrate the efficacy of the proposed study approach for a 100-MW offshore WPP consisting of 20 units of 5-MW full-converter turbines, a realistic benchmark system adapted from a WPP under development in Korea, and discuss lessons learned through this research.

17 WIND ENERGY↗

Computation of incompressible viscous flows through artificial heart devices with moving boundaries

The extension of computational fluid dynamics techniques to artificial heart flow simulations is illustrated. Unsteady incompressible Navier-Stokes equations written in 3-D generalized curvilinear coordinates are solved iteratively at each physical time step until the incompressibility condition is satisfied. The solution method is based on the pseudo compressibility approach and uses an implicit upwind differencing scheme together with the Gauss-Seidel line relaxation method. The efficiency and robustness of the time accurate formulation of the algorithm are tested by computing the flow through model geometries. A channel flow with a moving indentation is computed and validated with experimental measurements and other numerical solutions. In order to handle the geometric complexity and the moving boundary problems, a zonal method and an overlapping grid embedding scheme are used, respectively. Steady state solutions for the flow through a tilting disk heart valve was compared against experimental measurements. Good agreement was obtained. The flow computation during the valve opening and closing is carried out to illustrate the moving boundary capability.

Kiris, Cetin↗

Experimental Hypersonic Aerodynamic Characteristics of the 2001 Mars Surveyor Precision Lander with Flap

Aerodynamic wind-tunnel screening tests were conducted on a 0.029 scale model of a proposed Mars Surveyor 2001 Precision Lander (70 deg half angle spherically blunted cone with a conical afterbody). The primary experimental objective was to determine the effectiveness of a single flap to trim the vehicle at incidence during a lifting hypersonic planetary entry. The laminar force and moment data, presented in the form of coefficients, and shock patterns from schlieren photography were obtained in the NASA Langley Aerothermodynamic Laboratory for post-normal shock Reynolds numbers (based on forebody diameter) ranging from 2,637 to 92,350, angles of attack ranging from 0 tip to 23 degrees at 0 and 2 degree sideslip, and normal-shock density ratios of 5 and 12. Based upon the proposed entry trajectory of the 2001 Lander, the blunt body heavy gas tests in CF, simulate a Mach number of approximately 12 based upon a normal shock density ratio of 12 in flight at Mars. The results from this experimental study suggest that when traditional means of providing aerodynamic trim for this class of planetary entry vehicle are not possible (e.g. offset c.g.), a single flap can provide similar aerodynamic performance. An assessment of blunt body aerodynamic effects attributed to a real gas were obtained by synergistic testing in Mach 6 ideal-air at a comparable Reynolds number. From an aerodynamic perspective, an appropriately sized flap was found to provide sufficient trim capability at the desired L/D for precision landing. Inviscid hypersonic flow computations using an unstructured grid were made to provide a quick assessment of the Lander aerodynamics. Navier-Stokes computational predictions were found to be in very good agreement with experimental measurement.

Horvath, Thomas J.↗

Modeling and Analysis of Stirling Power Convertors

Modeling and Analysis of Stirling Power Convertors Luis A. Rodriguez1 Steven M. Geng, Terry V. Reid, Scott D. Wilson NASA Glenn Research Center, Cleveland, OH, 44135, USA NASA Glenn’s Thermal Energy Conversion Branch is supporting the development of the next generation free-piston Stirling power convertors. American Superconductor (AMSC) and Sunpower Inc. are the two firms under contract to develop the Flexure Isotope Stirling Convertor (FISC) and the Sunpower Robust Stirling Convertor (SRSC), respectively. To comprehend and forecast convertor performance, Sage, ANSYS® Maxwell, and ANSYS® Fluent were used to model the Stirling thermodynamic cycle, alternator electromagnetics, and piston and displacer dynamics. I. Introduction Stirling convertors are being developed by NASA as a potential steady source of electrical power for NASA’s future scientific space missions. Currently, NASA Glenn Research Center has two corporations under contract, American Superconductor (AMSC) and Sunpower Inc., for the development of the next generation of free-piston Stirling convertors for dynamic radioisotope power systems. AMSC is developing the Flexure Isotope Stirling Convertor (FISC), which uses flexures to prevent side motion and rubbing of the piston. Similarly, Sunpower Inc, is developing the Sunpower Robust Stirling Convertor (SRSC). The SRSC uses gas bearings to prevent radial contact of the moving piston. As convertor development continues, it is increasingly important to understand and predict the interactions of components in the system, how they respond to one another, and how they perform as a response to changes in operating conditions. A suitable and enlightening way to demonstrate and foresee these interactions is with the use of accurate modeling software. Sage, ANSYS® Maxwell, and ANSYS® Fluent are the current modeling tools used by NASA to analytically determine convertor performance. Sage is a one-dimensional object-oriented commercial software package used for modeling and optimizing Stirling convertors for Dynamic Radioisotope Power Systems (DRPS) and it is one of the most accurate Stirling convertor codes in use by NASA. This code is the successor to GLIMPS (Globally-Implicit Stirling Cycle Simulation) and GLOP (GLIMPS Optimization) software created by Gedeon Associates [1]. Model input parameters are typically material/gas type, component physical dimensions, temperatures, frequency, charge pressure, and number of time/space nodes. Sage is used to model both the FISC’s and SRSC’s Stirling cycle thermodynamics and piston/displacer dynamics. Performance maps were created and analyzed for both power systems to better understand the relationship between the following conditions: cold-end temperature, hot-end temperature, piston/displacer amplitudes, pressure drop, and thermal input power. The synergy between these conditions will help determine parameter sensitivity. ANSYS® Maxwell was used to create a three-dimensional (3-D) axisymmetric model for both FISC and SRSC alternators. The significant physical components included in each model are the magnets, magnet carrier, outer/inner laminations, and the coil. Inputs to the model are piston amplitude, piston frequency, alternator load, coil resistance, tuning capacitance, and specific material properties. The alternator models calculate terminal voltage, current, piston/current phase, voltage/current phase, coil inductance, terminal power and efficiency. The RI2 losses, core (hysteresis and eddy) losses, and magnet/can eddy losses are also a part of the final results. ANSYS® Fluent is used to build 3-D computational fluid dynamic (CFD) models to examine the Stirling cycle thermodynamics for both the FISC and SRSC systems. Three-dimensional Computer Aided Design (CAD) models were used to create the physical components of each convertor. Steady-state simulations were conducted for hardware testing, prediction of environmental losses during testing, and generation of radiation look-up tables. The model inputs to the aforementioned analysis are the material properties and boundary thermal conditions. The steady-state model calculates temperature and heat flow distributions. Transient 3-D calculations were also part of the CFD analysis. In this study a physically reduced version of the FISC is used to obtain a prediction of available engine power. For the gas bearing SRSC, the transient effort is used to obtain a prediction of bearing pad performance and its sensitivity to micro-channel geometric variation. The model inputs to the transient simulations are the piston amplitude, displacer amplitude, frequency, displacer/piston phase angle, dynamic deforming CFD grid, temperature boundary conditions, and user defined files describing motion profile of piston/displacer. The results of the model are temperature distributions, heat distributions, and PV power produced at pre-determined conditions.

Luis A Rodriguez↗

Airfoil Table Generation and Validation for the VR-12 and SSC-A09 Airfoils and Quadrotor Performance Prediction

As part of NASA’s Urban Air Mobility (UAM) mission, the Revolutionary Vertical Lift Technology (RVLT) project is creating rotorcraft reference designs that can be used by the rotorcraft community to develop quiet, efficient, and safer air vehicles. At NASA, this effort is made possible by the use of the RVLT toolchain software. Airfoil tables are often a critical part of the conceptual design process as they inform both comprehensive analysis and Computational Fluid Dynamics (CFD) codes of 2D aerodynamic coefficients, such as lift, drag, and moment coefficients. These airfoil tables, typically generated during experimental testing, are often proprietary and thus not widely distributable; further, there is not always clear or available reference documentation that provides information on the test conditions. The airfoil tables from these experiments are referred to in this work as Legacy airfoil tables. The scope of this study is to generate airfoil tables with the NASA OVERset grid CFD FLOW solver (OVERFLOW) coupled with the AirFoil Table Generator (AFTGen) software and validate these airfoil tables based on performance predictions calculated in Comprehensive Hierarchical Aeromechanics Rotorcraft Model (CHARM). The overall goal is to provide CFD generated airfoil tables, validate them through comprehensive analysis, and publish the airfoil tables for public distribution. The RVLT quadrotor reference model was selected as the basis for all analyses. The quadrotor uses the VR-12 and SSC-A09 airfoils for its rotor blades. Simulations were run in OVERFLOW for the VR-12 and SSC-A09 airfoils for ±20° angles of attack, Mach numbers from 0.3 to 1.0, and a Reynolds-Mach proportionality constant of 1.023x107. The C81 format airfoil tables from AFTGen were then blended with Legacy NACA 0012 airfoil table data for ±20° to ±180° angle of attack range. CHARM was used to compute thrust coefficient, power coefficient, and figure of merit for an isolated rotor in hover using proprietary legacy airfoil tables for the VR-12 and SSC-A09, as well as with the OVERFLOW generated airfoil tables for the VR-12 and SSC-A09. Overall, the OVERFLOW and Legacy performance predictions agreed well, with the best correlation between simulated and Legacy data observed with the VR-12 airfoil tables.

Airfoil↗

Retrieval of energy spectra for all flavours of neutrinos from core-collapse supernova with multiple detectors

ABSTRACT We present a new method by which to retrieve energy spectrum for all falvours of neutrinos from core-collapse supernova (CCSN). In the retrieval process, we do not assume any analytic formulas to express the energy spectrum of neutrinos but rather take a direct way of spectrum reconstruction from the observed data; the singular value decomposition algorithm with a newly developed adaptive energy-gridding technique is adopted. We employ three independent reaction channels having different flavour sensitivity to neutrinos. Two reaction channels, inverse beta decay on proton and elastic scattering on electrons, from a water Cherenkov detector such as Super-Kamiokande (SK) and Hyper-Kamiokande (HK), and a charged current reaction channel with Argon from the Deep Underground Neutrino Experiment (DUNE) are adopted. Given neutrino oscillation models, we iteratively search the neutrino energy spectra at the CCSN source until they provide the consistent event counts in the three reaction channels. We test the capability of our method by demonstrating the spectrum retrieval to a theoretical neutrino data computed by our recent three-dimensional CCSN simulation. Although the energy spectrum with either electron-type or electron-type antineutrinos at the CCSN source has relatively large error compared to that of other species, the joint analysis with HK + DUNE or SK + DUNE will provide precise energy spectrum of all flavours of neutrinos at the source. Finally, we discuss perspectives for improvements of our method by using neutrino data of other detectors.

79 ASTRONOMY AND ASTROPHYSICS↗

Real-time Simulation of Coordinated Sub-Transmission Volt-Var Control Tool under High Distributed PV Penetration

With the rapid penetration of intermittent solar photovoltaic (PV) and other distributed energy resources (DER) into the grid, and subsequent changes in power flow patterns in both distribution and sub-transmission, voltage regulation is becoming a major challenge. It is prudent to leverage PV and DER to provide ancillary services to the grid, such as voltage regulation. Recently, a quasi-static, coordinated real-time sub transmission volt-var control algorithm (CReST-VCT) was developed for voltage regulation under high PV penetration by dispatching the reactive power settings of the shunt devices and PV inverters. The algorithm was validated offline, on a quasistatic study. Such quasi-static dispatching algorithms cannot guarantee its performance for highly nonlinear dynamical power systems. In this paper the performance of the algorithm is validated with a real time nonlinear dynamic simulation of a modified IEEE 118 bus system in Opal-RT. The real time simulation is used to emulate the actual system operation, providing a more realistic testing environment for CReST-VCT. The dispatch and control signals are communicated between the power system (Opal solver) and the control center (GAMS solver) in real-time with a MODBUS bridge. The results demonstrate i) system is stable with the new dispatch points ii) significant improvement in system-wide voltage profiles compared to an uncontrolled scenario. Another significant contribution of the work is developing a framework for dispatchable, long duration dynamic simulations that can be leveraged for market/dispatch studies.

CReST-VCT, Opal-RT, Volt/Var Control, renewable in↗

Impacts of stochastic coalescence variability on warm rain initiation using Lagrangian microphysics in box and large-eddy simulations

Various coalescence methods for Lagrangian microphysics schemes are tested in box and large-eddy simulation (LES) models, including the stochastic all-or-nothing super-droplet method (SDM) and a version of SDM (dSDM) that applies a fractional approach similar to the average impact method. In LES, variabilities driven by microphysics and by flow realizations are separated using the “piggybacking” technique. Rain initiation averaged over many realizations of the box model is delayed and rain variability increases as the number of super-drops per collision volume ( N SD ) is decreased using SDM. In contrast, rain initiation time using SDM in LES is insensitive to N SD for 32 ≤ N SD ≤ 512. This is explained through the interaction between LES grid boxes, each acting as a separate collision volume. Variability across the ensemble of LES collision volumes using SDM results in rain quickly initiating in some of the LES grid cells at low N SD and leading to a similar overall timing of rain initiation from the cloud compared to simulations with high N SD . There is a ~20% decrease in the total rain mass and mean rain flux as N SD is increased from 32 to 256, with little additional change as N SD is increased from 256 to 512. The fractional coalescence approach in dSDM leads to reduced microphysical variability and a 15-18 min delay in rain initiation compared to SDM. Finally, an additional LES ensemble with microphysical variability feeding back to the dynamics shows that flow variability dominates the impact of microphysical variability on rain properties. Thus, flow variability must be constrained to isolate impacts of microphysical variability.

54 ENVIRONMENTAL SCIENCES↗

Harmonized Automatic Relay Mitigation of Nefarious Intentional Events (HARMONIE) - Special Protection Scheme (SPS)

The harmonized automatic relay mitigation of nefarious intentional events (HARMONIE) special protection scheme (SPS) was developed to provide adaptive, cyber-physical response to unpredictable disturbances in the electric grid. The HARMONIE-SPS methodology includes a machine learning classification framework that analyzes real time cyber-physical data and determines if the system is in normal conditions, cyber disturbance, physical disturbance, or cyber-physical disturbance. This classification then informs response, if needed and/or suitable, and included cyber-physical corrective actions. Beyond standard power system mitigations, a few novel approaches were developed that included a consensus algorithm-based relay voting scheme, an automated power system triggering condition and corrective action pairing algorithm, and a cyber traffic routing optimization algorithm. Both the classification and response techniques were tested within a newly integrated emulation environment composed of a real-time digital simulator (RTDS) and SCEPTRE™. This report details the HARMONIE-SPS methodology, highlighting both the classification and response techniques, and the subsequent testing results from the emulation environment.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Diabatic forcing and intialization with assimilation of cloud water and rainwater in a forecast model

In this study, diabatic forcing, and liquid water assimilation techniques are tested in a semi-implicit hydrostatic regional forecast model containing explicit representations of grid-scale cloud water and rainwater. Diabatic forcing, in conjunction with diabatic contributions in the initialization, is found to help the forecast retain the diabatic signal found in the liquid water or heating rate data, consequently reducing the spinup time associated with grid-scale precipitation processes. Both observational Special Sensor Microwave/Imager (SSM/I) and model-generated data are used. A physical retrieval method incorporating SSM/I radiance data is utilized to estimate the 3D distribution of precipitating storms. In the retrieval method the relationship between precipitation distributions and upwelling microwave radiances is parameterized, based upon cloud ensemble-radiative model simulations. Regression formulae relating vertically integrated liquid and ice-phase precipitation amounts to latent heating rates are also derived from the cloud ensemble simulations. Thus, retrieved SSM/I precipitation structures can be used in conjunction with the regression-formulas to infer the 3D distribution of latent heating rates. These heating rates are used directly in the forecast model to help initiate Tropical Storm Emily (21 September 1987). The 14-h forecast of Emily's development yields atmospheric precipitation water contents that compare favorably with coincident SSM/I estimates.

Raymond, William H.↗

TREAT Neutronics Analysis of Water-Loop Concept Accommodating LWR 9-rod Bundle

Abstract. Simulation of a variety of transient conditions has been successfully achieved in the Transient Reactor Test (TREAT) facility during operation between 1959 and 1994 to support characterization and safety analysis of nuclear fuels and materials. A majority of previously conducted tests were focused on supporting sodium-cooled fast reactor (SFR) designs. Experiments evolved in complexity. Simulation of thermal-hydraulic conditions expected to be encountered by fuels and materials in a reactor environment was realized in the development of TREAT sodium loop experiment vehicles. These loops accommodated up to 7-pin fuel bundles and served to simulate more closely the reactor environment while safely delivering large quantities of energy into the test specimen. Some of the immediate TREAT restart operations will be focused on testing light water reactor (LWR) accident tolerant fuels (ATF). Similar to the sodium loop objectives, a water loop concept, developed and analyzed in the 1990’s, aimed at achieving thermal-hydraulic conditions encountered in commercial power reactors. The historic water loop concept has been analyzed in the context of a reactivity insertion accident (RIA) simulation for high burnup LWR 2-pin and 3-pin fuel bundles. Findings showed sufficient energy could be deposited into the specimens for evaluation. Similar results of experimental feasibility for the water loop concept (past and present) have recently been obtained using MCNP6.1 with ENDF/B-VII.1 nuclear data libraries. The old water loop concept required only two central TREAT core grid spaces. Preparation for future experiments has resulted in a modified water loop conceptual design designated the TREAT water environment recirculating loop (TWERL). The current TWERL design requires nine TREAT core grid spaces in order to place the water recirculating pump under the TREAT core. Due to the effectiveness of water moderation, neutronics analysis shows that removal of seven additional TREAT fuel elements to facilitate the experiment will not inhibit the ability to successfully simulate a RIA for the 2-pin or 3-pin bundle. This new water loop design leaves room for accommodating a larger fuel pin bundle than previously analyzed. The 7-pin fuel bundle in a hexagonal array with similar spacing of fuel pins in a SFR fuel assembly was considered the minimum needed for one central fuel pin to encounter the most correct thermal conditions. The 9-rod fuel bundle in a square array similar in spacing to pins in a LWR fuel assembly would be considered the LWR equivalent. MCNP analysis conducted on a preliminary LWR 9-rod bundle design shows that sufficient energy deposition into the central pin can be achieved well within range to investigate fuel and cladding performance in a simulated RIA. This is achieved by surrounding the flow channel with an additional annulus of water. Findings also show that a highly significant increase in TREAT to specimen power coupling factor (PCF) within the central pin can be achieved by surrounding the experiment with one to two rings of TREAT upgrade fuel assemblies. The experiment design holds promise for the performance evaluation of PWR fuel at extremely high burnup under similar reactor environment conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Loosely Conditioned Emulation of Global Climate Models With Generative Adversarial Networks

Climate models encapsulate our best understanding of the Earth system, allowing research to be conducted on its future under alternative assumptions of how human-driven climate forces are going to evolve. An important application of climate models is to provide metrics of mean and extreme climate changes, particularly under these alternative future scenarios, as these quantities drive the impacts of climate on society and natural systems. Because of the need to explore a wide range of alternative scenarios and other sources of uncertainties in a computationally efficient manner, climate models can only take us so far, as they require significant computational resources, especially when attempting to characterize extreme events, which are rare and thus demand long and numerous simulations in order to accurately represent their changing statistics. Here we use deep learning in a proof of concept that lays the foundation for emulating global climate model output for different scenarios. We train two "loosely conditioned" Generative Adversarial Networks (GANs) that emulate daily precipitation output from a fully coupled Earth system model: one GAN modeling Fall-Winter behavior and the other Spring-Summer. Our GANs are trained to produce spatiotemporal samples: 32 days of precipitation over a 64x128 regular grid discretizing the globe. We evaluate the generator with a set of related performance metrics based upon KL divergence, and find the generated samples to be nearly as well matched to the test data as the validation data is to test. We also find the generated samples to accurately estimate the mean number of dry days and mean longest dry spell in the 32 day samples. Our trained GANs can rapidly generate numerous realizations at a vastly reduced computational expense, compared to large ensembles of climate models, which greatly aids in estimating the statistics of extreme events.

climate emulation, extreme climate, impacts, machi↗

Low Cost High Efficiency Magnetrons for Accelerator Applications

High efficiency, low cost RF sources are required for proton, ion, and muon accelerators. These sources must be compact and provide precise control of the output power and phase. Currently only large, expensive klystrons can provide this performance. Calabazas Creek Research Inc. proposes to develop a phase locked magnetron-based system that has phase and amplitude control via a low power amplifier. The system uses phase modulation to control the RF power. The efficiency can exceed 80% and the cost is expected to be less than $3/Watt, significantly lower than other high power amplifiers. CCR originally proposed to develop a magnetron with a grid for control of phase and amplitude. Detailed simulations showed that this was not a promising approach. The grid was found to be viable only in magnetrons with low duty cycle, and amplitude control was demonstrated; but the grid was found to be not effective for phase locking. Attention was turned to a concept from Fermilab, which is very promising for both amplitude and phase control and does not require a grid. CCR will develop an RF system providing phase and amplitude control of a 100 kW peak, 10 kW average 1.3 GHz magnetron. Procedures and techniques will be applicable to design of systems at other frequencies and power levels. The program will include development of the magnetron and locking electronics. The system will be tested at Fermilab. Commercial Applications and Other Benefits Magnetrons can provide high levels of RF power at very high efficiency at low cost. As oscillators, however, their applications are limited. Magnetrons with more precise control of the output power and control of the phase could find wide application in accelerators and colliders. They would also provide improved performance for high resolution radar.

43 PARTICLE ACCELERATORS↗