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 721 records · Page 40

Machine Learning for Adaptive Model Refinement to Bridge Scales

This whitepaper is responsive to focal area (2) 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 comprising of a hierarchy of models (e.g., AI driven model/component/parameterization selection). Here we describe scale-aware ML models for adaptive model refinement that allow us to bridge the spatial and/or temporal scales in simulation models and observation data for capturing and predicting extreme water cycles.

58 GEOSCIENCES↗

Machine learning-enabled model-data integration for predicting subsurface water storage

Subsurface water storage (SWS) is a key variable of the climate system and a storage component for precipitation and radiation anomalies, inducing persistence in the climate system. It plays a critical role in climate-change projections and can mitigate the impacts of climate change on ecosystems. However, because of the difficult accessibility of the underground, hydrologic properties and dynamics of SWS are poorly known. Direct observations of SWS are limited, and accurate incorporation of SWS dynamics into Earth system land models remains challenging. We propose a machine learning-enabled model-data integration framework to improve the SWS prediction at local to conus scales in a changing climate by leveraging all the available observation and simulation resources, as well as to inform the model development and guide the observation collection. The accurate prediction will enable an optimal decision of water management and land use and improve the ecosystem's resilience to the climate change.

Lu, Dan↗

Secure State Estimation with Asynchronous Measurements for Coordinated Cyber Attack Detection in Active Distribution Systems

Coordinated cyber attacks tamper with measurement data to disrupt the situational awareness of active distribution systems. Various sensors report measurements asynchronously at different rates, which introduces challenges during state estimation. In addition, this forces cyber intruders to exert greater effort to compromise multiple communication channels and launch coordinated attacks. Therefore, multi-channel and asynchronous measurements could be harnessed to develop more secure cyber defense strategies. In this paper, a prediction-correction-based multi-rate observer is designed to exploit the value of asynchronous measurements for the detection of coordinated false data injection (FDI) attacks. First, a time-function-dependent prediction-correction strategy is proposed to adjust the sampling interval for each sensor’s measurement. Then, an observer is designed based on the trade-off between estimation error and the optimal period of the most recent sampling instant, with the convergence of estimation error with the maximum permitted sampling interval. Moreover, the conditions for exponential stability are developed using the Lyapunov–Krasovskii functional technique. Next, a coordinated FDI attack detection strategy is developed based on the dual nonlinear minimization problem. The proposed attack detection and secure state estimation strategies are tested on the IEEE 13-node system. Simulation results show that these schemes are effective in enhancing attack detection based on asynchronous measurements or compromised data.

asynchronous measurements↗

Comparison of Electromagnetic Transient and Phasor Dynamic Simulations: Implications for Inverter Dominated Systems

The simulation of very high shares of inverter-based resources in power systems has begun to draw into question the validity of phasor domain tools in capturing relevant dynamics. Electromagnetic transient simulators can capture the dynamics of power electronics with substantially smaller time steps, but are computationally expensive. This work contrasts the results of phasor domain and electromagnetic transient tools for simulations on a validated model of the Maui power system operating at very high inverter-based resource shares with near zero voltage forming devices. The results show that the phasor domain tool predicts optimistic stability with fewer voltage forming elements on the network, and loses computational stability before the electromagnetic transient tool. As the electromagnet transient model is of the entire system, and system-wide discrepancies are observed, this case study of a physical power system highlights the potential need for system-wide detailed modeling during periods of very high shares of inverter-based resources and few voltage forming devices.

electromagnetic transient-domain↗

Benthic Kinetic Energy in an Earth System Model and Observations

The benthic layer of the ocean is closely correlated with every upper layer above it. Close analysis of benthic layer stress can yield an accurate depiction of the state of the ocean, as well as the atmosphere above it. Through the use of the Department of Energy’s Energy Exascale Earth System Model (E3SM), we took a look at how these features changed over time in a simulation spanning the historical period from 1958 to 2018 in the Arctic Ocean north of 60° N. Making use of E3SM’s regional refinement capability to place higher resolution in regions of interest, the E3SM-Arctic configuration has ~10 km resolution in the ocean and sea ice components in this region, coarsening up to ~60 km outside of the Arctic (Veneziani et al. 2022). In particular, we narrowly focused on the benthic kinetic energy levels, and how trends changed over time. To validate the fidelity of the E3SM model, we compared model current speed outputs over a 369-day period in the Bering Strait to actual recorded observations during this time period. Our analysis of the benthic kinetic energy mean, standard deviation, and trend throughout this time period can be used to characterize the effects of environmental impacts in the Arctic Ocean.

54 ENVIRONMENTAL SCIENCES↗

Weather and climate extremes in a changing Arctic

Weather and climate extremes are increasingly occurring in the Arctic. Here, in this Review, we evaluate historical and projected changes in rare Arctic extremes across the atmosphere, cryosphere and ocean and elucidate their driving mechanisms. Clear shifts occur in mean and extreme distributions after ~2000. For instance, pre-2000 to post-2000 observational probabilities of 1.5 standard deviation events increase by 20% for atmospheric heat waves, 76.7% for Atlantic layer warm events, 83.5% for Arctic sea ice loss and 62.9% for Greenland Ice Sheet melt extent — in many cases, low probability, rare extreme events in the early period become the norm in the latter period. These observed changes can be explained using a ‘pushing and triggering’ concept, representing interplay between external forcing and internal variability: long-term warming destabilizes the climate system and ‘pushes’ it to a new state, allowing subsequent variability associated with large-scale atmosphere–ocean–ice interactions and synoptic systems to ‘trigger’ extreme events over different timescales. Ongoing anthropogenic warming is expected to further increase the frequency and magnitude of extremes, such that simulated probabilities of 1.5 standard deviation events increase by 72.6% for atmospheric heat waves, 68.7% for Atlantic layer warm events and 93.3% for Greenland Ice Sheet melt rate between historic (1984–2014) and future (2069–2099) periods under a very high emission scenario. Future research should prioritize the development of physically based metrics, enhance high-resolution observation and modelling capabilities and improve understanding of multiscale Arctic climate drivers.

atmospheric dynamics↗

Data for A Hybrid Biophysical-Machine Learning Framework for Diurnal Surface Energy Flux Estimation Using Proximal Sensing

Thermal infrared-based remote sensing of surface energy fluxes has traditionally relied on high spatial resolution satellite data with revisit frequencies on the order of weeks. In this study, we evaluate a biophysics-based analytical surface energy balance model for predicting latent energy (LE) and sensible heat (H) fluxes using proximal sensing observations. The Surface Temperature Initiated Closure (STIC1.2) model has been extensively validated across a wide range of spatial and temporal scales using various satellite-derived thermal infrared data sets. Here we extend this validation by applying STIC at sub-hourly temporal resolution over multiple growing seasons for four distinct agricultural systems. We further develop and evaluate novel STIC variants that incorporate machine learning (ML) techniques to eliminate the need for surface energy balance observations, specifically net radiation and soil heat flux, thereby enhancing model applicability in data-sparse settings. The integration of a ML component to estimate surface available energy is shown to have strong predictive performance for both LE (R2 = 0.81–0.94) and H (R2 = 0.46–0.72) across all agricultural systems examined here, demonstrating the potential of hybrid biophysical-machine learning approaches for surface energy balance modeling with minimal data requirements. This study concludes with a novel application of explainable machine learning (exML) to diagnose sources of model error. This exML framework attributes residual prediction errors to both model input variables and environmental drivers not explicitly included in the simulation experiments. This approach provides a new pathway for improving model design and integrating previously overlooked yet influential variables into future model iterations.

AI/ML↗

A compartmentalized model of multiphase chemical kinetics

There are significant challenges in predicting multiphase chemical kinetics due to the complex coupling of reaction and mass transport across a phase boundary (i.e., interface). Here, we describe a framework for predicting multiphase kinetics that embeds the elementary kinetic steps of reaction, solvation, and diffusion into a coarse grain spatial description of two phases. The model is constructed to bridge the short-timescale interfacial dynamics observed in molecular simulations with the longer timescales observed in kinetic experiments. A simple set of governing differential equations is derived, which, when solved numerically or analytically, yield accurate predictions of multiphase kinetics in microdroplets. Although the equations are formulated for gas-liquid reactions, the underlying conceptual framework is general and can be applied to transformations in other two-phase systems (solid-liquid, liquid-liquid, etc.).

Chemical kinetics and dynamics↗

Creation of Speed-of-Sound Inject Data for Technical Nuclear Forensics

Technical Nuclear Forensics (TNF) exercises simulate a detonation of a nuclear device within the United States, usually in a city. Speed-of-sound (SOS) phenomenology are the atmospheric overpressure and ground shock mechanical motions that are observed at over-pressure (air-blast and infrasound) and seismic sensors, respectively (Figure 1). Amplitudes of SOS data are related to the explosive yield (Kinney and Graham, 1985; Koper et al., 2002; Bonner et al., 2013ab; Ford et al., 2014, 2021; Templeton et al., 2018; Schnurr et al., 2020). Within the country, the United States Prompt Diagnostics System (USPDS) includes a network of geophysical sensors to capture such SOS signals. During the TNF exercise, the event data will be analyzed by the players who know nothing about the technical details of the source (e.g. location, yield, explosive). Specifically, they will use the Integrated Yield Determination Tool (IYDT) to estimate the yield and height-of-burst or depth-of-burial (HOB/DOB). The IYDT allows the user to measure features on the overpressure and seismic channels, evaluate the consistency of features and jointly estimate yield and HOB/DOB. For these exercises, SOS data are simulated for the location, emplacement and yield of the device and signals are propagated to the observing stations. This document describes the steps undertaken in the simulation, validation, preparation and verification of SOS signals for TNF exercises.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Simulations of ocean deoxygenation in the historical era: insights from forced and coupled models

Ocean deoxygenation due to anthropogenic warming represents a major threat to marine ecosystems and fisheries. Challenges remain in simulating the modern observed changes in the dissolved oxygen (O 2 ). Here, we present an analysis of upper ocean (0-700m) deoxygenation in recent decades from a suite of the Coupled Model Intercomparison Project phase 6 (CMIP6) ocean biogeochemical simulations. The physics and biogeochemical simulations include both ocean-only (the Ocean Model Intercomparison Project Phase 1 and 2, OMIP1 and OMIP2) and coupled Earth system (CMIP6 Historical) configurations. We examine simulated changes in the O 2 inventory and ocean heat content (OHC) over the past 5 decades across models. The models simulate spatially divergent evolution of O 2 trends over the past 5 decades. The trend (multi-model mean and spread) for upper ocean global O 2 inventory for each of the MIP simulations over the past 5 decades is 0.03 ± 0.39×10 14 [mol/decade] for OMIP1, -0.37 ± 0.15×10 14 [mol/decade] for OMIP2, and -1.06 ± 0.68×10 14 [mol/decade] for CMIP6 Historical, respectively. The trend in the upper ocean global O2 inventory for the latest observations based on the World Ocean Database 2018 is -0.98×10 14 [mol/decade], in line with the CMIP6 Historical multi-model mean, though this recent observations-based trend estimate is weaker than previously reported trends. A comparison across ocean-only simulations from OMIP1 and OMIP2 suggests that differences in atmospheric forcing such as surface wind explain the simulated divergence across configurations in O 2 inventory changes. Additionally, a comparison of coupled model simulations from the CMIP6 Historical configuration indicates that differences in background mean states due to differences in spin-up duration and equilibrium states result in substantial differences in the climate change response of O 2 . Finally, we discuss gaps and uncertainties in both ocean biogeochemical simulations and observations and explore possible future coordinated ocean biogeochemistry simulations to fill in gaps and unravel the mechanisms controlling the O 2 changes.

54 ENVIRONMENTAL SCIENCES↗

Molecular dynamics analysis of particle number fluctuations in the mixed phase of a first-order phase transition

Molecular dynamics simulations are performed for a finite nonrelativistic system of particles with Lennard-Jones potential. We study the effect of liquid-gas mixed phase on particle number fluctuations in coordinate subspace. A metastable region of the mixed phase, the so-called nucleation region, is analyzed in terms of a noninteracting cluster model. Large fluctuations due to spinodal decomposition are observed. They arise due to the interplay between the size of the acceptance region and that of the liquid phase. These effects are studied with a simple geometric model. The model results for the scaled variance of particle number distribution are compared with those obtained from the direct molecular dynamic simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modeling the Joint Effects of Vegetation Characteristics and Soil Properties on Ecosystem Dynamics in a Panama Tropical Forest

Abstract In tropical forests, both vegetation characteristics and soil properties are important not only for controlling energy, water, and gas exchanges directly but also determining the competition among species, successional dynamics, forest structure and composition. However, the joint effects of the two factors have received limited attention in Earth system model development. Here we use a vegetation demographic model, the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM), ELM‐FATES, to explore how plant traits and soil properties affect tropical forest growth and composition concurrently. A large ensemble of simulations with perturbed vegetation and soil hydrological parameters is conducted at the Barro Colorado Island, Panama. The simulations are compared against observed carbon, energy, and water fluxes. We find that soil hydrological parameters, particularly the scaling exponent of the soil retention curve ( B sw ), play crucial roles in controlling forest diversity, with higher B sw values (>7) favoring late successional species in competition, and lower B sw values (1 ∼ 7) promoting the coexistence of early and late successional plants. Considering the additional impact of soil properties resolves a systematic bias of FATES in simulating sensible/latent heat partitioning with repercussion on water budget and plant coexistence. A greater fraction of deeper tree roots can help maintain the dry‐season soil moisture and plant gas exchange. As soil properties are as important as vegetation parameters in predicting tropical forest dynamics, more efforts are needed to improve parameterizations of soil functions and belowground processes and their interactions with aboveground vegetation dynamics.

54 ENVIRONMENTAL SCIENCES↗

Spin Dynamics of Triple-𝐐 Magnetic Orderings in a Triangular Lattice: Implications for Multi-𝐐 Orderings in General Two-Dimensional Lattices

Multi-𝐐 magnetic structures on two-dimensional (2D) lattices provide a key route to realizing topological physics in 2D magnetism. A major experimental challenge is to unambiguously confirm their formation by excluding the possibility of topologically trivial multidomain single- or double-𝐐 magnetic orders, which cannot be distinguished using conventional diffraction techniques. Here, we propose that long-wavelength spin dynamics offers a universal diagnostic for triangular lattices: Triple-𝐐 orders that preserve rotational symmetry and single- or double-𝐐 orders that break it exhibit qualitatively distinct anisotropies in their Goldstone-mode velocities, stemming from fundamental differences in their underlying spin configurations. We validate this concept using the metallic triangular-lattice antiferromagnet Co 0.325 ⁢TaS 2 , which hosts both a stripe-type single-𝐐 state and a triple-𝐐 tetrahedral ordering at different temperatures. Using inelastic neutron-scattering and spin dynamics simulations, we first refine the spin Hamiltonian by fitting the paramagnetic excitation spectra, allowing us to develop an unbiased model independent of magnetic ordering. We then show that the observed velocity profiles of the Goldstone modes agree with the high-temperature model’s predictions: markedly anisotropic for the single-𝐐 phase and near isotropic for the triple-𝐐 phase. Importantly, this contrast persists across various exchange parameters, highlighting its model-independent nature and suggesting potential applicability to other 2D lattice systems. Beyond the long-wavelength regime, we present a substantial discrepancy between the measured and simulated magnon spectra exclusively in the triple-𝐐 phase. We attribute this discrepancy to magnon energy renormalization arising from order-of-magnitude-enhanced magnon-magnon interactions in the triple-𝐐 phase, due to its noncollinear configuration. This work provides universal insight into the dynamical properties of topological multi-𝐐 magnetic orderings in 2D lattice structures, offering a broadly applicable diagnostic to distinguishing them from topologically trivial single- or double-𝐐 counterparts. The unequivocal confirmation of the triple-𝐐 structure in Co 0.325 ⁢TaS 2 further establishes it as a prominent material platform for exploring topological spin textures in the genuine 2D limit.

Park, Pyeongjae [Oak Ridge National Laboratory (OR↗

Visually resolving the direct Z-scheme heterojunction in CdS@ZnIn 2 S 4 hollow cubes for photocatalytic evolution of H 2 and H 2 O 2 from pure water

We report the direct Z-scheme heterojunction has been recently emerging as an appealing architecture for photocatalysts design. Its efficiency depends on the interfacial and structural features of the photocatalysts. Herein, the two-dimensional ZnIn 2 S 4 nanosheets are grown on the surface of CdS hollow cubes to construct the CdS@ZnIn 2 S 4 hierarchical hollow photocatalysts with chemically bonded interface. The visualized measurements based on spatial-resolved surface photovoltage spectroscopy, combined with other spectroscopic and simulation investigations, clearly disclose that the CdS@ZnIn 2 S 4 hollow cubes constitute a highly efficient direct Z-scheme system. This accounts for the stoichiometric generation of H 2 and H 2 O 2 from pure water observed for the CdS@ZnIn 2 S 4 sulfide-only photocatalysts under visible light irradiation with an apparent quantum efficiency of 1.63 % at 400 nm. The present work demonstrates an effective protocol to achieve comprehensive insights into the charge transfer route at semiconductor heterojunction, and offers a viable way for constructing efficient sulfide-only photocatalysts for driving water splitting reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

First-principle investigation on catalytic hydrogenation of benzaldehyde over Pt-group metals

We report that understanding the hydrogenation of organic compounds in the aqueous phase has always been fundamentally important for improving carbon neutral pathways to fuels and value-added chemicals. In this study, we investigated both thermodynamic and kinetic profiles of benzaldehyde hydrogenation over the Pd(111) and Pt(111) metal surfaces using density functional theory (DFT) and ab initio molecular dynamic (AIMD) simulations. The adsorption of H 2 shows the mixed preference of H adsorption sites on the Pt(111), while the fcc adsorption site is dominant for H on the Pd(111). When benzaldehyde is added to the systems, we observe a strong reduction of benzaldehyde on charged Pd (111) surface compared with that on neutral surface. In contrast, charged state of the Pt(111) surface does not change their interaction. Subsequent hydrogenation reaction of benzaldehyde over Pd(111), proceeding via Langmuir-Hinshelwood mechanism, is affected by two major factors: the presence of H 2 O solvent and surface charge. The presence of H 2 O solvent greatly reduces the activation energy of C-H and O-H bond formation during the hydrogenation process. Furthermore, the hydrogenation step via C-H bond formation is preferred thermodynamically and kinetically over O-H bond formation during thermocatalytic hydrogenation, while the opposite trend holds true during electrocatalytic hydrogenation.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

The Imprint of Southern Ocean Storms on Modeled Surface Chlorophyll, Their Drivers and Satellite Biases

Southern Ocean (SO) phytoplankton chlorophyll is highly variable on sub-seasonal time scales. Although the SO is the windiest ocean basin globally, it is not conclusively understood how storms impact SO phytoplankton dynamics. Much of our existing knowledge stems from satellites, but biases due to data gaps from cloud cover and low solar angles remain unquantified. Here, we use ocean–sea-ice simulations with the Community Earth System Model to quantify the climatological 1997–2018 imprint of storms on chlorophyll and phytoplankton dynamics in the ice-free SO. Additionally, by comparing the full-field model output to synthetic satellite observations, we quantify sampling biases in satellite-derived estimates. We find that both the sign and the magnitude of the average surface chlorophyll imprint vary substantially across storms but last for at least 4 days after the storm passing. Based on our analysis, more than one third of the storms explain the majority of local non-seasonal chlorophyll variability, but satellite-derived storm imprints are often too large in magnitude. On the day of the storm passing, changes in vertical mixing predominantly cause surface chlorophyll anomalies, and reduced light availability due to enhanced cloud cover outweighs the enhanced nutrient availability due to entrainment. Interestingly, storms imprint differently on total net primary production than on surface chlorophyll, demonstrating the difficulty to derive carbon-cycle impacts from a surface-chlorophyll assessment. With SO future storm activity projected to increase, complementing satellite observations with other observing technologies, for example, profiling floats, is necessary to better constrain how storms impact biological carbon cycling in the SO.

Nissen, Cara [University of Colorado, Boulder, CO ↗

Mass Changes the Diffusion Coefficient of Particles with Ligand-Receptor Contacts in the Overdamped Limit

Inertia does not generally affect the long-time diffusion of passive overdamped particles in fluids. Yet a model starting from the Langevin equation predicts a surprising property of particles coated with ligands that bind reversibly to surface receptors: heavy particles diffuse more slowly than light ones of the same size. We show this by simulation and by deriving an analytic formula for the mass-dependent diffusion coefficient in the overdamped limit. We estimate the magnitude of this effect for a range of biophysical ligand-receptor systems, and find it is potentially observable for tailored micronscale DNA-coated colloids.

97 MATHEMATICS AND COMPUTING↗

Generation of polarization entanglement via the quantum Zeno effect

The quantum Zeno effect reveals that continuous observation of a quantum system can significantly alter its evolution. Here, we present a method for establishing polarization entanglement between two initially unentangled photons in coupled waveguides via the quantum Zeno effect. We support our analytical investigation with numerical simulations of the underlying Schrodinger equation describing the system. Further, we extend our technique to three coupled waveguides in a planar configuration and determine the parameters required to generate three-qubit W-states. In contrast to existing schemes based on a vacuum and single-photon encoding, the polarization encoding in our approach is compatible with quantum information protocols that remove photon loss through post-selection. Our findings offer a powerful quantum state engineering approach for photonic quantum information technologies.

Nodurft, Ian C.↗