SEARCH · Search NASA

Results for “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.

563 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Quantum Circuits for the Preparation of Spin Eigenfunctions on Quantum Computers

The application of quantum algorithms to the study of many-particle quantum systems requires the ability to prepare wave functions that are relevant in the behavior of the system under study. Hamiltonian symmetries are important instruments used to classify relevant many-particle wave functions and to improve the efficiency of numerical simulations. In this work, quantum circuits for the exact and approximate preparation of total spin eigenfunctions on quantum computers are presented. Two different strategies are discussed and compared: exact recursive construction of total spin eigenfunctions based on the addition theorem of angular momentum, and heuristic approximation of total spin eigenfunctions based on the variational optimization of a suitable cost function. The construction of these quantum circuits is illustrated in detail, and the preparation of total spin eigenfunctions is demonstrated on IBM quantum devices, focusing on three- and five-spin systems on graphs with triangle connectivity.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Improved constraints on hematite refractive index for estimating climatic effects of dust aerosols

Abstract Uncertainty in desert dust composition poses a big challenge to understanding Earth’s climate across different epochs. Of particular concern is hematite, an iron-oxide mineral dominating the solar absorption by dust particles, for which current estimates of absorption capacity vary by over two orders of magnitude. Here, we show that laboratory measurements of dust composition, absorption, and scattering provide valuable constraints on the absorption potential of hematite, substantially narrowing its range of plausible values. The success of this constraint is supported by results from an atmospheric transport model compared with station-based measurements. Additionally, we identify substantial bias in simulating hematite abundance in dust aerosols with current soil mineralogy descriptions, underscoring the necessity for improved data sources. Encouragingly, the next-generation imaging spectroscopy remote sensing data hold promise for capturing the spatial variability of hematite. These insights have implications for enhancing dust modeling, thus contributing to efforts in climate change mitigation and adaptation.

Environmental Sciences & Ecology↗

DOE OSTI2024

Three-chamber electrochemical reactor for selective lithium extraction from brine

Efficient lithium recovery from geothermal brines is crucial for the battery industry. Current electrochemical separation methods struggle with the simultaneous presence of Na + , K + , Mg 2+ , and Ca 2+ because these cations are similar to Li + , making it challenging to separate effectively. We address these challenges with a three-chamber reactor featuring a polymer porous solid electrolyte in the middle layer. This design improves the transference number of Li + (t Li+ ) by 2.1 times compared to the two-chamber reactor and also reduces the chlorine evolution reaction, a common side reaction in electrochemical lithium extraction, to only 6.4% in Faradaic Efficiency. Employing a lithium-ion conductive glass ceramic (LICGC) membrane, the reactor achieved high t Li+ of 97.5% in LiOH production from simulated brine, while the concentrations of Na + K + , Mg 2+ , and Ca 2+ are below the detection limit. Electrochemical experiments and surface analysis elucidated the cation transport mechanism, highlighting the impact of Na + on Li + migration at the LICGC interface.

Science & Technology - Other Topics↗

DOE OSTI2025

Phosphorylation toggles the SARS-CoV-2 nucleocapsid protein between two membrane-associated condensate states

Abstract The Nucleocapsid protein (N) of SARS-CoV-2 plays a critical role in the viral lifecycle by regulating RNA replication and by packaging the viral genome. N and RNA phase separate to form condensates that may be important for these functions. Both functions occur at membrane surfaces, but how N toggles between these two membrane-associated functional states is unclear. Here, we reveal that phosphorylation switches how N condensates interact with membranes, in part by modulating condensate material properties. Our studies also show that phosphorylation alters N’s interaction with viral membrane proteins. We gain mechanistic insight through structural analysis and molecular simulations, which suggest phosphorylation induces a conformational change in N that softens condensate material properties. Together, our findings identify membrane association as a key feature of N condensates and provide mechanistic insights into the regulatory role of phosphorylation. Understanding this mechanism suggests potential therapeutic targets for COVID infection.

Science & Technology - Other Topics↗

NASA NTRS2023

CFD Modeling of Bi-Directional PMD inside Cryogenic Propellant Tanks Onboard Parabolic Flights

Future cryogenic propulsion systems will require efficient methods with which to transfer cryogenic propellants from a depot storage tank to a customer receiver tank to minimize cost and maximize reusability. The Reduced Gravity Cryogenic Transfer project is currently developing advanced cryogenic fluid management technology and developing and validating new numerical models for three phases of transfer: line chilldown, tank chilldown, and tank fill. Additionally, multiple liquid nitrogen (LN 2 ) parabolic flight transfer rigs are being designed by universities and NASA to investigate the gravitational sensitivities that exist in these three technologies. In order to maximize the collection of low-g data during flights, it is required to extract as much (LN 2 as possible from the supply tank, despite variable gravity levels. The purpose of this paper is to present computational fluid dynamics (CFD) volume of fluid simulations of (LN 2 behavior in the supply tank onboard parabolic flights to validate the optimal design of a bi-directional propellant management device (PMD) using the commercial software FLOW-3D. A parametric study is conducted on the effects of gravity level, fill level, pore size, open area, thickness, and type of baffle on PMD performance. Based on results, the PMD as designed exceeds the targeted expulsion efficiency.

Jason Hartwig↗

DOE OSTI2026

The CP-PAW Code Package for First-Principles Calculations from a User’s Perspective

CP-PAW is a combined electronic structure and ab initio molecular dynamics code to perform mixed quantum and classical simulations of atomistic condensed phase systems, such as solids, liquids, and molecular systems. As the name suggests, the CP-PAW code unifies the all-electron projector augmented-wave (PAW) method with the Car–Parrinello (CP) approach to determine not only the electronic and nuclear ground states of condensed matter but also to study their properties and dynamics. In addition to briefly outlining the underlying theory, the focus will be on the unique aspects of CP-PAW and how to correctly employ them as a user. How to install CP-PAW using the new build system will also be briefly mentioned.

Blöchl, Peter E [Institute for Theoretical Physic↗

DOE OSTI2025

Entanglement Structure of Non-Gaussian States and How to Measure It

Rapidly growing capabilities of quantum simulators to probe quantum many-body phenomena require new methods to characterize increasingly complex states. Here, we present a protocol that constrains quantum states using experimentally measured correlation functions. This method enables measurement of a quantum state’s entanglement structure, opening a new route to study entanglement-related phenomena. Our approach extends Gaussian state parameterizations by systematically incorporating higher-order correlations. We show the protocol’s usefulness in conjunction with current and forthcoming experimental capabilities, focusing on weakly interacting fermions as a proof of concept. Here, the lowest nontrivial expansion quantitatively predicts early time thermalization dynamics, including signaling the onset of quantum chaos indicated by the entanglement Hamiltonian.

Fermi gases↗

DOE OSTI2026

Conformational Isomerization of Imide Anions Governs Solvation and Transport in Water-in-Salt Electrolytes

The behavior of highly concentrated electrolytes departs radically from the dilute-solution theory, yet the molecular origin of this transformation remains unresolved. Here, we identify the conformational isomerization of molecular ions as a decisive, previously unrecognized control parameter governing structure and transport in crowded aqueous electrolytes. Across a series of fluorosulfonimide anions, we show that increasing concentration drives a collective shift from extended transoid to compact cisoid conformers, revealed by small-angle X-ray scattering, vibrational spectroscopy, pulsed-field gradient NMR, and molecular dynamics simulations. This conformational transition triggers a collapse of the hydrogen-bonded water network and the emergence of densely packed ionic domains with confined water, producing a qualitative change in Li+ transport from solvent-mediated diffusion to network-confined hopping. Anion size and asymmetry systematically tune the onset of this transition, demonstrating that molecular geometry dictates mesoscale organization and dynamics in the ion-rich regime. Our results establish ion conformation, not merely composition or coordination, as a fundamental thermodynamic variable in concentrated solutions, providing a chemical framework that unifies solvation structure and transport in water-in-salt electrolytes and suggesting new principles for designing dense ionic media.

Nguyen, Huong TD↗

DOE OSTI2024

Impact of hydrogenation on the stability and mechanical properties of amorphous boron nitride

Abstract Interconnect materials with ultralow dielectric constant, and good thermal and mechanical properties are crucial for the further miniaturization of electronic devices. Recently, it has been demonstrated that ultrathin amorphous boron nitride (aBN) films have a very low dielectric constant, high density (above 2.1 g cm −3 ), high thermal stability, and mechanical properties. The excellent properties of aBN derive from the nature and degree of disorder, which can be controlled at fabrication, allowing tuning of the physical properties for desired applications. Here, we report an improvement in the stability and mechanical properties of aBN upon hydrogen doping. With the introduction of a Gaussian approximation potential for atomistic simulations, we investigate the changing morphology of aBN with varying H doping concentrations. We found that for 8 at% of H doping, the concentration ofsp 3 -hybridized atoms reaches to a maximum which leads to an improvement of thermal stability and mechanical properties by 20%. These results will be a guideline for experimentalists and process engineers to tune the growth conditions of aBN films for numerous applications.

Materials Science↗

DOE OSTI2026

Stoichiometry dependent properties of cerium hydride: An active learning developed interatomic potential study

Cerium hydride has a variety of interesting properties, including a known lattice contraction and densification with increasing hydrogen content. However, precise stoichiometric control is not experimentally straightforward and ab initio approaches are not computationally feasible for many properties such as melting and low temperature diffusion. Therefore, we develop a machine-learned interatomic potential for cerium hydride that is valid for H to Ce ratios from 2.0 to 3.0. A query-by-committee active learning approach is used to develop the training set. Leveraging classical molecular dynamics simulations, we assess a range of properties and provide fundamental mechanisms for the trends with stoichiometry. Finally, a majority of the properties follow the trend of lattice contraction, being governed by the stronger lattice binding induced by adding octahedral atoms.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Terahertz conductivity of two-dimensional materials: a review

Two-dimensional (2D) van der Waals materials are shaping the landscape of next-generation devices, offering significant technological value thanks to their unique, tunable, and layer-dependent electronic and optoelectronic properties. Time-domain spectroscopic techniques at terahertz (THz) frequencies offer noninvasive, contact-free methods for characterizing the dynamics of carriers in 2D materials. They also pave the path toward the applications of 2D materials in detection, imaging, manufacturing, and communication within the increasingly important THz frequency range. In this paper, we overview the synthesis of 2D materials and the prominent THz spectroscopy techniques: THz time-domain spectroscopy, optical-pump THz-probe technique, and optical pump–probe THz spectroscopy. Through a confluence of experimental findings, numerical simulation, and theoretical analysis, we present the current understanding of the rich ultrafast physics of technologically significant 2D materials: graphene, transition metal dichalcogenides, MXenes, perovskites, topological 2D materials, and 2D heterostructures. Finally, we offer a perspective on the role of THz characterization in guiding future research and in the quest for ideal 2D materials for new applications.

2D materials↗

NASA NTRS2025

On the Stability Analysis of Astrophysical Cooling Functions

To model the temperature evolution of optically thin astrophysical environments at MHD scales, radiative and collisional cooling rates are typically either pretabulated or fit into a functional form and then input into MHD codes as a radiative loss function. Thermal balance requires estimates of the analogous heating rates, which are harder to calculate, and due to uncertainties in the underlying dissipative heating processes these rates are often simply parameterized. The resulting net cooling function defines an equilibrium curve that varies with density and temperature. Such cooling functions can make the gas prone to thermal instability (TI), which will cause departures from equilibrium. There has been no systematic study of thermally unstable parameter space for nonequilibrium states. Motivated by our recent finding that there is a related linear instability, catastrophic cooling instability, that can dominate over TI, here we carry out such a study. We show that Balbus instability criteria for TI can be used to define a critical cooling rate, Λc, that permits a nonequilibrium analysis of cooling functions through the mapping of TI zones. We furthermore illustrate how thermal conduction modifies the shape of TI zones. Upon applying a Λc-based stability analysis to coronal loop simulations, we find that loops undergoing periodic episodes of coronal rain formation are linearly unstable to catastrophic cooling instability, while TI is stabilized by thermal conduction.

Amanda Stricklan↗

DOE OSTI2025

Robust Spectral Anomaly Detection in EELS Spectral Images via 3D Convolutional Variational Autoencoders

Abstract A 3D Convolutional Variational Autoencoder (3D‐CVAE) is introduced for automated anomaly detection in electron energy‐loss spectroscopy spectrum imaging (EELS‐SI) data. This approach leverages the full 3D structure of EELS‐SI data to detect subtle spectral anomalies while preserving both spatial and spectral correlations across the datacube. By employing cross‐entropy loss and training on bulk spectra, the model learns to reconstruct bulk features characteristic of the defect‐free material. In exploring methods for anomaly detection, both the 3D‐CVAE approach and principal component analysis (PCA) are evaluated, testing their performance using FeL‐edge ΔEpeak shifts designed to simulate material defects. These results show that 3D‐CVAE achieves superior anomaly detection and maintains consistent performance across various shift magnitudes. The method demonstrates clear bimodal separation between bulk and anomalous spectra, enabling reliable classification. Further analysis verifies that lower‐dimensional representations are robust to anomalies in the data. While performance advantages over PCA diminish with decreasing anomaly concentration, our method maintains high reconstruction quality even in challenging, noise‐dominated spectral regions. This approach provides a robust framework for unsupervised automated detection of spectral anomalies in EELS‐SI data, particularly valuable for analyzing complex material systems.

Chemistry↗

DOE OSTI2026

Molecular dynamics based study on the effects of cation size on the local structure and diffusion in polymerized ionic liquids

We have used coarse-grained molecular dynamics simulations to understand the effect of cation size on ion diffusion in polymerized ionic liquids at temperatures well above the glass transition temperature of the polymers. We investigated dependencies of the diffusion constant on the cation radius and static dielectric constant and interpreted these results in terms of underlying structural changes, and decoupling phenomenon. We have found non-monotonic effects of the cation radius on the diffusion constant, with a maximum at an intermediate radius resulting from two different size dependent effects. Changes in the radial distribution functions characterizing spatial distribution of cations with respect to anions result in a monotonic increase in the coordination number. This yields an increase in the effective interaction energy as a function of cation radius, which can explain the decrease in the diffusion constant. This increase is counteracted at small cation radii by a corresponding increase in decoupling of cation and anion dynamics. In addition, the diffusion constant is found to increase on increasing the static dielectric constant for all radii of the cations. These results highlight the significant role of decoupling in designing polymerized ionic liquids with an enhanced diffusion constant of cations.

Gillespie, Colin [ORNL] (ORCID:0000000320287080)↗

DOE OSTI2026

UCB-GLOBES: An open-access mass spectral database of identified and unidentified atmospheric organic compounds

Chemical characterization of atmospheric organic aerosols using gas chromatography with 70 eV electron ionization mass spectrometry (GC/EI-MS) has been used for decades in advancing molecular marker detection and identification, though primarily through suspect screening and/or targeted analyses. To advance non-targeted analyses of environmental samples, we have catalogued approximately 27 000 mass spectra (MS) of the trimethylsilyl derivatives of semi-volatile organic aerosol (OA) analytes in the open-access University of California Berkeley Goldstein Library of Organic Biogenic Environmental Spectra (UCB-GLOBES). Analytes were observed in ambient samples from the U.S. and the Central Amazon and/or laboratory simulations of secondary OA (SOA) formation. These samples are representative of OA under urban and biomass burning influences as well as SOA derived from biogenic precursors (e.g., isoprene, monoterpenes, sesquiterpenes) and biomass burning intermediates. MS are documented in UCB-GLOBES without regard to known chemical identity, annotated with extensive metadata such as sample source/experimental conditions, any structural information gained from MS analyses, and predicted chemical properties such as average carbon oxidation state and carbon number. UCB-GLOBES MS are compatible for importing into the NIST MS Search program, and we have also provided a Jupyter Notebook for MS visualization and comparisons. We demonstrate the utility of UCB-GLOBES through MS reanalyses of prior analytes observed in ambient data, finding a 20 % reduction in the number of analytes assigned to OA source categories reliant solely on time series correlation and an overall 11 % increase in new MS-based OA source categorization for the Southeast U.S. For 1513 analytes observed previously in the Central Amazon, we found 375 MS matches using UCB-GLOBES vs. 136 MS matches during prior analyses, representing a 14 % gain in newly confirmed or newly categorized OA species. While OA from laboratory oxidation experiments in UCB-GLOBES are highly diverse chemically, on average only 29 % of UCB-GLOBES MS have a mass spectral match to another MS entry in UCB-GLOBES and/or in databases of known compounds (i.e. NIST MS Database, Adams Essential Oil, MANE Flavor and Fragrance Company). This indicates that roughly 70 % of UCB-GLOBES MS are unique thus far, not observed more than once among the laboratory oxidation samples and ambient data in UCB-GLOBES MS. Further, only 18 % can be positively identified using these databases or known authentic standards. This points to a large gap between these laboratory simulations and ambient OA. Overall, the UCB-GLOBES database can be utilized for improving confidence in OA source categorization and/or identification, novel chemical marker discovery, tracking chemical diversity, de novo structure and properties prediction, and improving MS search and matching algorithms. This can ultimately inform future research priorities for the chemical characterization of atmospheric organic samples.

Mass spectrometry↗

DOE OSTI2026

Energy-Optimal Vehicle Longitudinal Motion Control via Pontryagin’s Minimum Principle and Ultra-Local Model

Longitudinal vehicle motion control is essential for enhancing performance and optimizing a vehicle’s energy usage. However, it remains a challenging task due to the nonlinear and uncertain nature of vehicle dynamics, along with varying driving conditions. This paper presents a novel ultra-local optimal control approach based on Pontryagin’s Minimum Principle (PMP) that circumvents the need for detailed system identification by employing an ultra-local model. The control objective is to minimize the total energy consumption under boundary conditions while ensuring smooth traction force generation. The proposed approach is evaluated using a high-fidelity vehicle model in three representative scenarios: (i) nominal driving, (ii) a change in tire road friction coefficient (TRFC) from 0.5 to 0.65 and road slope from 0% to 5% during the maneuver, with target velocity unchanged, and (iii) a change in target velocity from 20 m/s to 0 m/s during the maneuver, while maintaining nominal TRFC and slope conditions. The simulation results demonstrate that the proposed method delivers robust performance, effectively balancing consumption and tracking accuracy in all tested scenarios.

Waleed khan, Muhammad [The University of Texas at ↗

NASA NTRS2018

Experimental Effervescence and Freezing Point Depression Measurements of Nitrogen in Liquid Methane-Ethane Mixtures

NASA is designing an unmanned submarine to explore the depths of the hydrocarbon-rich seas on Saturn’s moon Titan. Data from Cassini indicates that the Titan north polar environment sustains stable seas of variable concentrations of ethane, methane, and nitrogen, with a surface temperature near 93 K. The submarine must operate autonomously, study atmosphere/sea exchange, interact with the seabed, hover at the surface or any depth within the sea, and be capable of tolerating variable hydrocarbon compositions. Currently, the main thermal design concern is the effect of effervescence on submarine operation, which affects the ballast system, science instruments, and propellers. Twelve effervescence measurements on various liquid methane-ethane compositions with dissolved gaseous nitrogen are thus presented from 1.5 bar to 4.5 bar at temperatures from 92 K to 96 K to simulate the conditions of the seas. After conducting effervescence measurements, two freezing point depression measurements were conducted. The freezing liquid line was depressed more than 15 K below the triple point temperatures of pure ethane (90.4 K) and pure methane (90.7 K). Experimental effervescence measurements will be used to compare directly with effervescence modeling to determine if changes are required in the design of the thermal management system as well as the propellers.

Freezing Point Depression↗

DOE OSTI2026

Surrogate modeling of Monte Carlo radiation transport with convolutional neural networks for shielding optimization

Here, we present a machine learning (ML)-based surrogate model using convolutional neural networks (CNN) designed to emulate the attenuation of neutron fields as they pass through various shielding materials. This model can compute the outgoing neutron flux almost instantaneously and achieves reasonable accuracy compared to traditional Monte Carlo (MC)-based codes, which are computationally intensive. This emulator alleviates the complexity of neutron radiation transport through shielding materials by reducing the dimensionality and enables shielding optimization for a known radiation environment. This optimization process, which would have taken an unrealistic timeline due to several complex radiation transport simulations, can now be achieved in minutes, thus increasing computational capabilities in radiation shielding assessment. We demonstrate the applications of this emulator in computing effective dose rates and optimizing shielding solutions for a heavy-ion accelerator facility, such as the Facility for Rare Isotope Beams, where secondary neutrons produced via beam interactions dominate the radiation environment.

accelerator shielding↗