Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “power modeling and estimation”

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 1,261 records · Page 70

Optimization of Energy Storage System Economics and Controls by Incorporating Battery Degradation Costs in REopt

The use of stationary electrochemical energy storage systems utilizing lithium-ion batteries has increased rapidly as the production scale and price for lithium-ion batteries has decreased. These energy storage systems are crucial for maintaining grid resiliency, especially for grids operating with high penetration of renewable energy generation assets or for with a variety of distributed energy generation and storage systems. One challenging factor for the development of battery energy storage systems is estimating the proper sizing, in terms of both power and energy, that minimizes total costs over the lifetime of the systems; this calculation is difficult in simple cases, where a battery is costed independently, but is extremely challenging when building loads and electrical generation by photovoltaic resources are also considered. REopt is a techoeconomic optimization tool developed by NREL to address these challenges. Previously, battery degradation has been priced by simply assuming a 10-year replacement schedule for battery systems. However, this does not account for varying degradation trends observed across real-world batteries, or allow for batteries to be operated in a degradation-aware manner that optimizes battery dispatch based on operating costs. This work incorporates a battery life model into REopt. This battery life model is simple, so that it may be solvable within the constrains of a mixed-integer linear optimization problem, but is fit to accelerated aging data recorded in the lab. To achieve the best possible accuracy for lifetime estimates given these constraints, parameters for the battery life model in REopt are estimated by fitting 20-year simulations of battery life after identifying state-space battery degradation model from accelerated aging data. Comparisons of battery life predicted in REopt and from the state-space battery degradation model to ensure validity of lifetime estimates made by REopt. Battery life and cost is optimized by controlling three decision to minimize system life cost: battery sizing, daily state-of-charge, and daily energy-throughput. The cost of battery degradation as a function of these control variables is then estimated assuming two possible maintenance strategies: replacement, where the entire battery system is replaced if cell reach an end-of-life capacity threshold; and augmentation, which establishes a fund to pay for continual purchase of new batteries to maintain the initial energy capacity of the system. These two strategies offer conservative (for replacement) and optimistic (for augmentation) bounds for total system cost. The degradation cost incurred by these strategies is then used to control battery dispatch decisions, operating the battery in a degradation-aware manner that maximizes battery lifetime while also providing energy when favorable. Because the mixed-integer linear program has perfect foresight of future energy needs, batteries with degradation costs are always operated using 'just-in-time' charging, which is unrealistic, as no energy is left in the storage system to perform other energy services or to serve as emergency back-up power. To combat this, an inequality constraint on the average annual state-of-charge is imposed, and the sensitivity of system cost to average stored energy, e.g., the cost of system resiliency, can be quantified. Analysis of results has several conclusions, for instance, oversizing of battery storage systems is not a cost burden when battery storage is an optimal solution, as any additional battery capacity can simply be utilized to avoid costs of purchasing energy from a utility.

battery↗

100% Carbon-Free Electricity for Sandia NM and KAFB Using Concentrating Solar Power (CSP) (SAND Report)

This report provides a design study to produce 100% carbon-free electricity for Sandia NM and Kirtland Air Force Base (KAFB) using concentrating solar power (CSP). Annual electricity requirements for both Sandia and KAFB are presented, along with specific load centers that consume a significant and continuous amount of energy. CSP plant designs of 50 MW and 100 MW are then discussed to meet the needs of Sandia NM and the combined electrical needs of both Sandia NM and KAFB. Probabilistic modeling is performed to evaluate inherent uncertainties in performance and cost parameters on total construction costs and the levelized cost of electricity. Total overnight construction costs are expected to range between ~$300M - $400M for the 50 MW CSP plant and between ~$500M - $800M for the 100 MW plant. Annual operations and maintenance (O&M) costs are estimated together with potential offsets in electrical costs and CO 2 emissions. Other considerations such as interconnections, land use and permitting, funding options, and potential agreements and partnerships with Public Service Company of New Mexico (PNM), Western Area Power Administration (WAPA), and other entities are also discussed.

14 SOLAR ENERGY↗

Nonlinear burn control in ITER using adaptive allocation of actuators with uncertain dynamics

Abstract ITER will be the first tokamak to sustain a fusion-producing, or burning, plasma. If the plasma temperature were to inadvertently rise in this burning regime, the positive correlation between temperature and the fusion reaction rate would establish a destabilizing positive feedback loop. Careful regulation of the plasma’s temperature and density, or burn control, is required to prevent these potentially reactor-damaging thermal excursions, neutralize disturbances and improve performance. In this work, a Lyapunov-based burn controller is designed using a full zero-dimensional nonlinear model. An adaptive estimator manages destabilizing uncertainties in the plasma confinement properties and the particle recycling conditions (caused by plasma–wall interactions). The controller regulates the plasma density with requests for deuterium and tritium particle injections. In ITER-like plasmas, the fusion-born alpha particles will primarily heat the plasma electrons, resulting in different electron and ion temperatures in the core. By considering separate response models for the electron and ion energies, the proposed controller can independently regulate the electron and ion temperatures by requesting that different amounts of auxiliary power be delivered to the electrons and ions. These two commands for a specific control effort (electron and ion heating) are sent to an actuator allocation module that optimally maps them to the heating actuators available to ITER: an electron cyclotron heating system (20 MW), an ion cyclotron heating system (20 MW), and two neutral beam injectors (16.5 MW each). Two different actuator allocators are presented in this work. The first actuator allocator finds the optimal mapping by solving a convex quadratic program that includes actuator saturation and rate limits. It is nonadaptive and assumes that the mapping between the commanded control efforts and the allocated actuators (i.e. the effector model) contains no uncertainties. The second actuator allocation module has an adaptive estimator to handle uncertainties in the effector model. This uncertainty includes actuator efficiencies, the fractions of neutral beam heating that are deposited into the plasma electrons and ions, and the tritium concentration of the fueling pellets. Furthermore, the adaptive allocator considers actuator dynamics (actuation lag) that contain uncertainty. This adaptive allocation algorithm is more computationally efficient than the aforementioned nonadaptive allocator because it is computed using dynamic update laws so that finding the solution to a static optimization problem is not required at every time step. A simulation study assesses the performance of the proposed adaptive burn controller augmented with each of the actuator allocation modules.

Physics↗

Man‐in‐the‐middle attacks and defence in a power system cyber‐physical testbed

Abstract Man‐in‐The‐Middle (MiTM) attacks present numerous threats to a smart grid. In a MiTM attack, an intruder embeds itself within a conversation between two devices to either eavesdrop or impersonate one of the devices, making it appear to be a normal exchange of information. Thus, the intruder can perform false data injection (FDI) and false command injection (FCI) attacks that can compromise power system operations, such as state estimation, economic dispatch, and automatic generation control (AGC). Very few researchers have focused on MiTM methods that are difficult to detect within a smart grid. To address this, we are designing and implementing multi‐stage MiTM intrusions in an emulation‐based cyber‐physical power system testbed against a large‐scale synthetic grid model to demonstrate how such attacks can cause physical contingencies such as misguided operation and false measurements. MiTM intrusions create FCI, FDI, and replay attacks in this synthetic power grid. This work enables stakeholders to defend against these stealthy attacks, and we present detection mechanisms that are developed using multiple alerts from intrusion detection systems and network monitoring tools. Our contribution will enable other smart grid security researchers and industry to develop further detection mechanisms for inconspicuous MiTM attacks.

Wlazlo, Patrick↗

The helium chromosphere, coronal holes, and stellar X-rays

The solar chromosphere at the limb seen in D3 is an irregular bright band 1000 km thick with a dark band 1000 km thick beneath. The D3 chromosphere disappears in coronal holes. It is shown that the D3 emission, as well as the other He I and He II lines, can be explained quantitatively by photoionization by coronal back-radiation. A Chapman layer with N(He)H = 5 times 10 to the 17th power is formed near tau = 1 in the He I and He II continua. The chromospheric He emission or absorption is weak in coronal holes because there is no coronal back-radiation. Based on this model, the soft X-ray flux from stars with He 10830-A absorption lines is estimated as proportional to the 10830-A equivalent width and the apparent area.

Zirin, H.↗

Flow-field surveys on the windward side of the NASA 040A space shuttle orbiter at 31 deg angle of attack and Mach 20 in helium

Pitot pressure and flow angle distributions in the windward flow field of the NASA 040A space shuttle orbiter configuration and surface pressures were measured, at a Mach number of 20 and an angle of attack of 31 deg. The free stream Reynolds number, based on model length, was 5.39 x 10 to the 6th power. Results show that cores of high pitot pressure, which are related to the body-shock-wing-shock intersections, occur on the windward plane of symmetry in the vicinity of the wing-body junction and near midspan on the wing. Theoretical estimates of the flow field pitot pressures show that conical flow values for the windward plane of symmetry surface are representative of the average level over the entire lower surface.

Ashby, G. C., Jr.↗

Method for retrieving the true backscattering coefficient from measurements with a real antenna

Coherent and incoherent scattering produce contributions to the power backscattered by a surface at angles near normal incidence, and these contributions are further weighed by the antenna pattern. On the basis of a theoretical model describing the backscattering from a rough soil surface, a procedure is developed for retrieving the true angular pattern of the backscattering coefficient from its measured estimates, where these estimates are based on the usual assumption that the coefficient is approximately constant over the angular extent of the antenna beam for narrow beam systems. The retrieved backscattering coefficient patterns were then used to evaluate the dependence of the coefficient on soil surface roughness at 1.5, 4.25, and 7.25 GHz.

Ulaby, F. T.↗

Radar Reconnaissance of Near-Earth Asteroids

Groundbased radar observations of NEAs can help identify space resources with commercial potential and can dramatically reduce the cost and risk of the initial spacecraft missions to those objects. The near-Earth asteroid (NEA) population is thought to contain approximately 1500 objects as large as a kilometer, approximately 300,000 as large as 100 meters, and more than 100,000,000 as large as 10 meters. More than 10% of the NEAs are more accessible in terms of mission delta-V (i.e., fuel required) than the Moon, Mars, or the moons of Mars. Fewer than 1000 NEAs have been found, but the discovery rate is increasing rapidly. Once an asteroid is discovered (necessarily by wide-field optical telescopes), radar can provide otherwise unavailable information about its size, shape, spin state, and surface properties if it approaches within the range of the Goldstone (California) or Arecibo (Puerto Rico) radar telescopes. Asteroids generally appear as unresolved points through groundbased optical telescopes, but radar measurements of the distribution of echo power in time delay (range) and Doppler frequency (radial velocity) can yield images with resolution as fine as a decameter. Image sequences that furnish adequate orientational coverage can be inverted to construct geologically detailed 3-D models, to define the rotation state precisely, and to constrain the object's internal density distribution. Estimates of radar scattering properties characterize the surface's cm-to-m-scale roughness as well as its bulk density, which for asteroids depends primarily on porosity and metal concentration. A useful spin-off of radar detection is orbit refinement that simplifies navigation of flyby and rendezvous spacecraft.

Ostro, S. J.↗

Radar Reconnaissance of Near-Earth Asteroids

Groundbased radar observations of NEAs can help identify space resources with commercial potential and can dramatically reduce the cost and risk of the initial spacecraft missions to those objects. The near Earth asteroid (NEA) population is thought to contain approx. 1500 objects as large as a kilometer, approx. 300,000 as large as 100 meters, and more than 100,000,000 as large as 10 meters. More than 10% of the NEAs are more accessible in terms of mission delta-V (i.e., fuel required) than the Moon, Mars, or the moons of Mars. Fewer than 1000 NEAs have been found, but the discovery rate is increasing rapidly. Once an asteroid is discovered (necessarily by wide-field optical telescopes), radar can provide otherwise unavailable information about its size, shape, spin state, and surface properties if it approaches within the range of the Goldstone (California) or Arecibo (Puerto Rico) radar telescopes. Asteroids generally appear as unresolved points through groundbased optical telescopes, but radar measurements of the distribution of echo power in time delay (range) and Doppler frequency (radial velocity) can yield images with resolution as fine as a decameter. Image sequences that furnish adequate orientational coverage can be inverted to construct geologically detailed 3-D models, to define the rotation state precisely, and to constrain the object's internal density distribution. Estimates of radar scattering properties characterize the surface's cm-to-m-scale roughness as well as its bulk density, which for asteroids depends primarily on porosity and metal concentration. A useful spin-off of radar detection is orbit refinement that simplifies navigation of flyby and rendezvous spacecraft.

Ostro, S. J.↗

James Webb Space Telescope Deployment Brushless DC Motor Characteristics Analysis

A DC motor's performance is usually characterized by a series of tests, which are conducted by pass/fail criteria. In most cases, these tests are adequate to address the performance characteristics under environmental and loading effects with some uncertainties and decent power/torque margins. However, if the motor performance requirement is very stringent, a better understanding of the motor characteristics is required. The purpose of this paper is to establish a standard way to extract the torque components of the brushless motor and gear box characteristics of a high gear ratio geared motor from the composite geared motor testing and motor parameter measurement. These torque components include motor magnetic detent torque, Coulomb torque, viscous torque, windage torque, and gear tooth sliding torque. The Aerospace Corp bearing torque model and MPB torque models are used to predict the Coulomb torque of the motor rotor bearings and to model the viscous components. Gear tooth sliding friction torque is derived from the dynamo geared motor test data. With these torque data, the geared motor mechanical efficiency can be estimated and provide the overall performance of the geared motor versus several motor operating parameters such as speed, temperature, applied current, and transmitted power.

Tran, Ahn N.↗

Methods for estimating X-ray machine output through measurement and simulation

We report ball grid arrays are increasingly being applied in the electronics industry and may require X-ray inspection to ensure the integrity and correct placement of solder pins. However, as the architecture of integrated circuits continues to narrow while simultaneously growing more complex, the risk of electronic failure due to radiation damage increases. While medical X-ray devices have been held to high standards and are repeatedly shown to be well characterized, devices used for electronic inspection are often lacking detailed characterization. This study presents unique methods to solve for important properties in X-ray inspection devices such as source to object distance and energy spectrum. This information can then be applied to Monte Carlo models to achieve better overall dose estimates to electronics, which will lead to superior manufactured products. Since X-ray devices can vary greatly in source characteristics, this work investigates spectral measurement and Monte Carlo representation of three X-ray devices. For a Philips SRO 33 100 medical diagnostic device, the spectral output followed expected trends given by the prediction software SpekCalc and Spektr. For the Dage XD7500NT, direct measurement showed a spectral artifact that through the use of Gafchromic films, was shown to be a contributing effect in the dose output. For the Rad Source RS1800, a high powered irradiation device, direct spectral measurement was not achieved. However, a Monte Carlo model using an assumed spectra was found to match ion chamber measurements to a high degree.

47 OTHER INSTRUMENTATION↗

Two-Fluid and Discrete Element Modeling of a Parallel Plate Fluidized Bed Heat Exchanger for Concentrating Solar Power

A novel high-temperature particle solar receiver is developed using a light trapping planar cavity configuration. As particles fall through the cavity, the concentrated solar radiation warms the boundaries of the receiver and in turn heats the particles. Particles flow through the system, forming a fluidized bed at the lower section, leaving the system from the bottom at a constant flowrate. Air is introduced to the system as the fluidizing medium to improve particle heat transfer and mixing. A laboratory scale cavity receiver is built by collaborators at the Colorado School of Mines and their data are used for model validation. In this experimental setup, near IR quartz lamp is used to provide flux to the vertical wall of the heat exchanger. The system is modeled using the discrete element method and a continuum two-fluid method. The computational model matches the experimental system size and the particle size distribution is assumed monodisperse. A new continuum conduction model that accounts for the effects of solid concentration is implemented, and the heat flux boundary condition matches the experimental setup. Radiative heat transfer is estimated using a widely used correlation during the post-processing step to determine an overall heat transfer coefficient. The model is validated against testing data and achieves less than 30% discrepancy and a heat transfer coefficient greater than 1000 W/m2 K.

concentrating solar power↗

Supporting Energy-Related Societal Applications Using NASA's Satellite and Modeling Data

Improvements to NASA Surface Meteorology and Solar Energy (SSE) web site are now being made through the Prediction of Worldwide Energy Resource (POWER) project under NASA Science Mission Directorate Applied Science Energy Management Program. The purpose of this project is to tailor NASA Science Mission results for energy sector applications and decision support systems. The current status of SSE and research towards upgrading estimates of total, direct and diffuse solar irradiance from NASA satellite measurements and analysis are discussed. Part of this work involves collaborating with partners such as the National Renewable Energy Laboratory (NREL) and the Natural Resources Canada (NRCan). Energy Management and POWER plans including historic, near-term and forecast datasets are also overviewed.

Stackhouse, Paul W., Jr.↗

Quantifying agricultural productive use of energy load in Sub-Saharan Africa and its impact on microgrid configurations and costs

The use of advanced energy technologies for agricultural purposes—such as irrigation, refrigeration, crop processing, and egg incubation—has the potential to increase crop yield, reduce vulnerability to changing precipitation patterns, increase shelf life, strengthen income and employment opportunities in rural areas, and reduce emissions by displacing fossil fuel-based technologies. These productive uses of energy (PUE) in remote areas could potentially be powered by microgrids that additionally serve otherwise unelectrified communities, most of which are located in rural Sub-Saharan Africa. Here, in this paper, we use high-resolution geospatial data to estimate the end-use electricity demand for a range of agricultural PUE across Sub-Saharan Africa, and we share these data in an open-access mapping tool. Next, we use REopt®, a techno-economic optimization model of energy systems, to determine the cost and system sizing implications of incorporating agricultural PUE into microgrid designs in Kenya and Zambia. We estimate the upper bound of agricultural PUE demand for irrigation, milling, shelling, refrigeration, and egg incubation across Sub-Saharan to be 16.8 TWh/yr. We find that incorporating local agricultural PUE into microgrid system designs increases the required system sizing while having minimal impact on the levelized cost of energy of these systems. Our analysis is the first to demonstrate the PUE potential in the agricultural sector at a 10x10-kilometer resolution across Sub-Saharan Africa and to show, at scale, how site-specific PUE can impact the cost and sizing of microgrids that are otherwise deployed to serve local household and community load.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Valuing EV Managed Charging for Bulk Power Systems

When and where electric vehicle (EV) charging occurs has significant implications for power systems supporting widespread EV adoption, especially with high shares of wind and solar generation. This study extends previous works by leveraging detailed simulation models for EV adoption, EV use, EV charging, and bulk power system operations, and by linking them with methods for describing charging flexibility at both the individual vehicle and aggregate levels. This technical potential study focuses on how the value of EV managed charging (EVMC) changes depending on charging flexibility type (within-charging session or within-week scheduling), dispatch mechanism (direct load control or one of several price-based mechanisms), and managed charging participation rate. We show that naively aggregating EV charging flexibility from individual vehicles into megawatt-scale resources grossly overestimates the flexibility of the fleet, because such aggregate models can unrealistically pair, e.g., one already-fully-charged vehicle's ability to increase load with another already-charging vehicle's ability to accept more charge, effectively requesting a charging rate that is infeasible for the latter vehicle. We find per-vehicle bulk system value is highest at low participation rates for all dispatch mechanisms. Factoring in production cost savings, avoided firm capacity savings, and combustion-related power sector emissions savings, we estimate the value of EVMC at low participation rates (5%) to be $33/vehicle-year to $69/vehicle-yr for within-session charging flexibility and $40/vehicle-yr to $120/vehicle-yr for within-week charging flexibility in an envisioned 2038 New England power system and monetary value reported in 2016 U.S. dollars. At 100% participation, per-vehicle value declines to $25/vehicle-yr to $31/vehicle-yr for within-session charging flexibility and to $29/vehicle-yr to $36/vehicle-yr for within-week charging flexibility; however, 100% participation yields the highest total system savings.

ADVANCED PROPULSION SYSTEMS↗

Estimates of the wavenumber wavelet power spectrum of magnetic fluctuations during magnetic reconnection

Fluctuation analyses of experimental observations generally lack high temporal resolution and are in frequency-space f, contrary to theoretical efforts in wavenumber-space k. This is due to the inherent limits of the Fourier transform, though it is prominent due to the ease of diagnostic implementation. Advances in wavelet-based analysis have provided relief due to its temporal resolution, but in its common use, is still hard to compare to theoretical models. By using the two-point correlation technique in conjunction with large data sets, a wavelet power spectrum in wavenumber-space can be created. Dubbed the wavenumber wavelet power spectrum, this spectrum relates wavenumber to power in time. Further, this analysis technique more closely connects characterizations of experimentally observed fluctuations with other system parameters and theoretical predictions. In this article, we develop the wavenumber wavelet power spectrum using magnetic fluctuations caused by tearing instability driven magnetic reconnection in reproducible, high temperature laboratory plasmas. These dynamic magnetic fluctuations generated in reversed field pinch plasmas are broadband, ranging from the low frequency, 10's of kHz, up to the ion gyroradii frequencies, 100's of kHz. The dominant fluctuations have poloidal and toroidal mode numbers (m,n)=(1,6−10) and can grow to 2%–3% of the mean magnetic field. During these reconnection events, ions, and electrons are energized, magnetic fluctuation amplitudes increase, plasma flow is halted, and the toroidal magnetic flux increases, all on a semi-periodic basis. The newly developed spectrum provides better temporal resolution of spectrum characteristics to correlate with these particle energization phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

U.S. Agrivoltaics Irradiance Database

This is a foundational data set for research and deployment of agrivoltaics, which is the co-location of agriculture and solar power plants on the same land. This irradiance and shading dataset can be utilized to determine the suitability of agrivoltaics configurations for a given region and crop-type. The data is hourly, 4x4 km resolution across the contiguous United States and Hawaii. It is calculated from the National Solar Radiation Database sites, using the System Advisor Model (SAM) to simulate the shading patterns for 10 common agrivoltaics configurations. Sunlight availability data is reported for 10 locations on the ground between adjacent rows of solar panels, as well as averaged across areas of interest such as the average irradiance in the edge-to-edge open area or across 3-6 planting beds. Other available metrics include the input meteorological data from the NSRDB (e.g. global horizontal irradiance, wind speed, etc.) and estimates for comparing energy and agricultural characteristic across the 10 configurations, including power output per acre or per kW installed capacity and farmable land area per acre.

14 SOLAR ENERGY↗

Machine learning and process-based modeling of spatiotemporal changes in active layer thickness across Alaska

Permafrost degradation poses a growing threat to infrastructure stability and ecosystem resilience in the rapidly warming Arctic. We investigated the spatiotemporal dynamics of active layer thickness (ALT) across Alaska by integrating field observations, environmental datasets, a physically based Stefan model, and machine learning (ML) techniques. Using weather projections from the Coupled Model Intercomparison Project Phase 6 under two Shared Socioeconomic Pathways (SSP 2-4.5 and SSP 5-8.5), we assessed ALT sensitivity to projected future weather conditions. The random forest (RF) model outperformed the Stefan approach in predicting ALT on the training dataset (R² = 0.84 vs. 0.53) but demonstrated lower generalizability on the test dataset (R² = 0.24 vs. 0.54). The root mean square error (RMSE) for the RF model for training and testing ranged from 14 to 22 cm, compared to 17 and 18 cm for the Stefan model. Variable importance analysis revealed that mean annual temperature and slope angle were the strongest predictors of ALT, accounting for 19% and 18% of the variance, respectively, followed by sediment transport index (14%) and stream power index (11%). Comparative analysis of baseline ALT predictions showed the Stefan model tended to project a thicker active layer (mean ± SD: 65 ± 16 cm), compared to the RF model (mean ± SD: 59 ± 8.8) cm). Both models indicated a latitudinal gradient in ALT, with shallower depths at higher latitudes. Projected ALT increases by 2100 were estimated at 3.3 ± 2.2 cm under SSP 2-4.5 and 5.9 ± 4.0 cm under SSP 5-8.5 for the ML model, whereas the Stefan model projected substantially larger increases of 13 ± 2.6 cm (SSP 2-4.5) and 28 ± 4.4 cm (SSP5-8.5). Spatial analysis showed the greatest ALT increases in northern Alaska, with relatively smaller changes in southern regions. These findings highlight the complex, multifactorial nature of ALT dynamics and the value of hybrid modeling approaches. As rising temperatures accelerate permafrost thaw, changes in ALT can disrupt ecosystems, damage infrastructures, and enhance the release of stored soil carbon, highlighting the urgent need for improved predictive capabilities to inform adaptation strategies in the Arctic.

Climate sciences↗