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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 973 records · Page 54

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Optimization using pathwise algorithmic derivatives of electromagnetic shower simulations

Among the well-known methods to approximate derivatives of expectancies computed by Monte-Carlo simulations, averages of pathwise derivatives are often the easiest one to apply. Computing them via algorithmic differentiation typically does not require major manual analysis and rewriting of the code, even for very complex programs like simulations of particle-detector interactions in high-energy physics. However, the pathwise derivative estimator can be biased if there are discontinuities in the program, which may diminish its value for applications. This work integrates algorithmic differentiation into the electromagnetic shower simulation code HepEmShow based on G4HepEm, allowing us to study how well pathwise derivatives approximate derivatives of energy depositions in a sampling calorimeter with respect to parameters of the beam and geometry. We found that when multiple scattering is disabled in the simulation, means of pathwise derivatives converge quickly to their expected values, and these are close to the actual derivatives of the energy deposition. Additionally, we demonstrate the applicability of this novel gradient estimator for stochastic gradient-based optimization in a model example.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Risk-Averse Control of Fractional Diffusion with Uncertain Exponent

In this paper, we introduce and analyze a new class of optimal control problems constrained by elliptic equations with uncertain fractional exponents. We utilize risk measures to formulate the resulting optimization problem. We develop a functional analytic framework, study the existence of solution and rigorously derive the first-order optimality conditions. Additionally, we employ a sample-based approximation for the uncertain exponent and the finite element method to discretize in space. Further, we prove the rate of convergence for the optimal risk neutral controls when using quadrature approximation for the uncertain exponent and conclude with illustrative examples.

finite element method↗

NASA NTRS2021

NASA High Efficiency, High OPR Capable Small Core Compressor

As future aircraft become lighter and more aerodynamically efficient, thrust requirements will decrease, reducing the core size of the engine. Furthermore, in the pursuit of improved fuel burn, engine overall pressure ratio and bypass ratio will increase, further driving down engine core size (core size being defined as high-pressure compressor [HPC] exit-corrected flow). These drivers together mean that the core size for future single-aisle aircraft applications will shrink below 3.0 lb/s. Traditionally, this small core compressor size is in the domain of axi-centrifugal designs, machines that are typically less efficient and limited to pressure ratios of ~25 due to stress and thermomechanical fatigue in the centrifugal impeller. In this light, NASA and Pratt & Whtiney (P&W) embarked upon a program to develop technologies to enable an all-axial high-pressure compressor with a core size below 3.0 lb/s and an overall pressure ratio greater than 50. The challenge with an all-axial high-pressure compressor at this core size is the small span at the rear of the compressor. As core size is scaled down, the rotor tip clearances, stator hub seal clearances, fillet sizes and leading edge thicknesses do not scale, leading to significant efficiency penalties. The goal of this program is to recover this lapse and realize the cycle benefits of small core size and high overall pressure ratio. The small core challenges described are mitigated through design optimization and technology insertion, enabling an estimated 5 to 10% fuel burn reduction relative to 2020 best-in-class. Three test rigs run at NASA Glenn Research Center evolved the small core design: a low-speed rig to vet technology and validate tools, and two high speeds rigs, the first to demonstrate an optimized meanline design and the second to validate technology to manage large rotor tip gaps. The efficiency improvement validated with these rigs has unlocked the small core design space, demonstrating that small core compressors can maintain a similar efficiency to current best-in-class large core size compressors. In addition to advancing the state-of-the-art of technology, the program has also advanced the modeling standards for multistage compressors with large clearance-to-span ratios. A best practice modeling standard was developed over the course of the program, incorporating learning from all three rig programs.

Axial Compressors↗

NASA NTRS2001

Error Estimation of An Ensemble Statistical Seasonal Precipitation Prediction Model

This NASA Technical Memorandum describes an optimal ensemble canonical correlation forecasting model for seasonal precipitation. Each individual forecast is based on the canonical correlation analysis (CCA) in the spectral spaces whose bases are empirical orthogonal functions (EOF). The optimal weights in the ensemble forecasting crucially depend on the mean square error of each individual forecast. An estimate of the mean square error of a CCA prediction is made also using the spectral method. The error is decomposed onto EOFs of the predictand and decreases linearly according to the correlation between the predictor and predictand. Since new CCA scheme is derived for continuous fields of predictor and predictand, an area-factor is automatically included. Thus our model is an improvement of the spectral CCA scheme of Barnett and Preisendorfer. The improvements include (1) the use of area-factor, (2) the estimation of prediction error, and (3) the optimal ensemble of multiple forecasts. The new CCA model is applied to the seasonal forecasting of the United States (US) precipitation field. The predictor is the sea surface temperature (SST). The US Climate Prediction Center's reconstructed SST is used as the predictor's historical data. The US National Center for Environmental Prediction's optimally interpolated precipitation (1951-2000) is used as the predictand's historical data. Our forecast experiments show that the new ensemble canonical correlation scheme renders a reasonable forecasting skill. For example, when using September-October-November SST to predict the next season December-January-February precipitation, the spatial pattern correlation between the observed and predicted are positive in 46 years among the 50 years of experiments. The positive correlations are close to or greater than 0.4 in 29 years, which indicates excellent performance of the forecasting model. The forecasting skill can be further enhanced when several predictors are used.

Shen, Samuel S. P.↗

DOE OSTI2026

Intrusive Uncertainty Quantification and Optimal Experiment Design in the Open-Source Pyomo Ecosystem

This contribution describes ParmEst and Pyomo.DoE, two pillars of the open-source Python-based Pyomo ecosystem for computational optimization with (partial differential) algebraic equation mathematical models. Specifically, ParmEst facilitates intrusive frequentist parameter estimation (PE) and uncertainty quantification (UQ) through built-in features, such as covariance matrix estimation, bootstrapping, and likelihood ratio tests. Complementary, Pyomo.DoE enables optimal experiment design by maximizing various metrics of the Fisher information matrix, such as A-optimality (trace), D-optimality (determinant), E-optimality (minimum eigenvalue), and ME-optimality (condition number). ParmEst and Pyomo.DoE can solve high-dimensional optimization problems by leveraging the model structure and exact derivative information. Finally, we will discuss future opportunities to integrate PE and UQ capabilities with optimization under uncertainty, including robust optimization with non-convex models via PyROS.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1983

Optimal design and use of retry in fault tolerant real-time computer systems

A new method to determin an optimal retry policy and for use in retry of fault characterization is presented. An optimal retry policy for a given fault characteristic, which determines the maximum allowable retry durations to minimize the total task completion time was derived. The combined fault characterization and retry decision, in which the characteristics of fault are estimated simultaneously with the determination of the optimal retry policy were carried out. Two solution approaches were developed, one based on the point estimation and the other on the Bayes sequential decision. The maximum likelihood estimators are used for the first approach, and the backward induction for testing hypotheses in the second approach. Numerical examples in which all the durations associated with faults have monotone hazard functions, e.g., exponential, Weibull and gamma distributions are presented. These are standard distributions commonly used for modeling analysis and faults.

Lee, Y. H.↗

NASA NTRS1962

A STATISTICAL OPTIMIZING NAVIGATION PROCEDURE FOR SPACE FLIGHT

In a typical self-contained space navigation system, celestial observation data are gathered and processed to produce estimated velocity corrections. The results of this paper provide a basis for determining the best celestial measurements and the proper time s to implement velocity corrections. Fundamental to the navigation system is a procedure for processing celestial measurement data which permits incorporation of each individual measurement as it is made in order to provide an improved estimate of position and velocity. In order to "optimize" the navigation, a statistical evaluation of a number of alternative courses of action is made. The various alternatives, which form the basis of a decision process, concern the following: 1) which star and planet combination provides the "best" available observation; 2) whether the best observation gives a sufficient reduction in the predicted target error to warrant making the measurement; and 3) whether the uncertainty in the indicated velocity correction is a small enough percentage of the correction itself to justify an engine re-start and propellant expenditure. Numerical results are presented which illustrate the effectiveness of this approach to the space navigation problem.

NAVIGATION↗

NASA NTRS1973

Continued development of the Combined Pulsed Neutron Experiment (CPNE) for lunar and planetary surfaces

Current progress is reported on the inelastic scattering, capture, and activation gamma ray portions of the Combined Pulsed Neutron Experiment (CPNE). Experiments are described which have enabled a reduction in weight of the experimental probe to 7.3 kg. Parametric studies are described which enabled the optimization of experimental parameters (e.g., gate time settings, neutron pulse rate, etc.). Estimated detection sensitivities using this light weight probe and the optimized experimental parameters are discussed.

Mandler, J. W.↗

DOE OSTI2022

Continuous-Energy ENDF/B-VIII.0 Cross Section and SCALE 6.2.4 Performance for Nuclear Criticality Safety Applications: 1 H, C, 58,60 Ni, 182,183,184,186 W, 235,238 U, 239 Pu

As part of the nuclear data evaluation and validation cycle, the ENDF/B-VIII.0 cross-section library released in 2018 requires testing to determine areas of improvement and deterioration. Previous work by the authors investigated the performance of 16 O, 56 Fe, and 63,65 Cu cross sections, with this study acting as an extension of the prior work. In addition to the isotopes and nuclear criticality safety benchmarks of interest to the prior work, benchmarks from the International Criticality Safety Benchmark Evaluation Project Handbook were selected for their k eff sensitivity to 1 H, C, 58,60 Ni, 182,183,184,186 W, 235,238U, or 239 Pu cross sections and were modeled in the SCALE code system maintained by Oak Ridge National Laboratory. In total, 253 benchmark configurations were selected for their sensitivities and modeled using SCALE 6.2.4 Criticality Safety Analysis Sequences (CSAS) continuous-energy Monte Carlo k eff calculations. This collection includes and expands upon the 99 benchmarks in the prior work. The AMPX-processed ENDF/B-VIII.0 library was decomposed into individual ENDF/B-VIII.0 datum libraries for each isotope of interest. Doing so allowed for the individual substitution of an ENDF/B-VIII.0 cross section in the place of ENDF/B-VII.1, determining isotope-specific effects of ENDF/B-VIII.0 relative to ENDF/B-VII.1. Full library calculations with entirely ENDF/B-VII.1 data or entirely ENDF/B-VIII.0 data were also executed. As a measure of performance, the average relative deviation was determined as the ratio of the deviation between calculated and experimental keff to the propagated calculational and experimental uncertainty. With calculated full library and isotope-specific ENDF/B-VIII.0 k eff ’s, an optimized combination of data libraries was estimated and confirmed with SCALE calculations. This showed that reverting 239 Pu, 58 Ni, 16 O, and 65 Cu cross sections to ENDF/B-VII.1 resulted in improved performance relative to the full ENDF/B-VIII.0 library. Across all 253 benchmarks, the average relative deviation was 1.29σ for the full ENDF/B-VII.1 library, 1.17σ for the full ENDF/B-VIII.0 library, and 0.97σ for the optimized combination. The reversion of 239 Pu, 58 Ni, 16 O, and 65 Cu cross sections to ENDF/B-VII.1 in the 99 benchmarks of the prior work resulted in further improved experimental agreement compared to the previously reported improvement from 16 O and 65 Cu alone. Therefore, it is suggested that applications with significant sensitivities to 239 Pu, 58 Ni, 16 O, and 65 Cu consider their choice of nuclear data library.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

NASA NTRS1990

Geometric programming prediction of design trends for OMV protective structures

The global optimization trends of protective honeycomb structural designs for spacecraft subject to hypervelocity meteroid and space debris are presented. This nonlinear problem is first formulated for weight minimization of the orbital maneuvering vehicle (OMV) using a generic monomial predictor. Five problem formulations are considered, each dependent on the selection of independent design variables. Each case is optimized by considering the dual geometric programming problem. The dual variables are solved for in terms of the generic estimated exponents of the monomial predictor. The primal variables are then solved for by conversion. Finally, parametric design trends are developed for ranges of the estimated regression parameters. Results specify nonmonotonic relationships for the optimal first and second sheet mass per unit areas in terms of the estimated exponents.

Mog, R. A.↗

NASA NTRS2009

Hard X-Ray Flare Source Sizes Measured with the Ramaty High Energy Solar Spectroscopic Imager

Ramaty High Energy Solar Spectroscopic Imager (RHESSI) observations of 18 double hard X-ray sources seen at energies above 25 keV are analyzed to determine the spatial extent of the most compact structures evident in each case. The following four image reconstruction algorithms were used: Clean, Pixon, and two routines using visibilities maximum entropy and forward fit (VFF). All have been adapted for this study to optimize their ability to provide reliable estimates of the sizes of the more compact sources. The source fluxes, sizes, and morphologies obtained with each method are cross-correlated and the similarities and disagreements are discussed. The full width at half-maximum (FWHM) of the major axes of the sources with assumed elliptical Gaussian shapes are generally well correlated between the four image reconstruction routines and vary between the RHESSI resolution limit of approximately 2" up to approximately 20" with most below 10". The FWHM of the minor axes are generally at or just above the RHESSI limit and hence should be considered as unresolved in most cases. The orientation angles of the elliptical sources are also well correlated. These results suggest that the elongated sources are generally aligned along a flare ribbon with the minor axis perpendicular to the ribbon. This is verified for the one flare in our list with coincident Transition Region and Coronal Explorer (TRACE) images. There is evidence for significant extra flux in many of the flares in addition to the two identified compact sources, thus rendering the VFF assumption of just two Gaussians inadequate. A more realistic approximation in many cases would be of two line sources with unresolved widths. Recommendations are given for optimizing the RHESSI imaging reconstruction process to ensure that the finest possible details of the source morphology become evident and that reliable estimates can be made of the source dimensions.

Dennis, Brian R.↗

DOE OSTI2021

Isopod (inverse Optimization And Design): Moose-based Application For Performing Inverse Optimization On Multiphysics Modeling And Simulation

ISOPOD (InverSe OPtimizatiOn and Design) will be a MOOSE (Multiphysics Object Oriented Simulation Environment)-based framework which will iteratively estimate the parameters by minimizing the difference between simulated and experimental observables. While inverse optimization has long been used for parameter estimation, the unique feature of ISOPOD is rapid solution of new inverse problems, leading to accelerated Research, Development and Demonstration (RD&D) for complex parameter estimation problems. This will be achieved by extending the software structure of MOOSE to inverse optimization, as well as building on the unique automatic differentiation feature of MOOSE.

Munday, LynnB↗

DOE OSTI2025

Construction of the damped Ly⁢𝛼 absorber catalog for DESI DR2 Ly⁢𝛼 BAO

We present the Damped Ly⁢𝛼 Toolkit for automated detection and characterization of damped Ly⁢𝛼 absorbers (DLAs) in quasar spectra. Our method uses quasar spectral templates with and without absorption from intervening DLAs to reconstruct observed quasar forest regions. The best-fitting model determines whether a DLA is present while estimating the redshift and HI column density. With an optimized quality cut on detection significance (Δ⁢𝜒$^{2}_{𝑟}$ >0.03), the technique achieves an estimated 80% purity and 79% completeness when evaluated on simulated spectra with S/N>2 that are free of broad absorption lines (BALs). We provide a catalog containing candidate DLAs from the DLA Toolkit detected in DESI DR1 quasar spectra, of which 21 719 were found in S/N>2 spectra with predicted log 10 ⁡(𝑁 𝙷𝙸 )>20.3 and detection significance Δ⁢𝜒$^{2}_{𝑟}$ >0.03. We compare the Damped Ly⁢𝛼 Toolkit to two alternative DLA finders based on a convolutional neural network and Gaussian process models. We present a strategy for combining these three techniques to produce a high-fidelity DLA catalog from DESI DR2 for the Ly⁢𝛼 forest baryon acoustic oscillation measurement. The combined catalog contains 41 152 candidate DLAs with log 10 ⁡(𝑁 𝙷𝙸 )>20.3 from quasar spectra with S/N>2. We estimate this sample to be approximately 85% pure and 79% complete when BAL quasars are excluded.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRSDate not supplied

Surrogate Optimization for Quantum Circuits

Variational quantum Eigensolvers are touted as a near-term algorithm capable of impacting many applications. However, the potential has yet to be realized with few claims of quantum advantage and high resource estimates mainly due to the need for optimization in the presence of noise. Finding algorithms and methods to improve the convergence is essential to accelerate the capabilities of near-term hardware for VQE or more broad applications of hybrid methods in which optimization is required. To this goal we look to use modern approaches recently developed in circuit simulations and stochastic classical optimization that can be combined in a surrogate optimization approach to classical circuits. Using an approximate state vector simulator, we efficiently calculate an approximate Hessian, fed as an input for a detailed quantum circuit simulator. We demonstrate the capabilities of such an approach with and without sampling noise. We also show that this method outperforms Powell in the presence of quantum circuit shot noise by a factor of 2-4

quantum computing↗

DOE OSTI2025

Impact and mitigation of polarized extragalactic foregrounds on Bayesian cosmic microwave background lensing

Future low-noise cosmic microwave background (CMB) lensing measurements from e.g., CMB-S4 will be polarization dominated, rather than temperature dominated. In this new regime, statistically optimal lensing reconstructions outperform the standard quadratic estimator, but their sensitivity to extragalactic polarized foregrounds has not been quantified. Using realistic simulations of polarized radio and infrared point sources, we show for the first time that optimal Bayesian lensing from a CMB-S4-like experiment is insensitive to the expected level of polarized extragalactic foregrounds after masking, as long as an accurate foreground power spectrum is included in the analysis. For more futuristic experiments where these foregrounds could cause a detectable bias, we propose a new method to jointly fit for lensing and the Poisson foregrounds, generalizing the bias hardening from the standard quadratic estimator to Bayesian lensing.

cosmic microwave background↗

NASA NTRS1989

Maximum likelihood techniques in QELS

A framework for the analysis of Quasi Elastic Light Scattering (QELS) experiments designed to be used in microgravity environment is derived. Example calculations of the type to be used to design the QELS system are given. The framework for the analysis is based on the concepts of parameter estimation typified by Maximum Likelihood Estimation methods. These methods not only serve as the template for parameter estimation algorithms, but can also be used for optimal design of the experiments. Optimal design of experiments is facilitated by the fact that these methods not only give procedures for parameter estimation, but also estimates of the errors associated with the parameter estimation.

Edwards, Robert V.↗

DOE OSTI2025

Surrogate optimization of variational quantum circuits

Variational quantum eigensolvers are touted as a near-term algorithm capable of impacting many applications. However, the potential has not yet been realized, with few claims of quantum advantage and high resource estimates, especially due to the need for optimization in the presence of noise. Finding algorithms and methods to improve convergence is important to accelerate the capabilities of near-term hardware for VQE or more broad applications of hybrid methods in which optimization is required. To this goal, we look to use modern approaches developed in circuit simulations and stochastic classical optimization, which can be combined to form a surrogate optimization approach to quantum circuits. Using an approximate (classical CPU/GPU) state vector simulator as a surrogate model, we efficiently calculate an approximate Hessian, passed as an input for a quantum processing unit or exact circuit simulator. This method will lend itself well to parallelization across quantum processing units. We demonstrate the capabilities of such an approach with and without sampling noise and a proof-of-principle demonstration on a quantum processing unit utilizing 40 qubits.

Gustafson, Erik J. [RIACS, Mtn. View] (ORCID:00000↗

DOE OSTI2022

Designing Harvesting and Hauling Cost Models for Energy Cane Production for Biorefineries

The harvesting and hauling operations of bioenergy feedstock is an important area in biofuel production. Production costs can be minimized by maintaining optimal machinery units for these operations. The objective of this study is to design an optimal harvesting unit for bioenergy refinery and estimate harvesting and hauling costs of energy cane. A biorefinery with the annual capacity of processing twenty-five million imp. gallons of ethanol were considered. Given the efficiency of harvesting, a two-row soldier system was considered. Considering the year-round supply of energy cane to the refinery, the optimal machinery unit was designed, and the combined operation costs were derived. The average estimated ownership, repair, labor and fuel and lubricant costs of biomass harvest unit were calculated to be $\$$0.50, $\$$0.54, $\$$1.78 and $\$$1.51/mt, respectively. The costs distribution generated showed harvesting and hauling costs could range between $\$$5.47–$\$$9.23/mt of energy cane. The methodology and the research output will provide guidelines for investors in designing harvesting and hauling units and estimating costs for different scales of operation.

09 BIOMASS FUELS↗