Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time step”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

ETF Grout Expansion: Effects of Different Slags and Struvite Precipitation Protocols (Stirring Time)

The ETF ammonia tolerant grout waste form was developed by the Vitreous State Laboratory (VSL) to stabilize ammonium and thereby prevent emission of ammonia vapor during solidification of an ammonium-rich, concentrated sodium sulfate aqueous waste stream generated at the Hanford Effluent Treatment Facility (ETF). VSL personnel did not observe expansion in samples prepared during the waste form development work. However, Savannah River National Laboratory (SRNL) personnel detected expansion, in a few ETF grout samples, up to ~ 20 percent vertical expansion, while preforming work scope to evaluate the effect of ETF brine compositional ranges on the precipitation and solidification processes. The initial work requests can be found in Washington River Protection Services (WRPS) Statement of Work (SOW, Requisitions #: 339922 Revisions 0 and 1, March 11, 2021, and November 2021, respectively. As a result of the waste form expansion observed by SRNL, additional evaluation was requested by WRPS SOW Requisitions #: 356730 Revision 0 and Revision 1, February 1, 2022, and July 11, 2022, respectively. The goal of these requests was to determine whether differences in the slags and / or in sample preparation protocols used by VSL and SRNL were the cause of the expansion observed in the SRNL samples. G. Chen, WRPS, arranged a materials exchange between VSL and SRNL and coordinated the VSL protocol(s) for grouting and curing 1L batches of Base Case ETF base case simulant. Materials were exchanged between VSL and SRNL. Both laboratories performed the VSL struvite precipitation-solidification protocol and cured samples for at least 28 days. Differences in VSL and SRNL slag compositions, mineralogy and particle size results in themselves or combined were determined to not be the cause of the expansion because all samples expanded regardless of whether the VSL or SRNL slags were used. The fundamental cause of the observed expansion could not be attributed to minor variations in the following: 1) chemistry of the slag and other reagents, 2) container wall thickness, 3) ambient bench top curing conditions at VSL and SRNL, 4) pH adjustments to 7 ± a few tenths of unit, nor 5) stirring time between the precipitation and solidification process steps. Additional observations include: grout samples prepared with VSL slag expanded less than those prepared with the SRNL slag and, based on a very limited sample set, some samples stirred for longer times expanded less than those stirred for 20 minutes.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Direct Study of Changes in Catalyst Structure-Kinetic Properties During Redox Transitions

The Temporal Analysis of Products (TAP) pulse response methodology is a transient technique that provides the time resolution needed to deconvolve reaction steps from the complex networks typical in industrial catalytic processes. Traditionally, TAP measurements observe gas phase dynamics at the reactor exit but lack direct measurements of changes in the catalyst itself. Recently, a new operando technique was developed that couples gas phase transients to dynamic changes in metal centers with the precise TAP methodology for nanomole titration. Using an industrial CrOx/Al2O3 catalyst used for propane dehydrogenation, we demonstrate the capabilities of this unique device to reveal key catalytic processes: 1) total propane oxidation not associated with chromia centers, 2) reduction of Cr6+ to Cr3+ correlated with selective product formation, and 3) subsequent carbon accumulation. By utilizing incremental pulsing in a diffusion-only transport regime, the spectrokinetic device allows us to resolve detailed changes in catalyst structure, composition, and kinetic function that are otherwise indistinguishable in conventional operando devices. The unification of the TAP methodology with time-resolved spectroscopic measurements offers new and unique insights into the complex kinetic phenomena regulated by solid catalyst surfaces.

03 - NATURAL GAS↗
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 OSTI2021

Resilience Metrics and Framework for Distributed Wind Presentation

This presentation communicates information about the MIRACL project Resilience Metrics report and Resilience Framework report. It was created for the 2021 MIRACL advisory board meeting. We propose a three-tiered approach for the resilience framework. At the top level, we consider the time horizons on which resilience will be evaluated and executed. At the middle level, we consider the core functions of resilience, which span across the time horizons. At the lower level, we consider the process steps that correspond to implementing practices for resilience in each of the core functions. The framework considers three time horizons in order to enable organizations to assess and improve their system’s resilience throughout its lifecycle. We call these time horizons the planning, operational, and future stages. The planning stage uses organizational needs and current system evaluation to prepare for potential hazards. The operational stage seeks to execute responses to hazards as prudently and efficiently as possible to maintain system resilience. The future stage seeks to improve on current system resilience and feeds back into the planning stage to promote continuous improvement. While all three time horizons are important when considering a specific topic, the planning and evaluation phase (i.e., what is done in advance of the event) is critical in defining a system’s resilience characteristics and in outlining how a system responds to an event. This framework intentionally emphasizes the planning stage to highlight the overarching emphasis of this effort, not to imply that the other two time-related horizons (i.e., operational and change the future) are less important. The core functions in the framework are identify, prepare, detect, adapt, and recover. These five functions stem from a rigorous analysis of definitions used across the industry, and they represent the core capabilities that an organization must have to enable lifecycle resilience. Within each core function, process steps are described that help walk an organization through the information gathering, evaluation, decision-making, and implementation processes they will need to ensure their resilience goals are maintained throughout the system and the system lifecycle. Also highlighted in the figure is the concept that a resilience framework should be cyclical in nature. Because a system’s resilience is based on finite resources and time, it must continually evolve through this framework’s risk management and capital investment steps at an appropriate level of scope and pace.

17 WIND ENERGY↗
DOE OSTI2021

Toward active disruption avoidance via real-time estimation of the safe operating region and disruption proximity in tokamaks

This paper describes a real-time capable algorithm for identifying the safe operating region around a tokamak operating point. The region is defined by a convex set of linear constraints, from which the distance of a point from a disruptive boundary can be calculated. The disruptivity of points is calculated from an empirical machine learning predictor that generates the likelihood of disruption. While the likelihood generated by such empirical models can be compared to a threshold to trigger a disruption mitigation system, the safe operating region calculation enables active optimization of the operating point to maintain a safe margin from disruptive boundaries. The proposed algorithm is tested using a random forest disruption predictor fit on data from DIII-D. The safe operating region identification algorithm is applied to historical data from DIII-D showing the evolution of disruptive boundaries and the potential impact of optimization of the operating point. Real-time relevant execution times are made possible by parallelizing many of the calculation steps and implementing the algorithm on a graphics processing unit. Lastly, a real-time capable algorithm for optimizing the target operating point within the identified constraints is also proposed and simulated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Transmon Qubit Using Sn as a Junction Superconductor

Superconductor qubits typically use aluminum–aluminum oxide tunnel junctions to provide the nonlinear inductance. Junctions with semiconductor barriers make it possible to vary the superconductor material and explore beyond aluminum. We use InAs semiconductor nanowires coated with thin superconducting shells of β-Sn to realize transmon qubits. By tuning the Josephson energy with a gate voltage, we adjust the qubit frequency over a range of 3 GHz. The longest energy relaxation time, T 1 = 27 μs, is obtained at the lowest qubit frequencies, while the longest echo dephasing time, T 2 = 1.8 μs, is achieved at higher frequencies. Finally, we assess the possible factors limiting coherence times in these devices and discuss steps to enhance performance through improvements in materials fabrication and circuit design.

B-Sn↗
DOE OSTI2021

Development of a cosmic ray oriented trigger for the fluorescence telescope on EUSO-SPB2

The Extreme Universe Space Observatory on a Super Pressure Balloon 2 (EUSO-SPB2), in preparation, aims to make the first observations of Ultra-High Energy Cosmic Rays (UHECRs) from near space using optical techniques. EUSO-SPB2 will prototype instrumentation for future satellite-based missions, including the Probe of Extreme Multi-Messenger Astrophysics (POEMMA) and K-EUSO. The payload will consist of two telescopes. The first is a Cherenkov telescope (CT) being developed to quantify the background for future below-the-limb very high energy (E 10 PeV) astrophysical neutrino observations, and the second is a fluorescence telescope (FT) being developed for detection of UHECRs. The FT will consist of a Schmidt telescope, and a 6192 pixel ultraviolet camera with an integration time of 1.05 s. The first step in the data acquisition process for the FT is a hardware level trigger in order to decide which data to record. In order to maximize the number of UHECR induced extensive air showers (EASs) which can be detected, a novel trigger algorithm has been developed based on the intricacies and limitations of the detector. Finally, the expected performance of the trigger has been characterized by simulations and, pending hardware verification, shows that EUSO-SPB2 is well positioned to attempt the first near-space observation of UHECRs via optical techniques.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Adaptive workflow for simulation of RF heaters

Accurate RF (Radio Frequency) simulations of fusion systems like ITER require the definition of high-fidelity analysis geometries that include detailed antenna, reactor wall, and physics regions. Here, this paper will describe a workflow for the execution of adaptive high-performance simulations of RF fusion systems. In this workflow, the simulation input consists of a CAD model attributed with the needed analysis attributes. The analysis mesh is automatically generated and the analysis steps are executed using the time-harmonic Maxwell's equations solved using high-order Nédélec finite elements. A patch recovery-based error estimator is used to drive a conforming mesh adaptation procedure.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Peridynamics and surrogate modeling of pressure-driven well stimulation

In this work we use the peridynamics theory of solid mechanics to simulate fracture in an annular rock domain subject to an in-situ stress and create surrogate models that predict the area of the resulting cracks. Peridynamics is a non-local formulation of continuum mechanics that naturally accommodates material discontinuities. Furthermore, unlike other fracture modeling techniques there is no need to provide information about the crack path. We utilize the peridynamics code Peridigm and take a two-stage approach to fracture modeling. First an implicit solve is performed to compute the in-situ stress state. We then execute an explicit solve where a pressure loading designed to emulate fluid-driven hydraulic fracture is applied at the borehole and transmitted to the pre-stressed rock. We present results from polynomial and single and multi-level Gaussian process surrogate models constructed from a sampling study of the peridynamics model. The surrogates predict crack area given a measure of the in-situ stress anisotropy and rise time and amplitude of the pressure loading. These surrogates take a minuscule fraction of peridynamics model's running time to evaluate and are a step towards enabling advanced optimization and uncertainty quantification workflows that require many model evaluations.

42 ENGINEERING↗
DOE OSTI2025

Sequential Infiltration Synthesis of Multinary Metal Oxides: The Role of Precursor Diffusion and Reactivity

Understanding volatile metal–organic precursor transport in polymer matrices during sequential infiltration synthesis (SIS) is critical to the design of hybrid materials of tailored composition and uniform structure. Here, we investigate the diffusion and reaction behavior of diethylzinc (DEZ) in poly(methyl methacrylate) (PMMA) matrices with variable indium and zinc incorporation. Unlike conventional SIS where precursors adduct or react directly with the polymer backbone, DEZ shows minimal interaction with PMMA. However, DEZ may interact strongly with In–OH or Zn–OH moieties that are previously incorporated into a PMMA thick film. Using spectroscopic ellipsometry, X-ray photoelectron spectroscopy depth profiling, and EDX line scans, we quantify Zn incorporation as a function of exposure time, polymer thickness, and previous infiltration steps. The Zn distribution evolves from surface-localized to deep within the ∼1 μm thick film with increasing DEZ exposures up to 30 min. A reaction-diffusion model incorporating a hindering factor to capture diffusivity decay due to product accumulation reproduces the Zn depth profile and yields physically meaningful parameters to guide future growth. Furthermore, this study provides a framework for interpreting metal–organic precursor infiltration behavior and highlights the utility of model-guided SIS process design.

diethylzinc↗
DOE OSTI2020

Lightning Interferometry Uncertainty, Beam Steering Interferometry, and Evidence of Lightning Being Ignited by a Cosmic Ray Shower

Abstract We present an uncertainty analysis on correlation‐based time delay estimate, the basis for broadband lightning interferometry. A normal interferometry would yield much higher uncertainties than the theoretically predicted lower bound due to finite correlation window lengths. However, if the signals are aligned in time, the uncertainties approach the lower bound and can be used to indicate the interferometric uncertainties. Based on this, we introduce a beam steering interferometry technique. It first estimates a direction centroid for the lightning sources with a normal interferometry and computes the time delays among the signals. It then shifts the raw data with the time delays to align the signals roughly in time and reprocess the aligned signals. It finally pushes the reprocessed results to their correct positions, based on the time delays estimated in the first step. We apply this technique on a fast positive breakdown (FPB) process that started a normal intracloud lightning. The FPB process is shown to have a much more complex structure and development than a normal interferometry would provide. More importantly, from both interferometric and polarization analyses, the FPB appears to be ignited by a cosmic ray shower (CRS). We estimate the radio frequency strength and frequency content related to a presumed CRS in a thunderstorm electric field and find that they are in agreement with the observations. We examine the electric field effect of the CRS front and find that it could raise the field above the threshold for positive breakdown and is capable of igniting the FPB (and the lightning).

Shao, Xuan‐Min↗
DOE OSTI2022

Efficient Loads Surrogates for Waked Turbines in an Array

Accurately and efficiently predicting wind turbine structural loading is a crucial step in wind farm design. Without considering structural loading, wind farm optimization could negatively impact turbine fatigue and ultimate loads, especially for waked and partially waked turbines, which could result in higher maintenance costs and reduced turbine lifetime. However, predicting turbine loads throughout an array is a costly step, as these quantities require time-accurate results across long time histories, which is often intractable for large array optimization. Therefore, surrogate models that link array spacing to load outputs are often used, but the surrogates are then unique to the inflow conditions and array configurations in the training library. This work develops surrogate models for many wind turbine load outputs based solely on rotor plane velocity measurements, with no required input about array configuration or freestream inflow parameters. Surrogate models were constructed for many turbine quantities of interest (QoI), considering mean, standard deviation, ultimate, and fatigue loads. In general, most QoI statistics were accurately captured, as measured by predicted vs. actual correlation coefficient, confirming the suitability of the approach. Temporal mean values of the QoI required only temporal mean measurements of the rotor plane inflow velocity. However, accurate prediction of temporal standard deviation, ultimate, and fatigue values of QoI also required temporal standard deviations of the rotor plane velocity field. Poor surrogate performance was observed when too many correlated inputs were used, such as multiple velocity components. If the fewest inflow parameters are used to construct the surrogates, the average correlation coefficient value for all output QoI statistics is 0.89. Surrogates for standard deviations and damage equivalent loads (DELs) of turbine QoIs generally had lower accuracy and tower-base and shaft load channels posed the most difficult to capture accurately. The results suggest that these surrogates could be easily paired with analytic wake models, which are frequently used for pre-construction wind farm array optimization, to account for turbine loading in addition to power production. By including the optimal inflow conditions, the surrogate accuracy can improve to an average correlation coefficient value for all output QoI statistics of 0.92. This work has established the ability to build accurate surrogates for mean, standard deviation, ultimate load, and DEL turbine QoI values based on the rotor plane inflow velocity, and identified which inflow conditions lead to greater surrogate accuracy.

17 WIND ENERGY↗
DOE OSTI2023

Collective neutrino oscillations on a quantum computer with hybrid quantum-classical algorithm

We simulate the time evolution of collective neutrino oscillations in two-flavor settings on a quantum computer. We explore the generalization of Trotter-Suzuki approximation to time-dependent Hamiltonian dynamics. The trotterization steps are further optimized using the Cartan decomposition of two-qubit unitary gates U ϵ SU(4) in the minimum number of controlled-NOT (CNOT) gates making the algorithm more resilient to the hardware noise. As a result, a more efficient hybrid quantum-classical algorithm is also explored to solve the problem on noisy intermediate-scale quantum devices.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Mechanistic Insights to Degradation of PERC Minimodules with Differentiated Packaging Materials & Module Architectures

In this paper, we study the degradation behavior of glass-backsheet (GB) and double glass (DG) multicrystalline sil- icon monofacial PERC minimodules, fabricated using polyolefin elastomer (POE) and ethylene vinyl acetate (EVA) encapsulants, undergoing exposures of modified damp-heat (80 °C / 85% RH), with and without full spectrum light. The completed exposure time, at this time, is 1512 hours with each step of 504 hours. Step- wise measurements are conducted throughout exposures including current-voltage (I-V ) curves, Suns-V oc , electroluminescence (EL) images, and four-point proof loading with EL for tracking changes in mechanical properties in encapsulants by cell fracture probability. The results show that GB minimodules of both encapsulants, compared to DG minimodules, experience more changes in features of I-V and Suns-V oc measurements. There is no development of PV cell cracking during proof loading as seen from EL images and the load versus load line displacement plots currently.

14 SOLAR ENERGY↗
DOE OSTI2021

Wind Plant Performance Prediction Benchmark Phase 1 (Technical Report)

Financial risk resulting from the uncertainty associated with developing, owning, and operating wind power plants remains a barrier to reducing the levelized cost of energy (LCOE). On average, modern wind power plants in the U.S. underperform their expected annual energy output by 3.5-4.5% , with many underperforming by over 10%. To compensate for this uncertainty, investors require a larger return on investment (ROI) and apply "knock-down" factors that mask much of the underlying sources of uncertainty. Wind energy projects thus have reduced access to low-cost capital. Furthermore, operating wind plants often take a simple approach to estimating operations & maintenance (O&M) costs (e.g. straight-line estimates based on similar plants), which can eat into profits. To overcome these issues, the wind industry must improve the models they use for estimating wind plant performance and operations. An industry consortium (IC) requested that the National Renewable Energy Laboratory (NREL) lead a Department of Energy (DOE) working group to benchmark the accuracy of wind power plant energy predictions against real operational data. The IC was also motivated by DOE and NREL's potential to characterize systematic energy underperformance, identify sources of uncertainty, and explore root causes. The Wind Plant Performance Prediction (WP3) project was created out of this request, and this report represents the successful completion of Phase 1 of the WP3 project. During the project, wind plant owners provided both pre-construction and operational data to NREL. The pre-construction data was provided to wind resource assessment (WRA) consultants so they could conduct energy yield assessments (EYA). NREL took all of the completed EYAs, along with the operational data, and conducted an operational assessment to benchmark the EYA results against actual operational data. Given the large amounts of sensitive data required for this effort, as well as historical opposition to sharing data within industry, successful completion of Phase 1 represents an unprecedented milestone for industry data sharing. To improve the accuracy and confidence of pre-construction EYAs, wind plant owners and investors need better, more certain, energy yield predictions. The WP3 Benchmark Project is an industry-driven response to this reality. For the first time, industry has taken the important step of working together at scale, sharing valuable operational data with DOE and NREL in order to investigate the sources of bias and uncertainty in these energy estimates. This IC provides wind plant preconstruction and operational data to NREL in an organized and documented fashion and provides guidance and feedback as needed. The IC also provides introspection of the design of experiment, key metrics of success, data challenges, analysis best practices, and quality of results.

17 WIND ENERGY↗
DOE OSTI2022

Final Technical Report for "Cloud-based Low-Scaling Quantum Chemistry Simulations for Materials"

The objective of phase I was to develop the basic infrastructure for performing efficient hybrid DFT calculations on solid materials with a turnaround time amenable to the use in large-scale data-driven approaches. To fulfill the goal we have developed a novel algorithm that significantly reduces the cost of exchange matrix formation. The algorithm does so by making use of three ingredients (a) interpolative decomposition of the electron integrals (b) robust pseudospectral method and (c) occ-RI approach. Our published work demonstrates that the algorithm is orders of magnitude faster than any other hybrid-DFT method and is on the order of only 3-4 times slower than pure DFT. All steps of the algorithm developed during phase I can be accelerated to the point that, for large enough systems, the diagonalization of the Fock matrix will become the most expensive step. To treat such large systems we have created an interface to the ASCR-funded PEXSI library. The PEXSI method is used to compute the density matrix directly from the Fock matrix in a manner that preserves the sparsity of the local representation.

Shiozaki, Toru↗
DOE OSTI2022

Gamma Driven Catalysis of Ammonia

Experiments were conducted to investigate a passive production mechanism for the world’s most energy intensive commodity, ammonia. A novel method, gamma catalyzed ammonia production at ambient conditions, was investigated. Ammonia is currently produced through the highly energy intensive Haber-Fritz process, which requires an operation pressure of 400 atmosphere and 600 degrees Celsius. Due to the high demand and need for ammonia, the Haber-Bosch process consumes 25% all energy produced globally. Reported herein was an attempt to produce ammonia at ambient temperature (20 C) and ambient pressure (1 atm), through a novel process developed at PNNL, gamma driven catalysis of ammonia. Although the measurements of the ammonia production suggest wild success, reports in the literature by Gao et.al. suggest an experimental positive bias in the results. To rule out the potential positive bias, multiple additional production campaigns would be needed to with an alternate analysis technique such as ion chromatography, as suggested by Gao et.al. Unfortunately, due to this late determination of potential positive bias, the results of this study remain inconclusive to the feasibility of gamma driven catalysis of ammonia and more work is needed to describe the chemical evolution with time. The results are a first step and demonstrate that gamma-catalyst mediated reactions are possible. This represents a key opportunity to explore the fundamental chemistry of high band gap catalysts that can change the paradigm of radiation, transforming it from a waste to a valuable energy source.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗