Search NASA⌕ Search

SEARCH · Search NASA

Results for “linear 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,225 records · Page 68

Validation of the Pulmonary Function System for Use on the International Space Station

Aerobic deconditioning occurs during long duration space flight despite the use of exercise countermeasures (Convertino, 1996). As a part of International Space Station (ISS) medical operations, periodic tests designed to estimate aerobic capacity are performed to track changes in aerobic fitness and to determine the effectiveness of exercise countermeasures. These tests are performed prior to, during, and after missions of greater than 30 days in duration. Crewmembers selected for missions aboard the ISS perform a graded exercise test on a cycle ergometer approximately 270 days prior to their scheduled launch date in order to measure peak oxygen consumption (VO2PK) and peak heart rate (HRpk). Approximately 30 to 45 days prior to launch, crewmembers perform a submaximal cycle ergometer test at work rates set to elicit 25, 50 and 75% of their pre-flight VO2PK. This test, known as the Periodic Fitness Evaluation (PFE), serves as a baseline measure to which subsequent in-and post-flight exercise tests are compared. While onboard the ISS, crewmembers are normally scheduled to perform the PFE beginning with flight day (FD) 14 and every 30 days thereafter. The PFE is also conducted 5 and 30 days following flight. Using PFE data, aerobic fitness is estimated by quantifying the VO2 vs. HR relationship using linear regression and calculating the VO2 that would occur at the crewmember s previously measured HRpk. Currently, for data collected during flight, this technique assumes that the pre- vs. in-flight oxygen consumption per given cycle workload is similar. However, the validity of this assumption is based upon a sparse amount of data collected during the Skylab era (Michel, et al. 1977). The method of using heart rate and cycle ergometer work rates has been used to estimate aerobic fitness in normal gravity (Astrand and Ryhming, 1954; Lee, 1993). Due to spaceflight induced physiological alterations, such as shifts in extracellular fluid (e.g. plasma) volume, this method may not be valid during space flight. In addition, the ergometer onboard ISS is vibration-isolated and moves with the astronaut s application of force into the pedals. The effect of this movement on the VO2 of cycle exercise on ISS has not been quantified.

McCleary, Frank A.↗

NaCl nucleation from brine in seeded simulations: Sources of uncertainty in rate estimates

This work reexamines seeded simulation results for NaCl nucleation from a supersaturated aqueous solution at 298.15 K and 1 bar pressure. We present a linear regression approach for analyzing seeded simulation data that provides both nucleation rates and uncertainty estimates. Our results show that rates obtained from seeded simulations rely critically on a precise driving force for the model system. The driving force vs. solute concentration curve need not exactly reproduce that of the real system, but it should accurately describe the thermodynamic properties of the model system. We also show that rate estimates depend strongly on the nucleus size metric. We show that the rate estimates systematically increase as more stringent local order parameters are used to count members of a cluster and provide tentative suggestions for appropriate clustering criteria.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Estimating uncertainty in ambient and saturation nutrient uptake metrics from nutrient pulse releases in stream ecosystems

Nutrient spiraling is an important ecosystem process characterizing nutrient transport and uptake in streams. Various nutrient addition methods are used to estimate uptake metrics; however, uncertainty in the metrics is not often evaluated. A method was developed to quantify uncertainty in ambient and saturation nutrient uptake metrics estimated from saturating pulse nutrient additions (Tracer Additions for Spiraling Curve Characterization; TASCC). Using a Monte Carlo (MC) approach, the 95% confidence interval (CI) was estimated for ambient uptake lengths (S w-amb ) and maximum areal uptake rates (U max ) based on 100,000 datasets generated from each of four nitrogen and five phosphorous TASCC experiments conducted seasonally in a forest stream in eastern Tennessee, U.S.A. Uncertainty estimates from the MC approach were compared to the CIs estimated from ordinary least squares (OLS) and non-linear least squares (NLS) models used to calculate S w-amb and U max , respectively, from the TASCC method. The CIs for Sw-amb and Umax were large, but were not consistently larger using the MC method. Despite the large CIs, significant differences (based on nonoverlapping CIs) in nutrient metrics among seasons were found with more significant differences using the OLS/NLS vs. the MC method. Lastly, we suggest that the MC approach is a robust way to estimate uncertainty, as the calculation of S w-amb and U max violates assumptions of OLS/NLS while the MC approach is free of these assumptions. The MC approach can be applied to other ecosystem metrics that are calculated from multiple parameters, providing a more robust estimate of these metrics and their associated uncertainties.

54 ENVIRONMENTAL SCIENCES↗

COPSS II: the molecular gas content of ten million cubic megaparsecs at redshift z ∼ 3

We present a measurement of the abundance of carbon monoxide in the early universe, utilizing the final results from the CO Power Spectrum Survey (COPSS). Between 2013 and 2015, we performed observations with the Sunyaev–Zel’dovich Array to measure aggregate CO emission from z∼3 galaxies with the intensity mapping technique. Data were collected on 19 fields, covering an area of 0.7 square degrees, over the frequency range 27--35 GHz. With these data, along with data analyzed in COPSS I, we are able to observe the CO(1–0) transition within the redshift range z=2.3--3.3 for spatial frequencies between k=0.5--10 h Mpc{sup −1}, spanning a comoving volume of 4.9×10{sup 6} h{sup −3} Mpc{sup 3}. We present estimates of contributions from continuum sources and ground illumination within our measurement. We constrain the amplitude of the CO power spectrum to P{sub CO}=3.0{sub −1.3}{sup +1.3}×10{sup 3} μK{sup 2}(h{sup −1}Mpc){sup 3}, or Δ{sub CO}{sup 2}(k=1 h Mpc{sup −1})=1.5{sub −0.7}{sup +0.7}×10{sup 3} μK{sup 2}, at 68% confidence, and P{sub CO}>0 at 98.9% confidence. These results are a factor of 10 improvement in sensitivity compared to those of COPSS I. With this measurement, we constrain on the CO(1–0) galaxy luminosity function at z∼3. Assuming that CO emission is proportional to halo mass and using theoretical estimates of the scatter in this relationship, we constrain the ratio of CO(1--0) luminosity to halo mass to A{sub CO}=6.3{sub −2.1}{sup +1.4}×10{sup −7} L{sub ⊙} M{sub ⊙}{sup −1}. Assuming a Milky Way-like linear relationship between CO luminosity and molecular gas mass, we estimate a mass fraction of molecular gas of f{sub H{sub 2}}=5.5{sub −2.2}{sup +3.4}×10{sup −2} for halos with masses of ∼10{sup 12}M{sub ⊙}. Using theoretical estimates for the scaling of molecular gas mass fraction and halo mass, we estimate the cosmic molecular gas density to be ρ{sub z∼3}(H{sub 2})=1.1{sub −0.4}{sup +0.7}×10{sup 8} M{sub ⊙} Mpc{sup −3}.

79 ASTRONOMY AND ASTROPHYSICS↗

Computational Intelligence Based Data Fusion Algorithm for Dynamic sEMG and Skeletal Muscle Force Modelling

In this work, an array of three surface Electrography (sEMG) sensors are used to acquired muscle extension and contraction signals for 18 healthy test subjects. The skeletal muscle force is estimated using the acquired sEMG signals and a Non-linear Wiener Hammerstein model, relating the two signals in a dynamic fashion. The model is obtained from using System Identification (SI) algorithm. The obtained force models for each sensor are fused using a proposed fuzzy logic concept with the intent to improve the force estimation accuracy and resilience to sensor failure or misalignment. For the fuzzy logic inference system, the sEMG entropy, the relative error, and the correlation of the force signals are considered for defining the membership functions. The proposed fusion algorithm yields an average of 92.49% correlation between the actual force and the overall estimated force output. In addition, the proposed fusionbased approach is implemented on a test platform. Experiments indicate an improvement in finger/hand force estimation.

42 ENGINEERING↗

Fault Detection for Dynamical Systems using Differential Geometric and Concurrent Learning Approach

This paper presents a fault identification and estimation approach that brings together the differential geometric concept of observability codistribution with data driven concurrent learning. In order to identify faults in presence of unknown disturbances, we utilize the differential geometric approach to design a coordinate transformation, to find a subspace in which the effect of disturbances and system faults can be segregated. We then use concurrent learning to estimate magnitude of the constant fault. We illustrate the approach to fault isolation for a linear helicopter dynamics and a spherical pendulum dynamics. We use Lyapunov stability analysis to show that the fault estimate by concurrent learning converges to the actual fault value, and then illustrate the design of a recovery controller.

Fault detection and diagnosis↗

Anisotropic parameter estimation using velocity variation with offset analysis

Seismic anisotropy is defined as velocity dependent upon angle or offset. Knowledge about anisotropy effect on seismic data is important in amplitude analysis, stacking process and time to depth conversion. Due to this anisotropic effect, reflector can not be flattened using single velocity based on hyperbolic moveout equation. Therefore, after normal moveout correction, there will still be residual moveout that relates to velocity information. This research aims to obtain anisotropic parameters, ε and δ, using two proposed methods. The first method is called velocity variation with offset (VVO) which is based on simplification of weak anisotropy equation. In VVO method, velocity at each offset is calculated and plotted to obtain vertical velocity and parameter δ. The second method is inversion method using linear approach where vertical velocity, δ, and ε is estimated simultaneously. Both methods are tested on synthetic models using ray-tracing forward modelling. Results show that δ value can be estimated appropriately using both methods. Meanwhile, inversion based method give better estimation for obtaining ε value. This study shows that estimation on anisotropic parameters rely on the accuracy of normal moveout velocity, residual moveout and offset to angle transformation.

58 GEOSCIENCES↗

The Conductor-Dielectric Junctions in a Low Density Plasma

A conductor-dielectric junction exposed to the space environment is a frequent spacecraft design feature. Due to spacecraft charging and/or solar array operation, the conductor can acquire a high potential with respect to the surrounding plasma. If this potential is positive the insulators adjacent to exposed conductors can collect current as if they were conductors themselves. This phenomenon, called snapover, results in a substantial increase in current collection, and may even result in a glow discharge if the potential is high enough. If a conductor has a negative potential, arcing can occur at the site of a junction. Both of these phenomena negatively affect spacecraft operation. To prevent negative consequences, the physical mechanisms of snapover and arc inception require investigation. In this paper, results are presented of an experimental and theoretical study of snapover, glow discharge, and arc phenomena for different materials immersed in argon or xenon plasmas. The effect of snapover is investigated for several metal-dielectric junctions: copper-teflon, copper-Kapton, copper-glass, aluminum-teflon, aluminum-Kapton, steel-teflon, anodized aluminum with pinholes, and copper-ceramics. I-V curves are measured and snapover inception voltages, essential parameters (increase in current and collection area due to secondary electrons), and glow discharge inception thresholds are determined. Optical spectra are obtained for glow discharges in both argon and xenon plasmas. These spectra provide information regarding atomic species entrapped in the glow region. Some spectral lines can be used to estimate plasma parameters in the discharge area. A video-camera and linear array were used to confirm that snapover inception is accompanied by very low intensity visible light emission. This result seems to be important for the estimate of the light pollution around spacecraft. Optical spectra (wavelengths 380-650 nm) of arcs are also obtained on a negatively biased chromic acid anodized aluminum plate immersed in low density argon and xenon plasmas. Analysis of these spectra confirms our earlier findings that aluminum atoms are ejected from the arc site. Moreover, it is found that chromium atoms are also quite abundant in the arc plasma. It is believed that the latter results contribute considerably to the understanding of processes of plasma contamination caused by arcing.

Vayner, Boris↗

Reactive scattering calculations for {sup 87}Rb+{sup 87}RbHe→Rb{sub 2}({sup 3}Σ{sub u}{sup +},v)+He from ultralow to intermediate energies

We investigate atom-diatom reactive collisions, as a preliminary step, in order to assess the possibility of forming Rb{sub 2} molecules in their lowest triplet electronic state by cold collisions of rubidium atoms on the surface of helium nanodroplets. A simple model related to the well-known Rosen treatment of linear triatomic molecules [N. Rosen, J. Chem. Phys. 1, 319 (1933)] in relative coordinates is used, allowing to estimate reactive probabilities for different values of the total angular momentum. The best available full dimensional potential energy surface [Guillon et al., J. Chem. Phys. 136, 174307 (2012)] is employed through the calculations. Noticeable values of the probabilities in the ultracold regime, which numerically fulfill the Wigner threshold law, support the feasibility of the process. The rubidium dimer is mainly produced at high vibrational states, and the reactivity is more efficient for a bosonic helium partner than when the fermion species is considered.

74 ATOMIC AND MOLECULAR PHYSICS↗

Linear Correlations of Gibbs Free Energy of REE Phosphates (Monazite, Xenotime, and Rhabdophane) and Internally Consistent Binary Mixing Properties

Rare Earth Elements (REE) phosphates (monazite, xenotime, and rhabdophane) are critical REE-bearing minerals typically formed in hydrothermal and magmatic ore deposits. The thermodynamic properties of those REE minerals are crucial to understanding the solubility, speciation, and transport of REE complexes. However, the standard-state Gibbs free energy of formation (∆G°f) values reported for these minerals in the literature vary by up to 25 kJ mol−1. Here, we present linear free energy relationships that allow the evaluation and estimation of the ∆G°f values at 25 °C and 1 bar for the three minerals from the ionic radius (rREE3+) and the non-solvation Gibbs free energy contribution to the REE3+ aqua ion (∆G°n, REE3+): ∆G°f,monazite − 399.71 rREE3+ = 1.0059 ∆G°n,REE3+ − 2522.51; ∆G°f,xenotime − 344.08 rREE3+ = 0.9909 ∆G°n,REE3+ − 2451.53; and ∆G°f,rhabdophane − 416.17 rREE3+ = 1.0067 ∆G°n, REE3+ − 2688.86. Moreover, based on the new dataset derived for REE end-members, we re-fitted the binary Margules parameter (W) from previous theoretical calculations into linear correlations: W + 0.00204 ∆G°’n,monazite = 39.3549 ∆V + 0.0641; W + 0.00255 ∆G°’n,xenotime = 25.4885 ∆V − 0.0062. The internally consistent thermodynamic properties of these REE phosphates are incorporated into the computer program Supcrtbl, which is available online at Zhu’s research website.

Pan, Ruiguang (ORCID:0000000324896545)↗

Computation of output feedback gains for linear stochastic systems using the Zangnill-Powell Method

Because conventional optimal linear regulator theory results in a controller which 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. To this effect a stochastic linear model has been developed that accounts for process parameter and initial uncertainty, measurement noise, and a restricted number of measurable outputs. Optimization with respect to the corresponding output feedback gains was then 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.↗

Computation of output feedback gains for linear stochastic systems using the Zangwill-Powell method

Because conventional optimal linear regulator theory results in a controller which 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. To this effect a stochastic linear model has been developed that accounts for process parameter and initial uncertainty, measurement noise, and a restricted number of measurable outputs. Optimization with respect to the corresponding output feedback gains was then performed for both finite and infinite time performance indices without gradient computation by using Zangwill's modification of a procedure originally proposed by Powell.

Kaufman, H.↗

Optical properties of intact leaves for estimating chlorophyll concentration

Changes in leaf chlorophyll content can serve as relative indicators of plant vigor and environmental quality. This study identified reflectance, transmittance, and absorptance wavebands and band ratios within the 400- to 850-nm range for intact leaves that could be used to estimate extracted leaf chlorophyll per unit leaf area (areal concentration) with minimal error. Leaf optical properties along with chlorophyll a, b, and a + b concentrations were measured for the planar-leaved sweetgum (Liquidambar styraciflua L.), red maple (Acer rubrum L.), wild grape (Vitis rotundifolia Michx.), and switchcane [Arundinaria gigantea (Walter) Muhl.], and for needles of longleaf pine (Pinus palustris Miller). Generally, reflectance, transmittance, and absorptance corresponded most precisely with chlorophyll concentrations at wavelengths near 700 nm, although regressions were also strong in the 550- to 625-nm range. A power function was superior to a simple linear function in yielding low standard deviations of the estimate (s). When data were combined among the planar-leaved species, s values were low at approximately 50 mumol/m2 out of a 940 mumol/m2 range in chlorophyll a + b at best-fit wavelengths of 707 to 709 nm. Minimal s values for chlorophyll a + b ranged from 32 to 62 mumol/m2 across species when band ratios having numerator wavelengths of 693 to 720 nm were used with the application of a power function. Optimal denominator wavelengths for the band ratios were 850 nm for reflectance and transmittance and 400 nm for absorptance. This information can be applied in designing field portable chlorophyll meters and in the landscape-scale remote sensing of plant responses to the environment.

Chlorophyll/analysis↗

Error Estimates for the Optimal Control of a Parabolic Fractional PDE

In this work, we consider the integral definition of the fractional Laplacian and analyze a linear-quadratic optimal control problem for the so-called fractional heat equation; control constraints are also considered. We derive existence and uniqueness results, first order optimality conditions, and regularity estimates for the optimal variables. To discretize the state equation we propose a fully discrete scheme that relies on an implicit finite difference discretization in time combined with a piecewise linear finite element discretization in space. We derive stability results and a novel $L^2(0,T;L^2(\Omega))$ a priori error estimate. On the basis of the aforementioned solution technique, we propose a fully discrete scheme for our optimal control problem that discretizes the control variable with piecewise constant functions, and we derive a priori error estimates for it. We illustrate the theory with one- and two-dimensional numerical experiments.

97 MATHEMATICS AND COMPUTING↗

Microstructures in type III events in the solar wind

A model is presented for the generation and evolution of bump-in-tail driven Langmuir waves in the solar wind during type III emission which removes a number of apparent inconsistencies between theory and observations. Growth rates and energy densities of Langmuir waves are considerably enhanced, permitting growth to overcome linear scattering losses, and also allowing nonlinear decay into ion-acoustic waves, in line with observations. Estimates are made of the probability distribution p(E), of wave field strengths E, based on linear and nonlinear wave-packet evolution, yielding p(E) approximately equal to E exp -alpha, alpha greater than or equal to 3. This helps explain why very high values of E are rarely found in the measured spiky wave turbulence.

Melrose, D. B.↗

Comparing hard and soft prior bounds in geophysical inverse problems

In linear inversion of a finite-dimensional data vector y to estimate a finite-dimensional prediction vector z, prior information about X sub E is essential if y is to supply useful limits for z. The one exception occurs when all the prediction functionals are linear combinations of the data functionals. Two forms of prior information are compared: a soft bound on X sub E is a probability distribution p sub x on X which describeds the observer's opinion about where X sub E is likely to be in X; a hard bound on X sub E is an inequality Q sub x(X sub E, X sub E) is equal to or less than 1, where Q sub x is a positive definite quadratic form on X. A hard bound Q sub x can be softened to many different probability distributions p sub x, but all these p sub x's carry much new information about X sub E which is absent from Q sub x, and some information which contradicts Q sub x. Both stochastic inversion (SI) and Bayesian inference (BI) estimate z from y and a soft prior bound p sub x. If that probability distribution was obtained by softening a hard prior bound Q sub x, rather than by objective statistical inference independent of y, then p sub x contains so much unsupported new information absent from Q sub x that conclusions about z obtained with SI or BI would seen to be suspect.

Backus, George E.↗

Comparing hard and soft prior bounds in geophysical inverse problems

In linear inversion of a finite-dimensional data vector y to estimate a finite-dimensional prediction vector z, prior information about X sub E is essential if y is to supply useful limits for z. The one exception occurs when all the prediction functionals are linear combinations of the data functionals. Two forms of prior information are compared: a soft bound on X sub E is a probability distribution p sub x on X which describes the observer's opinion about where X sub E is likely to be in X; a hard bound on X sub E is an inequality Q sub x(X sub E, X sub E) is equal to or less than 1, where Q sub x is a positive definite quadratic form on X. A hard bound Q sub x can be softened to many different probability distributions p sub x, but all these p sub x's carry much new information about X sub E which is absent from Q sub x, and some information which contradicts Q sub x. Both stochastic inversion (SI) and Bayesian inference (BI) estimate z from y and a soft prior bound p sub x. If that probability distribution was obtained by softening a hard prior bound Q sub x, rather than by objective statistical inference independent of y, then p sub x contains so much unsupported new information absent from Q sub x that conclusions about z obtained with SI or BI would seen to be suspect.

Backus, George E.↗

A spectral reflectance estimation technique using multispectral data from the Viking lander camera

A technique is formulated for constructing spectral reflectance curve estimates from multispectral data obtained with the Viking lander camera. The multispectral data are limited to six spectral channels in the wavelength range from 0.4 to 1.1 micrometers and most of these channels exhibit appreciable out-of-band response. The output of each channel is expressed as a linear (integral) function of the (known) solar irradiance, atmospheric transmittance, and camera spectral responsivity and the (unknown) spectral responsivity and the (unknown) spectral reflectance. This produces six equations which are used to determine the coefficients in a representation of the spectral reflectance as a linear combination of known basis functions. Natural cubic spline reflectance estimates are produced for a variety of materials that can be reasonably expected to occur on Mars. In each case the dominant reflectance features are accurately reproduced, but small period features are lost due to the limited number of channels. This technique may be a valuable aid in selecting the number of spectral channels and their responsivity shapes when designing a multispectral imaging system.

Park, S. K.↗