Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation output analysis”

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 991 records · Page 55

Holistic Numerical Calibration of the iMESA Electrostatic Analyzer

A drifted Maxwellian velocity distribution is the most common model used to interpret the data from low-energy charged-particle instruments onboard spacecraft that are used to investigate the ambient plasma environment in the low Earth orbit (LEO). An original method is presented for determining the flow parameters (density, temperature, and flow energy) of such a distribution from the output of the integrated miniaturized electrostatic analyzer, which has been successfully flown on several LEO missions. Rather than attempting to deconvolve from the on-orbit data the analyzer’s response to an ideal, monoenergetic input, numerical simulation is used to predict and parameterize the response of the device to an input distribution that includes an isotropic, non-zero temperature, yielding a straightforward method for extracting the flow parameters from the spacecraft data. The method is computationally simple enough to be incorporated into a robust algorithm suitable for rapid batch processing or real-time analysis of data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Datasets for Widespread Residential Space Heating Electrification in Texas

In this experiment, we explore long term patterns in electricity demand driven by the dual effects of full electrification of space heating in Texas (by adoption of electric heat pumps), and climate change. We use a predictive model of electricity demand, climate projections, and an open source nodal power system (DC Optimal Power Flow) model of the Electric Reliability Council of Texas (ERCOT) system. Heat pumps are a more energy efficient way of providing space heating and cooling in homes. We attempt to exhaustively investigate the impacts of full residential space heating electrification by adoption of heat pumps for the segment of Texas households that currently rely on fossil fuels (about 40%), while simultaneously incorporating climate change meteorological variables. We explore a range of scenarios of heat pump efficiency and climate uncertainty over a long period of future years (2020-2099). In total, the simulation experiment generates 1,280 simulation years of hourly data. We report and analyze results in form of impacts on residential load, total load, peak load, seasonality of peaking, and reliability measured by occurrence and frequency of loss of load events. While the experiment is for ERCOT, the insights and approach can be applied to other regions. The results from the analysis can inform system planners on a range of potential capacity requirements/ reliability implications and/or risks of full space heating electrification via the adoption of electric heat pumps, given the uncertainty in the scenarios/ climate futures. The dataset includes model output for residential, non residential and total load, and the results from the GO ERCOT model runs for 4 RCP Scenarios (RCP 4.5 Cooler, RCP 4.5 Hotter, RCP 8.5 Cooler, RCP 8.5 Hotter), 4 Heating electrification Scenarios (Base , Standard Efficiency HP, High Efficiency HP, Ultra-High Efficiency HP) over 80 years (2020-2099). The metrological variables at BA scale were weighted weighted using population projections consistent with the SSP3 scenario.

Climate Change↗

Drive turbine systems for 20-inch turbofan simulators. 1: Duct turbine design

A study was made to evolve the turbine drive systems for 20-inch turbofan engine simulators. The fan designs used in the simulators included single-stage and two-stage configurations that covered a wide range of rotative speed and power requirement. The objective assumed for the study was to evolve one core turbine design that could drive all of the single-stage fans and, when operated in combination with one duct turbine design, drive all of the two-stage fans. The duct turbine power output is then needed to determine the make-up power required of the core turbine over the range of two-stage fan operating conditions. The duct turbine design analysis is reported and includes the selection of the duct turbine velocity diagram, a description of the blade design, and a determination of its off-design performance. Adjustable stators were found to be quite advantageous to the duct turbine off-design operation. The use of adjustable stators enabled the duct turbine to accommodate fan mass flow at all operating points and caused the duct turbine power output to increase as the total power requirement increased. This in turn resulted in a core turbine make-up power requirement that was not significantly greater than that required for driving the single-stage fans.

Whitney, W. J.↗

Control theory analysis of a three-axis VTOL flight director

A control theory analysis of a VTOL flight director and the results of a fixed-based simulator evaluation of the flight-director commands are discussed. The VTOL configuration selected for this study is a helicopter-type VTOL which controls the direction of the thrust vector by means of vehicle-attitude changes and, furthermore, employs high-gain attitude stabilization. This configuration is the same as one which was simulated in actual instrument flight tests with a variable stability helicopter. Stability analyses are made for each of the flight-director commands, assuming a single input-output, multi-loop system model for each control axis. The analyses proceed from the inner-loops to the outer-loops, using an analytical pilot model selected on the basis of the innermost-loop dynamics. The time response of the analytical model of the system is primarily used to adjust system gains, while root locus plots are used to identify dominant modes and mode interactions.

Niessen, F. R.↗

ASTEP user's guide and software documentation

The Algorithm Simulation Test and Evaluation Program (ASTEP) is a modular computer program developed for the purpose of testing and evaluating methods of processing remotely sensed multispectral scanner earth resources data. ASTEP is written in FORTRAND V on the UNIVAC 1110 under the EXEC 8 operating system and may be operated in either a batch or interactive mode. The program currently contains over one hundred subroutines consisting of data classification and display algorithms, statistical analysis algorithms, utility support routines, and feature selection capability. The current program can accept data in LARSC1, LARSC2, ERTS, and Universal formats, and can output processed image or data tapes in Universal format.

Gliniewicz, A. S.↗

Probabilistic Look-ahead Contingency Analysis Integration with Commercial Tool and Practical Data

This paper presents an initial effort of integrating a smart sampling-based probabilistic look-ahead contingency analysis algorithm with General Electric (GE) Grid Solutions’ commercial energy management system (EMS) tool as a proof-of-concept for a seamless research tool integration using real world large-scale grid data. With the increasing impact of random forces such as variable generation and load, their stochastic behaviors cannot be ignored. However, the current practices are still dominated by deterministic tools. They are becoming increasingly inadequate for the future grid. The developed look-ahead contingency analysis algorithm incorporates forecast errors of variable energy and load to address the challenges brought by the increasing uncertainty of power system. The algorithm can reveal the potential violations caused by the variance of variable energy and load that are not normally detected by traditional deterministic approaches. To test its performance under practical environments ( real data with real commercial tool), significant efforts have been made to prepare test cases, modify GE EMS tool, and adapt an extreme value distribution algorithm to analyze the GE EMS’s violation-only outputs. The test results clearly demonstrate the effectiveness of the developed algorithm as new transformer violations that were not previously detected have been identified. This performance provides better situational awareness to engineers for their decision-making process under uncertainty. Moreover, with the discussion of computational performance and future work, this paper has shown a clear path for integrating the probabilistic algorithm with commercial tools to make us better equipped for the changing power system.

Modeling and simulation of power systems, constrai↗

Radiation Transport in Random Media With Large Fluctuations

Strongly heterogeneous media arise in several applications that include radiation shields, nuclear fuel, BWR moderators, clouds, planetary and stellar atmospheres, turbulent gases and plasmas. Neutral and charged particle transport computations in such media have relied heavily on formulating transport equations with spatially random coefficients (physical data) and developing solution methods to deal with the additional stochastic dimensions. Attempts at developing approximate closures that yield only low order statistical information (e.g., mean and variance of the flux) have proved to be highly restrictive under real physics conditions or rely on techniques that require fluctuation amplitudes to be small for robustness. Recently, stochastic spectral methods such as polynomial chaos and stochastic collocation have been developed for aleatoric uncertainty quantification and sensitivity analysis, and successfully applied in radiation transport work. Advances in UQ techniques have tended to focus on efficiently handling large numbers of uncertain variables but the rigorous stochastic basis of the approach also promotes its use in situations where stochasticity is due to spatial heterogeneity and the associated uncertainty is large. Here we apply these techniques to radiation transport in media with spatially randomly varying cross sections without restriction on fluctuation amplitudes. Specifically, we represent the cross sections as a lognormal spatial random process with specified mean, variance and covariance function and use a Karhunen-Loeve (KL) decomposition to generate cross section realizations that are strictly positive. Woodcock Monte Carlo (WMC) is then used to simulate transport using random sampling of cross sections and deterministic sampling based on a stochastic collocation technique. Numerical results for the mean and variance in the scalar flux and leakage currents are obtained for weak and strong random variations. We focus in this work on demonstrating the approach and defer application to specific problems to a future investigation. We have demonstrated that a combination of stochastic spectral representation and Woodcock Monte Carlo simulation allows efficient computation of radiation transport in strongly random media. The use of a lognormal distribution to represent the fluctuations in medium properties facilitates enforcement of strict positivity of cross section realizations, a condition that has challenged standard approaches that assume Gaussian fluctuations. The fact that the lognormal random process is a memoryless nonlinear transformation of a Gaussian process further enables efficient and accurate reconstruction of the random cross section realizations using a Karhunen-Loeve representation of the cross section. Finally, the robustness of this model enables fluctuations of arbitrarily large amplitude to be studied, using both random sampling Monte Carlo and more efficient stochastic collocation techniques. The present work is being extended to allow reconstruction of probability density functions of output variables using polynomial chaos representation, and to accommodate general cross section covariance functions. Also, techniques are being explored to more judiciously select ceiling cross sections to enhance the efficiency of Woodcock Monte Carlo for transport in random media. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Representing Socio‐Economic Uncertainty in Human System Models

Abstract Socio‐economic development pathways and their implications for the environment are highly uncertain, and energy transitions will involve complex interactions among sectors. Here, traditional Monte Carlo analysis is paired with scenario discovery techniques to provide a richer portrait of these complexities. Modeled uncertain input variables include costs of advanced energy technologies, energy efficiency trends, fossil fuel resource availability, elasticities of substitution for labor, capital, and energy across economic sectors, population growth, and labor and capital productivity. The sampled values are simulated through a multi‐sector, multi‐region, recursively dynamic model of the world economy to explore a range of possible future outcomes. We find that many patterns of energy and technology development are possible for various long‐term environmental pathways and that sectoral output for most sectors is little affected through 2050 by the long‐term temperature target, but with tight constraints on emissions, emission intensities must fall much more rapidly. Scenario discovery techniques are applied to the large uncertainty ensembles to explore if there are prevailing storylines behind outcomes of interest. An illustrative investigation focused on different levels of economic growth shows many combinations of pathways and no single storyline emerging for a given economic outcome. This method can be extended to other outcomes of interest, exploring the nature of scenarios with both tail and median outcomes. Sampling from a Monte Carlo generated ensemble provides a rich set of scenarios to investigate, and potentially aids in avoiding heuristic biases in less structured scenario approaches.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

System Identification of a DC-DC Buck Converter Based on Two-Channel Relay Method

This work presents a system identification approach based on a two-channel relay applied to a DC-DC buck converter operating in closed-loop with a PI controller to regulate its output voltage. The implemented two-relay method allows system identification while the system is running online. Moreover, the algorithm works in a closed-loop configuration while introducing a minimal perturbation in the converter's operation. The presented algorithm was used to identify the frequency response of the converter with no-prior knowledge about the system structure. Furthermore, the identified transfer function parameters has a good agreement with the mathematical model of the converter. Simulation results of the comparison between the switching model, mathematical model and the identified model are provided to validate the theoretical analysis.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Data-Enabled Predictive Control for Building HVAC Systems

Model predictive control is widely used as a control technology for the computation of optimal control inputs of building heating, ventilating, and air conditioning (HVAC) systems. However, both the benefits and widespread adoption of model predictive control (MPC) are hindered by the effort of model creation, calibration, and accuracy of the predictions. In this paper, we apply the data-enabled predictive control (DeePC) algorithm for designing controls for building HVAC systems. The algorithm solely depends on input/output data from the system to predict future state trajectories without the need for system identification. The algorithm relies on the idea that a vector space of all input–output trajectories of a discrete-time linear time-invariant (LTI) system is spanned by time-shifts of a single measured trajectory, given the input signal is persistently exciting. Closed-loop simulations using EnergyPlus are performed to demonstrate the approach. The simulated building modeled in EnergyPlus is a modified commercial large office prototype building served by an air handling unit-variable air volume HVAC system. Temperature setpoints of zones are used as control variables to minimize the HVAC energy cost of the building considering a time-of-use electricity rate structure. Furthermore, sensitivity analysis is conducted to gain insights into the effect of parameter tuning on DeePC performance. Simulation results are used to illustrate the performance of the algorithm and compare the algorithm with model-based MPC and occupancy-based setpoint controller. Overall, DeePC achieves similar performance compared to MPC for lower engineering effort.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Forcing data (CESM2/CMIP6) for projection of drought impacts (2015-2100) at the K34 site in Manaus, Brazil

Historical and projected output data variables extracted and derived from Community Earth System Model 2 (CESM2) runs from the Coupled Model Intercomparison Project Phase 6 (CMIP6) archive. CESM2 is a fully coupled Earth system model used in simulations of Earth's past, present, and future climates (Danabasoglu et al., 2020). Variables in this dataset are in six-hourly resolution, and include air temperature (both in K and ℃), specific humidity (kg kg-1), air pressure (Pa), relative humidity (%), and vapor pressure deficit (kPa). CESM2 was the only CMIP6 model that provided VPD at the high temporal resolution required for this analysis. Data are included in .csv files, and the text file CESM2-CMIP6_forcing_K34-Manaus_headers.txt provides descriptions of data file headers.

54 ENVIRONMENTAL SCIENCES↗

Ex-Core Radiation Transport Modeling with VERA User Manual

It is important to be able to accurately predict the neutron flux outside the immediate reactor core for a variety of safety and material analyses. Monte Carlo radiation transport calculations are required to produce these high–fidelity ex-core responses. The Virtual Environment for Reactor Applications (VERA) provides the automated capability to launch independent Shift fixed-source and eigenvalue Monte Carlo (MC) calculations for user-specified state points during a standard Virtual Environment for Reactor Applications (VERA) calculation. VERA couples MPACT with COBRA-TF (CTF) to Shift to perform ex-core tallies for multiple state points concurrently, with each component capable of parallel execution on independent processor domains. In these ex-core calculations, MPACT is coupled to CTF and performs the in-core depletion and heat transfer calculation, followed by a fixed-source Shift transport calculation including ex-core regions to produce ex-core responses. The fission source, fuel pin temperatures, moderator temperature and density, boron concentration, and fuel pin depleted isotopic compositions can be transferred to Shift from the MPACT calculation. Specifically, VERA performs fluence calculations in the core barrel outward to the end of the pressure vessel and detector response calculations in ex-core detectors. It also performs the requested tallies in any user-defined ex-core regions. VERA takes advantage of the General Geometry (GG) package in Shift. This gives VERA the flexibility to explicitly model features outside the core barrel, including detailed vessel models, detectors, and power plant details. A very limited set of experimental and numerical benchmarks is available for ex-core simulation comparison. The Consortium for Advanced Simulation of Light Water Reactors has developed a set of ex-core benchmark problems to include as part of the VERA verification and validation set of problems. The ex-core capability in VERA has been tested on small representative assembly problems, multi-assembly problems, as well as quarter-core and full-core problems. VERA View has also been extended to visualize these vessel fluence results from VERA. This manual serves to present a guide to VERA users about the methodology behind ex-core calculations and the details of input, output, and analysis of results from these calculations. Details in this version of the manual are based on features in VERA 4.0.1.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Orbit determination error analysis and comparison of station-keeping costs for Lissajous and halo-type libration point orbits and sensitivity analysis using experimental design techniques

Spacecraft in orbit near libration point L1 in the Sun-Earth system are excellent platforms for research concerning solar effects on the terrestrial environment. One spacecraft mission launched in 1978 used an L1 orbit for nearly 4 years, and future L1 orbital missions are also being planned. Orbit determination and station-keeping are, however, required for these orbits. In particular, orbit determination error analysis may be used to compute the state uncertainty after a predetermined tracking period; the predicted state uncertainty levels then will impact the control costs computed in station-keeping simulations. Error sources, such as solar radiation pressure and planetary mass uncertainties, are also incorporated. For future missions, there may be some flexibility in the type and size of the spacecraft's nominal trajectory, but different orbits may produce varying error analysis and station-keeping results. The nominal path, for instance, can be (nearly) periodic or distinctly quasi-periodic. A periodic 'halo' orbit may be constructed to be significantly larger than a quasi-periodic 'Lissajous' path; both may meet mission requirements, but perhaps the required control costs for these orbits are probably different. Also for this spacecraft tracking and control simulation problem, experimental design methods can be used to determine the most significant uncertainties. That is, these methods can determine the error sources in the tracking and control problem that most impact the control cost (output); it also produces an equation that gives the approximate functional relationship between the error inputs and the output.

Gordon, Steven C.↗

Fault Diagnosis in HVAC Chillers

Modern buildings are being equipped with increasingly sophisticated power and control systems with substantial capabilities for monitoring and controlling the amenities. Operational problems associated with heating, ventilation, and air-conditioning (HVAC) systems plague many commercial buildings, often the result of degraded equipment, failed sensors, improper installation, poor maintenance, and improperly implemented controls. Most existing HVAC fault-diagnostic schemes are based on analytical models and knowledge bases. These schemes are adequate for generic systems. However, real-world systems significantly differ from the generic ones and necessitate modifications of the models and/or customization of the standard knowledge bases, which can be labor intensive. Data-driven techniques for fault detection and isolation (FDI) have a close relationship with pattern recognition, wherein one seeks to categorize the input-output data into normal or faulty classes. Owing to the simplicity and adaptability, customization of a data-driven FDI approach does not require in-depth knowledge of the HVAC system. It enables the building system operators to improve energy efficiency and maintain the desired comfort level at a reduced cost. In this article, we consider a data-driven approach for FDI of chillers in HVAC systems. To diagnose the faults of interest in the chiller, we employ multiway dynamic principal component analysis (MPCA), multiway partial least squares (MPLS), and support vector machines (SVMs). The simulation of a chiller under various fault conditions is conducted using a standard chiller simulator from the American Society of Heating, Refrigerating, and Air-conditioning Engineers (ASHRAE). We validated our FDI scheme using experimental data obtained from different types of chiller faults.

Choi, Kihoon↗

End-to-end Modeling of Lyot Coronagraph for LUVOIR: Sensitivity to Segment-Level Aberrations

Direct imaging and spectroscopic characterization of atmospheres of Earth-like exoplanets require starlight suppression to detect a faint planetary signal at small angular separations. Future space telescopes will require segmentation to survey a statistically significant population of nearby stars systems by increasing aperture size beyond the limits set by traditional monolithic apertures; however, segmentation also represents an additional challenge due to the introduction of temporal pupil discontinuities in terms of both phase and amplitude. For this study, a baseline Lyot coronagraph is used with the Large UV / Optical /Infrared (LUVOIR) mission’s 15-m on-axis, segmented optical telescope assembly (OTA). We present a high-fidelity, end-to-end modeling approach developed as part of the Segmented Coronagraph Design and Analysis (SCDA) group that includes diffractive optical propagations at the 1e-10 contrast level across the entire optical train including the OTA and intermediate re-imaging optics inside the coronagraph instrument. We perform a tolerancing analysis to assess the contrast sensitivity to wavefront perturbations for a set of key telescope terms including segment-level low-order and mid-spatial frequency aberrations. As an input to the end-to-end telescope model, we use two industry-provided segment-level wavefront aberration datasets pertaining: (1) thermal gradients and accelerations (courtesy of Ball-Aerospace and L3-Harris), and (2) telescope dynamics including line-of-sight and segment jitter (courtesy of Lockheed-Martin). Contrast performance outputs of the coronagraph are fed into a mission design reference simulator to assess the scientific yield sensitivity. End-to-end modeling will incorporate performance improvements due to active wavefront sensing and control and/or segment-level metrology to improve system performance.

Optics↗

Analysis of Operating Strategies Using Different Target Designs For 238Pu Production

An engineering effort is underway to re-establish capability to produce 238Pu oxide at the kilogram scale in the United States. A multi-step batch process is being developed to produce this important material. Recently, a portion of this process was studied using discrete-event simulation tools to determine whether the conceptual process might achieve its yearly production goal. The study showed the conceptual process can meet the yearly production goal under some circumstances, but process improvements would be needed to ensure greater likelihood of success. This study extends the work performed previously by examining the effects of changing the reactor target design on the yearly process output. Two new reactor target configurations are considered an aluminum-clad reactor target containing 50% greater 237Np oxide content than the original target, and a zirconium alloy-clad target using no aluminum. The results indicate that use of the new aluminum-clad target configuration may allow the process to achieve the same yearly production goal in less time using fewer targets. If the zirconium alloy-clad target is used, then even fewer targets would be needed to reach the production goal, but some process changes would be required to handle the zirconium cladding. The number of days needed to process a target batch to completion, and the steady state 238Pu oxide production rate, for each configuration are compared to the results from the initial simulation study.

campaign↗

Light Output and Thermal Blast Analysis of Nuclear Fireballs

During the 1940s, 50s, and 60s, there were 210 atmospheric nuclear tests. Every test was extensively recorded with high quality motion picture films and photographs, providing a unique historic record, containing ~10,000 films, with a total of ~22 million frames of information. Much of what is known about the effects of nuclear weapons has been obtained from these tests. There are approximately 6500 scientific films that are in the process of being scanned by Lawrence Livermore National Laboratory and re-analyzed using modern computer techniques. The four major effects of a nuclear detonation are prompt radiation, shockwave, thermal blast, and fallout. The objective of this interim report is to describe recent efforts to extract benchmark data from the existing streak films using modern Figure-processing techniques. These new and, hopefully, greatly-improved data will provide our code developers with the benchmark data necessary to validate our modern computer codes that simulate nuclear weapon effects.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Initial Performance Evaluation and Ranking of Thermal Energy Storage Options for Light Water Reactor Integration to Support Modeling and Simulation

This report provides an in-depth analysis of current thermal storage technologies in the marketplace as of 2019 and develops a phenomenological identification ranking table (PIRT) and Figure of Merit (FOM) study for near-term thermal energy storage technologies with light water reactor technology. Particular focus was given the NuScale reactor as part of the Joint Use Modular Plant (JUMP) program between INL, NuScale, and the Utah Associated Municipal Project (UAMPs). Ancillary service industries and technologies are explored in detail to ascertain thermodynamic requirements for integration with thermal energy storage technologies. These industries include steel manufacturing, hydrogen production, desalination, pulp and paper manufacturing. It was determined that a large percentage of these industries can make use of existing nuclear output temperatures and pressures. Storage technologies that can maintain these values are advantageous for coupling purposes. The report then analyzes ten thermal energy storage technologies, with thirteen specific systems discussed in detail regarding their potential to integrate with LWR technology. A ranking tool identified important characteristics of thermal storage and ranked each of the technologies. Concrete, molten salt, and thermal oil sensible heat storages, and steam accumulators evaluated highly. After including cost estimations and qualitatively comparing these technologies, it is recommended that either molten salt or thermal oil is used for implementation with the JUMP module. The TES project cost will likely cost between $14,000,000 and $50,000,000. The potential flexibility of a two-tank sensible heat storage system should be able to demonstrate various new uses for nuclear heat. From power peaking via steam or liquid air heat topping to process heat applications, a two-tank system should be able to best provide a constant power transfer capability to the auxiliary energy consumer. Through commencement of this work, a more refined thermal storage selection process for Nuscale and existing nuclear power plants can be conducted. The ranking system developed here can be used to re-evaluate the energy storage options periodically depending on the projected installation dates. Current results suggest two-tank thermal energy storage is the most relevant technology in the near-term.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗