Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “load data”

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 OSTI2023

Experimental investigation of advanced turbine control strategies and load-mitigation measures with a model-scale floating offshore wind turbine system

To advance the control co-design of offshore wind energy systems, the authors perform basin-scale experiments with a fully instrumented and actuated floating offshore wind turbine model. The model consists of a 1:70 scale performance-matched model of the International Energy Agency Wind Technology Collaboration Programme 15-MW reference turbine atop the VolturnUS-S semisubmersible platform. The Reference OpenSource Controller provides real-time blade pitch and generator torque control. Here we aim to develop an open data set for the validation of numerical models in predicting the influence of turbine control and load-mitigation measures. For this purpose, we measure the effects of advanced turbine control features, including peak shaving and floating feedback as well as hull-based control using tuned mass dampers on the system. Overall, the results demonstrate measurable and consistent influences from the control and load-mitigation measures, thus confirming the usefulness as a validation data set. Peak shaving attenuates the response to wind turbulence at near-rated wind speed. Floating feedback reduces the load and platform pitch motion associated with the negative damping induced by blade pitch control. The tuned mass dampers also attenuate the system response near the targeted frequencies under suitable conditions. We also identify detrimental side effects of each load-mitigation measure.

17 WIND ENERGY↗
DOE OSTI2021

Domain decomposition in the GPU-accelerated Shift Monte Carlo code

The GPU solver within the Shift continuous-energy Monte Carlo neutron transport code has been extended to provide domain decomposition in addition to domain replication to enable the solution of problems with memory requirements exceeding the capacity of a single GPU. The strategy follows the Multiple Set, Overlapping Domain (MSOD) approach that is used in Shift’s CPU solver and integrates into the event-based algorithm used for Shift’s GPU solver. Furthermore, the ability to assign processors to spatial domains non-uniformly has been maintained. In this work, two different approaches for communicating particle data between domains are considered, and multiple criteria for load balancing problems have been investigated. Numerical results are presented for both fresh and depleted small modular nuclear reactor (SMR) cores. A parallel efficiency of approximately 80% was achieved with up to 16 spatial domains measured relative to full domain replication. A scaling study on the Summit supercomputer demonstrates a weak scaling parallel efficiency of over 90% on over 24000 GPUs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A hybrid architecture for volt-var control in active distribution grids

Modern active distribution grids are characterized by the increasing penetration of distributed energy resources (DERs). The proper coordination and scheduling of a large numbers of these small-scale and spatially distributed DERs is necessary, and warrants the use of novel distributed approaches. In this paper, we propose a hybrid volt-var control architecture for the distribution grid, which leverages existing centralized and local approaches to planning, decision making, and control, and augments it with distributed optimization and distributed control for DER management. First, we propose a convex model to describe the power physics of distribution grids of meshed topology and unbalanced structure, based on current injection and McCormick Envelopes. Second, we employ the distributed proximal atomic coordination (PAC) algorithm to coordinate DERs to provide voltage support. We implement volt-var optimization by optimally coordinating DERs including PV smart inverters and demand response. We present results using the IEEE-34 bus network, using real data from a distribution feeder in Hawaii, to model load and PV generation. Different levels of DER penetration and objective functions are simulated. Finally, our results show the need for the coordination of DERs to improve voltage profiles, even in networks with existing voltage control devices. Further, we show the need for flexible reactive power capabilities to achieve desired grid performance.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

A nuclear-hydrogen hybrid energy system with large-scale storage: A study in optimal dispatch and economic performance in a real-world market

Here a study in optimal dispatch and economic analysis of a novel nuclear hybrid energy system (NHES) with large-scale hydrogen storage and a novel control scheme is conducted. The control scheme makes charge and discharge decisions by using the real-time electricity price relative to the historical price distribution. The pricing data is taken from the ERCOT Houston and West load zones with different price oscillations. Six case control studies with different levels of aggressiveness were chosen for each load zone for a total of 12 case studies. It is found that in a high price volatility load zone, a more aggressive control strategy maximizes revenue, while in a less oscillatory load zone, a moderate control strategy performs better. The NHES can store energy when the price is low and sell additional electricity when the price is high, as well as participate in the ancillary services market. Due to these advantages, the NHES can generate 10–40 million USD more revenue annually than a stand-alone nuclear plant. Due to higher costs, the NHES has a higher LCOE (81.36 USD/MWh) than a stand-alone plant (68.66 USD/MWh); however, the LCOE of the NHES is still comparable to other new forms of electricity generation. It is found that both the NHES and the stand-alone plant are not economical in the existing electricity market. Due to the extra revenue generated, the hybrid system has a positive annual cashflow whereas the stand-alone plant does not. It is found that the hybrid plant needs less cost subsidy than a stand-alone plant to become economical.

08 HYDROGEN↗
DOE OSTI2026

CO 2 Sorption in Moisture Swing Anion Exchange Resins for Direct Air Capture: Experimental Isotherm Determination and Modeling

Moisture swing (MS) sorption is a promising direct air capture (DAC) technology to achieve negative CO 2 emissions and counter global warming. Here, in this study, MS CO 2 sorption in a model MS sorbent, IRA900, was investigated. IRA900 is a macroporous commercial strong-base anion exchange resin (AER) with quaternary ammonium functional groups. A rigorous CO 2 desorption process was developed to desorb all CO 2 from the sample and obtain CO 2 sorption isotherms as a function of relative humidity, CO 2 partial pressure, and temperature. CO 2 sorption increased with decreasing relative humidity, decreasing temperature, and increasing CO 2 partial pressure. Remarkably, the total CO 2 uptake from the sorption isotherms matched the ion exchange capacity (IEC) of the material, suggesting a stoichiometry of one CO 2 molecule reacting per active site. An isotherm model for CO 2 sorption starting with empty sorption sites (i.e., with the AER in the OH – form) described the experimental data well, supporting the hypothesis that the bicarbonate-loaded AER fully unloaded to the OH – state following rigorous CO 2 desorption. CO 2 sorption kinetics were studied, and carbon diffusion coefficients were estimated as a function of %RH and CO 2 partial pressure. Overall, this study provides a roadmap for the systematic evaluation of MS CO 2 DAC materials.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Automation-Accelerated Electrolyte Design Mitigates Solubility Competition between Redox-Active Molecules and Supporting Salts

In nonaqueous redox-flow batteries (NRFBs), redox-active organic molecules (ROMs) and supporting salts compete for solvation sites, limiting achievable energy density. We combine automated high-throughput experimentation (HTE) with camera-based saturation monitoring and quantitative NMR to measure paired (ROM, salt) solubilities across single and mixed organic solvents. Using 2,1,3-benzothiadiazole (BTZ) with lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) as a model system, we find that a binary m-xylene/acetonitrile mixture dissolves ≈3 M of both BTZ and LiTFSI─surpassing the previously reported 2 M ceiling for neat acetonitrile─by leveraging complementary solvation (MX is BTZ-philic and salt-phobic; ACN stabilizes LiTFSI). A random-forest model (RMSE ≈ 0.24) trained on solvent descriptors highlights log P and salt concentration as dominant predictors and predicts MX/ACN ≈0.3/0.7 (v/v) to be near-optimal. These formulations retain practical viscosity and ∼5 mS·cm –1 conductivity at high loading. In conclusion, the workflow provides a reproducible, data-centric route to NRFB electrolyte design and motivates an open, standardized dual-solute solubility resource for accelerated electrolyte discovery.

Electrolytes↗
DOE OSTI2021

Frequency–Magnitude Statistics of Laboratory Foreshocks Vary With Shear Velocity, Fault Slip Rate, and Shear Stress

Understanding the nucleation phase of earthquakes is key for advancing earthquake hazard assessment and improving earthquake early warning systems. However, little progress has been made in this area due to a poor understanding of nucleation processes and incomplete seismic and fault zone measurements. The ability to integrate measured fault zone properties with seismic data could significantly improve our understanding of how earthquakes begin and whether there are systematic variations in seismic properties preceding failure. In this work, we use high-resolution laboratory measurements of fault zone properties along with acoustic emission data to document temporal variations of foreshock properties. Our data show that foreshock size increases with shear stress, loading rate, and fault slip rate. We propose that the preseismic fault slip rate and fault zone thickness (i.e., porosity) work in concert to modulate foreshock properties.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Interplay of Three-Dimensional Instabilities and Magnetic Reconnection in the Explosive Onset of Magnetospheric Substorms

Magnetospheric substorms are preceded by a slow growth phase of magnetic flux loading and current sheet thinning in the tail. Extensive data sets have provided evidence of the triggering of instabilities at substorm onset, including magnetic reconnection and ballooning instabilities. Using an exact kinetic magnetotail equilibrium we present particle-in-cell simulations which capture the explosive nature of substorms through a disruption of the dipolarization front by the ballooning instability. We use self-consistent particle tracking to determine the nonthermal particle acceleration mechanisms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Measurement of crack length in width tapered beam experiments

The width tapered beam method for measuring fracture toughness has been contemplated for use in measuring the adhesion of photovoltaic materials, is promoted by it being a viable method for use directly on modules rather than to engineered test specimens. However, precise determination of crack length is often difficult as failure often involves tendril formation, cavitation voids, and other plastic deformation at the crack front in addition to obfuscation by opaque materials. Because fracture energy varies as the square of crack length, imprecise crack tip length results in large measurement uncertainties. In this work, we develop a method to estimate the crack tip length using the relationship between the applied force and the displacement of the load frame. The crack tip is determined for every data point and can be defined from a calibration curve related to the amount of energy put into the beam, or from a theoretically idealized estimation. Once a calibration curve is set up and the equations programmed, the fracture energy as a function of position on the newly formed crack interface can be accurately estimated. The test additionally provides much easier analysis and insight into the joint failure through more precise knowledge of fracture toughness at each de-bonded area.

42 ENGINEERING↗
DOE OSTI2025

Compression–tension cell with sample manipulator for in situ X‐ray nanotomography experiments

In situ X-ray nanotomography experiments where tensile or compressive force is applied on the sample require specialized equipment. A compression-tension device with fluid flow-through capability has been designed for X-ray nanotomography beamlines. The compression-tension cell is equipped with a triaxial stage for sample alignment and a high sensitivity loadcell for measurement of applied force. To handle the <100 µm samples used for X-ray nanotomography imaging and for loading samples on the compression-tension cell a sample manipulator has been built. The sample manipulator is capable of selecting a single <100 µm particle for nanotomography scanning while viewing multiple samples under an optical microscope. To test the functionality of these two devices an initial compression experiment involving two glass beads was performed. To demonstrate instrument stability two spherical glass beads were compressed from a no load condition until one of the beads fractured. Nanotomography data were collected at each step of increasing compressive force. The experimentally observed contact area of the spherical glass beads was compared with the theoretical estimate using the Hertz analysis. To demonstrate the fluid flow capability, two calcite grains were compressed against each other under a calcite saturated solution. Surface topological changes were observed for the stressed grain contact area.

X-ray tomography↗
DOE OSTI2026

qSIEVE: Efficient qLDPC Memory via Systolic Movement in Atom Arrays

As quantum machines have scaled up in their number of qubits, significant research has turned towards increasing their fidelity with quantum error correction codes. Although promising results have been shown with the surface code, which only requires near-neighbor connections between qubits, the high qubit overhead of such local codes promises to be problematic. Consequently, recent work has explored non-local quantum LDPC (qLDPC) codes, which have good asymptotic encoding rates. Despite theoretical progress, hardware implementations of these codes have been a longstanding challenge. At the experimental level, demonstrations of movement based communication on atom arrays suggest this is a powerful new primitive to achieve non-local connectivity. Leveraging this, we present a protocol for implementing non-local qLDPC codes in hardware. Our protocol, qSIEVE, is a co-design of such codes with movement in atom arrays. qSIEVE defines a restricted family of qLDPC codes that can be implemented efficiently with systolic movement. We then quantify the utility of qSIEVE in the context of a complete fault tolerant architecture. We compare the cost of implementing benchmark programs in a standard, surface code only architecture and a mixed architecture where data is stored in qLDPC memory with qSIEVE and loaded to surface codes for computation.

Quantum error correction↗
DOE OSTI2020

Neutron Spectrum and 10B(n,α) Measurements in ZPPR-15

Neutron spectra were measured at the core center with proton recoil spectrometers in ZPPR15A loading 66, ZPPR-15B loading 90 and ZPPR loading 186. In addition, neutron spectra were measured in the SSNA shield configuration in ZPPR-15B loading 121 and in the BCNA shield configuration in ZPPR-15B loading 126. The ZPPR-15 neutron spectrum measurement and data processing techniques were described and analyzed. Total uncertainties were derived for all five spectrum measurements. A back-to-back proportional counter was used to measure the 10 B(n,α) reaction rate at six locations in the simulated shield, radial reflector and outer core in the ZPPR-15B BCNA shield experiment. The measurements were described and analyzed. A total uncertainty was derived for the 10 B(n,α) reaction rate measurements. The ZPPR-15 experimental records were used to develop detailed as-built Monte Carlo models for the five neutron spectrum measurement configurations and for the six 10 B(n,α) measurement configurations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Theoretical Yields for 188 W Production in the High Flux Isotope Reactor

Accurately predicting the production yield of 188 W and 191 Os from enriched 186W metal rings irradiated in HFIR is not straightforward. Two key factors can lead to significant uncertainty or inaccuracy in yield predictions: The effect of neutron flux depression in the dense tungsten metal target. Inaccuracy of the currently published 187 W thermal neutron absorption cross section. This manuscript reports the calculated yields for 188 W, 191 Os, and total osmium mass for the NM-888 target that was irradiated in High Flux Isotope Reactor (HFIR) Cycles 487 and 488. Three different yield calculation methods were used and are described herein. The results are presented in terms of total yield and in terms of yield per target mass. This allows for an approximate estimation of yields for future targets of varying mass loading but with similar irradiation profiles. For more accurate prediction yields, target-specific models that account for target loading, target positioning within HFIR, and accurate HFIR cycle data are required.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

High Solids Performance Testing in a Scaled TSCR System

The Tank Side Cesium Removal (TSCR) project is a technology demonstration that will pretreat Hanford tank waste supernatant in support of the Direct Feed Low-Activity Waste (DFLAW) mission. The TSCR system employs two key separation technologies: dead-end filtration (DEF) and ion exchange (IX) using crystalline silicotitanate (CST) media. DEF will be used to remove undissolved solids from tank waste to protect the functionality of the IX columns and the IX system will remove Cs-137 from tank waste. The separation technologies (DEF and IX) used in TSCR are technically mature and have also been successfully deployed at the Savannah River site in a similar facility known as the Tank Closure Cesium Removal (TCCR) system. While testing with simulants and real waste has been successfully performed under conditions expected during the initial operation of TSCR, test data is absent for assessing off normal high solids loading that may be in the TSCR waste feed. Normal TSCR treatment operations are expected to handle wastes with solids content on the order of 200 ppm, and off normal solids loading could be much larger than the nominal level. The testing program described in this report was conducted to understand the consequence of operating the TSCR system at elevated solids loadings up to the high-solids limit of 15,000 ppm [i.e., 1.5-wt%] identified in the TSCR design basis. Although the system is not required to make throughput above the nominal solids loading, the testing was intended to provide important information related to potential off normal operations. At off normal levels near the high-solids limit, there are potential implications for TSCR performance in the areas of throughput, DEF pressure drop, filter backflush frequency, and IX column pressure drop. In addition, intrusion of solids into the IX column was postulated to impact the Cs-137 loading behavior by promoting channeling or flow maldistribution in the column; since the magnitude of the postulated effect was unknown, assessing it was also of interest. The testing was performed using representative waste simulants and a prototypic, integrated TSCR system designed and assembled specifically to conduct the high solids performance assessment. Overall, the scaled TSCR testing demonstrated that full-scale unit operations can succeed in fulfilling their processing objectives in the presence of solids up to 3,000 ppm, but there are potential performance challenges to filtration operations at solids loadings as low as ~500 ppm. The severity of the challenge is likely to be dependent on the type and size distribution of solids, of which the current testing only examined a single type and size distribution. To provide some flexibility for future full-scale operations, the results of the testing suggest two possible risk reduction strategies that can be implemented without any changes in TSCR design or configuration. One option would be to enact an administrative limit on the solids loading to protect TSCR from feeds that are likely to require a high DEF swap frequency. Another option is to permit operation of the DEFs at differential pressures greater than 2 psid before swapping filters. The selection of a higher differential pressure target is not anticipated to adversely impact DEF backflushing efficacy and would reduce both swap frequency and the amount of waste sent to AP-108.

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

NREL Price Series Developed for the ARPA-E FLECCS Program

The price data for four regions (CAISO, ERCOT, MISO-W, and PJM-W) are developed using the ReEDS to PLEXOS conversion as described in (Gagnon et al. 2020). The reference ReEDS case chosen is based on the 2020 Standard Scenario Mid-case, which uses the 2020 ReEDS model version (Cole et al. 2020; Ho et al. 2021). All ReEDS model inputs use 2020 Standard Scenarios Mid-case assumptions except for CO2 prices, which are implemented as linearly increasing CO2 price trajectories beginning at $0/tCO2 in 2020 and ending at either $100/tCO2 or $150/tCO2 in 2035 to dive capacity expansion towards a low-carbon system that could support CCS deployment. However, these scenarios also prohibit CCS deployment in this time frame so that resulting price data are not influenced by the deployment and operation of CCS itself. Implementing the ReEDS to PLEXOS conversion tool, PLEXOS is then simulated using the 2035 ReEDS infrastructure for both CO2 price scenarios, with the following model version and setup: PLEXOS Version: 8.2 Solver: Xpress-MP 35.01.01 Mixed integer optimization relative gap 1% System configuration: • Total number of nodes: 134 (consistent with ReEDS balancing areas) • Line losses enforced using piecewise linear approximation • Energy dump was enabled Price data is aggregated to the ISO/RTO level using load-weighted averages. References: Cole, Wesley, Sean Corcoran, Nathaniel Gates, Daniel Mai, Trieu, and Paritosh Das. 2020. “2020 Standard Scenarios Report: A U.S. Electricity Sector Outlook.” NREL/TP-6A20-77442. Golden, CO: National Renewable Energy Laboratory. https://www.nrel.gov/docs/fy21osti/77442.pdf. Gagnon, Pieter, Will Frazier, Elaine Hale, and Wesley Cole. 2020. “Cambium Documentation: Version 2020.” NREL/TP-6A20-78239. National Renewable Energy Lab. (NREL), Golden, CO (United States). https://doi.org/10.2172/1734551. Ho, Jonathan, Jonathon Becker, Maxwell Brown, Patrick Brown, Ilya (ORCID:0000000284917814) Chernyakhovskiy, Stuart Cohen, Wesley (ORCID:000000029194065X) Cole, et al. 2021. “Regional Energy Deployment System (ReEDS) Model Documentation: Version 2020.” NREL/TP-6A20-78195. Golden, CO: National Renewable Energy Laboratory. https://doi.org/10.2172/1788425.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Getting to the Core of PARAFAC2, A Nonnegative Approach

In this paper, the authors present a novel method of performing PARAFAC2 factorization of three-way data using a compact representation of that data. In the standard PARAFAC2 algorithm, two modes of the data are recovered directly during the decomposition while the third mode is returned as a transformation matrix, which is then used to rotate sets of orthogonal third-mode basis factors into interpretable factors. In our new method, the data are first decomposed into a core matrix and orthogonal factor loading matrices in the first two modes as well as sets of orthogonal factors in the third mode (as in standard PARAFAC2). The core matrix is then decomposed using a the standard PARAFAC2 strategy to produce transformation matrices in all three modes. The algorithm is particularly useful for very large data sets and essentially permits imposition of nonnegativity in all three modes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Local measurements of the spatial magnetic field distribution in a z-pinch plasma during and near stagnation using polarization spectroscopy

We present here the detailed measurements of radial distribution of the magnetic field in a gas-puff z-pinch plasma at the final stages of the implosion phase and at stagnation. While the measurements are chordal, the radial distribution of different charge states was utilized to measure the magnetic field locally for certain radii, so that unlike chordal measurements in general, the magnetic field radial distribution was obtained with no need for the Abel inversion of the data. The distribution was measured using the Zeeman effect via a novel spectroscopic technique, at several axial locations, and demonstrates striking features such as the peak field remaining at a radius much larger than the stagnation radius at all times. Furthermore, while the distribution observed is sometimes monotonic with respect to the radius, it is often not, a behavior that can be linked to 2D features in the plasma column resulting from the Rayleigh–Taylor instability. The current flowing through the stagnating plasma was found to be a small fraction of the total current, resulting in clearly insufficient magnetic pressure to balance the plasma pressure at stagnation. The magnetic field data, taken over several axial positions, are used to obtain the true inductance in the imploding plasma for the first time; it is found that the data cannot explain the current turnover at stagnation. A simulation with the MACH2-Tabular Collisional-Radiative Equilibrium magnetohydrodynamics code in the r–z plane shows that the peak of the magnetic field pinches to a much smaller radius than is observed in the spectroscopic data. Furthermore, the depth of the computed current turnover at stagnation is smaller than the measured one. The two observed features of a radially extended magnetic field at stagnation together with a deep current turnover are a challenge to match in simulations. Various calculations and estimates of the inductive and resistive load voltages are examined to ascertain if they are responsible for the observed current notch. Lastly, the results demonstrate that the knowledge of the true inductance in the driven load requires such magnetic-field-distribution measurements and that imaging data or electrical measurements are insufficient.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗