Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Function 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,099 records · Page 61

Performance of a 2-micrometer coherent Doppler lidar for wind measurements

Measurements of boundary layer winds are presented using a 2-micrometer coherent Doppler lidar and the optimal performance of the maximum likelihood estimator. The systematic error for single-shot estimates was estimated as 3.6 cm/s using measurements from a stationary hard target. The estimation error for measurements of the radial component of the wind field was determined, as well as the fraction of the estimates that are randomly distributed over the velocity search space, when the signal power is low and speckle fading is important. The results from actual data are compared with the results from ideal simulations. The first direct estimation of the spatial structure function of the radial wind field and of the energy dissipation rate is presented for both horizontal and vertical directions of propagation. The rms estimation error of the velocity estimates is found to be within 30% of ideal performance based on simulation.

Frehlich, Rod↗

Polarimetric, Two-Color, Photon-Counting Laser Altimeter Measurements of Forest Canopy Structure

Laser altimeter measurements of forest stands with distinct structures and compositions have been acquired at 532 nm (green) and 1064 nm (near-infrared) wavelengths and parallel and perpendicular polarization states using the Slope Imaging Multi-polarization Photon Counting Lidar (SIMPL). The micropulse, single photon ranging measurement approach employed by SIMPL provides canopy structure measurements with high vertical and spatial resolution. Using a height distribution analysis method adapted from conventional, 1064 nm, full-waveform lidar remote sensing, the sensitivity of two parameters commonly used for above-ground biomass estimation are compared as a function of wavelength. The results for the height of median energy (HOME) and canopy cover are for the most part very similar, indicating biomass estimations using lidars operating at green and near-infrared wavelengths will yield comparable estimates. The expected detection of increasing depolarization with depth into the canopies due to volume multiple-scattering was not observed, possibly due to the small laser footprint and the small detector field of view used in the SIMPL instrument. The results of this work provide pathfinder information for NASA's ICESat-2 mission that will employ a 532 nm, micropulse, photon counting laser altimeter.

Harding, David J.↗

An investigation of using an RQP based method to calculate parameter sensitivity derivatives

Estimation of the sensitivity of problem functions with respect to problem variables forms the basis for many of our modern day algorithms for engineering optimization. The most common application of problem sensitivities has been in the calculation of objective function and constraint partial derivatives for determining search directions and optimality conditions. A second form of sensitivity analysis, parameter sensitivity, has also become an important topic in recent years. By parameter sensitivity, researchers refer to the estimation of changes in the modeling functions and current design point due to small changes in the fixed parameters of the formulation. Methods for calculating these derivatives have been proposed by several authors (Armacost and Fiacco 1974, Sobieski et al 1981, Schmit and Chang 1984, and Vanderplaats and Yoshida 1985). Two drawbacks to estimating parameter sensitivities by current methods have been: (1) the need for second order information about the Lagrangian at the current point, and (2) the estimates assume no change in the active set of constraints. The first of these two problems is addressed here and a new algorithm is proposed that does not require explicit calculation of second order information.

Beltracchi, Todd J.↗

Noncoherent DTTLs for Symbol Synchronization

Noncoherent data-transition tracking loops (DTTLs) have been proposed for use as symbol synchronizers in digital communication receivers. [Communication- receiver subsystems that can perform their assigned functions in the absence of synchronization with the phases of their carrier signals ( carrier synchronization ) are denoted by the term noncoherent, while receiver subsystems that cannot function without carrier synchronization are said to be coherent. ] The proposal applies, more specifically, to receivers of binary phase-shift-keying (BPSK) signals generated by directly phase-modulating binary non-return-to-zero (NRZ) data streams onto carrier signals having known frequencies but unknown phases. The proposed noncoherent DTTLs would be modified versions of traditional DTTLs, which are coherent. The symbol-synchronization problem is essentially the problem of recovering symbol timing from a received signal. In the traditional, coherent approach to symbol synchronization, it is necessary to establish carrier synchronization in order to recover symbol timing. A traditional DTTL effects an iterative process in which it first generates an estimate of the carrier phase in the absence of symbol-synchronization information, then uses the carrier-phase estimate to obtain an estimate of the symbol-synchronization information, then feeds the symbol-synchronization estimate back to the carrier-phase-estimation subprocess. In a noncoherent symbol-synchronization process, there is no need for carrier synchronization and, hence, no need for iteration between carrier-synchronization and symbol- synchronization subprocesses. The proposed noncoherent symbolsynchronization process is justified theoretically by a mathematical derivation that starts from a maximum a posteriori (MAP) method of estimation of symbol timing utilized in traditional, coherent DTTLs. In that MAP method, one chooses the value of a variable of interest (in this case, the offset in the estimated symbol timing) that causes a likelihood function of symbol estimates over some number of symbol periods to assume a maximum value. In terms that are necessarily oversimplified to fit within the space available for this article, it can be said that the mathematical derivation involves a modified interpretation of the likelihood function that lends itself to noncoherent DTTLs. The proposal encompasses both linear and nonlinear noncoherent DTTLs. The performances of both have been computationally simulated; for comparison, the performances of linear and nonlinear coherent DTTLs have also been computationally simulated. The results of these simulations show that, among other things, the expected mean-square timing errors of coherent and noncoherent DTTLs are relatively insensitive to window width. The results also show that at high signal-to-noise ratios (SNRs), the performances of the noncoherent DTTLs approach those of their coherent counterparts at, while at low SNRs, the noncoherent DTTLs incur penalties of the order of 1.5 to 2 dB.

Simon, Marvin↗

Mode Shape Estimation Algorithms Under Ambient Conditions: A Comparative Review

Abstract—This paper provides a comparative review of five existing ambient electromechanical mode shape estimation algorithms, i.e., the Transfer Function (TF), Spectral, Frequency Domain Decomposition (FDD), Channel Matching, and Subspace Methods. It is also shown that the TF Method is a general approach to estimating mode shape and that the Spectral, FDD, and Channel Matching Methods are actually special cases of it. Additionally, some of the variations of the Subspace Method are reviewed and the Numerical algorithm for Subspace State Space System IDentification (N4SID) is implemented. The five algorithms are then compared using data simulated from a 17-machine model of the Western Electricity Coordinating Council (WECC) under ambient conditions with both low and high damping, as well as during the case where ambient data is disrupted by an oscillatory ringdown. The performance of the algorithms is compared using the statistics from Monte Carlo Simulations and results from measured WECC data, and a discussion of the practical issues surrounding their implementation, including cases where power system probing is an option, is provided. The paper concludes with some recommendations as to the appropriate use of the various techniques. Index Terms—Electromechanical mode shape, small-signal stability, phasor measurement units (PMU), system identification, N4SID, subspace.

Dosiek, Luke↗

Inverse regression for ridge recovery: a data-driven approach for parameter reduction in computer experiments

Parameter reduction can enable otherwise infeasible design and uncertainty studies with modern computational science models that contain several input parameters. In statistical regression, techniques for sufficient dimension reduction (SDR) use data to reduce the predictor dimension of a regression problem. A computational scientist hoping to use SDR for parameter reduction encounters a problem: a computer prediction is best represented by a deterministic function of the inputs, so data comprised of computer simulation queries fail to satisfy the SDR assumptions. To address this problem, here we interpret SDR methods sliced inverse regression (SIR) and sliced average variance estimation (SAVE) as estimating the directions of a ridge function, which is a composition of a low-dimensional linear transformation with a nonlinear function. Within this interpretation, SIR and SAVE estimate matrices of integrals whose column spaces are contained in the ridge directions’ span; we analyze and numerically verify convergence of these column spaces as the number of computer model queries increases. Moreover, we show example functions that are not ridge functions but whose inverse conditional moment matrices are low-rank. Consequently, the computational scientist should beware when using SIR and SAVE for parameter reduction, since SIR and SAVE may mistakenly suggest that truly important directions are unimportant.

97 MATHEMATICS AND COMPUTING↗

Maximum Likelihood Estimation: Some Basics

The maximum likelihood estimation is a general estimation procedure. It is often compared to estimation procedures like the ordinary least squares regression or generalized method of moments, to name a few. We discuss some basics about the maximum likelihood estimation, its advantages and disadvantages, and provide an example application to a gamma distribution function.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Maximum Likelihood Estimation: Some Basics

The maximum likelihood estimation is a general estimation procedure. It is often compared to estimation procedures like the ordinary least squares regression or generalized method of moments, to name a few. We discuss some basics about the maximum likelihood estimation, its advantages and disadvantages, and provide an example application to a gamma distribution function.

97 MATHEMATICS AND COMPUTING↗

Evaluation of high temperature superconductive thermal bridges for space borne cryogenic detectors

Infrared sensor satellites are used to monitor the conditions in the earth's upper atmosphere. In these systems, the electronic links connecting the cryogenically cooled infrared detectors to the significantly warmer amplification electronics act as thermal bridges and, consequently, the mission lifetimes of the satellites are limited due to cryogenic evaporation. High-temperature superconductor (HTS) materials have been proposed by researchers at the National Aeronautics and Space Administration Langley's Research Center (NASA-LaRC) as an alternative to the currently used manganin wires for electrical connection. The potential for using HTS films as thermal bridges has provided the motivation for the design and the analysis of a spaceflight experiment to evaluate the performance of this superconductive technology in the space environment. The initial efforts were focused on the preliminary design of the experimental system which allows for the quantitative comparison of superconductive leads with manganin leads, and on the thermal conduction modeling of the proposed system. Most of the HTS materials were indicated to be potential replacements for the manganin wires. In the continuation of this multi-year research, the objectives of this study were to evaluate the sources of heat transfer on the thermal bridges that have been neglected in the preliminary conductive model and then to develop a methodology for the estimation of the thermal conductivities of the HTS thermal bridges in space. The Joule heating created by the electrical current through the manganin wires was incorporated as a volumetric heat source into the manganin conductive model. The radiative heat source on the HTS thermal bridges was determined by performing a separate radiant interchange analysis within a high-T(sub c) superconductor housing area. Both heat sources indicated no significant contribution on the cryogenic heat load, which validates the results obtained in the preliminary conduction model. A methodology was presented for the estimation of the thermal conductivities of the individual HTS thermal bridge materials and the effective thermal conductivities of the composite HTS thermal bridges as functions of temperature. This methodology included a sensitivity analysis and the demonstration of the estimation procedure using simulated data with added random errors. The thermal conductivities could not be estimated as functions of temperature; thus the effective thermal conductivities of the HTS thermal bridges were analyzed as constants.

Scott, Elaine P.↗

The most luminous z ∼ 9-10 galaxy candidates yet found: The luminosity function, cosmic star-formation rate, and the first mass density estimate at 500 Myr

We present the discovery of four surprisingly bright (H {sub 160} ∼ 26-27 mag AB) galaxy candidates at z ∼ 9-10 in the complete HST CANDELS WFC3/IR GOODS-N imaging data, doubling the number of z ∼ 10 galaxy candidates that are known, just ∼500 Myr after the big bang. Two similarly bright sources are also detected in a reanalysis of the GOODS-S data set. Three of the four galaxies in GOODS-N are significantly detected at 4.5σ-6.2σ in the very deep Spitzer/IRAC 4.5 μm data, as is one of the GOODS-S candidates. Furthermore, the brightest of our candidates (at z = 10.2 ± 0.4) is robustly detected also at 3.6 μm (6.9σ), revealing a flat UV spectral energy distribution with a slope β = –2.0 ± 0.2, consistent with demonstrated trends with luminosity at high redshift. Thorough testing and use of grism data excludes known low-redshift contamination at high significance, including single emission-line sources, but as-yet unknown low redshift sources could provide an alternative solution given the surprising luminosity of these candidates. Finding such bright galaxies at z ∼ 9-10 suggests that the luminosity function for luminous galaxies might evolve in a complex way at z > 8. The cosmic star formation rate density still shows, however, an order-of-magnitude increase from z ∼ 10 to z ∼ 8 since the dominant contribution comes from low-luminosity sources. Based on the IRAC detections, we derive galaxy stellar masses at z ∼ 10, finding that these luminous objects are typically 10{sup 9} M {sub ☉}. This allows for a first estimate of the cosmic stellar mass density at z ∼ 10 resulting in log{sub 10} ρ{sub ∗}=4.7{sub −0.8}{sup +0.5} M {sub ☉} Mpc{sup –3} for galaxies brighter than M {sub UV} ∼ –18. The remarkable brightness, and hence luminosity, of these z ∼ 9-10 candidates will enable deep spectroscopy to determine their redshift and nature, and highlights the opportunity for the James Webb Space Telescope to map the buildup of galaxies at redshifts much earlier than z ∼ 10.

79 ASTRONOMY AND ASTROPHYSICS↗

Kidney function and blood pressure in preschool-aged children exposed to cadmium and arsenic - potential alleviation by selenium

Background: Early-life exposure to toxic compounds may cause long-lasting health effects, but few studies have investigated effects of childhood exposure to nephrotoxic metals on kidney and cardiovascular function. Objectives: To assess effects of exposure to arsenic and cadmium on kidney function and blood pressure in pre-school-aged children, and potential protection by selenium. Methods: This cross-sectional study was part of the 4.5 years of age (range: 4.4–5.4 years) follow-up of the children from a supplementation trial in pregnancy (MINIMat) in rural Bangladesh, and nested studies on early-life metal exposures. Exposure to arsenic, cadmium and selenium from food and drinking water was assessed by concentrations in children's urine, measured by ICP-MS. Kidney function was assessed by the estimated glomerular filtration rate (eGFR, n=1106), calculated from serum cystatin C, and by kidney volume, measured by ultrasound (n=375). Systolic and diastolic blood pressure was measured (n=1356) after five minutes rest. Results: Multivariable-adjusted regression analyzes showed that exposure to cadmium, but not arsenic, was inversely associated with eGFR, particularly in girls. A 0.5 µg/L increase in urinary cadmium among the girls (above spline knot at 0.12) was associated with a decrease in eGFR of 2.6 ml/min/1.73 m{sup 2}, corresponding to 0.2SD (p=0.022). A slightly weaker inverse association with cadmium was also indicated for kidney volume, but no significant associations were found with blood pressure. Stratifying on children's urinary selenium (below or above median of 12.6 µg/L) showed a three times stronger inverse association of U-Cd with eGFR (all children) in the lower selenium stratum (B=−2.8; 95% CI: −5.5, −0.20; p=0.035), compared to those with higher selenium (B=−0.79; 95% CI: −3.0, 1.4; p=0.49). Conclusions: Childhood cadmium exposure seems to adversely affect kidney function, but not blood pressure, in this population of young children in rural Bangladesh. Better selenium status appears to be protective. However, it is important to follow up these children to assess potential long-term consequences of these findings. - Highlights: • Bangladesh has high levels of arsenic and cadmium in drinking water and food • We assessed toxic exposure, kidney and cardiovascular function in Bangladeshi children • Cadmium appeared to decrease estimated glomerular filtration rate in girls • These effects were somewhat alleviated by selenium.

54 ENVIRONMENTAL SCIENCES↗

Determination of lift and drag characteristics of Space Shuttle Orbiter using maximum likelihood estimation technique

This paper presents the technique and results of maximum likelihood estimation used to determine lift and drag characteristics of the Space Shuttle Orbiter. Maximum likelihood estimation uses measurable parameters to estimate nonmeasurable parameters. The nonmeasurable parameters for this case are elements of a nonlinear, dynamic model of the orbiter. The estimated parameters are used to evaluate a cost function that computes the differences between the measured and estimated longitudinal parameters. The case presented is a dynamic analysis. This places less restriction on pitching motion and can provide additional information about the orbiter such as lift and drag characteristics at conditions other than trim, instrument biases, and pitching moment characteristics. In addition, an output of the analysis is an estimate of the values for the individual components of lift and drag that contribute to the total lift and drag. The results show that maximum likelihood estimation is a useful tool for analysis of Space Shuttle Orbiter performance and is also applicable to parameter analysis of other types of aircraft.

Trujillo, B. M.↗

Transfer function analysis of dynamic cerebral autoregulation in humans

To test the hypothesis that spontaneous changes in cerebral blood flow are primarily induced by changes in arterial pressure and that cerebral autoregulation is a frequency-dependent phenomenon, we measured mean arterial pressure in the finger and mean blood flow velocity in the middle cerebral artery (VMCA) during supine rest and acute hypotension induced by thigh cuff deflation in 10 healthy subjects. Transfer function gain, phase, and coherence function between changes in arterial pressure and VMCA were estimated using the Welch method. The impulse response function, calculated as the inverse Fourier transform of this transfer function, enabled the calculation of transient changes in VMCA during acute hypotension, which was compared with the directly measured change in VMCA during thigh cuff deflation. Beat-to-beat changes in VMCA occurred simultaneously with changes in arterial pressure, and the autospectrum of VMCA showed characteristics similar to arterial pressure. Transfer gain increased substantially with increasing frequency from 0.07 to 0.20 Hz in association with a gradual decrease in phase. The coherence function was > 0.5 in the frequency range of 0.07-0.30 Hz and < 0.5 at < 0.07 Hz. Furthermore, the predicted change in VMCA was similar to the measured VMCA during thigh cuff deflation. These data suggest that spontaneous changes in VMCA that occur at the frequency range of 0.07-0.30 Hz are related strongly to changes in arterial pressure and, furthermore, that short-term regulation of cerebral blood flow in response to changes in arterial pressure can be modeled by a transfer function with the quality of a high-pass filter in the frequency range of 0.07-0.30 Hz.

NASA Discipline Cardiopulmonary↗

Development of a Portfolio Management Approach with Case Study of the NASA Airspace Systems Program

A portfolio management approach was developed for the National Aeronautics and Space Administration s (NASA s) Airspace Systems Program (ASP). The purpose was to help inform ASP leadership regarding future investment decisions related to its existing portfolio of advanced technology concepts and capabilities (C/Cs) currently under development and to potentially identify new opportunities. The portfolio management approach is general in form and is extensible to other advanced technology development programs. It focuses on individual C/Cs and consists of three parts: 1) concept of operations (con-ops) development, 2) safety impact assessment, and 3) benefit-cost-risk (B-C-R) assessment. The first two parts are recommendations to ASP leaders and will be discussed only briefly, while the B-C-R part relates to the development of an assessment capability and will be discussed in greater detail. The B-C-R assessment capability enables estimation of the relative value of each C/C as compared with all other C/Cs in the ASP portfolio. Value is expressed in terms of a composite weighted utility function (WUF) rating, based on estimated benefits, costs, and risks. Benefit utility is estimated relative to achieving key NAS performance objectives, which are outlined in the ASP Strategic Plan.1 Risk utility focuses on C/C development and implementation risk, while cost utility focuses on the development and implementation portions of overall C/C life-cycle costs. Initial composite ratings of the ASP C/Cs were successfully generated; however, the limited availability of B-C-R information, which is used as inputs to the WUF model, reduced the meaningfulness of these initial investment ratings. Development of this approach, however, defined specific information-generation requirements for ASP C/C developers that will increase the meaningfulness of future B-C-R ratings.

Neitzke, Kurt W.↗

Lidar-Based Estimates of Above-Ground Biomass in the Continental US and Mexico Using Ground, Airborne, and Satellite Observations

Existing national forest inventory plots, an airborne lidar scanning (ALS) system, and a space profiling lidar system (ICESat-GLAS) are used to generate circa 2005 estimates of total aboveground dry biomass (AGB) in forest strata, by state, in the continental United States (CONUS) and Mexico. The airborne lidar is used to link ground observations of AGB to space lidar measurements. Two sets of models are generated, the first relating ground estimates of AGB to airborne laser scanning (ALS) measurements and the second set relating ALS estimates of AGB (generated using the first model set) to GLAS measurements. GLAS then, is used as a sampling tool within a hybrid estimation framework to generate stratum-, state-, and national-level AGB estimates. A two-phase variance estimator is employed to quantify GLAS sampling variability and, additively, ALS-GLAS model variability in this current, three-phase (ground-ALS-space lidar) study. The model variance component characterizes the variability of the regression coefficients used to predict ALS-based estimates of biomass as a function of GLAS measurements. Three different types of predictive models are considered in CONUS to determine which produced biomass totals closest to ground-based national forest inventory estimates - (1) linear (LIN), (2) linear-no-intercept (LNI), and (3) log-linear. For CONUS at the national level, the GLAS LNI model estimate (23.95 +/- 0.45 Gt AGB), agreed most closely with the US national forest inventory ground estimate, 24.17 +/- 0.06 Gt, i.e., within 1%. The national biomass total based on linear ground-ALS and ALS-GLAS models (25.87 +/- 0.49 Gt) overestimated the national ground-based estimate by 7.5%. The comparable log-linear model result (63.29 +/-1.36 Gt) overestimated ground results by 261%. All three national biomass GLAS estimates, LIN, LNI, and log-linear, are based on 241,718 pulses collected on 230 orbits. The US national forest inventory (ground) estimates are based on 119,414 ground plots. At the US state level, the average absolute value of the deviation of LNI GLAS estimates from the comparable ground estimate of total biomass was 18.8% (range: Oregon,−40.8% to North Dakota, 128.6%). Log-linear models produced gross overestimates in the continental US, i.e., N2.6x, and the use of this model to predict regional biomass using GLAS data in temperate, western hemisphere forests is not appropriate. The best model form, LNI, is used to produce biomass estimates in Mexico. The average biomass density in Mexican forests is 53.10 +/- 0.88 t/ha, and the total biomass for the country, given a total forest area of 688,096 sq km, is 3.65 +/- 0.06 Gt. In Mexico, our GLAS biomass total underestimated a 2005 FAO estimate (4.152 Gt) by 12% and overestimated a 2007/8 radar study's figure (3.06 Gt) by 19%.

Hybrid 3-Phase sampling↗

Methods of Constructing a Blended Performance Function Suitable for Formation Flight

Two methods for constructing performance functions for formation fight-for-drag-reduction suitable for use with an extreme-seeking control system are presented. The first method approximates an a prior measured or estimated drag-reduction performance function by combining real-time measurements of readily available parameters. The parameters are combined with weightings determined from a minimum squares optimization to form a blended performance function.

control systems↗

Isolating Added Mass Load Components of CPAS Main Clusters

The current simulation for the Capsule Parachute Assembly System (CPAS) lacks fidelity in representing added mass for the 116 ft Do ringsail Main parachute. The availability of 3-D models of inflating Main canopies allowed for better estimation the enclosed air volume as a function of time. This was combined with trajectory state information to estimate the components making up measured axial loads. A proof-of-concept for an alternate simulation algorithm was developed based on enclosed volume as the primary independent variable rather than drag area growth. Databases of volume growth and parachute drag area vs. volume were developed for several flight tests. Other state information was read directly from test data, rather than numerically propagated. The resulting simulated peak loads were close in timing and magnitude to the measured loads data. However, results are very sensitive to data curve fitting and may not be suitable for Monte Carlo simulations. It was assumed that apparent mass was either negligible or a small fraction of enclosed mass, with little difference in results.

Ray, Eric S.↗