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

An Open-Source Frequency-Domain Model for Floating Wind Turbine Design Optimization

A new frequency-domain dynamics model has been developed that uses open-source components to efficiently represent a complete floating wind turbine system. The model, called RAFT (Response Amplitudes of Floating Turbines), incorporates quasi-static mooring reactions, strip-theory and potential-flow hydrodynamics, blade-element-momentum aerodynamics, and linear turbine control. The formulation is compatible with a wide variety of support structure configurations and no manual or time-domain preprocessing steps are required, making RAFT very practical in design and optimization workflows. The model is applied to three reference floating wind turbine designs and its predictions are compared with results from time-domain OpenFAST simulations. There is good agreement in mean offsets as well the statistics and spectra of the dynamic response, verifying RAFT’s general suitability for floating wind analysis. Follow-on work will include verification of potential-flow and turbine-control features and application to optimization problems.

17 WIND ENERGY↗
DOE OSTI2020

Femtosecond Compression Dynamics and Timing Jitter Suppression in a THz-driven Electron Bunch Compressor

Herein we present the first demonstration of THz driven bunch compression and timing stabilization of a relativistic electron beam. Quasi-single-cycle strong field THz radiation is used in a shorted parallel-plate structure to compress a few-fC beam with 2.5 MeV kinetic energy by a factor of 2.7, producing a 39 fs rms bunch length and a reduction in timing jitter by more than a factor of 2 to 31 fs rms. This THz driven technique offers a significant improvement to beam performance for applications like ultrafast electron diffraction, providing a critical step towards unprecedented timing resolution in ultrafast sciences, and other accelerator applications using femtosecond-scale electron beams.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Free energy difference to create the M-OH* intermediate of the oxygen evolution reaction by time-resolved optical spectroscopy

Theoretical descriptors differentiate the catalytic activity of materials for the oxygen evolution reaction by the strength of oxygen binding in the reactive intermediate created upon electron transfer. Recently, time-resolved spectroscopy of a photo-electrochemically driven oxygen evolution reaction followed the vibrational and optical spectra of this intermediate, denoted M-OH*. However, these inherently kinetic experiments have not been connected to the relevant thermodynamic quantities. Here we discover that picosecond optical spectra of the Ti-OH* population on lightly doped SrTiO 3 are ordered by the surface hydroxylation. A Langmuir isotherm as a function of pH extracts an effective equilibrium constant relatable to the free energy difference of the first oxygen evolution reaction step. Thus, time-resolved spectroscopy of the catalytic surface reveals both kinetic and energetic information of elementary reaction steps, which provides a critical new connection between theory and experiment by which to tailor the pathway of water oxidation and other surface reactions.

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

Dynamics retrieval from stochastically weighted incomplete data by low-pass spectral analysis

Time-resolved serial femtosecond crystallography (TR-SFX) provides access to protein dynamics on sub-picosecond timescales, and with atomic resolution. Due to the nature of the experiment, these datasets are often highly incomplete and the measured diffracted intensities are affected by partiality. To tackle these issues, one established procedure is that of splitting the data into time bins, and averaging the multiple measurements of equivalent reflections within each bin. This binning and averaging often involve a loss of information. Here, we propose an alternative approach, which we call low-pass spectral analysis (LPSA). In this method, the data are projected onto the subspace defined by a set of trigonometric functions, with frequencies up to a certain cutoff. This approach attenuates undesirable high-frequency features and facilitates retrieving the underlying dynamics. A time-lagged embedding step can be included prior to subspace projection to improve the stability of the results with respect to the parameters involved. Subsequent modal decomposition allows to produce a low-rank description of the system's evolution. Using a synthetic time-evolving model with incomplete and partial observations, we analyze the LPSA results in terms of quality of the retrieved signal, as a function of the parameters involved. We compare the performance of LPSA to that of a range of other sophisticated data analysis techniques. We show that LPSA allows to achieve excellent dynamics reconstruction at modest computational cost. Finally, we demonstrate the superiority of dynamics retrieval by LPSA compared to time binning and merging, which is, to date, the most commonly used method to extract dynamical information from TR-SFX data.

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

Time-Varying Feedback Optimization for Quadratic Programs with Heterogeneous Gradient Step Sizes

Online feedback-based optimization has become a promising framework for real-time optimization and control of complex engineering systems. This tutorial paper surveys the recent advances in the field as well as provides novel convergence results for primal-dual online algorithms with heterogeneous step sizes for different elements of the gradient. The analysis is performed for quadratic programs and the approach is illustrated on applications for adaptive step-size and model-free online algorithms, in the context of optimal control of modern power systems.

gradient methods↗
DOE OSTI2023

Nanosecond time-resolved infrared spectroscopy for the study of electron transfer in photosystem I

Abstract Microsecond time-resolved step-scan FTIR difference spectroscopy was used to study photosystem I (PSI) from Thermosynechococcus vestitus BP-1 ( T. vestitus, formerly known as T. elongatus ) at 77 K. In addition, photoaccumulated (P700 + –P700) FTIR difference spectra were obtained at both 77 and 293 K. The FTIR difference spectra are presented here for the first time. To extend upon these FTIR studies nanosecond time-resolved infrared difference spectroscopy was also used to study PSI from T. vestitus at 296 K. Nanosecond infrared spectroscopy has never been used to study PSI samples at physiological temperatures, and here it is shown that such an approach has great value as it allows a direct probe of electron transfer down both branches in PSI. In PSI at 296 K, the infrared flash-induced absorption changes indicate electron transfer down the B- and A-branches is characterized by time constants of 33 and 364 ns, respectively, in good agreement with visible spectroscopy studies. These time constants are associated with forward electron transfer from A 1 – to F X on the B- and A-branches, respectively. At several infrared wavelengths flash-induced absorption changes at 296 K recover in tens to hundreds of milliseconds. The dominant decay phase is characterized by a lifetime of 128 ms. These millisecond changes are assigned to radical pair recombination reactions, with the changes being associated primarily with P700 + rereduction. This conclusion follows from the observation that the millisecond infrared spectrum is very similar to the photoaccumulated (P700 + –P700) FTIR difference spectrum.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Geochemical Impact of Acid Spearhead and Slickwater Stimulation on Wolfcamp Shale from the Hydraulic Fracturing Test Site

The Hydraulic Fracturing Test Site 1 (HFTS-1) was a field study performed in the Wolfcamp Formation in the West Texas Permian (Midland) Basin, USA, with a focus on improving the efficiency of hydraulic fracturing. Investigating site-specific rock-fluid geochemical interactions during hydraulic fracturing is an important step to understanding the impact on formation shale porosity, permeability, and long-term shale gas production. During field operations in this region, hydraulic fracturing fluid (HFF) injection usually starts with a concentrated acid spearhead for rapid rock dissolution, followed by the injection of near-neutral pH slickwater containing chemicals and proppants. A multistep sequential injection approach was used to investigate different stages of rock-fluid interactions. The carbonate content in the host rock is important when acid spearhead is considered, as carbonate mineral dissolution is rapid and can result in porosity and permeability changes in the shale matrix. Here, in this study, we designed flow-through experiments using fractured carbonate-rich and clay-rich Wolfcamp shale cores with (1) a short-time acid soaking step and (2) a long-term slickwater flow-through step to simulate the injection method used at HFTS-1. The fluid chemistry was analyzed. A thorough mineralogical progression [e.g., Calcium (Ca) dissolution and iron (Fe) redox progression] in the cores during HFF injection was also characterized and imaged by synchrotron microprobe. Reactive transport modeling was performed based on the experimental setup. The results showed that the acid spearhead is a crucial step in creating a reaction front by mineral dissolution, especially in carbonate-rich shales. A slight layer of ferrihydrite precipitated during the slickwater flow-through period. This study provides insights into potential geochemical impact due to hydraulic fracturing operations in the Permian Basin.

clastic rock↗
DOE OSTI2024

Parallel computing for power system climate resiliency: Solving a large-scale stochastic capacity expansion problem with mpi-sppy

Here we propose a nodal stochastic generation and transmission expansion planning model that incorporates the output from high-resolution global climate models through load and generation availability scenarios. We implement our model in Pyomo and perform computational studies on a realistically-sized test case of the California electric grid in a high performance computing environment. We propose model reformulations and algorithm tuning to efficiently solve this large problem using a variant of the Progressive Hedging Algorithm. We utilize the parallelization capabilities and overall versatility of mpi-sppy, exploiting its hub-and-spoke architecture to concurrently obtain inner and outer bounds on an optimal expansion plan. Initial results show that instances with 360 representative days on a system with over 8,000 buses can be solved to within 5% of optimality in under 4 h of wall clock time, a first step towards solving a large-scale power system expansion planning problem across a wide range of climate-informed operational scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Extrusion of Unhomogenized Castings of 7075 Aluminum via ShAPE

Extrusion of unhomogenized castings of 7075 aluminum (Al) alloy has been accomplished using shear assisted processing and extrusion (ShAPE). The simultaneous plastic deformation and heat generation during ShAPE rapidly fracture and dissolve interdendritic and intragranular secondary phases of Al-Zn-Mg-Cu, accomplishing homogenization in seconds rather than many hours in a furnace prior to extrusion. ShAPE thereby eliminates the energy-intensive and time-consuming homogenization step required to prepare as-cast microstructures for conventional extrusion. Concurrently, extensive grain refinement occurs due to gradient activation of dynamic recrystallization during ShAPE, which facilitates a threefold increase in extrusion speed compared to the conventional extrusion method. Evident enhancement of the mechanical properties of ShAPE+T6 samples is achieved, compared to the ASTM standard values for conventional extrusion products. The microstructural evolution pathway, involving the processes of grain refinement and secondary particle dissolution, is explained using scanning and transmission electron microscopy.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Vortical interactions in turbulent thermoacoustic systems

This study examines the dynamics of vortical interactions and their implications for mitigating thermoacoustic instability in a turbulent combustor. The regions of intense vortical interactions are identified as vortical communities in the network space of weighted directed vortical networks constructed from two-dimensional experimental velocity data. One can expect vortical interactions in the combustor to be strongest near the moment of vortex shedding, as the shed vortices gradually weaken due to dissipation while convecting downstream. However, we show that, during the state of thermoacoustic instability, there is a non-trivial consistent phase lag of approximately 52° between the shedding of the coherent structures from the backward-facing step and the time instant when the vortical interactions attain their local maximum value. We explain this phase lag by investigating the correlation between acoustic pressure fluctuations, spatio-temporal dynamics of coherent structures and vortical interactions in the reaction field of the combustor. We also show the aperiodic variation of vortical interactions during the states of combustion noise and aperiodic epochs of intermittency. Furthermore, the spatio-temporal evolution of pairs of vortical communities with the maximum inter-community interactions provides insight into explaining the critical regions detected in the reaction field during the states of intermittency and thermoacoustic instability, also identified in previous studies. As a result, we further show that the most efficient suppression of thermoacoustic instability via air microjet injection is achieved when steady air jets are introduced to disrupt the maximum inter-community interactions present during the state of thermoacoustic instability.

Sahay, Ankit [Indian Institute of Technology Madra↗
DOE OSTI2026

Structural Conservation of the A 1 Binding Site in Photosystem I across Cyanobacteria and Green Algae

Time-resolved step-scan Fourier transform infrared (FTIR) difference spectroscopy was used to obtain (A 1 − − A 1 ) FTIR difference spectra from photosystem I (PSI) samples isolated from eight phylogenetically diverse cyanobacterial strains and one green alga, totaling 13 PSI preparations. These included samples from cells grown under farred light and PSI in monomeric, dimeric, trimeric, and tetrameric states. Spectral profiles were shown to be independent of oligomeric state. Remarkably, all (A 1 − − A 1 ) FTIR difference spectra exhibited high similarity, underscoring the robustness of the technique and indicating minimal experimental variability. This congruence reveals a highly conserved environment for the phylloquinone cofactor at the A 1 binding site across diverse taxa. Conserved bands associated with the A 0 pigment further suggest structural continuity from A 0 to A 1 . To leverage this consistency, we constructed a composite (A 1 − − A 1 ) FTIR difference spectrum by averaging all 13 spectra. This composite spectrum provides enhanced resolution, enabling unambiguous identification of previously unresolved bands. The fact that a highly resolved composite spectrum can be obtained by averaging demonstrates the similarity in the spectra from the different types of samples. Band assignments were refined using prior studies, yielding an improved spectral framework for future investigations of PSI electron transfer cofactors.

Charge transfer↗
DOE OSTI2022

Superradiance and Directional Exciton Migration in Metal–Organic Frameworks

Crystalline metal–organic frameworks (MOFs) are promising synthetic analogues of photosynthetic light-harvesting complexes (LHCs). The precise assembly of linkers (organic chromophores) around the topology-defined pores offers the evolution of unique photophysical behaviors that are reminiscence of LHCs. These include MOF excited states with photoabsorbed energy that is spatially dispersed over multiple linkers defining the molecular excitons. The multilinker molecular excitons display superradiance–a hallmark of coupled oscillators seen in LHCs–with radiative rate constant (k rad ) exceeding that of a single linker. Our theoretical model and experimental results on three zirconium MOFs, namely, PCN-222(Zn), NU-1000, and SIU-100, with similar topology but varying linkers suggest that the size of such molecular excitons depends on the electronic symmetry of the linker. This multilinker exciton model effectively predicts the energy transfer rate constant; corresponding single-step exciton hopping time, ranging from a few picoseconds in SIU-100 and NU-1000 to a few hundreds of picoseconds in PCN-222(Zn), matches well with the experimental data. The model also predicts the anisotropy of exciton displacement with preferential migration along the crystallographic c-axis. Overall, these findings establish various missing links defining the exciton size and dynamics in MOF-assembled linkers. Furthermore, the understandings will provide design principles, especially, positioning the catalysts or electrode relative to the linker orientation for low-density solar energy conversion systems.

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

Evaluation of Waveform Profiles for Traveling Wave Ion Mobility Separations in Structures for Lossless Ion Manipulations

Structures for lossless ion manipulations (SLIM) have recently enabled a powerful implementation of traveling wave ion mobility spectrometry (TWIMS) for ultra-high resolution separations. However, optimization of experimental parameters remains challenging. Most TWIMS separations have utilized square-shaped waveforms applied by time-dependent voltage stepping across a repeating sets of electrodes; however, alternative waveforms may provide further improvements to resolution. In this work, we characterize five waveforms (including square and sine) in terms of their transmission efficiency, IMS resolution, and resolving power, and explore the effects of TW amplitude and speed on the performance of each. We found, consistent with previous work, separations were generally improved with higher TW amplitudes, moderately improved by lower speeds (limited by ion ‘surfing’ with the waves), and found decreases in signal intensity at the extremes of operating conditions. The triangle and asymmetric ‘ramp forward’ shaped profiles were found to provide modestly greater resolution and resolving power, an observation we tentatively attribute to their relatively uniform fields and minimal low-field regions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

AI Driven Experiment Calibration and Control

One critical step on the path from data taking to physics analysis is calibration. For many experiments this step is both time consuming and computationally expensive. The AI Experimental Calibration and Control project seeks to address these issues, starting first with the GlueX Central Drift Chamber (CDC). We demonstrate the ability of a Gaussian Process to estimate the gain correction factor (GCF) of the GlueX CDC accurately, and also the uncertainty of this estimate. Using the estimated GCF, the developed system infers a new high voltage (HV) setting that stabilizes the GCF in the face of changing environmental conditions. This happens in near real time during data taking and produces data which are already approximately gain-calibrated, eliminating the cost of performing those calibrations which vary ±15% with fixed HV. We also demonstrate an implementation of an uncertainty aware system which exploits a key feature of a Gaussian process.

Britton, Thomas↗
DOE OSTI2024

Drift kinetic electrostatic simulations of the edge localized mode heat pulse

In the present work, electrostatic drift kinetic simulations of parallel plasma transport within the tokamak scrape-off layer (SOL) are conducted using the COGENT code. The SOL configuration is represented in one-dimensional slab geometry, incorporating a heat source localized in the midplane. The heat source parameters correspond to those characterizing edge-localized modes observed in the Joint European Torus (JET) tokamak. The numerical model includes kinetic treatment of both ions and electrons, a simplified model for the gyrokinetic Poisson equation that allows one to step over short time scales associated with fast electrostatic shear Alfvèn waves, and the logical sheath boundary condition (LSBC) that enforces global system quasineutrality. A third-order accurate LSBC is derived to be consistent with the third-order accurate upwind advection scheme utilized in the code, and it was shown to noticeably impact the simulation results, especially parallel heat flux at the target plate. The findings of this study are in agreement with results from preceding fluid and kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Identification Uncertainty in Inverse Material Model Parameter Determination: A Sensitivity‐Based Decision Process for Load Path Selection

This research proposes a sensitivity-based framework for selecting the optimal prescribed loading path for a biaxial cruciform specimen. Optimality here is determined by the direction and magnitude of the prescribed displacement that minimizes the influence of random noise on the material model parameter identification. Using simulated experimental data based on finite element simulation, in this work, we identify the material model parameters of a Ludwik hardening model and plane stress implementation of the Hill-48 yield criterion using finite element model updating (FEMU). Our analysis reveals that the identification (or estimator) uncertainty of model parameters depends on the displacement boundary conditions (i.e., loading sequence) and the ground-truth value of the individual parameters. Optimal experimental design (OED) criteria based on the Fisher information matrix were investigated to mitigate indecision in the choice of optimal load path when the identification uncertainty of different material model parameters optimized at different load paths. The determinant of the Fisher information matrix was chosen here as the more useful metric due to its ability to capture uncertainty of the most influential material model parameters. The proposed framework demonstrates potential for real-time automated load step selection using scalar criteria derived prior to mechanical loading. The framework can be generalized to other geometries, boundary conditions and material models, allowing this procedure to be utilized for different experimental configurations and materials.

Fayad, Samuel S. [University of Illinois at Urbana↗
DOE OSTI2025

Tailoring optical response of MXene thin films

Due to their attractive optical properties, 2D MXenes have garnered interest in nanophotonic and optoelectronic applications. However, tuning their properties typically requires the iterative synthesis of MXenes with a specific set of properties, such as the absorption band position, electronic conductivity, and dielectric constant. We demonstrate how to tailor the optical properties of MXene thin films over a broad 1500-nm wavelength range by mixing different ratios of highly conductive Ti 3 C 2 T x with poorly conductive Nb 2 CT x . By changing the MXene film composition, the epsilon-near-zero (ENZ) point, where the optical properties transit from dielectric to metallic, was varied in the spectral range from 1.1 to 2.6 µm. Additionally, we observed a reduction in absorption in some compositions compared to the absorption of the pure MXene films. Compared to other methods, this approach enables simple and continual tuning of MXene optical properties without requiring multiple time-consuming synthesis steps.

2D materials↗
DOE OSTI2023

Higher-Order Corrections to Optimisers based on Newton's Method

The Newton, Gauss–Newton and Levenberg–Marquardt methods all use the first derivative of a vector function (the Jacobian) to minimise its sum of squares. When the Jacobian matrix is ill-conditioned, the function varies much faster in some directions than others and the space of possible improvement in sum of squares becomes a long narrow ellipsoid in the linear model. This means that even a small amount of nonlinearity in the problem parameters can cause a proposed point far down the long axis of the ellipsoid to fall outside of the actual curved valley of improved values, even though it is quite nearby. This paper presents a differential equation that ‘follows’ these valleys, based on the technique of geodesic acceleration, which itself provides a 2 nd order improvement to the Levenberg–Marquardt iteration step. Higher derivatives of this equation are computed that allow n th order improvements to the optimisation methods to be derived. These higher-order accelerated methods up to 4 th order are tested numerically and shown to provide substantial reduction of both number of steps and computation time.

43 PARTICLE ACCELERATORS↗