SEARCH · Engineering Papers

Results for “Simulation Results”

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.

SourceResearch connectionFrom yearThrough year

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

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Dynamics modeling of molten salt reactor with reduced and expanded representations of delayed neutron precursors

Molten salt reactors (MSRs) present unique challenges in dynamic behavior due to the mobility of their fuel. In these reactors, delayed neutron precursors (DNPs) drift with the fuel circulation through the primary loop. As a result, a fraction of DNPs decays outside the core, effectively reducing the available delayed neutron population for reactivity control. Consequently, precise modeling of the distribution and behavior of DNPs is critical for accurate reactor dynamics simulations. In this study, the System Dynamics Analysis Tool (SDAT) was used to simulate a thermal-spectrum MSR under steady-state conditions and following transients. The effects of using reduced and expanded representations of DNPs with fewer or more groups than the conventional 6-group model were investigated. Their impact on the simulated distribution of precursors in the primary loop, reactivity loss value, and reactor response to transients was analyzed. Simulation results showed that reduced models lead to the loss of the actual DNPs distribution data, resulting in less accurate estimates of reactivity loss. Reactor power predictions using these reduced models showed significant deviations compared to those using the conventional 6-group model in transient simulations. Expanded models offered a more accurate representation of the distribution of DNPs and reactivity loss estimates. Reactor power predictions using expanded models showed minimal deviation from the conventional 6-group model during the simulated transients.

analysis↗

DOE OSTI2020

WEC fault modelling and condition monitoring: A graph-theoretic approach

The nature of wave resources usually requires wave energy converter (WEC) components to handle peak loads (i.e., torques, forces, and powers) that are many times greater than their average loads, accelerating equipment degradation. Moreover, due to their isolated nature and harsh operating environment, WEC systems are projected to possess high operations and maintenance (O&M) cost, i.e., around 27% of their leveled cost of energy. As such, developing techniques to mitigate these costs through the application of condition monitoring and fault tolerant control will significantly impact the economic feasibility of grid connected WEC power. Toward this goal, models of faulty components are developed in the open source modeling platform, WEC-Sim, to estimate the performance and measurable states of a WEC operating with likely device and sensor failures. Two types of faulty component models are then applied to a point absorber WEC model with basic controller damping and spring forces. Resulting changes in device behavior are recorded as a benchmark, and a graph-theoretic approach is proposed for fault detection and identification utilizing multivariate time series. Simulation results demonstrate that these faults can greatly affect the WEC performance, and that the proposed method can effectively detect and classify different types of faults.

16 TIDAL AND WAVE POWER↗

DOE OSTI2026

Power-Hardware-in-the-Loop Experiments of a Microgrid with a Grid-Forming Battery Inverter

Microgrids continue to proliferate, and they are transitioning away from using conventional generating resources to increasingly relying on inverter-based resources (IBRs) as the voltage and frequency leaders. It is crucial to evaluate the capability of IBRs to provide microgrid stability and resilience. Hardware-in-the-loop (HIL) experiments were conducted to de-risk the field deployment of the San Diego Gas & Electric Company Borrego Springs Microgrid, where a battery inverter was upgraded with grid-forming (GFM) capability to serve as the island leader. This paper presents the HIL experimental results from an HIL test bed that uses a power-hardware-in-the-loop (PHIL) interface with a power inductor that was previously developed for PHIL simulations of microgrids where the inverters need to switch modes, i.e., between grid-following and GFM as the microgrid transitions between grid-connected and islanded operation. This paper presents more details on the interface and HIL simulation results of the planned islanding and load steps in islanded operation to show the effectiveness of the inverters in managing the voltage and frequency.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

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.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

PAMTRA 1.0: the Passive and Active Microwave radiative TRAnsfer tool for simulating radiometer and radar measurements of the cloudy atmosphere

Forward models are a key tool to generate synthetic observations given knowledge of the atmospheric state. In this way, they are an integral part of inversion algorithms that aim to retrieve geophysical variables from observations or in data assimilation. Their application for the exploitation of the full information content of remote sensing observations becomes increasingly important when these are used to evaluate the performance of cloud-resolving models (CRMs). Herein, CRM profiles or fields provide the input to the forward model whose simulation results are subsequently compared to the observations. This paper introduces the freely available comprehensive microwave forward model PAMTRA (Passive and Active Microwave TRAnsfer), demonstrates its capabilities to simulate passive and active measurements across the microwave spectral region for upward- and downward-looking geometries, and illustrates how the forward simulations can be used to evaluate CRMs and to interpret measurements to improve our understanding of cloud processes. PAMTRA is unique as it treats passive and active radiative transfer (RT) in a consistent way with the passive forward model providing upwelling and downwelling polarized brightness temperatures and radiances for arbitrary observation angles. The active part is capable of simulating the full radar Doppler spectrum and its moments. PAMTRA is designed to be flexible with respect to instrument specifications and interfaces to many different formats of input and output, especially CRMs, spanning the range from bin-resolved microphysical output to one- and two-moment schemes, and to in situ measured hydrometeor properties. A specific highlight is the incorporation of the self-similar Rayleigh–Gans approximation (SSRGA) for both active and passive applications, which becomes especially important for the investigation of frozen hydrometeors.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Power-Hardware-in-the-Loop Experiments of a Microgrid With a Grid-Forming Battery Inverter: Preprint

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2020

Biomass Combustion in a Circulating Fluidized Bed Combustor

Interest in circulating fluidized bed (CFB) boilers as a power generation technology has skyrocketed in recent years because of several advantages this technology offers over conventional boilers, such as increased gas-solid mixing resulting in higher combustion efficiency and the ability to use lower quality fuels. CFB combustors are operated at lower temperatures than conventional thermal power generation combustors, thus reducing NO x emissions. SO 2 emissions are conveniently controlled through the addition of Ca-based sulphur sorbents within the combustor. This report summarizes the current modeling effort on a 50 kW th CFB combustor with a diameter of 10 cm and a height of 5 m; designed, built, and operated at CanmetENERGY in Ottawa, Canada employing the multiphase particle-in-cell (PIC) approach in the open-source Multiphase Flow with Interphase eXchanges (MFiX) Software Suite. The MFiX-PIC model parameters for the simulation are tuned against cold-flow experiments from CanmetENERGY using 9 kg of olivine sand as the inert bed material. It is shown that for the relatively coarse fluid meshes and large parcel sizes necessitated by the scale of the simulation, filter size dependent corrections to the drag law must be incorporated to ensure accuracy of the simulation results. The validated cold flow model is extended to simulate reacting flow with torrefied hardwood as the feedstock and to validate the combustion reaction scheme. The species concentrations at the riser outlet are compared against CanmetENERGY’s experiments and show satisfactory agreement. The simulations demonstrate the ability of MFiX-PIC to accurately capture the physics and chemistry of a circulating fluidized bed combustor at bench scales, which can be further extended to pilot- and industrial-scale systems.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2023

Multiscale numerical investigation of ratchet growth damage effects in PBX 9502

This paper presents results of numerical experiments conducted on the high explosive PBX 9502 to investigate how recently observed grain-scale damage mechanisms of ratchet growth affect uniaxial compression measurements. Simulations are multiscale in the sense of directly resolving grains, pores, cracks, and grain-interfaces based upon scanning electron microscope (SEM) images of damaged and undamaged samples. The combined finite-discrete element method (FDEM) is utilized to resolve both grain-scale microfracture and elastoplastic deformation of solid grains. Pristine (undamaged) and damaged microstructures are compared in simulation of unconfined compression tests of the same material from the literature. Here, the simulation results show the observed microscale mechanisms of damage, specifically microfracture predominantly around and sometimes through grains and crack-associated pore growth, can well-explain the effective degradation of strength and stiffness observed in the laboratory measurements.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Optimal Droop Setting for Congestion Reduction in a 100% Grid-Forming Inverter-based Power System

he high penetration of inverter-based resources (IBRs) introduces new challenges to power systems due to the complex inverter control. However, IBRs can be configured to maximize their benefits to improve system resilience and reliability. This paper proposes a steady-state optimization model that aims to mitigate transmission congestion in a 100% grid- forming (GFM) IBR-based power system. This goal is achieved by determining the optimal droop settings for the GFM IBRs under different congestion conditions due to renewable energy and load variations. The numerical solution is rigorously verified by a high-fidelity model of the IEEE 39-bus test system with detailed GFM IBR control in the time-domain electromagnetic transient (EMT) simulation tool PSCAD. The numerical solution and simulation results show a significant congestion reduction while meeting all other operating requirements. It is also observed that the numerical solving time is substantially less compared to the EMT simulation time.

Nguyen, Quan H.↗

DOE OSTI2020

A robust dynamic state estimation approach against model errors caused by load changes

Dynamic state estimation (DSE) plays an important role in power system security monitoring and online control. In practice, there are two approaches to implementing DSE. The first approach is distributed DSE, which is based on the assumption that the terminal bus of each generator can be measured by PMUs (phasor measurement units). The assumption cannot be satisfied currently, however, because PMUs usually are installed at important high-voltage buses such as 500-kV buses installed in portions of the grid overseen by the Western Electricity Coordinating Council. Another issue of this approach is that performance of DSE is vulnerable to bad measurement data. The reason for this vulnerability is that DSE is performed separately through measurements at each terminal bus, and measurements at terminal buses are the only measurement upon which DSE can rely. Therefore, important redundant measurements are not included in this approach. The second approach is centralized DSE. This approach does not have the requirement for PMU location, and redundant measurements can be considered fully. However, load changes and grid topology changes impact centralized DSE. In this paper, we propose a new approach for handling the impact of load changes on DSE. We have developed a new algorithm that includes two sequential steps. In the first step, errors caused by load changes are detected by analyzing the difference between prediction results and measured results. In the second step, once model error is detected, a model optimization procedure is run to correct the error so the state estimation error can be mitigated. Simulation results from the IEEE 68 bus system show that the proposed approach can effectively handle model errors caused by load changes.

robust dynamic state estimation, load change, powe↗

DOE OSTI2023

Fuel Performance Analysis of Fast Flux Test Facility MFF-3 and -5 Fuel Pins Using BISON with Post Irradiation Examination Data

Using the BISON fuel-performance code, simulations were conducted of an automated process to read initial and operating conditions from the Pacific Northwest National Laboratory (PNNL) database and reports, which contain metallic-fuel data from the Fast Flux Test Facility (FFTF) MFF Experiments. This work builds on previous modeling efforts involving 1977 EBR-II metallic fuel pins from experiments. Coupling the FFTF PNNL reports to BISON allowed for all 338 pins from MFF-3 and MFF-5 campaigns to be simulated. Each BISON simulation contains unique power and flux histories, axial power and flux profiles, and coolant-channel flow rates. Fission-gas release (FGR), fuel axial swelling, cladding profilometry, and burnup were all simulated in BISON and compared to available post-irradiation examination (PIE) data. Cladding profilometry, FGR, and fuel axial swelling simulation results for full-length MFF metallic pins were found to be in agreement with PIE measurements using FFTF physics and models used previously for EBR-II simulations. The main two peaks observed within the cladding profilometry were able to be simulated, with fuel-cladding mechanical interaction (FCMI), fuel-cladding chemical interaction (FCCI), and thermal and irradiation-induced creep being the cause. A U-Pu-Zr hot-pressing model was included in this work to allow pore collapse within the fuel matrix. This allowed better agreement between BISON-simulated cladding profilometry and PIE measurements for the peak caused by FCMI. This work shows that metallic fuel models used to accurately represent fuel performance for smaller EBR-II pins may be used for full-length metallic fuel, such as FFTF MFF assemblies and the Versatile Test Reactor (VTR). As new material models and PIE measurements become available, FFTF MFF assessment cases will be reassessed to further BISON model development.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Numerical Analysis of Water Harvesting from Saturated Vapor with the Electrospray Water Droplet Injection in The Solar-Driven Desalination Applications

Solar-driven seawater desalination systems provide freshwater through environmentally friendly and carbon-neutral processes. Air humidification and dehumidification desalination (HDD) systems extract water directly from moist air. This paper focuses on a new solar-driven seawater desalination system that shares some common traits with HDD systems but in which the air is eliminated from the process. Seawater is vaporized by solar thermal radiation in high-performance solar panels. Water vapor flows to a new electro-condensation chamber and passes through a series of electrosprays, which inject small nuclei of freshwater droplets into the vapor clouds. Here, due to dielectrophoresis and electrohydrodynamic flows, the vapor molecules are captured by the charged droplets, and vapor is condensed at the droplet surface. The water vapor condensation was modeled by using the multiphase flow numerical model. The CFD model, which was also experimentally validated, used SprayFoam solver modules available in the Open FOAM® CFD open-source software. In preliminary results, 802 grams/(hr-m2) (grams per hour per unit area) of water was harvested directly from moist air at 24C and 90% relative humidity. The simulation results indicated that freshwater productivity could increase up to 4,000 grams/(hr-m2) for fully saturated air conditions. This represents about 5 times the water productivity rate of current conventional HDD systems today. In addition, energy consumption could be reduced by 80 percent compared to conventional HDD systems.

Electrospray, electro-condensation, solar, droplet↗

DOE OSTI2026

Classic and Quantum Task-Based Intelligent Runtime for QIRs Running on Multiple QPUs

High-performance computing systems are rapidly evolving into heterogeneous platforms that fuse quantum accelerators with traditional classical processing units (CPUs) and graphical processing units (GPUs). This convergence calls for runtimes capable of managing both classical and quantum workloads in a unified manner. We introduce an intelligent, task-based runtime that marries the Intelligent RuntIme System (IRIS) asynchronous scheduler with a quantum programming stack through the Quantum Intermediate Representation Execution Engine (QIR-EE). Our design allows programs written in the quantum intermediate representation (QIR) to be dispatched concurrently to a variety of back-ends, including multiple quantum simulators and nascent quantum processors, enabling genuine hybrid execution on a single node. To illustrate its practicality, we partition a 4-qubit and 20-qubit circuit into three sub-circuits using quantum circuit cutting via the QCut library. Each sub-circuit is simulated independently by the QIR-EE driver within IRIS, after which a classical post-processing step merges the simulation results to recover the outcome of the original full-circuit computation. This case study demonstrates how finer task granularity can enable the parallel execution and lower the simulation burden per quantum task while preserving overall accuracy, highlighting the feasibility of our hybrid approach.

Miniskar, Narasinga Rao [ORNL] (ORCID:000000018259↗

DOE OSTI2021

Propagating neutronic uncertainties for FFTF LOFWOS Test #13

The safety evaluation conducted for licensing a Sodium-cooled Fast Reactor (SFR) may require modeling transients with best-estimate calculation tools that must first be validated against real-world measurements. To provide the community with a valuable benchmarking opportunity for validating SFR analysis tools and methods, the International Atomic Energy Agency (IAEA) initiated a coordinated research project (CRP) in 2018 for the analysis of the Fast Flux Test Facility (FFTF) Loss of Flow Without Scram (LOFWOS) Test #13. The impact of nuclear data uncertainties on neutronics parameters was previously investigated based on the COMMARA-2.0 covariance matrix. Since the transient simulation results are very sensitive to certain reactivity coefficients, it was decided to employ rigorous uncertainty propagation methods to quantify the impact of nuclear data uncertainties on the best-estimate predication of FFTF LOFWOS Test #13. The DAKOTA code is used to propagate neutronic uncertainties through SAS4A/SASSYS-1 transient simulations, while taking into account spatial and reaction-wise correlations within these uncertainties. This study shows that the remaining discrepancies observed between the Argonne National Laboratory (ANL) best-estimate results and the experimental measurements can be partly explained by the uncertainty associated with Gas Expansion Module (GEM) worth, which contributes the majority of the overall nuclear data uncertainty on the output from the FFTF LOFWOS Test #13 transient simulation. This study also confirmed the importance of including spatial and reaction-wise correlations of nuclear data uncertainties on feedback coefficients in the uncertainty propagation to avoid under-estimating their impact during the transient simulations.

42 ENGINEERING↗

DOE OSTI2024

Prediction of The TREAT Facility Power and Temperature of Sirius-2c Nuclear Propulsion Fuel Experiment

This paper presents simulation results of the Sirius-2c full power experiment that was performed at the TREAT facility to support experiment modeling and transient analysis of the NASA-sponsored Sirius series of experiments for nuclear thermal propulsion applications. We presents a novel approach to predict the control rod motion and power of the reactor and the temperature of the specimen during transient experiments using the Griffin and Bison packages of the Multiphysics Object Oriented Simulation Environment (MOOSE). First, the model predicts the transient rod motion that results in the desired reactor power shape. Then, the temperature of the specimen is predicted based on the power deposited in the specimen. The prediction results are in good agreement with the measured reactor power and control rod position and over-predicted the specimen temperature which is mainly related to specimen surface emissivities. The reactor reaches a peak power of 80 MW and a specimen temperature of 2500 K and considering the right specimen emissivities the predicted temperature is in an excellent agreement with the measured values. This model can be utilized in helping and supporting the design and optimization a new experiments that are going to be irradiated in the TREAT facility.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Bayesian calibration of bubble size dynamics applied to CO 2 gas fermenters

To accelerate the scale-up of gaseous CO 2 fermentation reactors, computational models need to predict gas-to-liquid mass transfer which requires capturing the bubble size dynamics, i.e. bubble breakup and coalescence. However, the applicability of existing models beyond air–water mixtures remains to be established. Here, an inverse modeling approach, accelerated with a neural network surrogate, calibrates the breakup and coalescence closure models, that are used in class methods for population balance modeling (PBM). The calibration is performed based on experimental results obtained in a CO 2 -air–water-coflowing bubble column reactor. Bayesian inference is used to account for noise in the experimental dataset and bias in the simulation results. To accurately capture gas holdup and interphase mass transfer, the results show that the breakage rate needs to be increased by one order of magnitude. In conclusion, the inferred model parameters are then used on a separate configuration and shown to also improve bubble size distribution predictions.

09 BIOMASS FUELS↗

DOE OSTI2022

A phase-field model of quasi-brittle fracture for pressurized cracks: Application to UO 2 high-burnup microstructure fragmentation

Here, in this paper, we present a phase-field model of quasi-brittle fracture with pressurized cracks, with dedicated applications for polycrystalline materials. The model is formulated as a minimization problem within the variational framework. The external work done by pressure on the crack surfaces is included in the objective function. Several careful modeling choices lead to a regularization-length-independent critical strength. The model is constructed to give a softening response with an underlying linear traction-separation law. The pressure-dependent softening response and the regularization of the prescribed pressure are demonstrated with a (quasi) one-dimensional numerical analysis. In a two-dimensional numerical analysis under plane strain assumptions, the critical stress corresponding to crack propagation (as predicted by our quasi-brittle fracture model) is compared with linear elastic fracture mechanics (LEFM) analytical solutions. Our model is further utilized to simulate fission-gas-induced fragmentation of the UO 2 high-burnup structure (HBS). Simulation results show that pressurized bubbles can cause crack nucleation and propagation, and that the bubble size and the surrounding external pressure affect the critical pressure corresponding to crack nucleation. Simulations of a partial HBS at different recrystallization stages show that different grain structures (due to recrystallization) also influence crack paths and fragmentation morphology.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Evolution of the self-injection process in long wavelength infrared laser driven LWFA

In this work, long wavelength infrared laser-driven plasma wakefield accelerators are investigated in the self-modulated laser wakefield acceleration (SM-LWFA) and blowout regimes using 3D particle-in-cell simulations. The simulation results show that in the SM-LWFA regime, self-injection arises with wave breaking, whereas in the blowout regime, self-injection is not observed under the simulation conditions. The wave breaking process in the SM-LWFA regime occurs at a field strength that is significantly below the 1D wave-breaking threshold. This process intensifies at higher laser power and plasma density and is suppressed at low plasma densities (≤1×10 17 cm –3 here). The produced electrons show spatial modulations with a period matching that of the laser wavelength, which is a clear signature of direct laser acceleration.

97 MATHEMATICS AND COMPUTING↗