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 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Local models of two-temperature accretion disc coronae – II. Ion thermal conduction and the absence of disc evaporation

We use local stratified shearing-box simulations with magnetic field-aligned thermal conduction to study an idealized model of the coupling between a cold, radiatively efficient accretion disc, and an overlying, hot, two-temperature corona. Evaporation of a cold disc by conduction from the hot corona has been proposed as a means of mediating the soft-to-hard state transitions observed in X-ray binary systems. We model the coronal plasma in our local disc patch as an MHD fluid subject to both free-streaming ion conduction and a parametrized cooling function that captures the collisional transfer of energy from hot ions to colder, rapidly cooling leptons. In all of our models, independent of the initial net vertical magnetic flux (NF) threading the disc, we find no evidence of disc evaporation. The ion heat flux into the disc is radiated away before conduction can heat the disc’s surface layers. When an initial NF is present, steady-state temperature, density, and outflow velocities in our model coronae are unaffected by conduction. Instead of facilitating disc evaporation, thermal conduction is more likely to feed the disc with plasma condensing out of the corona, particularly in flows without NF. Our work indicates that uncertainties in the amount of NF threading the disc hold far greater influence over whether or not the disc will evaporate into a radiatively inefficient accretion flow compared to thermal conduction. We speculate that a change in net flux mediates disc truncation/evaporation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Shapiro steps and stability of skyrmions interacting with alternating anisotropy under the influence of ac and dc drives

Here we use atomistic simulations to examine the sliding dynamics of a skyrmion in a two-dimensional system containing a periodic one-dimensional stripe pattern of variations between low and high values of the perpendicular magnetic anisotropy. The skyrmion changes in size as it crosses the interface between two anisotropy regions. On applying combined dc and ac driving in either parallel or perpendicular directions, we observe a wide variety of Shapiro steps, Shapiro spikes, and phase-locking phenomena. The phase-locked orbits have two-dimensional dynamics due to the gyrotropic or Magnus dynamics of the skyrmions and are distinct from the phase-locked orbits found for strictly overdamped systems. Along a given Shapiro step when the ac drive is perpendicular to the dc drive, the velocity parallel to the ac drive is locked while the velocity in the perpendicular direction increases with increasing drive to form Shapiro spikes. At the transition between adjacent Shapiro steps, the parallel velocity jumps up to the next step value, and the perpendicular velocity drops. The skyrmion Hall angle shows a series of spikes as a function of increasing dc drive, where the jumps correspond to the transition between different phase-locked steps. At high drives, the Shapiro steps and Shapiro spikes are lost. When both the ac and dc drives are parallel to the stripe periodicity direction, Shapiro steps appear, while if the dc drive is parallel to the stripe periodicity direction and the ac drive is perpendicular to the stripe periodicity, then there are only two locked phases, and the skyrmion motion consists of a combination of sliding along the interfaces between the two anisotropy values and jumping across the interfaces.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Simulating lattice fermion doubling with a Floquet drive

We consider a recently discovered mathematical correspondence between the spectra of a naively discretized lattice fermion and that of a periodically driven (i.e., Floquet) quantum system and enhance it into an infrared equivalence between the two systems. The equivalence can be framed as a duality relation, allowing us to simulate a two-flavor discrete-time fermion theory on the lattice side, where the two flavors arise from time discretization, using a single-flavor fermion theory on the Floquet side. Our demonstration establishes an equivalence between (i) the fermion content, (ii) the correlation functions, and consequently (iii) observables of the two theories in the infrared, going substantially beyond the previously discovered spectral equivalence. We also show how interactions may be incorporated into this enhanced infrared equivalence.

Briceño, Raúl A. [University of California, Berkel↗
DOE OSTI2020

Battery Equivalent Model for Residential HVAC

Flexible loads, especially heating, ventilation, and air-conditioning (HVAC) systems can be used to provide a battery-like service to the power grid by varying their demand up and down over a baseline. Recent work has reported that the round-trip efficiency for providing grid services with HVAC systems can be close to 1. Inspired by this promising result, in this work, we rigorously develop a battery equivalent model (BEM) for residential HVAC systems. Relying only on typically available data, the proposed BEM will bridge the observability and controllability gaps between utility companies and customers. The proposed BEM will enable the utility companies to manage the flexible HVAC devices as easy as charging/discharging the conventional energy storage devices. Validated physics-based models are used to train the BEM, whose prediction accuracy is illustrated through a simulation study.

Dong, Jin↗
DOE OSTI2020

SSA Monte Carlo and Master Equation Modeling of Neutron Leakage Distributions

Strongly stochastic neutron populations in highly subcritical multiplying media containing weak spontaneous fission sources are unique indicators of the presence of SNM. For instance, in neutron multiplicity counting, an extended tail in the count number probability distribution is observed that is a characteristic of and sensitive to the amount of SNM present. Theoretical modeling and efficient computation of stochastic neutron distributions are complementary to experimental measurements and moreover shed light on the fundamental interactions and transport processes that underlie nuclear signatures. Modeling approaches that are intermediate between gross-system level models and high-fidelity but computationally demanding Monte Carlo simulations are of particular interest in this regard. Furthermore we present three advances on existing work in point kinetic or lumped geometry: (i) implementation of the Stochastic Simulation Algorithm (SSA), a system state-updating Monte Carlo method that is more e cient than single-event Monte Carlo in lumped geometries, (ii) application of the backward Master equation to obtain time-gated count moments, and (iii) reconstructing the count number distribution from the deterministically computed moments using a generalized Laguerre polynomial expansion.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

In Situ Measurements of the Reflectances of the LSSTCam Optics and Assessing the Impact of Optical Ghosts

Optical ghosts are image artifacts caused by successive reflections of light between optical surfaces such as lenses, filters, and detectors. These artifacts are unavoidable due to the nonzero reflectances of optical elements and are a major source of contamination for low-surface-brightness science. We use optical ray tracing simulations tuned to observations from LSST Commissioning to quantify the impact of optical ghosts on the LSST data. In particular, we find that ~0.57% of the LSSTCam focal plane is impacted by optical ghosts when averaged across all bands. We also use data from the Collimated Beam Projector to measure the reflectances of various optical elements, generally confirming estimates of ~2% from the systems engineering throughput predictions.

Pai, Aashay↗
DOE OSTI2022

Decline in Seasonal Snow during a Projected 20-Year Dry Spell

Snowpack loss in midlatitude mountains is ubiquitously projected by Earth system models, though the magnitudes, persistence, and time horizons of decline vary. Using daily downscaled hydroclimate and snow projections, we examine changes in snow seasonality across the U.S. Pacific Southwest region during a simulated severe 20-year dry spell in the 21st century (2051–2070) developed as part of the 4th California Climate Change Assessment to provide a “stress test” for water resources. Across California’s mountains, substantial declines (30–100% loss) in median peak annual snow water equivalent accompany changes in snow seasonality throughout the region compared to the historic period. We find that 80% of historic seasonal snowpacks transition to ephemeral conditions. Subsetting empirical-statistical wildfire projections for California by snow seasonality transition regions indicates a two-to-four-fold increase in the area burned, consistent with recent observations of high elevation wildfires following extended drought conditions. By analyzing six of the major California snow-fed river systems, we demonstrate snowpack reductions and seasonality transitions result in concomitant declines in annual runoff (47–58% of historical values). The negative impacts to statewide water supply reliability by the projected dry spell will likely be magnified by changes in snowpack seasonality and increased wildfire activity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Future Cost Benefits Analysis for Electrified Vehicles from Advances Due to U.S. Department of Energy Targets

The U.S. Department of Energy’s Vehicle Technologies Office (DOE-VTO) supports research and development (R&D), as well as deployment of efficient and sustainable transportation technologies, that will improve energy efficiency and fuel economy and enable America to use less petroleum. To accelerate the creation and adoption of new technologies, DOE-VTO has developed specific targets for a wide range of powertrain technologies (e.g., engine, battery, electric machine, lightweighting, etc.). This paper quantifies the impact of VTO R&D on vehicle energy consumption and cost compared to expected historical improvements across vehicle classes, powertrains, component technologies and timeframes. We have implemented a large scale simulation process to develop and simulate tens of thousands of vehicles on U.S. standard driving cycles using Autonomie, a vehicle simulation tool developed by Argonne National Laboratory. Results demonstrate significant additional reductions in both cost and energy consumption due to the existence of VTO R&D targets compared to predicted historical trends. It is observed that, over time, the fuel consumption of different electrified vehicles is expected to decrease by 40–50% and a reduction of 45–55% for vehicle manufacturing costs owing to significant improvements through various VTO R&D targets.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Electrostatic Interactions Control the Nanostructure of Conjugated Polyelectrolyte–Polymeric Ionic Liquid Blends

Polyelectrolyte complexation offers unique opportunities to compatibilize polymers with very different backbone chemistries and to control the morphology of the resulting blend via electrostatic manipulation. In this study, we demonstrate the ability to formulate homogeneous complexes of a conjugated polyelectrolyte with a polymeric ionic liquid, utilizing the electrostatic attraction among their oppositely charged side chains. Variation of electrostatic parameters, such as counterion concentration or polymer charge fraction, tunes the morphology of these polymer complexes from homogeneously disordered blend to weakly structured microemulsion where the local ordering arises from backbone-immiscibility-induced microphase segregation. Our experimental observations are in qualitative agreement with both field-theoretic simulation and random-phase approximation calculations. Simulated morphology snapshots suggest and experimental evidence also indicates that the microphase-segregated complex likely takes on a cocontinuous microemulsion structure. Our findings show that ionic interactions are an effective pathway to compatibilize polymers at macroscopic length scales while achieving controlled nanostructures in these ionic blends. Such systems have great potential for engineering the nanostructure of polymers to tailor applications such as nanofiltration, catalysis, and energy storage, where local ordering can enhance the physical properties of an otherwise macroscopically homogeneous structure.

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

Development and scaling of MeV x-ray radiography at NIF-ARC

We investigate the generation of MeV x-rays using the Advanced Radiography Capability laser system at the National Ignition Facility using 1, 10, and 38 ps pulse durations, with laser energies reaching up to 4 kJ and using compound parabolic concentrators. Hydrodynamic simulations using up-to-date measurements of the contrast of the ARC laser are conducted, allowing us to employ an electron scaling model that incorporates scale length and pulse duration, which aligns closely with the observed temperature distributions. Comparable x-ray sources, in terms of dose, are generated when using a 10 ps pulse duration with 2.4 kJ at ∼1.9 ± 0.4 × 10 18 W/cm 2 and when using a 38 ps pulse duration with 4 kJ at 9.9 ± 0.4 × 10 17 W/cm 2 , both achieving ∼16 rad in air at 1 m for x-ray energies > 0.5 MeV. Radiographs performed on the laser “line-of-sight” show significant improvements in image quality than radiographs performed at 65 degrees to the laser axis. We verify the radiography performance using Monte Carlo Simulations.

Physics - Plasma physics↗
DOE OSTI2025

Topological Phenomena in Artificial Quantum Materials Revealed by Local Chern Markers

A striking example of frustration in physics is Hofstadter’s butterfly, a fractal structure that emerges from the competition between a crystal’s lattice periodicity and the magnetic length of an applied field. Current methods for predicting the topological invariants associated with Hofstadter’s butterfly are challenging or impossible to apply to a range of materials, including those that are disordered or lack a bulk spectral gap. Here, in this work, we demonstrate a framework for predicting a material’s local Chern markers using its position-space description and validate it against experimental observations of quantum transport in artificial graphene in a semiconductor heterostructure, inherently accounting for fabrication disorder strong enough to close the bulk spectral gap. By resolving local changes in the system’s topology, we reveal the topological origins of antidot-localized states that appear in artificial graphene in the presence of a magnetic field. Moreover, we show the breadth of this framework by simulating how Hofstadter’s butterfly emerges from an initially unpatterned 2D electron gas as the system’s potential strength is increased and predict that artificial graphene becomes a topological insulator at the critical magnetic field. Overall, we anticipate that a position-space approach to determine a material’s Chern invariant without requiring prior knowledge of its occupied states or bulk spectral gaps will enable a broad array of fundamental inquiries and provide a novel route to material discovery, especially in metallic, aperiodic, and disordered systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Impact of QCD energy evolution on observables in heavy-ion collisions

Here, we study how the inclusion of energy dependence as dictated by quantum chromodynamic (QCD) small-x evolution equations affects key observables in ultra-relativistic heavy-ion collisions. Specifically, we incorporate JIMWLK evolution into the IP-Glasma framework, which serves as the initial condition for a simulation pipeline that includes viscous relativistic hydrodynamics and a hadronic afterburner. This approach enables a consistent modeling of highly energetic nuclei across varying Bjorken-x values, which are relevant for different collision energies and rapidity regions. In comparison to the standard IP-Glasma setup without small-x evolution, we observe pronounced changes in particle multiplicities and spectral distributions, especially in smaller systems and at the highest available energies. We further explore effects on anisotropic flow observables and correlations between mean transverse momentum and elliptic flow. Our findings underscore the critical role of nonlinear QCD evolution in accurately modeling the early stages of heavy-ion collisions, as well as its implications for extracting transport properties of the quark-gluon plasma.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

A Cyber Threat Mitigation Approach For Wide Area Control of SVCs using Stability Monitoring

We propose a stability monitoring approach for the mitigation of cyber threats directed at the wide area control (WAC) system used for coordinated control of Flexible AC Transmission Systems (FACTS) used for power oscillation damping (POD) of active power flow on inter-area tie lines. The approach involves monitoring the modes of the active power oscillation on an inter-area tie line using the Matrix Pencil (MP) method. We use the stability characteristics of the observed modes as a proxy for the presence of destabilizing cyber threats. We monitor the system modes to determine whether any destabilizing modes appear after the WAC system engages to control the POD. If the WAC signal exacerbates the POD performance, the FACTS falls back to POD using local measurements. The proposed approach does not require an expansive system-wide view of the network. We simulate replay, control integrity, and timing attacks for a test system and present results that demonstrate the performance of the SM approach for mitigation.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Data-Driven Closures and Assimilation for Stiff Multiscale Random Dynamics

Here, we introduce a data-driven and physics-informed framework for propagating uncertainty in stiff, multiscale random ordinary differential equations (RODEs) driven by correlated (colored) noise. Unlike systems subjected to Gaussian white noise, a deterministic equation for the joint probability density function (PDF) of RODE state variables does not exist in closed form. Moreover, such an equation would require as many phase-space variables as there are states in the RODE system. To alleviate this curse of dimensionality, we instead derive exact, albeit unclosed, reduced-order PDF (RoPDF) equations for low-dimensional observables/quantities of interest. The unclosed terms take the form of state-dependent conditional expectations, which are directly estimated from data at sparse observation times. However, for systems exhibiting stiff, multiscale dynamics, data sparsity introduces regression discrepancies that compound during RoPDF evolution. This is overcome by introducing a kinetic-like defect term to the RoPDF equation, which is learned by assimilating in sparse, low-fidelity RoPDF estimates. Two assimilation methods are considered, namely nudging and deep neural networks, which are successfully tested against Monte Carlo simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Addressing GPU memory limitations for Graph Neural Networks in High-Energy Physics applications

Introduction Reconstructing low-level particle tracks in neutrino physics can address some of the most fundamental questions about the universe. However, processing petabytes of raw data using deep learning techniques poses a challenging problem in the field of High Energy Physics (HEP). In the Exa.TrkX Project, an illustrative HEP application, preprocessed simulation data is fed into a state-of-art Graph Neural Network (GNN) model, accelerated by GPUs. However, limited GPU memory often leads to Out-of-Memory (OOM) exceptions during training, due to the large size of models and datasets. This problem is exacerbated when deploying models on High-Performance Computing (HPC) systems designed for large-scale applications. Methods We observe a high workload imbalance issue during GNN model training caused by the irregular sizes of input graph samples in HEP datasets, contributing to OOM exceptions. We aim to scale GNNs on HPC systems, by prioritizing workload balance in graph inputs while maintaining model accuracy. Our paper introduces diverse balancing strategies aimed at decreasing the maximum GPU memory footprint and avoiding the OOM exception, across various datasets. Results Our experiments showcase memory reduction of up to 32.14% compared to the baseline. We also demonstrate the proposed strategies can avoid OOM in application. Additionally, we create a distributed multi-GPU implementation using these samplers to demonstrate the scalability of these techniques on the HEP dataset. Discussion By assessing the performance of these strategies as data loading samplers across multiple datasets, we can gauge their effectiveness in both single-GPU and distributed environments. Our experiments, conducted on datasets of varying sizes and across multiple GPUs, broaden the applicability of our work to various GNN applications that handle input datasets with irregular graph sizes.

Lee, Claire Songhyun↗
DOE OSTI2022

Importance of different parameterization changes for the updated dust cycle modeling in the Community Atmosphere Model (version 6.1)

Abstract. The Community Atmosphere Model (CAM6.1), the atmospheric component of the Community Earth System Model (CESM; version 2.1), simulates the life cycle (emission, transport, and deposition) of mineral dust and its interactions with physio-chemical components to quantify the impacts of dust on climate and the Earth system. The accuracy of such quantifications relies on how well dust-related processes are represented in the model. Here we update the parameterizations for the dust module, including those on the dust emission scheme, the aerosol dry deposition scheme, the size distribution of transported dust, and the treatment of dust particle shape. Multiple simulations were undertaken to evaluate the model performance against diverse observations, and to understand how each update alters the modeled dust cycle and the simulated dust direct radiative effect. The model–observation comparisons suggest that substantially improved model representations of the dust cycle are achieved primarily through the new more physically-based dust emission scheme. In comparison, the other modifications induced small changes to the modeled dust cycle and model–observation comparisons, except the size distribution of dust in the coarse mode, which can be even more influential than that of replacing the dust emission scheme. We highlight which changes introduced here are important for which regions, shedding light on further dust model developments required for more accurately estimating interactions between dust and climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Quantum computing universal thermalization dynamics in a (2 + 1)D Lattice Gauge Theory

Simulating non-equilibrium phenomena in strongly-interacting quantum many-body systems, including thermalization, is a promising application of near-term and future quantum computation. By performing experiments on a digital quantum computer consisting of fully-connected optically-controlled trapped ions, we study the role of entanglement in the thermalization dynamics of a Z 2 lattice gauge theory in 2+1 spacetime dimensions. Using randomized-measurement protocols, we efficiently learn a classical approximation of non-equilibrium states that yields the gap-ratio distribution and the spectral form factor of the entanglement Hamiltonian. These observables exhibit universal early-time signals for quantum chaos, a prerequisite for thermalization. Our work, therefore, establishes quantum computers as robust tools for studying universal features of thermalization in complex many-body systems, including in gauge theories.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Achieving record hot spot energies with large HDC implosions on NIF in HYBRID-E

HYBRID-E is an inertial confinement fusion implosion design that increases energy coupled to the hot spot by increasing the capsule scale in cylindrical hohlraums while operating within the current experimental limits of the National Ignition Facility. HYBRID-E reduces the hohlraum scale at a fixed capsule size compared to previous HYBRID designs, thereby increasing the hohlraum efficiency and energy coupled to the capsule, and uses the cross-beam energy transfer (CBET) to control the implosion symmetry by operating the inner (23° and 30°) and outer (44° and 50°) laser beams at different wavelengths (Δλ> 0). Small case to capsule ratio designs can suffer from insufficient drive at the waist of the hohlraum. We show that only a small amount of wavelength separation between the inner and outer beams (Δλ1-2 Å) is required to control the symmetry in low-gas-filled hohlraums (0.3 mg/cm3 He) with enough drive at the waist of the hohlraum to symmetrically drive capsules 1180 μm in outer radius. This campaign is the first to use the CBET to control the symmetry in 0.3 mg/cm3 He-filled hohlraums, the lowest gas fill density yet fielded with Δλ> 0. We find a stronger sensitivity of hot spot P2 in μm per Angstrom (40–50 μm/Å wavelength separation) than observed in high-gas-filled hohlraums and previous longer pulse designs that used a hohlraum gas fill density of 0.6 mg/cm3. There is currently no indication of transfer roll-off with increasing Δλ, indicating that even longer pulses or larger capsules could be driven using the CBET in cylindrical hohlraums. We show that the radiation flux symmetry is well controlled during the foot of the pulse, and that the entire implosion can be tuned symmetrically in the presence of the CBET in this system, with low levels of laser backscatter out of the hohlraum and low levels of hot electron production from intense laser–plasma interactions. Radiation hydrodynamic simulations can accurately represent the early shock symmetry and be used as a design tool, but cannot predict the late-time radiation flux symmetry during the peak of the pulse, and semi-empirical models are used to design the experiments. Deuterium–tritium (DT)-layered tests of 1100 μm inner radius implosions showed performance close to expectations from simulations at velocities up to ~360 km/s, and record yields at this velocity, when increasing the DT fuel layer thickness to mitigate hydrodynamic mixing of the ablator into the hot spot as a result of defects in the ablator. However, when the implosion velocity was increased, mixing due to these defects impacted performance. The ratio of measured to simulated yield for these experiments was directly correlated with the level of observed mixing. These simulations suggest that reducing the mixing, e.g., by improving the capsule defects, could result in higher performance. In addition, future experiments are planned to reduce the coast time at this scale, delay between the peak compression and the end of the laser, to increase the hot spot convergence and pressure. To reduce the coast time by several hundred ps compared to the 1100 μm inner radius implosions, HYBRID-E has also fielded 1050 μm inner radius capsules, which resulted in higher hot spot pressure and a fusion energy yield of ~170 kJ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗