Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “criticality simulations”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Measurement Techniques for the High-Contrast and in-Field Stray-Light Performance of OCI

Scheduled for launch in January 2024, the PACE mission represents NASA’s next investment in ocean biology, clouds, and aerosol data records [1]. A key feature of PACE is the inclusion of an advanced satellite radiometer known as the Ocean Color Instrument (OCI), a global mapping radiometer that combines multispectral and hyperspectral remote sensing. A critical requirement for OCI is the high-contrast or spatial cross-talk specification (also referred to as in-field stray-light response). The requirement states that for global top-of-atmosphere radiances based on measured MODIS radiances, the global average residual contamination shall be less than 0.15% for any multispectral band. Accurate resolution of high contrast in TOA radiance images is important to estimate stray light contamination due to clouds, for studying small scale features like ocean fronts and for working in coastal and estuarine areas where the scales are 1km. This occurs in all wavelengths in the spatial direction. Knowledge of high contrast resolution makes up part of the artifact budget. Accurate measurement of the high-contrast performance of OCI requires laboratory Ground Support Equipment (GSE) that projects a scene of sufficient quality that the unwanted stray light of the GSE itself is not confused with the stray light response of the telescope. This paper concerns the development, analyses and test of the GSE to ensure the quality of the projected image is sufficient to verify the OCI requirements. Optical models were developed for both the instrument as well as the GSE and laboratory environment. Simulation of various non-ideal parameters were critical to accurately predict performance. Measurements using COTS cameras and lenses were also made of the projected GSE image to reasonably verify the optical model predictions. Measured and modelled results from OCI are discussed.

OCI, PACE, LSF, PSF, Stray-light, Contrast, Cross-↗
DOE OSTI2024

Performance Improvements of Poincaré Analysis for Exascale Fusion Simulations

Understanding the time-varying magnetic field in a fusion device is critical for the successful design and construction of clean-burning fusion power plants. Poincaré analysis provides a powerful method for the visualization of magnetic fields in fusion devices. However, Poincaré plots can be very computationally expensive making it impractical, for example, to generate these plots in situ during a simulation. In this short paper, we describe a collaboration among computer science and physics researchers to develop a new Poincaré tool that provides a significant reduction in the time to generate analysis results.

Pugmire, Dave↗
DOE OSTI2020

Parallel-in-Time Simulation of Lindblad's Equation

Constructing fast quantum logic gates is critical to building a scalable quantum computer. We consider a qudit, a quantum version of a bit that can take an arbitrary number of states, coupled with a cavity. In this project, we wish to force the qudit to reach the 0-state, for any possible initial state. The coupled system changes in time according to Lindblad’s equation, an ordinary differential equation on the density matrix of the quantum system. Lindblad’s equation contains some parameters that we can control, so-called control functions. We seek control functions which force the qudit to the 0-state within 2 microseconds, which is much faster than what is currently done in practice. The search method is gradient descent, a numerical optimization method that uses gradient information to iteratively improve the control parameters. My contribution to this project is an attempt to speed up the computation of the gradient. It currently takes about 40 seconds to compute the gradient which involves solving a set of ODEs sequentially. Current supercomputers have thousands of cores, but sequential computations can only make use of 1 core at a time. We wish to divide up the work better, so that we can use many more cores at once. To this end, we have implemented the Multigrid Reduction in Time (MGRIT) algorithm. We perform a systematic parameter search on how to best apply this algorithm. Results indicate a 25 percent speed up for solving Lindblad’s equation and determining how close the final state is the 0-state.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Equilibrium expectations for non-Gaussian fluctuations near a QCD critical point

With the highly anticipated results from the Beam Energy Scan II program at RHIC being recently revealed, an understanding of particle-number fluctuations and their significance as a potential signature of a possible QCD critical point is crucial. Early works that embarked on this endeavor sought to estimate the fluctuations due to the presence of a critical point assuming they stay in equilibrium. From these results came the proposal to focus efforts on higher, non-Gaussian, moments of the event-by-event distributions, in particular of the number of protons. These non-Gaussian moments are especially sensitive to critical fluctuations, as their magnitudes are proportional to high powers of the critical correlation length. As the equation of state provides key input for hydrodynamical simulations of heavy-ion collisions, we estimate equilibrium fluctuations from the BEST equation of state (EoS) that includes critical features from the 3D Ising Model. In particular, the proton factorial cumulants and their dependence on non-universal mapping parameters is investigated within the BEST EoS. Furthermore, the correlation length, as a central quantity for the assessment of fluctuations in the vicinity of a critical point, is also calculated in a consistent manner with the scaling equation of state. An understanding of the equilibrium estimates of proton factorial cumulants will be useful for further comparison to estimates of out-of-equilibrium fluctuations in order to determine the magnitude of the observable fluctuations to be expected in heavyion collision experiments, in which the time spent near a critical point is short.

Karthein, Jamie M. [Massachusetts Institute of Tec↗
NASA NTRS2001

An Improved Model for the Turbulent PBL

Second order turbulence models of the Mellor and Yamada type have been widely used to simulate the PBL. It is however known that these models have several deficiencies. For example, they all predict a critical Richardson number which is about four times smaller than the Large Eddy Simulation (LES) data, they are unable to match the surface data, and they predict a boundary layer height lower than expected. In the present model, we show that these difficulties are all overcome by a single new physical input: the use of the most complete expression for both the pressure-velocity and the pressure-temperature correlations presently available. Each of the new terms represents a physical process that, was not accounted for by previous models. The new model is presented in three different levels according to Mellor and Yamada's terminology, with new, ready-to-use expressions for the turbulent, moments. We show that the new model reproduces several experimental and LES data better than previous models. As far as the PBL is concerned, we show that the model reproduces both the Kansas data as analyzed by Businger et al. in the context of Monin-Obukhov similarity theory for smaller Richardson numbers, as well as the LES and laboratory data up to Richardson numbers of order unity. We also show that the model yields a higher PBL height than the previous models.

Cheng, Y.↗
NASA NTRSDate not supplied

Scalable Variable Charge Molecular Dynamics Simulations of Metal-Oxide Systems

Interfaces between metals and oxides are an important feature in many technologically relevant materials, e.g., oxidation of metal surfaces, oxide-dispersion strengthened (ODS) alloys, dielectric components, and thermal barrier coatings among others. Experimental studies of such interfaces are challenging since the majority are buried within the bulk, making computational modeling an attractive alternative. Molecular dynamics (MD) simulations operate at the length scales relevant to many interface-mediated mechanisms, but the requisite interatomic potentials for metal-oxide systems require computationally expensive variable charge schemes to account for the disparate bonding types, thus often limiting their effectiveness. Here we introduce several improvements to the charge transfer interatomic potential (CTIP) model which enable greater computational efficiency for large scale MD simulations. Then, using a new CTIP parametrization for the Ni-O system, we demonstrate its capabilities to capture critical atomic scale mechanisms associated with metal-oxide interfaces. Long time scale simulations (>10 ns) are used to investigate high temperature oxidation and oxide precipitation from the melt, and large length scale simulations (> 1 million atoms) are used to study the interaction of dislocations with oxide particles. We have implemented the new CTIP model in the widely used, open-source MD code LAMMPS.

Gabriel Plummer↗
NASA NTRS2004

Reducing the Time and Cost of Testing Engines

Producing a new aircraft engine currently costs approximately $1 billion, with 3 years of development time for a commercial engine and 10 years for a military engine. The high development time and cost make it extremely difficult to transition advanced technologies for cleaner, quieter, and more efficient new engines. To reduce this time and cost, NASA created a vision for the future where designers would use high-fidelity computer simulations early in the design process in order to resolve critical design issues before building the expensive engine hardware. To accomplish this vision, NASA's Glenn Research Center initiated a collaborative effort with the aerospace industry and academia to develop its Numerical Propulsion System Simulation (NPSS), an advanced engineering environment for the analysis and design of aerospace propulsion systems and components. Partners estimate that using NPSS has the potential to dramatically reduce the time, effort, and expense necessary to design and test jet engines by generating sophisticated computer simulations of an aerospace object or system. These simulations will permit an engineer to test various design options without having to conduct costly and time-consuming real-life tests. By accelerating and streamlining the engine system design analysis and test phases, NPSS facilitates bringing the final product to market faster. NASA's NPSS Version (V)1.X effort was a task within the Agency s Computational Aerospace Sciences project of the High Performance Computing and Communication program, which had a mission to accelerate the availability of high-performance computing hardware and software to the U.S. aerospace community for its use in design processes. The technology brings value back to NASA by improving methods of analyzing and testing space transportation components.

Source record↗
DOE OSTI2024

Deep Multi-Agent Reinforcement Learning for Real-World Signalized Traffic Corridor Control

Signalized traffic control problem has been addressed recently with deep Reinforcement Learning (RL) approaches involving diverse state, action, and reward structures. While significant progress has been noted in the literature, open challenges still remain in the areas of adaptive signal phase timing, coordination in a multi-intersection corridor setting, and consideration of real-world traffic conditions. In the context of deep RL-based problem framing, extensions are needed that enable adaptive signal phase timings in an intersection agent's action space, computationally efficient information sharing among neighboring signalized intersection agents along a corridor, and experimentation in realistic simulation environments. In this paper, we develop a deep Advantage Actor Critic (A2C) multi-agent RL (MARL) approach capturing the research extensions above and apply it within a real-world calibrated Aimsun Next traffic corridor simulation model based on traffic data from the City of Coral Gables, Florida. For a multi-intersection corridor control setting, our numerical simulation experiments with a decentralized A2C MARL algorithm applied at different time periods led to a total average corridor travel delay reduction (expressed in seconds/mile averaged over vehicles) from 4.9% to 19.9% compared to state-of-the-art actuated control.

Shuvo, Salman S. [BATTELLE (PACIFIC NW LAB)]↗
NASA NTRS1979

Effects of space environment on composites: An analytical study of critical experimental parameters

A generalized methodology currently employed at JPL, was used to develop an analytical model for effects of high-energy electrons and interactions between electron and ultraviolet effects. Chemical kinetic concepts were applied in defining quantifiable parameters; the need for determining short-lived transient species and their concentration was demonstrated. The results demonstrates a systematic and cost-effective means of addressing the issues and show qualitative and quantitative, applicable relationships between space radiation and simulation parameters. An equally important result is identification of critical initial experiments necessary to further clarify the relationships. Topics discussed include facility and test design; rastered vs. diffuse continuous e-beam; valid acceleration level; simultaneous vs. sequential exposure to different types of radiation; and interruption of test continuity.

Gupta, A.↗
DOE OSTI2021

Manifestations of metastable criticality in the long-range structure of model water glasses

Abstract Much attention has been devoted to water’s metastable phase behavior, including polyamorphism (multiple amorphous solid phases), and the hypothesized liquid-liquid transition and associated critical point. However, the possible relationship between these phenomena remains incompletely understood. Using molecular dynamics simulations of the realistic TIP4P/2005 model, we found a striking signature of the liquid-liquid critical point in the structure of water glasses, manifested as a pronounced increase in long-range density fluctuations at pressures proximate to the critical pressure. By contrast, these signatures were absent in glasses of two model systems that lack a critical point. We also characterized the departure from equilibrium upon vitrification via the non-equilibrium index; water-like systems exhibited a strong pressure dependence in this metric, whereas simple liquids did not. These results reflect a surprising relationship between the metastable equilibrium phenomenon of liquid-liquid criticality and the non-equilibrium structure of glassy water, with implications for our understanding of water phase behavior and glass physics. Our calculations suggest a possible experimental route to probing the existence of the liquid-liquid transition in water and other fluids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Exploring Wildfire & Energy data toward State Prioritization Index (WESPI)

Energy infrastructure can both induce and suffer risks from wildfires ranging from direct damage to energy assets such as substations and power lines to Public Safety Power Shutoffs. Recent wildfire events underscore the need for data-driven approaches that help states and utilities proactively plan for wildfire risk. Existing national tools such as Federal Emergency Management Agency (FEMA)’s National Risk Index (NRI) are valuable for community hazard planning. However, they are less suited for energy infrastructure, as they emphasize population and building exposure rather than system vulnerabilities. In this paper, we explore relationships between energy and wildfire data and present a Wildfire-Energy State Prioritization Index (WESPI). Our methodology combines data from the US Forest Service’s Fire Simulation (FSIM) dataset with energy resilience metrics, historical fire incidents, and geospatial data on transmission lines and fire stations. Correlation analyses suggest that FSIM burn probability is more strongly associated with power outage metrics (ρ = 0.32) than NRI wildfire frequency, and counties with a greater density of fire stations experience more frequent, but less intense wildfires. We further leverage data for burn probability, transmission line density, and fire station density to develop a Wildfire-Energy State Prioritization Index (WESPI) to highlight counties where wildfire hazard, infrastructure exposure, and limited suppression capacity converge. The index provides a consistent, scalable framework for state energy offices and utilities to screen counties for vegetation management, optimization of outage management system deployment, and to inform wildfire mitigation plans.

Critical infrastructure↗
NASA NTRS1990

Cryogenic liquid hydrogen reorientation activated by constant reverse gravity acceleration of geyser initiation

A key objective for cryogenic fluid management in a spacecraft propulsion system is development of the technology necessary for acquisition or positioning of liquid outflow or vapor venting. Numerical simulation of positive liquid acquisition is attempted by introducing a reverse gravity acceleration from the propulsive thrust of auxiliary engines which exceeds critical value for the initiation of a geyser. Based on the computer simulation of flow fields during the course of fluid reorientation, six dimensionless parameters resulted. These parameters hold near-constant values through the entire range of liquid filled levels, from 30 to 80 percent, during the course of fluid reorientation.

Hung, R. J.↗
DOE OSTI2023

Fluid-phase helium: Shock-compression experiments, quantum molecular dynamics simulations, and development of an equation of state

Helium (He) plays a critical role in numerous areas ranging from the study of celestial objects like brown dwarfs and gas giants to modern-day technologies like nuclear energy and rocket propulsion. For many of these applications, it is essential to have a reliable equation of state (EOS) for He that yields an accurate representation of its thermodynamic behavior. To help constrain and develop such EOS models, we have performed a series of shock-compression experiments on cryogenic liquid He to pressures exceeding 100 GPa using a magnetically accelerated flyer plate on Sandia National Laboratories' Z-machine. We have also performed quantum molecular dynamics simulations that are consistent with our shock measurements. None of the previously available EOSs agree with our experimental and simulation results, motivating the development of a fluid-phase He EOS that we present in this study. Here, we show that our EOS yields good agreement with published data that span temperatures and pressures encountered across a diverse array of applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

The effects of pre-plasma scale length on the relativistic electron beam directionality

The effects of pre-plasmas on the electron beam directionality was experimentally and numerically investigated. Single material and layered targets made of Ti and/or CH were used to simultaneously measure high-energy (≥3 MeV) electrons along two directions, pre-pulse energy and pre-plasma density. The electron directionality is quantified by using a new parameter, the electron energy ratio of the total kinetic energies along the two directions. Measurements and radiation–hydrodynamic (RH) simulations show that a large (≥3.5 μm) plasma scale length at the critical surface enhances electrons along the laser axis, and such pre-plasma conditions could only be achieved with the CH targets. Particle-in-cell simulations were performed on the RH generated pre-plasmas from Ti and CH targets, and the results show that the CH target provided conditions for higher forward momentum gains by electrons. First, the CH target allowed longer distances for electrons to interact with laser. Second, the intense laser pulse modified the critical surface, but the resulting surface differed. The CH target resulted in a smooth surface where a retro-reflection was observed while the Ti target resulted in a rippled surface that scattered the reflected light. As results, the CH electrons gained higher forward momentum via a direct-laser-acceleration in the counter propagating laser fields. So the results presented in this article show a way of controlling the high-energy electron directionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

gdess: A framework for evaluating simulated atmospheric CO 2 in Earth System Models

Atmospheric carbon dioxide (CO 2 ) plays a key role in the global carbon cycle and global warming. Climate-carbon feedbacks are often studied and estimated using Earth System Models (ESMs), which couple together multiple model components—including the atmosphere, ocean, terrestrial biosphere, and cryosphere—to jointly simulate mass and energy exchanges within and between these components. Despite tremendous advances, model intercomparisons and benchmarking are aspects of ESMs that warrant further improvement (Fer et al., 2021; Smith et al., 2014). Such benchmarking is critical because comparing the value of state variables in these simulations against observed values provides evidence for appropriately refining model components; moreover, researchers can learn much about Earth system dynamics in the process (Randall et al., 2019). We introduce `gdess` (a.k.a., Greenhouse gas Diagnostics for Earth System Simulations), which parses observational datasets and ESM simulation output, combines them to be in a consistent structure, computes statistical metrics, and generates diagnostic visualizations. In its current incarnation, `gdess` facilitates evaluating a model's ability to reproduce observed temporal and spatial variations of atmospheric CO 2 . The diagnostics implemented modularly in `gdess` support more rapid assessment and improvement of model-simulated global CO 2 sources and sinks associated with land and ocean ecosystem processes. We intend for this set of automated diagnostics to form an extensible, open source framework for future comparisons of simulated and observed concentrations of various greenhouse gases across Earth system models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A quantitative risk assessment framework for fault reactivation in underground hydrogen storage: Coupled simulation and deep learning approach

Underground hydrogen storage (UHS) is emerging as a critical solution for large-scale energy storage. However, like all subsurface fluid injection activities, UHS poses the risk of injection-induced fault reactivation. Accurate risk assessment is essential to ensuring the safety and efficiency of UHS operations. This study presents the development of deep-learning surrogate models for fault reactivation prediction in UHS, trained on a comprehensive database of fully coupled fluid flow-geomechanics simulations. Our findings reveal that analytical models often yield unreliable estimates, with errors up to 54% in the allowable injection pressure, potentially leading to a 40% reduction in UHS operational capacity. The developed surrogate models were incorporated into a quantitative risk assessment (QRA) framework, enabling probabilistic evaluation of fault reactivation risk while accounting for uncertainties in the input variables. Site-specific features, such as horizontal stress gradients, fault’s dip and strike angles, and operational parameters like bottom-hole injection pressure and well-fault distance, were identified as the primary drivers of fault reactivation across various stress regimes. Whereas other hydraulic, geological, and poroelastic reservoir properties were found to have a secondary impact. Notably, we observed that the risk of fault reactivation for a critically oriented fault with a static friction coefficient greater than 0.55 remains below 10% in a normal faulting stress regime. However, the risk significantly increases as the stress regime transitions from normal to strike-slip and ultimately to reverse faulting conditions. These findings underscore the importance of rigorous site characterization and comprehensive QRA evaluations to optimize UHS performance and minimize geomechanical risks.

25 ENERGY STORAGE↗
DOE OSTI2024

Application of the Cohn-alpha Method on Bare Highly Enriched Uranium Using Organic Scintillators [Poster]

The prompt neutron decay constant (α), can be used to calculate the lifetime of prompt fission chains within a multiplying system and verify reactor startup and shutdown. We analyze a fast, bare system consisting of highly enriched uranium (93% 235 U) hemisphere shells known as the Measurements of Uranium Subcritical and Critical (MUSiC). We measured the system with a three by-four array of organic scintillators (OSCAR) to obtain the gamma-ray and fast neutron signals. Pulse shape discrimination was used to isolate the neutron time series data. The Cohn-α method was used to estimate α for 3 supercritical measurements of Configurations 8 MUSiC. We use these estimates to linearly extrapolate α at delayed critical. These values can be compared to highly detailed simulations in MCNP6.2.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1983

New perspectives for advanced automobile diesel engines

Computer simulation results are presented for advanced automobile diesel engine performance. Four critical factors for performance enhancement were identified: (1) part load preheating and exhaust gas energy recovery, (2) fast heat release combustion process, (3) reduction in friction, and (4) air handling system efficiency. Four different technology levels were considered in the analysis. Simulation results are compared in terms of brake specific fuel consumption and vehicle fuel economy in km/liter (miles per gallon). Major critical performance sensitivity areas are: (1) combustion process, (2) expander and compressor efficiency, and (3) part load preheating and compound system. When compared to the state of the art direct injection, cooled, automobile diesel engine, the advanced adiabatic compound engine concept showed the unique potential of doubling the fuel economy. Other important performance criteria such as acceleration, emissions, reliability, durability and multifuel capability are comparable to or better than current passenger car diesel engines.

Tozzi, L.↗