Search NASA⌕ Search

SEARCH · Search NASA

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

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

At least 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Demonstration of an Automated System for Vertical Profiles of Volatile Organic Compounds

Volatile organic compounds (VOCs) play important roles throughout the atmosphere, many of which are altitude dependent. This highlights the need for easily deployable devices to sample VOCs across different atmospheric layers. To address this, we present the design and initial application of a Time Resolved Automated Volatile organIc compounds Sampling system (TRAVIS). VOCs are collected on sorbent tubes, which are subsequently analyzed by a thermal desorption gas chromatography mass spectrometry pipeline. TRAVIS leverages a piezoelectric pump with an integrated pressure sensor for precise (0.1% flow rate relative standard deviation) and accurate (−3 ± 2% error in VOC quantitation) measurements. Via deployment on a tethered balloon system over an agricultural area, TRAVIS is used to show consistent vertically resolved VOC concentrations in a well-mixed (i.e., turbulent) atmosphere (e.g., 5% relative standard deviation for isoprene) and vertically dependent concentrations for a stratified atmosphere (e.g., prior to boundary layer development). Furthermore, we also show VOC information from an intermittent plume via both targeted and untargeted analysis, highlighting future applications for spurious events in agriculture, air quality monitoring, and environmental impact. Overall, the development of TRAVIS represents a lightweight, accurate, sensitive, and precise VOC sampling module for the scientific community.

Aerosols↗
DOE OSTI2025

Advancing Aerosol Chemical Characterization and Vertical Profiling over the Southern Great Plains Using Uncrewed Aerial Sampling and Offline Aerosol Mass Spectrometry

Recent advancements in uncrewed aerial systems (UASs) and particulate matter (PM) analytical techniques have provided opportunities for atmospheric research. In this study, we deployed the Department of Energy’s fixed-wing ArcticShark UAS to examine PM 2.5 composition at varying altitudes─within and above the planetary boundary layer (PBL)─over the Southern Great Plains atmospheric observatory (SGP). A total of 22 flights were conducted across March, June, and August 2023. Composite filter samples were collected during each flight and analyzed with offline aerosol mass spectrometry (AMS), complemented by on-board real-time sensors and ground-based instrumentation, to provide a comprehensive view of regional aerosol characteristics. Results show clear vertical and seasonal differences in the aerosol composition. Relative to ground-level measurements, aloft samples exhibited shifts in the distribution of organic and inorganic PM, with the organic composition varying distinctly across seasons. Particulate organic nitrogen (ON) was elevated, with bulk compositions similar in March and June but strongly altered in August, likely driven by biomass burning and enhanced photochemical activity. Combined AMS and chemical ionization mass spectrometry analyses detected amines, amides, and amino acids. PM above the planetary boundary layer was enriched in oxidized organic aerosols, while ground-level PM contained higher nitrate and sulfate. Seasonal differences in aqueous-phase processing were also observed, which were strongest in March during persistent cloud cover and weaker in the drier August period, suggesting a shift from aqueous- to gas-phase SOA formation. In conclusion, these findings highlight the value of UAS in advancing PM measurements and vertical profiling of aerosol composition.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Sparge Sampling of Molten Salts for Online Monitoring via Laser-Induced Breakdown Spectroscopy

A method was developed to sample molten salts by sparging to generate and transport aerosols to an isolated instrument for compositional analysis by laser-induced breakdown spectroscopy (LIBS). Real-time monitoring of molten salt composition is critical to developing molten salt nuclear reactors, which offer enhanced safety and efficiency. In this article, the sparge sampling method is described and compared with sampling using a Collison nebulizer. The size distribution and transport of aerosols produced from molten eutectic NaNO 3 –KNO 3 salt were compared for multiple gas flow rates (75–1200 mL min –1 ) and transport distances (0.68–2.61 m). Both methods produced aerosols ranging from 0.5 to 5.0 μm determined using a cascade impactor. Aerosols were effectively transported without pre- or trace-heating of gas lines, but transport efficiency was reduced by the formation of agglomerates. Sparge sampling was found to use less sample and less gas than a Collison nebulizer while producing a more concentrated aerosol stream (up to 5 μg L –1 ). The effects of laser energy and delay time on the signal quality of LIBS measurements of these aerosols were also studied. High energy and short delay times were found to enhance signal and repeatability, whereas signal-to-background and signal-to-noise ratios were highest at low energy and longer delay times. The capabilities of this system for online monitoring of molten salts were demonstrated with calibrations for Sr and Li with relative standard deviations of 2.6% and 1.5% and limits of detection of 380 and 180 μg g –1 , respectively.

Aerosols↗
DOE OSTI2025

Molecular Insights into Geochemical Reactions of Iron-Bearing Minerals: Implications for Hydrogen Geo-Storage

Here, this study investigates the reaction of hydrogen (H 2 ) with pyrite (FeS 2 ), focusing on how temperature and the presence of water influence the reaction pathways and kinetics. Utilizing computational molecular simulations and kinetic analyses, we explore the impact of these factors on the formation of hydrogen sulfide (H 2 S) and related species. First, grand canonical Monte Carlo/molecular dynamics (GCMC/MD) simulations reveal that physical H 2 adsorption occurs in distinct layers on the pyrite surface. In addition, increased temperatures reduce the absolute adsorption capacities. Reactive MD simulations demonstrate that H 2 interacts differently with pyrite under varying conditions. At 298 K, H 2 reacts with pyrite to form HS – , leading to the formation of HS – through covalent bonding with sulfur of pyrite. However, no H 2 S is produced at this temperature, suggesting that a kinetic barrier (i.e., activation energy) may prevent this reaction. At higher temperatures, H 2 S production significantly increases. The presence of water introduces additional complexity to the reaction mechanism. Unlike dry conditions, water enhances H 2 S generation, even at low temperatures. Water also facilitates the formation of additional products, such as SOH, indicating a more intricate chemical environment on the pyrite surface. Our findings identify the association of HS – ions to form H 2 S as the rate-limiting step, with temperature influencing this process. This finding suggests that while the presence of water can create a more dynamic reaction environment, the overall mechanisms leading to H 2 S formation remain consistent. These outcomes suggest the need for developing targeted strategies to manage and control H 2 S emissions within the context of underground hydrogen storage.

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

Basin-Size Mapping: Prediction of Metastable Polymorph Synthesizability Across TaC–TaN Alloys

The sizes of the basins of attraction on the potential energy surface are helpful indicators in determining the experimental synthesizability of metastable phases. In principle, these basins can be controlled with changes in thermodynamic conditions such as composition, pressure, and surface energy. Herein, we use random structure sampling to computationally study how alloying smoothly perturbs basin of attraction sizes. The TaC 1-x N x pseudobinary is an ideal test system given the structural and polymorphic contrast of its parent compounds and their technological relevance as epitaxial substrates for Al 1-x Ga x N. While we find limited thermodynamic stability across all computationally observed phases, random structure sampling shows a significant composition region where the rocksalt basin dominates. As such, we predict the potential for the nonequilibrium synthesis of metastable rocksalt TaC 1-x N x alloys as substrates for Al 1-x Ga x N. At higher nitrogen concentrations, other low-energy metastable polymorphs emerge that continue to retain the hexagonal close packing suitable for III-N growth. Confidence in these trends was established through uncertainty quantification of the basin sizes and energy distributions; such analysis utilized the Beta and Dirichlet distributions. In conclusion, we also find (a) polymorph basin sizes can be rationalized in terms of energetic preferences for different coordination environments; and (b) basin sizes universally shrink with increasing nitrogen content, making the system more prone to amorphous growth.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Methods for Carbon Mass Closure in Polyolefin Hydrocracking

Heterogeneous catalytic hydrocracking of polyolefins is a promising approach for the processing of postconsumer plastics, but product quantification methods remain inconsistent across the literature. In systems that generate a large fraction of vapor-phase products, typical product capture methods can result in large carbon balance deficits, exceeding 50%, compromising reported yields and selectivities. Here, we identify the major sources of product loss and develop enhanced capture methods to improve the quantification accuracy. Seven supplemental techniques were evaluated, targeting either increased vapor recovery (by increasing the volatility or system volume) or enhanced retention in the liquid phase (by decreasing volatility). Among these, a flow collection approach using a continuous helium sweep and downstream gas sampling bag capture yielded the highest recovery, achieving a 96 ± 9.2% carbon balance closure. We show that the efficacy of these methods is strongly dependent on product distribution. In general, solvent addition was most effective when condensable species dominate the product distribution, while flow collection was preferred when both condensable species and light gases are present in high concentrations. These results highlight the need for method-specific workup strategies and demonstrate that no single protocol is universally optimal. We provide general guidelines for selecting and implementing robust product capture techniques, enabling accurate yield and selectivity determinations in polyolefin hydrocracking systems.

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

Predicting Flow in Fracture Networks With Quantum Algorithms

Uncertainty quantification plays a crucial role in the modeling of subsurface flow. For instance, uncertainties in the properties of geologic fracture networks significantly impact flow, requiring numerous simulations to accurately estimate quantities of interest. However, each simulation is computationally expensive because it requires solving a large linear system to capture features that involve both small and large fractures. An example is in percolation, where the interaction of many small fractures (which cumulatively can have a large surface area) with the rock matrix must be modeled precisely. Quantum computing is an emerging tool with the potential to address this issue. Quantum algorithms offer a significant speedup in solving linear systems, achieving efficiencies that are challenging to match with classical approaches. These classical approaches include direct solvers, such as LU decomposition, and iterative methods, notably preconditioned conjugate gradient, commonly used in subsurface modeling to solve large sparse systems. However, applying quantum algorithms to geologic fracture flow requires careful attention to algorithmic and problem-specific constraints to fully realize this quantum advantage. In this work we describe a quantum algorithm for generalized Monte Carlo applications with a quadratic speedup over the classical approaches which can be combined with the quantum speedup, currently under investigation, for solving quantum linear systems for subsurface flow. We show that for quantum algorithms the computational cost of estimating a quantity of interest for a statistical ensemble of networks is roughly the same as that of a single realization, essentially implying that one can get uncertainty quantification for free.

58 GEOSCIENCES↗
DOE OSTI2024

Probing phonon focusing, thermomechanical behavior, and moiré patterns in van der Waals architectures using surface acoustic waves

Surface acoustic waves (SAWs) propagate along solid-air, solid-liquid, and solid-solid interfaces. Their characteristics depend on the elastic properties of the solid. Combining transmission electron microscopy (TEM) experiments with molecular dynamics (MD) simulations, we probe atomic environments around intrinsic defects that generate SAWs in vertically stacked two-dimensional (2D) bilayers of MoS 2 . Our joint experimental-simulation study provides insights into SAW-induced structural and dynamical changes and thermomechanical responses of MoS 2 bilayers. Using MD simulations, we compute mechanical properties from the SAW velocity and thermal conductivity from thermal diffusion of SAWs. The results for Young’s modulus and thermal conductivity of an MoS 2 monolayer are in good agreement with experiments. The presence of defects, such as nanopores which generate SAWs, reduces the thermal conductivity of 2D-MoS 2 by an order of magnitude. We also observe dramatic changes in moiré patterns, phonon focusing, and cuspidal structures on 2D-MoS 2 layers.

36 MATERIALS SCIENCE↗
DOE OSTI2026

High-fidelity entanglement and coherent multi-qubit mapping in an atom array

Neutral atoms in optical tweezer arrays possess broad applicability for quantum technologies, such as computing, analogue simulation and metrology. The atomic species ytterbium-171 is able to host multiple types of qubits, making it a strong candidate for bridging various applications of tweezer arrays. Realizing this potential requires high-fidelity generation and transfer of many-body entanglement between these distinct qubit degrees of freedom. Here, in this study, we demonstrate the creation and coherent mapping of entangled quantum states across multiple qubits in ytterbium-171 tweezer arrays. We map entangled states onto the optical clock qubit from the nuclear spin qubit or the Rydberg qubit. We coherently transfer Z 2 -ordered Greenberger–Horne–Zeilinger states of up to 20 atoms from the interacting Rydberg manifold to the metastable nuclear spin manifold. Furthermore, we find that clock-qubit-based spin detection, when applied to Rydberg and nuclear spin qubits, facilitates atom-loss-detectable qubit measurements and Rydberg decay detection. This enables delayed-erasure detection, yielding an error-detected two-qubit gate fidelity of 99.78(4)% in metastable qubits. These results establish a versatile architecture that advances multiple fields of quantum information science while also establishing bridges between them.

Senoo, Aruku [University of Colorado, Boulder, CO ↗
DOE OSTI2025

Continuous operation of a coherent 3,000-qubit system

Neutral atoms are a promising platform for quantum science, enabling advances in areas ranging from quantum simulations and computation to metrology, atomic clocks and quantum networking. Although atom losses typically limit these systems to a pulsed mode, continuous operation could substantially enhance cycle rates, remove bottlenecks in metrology and enable deep-circuit quantum evolution through quantum error correction. Here we demonstrate an experimental architecture for high-rate reloading and continuous operation of a large-scale atom-array system while realizing coherent storage and manipulation of quantum information. Our approach utilizes a series of two optical lattice conveyor belts to transport atom reservoirs into the science region, where atoms are repeatedly extracted into optical tweezers without affecting the coherence of qubits stored nearby. Using a reloading rate of 300,000 atoms in tweezers per second, we create over 30,000 initialized qubits per second, which we leverage to assemble and maintain an array of over 3,000 atoms for more than 2 hours. Furthermore, we demonstrate persistent refilling of the array with atomic qubits in either a spin-polarized or a coherent superposition state while preserving the quantum state of stored qubits. Our results pave the way for the realization of large-scale continuously operated atomic clocks, sensors and fault-tolerant quantum computers.

atomic and molecular physics↗
DOE OSTI2026

Molecular concentration field design using closed-form steady-state solutions

Control over spatial concentration fields represents a fundamental challenge in designing synthetic biological systems and programmable soft materials. While nature creates morphogen gradients that orchestrate complex developmental processes, synthetic approaches have largely relied on empirical optimization and computationally intensive simulations. Here, we present an analytical framework for steady-state concentration fields generated by finite-sized localized sources in diffusion–degradation systems and derive closed-form solutions for one-, two-, and three-dimensional geometries. By expressing these solutions in dimensionless form, we show that gradient steepness and spatial structure are organized by the Thiele modulus, which captures the competition between diffusion and degradation length scales. The analysis reveals distinct design regimes: in degradation-dominated systems, gradient shape is governed by exponential decay and becomes dimension-independent, whereas in diffusion-dominated systems, gradient magnitude and extent follow dimension-dependent power-law scaling. Building on these results, we introduce a quantitative design strategy that uses threshold-based criteria to program concentration ranges by tuning physically accessible parameters, most directly the production rate, while holding transport and degradation properties fixed. Comparisons with numerical solutions and reported experimental systems demonstrate consistency with the predicted scaling behavior. Together, this work provides a generalizable and physically transparent framework for designing steady-state concentration fields in synthetic biological and soft matter systems, enabling predictive control of gradient-mediated organization without reliance on extensive numerical optimization.

Kim, Dong Woo [Johns Hopkins University, Baltimore↗
DOE OSTI2024

A nonlinear journey from structural phase transitions to quantum annealing

Motivated by an exact mapping between equilibrium properties of a one-dimensional chain of quantum Ising spins in a transverse field (the transverse field Ising (TFI) model) and a two-dimensional classical array of particles in double-well potentials (the “$\phi$ 4 model”) with weak inter-chain coupling, we explore connections between the driven variants of the two systems. Here, we argue that coupling between the fundamental topological solitary waves in the form of kinks between neighboring chains in the classical $\phi$ 4 system is the analog of the competing effect of the transverse field on spin flips in the quantum TFI model. As an example application, we mimic simplified measurement protocols in a closed quantum model system by studying the classical $\phi$ 4 model subjected to periodic perturbations. This reveals memory/loss of memory and coherence/decoherence regimes, whose quantum analogs are essential in annealing phenomena. In particular, we examine regimes where the topological excitations control the thermal equilibration following perturbations. This paves the way for further explorations of the analogy between lower-dimensional linear quantum and higher-dimensional classical nonlinear systems.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Inference of the neutron down-scatter ratio using neutron time-of-flight data for DT-layered implosions on the OMEGA laser

First analysis of the neutron time of flight (nTOF) data is presented to infer the areal density using the neutron down-scatter ratio (DSR) for OMEGA deuterium–tritium (DT)-layered (cryogenic) implosions. The required very-high-dynamic range (nTOF) signal is constructed using an nTOF detector [Forrest et al., Rev. Sci. Ins. 83, 10D919 (2012)] with multiple time-gated photomultiplier tubes. The nTOF data are analyzed using a forward-fit technique to account for the detector responses and infer the DSR and areal density for DT-layered implosions. The areal densities inferred using this analysis are found to be in agreement with the areal densities measured using the magnetic recoil spectrometer detector that is located in a similar line-of-sight. This work establishes the feasibility of the DSR measurement using nTOF data on OMEGA and provides an additional areal-density measurement that will benefit the assessment of the implosion performance and 3D reconstruction of the imploded core.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

A physics-constrained deep learning treatment of runaway electron dynamics

An adjoint formulation leveraging a physics-informed neural network (PINN) is employed to advance the density moment of a runaway electron (RE) distribution forward in time. A distinguishing feature of this approach is that once the adjoint problem is solved, its solution can be used to project the RE density forward in time for an arbitrary initial momentum space distribution of REs. Furthermore, by employing a PINN, a parametric solution to the adjoint problem can be learned. Thus, once trained, this adjoint-deep learning framework is able to efficiently project the RE density forward in time across various plasma conditions while still including a fully kinetic description of RE dynamics. As an example application, the temporal evolution of the density of primary electrons is studied, with particular emphasis on evaluating the decay of a RE population when below threshold. Predictions from the adjoint-deep learning framework are found to be in good agreement with a traditional relativistic electron Fokker–Planck solver, for several distinct initial conditions, and across an array of physics parameters. Once trained, the PINN thus provides a means of generating RE density time histories with exceptionally low online execution time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Many roads to the seam: How conformational flexibility drives nonadiabatic relaxation in a prototypical tetrapyrrolic chromophore

Large and structurally flexible chromophores pose challenges for in silico modeling of photodeactivation due to the many vibrational modes that can funnel the system toward energy degeneracy. In this work, we examine how the multiple degrees of freedom in biliverdin, a prototypical tetrapyrrolic chromophore, cooperate to drive access to the S 1 /S 0 intersection seam in vacuo. We begin by mapping the ground-state potential energy surface to identify representative biliverdin conformers relevant to photoexcitation. We then use a CASSCF-based framework to map the excited-state landscape and characterize the intersection seam, identifying distinct conical-intersection types. Finally, we employ ab initio multiple spawning to resolve the dynamical pathways by which the system accesses these regions. DFT potential-energy and free-energy mappings indicate that, although several conformers are relevant, the “locked-helix” ZsZsZs conformer predominates in the ground state. The intersection seam comprises numerous geometrically distinct regions characterized by varying degrees and combinations of dihedral torsion, pyramidalization, and bond-length alternation. Yet only select regions lie within energetic reach, and moderate barriers separate them from the S 1 minimum. Nonadiabatic dynamics combined with multivariate analyses show that, despite extensive mode coupling during deactivation that guides the system toward multiple regions of the seam, a single dihedral torsion, together with bond-length alternation, predominantly drives energy degeneracy. This work offers new insight into biliverdin’s intrinsic photochemical response and underscores a general feature of flexible chromophores: many modes may participate during photorelaxation, but only a limited subset ultimately dictates seam accessibility.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗