Engineering Papers⌕ Search

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.

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

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2023

Coarse-grained explicit-solvent molecular dynamics simulations of semidilute unentangled polyelectrolyte solutions

In this study, we present results from explicit-solvent coarse-grained molecular dynamics (MD) simulations of fully charged, salt-free, and unentangled polyelectrolytes in semidilute solutions. The inclusion of a polar solvent in the model allows for a more physical representation of these solutions at concentrations, where the assumptions of a continuum dielectric medium and screened hydrodynamics break down. The collective dynamic structure factor of polyelectrolytes, S(q, t), showed that at q > q*, where q* = 2π/ξ is the polyelectrolyte peak in the structure factor S(q) and ξ is the correlation length, the relaxation time obtained from fits to stretched exponential was $\tau$ KWW ~ q -3 , which describes unscreened Zimm-like dynamics. This is in contrast to implicit-solvent simulations using a Langevin thermostat where $\tau$ KWW ~ q -2 . At q < q*, a crossover region was observed that eventually transitions to another inflection point $\tau$ KWW ~ q -2 at length scales larger than ξ for both implicit- and explicit-solvent simulations. The simulation results were also compared to scaling predictions for correlation length, ξ ~ c$-½\atop{p}$, specific viscosity, η sp ~ c$½\atop{p}$, and diffusion coefficient, D ~ c$0\atop{p}$, where c p is the polyelectrolyte concentration. The scaling prediction for ξ holds; however, deviations from the predictions for η sp and D were observed for systems at higher c p , which are in qualitative agreements with recent experimental results. This study highlights the importance of explicit-solvent effects in molecular dynamics simulations, particularly in semidilute solutions, for a better understanding of polyelectrolyte solution behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Bicrystallography-informed Frenkel–Kontorova model for interlayer dislocations in strained 2D heterostructures

In recent years, van der Waals (vdW) heterostructures and homostructures, which consist of stacks of two-dimensional (2D) materials, have risen to prominence due to their association with exotic quantum phenomena originating from correlated electronic states harbored by them. Atomistic scale relaxation effects play an extremely important role in the electronic scale quantum physics of these systems, providing means of manipulation of these materials and allowing them to be tailored for emergent technologies. We investigate such structural relaxation effects in this work using atomistic and mesoscale models, within the context of twisted bilayer graphene — a well-known heterostructure system that features moiré patterns arising from the lattices of the two graphene layers. For small twist angles, atomic relaxation effects in this system are associated with the natural emergence of interface dislocations or strain solitons, which result from the cyclic nature of the generalized stacking fault energy (GSFE), that measures the interface energy based on the relative movement of the two layers. Here, in this work, we first demonstrate using atomistic simulations that atomic reconstruction in bilayer graphene under a large twist also results from interface dislocations, although the Burgers vectors of such dislocations are considerably smaller than those observed in small-twist systems. To reveal the translational invariance of the heterointerface responsible for the formation of such dislocations, we derive the translational symmetry of the GSFE of a 2D heterostructure using the notions of coincident site lattices (CSLs) and displacement shift complete lattices (DSCLs). The workhorse for this exercise is a recently developed Smith normal form bicrystallography framework. Next, we construct a bicrystallography-informed and frame-invariant Frenkel–Kontorova model, which can predict the formation of strain solitons in arbitrary 2D heterostructures, and apply it to study a heterostrained, large-twist bilayer graphene system. Our mesoscale model is found to produce results consistent with atomistic simulations. We anticipate that the model will be invaluable in predicting structural relaxation and for providing insights into various heterostructure systems, especially in cases where the fundamental unit cell is large and therefore, atomistic simulations are computationally expensive.

2D heterostructures↗
DOE OSTI2024

Piston geometry and stroke optimization for high efficiency propane spark ignition engines

Propane has unique properties and offers interesting characteristics for high-efficiency spark ignition engines. Its high volatility reduces or completely eliminates fuel-wall wetting and facilitates fuel air mixing. Furthermore, propane has a research octane number of 112 and a high octane sensitivity of 15. Finally, its laminar flame speed is on the same order as that of conventional gasoline, and it exhibits high dilution tolerance. Modern spark ignition internal combustion engines rely on fast combustion rates and high dilution to achieve high brake thermal efficiencies. To accomplish this, high stroke-to-bore ratios and high geometric compression ratios have been used in new engine designs. Therefore, propane’s relatively high laminar flame speeds, high knock resistance, and dilution tolerance make it an excellent candidate fuel for modern spark ignition engines. The objective of this work is to co-optimize the piston geometry and the engine stroke to maximize the efficiency of a spark-ignition engine fueled with propane. 3D computational fluid dynamics (CFD) simulations employing the extended coherent flamelet model were used to study the parametric effects of piston shape and stroke length. A piston geometry based on high performing pistons was parameterized using four controlling parameters. The piston geometry and engine stroke design space was explored using deterministic and quasi-random sampling techniques. In conclusion, a Gaussian process regression model was built using the simulation data to explain the results observed.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Temperature response and brine availability to heated boreholes in bedded salt

Abstract There is a growing need for disposal of high‐level nuclear waste. To reduce uncertainty associated with brine availability to repository excavations in salt formations, a collaboration between Sandia, Los Alamos, and Lawrence Berkeley National Laboratories is performing a series of borehole‐scale coupled process tests. Here, we report on the first round of the Brine Availability Test in Salt (BATS) project, a “shakedown” experiment called Phase 1s. Experimental testing included placing a resistive heater, a 260‐W radiative heater, and a 750‐W radiative heater within previously drilled horizontal boreholes at the Waste Isolation Pilot Plant (WIPP) while monitoring temperature and water inflow. The experiments successfully achieved the targeted temperature of 120 °C when using the 750‐W radiative heater. Simulations using FEHM (Finite Element Heat and Mass transfer code) and TOUGH‐FLAC (Transport Of Unsaturated Groundwater and Heat–Fast Lagrangian Analysis of Continua) were able to accurately predict the coupled thermo–hydro–mechanical–chemical response of salt, matching the observed temperature and brine production. Due to the extremely low permeability of salt, these systems take many years to reach steady state when perturbed by mining activities. Long‐term numerical simulations are used to develop the initial pressure and saturation conditions. The inclusion of a damaged rock zone with higher permeability around the borehole also affects the saturation and pressure distributions and plays an important role in dissipating the potential for thermal pressurization. Knowledge gained from this round of experimentation and simulation will be used to conduct the next BATS project experiment in new boreholes at WIPP.

58 GEOSCIENCES↗
DOE OSTI2026

Digital quantum magnetism on a trapped-ion quantum computer

Digital quantum matter—realized when discrete quantum gates approximate continuous time evolution—is susceptible to heating into chaotic, structureless states. If digitization errors are adequately suppressed, a long-lived transient regime of approximately energy-conserving dynamics can be observed on gate-based quantum computers. Conservation of energy, in turn, enables the exploration of a wide variety of complex behaviours observed in equilibrium systems, ranging from the non-trivial microscopic origins of thermalization itself to the stabilization of effective models hosting exotic emergent properties. Here we use Quantinuum’s H2 quantum computer to simulate digitized dynamics of the quantum Ising model, suppressing digitization errors well enough to observe thermalization on timescales that severely challenge classical simulation methods. Relaxation of an inhomogeneous state reveals an emergent hydrodynamics owing to approximate energy conservation and we compute the associated diffusion constant. By reprogramming our simulations to take place on a triangular lattice with periodic boundary conditions, we observe thermalization consistent with emergent gauge and topological constraints resulting from lattice frustration. Furthermore, our results were enabled by continued advances in two-qubit gate quality (native partial entangler fidelities of 99.94(1)%) and establish digital quantum computers as powerful tools for studying (effectively) continuous-time dynamics.

Information theory and computation↗
DOE OSTI2022

Digital quantum simulation of Floquet symmetry-protected topological phases

Quantum many-body systems away from equilibrium host a rich variety of exotic phenomena that are forbidden by equilibrium thermodynamics. A prominent example is that of discrete time crystals, in which time-translational symmetry is spontaneously broken in periodically driven systems. Pioneering experiments have observed signatures of time crystalline phases with trapped ions, solid-state spin systems, ultracold atoms and superconducting qubits. Here we report the observation of a distinct type of non-equilibrium state of matter, Floquet symmetry-protected topological phases, which are implemented through digital quantum simulation with an array of programmable superconducting qubits. We observe robust long-lived temporal correlations and subharmonic temporal response for the edge spins over up to 40 driving cycles using a circuit of depth exceeding 240 and acting on 26 qubits. We demonstrate that the subharmonic response is independent of the initial state, and experimentally map out a phase boundary between the Floquet symmetry-protected topological and thermal phases. Our results establish a versatile digital simulation approach to exploring exotic non-equilibrium phases of matter with current noisy intermediate-scale quantum processors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Application of a Machine Learning Algorithm in Generating an Evapotranspiration Data Product From Coupled Thermal Infrared and Microwave Satellite Observations

Land surface evapotranspiration (ET) is one of the main energy sources for atmospheric dynamics and a critical component of the local, regional, and global water cycles. Consequently, accurate measurement or estimation of ET is one of the most active topics in hydro-climatology research. With massive and spatially distributed observational data sets of land surface properties and environmental conditions being collected from the ground, airborne or space-borne platforms daily over the past few decades, many research teams have started to use big data science to advance the ET estimation methods. The Geostationary satellite Evapotranspiration and Drought (GET-D) product system was developed at the National Oceanic and Atmospheric Administration (NOAA) in 2016 to generate daily ET and drought maps operationally. The primary inputs of the current GET-D system are the thermal infrared (TIR) observations from NOAA GOES satellite series. Because of the cloud contamination to the TIR observations, the spatial coverage of the daily GET-D ET product has been severely impacted. Based on the most recent advances, we have tested a machine learning algorithm to estimate all-weather land surface temperature (LST) from TIR and microwave (MW) combined satellite observations. With the regression tree machine learning approach, we can combine the high accuracy and high spatial resolution of GOES TIR data with the better spatial coverage of passive microwave observations and LST simulations from a land surface model (LSM). The regression tree model combines the three LST data sources for both clear and cloudy days, which enables the GET-D system to derive an all-weather ET product. This paper reports how the all-weather LST and ET are generated in the upgraded GET-D system and provides an evaluation of these LST and ET estimates with ground measurements. The results demonstrate that the regression tree machine learning method is feasible and effective for generating daily ET under all weather conditions with satisfactory accuracy from the big volume of satellite observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Advanced measurement techniques in quantum Monte Carlo: The permutation matrix representation approach

In a typical finite temperature quantum Monte Carlo (QMC) simulation, estimators for simple static observables such as specific heat and magnetization are known. With a great deal of system-specific manual labor, one can sometimes also derive more complicated non-local or even dynamic observable estimators. In contrast, we show that arbitrary static observables can be estimated within the permutation matrix representation (PMR) flavor for any Hamiltonian. We then generalize these results to general imaginary-time correlation functions and non-trivial integrated susceptibilities thereof. Finally, we demonstrate the practical versatility of our method by estimating various non-local, random observables for the transverse-field Ising model on a square lattice and a toy random model.

Permutation matrix representation↗
DOE OSTI2023

Synchrotron X-ray diffraction analysis of constituent phases in transition joint between nickel alloy 738LC and a MnFeCoNiCu alloy

Synchrotron x-ray diffraction (XRD) analysis was performed on transition joints between a single-phase MnFeCoNiCu alloy and Ni-base Alloy 738LC to efficiently identify the constituent phases across the interface, with different levels of material intermixing generated by laser-welding with variable power. Here, the intermixing extent was quantified by postmortem energy dispersive x-ray spectroscopy mapping. Diffraction-based analyses on complex composition spaces with limited prior data present significant challenges because elemental substitution in both disordered and ordered phases is more extensive than in conventional alloy spaces, which may effectuate relatively large shifts in the observed lattice parameters that convolute the analysis. Therefore, thermodynamic simulations and crystallographic literature data were employed to construct a system-specific diffraction library of twelve prospective phases for the composition space investigated. Subsequently, for predicted disordered cubic phases, statistical hard-sphere models were established to estimate the lattice parameters and predict diffraction peak positions for inclusion in the library. The library was then employed to analyze diffraction profiles measured from the variably intermixed transition joints, with focus on accounting for both high and low-intensity peaks. 99.0 % of diffraction peaks with relative intensity greater than 0.001 were assigned to phases from the system-specific library, exemplifying rigorous peak accounting and indicating that no unexpected phases were present. Up to six of the twelve library phases were experimentally found in the transition joints. The lattice parameters predicted by the statistical hard-sphere model based upon thermodynamic simulations agree reasonably well with the measured values for the disordered FCC matrix phase.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Quantum algorithms for disordered physics

We show how a quantum computer may efficiently simulate a disordered Hamiltonian, by incorporating a pseudo-random number generator directly into the time evolution circuit. This technique is applied to quantum simulation of few-body disordered systems in the large volume limit; in particular, Anderson localization. The method requires a number of (error corrected) qubits proportional to the logarithm of the volume of the system, and each time evolution step requires a number of gates polylogarithmic in the volume. We simulate the method to observe the metal-insulator transition on a three-dimensional lattice. Additionally, we demonstrate the algorithm on a one-dimensional lattice, using physical quantum processors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

SOHIP Abel Transform and Onion Peeling Model Module

This software provides tools for analyzing and modeling physical systems using mathematical transforms and layered models. It includes (1) functions for performing the Abel transform, which is used to relate measurements of bending angles to properties such as refractive index and radius in a medium. The code can compute bending angles from input profiles and also reconstruct these profiles from observed data; (2) the functions for modeling systems with multiple layers using an onion-peeling approach, allowing users to simulate and analyze the behavior of layered materials or structures. These capabilities are useful for researchers and engineers working in fields such as optics, atmospheric science, and materials analysis, enabling them to interpret and model data from experiments or simulations relates to refraction in spherical symmetric medium.

Xu, Shuang [Lawrence Livermore National Laboratory↗
DOE OSTI2024

Sensitivity of Arctic Clouds to Ice Microphysical Processes in the NorESM2 Climate Model

Abstract Ice formation remains one of the most poorly represented microphysical processes in climate models. While primary ice production (PIP) parameterizations are known to have a large influence on the modeled cloud properties, the representation of secondary ice production (SIP) is incomplete and its corresponding impact is therefore largely unquantified. Furthermore, ice aggregation is another important process for the total cloud ice budget, which also remains largely unconstrained. In this study, we examine the impact of PIP, SIP, and ice aggregation on Arctic clouds, using the Norwegian Earth System Model, version 2 (NorESM2). Simulations with both prognostic and diagnostic PIP show that heterogeneous freezing alone cannot reproduce the observed cloud ice content. The implementation of missing SIP mechanisms (collisional breakup, drop shattering, and sublimation breakup) in NorESM2 improves the modeled ice properties, while improvements in liquid content occur only in simulations with prognostic PIP. However, results are sensitive to the description of collisional breakup. This mechanism, which dominates SIP in the examined conditions, is very sensitive to the treatment of the sublimation correction factor, a poorly constrained parameter that is included in the utilized parameterization. Finally, variations in ice aggregation treatment can also significantly impact cloud properties, mainly through their impact on collisional breakup efficiency. Overall, enhancement in ice production through the addition of SIP mechanisms and the reduction in ice aggregation (in line with radar observations of shallow Arctic clouds) result in enhanced cloud cover and decreased TOA radiation biases, compared to satellite measurements, especially during the cold months. Significance Statement Arctic clouds remain a large source of uncertainty in projections of the future climate due to the poor representation of the microphysical processes that govern their life cycle. Ice formation is among the least understood processes. While it is widely recognized that better constraints on primary ice production (PIP) are needed to improve existing parameterizations, we show that secondary ice production (SIP) and ice aggregation can have also a significant impact on ice number concentrations. Constraining ice formation through the addition of missing SIP mechanisms and reducing ice aggregation can improve the representation of the cloud macrophysical properties and enhance total cloud cover in the Arctic region, which in turn contributes to decreased TOA radiation biases in the cold months.

Meteorology & Atmospheric Sciences↗
DOE OSTI2021

Representing the Unrepresented Impact of River Ice on Hydrology, Biogeochemistry, Vegetation, and Geomorphology: A Hybrid Physics-Machine Learning Approach

Focal areas include: Predictive modeling through the use of AI techniques and AI-derived model components; the use of AI and other tools to design a prediction system composed of a hierarchy of models (e.g., AI driven model/component/parameterization selection). Insight gleaned from complex data (both observed and simulated) using AI, big data analytics, and other advanced methods, including explainable AI and physics- or knowledge-guided AI.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Development of Laboratory Testbeds Simulating Real-World Greenhouse Gas Emissions

Climate change leads to drastic variations in the environment that can threaten food security. To protect the security and prosperity of the environment, new agricultural practices that mitigate greenhouse gas emissions need to be developed. Agriculture plays a large role in the emissions of greenhouse gases, thus a need for new practices. At present, testing landscape management practices for carbon sequestration and reduction of greenhouse gas emissions are expensive because such tests often require intricate field experiments. Laboratory tests on landscape management practices often fail to simulate and accurately represent the natural systems researchers intend to study. There are methods that can be used to make these laboratory tests more closely replicate natural systems for the improvement of agricultural lands used in biofuel and food production. The goal of this project is to develop a laboratory testbed replicating field conditions that will mimic field greenhouse gas emissions. To achieve this goal, we used soil cores taken from agricultural and marginal lands, then attached fiberglass wicks to improve the drainage system in the soil core to simulate a more natural hydraulic system. The soil cores with the wicks were used to compared to greenhouse gas emissions observed in the field. Our preliminary results suggest that the soil carbon dioxide and nitrous oxide flux from the agricultural land site were comparable to the soil cores equipped with the drainage wicks. The data suggests that laboratory testbeds can be used to mimic field greenhouse gas emissions if parameters, such as soil water content were in a set range. The laboratory testbed produces similar greenhouse gas emissions to the agricultural site when the soil water content was between 15-25%.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Ice‐Nucleating Particles That Impact Clouds and Climate: Observational and Modeling Research Needs

Abstract Atmospheric ice‐nucleating particles (INPs) play a critical role in cloud freezing processes, with important implications for precipitation formation and cloud radiative properties, and thus for weather and climate. Additionally, INP emissions respond to changes in the Earth System and climate, for example, desertification, agricultural practices, and fires, and therefore may introduce climate feedbacks that are still poorly understood. As knowledge of the nature and origins of INPs has advanced, regional and global weather, climate, and Earth system models have increasingly begun to link cloud ice processes to model‐simulated aerosol abundance and types. While these recent advances are exciting, coupling cloud processes to simulated aerosol also makes cloud physics simulations increasingly susceptible to uncertainties in simulation of INPs, which are still poorly constrained by observations. Advancing the predictability of INP abundance with reasonable spatiotemporal resolution will require an increased focus on research that bridges the measurement and modeling communities. This review summarizes the current state of knowledge and identifies critical knowledge gaps from both observational and modeling perspectives. In particular, we emphasize needs in two key areas: (a) observational closure between aerosol and INP quantities and (b) skillful simulation of INPs within existing weather and climate models. We discuss the state of knowledge on various INP particle types and briefly discuss the challenges faced in understanding the cloud impacts of INPs with present‐day models. Finally, we identify priority research directions for both observations and models to improve understanding of INPs and their interactions with the Earth System.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Reactive transport modeling of the Aquifer Thermal Energy Storage (ATES) system at Stockton University, New Jersey during seasonal operations

Hydrogeochemical processes associated with Aquifer Thermal Energy Storage (ATES) operations can often impact the system performance owing to mineral precipitation either at the wellbore or in the aquifer owing to changes in temperature and fluid disequilibria. Although failure of ATES systems due to mineral precipitation ("fouling") is common, predictive reactive-transport models have rarely been applied to plan their design and operation. Here, the objective of this study is to develop a reactive-transport model by coupling thermal, hydrological, and chemical (THC) processes to evaluate effects of introduced atmospheric oxygen on water chemistry, mineral precipitation/dissolution, porosity, and permeability changes associated with an ATES system at Stockton University (New Jersey, USA). The THC model builds on a Thermal-Hydrological-Mechanical (THM) model of the site that evaluated system failure owing to possible fracturing in the caprock or around the wellbore. The causes of the system failure are not known – potential causes include hydraulic fracturing owing to elevated pump pressures that took place, a flow pathway created by one of the boreholes, or a pre-existing natural hydrologic connection between the upper unconfined aquifer and the ATES aquifer, any of which could have led to oxygenated water entering the reservoir and causing the observed Fe-oxide fouling on well screens. The THC model is used to evaluate some of the hypotheses and observations regarding system failure owing to geochemical processes. The reactive-transport code TOUGHREACT V4 was used to model the THC processes during seasonal heating and cooling operations at the Stockton ATES site over 6 years of operation. In the THC simulations, the primary effects on geochemistry were observed when the injection water is saturated with atmospheric oxygen. Simulations show greater precipitation of goethite near the cold wells as compared to the warm wells. Although volume fractions of Fe-hydroxides were relatively small, the model was aimed at processes in the aquifer at the scale of meters and larger rather than at the scale of mm or cm (i.e., a well screen). Kaolinite is the dominant precipitating phase, also around the cold wells. Illite dissolves near the cold wells and precipitates near the warm wells. There is a net decrease in the porosity near the cold wells and increase near the warm wells, although a slight amount of thermal contraction near the cold wells and expansion near the warm wells is responsible for a significant proportion of the porosity change. Owing to the coarse discretization of the numerical grid near the wells (compared to the screen thickness) the magnitude of permeability changes at the wellbore are likely underestimated. The reactive transport model in this study can be used for characterization of aquifers, optimizing the operational parameters (temperature, pressure, pH etc.), and planning of mitigation strategies for ATES systems.

15 GEOTHERMAL ENERGY↗