SEARCH · Engineering Papers

Results for “Observing System Simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Bridging microscopy with molecular dynamics and quantum simulations: an atomAI based pipeline

Recent advances in (scanning) transmission electron microscopy have enabled a routine generation of large volumes of high-veracity structural data on 2D and 3D materials, naturally offering the challenge of using these as starting inputs for atomistic simulations. In this fashion, the theory will address experimentally emerging structures, as opposed to the full range of theoretically possible atomic configurations. However, this challenge is highly nontrivial due to the extreme disparity between intrinsic timescales accessible to modern simulations and microscopy, as well as latencies of microscopy and simulations per se. Addressing this issue requires as a first step bridging the instrumental data flow and physics-based simulation environment, to enable the selection of regions of interest and exploring them using physical simulations. Here we report the development of the machine learning workflow that directly bridges the instrument data stream into Python-based molecular dynamics and density functional theory environments using pre-trained neural networks to convert imaging data to physical descriptors. Additionally, the pathways to ensure structural stability and compensate for the observational biases universally present in the data are identified in the workflow. This approach is used for a graphene system to reconstruct optimized geometry and simulate temperature-dependent dynamics including adsorption of Cr as an ad-atom and graphene healing effects. However, it is universal and can be used for other material systems.

36 MATERIALS SCIENCE↗

DOE OSTI2021

A crystal plasticity model with an atomistically informed description of grain boundary sliding for improved predictions of deformation fields

Crystal plasticity (CP) is a powerful meso-scale technique for deformation modeling in polycrystalline materials. Crystal plasticity models typically do not include an explicit description of grain boundary (GB) sliding, which could lead to inaccurate predictions of strain distributions, especially in the vicinity of GBs. In the present study, a CP model is developed that includes GB sliding as an additional deformation mechanism and models the interaction between slip and GB sliding. Atomistic simulations are used to formulate the constitutive model for GB sliding and its interaction with incident slip. The deformation fields and structure of the GBs are obtained directly from experimental characterization and are faithfully reproduced as bicrystal systems for molecular dynamics simulations. The GBs are modeled as random-type (non-coincidence, asymmetrical) boundaries as observed from experimental data. The underlying atomistic-scale mechanism for pure sliding of random GBs was found to be analogous to fluid-flow. The interaction between incident slip and a sliding GB caused local increases in stress concentration, which further led to a momentary increase in the local sliding rate. The displacement profiles at the sliding GBs computed from the CP model with sliding accommodation are 38% more accurate than the baseline CP model as quantified by the mean squared error between the simulations and experiment. Here these results help improve the sophistication and accuracy of deformation modeling by including physics-based descriptions for GB sliding and its interaction with slip, eventually leading to more reliable predictions of micromechanical quantities.

36 MATERIALS SCIENCE↗

DOE OSTI2024

On the Sensitivity of Coastal Hypoxia to Its External Physical Forcings

Abstract The development of low‐oxygen zones threatening marine life (hypoxia) occurs annually in multiple coastal regions of the world. The largest estuary of the continental United States, the Chesapeake Bay, typically has ≈10 km 3 of water with dioxygen concentrations <3 mg L −1 in July. As numerical methods for refining model resolutions in targeted areas are becoming common, there is interest in assessing the feasibility of simulating coastal hazards such as hypoxia in Earth System Models (ESMs). These coupled models are typically not constrained by observations and thus likely to feature systematic biases in their land, atmosphere, or ocean components. This study relies on four numerical experiments to evaluate the sensitivity of Chesapeake Bay hypoxia to changes or biases in its external physical forcings. Hypoxia exhibits only a minor decrease (−1.6%) after reducing the Bay's terrestrial freshwater discharge by 9.5% (but keeping terrestrial nutrient loadings the same). Changes in freshwater discharges have their largest impact on hypoxia during one extreme event (−37% during 2011 tropical storm Lee). Similarly, changing oceanic conditions on the shelf or their temporal frequency impact hypoxia by only 5%–6%, indicating that the latter is predominantly dictated by local conditions. Although these results are promising from the perspective of ESMs, additional components of ESMs will need to be evaluated before general conclusions can be reached. We notably speculate that the Bay's hypoxia would exhibit higher sensitivity to other forcings not examined here, notably air temperatures and nutrient loadings.

Meteorology & Atmospheric Sciences↗

DOE OSTI2025

Evaluating Mean State Cloud Properties in the Simple Cloud‐Resolving E3SM Atmosphere Model (SCREAM)

Accurately simulating clouds remains a key challenge in global climate models, primarily because cloud formation involves sub-grid processes that are parameterized and crudely represented in models. This study examines the performance of DOE's Simple Cloud-Resolving Energy Exascale Earth System (E3SM) Atmosphere Model (SCREAM) in simulating cloud properties and their spatio-temporal distribution by comparing against satellite observations. Two horizontal resolutions of SCREAM (3 and 12 km) are examined, and both depict a realistic spatial structure of mean-state cloud cover but underestimate its global mean magnitude. SCREAM 3 km reasonably reproduces the distribution of mean-state cloud properties across various cloud optical thickness and cloud-top pressure regimes, with performance comparable to CMIP5 and CMIP6 ensemble and marginally outperforming SCREAM 12 km. Still, SCREAM 3 km tends to underpredict low clouds and optically thin clouds, highlighting the need for continued improvement in representing unresolved processes. This study provides a basis for confidence in the representation of clouds in SCREAM, as simulating mean-state clouds is a necessary prerequisite for trusting its cloud responses to changes in aerosols and greenhouse gases.

Geosciences↗

DOE OSTI2026

Predicting Open Quantum Dynamics with Data-Informed Quantum-Classical Dynamics

We introduce a data-informed quantum-classical dynamics (DIQCD) approach for predicting the evolution of an open quantum system. The equation of motion in DIQCD is a Lindblad equation with a flexible, time-dependent Hamiltonian that can be optimized to fit sparse and noisy data from local observations of an extensive open quantum system. We demonstrate the accuracy and efficiency of DIQCD for both experimental and simulated quantum devices. We show that DIQCD can predict entanglement dynamics of ultracold molecules (calcium fluoride) in optical tweezer arrays. DIQCD also successfully predicts carrier mobility in organic semiconductors (rubrene) with accuracy comparable to nearly exact numerical methods.

Lindblad equation↗

DOE OSTI2024

Random insights into the complexity of two-dimensional tensor network calculations

Projected entangled pair states (PEPS) offer memory-efficient representations of some quantum many-body states that obey an entanglement area law and are the basis for classical simulations of ground states in two-dimensional (2d) condensed matter systems. However, rigorous results show that exactly computing observables from a 2d PEPS state is generically a computationally hard problem. Yet approximation schemes for computing properties of 2d PEPS are regularly used, and empirically seen to succeed, for a large subclass of (“not too entangled”) condensed matter ground states. Adopting the philosophy of random matrix theory, in this work, we analyze the complexity of approximately contracting a 2d random PEPS by exploiting an analytic mapping to an effective replicated statistical mechanics model that permits a controlled analysis at a large bond dimension. Through this statistical-mechanics lens, we argue that (i) although approximately sampling wave-function amplitudes of random PEPS faces a computational-complexity phase transition above a critical bond dimension, and (ii) one can generically efficiently estimate the norm and correlation functions for any finite bond dimension. Furthermore, these results are supported numerically for various bond-dimension regimes. It is an important open question whether the above results for random PEPS apply more generally also to PEPS representing physically relevant ground states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Observing Quantum Measurement Collapse as a Learnability Phase Transition

During a quantum measurement, superpositions of states with different observable properties probabilistically collapse into one with a sharp value of the measured observable. In macroscopic quantum systems, this collapse arises via a continuous measurement-induced phase transition (MIPT) at a critical value of the strength of interaction with the measurement apparatus. MIPTs lie outside established paradigms for equilibrium or nonequilibrium critical phenomena and delineate distinct, stable dynamical and computational phases of matter. Quantum computers enable programmable simulation of the interaction of a measurement apparatus with a dynamical quantum system, to explore MIPT phenomena over a range of system sizes while retaining quantum coherence. Yet, existing experimental protocols rely on fundamentally nonscalable postselection techniques or direct classical simulation of quantum circuits. Here, we report the scalable observation of finite-size scaling evidence for an observable-sharpening MIPT in monitored quantum circuits in a chain of Yb + 171 ions in Quantinuum’s H1-1 trapped-ion quantum processor. By leveraging an equivalent description as a statistical physics problem, we implement scalable classical algorithms to infer the value of the measured observable from a single experimental shot. This technique enables a truly scalable protocol to observe observable-sharpening MIPTs in generic classes of circuits that cannot be directly classically simulated and also provides enhanced means to detect and suppress errors in the quantum simulation. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Simulation of African Easterly Waves in a Global Climate Model

Abstract African easterly waves (AEWs) exert significant influence on local and downstream high-impact weather including tropical cyclone (TC) genesis over the Atlantic. Accurate representation of AEWs in climate and weather prediction models therefore is necessary for skillful predictions. In this study, we examine simulated AEWs, including their evolution, vertical structure, and linkage to tropical cyclone genesis, in the NASA Goddard Earth Observing System Model, version 5 (GEOS-5) atmospheric global climate model. Identified by the leading empirical orthogonal function mode of time-filtered precipitation, the observed westward propagating AEWs along the southern track over the Atlantic are largely captured in GEOS-5, but with a slower phase speed and significantly weaker amplitude downstream off the West Africa coast. The weak downstream development of AEWs in GEOS-5 is accompanied by much reduced TC genesis over the main development region. Further analyses suggest that the slow westward propagation and weaker AEW amplitude downstream can be ascribed to a weak African easterly jet, while overestimated negative (positive) meridional potential vorticity (PV) gradients appear to the north (south) of 10°N in GEOS-5. The greatly overestimated positive meridional PV gradient to the south of 10°N is expected to generate strong horizontal stretching in the AEW wave pattern in the model, which hinders organization of convection and its feedback to sustain the AEW development. Persistent and vigorous AEW precipitation in the Guinea Highlands of the West Africa coast could also be responsible for reduced westward propagation of AEWs in the model.

Meteorology & Atmospheric Sciences↗

DOE OSTI2024

Heat Transport Hysteresis Generated Through Frequency Switching of a Time-Dependent Temperature Gradient

A stochastic energetics framework is applied to examine how periodically shifting the frequency of a time-dependent oscillating temperature gradient affects heat transport in a nanoscale molecular model. We specifically examine the effects that frequency switching, i.e., instantaneously changing the oscillation frequency of the temperature gradient, has on the shape of the heat transport hysteresis curves generated by a particle connected to two thermal baths, each with a temperature that is oscillating in time. Analytical expressions are derived for the energy fluxes in/out of the system and the baths, with excellent agreement observed between the analytical expressions and the results from nonequilibrium molecular dynamics simulations. We find that the shape of the heat transport hysteresis curves can be significantly altered by shifting the frequency between fast and slow oscillation regimes. We also observe the emergence of features in the hysteresis curves such as pinched loops and complex multi-loop patterns due to the frequency shifting. The presented results have implications in the design of thermal neuromorphic devices such as thermal memristors and thermal memcapacitors.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Regime-based aerosol–cloud interactions from CALIPSO-MODIS and the Energy Exascale Earth System Model version 2 (E3SMv2) over the Eastern North Atlantic

This study investigates aerosol-cloud interactions in marine boundary layer (MBL) clouds using an advanced deep-learning-driven synoptic-regime-based framework, combining satellite data (CALIPSO vertically resolved aerosol extinction and MODIS cloud properties) with 1° nudged Energy Exascale Earth System Model version 2 (E3SMv2) simulation over the Eastern North Atlantic (ENA; ∼10°×10°, 2006–2014). The E3SMv2 captures observed seasonal variations in cloud droplet number concentrations (N d ) and liquid water path (LWP), though it systematically underestimates N d . We then partition ENA meteorology into four synoptic regimes (Pre-Trough, Post-Trough, Ridge, Trough) via a deep-learning clustering of ERA5 reanalysis fields, enabling regime-dependent aerosol-cloud interactions analyses. Both satellite and E3SMv2 exhibit an inverted-V LWP-N d relationship. In Post-Trough and Ridge regimes, the satellite shows stronger negative LWP-N d sensitivities than in Pre-Trough regime. The Trough regime displays a muted satellite LWP response. In comparison, the model predicts more exaggerated LWP responses across regimes, with LWP increasing too quickly at low N d and decreasing more sharply at high N d , especially in Pre-Trough and Trough regimes. These exaggerated model LWP sensitivities may stem from uncertainties in representing drizzle processes, entrainment, and turbulent mixing. As for N d susceptibility to aerosols, N d increases with MBL aerosol extinction in both datasets, but the simulated aerosol-cloud interactions appear oversensitive to meteorological conditions. Overall, E3SMv2 better captures aerosol effects under regimes that favor stratiform clouds (Post-Trough, Ridge), but performance deteriorates for regimes with deeper, dynamically complex clouds (Trough), highlighting the need for improved representations of those cloud processes in climate models.

Environmental sciences↗

DOE OSTI2021

Large scale multi-node simulations of $\mathbb{Z}_2$ gauge theory quantum circuits using Google Cloud Platform

Simulating quantum field theories on a quantum computer is one of the most exciting fundamental physics applications of quantum information science. Dynamical time evolution of quantum fields is a challenge that is beyond the capabilities of classical computing, but it can teach us important lessons about the fundamental fabric of space and time. Whether we may answer scientific questions of interest using near-term quantum computing hardware is an open question that requires a detailed simulation study of quantum noise. Here we present a large scale simulation study powered by a multi-node implementation of qsim using the Google Cloud Platform. We additionally employ newly-developed GPU capabilities in qsim and show how Tensor Processing Units -- Application-specific Integrated Circuits (ASICs) specialized for Machine Learning -- may be used to dramatically speed up the simulation of large quantum circuits. We demonstrate the use of high performance cloud computing for simulating $\mathbb{Z}_2$ quantum field theories on system sizes up to 36 qubits. We find this lattice size is not able to simulate our problem and observable combination with sufficient accuracy, implying more challenging observables of interest for this theory are likely beyond the reach of classical computation using exact circuit simulation.

Gustafson, Erik↗

DOE OSTI2021

Laminar and turbulent plasmoid ejection in a laboratory Parker Spiral current sheet

Quasi-periodic plasmoid formation at the tip of magnetic streamer structures is observed to occur in experiments on the Big Red Ball as well as in simulations of these experiments performed with the extended magnetohydrodynamics code, NIMROD. This plasmoid formation is found to occur on a characteristic time scale dependent on pressure gradients and magnetic curvature in both experiment and simulation. Single mode, or laminar, plasmoids exist when the pressure gradient is modest, but give way to turbulent plasmoid ejection when the system drive is higher, which produces plasmoids of many sizes. However, a critical pressure gradient is also observed, below which plasmoids are never formed. A simple heuristic model of this plasmoid formation process is presented and suggested to be a consequence of a dynamic loss of equilibrium in the high- $\beta$ region of the helmet streamer. This model is capable of explaining the periodicity of plasmoids observed in the experiment and simulations, and produces plasmoid periods of 90 minutes when applied to two-dimensional models of solar streamers with a height of $3R_\odot$ . This is consistent with the location and frequency at which periodic plasma blobs have been observed to form by Large Angle and Spectrometric Coronograph and Sun Earth Connection Coronal and Heliospheric Investigation instruments.

Physics↗

DOE OSTI2023

WHOLESCALE: Coordinates of wells at San Emidio, Nevada

This dataset includes position coordinates and elevation information for wells at the WHOLESCALE San Emidio project location. Well positions in the attached file are characterized by UTM coordinates (Easting, Northing) in meters, and WHOLESCALE coordinates (Easting, Northing) relative in meters from a chosen reference point. The elevation of the top and bottom of open intervals within each well (representing locations of perforated or open-hole sections) are measured in meters (positive) above WGS84 geoid (mean sea level). The WHOLESCALE acronym stands for Water & Hole Observations Leverage Effective Stress Calculations and Lessen Expenses. The goal of the WHOLESCALE project is to simulate the spatial distribution and temporal evolution of stress in the geothermal system at San Emidio in Nevada, United States.

15 GEOTHERMAL ENERGY↗

DOE OSTI2026

Functionalized benzylamines from commercial kraft lignin

Benzylamines are key intermediates in pharmaceuticals, agrochemicals, and polymers, but their conventional production relies on benzyl chloride - a petroleum-derived compound with high toxicity and energy demands. Lignin, accounting for up to 30% of plant biomass, is the largest renewable source of aromatic carbon on Earth. However, its highly complex and recalcitrant structure poses a major barrier to efficient conversion into high-value chemicals. Here, in this study, we developed a catalytic approach to convert commercial kraft lignin into phenolic benzylamines through selective depolymerization and subsequent functionalization. We systematically evaluated the effects of three alcohol solvents, formic acid (FA), and a ruthenium-on‑carbon (Ru/C) catalyst on monophenol yield and selectivity. Up to 6.5 wt% monophenol yield was achieved using methanol (MeOH), FA, and Ru/C at 300 °C for 2 h. Quantum thermodynamic simulations based on the COSMO-RS model confirmed the superior solvation and reactivity of the MeOH + FA system, rationalizing observed product yield. The purified monophenolic products, primarily guaiacol and alkyl guaiacols, were then converted into functionalized benzylamines with >90% yield via a multicomponent Mannich reaction under mild conditions. Techno economic analysis (TEA) and life cycle assessment (LCA) underscore the importance of improving lignin depolymerization yields and expanding biorefinery scale. Solvent-only configurations outperform other options in both cost and emissions, with the methanol-only case performing the best ($\$$105 /kg and 26 kg CO 2e /kg) at a large-scale facility. This study establishes a scalable, bio-based pathway for producing benzylamines from commercial kraft lignin, advancing lignin valorization and offering a sustainable alternative to produce petrochemical-based benzylamines.

Catalytic reduction↗

DOE OSTI2024

Responses of summer mesoscale convective systems to irrigation over the North China Plain based on convection-permitting model simulations

Extensive irrigation activities in the North China Plain (NCP) significantly influence regional weather and climate. However, previous studies focusing on the NCP were primarily based on coarse-resolution models, which are unable to explicitly resolve convection systems, causing large uncertainty in precipitation simulations. In this study, a convection-permitting model coupled with a dynamic irrigation scheme is utilized to investigate the impacts of irrigation on summertime mesoscale convective systems (MCSs) over the NCP. Sensitivity experiments with irrigation off and on are conducted for 5 summers and an MCS identification and tracking algorithm is applied to both satellite observations and model simulations. We find that incorporating irrigation in the model increases MCS precipitation, which agrees more with observations. The probability distributions of MCS lifetime, area, propagation speed, and intensity are all better simulated with irrigation. Irrigation increases the occurrence frequency of MCSs throughout the entire day. The nighttime increase is partly because of more frequent local initiation of MCS developed from isolated deep convection, while the daytime increase is mainly attributed to the changes in MCSs initiating elsewhere and then propagating to the NCP. On average, irrigation induces additional moisture that is more thermodynamically favorable for precipitation, but this effect is partially offset by the weakened ascending air motion primarily caused by irrigation surface cooling. Compared to weak MCS precipitation events, strong MCS precipitation events experience greater enhancement in precipitation intensity when including irrigation because the offset effect from the change in large-scale ascending air motion is insignificant. In addition, irrigation makes the variation of MCS precipitation intensity more correlated with the variation in ascending motion but less correlated with that in atmospheric moisture. Our results suggest the pronounced impacts of irrigation on MCSs over the NCP which should be included in numerical models to improve regional precipitation simulation and prediction.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Interaction between cloud–radiation, atmospheric dynamics and thermodynamics based on observational data from GoAmazon 2014/15 and a cloud-resolving model

Observational meteorological data from the field experiment GoAmazon 2014/15 and data from numerical simulations with the cloud-resolving model (CRM) called the System for Atmospheric Modeling (SAM) are used to study the interaction between the cloudiness–radiation as well as the atmospheric dynamics and thermodynamics variables for a site located in the central Amazon region (–3.2°S, –60.6°W) during the wet and dry periods. The main aims are to (a) analyze the temporal series of the integrated cloud fraction, precipitation rate and downward shortwave flux as well as (b) to determine the relationship between the integrated cloud fraction, radiative fluxes and large-scale variable anomalies as a function of the previous day's average. The temporal series of the integrated cloud fraction, precipitation rate and downward shortwave flux from SAM simulations showed physical consistency with the observations from GoAmazon 2014/15. Shallow and deep convection clouds show to have a meaningful impact on radiation fluxes in the Amazon region during wet and dry periods. Anomalies of large-scale variables (relative to the previous day's average) are physically associated with cloud formation, evolution and dissipation. SAM consistently simulated these results, where the cloud fraction vertical profile shows a pattern very close to the observed data (cloud type). Additionally, the integrated cloud fraction and large-scale variable anomalies, as a function of the previous day's average, have a good correlation. These results suggest that the memory of the large-scale dynamics from the previous day can be used to estimate the cloud fraction as well as the water content, which is a variable of the cloud itself. In general, the SAM satisfactorily simulated the interaction between cloud–radiation as well as dynamic and thermodynamic variables of the atmosphere during the periods of this study, being able to obtain atmospheric variables that are impossible to obtain in an observational way.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Uncovering dark matter density profiles in dwarf galaxies with graph neural networks

Dwarf galaxies are small, dark-matter-dominated galaxies, some of which are embedded within the Milky Way. Their lack of baryonic matter (e.g., stars and gas) makes them perfect test beds for probing the properties of dark matter—understanding the spatial dark matter distribution in these systems can be used to constrain microphysical dark matter interactions that influence the formation and evolution of structures in our Universe. Here, we introduce a new method that leverages simulation-based inference and graph-based machine learning in order to infer the dark matter density profiles of dwarf galaxies from observable kinematics of stars gravitationally bound to these systems. Our approach aims to address some of the limitations of established methods based on dynamical Jeans modeling. We show that this novel method can place stronger constraints on dark matter profiles and, consequently, has the potential to weigh in on some of the ongoing puzzles associated with the small-scale structure of dark matter halos, such as the core-cusp discrepancy.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Sachdev-Ye-Kitaev model on a noisy quantum computer

Here we study the SYK model -- an important toy model for quantum gravity on IBM's superconducting qubit quantum computers. By using a graph-coloring algorithm to minimize the number of commuting clusters of terms in the qubitized Hamiltonian, we find the gate complexity of the time evolution using the first-order product formula for N Majorana fermions is $\mathscr{O}$(N 5 J 2 t 2 /ε) where J is the dimensionful coupling parameter, t is the evolution time, and ε is the desired precision. With this improved resource requirement, we perform the time evolution for N=6,8 with maximum two-qubit circuit depth of 343. We perform different error mitigation schemes on the noisy hardware results and find good agreement with the exact diagonalization results on classical computers and noiseless simulators. In particular, we compute return probability after time t and out-of-time order correlators (OTOC) which is a standard observable of quantifying the chaotic nature of quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗