Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Optimal 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,585 records · Page 88

Estimating Switchgrass Biomass Yield and Lignocellulose Composition from UAV-Based Indices

Innovative methods for estimating commercial-scale switchgrass yields and feedstock quality are essential to optimize harvest logistics and biorefinery efficiency for sustainable aviation fuel production. This study utilized vegetation indices (VIs) derived from multispectral images to predict biomass yield and lignocellulose concentrations of advanced bioenergy-type switchgrass cultivars (“Liberty” and “Independence”) under two N rates (28 and 56 kg N ha –1 ). Field-scale plots were arranged in a randomized complete block design (RCBD) and replicated three times at Urbana, IL. Multispectral images captured during the 2021–2023 growing seasons were used to extract VIs. The results show that linear and exponential models outperformed partial least square and random forest models, with mid-August imagery providing the best predictions for biomass, cellulose, and hemicellulose. The green normalized difference vegetation index (GNDVI) was the best univariate predictor for biomass yield (R 2 = 0.86), while a multivariate combination of the GNDVI and normalized difference red-edge index (NDRE) enhanced prediction accuracy (R 2 = 0.88). Cellulose was best predicted using the NDRE (R 2 = 0.53), whereas hemicellulose prediction was most effective with a multivariate model combining the GNDVI, NDRE, NDVI, and green ratio vegetation index (GRVI) (R 2 = 0.44). These findings demonstrate the potential of UAV-based VIs for the in-season estimation of biomass yield and cellulose concentration.

09 BIOMASS FUELS↗

CalWave's xWave Design for PacWave (Final Technical Report)

CalWave Inc. (CalWave) is developing a wave energy converter (WEC) technology that can generate electricity from ocean waves. CalWave’s design offers a unique approach to wave energy conversion that operates fully submerged and can actively adjust the wave excitation. This capability gives the architecture enhanced survivability in ocean storms without adding significant costs. Prior to this project, CalWave had completed a demonstration of a fully functional WEC system in an open ocean demonstration at nominal 1:5 scale under FOA 1663. The goal of this project was the detailed design, following relevant standards and industry best-practices, of a variant of the xWave WEC technology that can safely and efficiently operate at the DOE’s PacWave South test site for a targeted deployment of up two years. The WEC design and associated review processes proceeded in two distinct project phases: a ‘Preliminary’ and a ‘Final’ design phase. The first phase of the project consisted of the systematic design of the WEC’s key features with regards to appropriate IEC standards. The work resulted in a preliminary design of the xWave hull including structural and Power Take-Off (PTO) load estimates, as well as performance estimates for all ocean conditions the WEC would operate in at PacWave South. Following the first open-water demonstration of CalWave’s small-scale “x1” device under FOA 1663, lessons learned were fed directly into a comprehensive review of the xWave design in the second design phase of this FOA project. CalWave’s work was supported by Sandia National Lab (SNL) and the National Renewable Energy Lab (NREL) on the holistic WEC design, and detailed feedback from specialized partners on hull design, mooring and anchoring specification, and electrical grid interconnection. Optimization of the WEC system was performed using a novel numerical optimization tool developed by Sandia and optimization trends were confirmed via an experimental model scale tank test campaign. Performance estimates for PacWave and a detailed xWave design including integration of all relevant system components were concluded. The mooring design was also concluded in the Final design phase using the most up to date sea floor characterization (CPT) data.

16 TIDAL AND WAVE POWER↗

Refining Lucy Mission Delta-V During Spacecraft Design Using Trajectory Optimization Within High-Fidelity Monte Carlo Maneuver Analysis

Recent advances linking medium-fidelity trajectory optimization and high-fidelity trajectory propagation/maneuver design software with Monte Carlo maneuver analysis and parallel processing enabled realistic statistical delta-V estimation well before launch. Completing this high-confidence, refined statistical maneuver analysis early enabled release of excess delta-V margin for increased dry mass margin for the Lucy Jupiter Trojan flyby mission. By 3.3 years before launch, 16 of 34 TCMs had 1000 re-optimized trajectory design samples, yielding tens of m/s lower 99%-probability delta-V versus targeting maneuvers to one optimal trajectory. One year later, 1000 re-optimized samples of all deterministic maneuvers and subsequent flybys further lowered estimated delta-V.

Trojan↗

An optimization approach for cryogenic distillation-based protium removal systems in magnetic fusion energy fuel cycles

Tritium inventory reduction for fusion power plants is central to the successful adoption of fusion energy. The advent of direct internal recycling in deuterium-tritium fusion fuel cycle design has led to significant reduction in startup tritium inventory estimates for fusion power plants but requires an isotope rebalancing and protium removal (IRPR) system to ensure acceptable isotopic fuel composition. Cryogenic distillation is a potential solution for near-term deployment in an IRPR system due to its demonstrated performance in similar operating regimes. Using protium removal in the fuel cycle as the primary performance metric, this paper presents an optimization methodology for a single-column continuous cryogenic distillation-based IRPR system for a 500 MW fus magnetic fusion device. Distillation column optimization was performed using the CryOgenic Distillation For Isotopic Separation of Hydrogen (CODFISH) code developed at Savannah River National Laboratory. The distillation column design presented maximizes direct recycling of hydrogen isotopes from the fusion chamber exhaust to the fueling system while minimizing IRPR system steady-state tritium inventory. The optimized IRPR distillation column presented achieves direct recycling of 60 % of the hydrogen isotopes in the fusion chamber exhaust with an estimated steady-state system tritium inventory <30 g. In conclusion, the optimized IRPR distillation column operation was then used to estimate the design requirements for a detritation column to treat the IRPR system effluent stream.

CryOgenic distillation↗

Cloud Thickness from Diffusion of Lidar Pulses in Clouds

Measurements of the distribution of reflected light from a laser beam incident on an aqueous suspension of particles or "cloud" with known thickness and particle size distribution are reported. The distribution is referred to as the "cloud radiative Green's function", G. In the diffusion domain, G is sensitive to cloud thickness, allowing that important quantity to be retrieved. The goal of the laboratory simulation is to provide preliminary estimates of sensitivity of G to cloud thickness,for use in the optimal design of an offbeam Lidar instrument for remote sensing of cloud thickness (THOR, Thickness from Offbeam Returns). These clouds of polystyrene microspheres suspended in water are analogous to real clouds of water droplets suspended in air. The microsphere size distribution is roughly lognormal, from 0.5 microns to 25 microns, similar to real clouds. Density of suspended spheres is adjusted so mean-free-path of visible photons is about 10 cm, approximately 1000 times smaller than in real clouds. The light source is a ND:YAG laser at 530 nm. Detectors are flux and photon-counting Photomultiplier Tube (PMTS), with a glass probe for precise positioning. A Labview 5 VI controls positioning, and data acquisition, via an NI Motion Control board connected to a stepper motor driving an Edmund linear slider, and a 16-channel 16-bit NI-DAQ board. The stepper motor is accurate to 10 microns, and step size is selectable from the VI software. Far from the incident beam, the rate of exponential increase as the direction of the incident beam is approached scales as expected from diffusion theory, linearly with the cloud thickness, and inversely as the square root of the reduced optical thickness, and is independent of particle size. Near the beam the signal begins to increase faster than exponential, due to single and low-order scattering near the backward direction, and here the distribution depends on particle size. Results are being used to verify 3D Monte Carlo radiative transfer simulations, used to estimate signal-to-noise ratios for remotely sensed off beam returns, for both homogeneous and inhomogeneous clouds. Signal-to-noise estimates show that unfiltered observations are straight forward at night, while narrow band pass filters are being studied for day.

Cahalan, Robert F.↗

Transient Efficiency, Flexibility, and Reliability Optimization of Coal-Fired Power Plants - Final Program Review

This program developed an advanced model-based monitoring and model-predictive control algorithms for a coal fired power plant (CFPP), and deployed these algorithms in a real-time platform to demonstrate performance benefits for transient flexibility and plant operation efficiency. More specifically, the objectives were successfully achieved through a combination of (i) developing a high-fidelity transient plant model in Apros, which was used as a high-fidelity plant simulation between $100-50\% TMCR$ where TMCR denotes the turbine maximum continuous rating, i.e., baseload, (ii) developing a very fast physics-based reduced-order model (ROM) of the plant, which ran more than $100\times$ faster than real-time, enabling its use as real-time embedded model for model-based estimation (MBE) and model predictive control (MPC) (iii) implementing a real-time MBE based on ROM using a robust unscented Kalman filter (UKF) to continuously tune the ROM to match the measurements from high-fidelity Apros plant model despite significant plant-model mismatch, and thus, obtain a Digital Twin of the plant (iv) designing and implementing a real-time MPC with dual objectives of transient plant load tracking with high ramp rates and minimizing coal consumption, i.e., improving plant efficiency in the baseload-partload operation range of $100-50\% TMCR$. Each key element above was developed and tested individually, and has been reported in corresponding Topical Reports. Finally, all the individual elements were integrated in an overall closed-loop system, that was successfully tested in desktop Simulink test harness simulations with ROM or high-fidelity model as the plant. Thereafter, the Simulink implementation was used to auto-generate C-code and deploy as real-time Docker microservice containers in Linux, and validate that they can run in real-time in the hardware-in-the loop (HIL) setup and produce the same results as in Simulink. The results of the integrated simulation tests in Simulink as well as the real-time HIL deployment are documented in this final report, showing good load tracking for load ramps at $3-4\%/min$ ramp rates, and achieving up to $5.5\%$ reduction in coal relative to baseline operation at $50\% TMCR$ load. The desktop and HIL simulations show successful performance of the overall model based estimation and control solution and achieve the key objectives of the program for flexible, efficient and reliable operation of sub-critical coal fired power plants.

20 FOSSIL-FUELED POWER PLANTS↗

Optimal experiment design for identification of large space structures

The optimal experiment design for on-orbit identification of modal frequency and damping parameters in large flexible space structures is discussed. The main result is a separation principle for D-optimal design which states that under certain conditions the sensor placement problem is decoupled from the input design problem. This decoupling effect significantly simplifies the overall optimal experiment design determination for large MIMO structural systems with many unknown modal parameters. The error from using the uncoupled design is estimated in terms of the inherent damping of the structure. A numerical example is given, demonstrating the usefulness of the simplified criteria in determining optimal designs for on-orbit Space Station identification experiments.

Bayard, D. S.↗

A Machine Learning-Based Method to Estimate Transformer Primary-Side Voltages with Limited Customer-Side AMI Measurements

Distribution control applications such as volt/var optimization, network reconfiguration, and distribution automation require accurate knowledge of the distribution system state. The lack of sufficient sensors on the primary side of distribution networks often limits the accuracy of the control decisions by these applications. The deployment of advanced metering infrastructure (AMI) provides utilities an opportunity to translate the AMI data on the secondary onto the primary so that it can be used as pseudo-measurements to augment the limited existing measurements on the primary. This paper develops an approach for estimating service transformer primary-side voltages by using limited secondary-side AMI measurements. The estimated primary-side voltages can be used by utilities as pseudo-measurements for distribution control applications. The detailed secondary model topology, which is an essential input data for many existing algorithms, is not required for the proposed method. The performance of the proposed method is validated by using AMI measurements from the field and an actual distribution feeder model of San Diego Gas & Electric Company.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Locally adaptive artificial viscosity strategies for Lagrangian hydrodynamics

To accurately model inviscid flow with shock waves using staggered-grid Lagrangian hydrodynamics, artificial viscosity is introduced to convert kinetic energy into internal energy, thereby providing a mechanism to generate the required entropy increase across shocks. In this paper, we propose a new method for constructing an adaptive, artificial viscosity in the context of one-dimensional, staggered-grid Lagrangian hydrodynamics. Our adaptive, artificial viscosity is defined in terms of two parameters that depend on density locally in a neighborhood around a shock, and hence, vary cell-by-cell. Our methodology is based on building a reference set of pre-computed optimal, globally constant artificial viscosity coefficients for a family of isolated shock test problems. For arbitrary flows, the evaluation of the unknown coefficients is automated by first estimating shock intensity locally, and second computing the corresponding adaptive parameter value via interpolation over the data from the pre-computed reference set of optimal, globally constant coefficient values. To illustrate the performance of our new approach, we compare our results against two existing methods. The first method is a limiter-based approach, which relies on estimating velocity gradients of the flow, and the second method is utilized in several commercial codes. Finally, we demonstrate that our new adaptive methodology produces more accurate results for a variety of tests with propagating shock waves, as well as for the aforementioned family of isolated shock problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On-line estimation of error covariance parameters for atmospheric data assimilation

A simple scheme is presented for on-line estimation of covariance parameters in statistical data assimilation systems. The scheme is based on a maximum-likelihood approach in which estimates are produced on the basis of a single batch of simultaneous observations. Simple-sample covariance estimation is reasonable as long as the number of available observations exceeds the number of tunable parameters by two or three orders of magnitude. Not much is known at present about model error associated with actual forecast systems. Our scheme can be used to estimate some important statistical model error parameters such as regionally averaged variances or characteristic correlation length scales. The advantage of the single-sample approach is that it does not rely on any assumptions about the temporal behavior of the covariance parameters: time-dependent parameter estimates can be continuously adjusted on the basis of current observations. This is of practical importance since it is likely to be the case that both model error and observation error strongly depend on the actual state of the atmosphere. The single-sample estimation scheme can be incorporated into any four-dimensional statistical data assimilation system that involves explicit calculation of forecast error covariances, including optimal interpolation (OI) and the simplified Kalman filter (SKF). The computational cost of the scheme is high but not prohibitive; on-line estimation of one or two covariance parameters in each analysis box of an operational bozed-OI system is currently feasible. A number of numerical experiments performed with an adaptive SKF and an adaptive version of OI, using a linear two-dimensional shallow-water model and artificially generated model error are described. The performance of the nonadaptive versions of these methods turns out to depend rather strongly on correct specification of model error parameters. These parameters are estimated under a variety of conditions, including uniformly distributed model error and time-dependent model error statistics.

Dee, Dick P.↗

A Machine Learning-Based Method to Estimate Transformer Primary-Side Voltages with Limited Customer-Side AMI Measurements

Distribution control applications such as volt/var optimization, network reconfiguration, and distribution automation require accurate knowledge of the distribution system state. The lack of sufficient sensors on the primary side of distribution networks often limits the accuracy of the control decisions by these applications. The deployment of advanced metering infrastructure (AMI) provides utilities an opportunity to translate the AMI data on the secondary onto the primary so that it can be used as pseudo-measurements to augment the limited existing measurements on the primary. This paper develops a machine learning based approach for estimating service transformer primary-side voltages by using limited secondary-side AMI measurement. The machine learning model is developed by using random forest algorithm. The estimated primary-side voltages can be used by utilities as pseudo-measurements for distribution control applications. The detailed secondary model topology, which is an essential input data for many existing algorithms, is not required for the proposed method. The performance of the proposed method is validated by using AMI measurements from the field and an actual distribution feeder model of San Diego Gas & Electric Company.

advanced metering infrastructure↗

Information content of spatially distributed ground-based measurements for hydrologic-parameter calibration in mixed rain-snow mountain headwaters

Parameters in hydrologic models used in mixed rain-snow regions are often uncertain to calibrate and overfitted on streamflow. To contribute addressing these challenges, we used an algorithm that assesses modeling performances through time (Dynamic Identifiability Analysis) to quantify the information content of spatially distributed ground-based measurements for identifying optimal parameter values in the Precipitation Runoff Modeling System (PRMS) model. Including spatially distributed ground-based measurements in Identifiability Analysis allowed us to unambiguously estimate more parameter values than only using streamflow (seven parameters instead of two out of a pool of thirty-three). Peaks in information gain were obtained when using dew-point temperature to identify precipitation phase-partitioning parameters. Multi-attribute identifiability analysis also yielded optimal parameter values that were temporally less variable than those estimated using streamflow alone. Overall, identifying parameter values using ground-based measurements improved the simulation of key drivers of the surface-water budget, such as air temperature and precipitation-phase partitioning. However, parameters simulating surface-to-subsurface mass fluxes like snow accumulation and melt or evapotranspiration were poorly identified by any attribute and so emerged as key sources of predictive uncertainty for this distributed-parameter hydrologic model. In conclusion, this work demonstrates the value of expanded ground-based measurements for identifying parameters in distributed-parameter hydrologic models and so diagnosing their conceptual uncertainty across the water budget.

54 ENVIRONMENTAL SCIENCES↗

Optimization performance, fidelity, and cost: SIAM VQE

This dataset contains files storing results from classically-simulated quantum subroutines within a dynamical mean-field theory workflow, and jupyter notebooks processing the data in these files to generate plots. The files store: (1) Results from variational quantum eigensolver (VQE) simulations searching for optimal parameters allowing parametrized quantum circuits to prepare approximations to ground states of different Anderson impurity models (AIMs) (2) Results from simulations of a quantum Lanczos algorithm (QLA) estimating the Lanczos coefficients defining the continued-fraction representation of an (AIM) Green’s function Description: Any file named vqe_gs_results* stores approximations to the ground state and energy of a given AIM estimated using three different methods: (1) Numerical diagonalization (2) Ideal VQE simulation (3) VQE simulation with sampling noise For each VQE simulations metadata about the optimization (optimization results plus number of quantum circuits that would have been executed on real hardware) is also stored. Any file named qla_dos_results* estimations for the Lanczos coefficients defining the Green’s function of an AIM. The stored estimations are achieved using different methods: (1) Numerical Lanczos algorithm from initial states obtained from numerical diagonalization (2) Simulated quantum Lanczos algorithm from initial states prepared from parametrized quantum circuits yielded by corresponding ideal and noisy VQE subroutines. The dataset is used and described in M. Karabin et al., "Quantum solver for single-impurity Anderson models with particle-hole symmetry", Phys. Rev. Research 8, 033066 (2026). DOI: https://doi.org/10.1103/7ys3-tl4l

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Global Radius of Curvature Estimation and Control for the Hobby-Eberly Telescope

A system, which estimates the global radius of curvature (GroC) and corrects for changes in GroC on a segmented primary mirror has been developed for and verified on McDonald Observatory's Hobby Eberly Telescope (HET). The GroC estimation and control system utilizes HET's primary mirror control (PMC) system and the Segment Alignment Maintenance System (SAMS), an inductive edge sensor system. A special set of boundary conditions is applied to the derivation of the optimal edge match control. The special boundary conditions allow the further derivation of an observer, which enables estimation and control of the Groc mode to within HET's specification. The magnitude of the GroC mode can then be controlled despite the inability of the SAMS edge sensor system, by itself, to observe or control the GroC mode. The observer can be extended to any segmented mirror telescope. It will be shown that the observer improves with accuracy as the number of segments increases. This paper presents the mathematical theory of the observer. Simulation results will demonstrate the inherent accuracy and robustness of the system. Performance verification data from the HET will be presented.

Rakoczy, John↗

A Dynamic Amplitude-Correcting Gradient Estimation Technique to Align X-ray Focusing Optics

High-brightness X-rays, as produced at synchrotrons and X-ray free electron laser (XFEL) facilities, are used to characterize materials in a variety of scientific experiments. In most cases, effective use of the high-energy light requires precisely-aligned focusing optics; one example being a compound refractive lens (CRL). To align a CRL, the position and rotation must be optimized along four axes. In practice, this is a labor-intensive, time-consuming manual process that can monopolize scarce experimental time at the necessary X-ray facilities. Models of the expected Xray transmission function suggest that this task can be automated; however, the temporally-varying intensity at X-ray free electron laser facilities preclude the direct use of standard implementations of optimization solvers such as steepest descent algorithms. In this paper, we propose a novel technique to estimate the gradient of noisy functions with temporally-varying amplitudes. We construct this dynamicamplitude correction by systematically sampling a fixed central location within the standard finite difference stencil, accounting for the observed changes in time, and normalizing the difference quotients against those fluctuations. In addition to a rigorous error analysis of the sampling technique, we demonstrate its efficacy in stochastic descent optimization methods. Further, we demonstrate how this approach may be implemented to optimize X-ray focusing optics at synchrotrons or XFEL facilities

97 MATHEMATICS AND COMPUTING↗

Thermal Property Estimation of Fibrous Insulation: Heat Transfer Modeling and the Continuous Genetic Algorithm

Thermal properties of high-temperature fibrous insulation materials were estimated from transient thermal tests using an inverse heat transfer technique. Transient temperature data from an experimental set up was collected and a simple, one-dimensional numerical model was constructed to replicate the temperatures within the test assembly. The Continuous Genetic Algorithm optimization technique in conjunction with a numerical thermal model and transient test data was used to estimate coefficients of a functional representation of the thermal property. The thermal properties, i.e., thermal conductivity and specific heat, of an alumina insulation felt were estimated over the temperature range of 300 K to 1700 K at various constant static pressures in nitrogen gas and compared with published data. The resulting thermal property estimates were within 10% of published values over the entire temperature range at various pressures. The methodology, application, and results are presented.

Frye, Elora↗

Adaptive sampling quasi-Newton methods for zeroth-order stochastic optimization

Here, we consider unconstrained stochastic optimization problems with no available gradient information. Such problems arise in settings from derivative-free simulation optimization to reinforcement learning. We propose an adaptive sampling quasi-Newton method where we estimate the gradients using finite differences of stochastic function evaluations within a common random number framework. We develop modified versions of a norm test and an inner product quasi-Newton test to control the sample sizes used in the stochastic approximations and provide global convergence results to the neighborhood of a locally optimal solution. We present numerical experiments on simulation optimization problems to illustrate the performance of the proposed algorithm. When compared with classical zeroth-order stochastic gradient methods, we observe that our strategies of adapting the sample sizes significantly improve performance in terms of the number of stochastic function evaluations required.

97 MATHEMATICS AND COMPUTING↗

Reduced state feedback gain computation

Because application of conventional optimal linear regulator theory to flight controller design requires the capability of measuring and/or estimating the entire state vector, it is of interest to consider procedures for computing controls which are restricted to be linear feedback functions of a lower dimensional output vector and which take into account the presence of measurement noise and process uncertainty. Therefore, a stochastic linear model that was developed is presented which accounts for aircraft parameter and initial uncertainty, measurement noise, turbulence, pilot command and a restricted number of measurable outputs. Optimization with respect to the corresponding output feedback gains was performed for both finite and infinite time performance indices without gradient computation by using Zangwill's modification of a procedure originally proposed by Powell. Results using a seventh order process show the proposed procedures to be very effective.

Kaufman, H.↗