SEARCH · Search NASA
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.
popclass is a lightweight python package that allows fast, probabilistic classification of the lens of a microlensing event given the event's posterior distribution and a model of the Galaxy. popclass provides the bridge between Galactic simulation and lens classification, an interface to common Bayesian inference libraries, and the ability for users to flexibly specify their own Galactic model and classification parameters.
The Fermilab Linac, a pivotal and historic accelerator at Fermilab, is crucial to the laboratory's operations, supplying a 400 MeV beam to various acceleration facilities. Due to daily variations in the Linac's output energy, precise monitoring and machine tuning are important to ensure the exiting energy meets the Booster's acceptance criteria. To address this, Beam Position Monitors (BPMs) are employed to assess and adjust the beam's energy, providing essential data on both the transverse position and longitudinal phase of the beam. We also aim to regulate the longitudinal phase profile to match our simulation predictions. By utilizing a Python-based simulation and beam data from three BPMs, we can determine the optimal phasing correction required for the final RF stage to achieve the desired energy. Currently, the calculations of longitudinal phase profile are based on simulations, so additional research is needed to verify if the predicted longitudinal phase distribution aligns with actual real-world data.
Methane (CH 4 ) emissions from wetland ecosystems are controlled by redox conditions in the soil, which are currently underrepresented in Earth system models. Plant-mediated radial oxygen loss (ROL) can increase soil O 2 availability, affect local redox conditions, and cause heterogeneous distribution of redox-sensitive chemical species at the root scale, which would affect CH 4 emissions integrated over larger scales. In this study, we used a subsurface geochemical simulator (PFLOTRAN) to quantify the effects of incorporating either spatially homogeneous ROL or more complex heterogeneous ROL on model predictions of porewater solute concentration depth profiles (dissolved organic carbon, methane, sulfate, sulfide) and column integrated CH 4 fluxes for a tidal coastal wetland. From the heterogeneous ROL simulation, we obtained 18% higher column averaged CH 4 concentration at the rooting zone but 5% lower total CH 4 flux compared to simulations of the homogeneous ROL or without ROL. This difference is because lower CH 4 concentrations occurred in the same rhizosphere volume that was directly connected with plant-mediated transport of CH 4 from the rooting zone to the atmosphere. Sensitivity analysis indicated that the impacts of heterogeneous ROL on model predictions of porewater oxygen and sulfide concentrations will be more important under conditions of higher ROL fluxes or more heterogeneous root distribution (lower root densities). Despite the small impact on predicted CH 4 emissions, the simulated ROL drastically reduced porewater concentrations of sulfide, an effective phytotoxin, indicating that incorporating ROL combined with sulfur cycling into ecosystem models could potentially improve predictions of plant productivity in coastal wetland ecosystems.
Low-rank methods have shown success in accelerating simulations of a collisionless plasma described by the Vlasov equation, but still rely on computationally costly linear algebra every time step. We propose a data-driven factorization method using artificial neural networks, specifically with convolutional layer architecture, that trains on existing simulation data. At inference time, the model outputs a low-rank decomposition of the distribution field of the charged particles, and we demonstrate that this step is faster than the standard linear algebra technique. Numerical experiments show that the method achieves comparable reconstruction accuracy for interpolation tasks, generalizing to unseen test data in a manner beyond just memorizing training data; patterns in factorization also inherently followed the same numerical trend as those within algebraic methods (e.g., truncated singular-value decomposition). However, when training on the first 70% of a time-series data and testing on the remaining 30%, the method fails to meaningfully extrapolate. Despite this limiting result, the technique may have benefits for simulations in a statistical steady-state or otherwise showing temporal stability. These results suggest that while the model offers a computationally efficient alternative for datasets with temporal stability, its current formulation is best suited for interpolation rather than for predicting future states in time-evolving systems. This study thus lays the groundwork for further refinement of neural network-based approaches to low-rank matrix factorization in high-dimensional plasma simulations.
Stripped-envelope supernovae (SESNe) represent a significant fraction of core-collapse supernovae, arising from massive stars that have shed their hydrogen and, in some cases, helium envelopes. The origins and explosion mechanisms of SESNe remain a topic of active investigation. In this work, we employ radiative-transfer simulations to model the light curves and spectra of a set of explosions of single, solar-metallicity, massive Wolf–Rayet stars with ejecta masses ranging from 4 to 11 M ⊙ , which were computed from a turbulence-aided and neutrino-driven explosion mechanism. We analyze these synthetic observables to explore the impact of varying ejecta mass and helium content on observable features. We find that the light curve shape of these progenitors with high ejecta masses is consistent with observed SESNe with broad light curves but not the peak luminosities. The commonly used analytic formula based on rising bolometric light curves overestimates the ejecta mass of these high-initial-mass progenitor explosions by a factor of up to 2.6. In contrast, the calibrated method by Haynie et al., which relies on late-time decay tails, reduces uncertainties to an average of 20% within the calibrated ejecta mass range. Spectroscopically, the He I 1.083 μm line remains prominent even in models with as little as 0.02 M ⊙ of helium. However, the strength of the optical He I lines is not directly proportional to the helium mass but instead depends on a complex interplay of factors such as the 56 Ni distribution, composition, and radiation field. Thus, producing realistic helium features requires detailed radiative transfer simulations for each new hydrodynamic model.
Large-scale molecular dynamics of deuterium under shock can elucidate kinetic processes vital to the target design in inertial confinement fusion and high-energy-density experiments. However, modeling the complex evolution of this material from an insulating molecular state at ambient pressure to an ionized, atomic fluid under strong shock is beyond the capability of simple pair and even bond order potentials. We thus train a quantum-accurate and broadly transferable machine-learning interatomic potential for deuterium using the Chebyshev Interaction Model for Efficient Simulations framework. We show that due to an improved description of the molecular-to-atomic transition, our model is able to better reproduce the ab initio equation of state, radial distribution functions, and principal Hugoniot than bond order potentials. This represents an important step toward large-scale quantum-accurate and nonequilibrium simulations of complicated systems under dynamic changes including phase transitions.
The ever-increasing volume of data produced by HPC simulations necessitates scalable methods for data exploration and knowledge extraction. Scientific data analysis often involves complex queries across distributed datasets, requiring manipulation of multiple primary variables and generating derived data that needs to be handled efficiently, creating challenges for applications that need to parse many large datasets. Relying on individual applications to handle all intermediate data generally leads to redundant computations across studies and unnecessary data transfers. In this paper, we investigate the performance of different approaches where applications define derived variables as quantities of interest (QoIs) and offload the computation and transfer of these QoIs to the I/O library. This significantly reduces redundancy and optimizes data movement across the distributed storage and processing infrastructure by allowing control over when and where derived variables are computed. We present a detailed analysis of the performance-storage trade-offs associated with different solutions and showcase results for our study on two large-scale datasets created from climate and combustion simulations.
This work will compare the Allen method and Walker method of accounting for subhourly inverter clipping power losses in hourly PV performance models. The Allen method uses a matrix lookup based on DNI clearness and clipping potential to assign a clipping correction loss at each simulation timestep. The Walker method models the PV DC power input to the inverter as a distribution over the hourly timestep and uses integration over the timestep to determine the amount of clipping that occurs within the timestep. Both these models have been recently implemented in the System Advisor Model's (SAM) open-source code, and will be applied to hourly SURFRAD datasets to analyze the subhourly clipping loss predicted by each model for different system designs and inverter loading conditions. Both models will be compared to "true" 1-minute SURFRAD data simulations to see their accuracy against more accurate 1-minute clipping correction loss predictions. This model comparisons will be investigated in more detail at the PVSC conference in Seattle, Washington June 2024.
Inverter-based distributed energy resources (DERs) connected to distribution networks (DNs) can provide fast frequency support, but their reserve deliverability depends on feeder constraints and differs from synchronous primary frequency response (PFR). Existing transmission–distribution coordination studies usually treat reserve generically or neglect feeder-level feasibility, while frequency-security scheduling studies rarely represent distribution feeders explicitly. This paper develops a bi-level day-ahead scheduling framework for integrated transmission–distribution networks that jointly clears energy, transmission-side PFR, and distribution-side fast frequency response (FFR) under exogenous hourly inertia and largest-loss inputs from an external unit commitment (UC) schedule. The transmission problem is modeled with DC-optimal power flow (OPF) and closed-form second-order cone (SOC) frequency-security constraints, whereas each DN is represented by a reserve-aware branch-flow AC-OPF so that scheduled fast reserves remain deliverable during activation. The bi-level problem is reformulated through Karush–Kuhn–Tucker (KKT) conditions into a mixed-integer SOC program, and a penalty term is used to tighten the distribution-network relaxation. In the reduced test system, lower exogenous inertia increased the required primary response from 179.64 MW to 191.08 MW, distribution-side fast response reduced total frequency-response procurement by up to 4.9%, and neglecting distribution constraints overstated the combined distribution-side energy and reserve award by up to 18%. In the expanded study, the largest case was solved in 2.02 s with a 0.00% optimality gap. Time-domain simulations kept the frequency nadir above 59.0 Hz in all tested hours. These results demonstrate the value of fast-response modeling and distribution-feasible reserve delivery in coordinated market clearing.
Evaluating the mechanical response of fiber-reinforced composites can be extremely time-consuming and expensive. Machine learning (ML) techniques offer a means for faster predictions via models trained on existing input–output pairs and have exhibited success in composite research. This paper explores a fully convolutional neural network modified from StressNet, which was originally used for linear elastic materials, and extended here for a non-linear finite element (FE) simulation to predict the stress field in 2D slices of segmented tomography images of a fiber-reinforced polymer specimen. The network was trained and evaluated on data generated from the FE simulations of the exact microstructure. The testing results show that the trained network accurately captures the characteristics of the stress distribution, especially on fibers, solely from the segmented microstructure images. The trained model can make predictions within seconds in a single forward pass on an ordinary laptop, given the input microstructure, compared to 92.5 h to run the full FE simulation on a high-performance computing cluster. These results show promise in using ML techniques to conduct fast structural analysis for fiber-reinforced composites and suggest a corollary that the trained model can be used to identify the location of potential damage sites in fiber-reinforced polymers.
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.
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.
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.
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.
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.
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.
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.
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.