SEARCH · Search NASA

Results for “distributed simulation”

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 577 records · Page 32

All sources NASA NTRS DOE OSTI

DOE OSTI2024

The Monolithic Heat Pipe Microreactor Reference Plant Model

This work introduces a reference plant model for a generic monolithic heat-pipe-cooled microreactor. The model will serve as a springboard to develop future evaluation models in the licensing process of similar microreactor designs at the U.S. Nuclear Regulatory Commission. This model has been developed with the Comprehensive Reactor Analysis Bundle and its specifications are based on open literature publications for the eVinci TM design. BlueCRAB is the U.S. Nu- clear Regulatory Commission non-light-water reactor analysis system based on the Multiphysics Object-Oriented Simulation Environment framework, which can couple the Griffin, BISON, and Sockeye applications to resolve the various physics that are essential for the safety analysis of this type of reactor system. The core specifications includes tristructural isotropic fuel, graphite monolith, graphite reflectors, and drums composed of graphite and B 4 C. No moderator or burnable poison pins are used in the design. The fuel enrichment is reduced to control excess reactivity in the core. This core design is not optimized and only serves for testing purposes, since the primary objective of this work is to exercise the multiphysics coupling for this type of reactor system. A three dimensional (3D) core heterogeneous Griffin discrete ordinates (SN) transport model allows the precise calculation of the flux distribution and pin powers. Griffin transfers the power density distribution and obtains a temperature distribution to and from BISON. The BISON model com- putes the 3D core temperature distribution and is coupled to 876 Sockeye subapplications running a heat pipe model. This 3D conduction model is coupled to the various heat pipes via heat flux boundary conditions. The model includes a small gap between the heat pipe and the monolith. Convective heat transfer boundaries with either ambient temperature or condenser temperature as heat sinks are imposed at the model boundaries. The 2D Sockeye heat pipe model uses a vapor- only methodology, which provides the needed resolution for transient calculations and allows the determination of various heat pipe limits. This approach is superior to the superconductor model traditionally used in steady-state calculations. BlueCRAB computes steady-state power and temperature distributions that serve as the initial condition for a loss-of-heat-sink transient simulation. The steady-state results show significant peaking due to the position of the control drum, but this is a characteristic of the particular design used, which is not optimized at this stage. The transient results show the reactor power slowly stabilizing towards a 3% power level after the partial loss of secondary heat removal. Several recriticalities are observed due to cooling through the secondary system but the reactor is self-stabilizing and behaves as expected.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2024

Simulation Evaluation of a Large-Scale Implementation of Virtual-Phase Link-Based Model Predictive Control

Traffic congestion is a serious problem in the US, and traffic signal control is one of the effective solutions to congestion. Previous research on model predictive control (MPC)-based traffic signal control showed substantial benefits over conventional methods. This study focused on implementing MPC over a large-scale network with complex intersections and the impact of cycle length, network size, and imperfect state estimation on performances. This study implemented a virtual phase link (VPL)-based model predictive control method which used the number of vehicles in each VPL as input state variables and was suitable for National Electrical Manufacturing Association (NEMA) ring-barrier control. To test the impact of network size, the performance of distributed MPC (36 intersections in the network are divided into five subnetworks) was compared with that of MPC over the full network for a set of cycle lengths. To test the impact of imperfect state estimation, we synthetically infused estimation error and developed two scenarios, MPC-error and MPC-error narrow, which had higher and lower estimation errors, respectively. The performance of these MPC methods was compared with that of the existing time-of-day (TOD) method and an offline method that used Webster's method for split and MULTIBAND for cycle length and offset optimization. Trajectory and linkwise signal performance measures were collected from the simulation to evaluate performance. The distributed MPC method with perfect state estimation had the lowest delay and highest energy efficiency of all the methods. The performance of MPC decreased as the prediction inaccuracy increased. MPC-error had 7% and 11% more delay than MPC-error narrow in the morning and evening peaks, respectively. Overall, simulation results suggest that even with imperfect state estimation, MPC methods will outperform offline methods significantly.

large-scale simulation↗

DOE OSTI2024

Data assimilation for combustion ignition delay time simulation using schlieren image velocimetry

This study sought to improve the accuracy of simulating spray penetration and combustion ignition delay by means of data assimilation (DA). The simulations were conducted using the Reynolds-averaged Navier–Stokes (RANS) equations and assimilating the schlieren image data. In DA, an ensemble square root filter (EnSRF) was used to build the statistical model, making the simulation results more accurate without any change in the governing equations. Recognizing that the spray-cone injection angle has a large effect on penetration, we created ensemble members with different injection angles. And we applied the two-component velocity distribution calculated via SIV and updated both velocity and temperature by using a DA statistical model derived from RANS ensemble simulations. The ignition delay time is generally known to vary even under the same experimental conditions because it is influenced by many factors. In this study, we attempted the transient DA-assisted RANS simulation to predict the ignition delay time even when the temporal resolution and accuracy of the observation data ware insufficient. Our trials offer an example of how a combination of techniques can be effectively used to assimilate experimental data obtained under restricted conditions.

Combustion simulation↗

DOE OSTI2025

Numerical framework for integrated additive manufacturing-compression molding (AM-CM) of thermoplastic composites

Additive manufacturing-compression molding (AM-CM) has emerged as a transformative technology in advanced composite manufacturing. Additive manufacturing (AM) offers high design flexibility and the ability to produce complex geometries with precisely aligned fibers in the preferred orientation. Compression molding (CM) enhances composite materials by providing excellent dimensional stability, reduced porosity, high production rates, and a smooth surface finish. Despite these advantages, extensive integrated analysis is required to optimize processing conditions for improved fiber orientation distribution (FOD) and porosity control. Here, this study develops a comprehensive numerical model to simulate the AM-CM manufacturing process. The model isolates the effects of both the AM and CM phases while also capturing their integration. Additionally, it accounts for heat transfer, temperature-dependent viscosity, and fiber orientation in the extruded fiber-filled polymer, accurately representing material behavior during processing. This approach enables the analysis of interactions between deposited beads of complex strand shapes and their interface regions after full compression. Moreover, the model predicts key parameters such as polymer flowability, fiber orientation, and temperature evolution in AM-CM parts. By optimizing processing conditions, it facilitates a controlled and predictable microstructure.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Study on Conversion from an Existing Alternating Current Medium-Voltage Line to a Direct Current Line: A Case of Load or Source Expansion

Simultaneous increase of power generation through renewable systems and load due owing to the rapid growth in electrified road transportation systems, leads to new challenges in the design, operation, and protection of distribution-level electrical power systems. This paper addresses challenges related to the integration of new sources or loads to a distribution grid through long alternating current lines. The International Council on Large Electric Systems medium-voltage benchmark model developed for North America was used for simulation studies. Use cases were developed to illuminate the issues related to long alternating current lines and the usefulness of medium-voltage direct current lines in a distribution network.

Jaldanki, Sreenivasa↗

DOE OSTI2025

Virtual Agents-Based Attack-Resilient Distributed Control for Islanded AC Microgrid

Due to its dependence on a communication network, distributed secondary control of microgrids is susceptible to denial-of-service (DoS) attacks in channel shutdown mode, which may negatively impact the network connectivity and thus deteriorate the coordination and power sharing among distributed generators (DGs). Honeypot is a common method for cyber deception by introducing fake targets. However, in the context of microgrid, the misleading information spread by honeypots will also impact the system performance. This paper proposes an attack-resilient distributed control for AC microgrids utilizing virtual agents (VAs) to counteract both DoS edge and node attacks. The VAs are designed to not impact the system’s steady state during normal operation but to share information among neighboring real agents and serve as dummy targets for DoS attacks. The control with VAs is implemented by a primal-dual gradient based distributed algorithm to efficiently obtain a practical solution for voltage/frequency regulation and power sharing. The simulation results on a 4-DG test system and a modified IEEE 34-bus system show that 1) VAs do not impact the normal functionality of the test system, and 2) deploying VAs can enhance the resilience of the microgrid control against DoS edge and node attacks.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2025

The phase diagram of quantum chromodynamics in one dimension on a quantum computer

The quantum chromodynamics (QCD) phase diagram, which reveals the state of strongly interacting matter at different temperatures and densities, is key to answering open questions in physics, ranging from the behaviour of particles in neutron stars to the conditions of the early universe. However, classical simulations of QCD face significant computational barriers, such as the sign problem at finite matter densities. Quantum computing offers a promising solution to overcome these challenges. Here, we take an important step toward exploring the QCD phase diagram with quantum devices by preparing thermal states in one-dimensional non-Abelian gauge theories. We experimentally simulate the thermal states of SU(2) and SU(3) gauge theories at finite densities on a trapped-ion quantum computer using a variational method. This is achieved by introducing two features: Firstly, we add motional ancillae to the existing qubit register to efficiently prepare thermal probability distributions. Secondly, we introduce charge-singlet measurements to enforce colour-neutrality constraints. This work pioneers the quantum simulation of QCD at finite density and temperature for two and three colours, laying the foundation to explore QCD phenomena on quantum platforms.

Quantum information↗

DOE OSTI2024

Learning-Based Building Flexibility Estimation and Control to Improve Microgrid Economics and Resilience: Preprint

This paper proposes a learning-based building flexibility estimation and control framework to improve system economics and resilience. A data-driven building load flexibility model consisting of weather forecasting and estimating load consumption is proposed to quantify building heating, ventilation, and air conditioning (HVAC) load flexibility. A reinforcement learning-based microgrid controller is proposed to dispatch distributed generators, distributed energy resources, and build HVAC loads while taking flexibility information as one of the inputs. Simulation analysis is conducted on the model of a real microgrid in California. The effectiveness of the proposed learning-based building flexibility estimation and control in reducing microgrid energy costs and improving the sustainability of critical loads is demonstrated.

building load flexibility↗

DOE OSTI2025

Porous and Amorphous Mn x Mo 3 S 13 Chalcogel Electrode for High-Capacity Conversion-Based Lithium-Ion Batteries

While Li-ion batteries (LIBs) are a leading energy storage technology, their energy densities are limited by the low capacity of conventional intercalation cathodes, driving interest in high energy-density Li–S batteries that make use of conversion chemistry. Achieving high capacity, reversibility, and cycle stability, and controlling volume changes in conversion batteries during the charge–discharge process, however, remains challenging. Here, in this study, we present a porous, amorphous, sulfide-based Mn x Mo 3 S 13 chalcogel, which concurrently offers high capacity and cycle stability. The solution-processable room temperature synthesized Mn x Mo 3 S 13 (x = 0.25) chalcogel exhibits a local structure that resembles the Mo 3 S 13 cluster with Mn 2+ distributed across the Mo 3 S 13 matrix, as determined by synchrotron X-ray pair distribution function (PDF) and extended X-ray absorption fine structure (EXAFS). Ab initio molecular dynamics (AIMD) simulations reveal that Mn 2+ incorporation shortens the polysulfide chain in the gel matrix compared to the Mo 3 S 13 chalcogel, while forming a coordination environment with disulfide groups, analogous to the experimental findings. A Li/Mn 0.25 Mo 3 S 13 half-cell delivers 897 mAh g –1 capacity during the first discharge and retains 571 mAh g –1 capacity after 100 cycles at a C/3 rate. Distribution of relaxation time (DRT) unveils a stable solid–electrolyte interphase (SEI) formation upon cycling that enables charge–discharge reversibility. Here, the enhanced capacity retention and cycle stability compared to those of the Li/Mo 3 S 13 cell are attributed to the reduced dissolution of active mass into the electrolyte, facilitated by the formation of shorter polysulfide chains within the Mn 0.25 Mo 3 S 13 structure and the strong affinity of Lewis-acidic Mn 2+ for polysulfide anions generated during the charge–discharge process of the Li/Mn 0.25 Mo 3 S 13 cell. Thus, this work illustrates a design principle of material for high-capacity and cycle-stable Li-metal sulfide batteries.

25 ENERGY STORAGE↗

DOE OSTI2024

How Does Vertical Wind Shear Influence Updraft Characteristics and Hydrometeor Distributions in Supercell Thunderstorms?

Abstract Vertical wind shear is known to affect supercell thunderstorms by displacing updraft hydrometeor mass downshear, thereby facilitating the storms’ longevity. Shear also impacts the size of supercell updrafts, with stronger shear leading to wider, less dilute, and stronger updrafts with likely greater hydrometeor production. To more clearly define the role of shear across different vertical layers on hydrometeor concentrations and displacements relative to supercell updrafts, a suite of idealized numerical model simulations of supercells was conducted. Shear magnitudes were systematically varied across the 0–1, 1–6, and 6–12 km AGL layers, while the thermodynamic environment was held fixed. Simulations show that as shear magnitude increases, especially from 1 to 6 km, updrafts become wider and less dilute with an increase in hydrometeor loading, along with an increase in the low-level precipitation area/rate and total precipitation accumulation. Even with greater updraft hydrometeor loading amid stronger shear, updrafts are more intense in stronger shear simulations due to larger thermal buoyancy owing to wider, less dilute updraft cores. Furthermore, downshear hydrometeor displacements are larger in environments with stronger 1–6-km shear. In contrast, there is relatively less sensitivity of hydrometeor concentrations and displacements to variations in either 0–1- or 6–12-km shear. Results are consistent across free tropospheric relative humidity sensitivity simulations, which show an increase in updraft size and hydrometeor mass with increasing free tropospheric relative humidity owing to a reduction in entrainment-driven dilution for wider updrafts in moister environments. Significance Statement Rotating thunderstorms, known as supercells, are able to persist for multiple hours. One common explanation is that large changes in wind speed and/or direction with height, or shear, transport rain/hail away from supercell updrafts, supporting their maintenance. The strong shear within supercell environments, however, may also lead to greater rail/hail amounts, thereby leading to weaker storms due to this extra mass of water/ice within updrafts. Furthermore, the impact of shear across different height layers on supercell rain/hail characteristics has not been thoroughly investigated. In this study, computer simulations of supercells were conducted to determine that shear occurring between 1 and 6 km above ground level has a large impact on rain/hail distribution in supercells and that stronger shear in this layer leads to wider/stronger supercells with greater rain/hail accumulations at the surface. Additionally, some of the extra mass of water/ice is transported farther away from updrafts due to the stronger environmental storm-relative winds.

Meteorology & Atmospheric Sciences↗

DOE OSTI2026

Thermophysical Properties of Liquid Tritium: A Path Integral Monte Carlo Study

Here, we present worm-algorithm, path integral Monte Carlo simulations of bulk liquid tritium. The simulations are benchmarked against empirically known thermophysical properties of liquid deuterium and liquid tritium. Results for the pair correlation function, chemical potential, isothermal compressibility, isochoric heat capacity, and single-particle momentum distributions are reported. Given the benchmark comparisons, our predictions of liquid tritium properties are expected to be accurate to within a few percent. Our simulations unambiguously demonstrate the significance of nuclear quantum effects to the properties of liquid tritium. In particular, under saturated vapor pressure, the average molecular kinetic energy of the liquid is found to be more than 60% higher than the value expected from the classical equipartition theorem.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Xyce™ Parallel Electronic Simulator Users' Guide (V.7.9)

This manual describes the use of the Xyce Parallel Electronic Simulator. Xyce has been designed as a SPICE-compatible, high-performance analog circuit simulator, and has been written to support the simulation needs of the Sandia National Laboratories electrical designers. This development has focused on improving capability over the current state-of-the-art in the following areas: • Capability to solve extremely large circuit problems by supporting large-scale parallel computing platforms (up to thousands of processors). This includes support for most popular parallel and serial computers. • A differential-algebraic-equation (DAE) formulation, which better isolates the device model package from solver algorithms. This allows one to develop new types of analysis without requiring the implementation of analysis-specific device models. • Device models that are specifically tailored to meet Sandia’s needs, including some radiation-aware devices (for Sandia users only). • Object-oriented code design and implementation using modern coding practices. Xyce is a parallel code in the most general sense of the phrase — a message passing parallel implementation — which allows it to run efficiently a wide range of computing platforms. These include serial, shared-memory and distributed-memory parallel platforms. Attention has been paid to the specific nature of circuit-simulation problems to ensure that optimal parallel efficiency is achieved as the number of processors grows.

42 ENGINEERING↗

DOE OSTI2025

Xyce™ Parallel Electronic Simulator Users’ Guide, Version 7.10

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

A new approach to the evaluation and solution of the relativistic kinetic dispersion relation and verification with continuum kinetic simulation

Here, the present work describes a new approach to evaluation and root finding for the kinetic dispersion relation of Langmuir waves, which is central to the analytical understanding of collisionless damping in plasmas. The plasma dispersion function is solved to machine precision using direct integration in the complex plane in combination with an analytic evaluation of the residue to account for the deformation along the Landau contour. To efficiently attain machine precision, the contour is displaced in the complex plane prior to integration, and numerical subtleties related to the placement of the contour are discussed. The approach is generic in that it applies to arbitrary distribution functions, with the present manuscript focused on relativistic cases. Detailed verification of results via direct kinetic simulation in a variety of configuration space dimensions is also presented. Finally, the technique is applied to the challenging case of highly relativistic (i.e. extremely hot) plasmas. Here we show both qualitative agreement with prior work, as well as the disappearance of the Landau root which would have significant implication for real-life observation or experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Q-Cluster: Quantum Error Mitigation Through Noise-Aware Unsupervised Learning

Quantum error mitigation (QEM) is critical in reducing the impact of noise in the pre-fault-tolerant era, and is expected to complement error correction in fault-tolerant quantum computing (FTQC). In this work, we propose a novel QEM approach, Q-Cluster, that uses unsupervised learning (clustering) to reshape the measured bit-string distribution. Our approach starts with a simplified bit-flip noise model. It first performs clustering on noisy measurement results, i.e., bit-strings, based on the Hamming distance. The centroid of each cluster is calculated using a qubit-wise majority vote. Next, the noisy distribution is adjusted with the clustering outcomes and the bitflip error rates using Bayesian inference. Our simulation results show that Q-Cluster can mitigate high noise rates (up to 40% per qubit) with the simple bit-flip noise model. However, real quantum computers do not fit such a simple noise model. To address the problem, we (a) apply Pauli twirling to tailor the complex noise channels to Pauli errors, and (b) employ a machine learning model, ExtraTrees regressor, to estimate an effective bit-flip error rate using a feature vector consisting of machine calibration data (gate & measurement error rates), circuit features (number of qubits, numbers of different types of gates, etc.) and the shape of the noisy distribution (entropy). Our experimental results show that our proposed Q-Cluster scheme improves the fidelity by a factor of 1.46x, on average, compared to the unmitigated output distribution, for a set of low-entropy benchmarks on five different IBM quantum machines. Our approach outperforms the state-of-art QEM approaches RZNE [28], M3 [24], Hammer [35], and QBEEP [33] by 1.26x,1.29x,1.47x, and 2.65 x, respectively.

42 ENGINEERING↗

DOE OSTI2025

Impact of super-Gaussian electron distributions on plasma K-shell emission

Electron distributions in laser-produced plasmas will be driven toward a super-Gaussian distribution due to inverse bremsstrahlung absorption [Langdon, Phys. Rev. Lett. 44, 575 (1980)]. Both theoretical and experimental evidence suggest that fundamental plasma properties are altered by the super-Gaussian distribution. Here, this paper examines how the super-Gaussian distribution affects the ionization balance and K-shell emission of atomic plasmas, utilizing approximate formulas and detailed collisional-radiative simulations. While the impact on plasma ionization is small, K-shell spectra can be significantly modified. Based on these findings, we demonstrate that K-shell spectroscopy can be used to infer super-Gaussian or other similar nonequilibrium electron distributions.

Atomic spectra↗

DOE OSTI2025

Analysis of Screening Current Effects in a Hybrid Nb$_{3}$Sn/REBCO Superconducting Accelerator Magnet Using a T - A Formulation

To explore the feasibility of using high-temperature superconducting (HTS) REBCO coated conductors in future accelerator magnets, two REBCO flat racetrack coils were fabricated using 4-mm wide EuBCO tapes at the High Energy Accelerator Research Organization (KEK). These coils were tested as an insert inside a Nb$_{3}$Sn common-coil dipole magnet, which provides a background field of up to $\sim$ 9.5T, at the Brookhaven National Laboratory (BNL). REBCO tapes offer exceptionally high critical current density under strong magnetic fields; however, they also exhibit significant magnetization due to screening currents, leading to magnetic field errors. Here, this study presents a 2D finite element model of screening current-induced fields (SCIF) in REBCO coils using the T-A formulation, along with the results obtained. Simulations were then performed for two KEK test cases: one where the REBCO conductors were oriented with the HTS tapes parallel to the background field, and another where the tapes were perpendicular to it. Since screening currents also influence the stress distribution and increase the peak stress in the coils, the mechanical effects of these currents were analyzed. The implications of these simulation and test results for the design of Nb$_{3}$Sn/REBCO superconducting accelerator magnets are discussed.

Accelerator magnets↗

DOE OSTI2024

Modeling the Effect of Material Properties on Liquid-Alkaline Water Electrolysis

Liquid-alkaline water electrolyzers (LAWEs) use electricity to drive the conversion of water to H 2 and O 2 gas. These devices benefit from the use of low-cost nickel electrodes and metal-oxide separators, but suffer from lower current densities and higher cell voltages than proton-exchange-membrane water electrolyzers. Identifying the inefficiencies that result in this poor performance is key to mitigating losses and optimizing LAWEs. Here, we report an experimentally-validated 1-D continuum model of a LAWE that elucidates the gradients within the cell, simulates H 2 crossover, and projects the energy improvements made possible by modulating the properties of the electrodes and separator. The model captures the Nernstian polarization losses and the distribution of gas- and liquid-phases within the electrodes, enabling quantification of energy losses associated with kinetic, ohmic, and bubble-induced (mass-transport) resistances. Simulations demonstrate that LAWE can achieve energy intensities of 50 kWh kg −1 of H 2 at 1 A cm −2 using improved electrode and separator properties.

08 HYDROGEN↗