Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “input 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,297 records · Page 72

Phase-locked loop FM demodulator

A conventional phase-locked loop is improved by replacing its phase detector with one comprising a linear ramp generator and a sample-and-hold circuit, thus eliminating the need for a lowpass loop filter, although the output of the sample-and-hold circuit may be filtered in the case of a very low level modulating signal on the incoming FM signal, but then filtering is not a difficult problem as in a conventional phase-locked loop. The result is FM demodulation by zero-order estimation. For FM demodulation by first-order estimation, the arithmetic difference between adjacent samples is formed, and using a second sample-and-hold circuit an arithmetic difference signal is produced as an input to a second ramp generator that is reset after each sampling cycle to generate a ramp the slope of which is a function of the arithmetic difference signal stored in the second sample-and-hold circuit. The ramp thus generated by the second ramp generator is arithmetically summed with the zero-estimation signal from the first sample-and-hold circuit to form a first-order estimation signal. Filtering such a first-order estimation signal is less of a problem than filtering a zero-order estimation signal.

Kirkham, Harold↗

A hybrid calibration approach to Hertz-type contact parameters for discrete element models

This study aims at providing a hybrid calibration framework to estimate Hertz-type contact parameters (particle-scale shear modulus and Poisson ratio) for both two-dimensional and three-dimensional discrete element modelling (DEM). On the basis of statistically isotropic granular packings, a set of analytical formulae between macroscopic material parameters (Young modulus and Poisson ratio) and particle-scale Hertz-type contact parameters for granular systems are derived under small-strain isotropic stress conditions. However, the derived analytical solutions are only estimated values for general models. By viewing each DEM modelling as an implicit mathematical function taking the particle-level parameters as independent variables and employing the derived analytical solutions as the initial input parameters, an automatic iterative scheme is proposed to obtain the calibrated parameters with higher accuracies. Considering highly nonlinear features and discontinuities of the macro-micro relationship in Hertz-based discrete element models, the adaptive moment estimation algorithm is adopted in this study because of its capacity of dealing with noise gradients of cost functions. Here, the proposed method is validated with several numerical cases including randomly distributed monodisperse and polydisperse packings. Noticeable improvements in terms of calibration efficiency and accuracy have been made.

Constitutive law↗

The PAWS and STEM reliability analysis programs

The PAWS and STEM programs are new design/validation tools. These programs provide a flexible, user-friendly, language-based interface for the input of Markov models describing the behavior of fault-tolerant computer systems. These programs produce exact solutions of the probability of system failure and provide a conservative estimate of the number of significant digits in the solution. PAWS uses a Pade approximation as a solution technique; STEM uses a Taylor series as a solution technique. Both programs have the capability to solve numerically stiff models. PAWS and STEM possess complementary properties with regard to their input space; and, an additional strength of these programs is that they accept input compatible with the SURE program. If used in conjunction with SURE, PAWS and STEM provide a powerful suite of programs to analyze the reliability of fault-tolerant computer systems.

Butler, Ricky W.↗

Uncertainty Quantification of Global Net Methane Emissions From Terrestrial Ecosystems Using a Mechanistically Based Biogeochemistry Model

Quantification of methane (CH4) emissions from wetlands and its sinks from uplands is still fraught with large uncertainties. Here, a methane biogeochemistry model was revised, parameterized, and verified for various wetland ecosystems across the globe. The model was then extrapolated to the global scale to quantify the uncertainty induced from four different types of uncertainty sources including parameterization, wetland type distribution, wetland area distribution, and meteorological input. We found that global wetland emissions are 212 ± 62 and 212 ± 32 Tg CH4 year−1 (1Tg = 1012 g) due to uncertain parameters and wetland type distribution, respectively, during 2000–2012. Using two wetland distribution data sets and three sets of climate data, the model simulations indicated that the global wetland emissions range from 186 to 212 CH4 year−1 for the same period. The parameters were the most significant uncertainty source. After combining the global methane consumption in the range of −34 to −46 Tg CH4 year−1, we estimated that the global net land methane emissions are 149–176 Tg CH4 year−1 due to uncertain wetland distribution and meteorological input. Spatially, the northeast United States and Amazon were two hotspots of methane emission, while consumption hotspots were in the Eastern United States and eastern China. During 1950–2016, both wetland emissions and upland consumption increased during El Niño events and decreased during La Niña events. This study highlights the need for more in situ methane flux data, more accurate wetland type, and area distribution information to better constrain the model uncertainty.

wetland methane emission↗

The use of multispectral thermal infrared image data to estimate the sulfur dioxide flux from volcanoes: A case study from Mount Etna, Sicily, July 29, 1986

We have found that image data acquired with NASA's airborne Thermal Infrared Multispectral Scanner (TIMS) can be used to make estimates of the SO2 content of volcanic plumes. TIMS image data are most applicable to the study of partially transparent SO2 plumes, such as those released during quiescent periods or nonexplosive eruptions. The estimation procedure is based on the LOWTRAN 7 radiative transfer code, which we use to model the radiance perceived by TIMS as it views the ground through an SO2 plume. The input to the procedure includes the altitudes of the aircraft and ground, the altitude and thickness of the SO2 plume, the emissivity of the ground, and altitude profiles of the atmospheric pressure, temperature, and relative humidity. We use the TIMS data to estimate both ground temperatures beneath a plume and SO2 concentrations within a plume. Applying our procedure to TIMS data acquired over Mount Etna, Sicily, on July 29, 1986, we estimate that the SO2 flux from the volcano was approximately 6700 t d(exp -1). The use of TIMS to study SO2 plumes represents a bridge between highly localized methods, such as correlation spectroscopy or direct sampling, and small-scale mapping techniques involving satellite instruments such as the Total Ozone Mapping Spectrometer or Microwave Limb Sounder. We require further airborne experiments to refine our estimation procedure. This refinement is a necessary preparation for the schedueled 1998 launch of the Advanced Spaceborne Thermal Emission and Reflectance Radiometer, which will allow large-scale multispectral thermal infrared image data to be collected over virtually any volcano on Earth at least once every 16 days.

Realmuto, Vincent J.↗

Temperature-moisture biases in ECMWF analyses based on clear sky longwave simulations constrained by SSMI and MSU measurements and comparisons to ERBE estimates

Clear sky longwave radiation fluxes for the summer of 1988 and winter of 1989 have been simulated with a radiative transfer model that includes detailed treatment of atmospheric gas absorption. The input data to the model are humidity and temperature profiles from European Center for Medium-Range Weather Forecasts (ECMWF) analyses and surface temperature measurements from the International Satellite Cloud Climatology Project. To reduce the inherent uncertainties in humidity profiles, the ECMWF moisture fields have been adjusted, based on the climatological relationship between the moisture profile and total precipitable water (PW), so that the ECMWF PW is equal to that derived from Special Sensor Microwave Imager (SSMI) data. SSMI and ECMWF PW patterns are generally similar, but significant differences in magnitude are found in the low latitudes, particularly over the subtropical oceans in the Southern Hemisphere. At the same time the tropospheric brightness temperature (T(sub B)) at 53.74 GHz has been computed using a microwave transfer model with ECMWF temperature and humidity fields as inputs. The comparison of computed T(sub B) with the Microwave Sounding Unit (MSU) channel 2 T(sub B) revealed differences larger than 1.5 K over most of the oceans, suggesting that the ECMWF model atmosphere has a cold bias. Maximum biases greater than 3 K are found in low latitudes. After removing the bias from the ECMWF temperature field, a new set of temperature profiles yielding brightness temperature equal to MSU channel 2 T(sub B) has been obtained. Simulation results clearly demonstrate that inclusion of satellite estimates of PW and T(sub B) enhances the accuracy of the clear sky LW flux simulation, substantially reducing the differences from Earth Radiation Budget Experiment (ERBE) estimates. Global averages indicate that the model-derived top-of-atmosphere (TOA) clear sky fluxes from the adjusted moisture and temperature profiles are in very good agreement with satellite-measured ERBE values, differing by only 1.4 W/sq m. By contrast, the calculation using only ECMWF analyses without any adjustments gives a significant underestimation of TOA clear sky fluxes up to 10.7 W/sq m compared to the satellite measurements, suggesting that the constraint method can significantly improve the accuracy in the radiation budget simulation. The results also indicate that the ECMWF-based simulation errors are, at least for the two seasons studied, mainly attributed to the error in the temperature field rather than to errors in the moisture field.

Sohn, Byung-Ju↗

Estimation of distributions via multilevel Monte Carlo with stratified sampling

We design and implement a novel algorithm for computing a multilevel Monte Carlo (MLMC) estimator of the joint cumulative distribution function (CDF) of a vector-valued quantity of interest in problems with random input parameters and initial conditions. Our approach combines MLMC with stratified sampling of the input sample space by replacing standard Monte Carlo at each level with stratified Monte Carlo initialized with proportionally allocated samples. We show that the resulting stratified MLMC (sMLMC) algorithm is more efficient than its standard MLMC counterpart due to the additional variance reduction provided by the stratification of the random parameter's domain, especially at the coarsest levels. Additional computational cost savings are obtained by smoothing the indicator function with a Gaussian kernel, which proves to be an efficient and robust alternative to recently developed polynomial-based techniques.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Estimation of Kalman filter gain from output residuals

This paper presents a procedure for estimating the Kalman filter gain from output residuals. The system state space model is assumed to be known, but the process and noise covariance are unknown. The proposed procedure consists of three basic steps. First, the output residuals are computed from the given model and a given set of input-output data. Second, a linear regression model for this part of the response is computed by a least squares solution. Third, the Kalman filter gain is then estimated from the coefficients of this model. Numerical results using experimental data are presented to illustrate the validity of the developed procedure.

Juang, Jer-Nan↗

Protocol to detect dilution cycles in chemostat experiments and estimate growth rate slopes with linear modeling with R software chemostat_regression

Chemostat growth chambers measure optical density over time and require manual calculation of growth rates. Here, we present chemostat_regression, R software that enables users to automatically identify chemostat cycles and estimate growth rate using a linear regression approach. We describe steps for creating requisite software environment(s), formatting input data, executing the software via command line/RStudio/R-Shiny, interpreting results, assessing the validity of results, and modifying input parameters.

59 BASIC BIOLOGICAL SCIENCES↗

Alpert multi-wavelets for functional inverse problems: direct optimization and deep learning

Computational engineering models often contain unknown entities (e.g. parameters, initial and boundary conditions) that require estimation from other measured observable data. Estimating such unknown entities is challenging when they involve spatio-temporal fields because such functional variables often require an infinite-dimensional representation. Here, we address this problem by transforming an unknown functional field using Alpert wavelet bases and truncating the resulting spectrum. Hence the problem reduces to the estimation of few coefficients that can be performed using common optimization methods. We apply this method on a one-dimensional heat transfer problem where we estimate the heat source field varying in both time and space. The observable data is comprised of temperature measured at several thermocouples in the domain. This latter is composed of either copper or stainless steel. The optimization using our method based on wavelets is able to estimate the heat source with an error between 5% and 7%. We analyze the effect of the domain material and number of thermocouples as well as the sensitivity to the initial guess of the heat source. Finally, we estimate the unknown heat source using a different approach based on deep learning techniques where we consider the input and output of a multi-layer perceptron in wavelet form. We find that this deep learning approach is more accurate than the optimization approach with errors below 4%.

97 MATHEMATICS AND COMPUTING↗

Extended State Observer-Based Robust Model Predictive Velocity Control for Permanent Magnet Synchronous Motor

This article proposes an extended state observer based robust model predictive velocity control to decrease system prediction error under parameter uncertainties for permanent magnet synchronous motor (PMSM). We develop a new PMSM model that consists of velocity and acceleration to lump the system information and an external disturbance into a disturbance. The extended state observer (ESO) is designed to estimate the velocity, acceleration, and disturbance. By estimating the state variables and disturbance using the ESO, the model predictive control (MPC) finds the optimal control input by predicting future system behavior. Additionally, the direct current controller is designed so that the direct current converges to zero. Because the proposed method is not designed based on the cascade structure from the viewpoint of velocity control, the optimization control for the velocity and currents can be defined. Thus, the proposed method is robust against external disturbances and parameter uncertainties owing to feedback linearization, state feedback, and ESO-based MPC using the acceleration PMSM model. The proposed control algorithm was experimentally verified and it showed improved velocity tracking performance compared with ESO-based MPC using the conventional PMSM model.

42 ENGINEERING↗

Hyporheic-zone Processes and Stream Oxygen Dynamics: Insights from a Multiscale Reactive Transport Model: Modeling Archive

This archive contains the data and Python scripts required to reproduce the analyses and figures in the study: Gomez-Velez, J. D., Rathore, S. S., Cohen, M. J., & Painter, S. L. (2025). Hyporheic-zone Processes and Stream Oxygen Dynamics: Insights from a Multiscale Reactive Transport Model. Submitted to Water Resources Research. The analysis utilizes the subgrid model Advection Dispersion Equation with Lagrangian Subgrids (ADELS) implemented in the Advanced Terrestrial Simulator (ATS; https://amanzi.github.io/ats/stable/). In this case, the ATS and Amanzi versions are (1) ATS version 1.5.1_f5ba18f8 and (2) Amanzi version 1.6-dev_53444cca4. The repository includes a Jupyter Notebook and the necessary data (Pandas DataFrames stored as pickle files) to generate the figures for the manuscript. Additionally, it contains Python scripts to create ATS input files, run the ATS simulations, and post-process the results. Finally, it provides routines for parameter estimation using the Single-Station Metabolism (SSM) model with the Differential Evolution Adaptive Metropolis (DREAM) Markov Chain Monte Carlo (MCMC) algorithm with ZS enhancements (DREAM-ZS).

54 ENVIRONMENTAL SCIENCES↗

NUCLEI (SciDAC-4) effort at The University of Tennessee, Knoxville

This Report summarizes the accomplishments of NUCLEI researchers at The University of Tennessee, Knoxville (UTK) and those supported via subcontracts at the University of North Carolina at Chapel Hill (UNC) and Iowa State University (ISU). The SciDAC-4 NUCLEI collaboration had the goal to advance the computation of atomic nuclei and to make reliable predictions including uncertainty estimates for observables of interest. These goals were achieved by developing improved nuclear interactions and currents that are input to the computations, and by advancing the algorithms and techniques used for computing nuclear observables of interest. The computed results advanced our understanding of how nuclei work; they helped to guide and interpret experiments at existing low-energy nuclear physics facilities, the Facility for Rare Isotope Beams, Jefferson Laboratory, and neutrino experiments. In addition, predictions were made that await confrontation with experimental results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Evaluating the 238 U PFNS Including Chi-Nu Experimental Data

This report documents an evaluation of 238 U prompt fission neutron spectra (PFNS) which is a deliverable for a FY2024 Q4 NCSP (Nuclear Criticality Safety Program) milestone. This evaluation is new; its prior input is based on extended Los Alamos and exciton models implemented in the code CoH. Experimental covariances were estimated for five experimental data sets. One of these data sets that was measured by the Chi-Nu team of LANL and LLNL. It covers the 238 U PFNS for continuous incident-neutron energies of 1–20 MeV and outgoing-neutron energies from 10 keV– 10 MeV with high precision. Contrary to Chi-Nu data, previous data sets were measured in a limited energy range. The resulting evaluated data correspond well to the experimental PFNS taken into account for the evaluation. The evaluated PFNS also produce average mean energies in agreement with associated Chi-Nu data. If one uses the new evaluated data to predict the neutron multiplication factor, k eff , of the Flattop, Flattop-Pu and BigTen ICSBEP critical assemblies (which all have thick reflectors with high percentages of 238 U), the differences of simulated values compared to those using ENDF/B-VIII.1β3 is modest (less than 25 pcm). In addition to that, the new PFNS predict on average 238 U LLNL pulsed-sphere neutron-leakage spectra slightly better than ENDF/BVIIII.0 and ENDF/B-VIII.1β3 PFNS. The differences are, however, well within the experimental uncertainties.

238U↗

P-model v1.0: an optimality-based light use efficiency model for simulating ecosystem gross primary production

Terrestrial photosynthesis is the basis for vegetation growth and drives the land carbon cycle. Accurately simulating gross primary production (GPP, ecosystem-level apparent photosynthesis) is key for satellite monitoring and Earth system model predictions under climate change. While robust models exist for describing leaf-level photosynthesis, predictions diverge due to uncertain photosynthetic traits and parameters which vary on multiple spatial and temporal scales. Here, we describe and evaluate a GPP (photosynthesis per unit ground area) model, the P-model, that combines the Farquhar–von Caemmerer–Berry model for C 3 photosynthesis with an optimality principle for the carbon assimilation–transpiration trade-off, and predicts a multi-day average light use efficiency (LUE) for any climate and C 3 vegetation type. The model builds on the theory developed in Prentice et al. (2014) and Wang et al. (2017a) and is extended to include low temperature effects on the intrinsic quantum yield and an empirical soil moisture stress factor. The model is forced with site-level data of the fraction of absorbed photosynthetically active radiation (fAPAR) and meteorological data and is evaluated against GPP estimates from a globally distributed network of ecosystem flux measurements. Although the P-model requires relatively few inputs, the R 2 for predicted versus observed GPP based on the full model setup is 0.75 (8 d mean, 126 sites) – similar to comparable satellite-data-driven GPP models but without predefined vegetation-type-specific parameters. The R 2 is reduced to 0.70 when not accounting for the reduction in quantum yield at low temperatures and effects of low soil moisture on LUE. The R 2 for the P-model-predicted LUE is 0.32 (means by site) and 0.48 (means by vegetation type). Applying this model for global-scale simulations yields a total global GPP of 106–122 Pg C yr –1 (mean of 2001–2011), depending on the fAPAR forcing data. The P-model provides a simple but powerful method for predicting – rather than prescribing – light use efficiency and simulating terrestrial photosynthesis across a wide range of conditions. The model is available as an R package (rpmodel).

54 ENVIRONMENTAL SCIENCES↗

Prediction of overall and blade-element performance for axial-flow pump configurations

A method and a digital computer program for prediction of the distributions of fluid velocity and properties in axial flow pump configurations are described and evaluated. The method uses the blade-element flow model and an iterative numerical solution of the radial equilbrium and continuity conditions. Correlated experimental results are used to generate alternative methods for estimating blade-element turning and loss characteristics. Detailed descriptions of the computer program are included, with example input and typical computed results.

Serovy, G. K.↗

Flight test design for efficient extraction of aircraft parameters

This paper considers the design of flight control inputs which result in a minimization of the nonuniqueness problem of parameter identification due to statistical correlation. The performance index for the control design is taken as a linear weighted sum of the squares of the correlation coefficients of the aerodynamic stability and control derivatives as computed from the Cramer-Rao lower bound matrix. An optimal input design is demonstrated for an F-8 aircraft with supercritical wing. Reduced correlations are noted and corresponding estimates and confidence levels presented.

Wells, W. R.↗

Algorithms for adaptive stochastic control for a class of linear systems

Control of linear, discrete time, stochastic systems with unknown control gain parameters is discussed. Two suboptimal adaptive control schemes are derived: one is based on underestimating future control and the other is based on overestimating future control. Both schemes require little on-line computation and incorporate in their control laws some information on estimation errors. The performance of these laws is studied by Monte Carlo simulations on a computer. Two single input, third order systems are considered, one stable and the other unstable, and the performance of the two adaptive control schemes is compared with that of the scheme based on enforced certainty equivalence and the scheme where the control gain parameters are known.

Toda, M.↗