Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Local time stepping”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Adaptive sampling for accelerating neutron diffraction-based strain mapping *

Abstract Neutron diffraction is a useful technique for mapping residual strains in dense metal objects. The technique works by placing an object in the path of a neutron beam, measuring the diffracted signals and inferring the local lattice strain values from the measurement. In order to map the strains across the entire object, the object is stepped one position at a time in the path of the neutron beam, typically in raster order, and at each position a strain value is estimated. Typical dwell times at neutron diffraction instruments result in an overall measurement that can take several hours to map an object that is several tens of centimeters in each dimension at a resolution of a few millimeters, during which the end users do not have an estimate of the global strain features and are at risk of incomplete information in case of instruments outages. In this paper, we propose an object adaptive sampling strategy to measure the significant points first. We start with a small initial uniform set of measurement points across the object to be mapped, compute the strain in those positions and use a machine learning technique to predict the next position to measure in the object. Specifically, we use a Bayesian optimization based on a Gaussian process regression method to infer the underlying strain field from a sparse set of measurements and predict the next most informative positions to measure based on estimates of the mean and variance in the strain fields estimated from the previously measured points. We demonstrate our real-time measure-infer-predict workflow on additively manufactured steel parts—demonstrating that we can get an accurate strain estimate even with 30%–40% of the typical number of measurements—leading the path to faster strain mapping with useful real-time feedback. We emphasize that the proposed method is general and can be used for fast mapping of other material properties such as phase fractions from time-consuming point-wise neutron measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

On a nonlocal Cahn–Hilliard model permitting sharp interfaces

A nonlocal Cahn–Hilliard model with a non-smooth potential of double-well obstacle type that promotes sharp interfaces in the solution is presented. To capture long-range interactions between particles, a nonlocal Ginzburg–Landau energy functional is defined which recovers the classical (local) model as the extent of nonlocal interactions vanish. In contrast to the local Cahn–Hilliard problem that always leads to diffuse interfaces, the proposed nonlocal model can lead to a strict separation into pure phases of the substance. In this work, the lack of smoothness of the potential is essential to guarantee the aforementioned sharp-interface property. Mathematically, this introduces additional inequality constraints that, in a weak formulation, lead to a coupled system of variational inequalities which at each time instance can be restated as a constrained optimization problem. We prove the well-posedness and regularity of the semi-discrete and continuous in time weak solutions, and derive the conditions under which pure phases are admitted. Moreover, we develop discretizations of the problem based on finite element methods and implicit–explicit time-stepping methods that can be realized efficiently. Finally, we illustrate our theoretical findings through several numerical experiments in one and two spatial dimensions that highlight the differences in features of local and nonlocal solutions and also the sharp interface properties of the nonlocal model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Mineral Deposition on the Rough Walls of a Fracture

Modeling carbonate growth in fractures and pores is important for understanding carbon sequestration in the environment or when supersaturated solutions are injected into rocks. Here, we study the simple but nontrivial problem of calcite growth on fractures with rough walls of the same mineral using kinetic Monte Carlo simulations of attachment and detachment of molecules and scaling approaches. First, we consider wedge-shaped fracture walls whose upper terraces are in the same low-energy planes and show that the valleys are slowly filled by the propagation of parallel monolayer steps in the wedge sides. The growth ceases when the walls reach these low-energy configurations so that a gap between the walls may not be filled. Second, we consider fracture walls with equally separated monolayer steps (vicinal surfaces with roughness below 1 nm) and show that growth by step propagation will eventually clog the fracture gap. In both cases, scaling approaches predict the times to attain the final configurations as a function of the initial geometry and the step-propagation velocity, which is set by the saturation index. The same reasoning applied to a random wall geometry shows that step propagation leads to lateral filling of surface valleys until the wall reaches the low-energy crystalline plane that has the smallest initial density of molecules. Thus, the final configurations of the fracture walls are much more sensitive to the crystallography than to the roughness or the local curvature. The framework developed here may be used to determine those configurations, the times to reach them, and the mass of deposited mineral. Effects of transport limitations are discussed when the fracture gap is significantly narrowed.

calcite↗
DOE OSTI2023

A Review of the Evaluation of Building Energy Code Compliance in the United States

Building energy codes are essential tools for achieving energy efficiency in buildings. However, the full energy savings potential of these codes can only be realized if buildings are constructed in compliance with them. Therefore, evaluating building energy code compliance is crucial in bridging the gap between the energy efficiency requirements set by energy codes and the actualized energy savings achieved. An energy code compliance evaluation serves as a mechanism to assess construction practices, evaluate the effectiveness of code enforcement, identify gaps in compliance, and guide strategies for improvement through training and education. Conducting code compliance evaluation activities involves field studies that require careful design and significant resources. Historically, more emphasis has been placed on developing and adopting building energy codes, while efforts to evaluate compliance have been relatively limited and lacking consistent approaches. The passage of the 2009 American Recovery and Reinvestment Act (ARRA), which mandated that states create plans for achieving 90% compliance within eight years, stimulated the need for an energy code compliance evaluation. As a result, federal, state, and local governments, and utilities have invested in the development of methodologies and tools for code compliance evaluation studies. This paper reviews the code compliance evaluation studies conducted in the United States over the past three decades. It describes and compares the methodologies and metrics used to assess building energy code compliance, summarizes the general elements and steps involved in the evaluation process, and discusses common issues in these studies. Over time, code compliance evaluation methodologies have evolved from isolated development within individual states, regions, and utilities, to widely accepted protocols applicable across different states and local jurisdictions. There has been a transition in compliance metrics, shifting from historical compliance rates to energy-consumption-oriented approaches.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Real-Time Radiological Source Term Estimation for Multiple Sources in Cluttered Environments

A particle filter algorithm is presented to estimate the position, strength, and cardinality of an unknown number of radioactive point sources in an obstacle-rich environment using count measurements. The algorithm addresses gaps in the prior literature by incorporating two novel elements. The first is a precomputation step in which local terrain and obstacle data is processed to compute attenuation kernels throughout the search area. This enables rapid estimation performance in obstacle-rich environments as measurements are gathered. The second novel feature is a dynamic particle allocation technique in which the number of particles is adjusted in real time to meet convergence goals. This feature allows the algorithm to scale more efficiently to scenarios with a larger number of sources. Furthermore, a series of computational experiments using simulated data demonstrates the algorithm’s performance in a cluttered environment with up to eight sources.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Real-Time Radiological Source Term Estimation for Multiple Sources in Cluttered Environments

A particle filter algorithm is presented to estimate the position, strength, and cardinality of an unknown number of radioactive point sources in an obstacle-rich environment using count measurements. The algorithm addresses gaps in the prior literature by incorporating two novel elements. The first is a precomputation step in which local terrain and obstacle data is processed to compute attenuation kernels throughout the search area. This enables rapid estimation performance in obstacle-rich environments as measurements are gathered. The second novel feature is a dynamic particle allocation technique in which the number of particles is adjusted in real time to meet convergence goals. This feature allows the algorithm to scale more efficiently to scenarios with a larger number of sources. A series of computational experiments using simulated data demonstrates the algorithm’s performance in a cluttered environment with up to eight sources.

Kemp, Samuel↗
DOE OSTI2023

Establishing connectivity patterns of eastern oysters ( Crassostrea virginica ) on regional oceanographic scales

Localized depletion of marine and estuarine populations often results from large-scale natural and anthropogenic disturbances (e.g., hurricanes, oil spills) as well as overharvest of fisheries resources. Understanding how such localized depletions may affect populations at larger regional scales requires knowledge of connectivity among local populations within the larger regional landscape or metapopulation. Efforts to restore populations following such decreases require similar knowledge. During the 2010 Deepwater Horizon oil spill, dramatic declines of oysters occurred throughout estuaries west and east of the Mississippi River. We examined trajectories of particles, which were parameterized to mimic oyster larvae, using the ADvanced CIRCulation (ADCIRC; https://adcirc.org ) model to evaluate potential connectivity within and among embayments from Western Louisiana to Alabama. Patterns of larval settlement, which we defined as the intersection of a larval particle with known or expected oyster habitat at any point 13–21 days post-release, reflected much greater local contributions, with 50%–90% of settled particles originating within the same subdivision of embayments. Exchange among subdivision was much less (0%–40%) and settlement originating from outside the embayment of release was trivial under most scenarios (0%–14%). Connectivity between adjacent basins was greatest for larvae released in the southern portions of the embayment, whereas connectivity among nonadjacent basins was not predicted under the scenarios modeled. Because most local populations are relatively isolated on ecological timescales, areas suffering from extensive local depletion are likely to require extensive time to recover due to the lack of larval subsidy from the overall regional population. Restoration would require building stepping-stone populations or reefs within and among basin to restore a high degree of connectivity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Photocurrent in carbon nanotube bundle: Graded Seebeck coefficient phenomenon

Here, polarized photovoltage of a suspended aligned carbon nanotube (CNT) bundle under uniform optical irradiation is discovered without additional structural modification or bias voltage. Such a phenomenon is very surprising considering the metallic behavior of the overall bundle and zero temperature difference between ends. The photovoltage characteristic time is found similar to the thermal response time under step Joule heating and implies a relation to the thermal behavior of the CNT bundle. A similar thermoelectric voltage is also observed during step Joule heating. Localized laser heating and scanning along the axial direction of the bundle uncovers a linear spatial variation of the local Seebeck coefficient. The Seebeck coefficient linearly decreases from root to tip of the CNT bundle with a rate of a few - µV∙K –1 ·mm –1 . Deep investigation in both the microscopic and macroscopic structures of the CNT bundle reveals that the local alignment of CNT assemblies rather than the minor defects in individual CNTs brings about this linear distribution of Seebeck coefficient in space. The finding presents a new way for direct photon-to-electric energy conversion via Seebeck coefficient grading in CNT structures.

42 ENGINEERING↗
DOE OSTI2024

Quantum simulations of hadron dynamics in the Schwinger model using 112 qubits

Hadron wave packets are prepared and time evolved in the Schwinger model using 112 qubits of IBM’s 133-qubit Heron quantum computer ibm_torino. The initialization of the hadron wave packet is performed in two steps. First, the vacuum is prepared across the whole lattice using the recently developed SC-ADAPT-VQE algorithm and workflow. SC-ADAPT-VQE is then extended to the preparation of localized states, and used to establish a hadron wave packet on top of the vacuum. This is done by adaptively constructing low-depth circuits that maximize the overlap with an adiabatically prepared hadron wave packet. Due to the localized nature of the wavepacket, these circuits can be determined on a sequence of small lattices using classical computers, and then robustly scaled to prepare wave packets on large lattices for simulations using quantum computers. Time evolution is implemented with a second-order Trotterization. To reduce both the required qubit connectivity and circuit depth, an approximate quasilocal interaction is introduced. This approximation is made possible by the emergence of confinement at long distances, and converges exponentially with increasing distance of the interactions. Using multiple error-mitigation strategies, up to 14 Trotter steps of time evolution are performed, employing 13,858 two-qubit gates (with a CNOT depth of 370). The propagation of hadrons is clearly identified, with results that compare favorably with Matrix Product State simulations. Finally, prospects for a near-term quantum advantage in simulations of hadron scattering are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Nonoverlapping Localized Exponential Time Differencing Methods for Diffusion Problems

In this paper, we propose nonoverlapping localized exponential time differencing (ETD) methods for diffusion problems. The model time-dependent diffusion equation is first reformulated on subdomains based on the nonoverlapping domain decomposition, in which Neumann boundary conditions are imposed on the interfaces for the subdomain problems and Dirichlet type conditions are enforced to form a space-time interface problem. After spatial discretization by standard central finite differences and temporal integration with the first or second order ETD methods, the fully discrete interface problem is obtained. Such an interface problem is then solved iteratively either at each time step or over the whole time interval: the former involves the solution of stationary problems in each subdomain at each iteration while the latter involves the solution of time-dependent subdomain problems at each iteration. For both approaches, we prove that localized ETD solutions conserve mass exactly and converge in time to the exact space semidiscrete solution. Numerical experiments in two dimensions are also presented to illustrate the performance of the proposed methods.

Mass conservation↗
DOE OSTI2023

Verifying Ray Tracing Amplitude Methods for Global Magnetospheric Modeling

Abstract Ray tracing is a commonly used method for modeling the propagation of electromagnetic waves in Earth's magnetosphere. To apply ray tracing results to global models of wave‐particle interaction such as energetic electron scattering, it is useful to map the discrete rays to a volume filling mesh. However, some methods have inherent losses of energy from the wave source, or do not account for the full range of wave properties within a sample volume. We have developed and tested a 3D magnetospheric ray tracing code “MESHRAY” which resolves these issues. MESHRAY uses the conservation of Poynting flux through ray triplets with finite volume to determine the local field amplitudes. Electromagnetic wave energy density from all ray data points is mapped to a mesh and verified against the wave source power for energy conservation varying time step length, number of rays, and total time steps. We find that the method is self‐consistent and numerically robust. We further investigate whether the neglect of phase information and superposition has a significant impact on the accuracy of mapping wave intensity to a mesh. We find excellent agreement between the analytic solution for waves emitted by a line source in a plane‐stratified medium and an equivalent ray tracing solution. When phase information is excluded, ray tracing reproduces an average amplitude spread over regions of coherent constructive and destructive interference. This may be an important consideration for interpolating ray tracing results of longer wavelength waves such as magnetosonic, electromagnetic ion cyclotron, or ULF waves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Echo state network for coarsening dynamics of charge density waves

An echo state network (ESN) is a type of reservoir computer that uses a recurrent neural network with a sparsely connected hidden layer. Compared with other recurrent neural networks, one great advantage of ESN is the simplicity of its training process. Yet, despite the seemingly restricted learnable parameters, ESN has been shown to successfully capture the spatial-temporal dynamics of complex patterns. Here we build an ESN to model the coarsening dynamics of charge-density waves (CDWs) in a semiclassical Holstein model, which exhibits a checkerboard electron density modulation at half-filling stabilized by a commensurate lattice distortion. The inputs to the ESN are local CDW order parameters in a finite neighborhood centered around a given site, while the output is the predicted CDW order of the center site at the next time step. Special care is taken in the design of couplings between hidden layer and input nodes to ensure lattice symmetries are properly incorporated into the ESN model. Since the model predictions depend only on CDW configurations of a finite domain, the ESN is scalable and transferrable in the sense that a model trained on dataset from a small system can be directly applied to dynamical simulations on larger lattices. Furthermore, our work opens avenues for efficient dynamical modeling of pattern formations in functional electron materials.

2-dimensional systems↗
DOE OSTI2020

Signaling and scrambling with strongly long-range interactions

Strongly long-range interacting quantum systems—those with interactions decaying as a power law 1/$r^α$ in the distance $\textit{r}$ on a $\textit{D}$-dimensional lattice for $\textit{α ≤ D}$—have received significant interest in recent years. They are present in leading experimental platforms for quantum computation and simulation, as well as in theoretical models of quantum-information scrambling and fast entanglement creation. Since no notion of locality is expected in such systems, a general understanding of their dynamics is lacking. In a step towards rectifying this problem, we prove two Lieb-Robinson-type bounds that constrain the time for signaling and scrambling in strongly long-range interacting systems, for which no tight bounds were previously known. In this work, our first bound applies to systems mappable to free-particle Hamiltonians with long-range hopping, and is saturable for $\textit{α ≤ D/2}$. Our second bound pertains to generic long-range interacting spin Hamiltonians and gives a tight lower bound for the signaling time to extensive subsets of the system for all $\textit{α < D}$. This many-site signaling time lower bounds the scrambling time in strongly long-range interacting systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Distributed model predictive control for coordinated, grid-interactive buildings

Continued focus on reducing carbon emissions and improving energy efficiency requires buildings to become grid-interactive and not just behave as static consumers. A distributed model predictive control (DMPC) algorithm known as Limited-Communication (LC) DMPC is modified to enable grid-interactive buildings. A grid-aggregator subsystem is added that allows for a bulk grid power reference signal to be followed while the individual building subsystems also achieve their local comfort objectives. The LC-DMPC algorithm is applied for the first time to systems with multiple buildings. Adequate power tracking is shown for different simulation scenarios involving heterogeneous buildings, and next steps are discussed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

SOLAr Critical Infrastructure Energization (SOLACE): Leveraging Distributed Energy Resources to Provide Local Power

This document is a technical report based on a large-scale DOE-funded project conducted between 2019 and 2022. The project, called SOLAr Critical infrastructure Energization (SOLACE), is aimed at leveraging distributed energy resources (DERs) to provide local power to communities, feeders, or other regions. The project included a particular focus on high-value critical loads such as municipal water supply, telecommunication hubs, and disaster shelters. The key developments of the project were: (1) A comprehensive pre-event power system analysis methodology that identifies and characterizes the viable local power options. (2) A DMS-based example control system for activating and operating the local power solution during a time of crisis. (3) Grid-forming inverter technology to enable isolated local power operation and black-starting. (4) Cyber-security considerations for isolated systems when wider-area communication may be offline. This report describes the overall process that identifies and assesses the viable DER-based local power solutions for a given facility, region, or community. The report documents the process, including the key analysis steps, tools required, data requirements, and recommended pass/fail criteria for each step. It also provides a sample implementation of this process through a test case. Finally, it provides details of how a viable pre-event plan can be selected and activated at go-time when an event has occurred.

14 SOLAR ENERGY↗
DOE OSTI2025

Local chemical ordering of a neutron-irradiated CrFeMnNi compositionally complex alloy

While ion-irradiation studies are a critical first step in studying compositionally complex alloys (CCAs) for nuclear applications, they do not capture all the microstructural changes occurring under the low irradiation dose rates and different particles’ scattering patterns experienced in a nuclear reactor setting. To explore these phenomena in reactor-relevant conditions for the first time in CCA, the single-phase solid-solution Cr 10 Fe 30 Mn 30 Ni 30 was neutron irradiated up to 6.61 displacements per atom at 395 and 579 °C. Irradiation-enhanced local chemical ordering (LCO) well beyond the range of short range ordering was observed, and is predicted to be the precursor to the precipitation of a coherent Ni-Mn L1 0 phase and a Cr-rich α’ phase, though TEM analysis did not indicate the presence of either in any irradiation condition. The line density of faulted dislocation loops decreased from 6.47 to 1.69 ∙ 10 15 m -2 from 3.43 to 6.61 dpa at 579 °C despite no appreciable faulted loop content in the unirradiated material. LCO is expected to increase the complexity of the energy landscape within this alloy, restricting interstitial point defect mobility and creating local regions of greater stacking fault energy. These contribute to the negative correlation between irradiation dose and faulted dislocation loop density in this study, as well as the lack of void swelling observed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Operator Relaxation and the Optimal Depth of Classical Shadows

Classical shadows are a powerful method for learning many properties of quantum states in a sample-efficient manner, by making use of randomized measurements. Here we study the sample complexity of learning the expectation value of Pauli operators via “shallow shadows,” a recently proposed version of classical shadows in which the randomization step is effected by a local unitary circuit of variable depth t. Here we show that the shadow norm (the quantity controlling the sample complexity) is expressed in terms of properties of the Heisenberg time evolution of operators under the randomizing (“twirling”) circuit—namely the evolution of the weight distribution characterizing the number of sites on which an operator acts nontrivially. For spatially contiguous Pauli operators of weight k, this entails a competition between two processes: operator spreading (whereby the support of an operator grows over time, increasing its weight) and operator relaxation (whereby the bulk of the operator develops an equilibrium density of identity operators, decreasing its weight). From this simple picture we derive (i) an upper bound on the shadow norm which, for depth t~log⁡(k), guarantees an exponential gain in sample complexity over the t=0 protocol in any spatial dimension, and (ii) quantitative results in one dimension within a mean-field approximation, including a universal subleading correction to the optimal depth, found to be in excellent agreement with infinite matrix product state numerical simulations. Our Letter connects fundamental ideas in quantum many-body dynamics to applications in quantum information science, and paves the way to highly optimized protocols for learning different properties of quantum states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Quantifying Information without Entropy: Identifying Intermittent Disturbances in Dynamical Systems

A system’s response to disturbances in an internal or external driving signal can be characterized as performing an implicit computation, where the dynamics of the system are a manifestation of its new state holding some memory about those disturbances. Identifying small disturbances in the response signal requires detailed information about the dynamics of the inputs, which can be challenging. This paper presents a new method called the Information Impulse Function (IIF) for detecting and time-localizing small disturbances in system response data. The novelty of IIF is its ability to measure relative information content without using Boltzmann’s equation by modeling signal transmission as a series of dissipative steps. Since a detailed expression of the informational structure in the signal is achieved with IIF, it is ideal for detecting disturbances in the response signal, i.e., the system dynamics. Those findings are based on numerical studies of the topological structure of the dynamics of a nonlinear system due to perturbated driving signals. The IIF is compared to both the Permutation entropy and Shannon entropy to demonstrate its entropy-like relationship with system state and its degree of sensitivity to perturbations in a driving signal.

42 ENGINEERING↗