Search NASA⌕ Search

SEARCH · Search NASA

Results for “estimation steps”

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 307 records · Page 17

Precision in the perception of direction of a moving pattern

The precision of the model of pattern motion analysis put forth by Adelson and Movshon (1982) who proposed that humans determine the direction of a moving plaid (the sum of two sinusoidal gratings of different orientations) in two steps is qualitatively examined. The volocities of the grating components are first estimated, then combined using the intersection of constraints to determine the velocity of the plaid as a whole. Under the additional assumption that the noise sources for the component velocities are independent, an approximate expression can be derived for the precision in plaid direction as a function of the precision in the speed and direction of the components. Monte Carlo simulations verify that the expression is valid to within 5 percent over the natural range of the parameters. The expression is then used to predict human performance based on available estimates of human precision in the judgment of single component speed. Human performance is predicted to deteriorate by a factor of 3 as half the angle between the wavefronts (theta) decreases from 60 to 30 deg, but actual performance does not. The mean direction discrimination for three human observers was 4.3 plus or minus 0.9 deg (SD) for theta = 60 deg and 5.9 plus or minus 1.2 for theta = 30 deg. This discrepancy can be resolved in two ways. If the noises in the internal representations of the component speeds are smaller than the available estimates or if these noises are not independent, then the psychophysical results are consistent with the Adelson-Movshon hypothesis.

Stone, Leland S.↗

Galerkin/Runge-Kutta discretizations for parabolic equations with time-dependent coefficients

A new class of fully discrete Galerkin/Runge-Kutta methods is constructed and analyzed for linear parabolic initial boundary value problems with time dependent coefficients. Unlike any classical counterpart, this class offers arbitrarily high order convergence while significantly avoiding what has been called order reduction. In support of this claim, error estimates are proved, and computational results are presented. Additionally, since the time stepping equations involve coefficient matrices changing at each time step, a preconditioned iterative technique is used to solve the linear systems only approximately. Nevertheless, the resulting algorithm is shown to preserve the original convergence rate while using only the order of work required by the base scheme applied to a linear parabolic problem with time independent coefficients. Furthermore, it is noted that special Runge-Kutta methods allow computations to be performed in parallel so that the final execution time can be reduced to that of a low order method.

Keeling, Stephen L.↗

On adaptive weighted polynomial preconditioning for Hermitian positive definite matrices

The conjugate gradient algorithm for solving Hermitian positive definite linear systems is usually combined with preconditioning in order to speed up convergence. In recent years, there has been a revival of polynomial preconditioning, motivated by the attractive features of the method on modern architectures. Standard techniques for choosing the preconditioning polynomial are based only on bounds for the extreme eigenvalues. Here a different approach is proposed, which aims at adapting the preconditioner to the eigenvalue distribution of the coefficient matrix. The technique is based on the observation that good estimates for the eigenvalue distribution can be derived after only a few steps of the Lanczos process. This information is then used to construct a weight function for a suitable Chebyshev approximation problem. The solution of this problem yields the polynomial preconditioner. In particular, we investigate the use of Bernstein-Szego weights.

Fischer, Bernd↗

A Hierarchical Approach to Buckling Load Calculations

The advantages of using a hierarchical analysis approach to calculate the buckling load of an axially compressed composite cylindrical shell is demonstrated using an example taken from a recent experimental program. The Delft Interactive Shell DEsign COde (DISDECO) shell design code is used for this hierarchical analysis approach to provide an accurate prediction of the critical buckling load of the actual shell structure. DISDECO includes the influence of the boundary conditions, initial geometric imperfections, the effects of stiffener and load eccentricities, and the effects of prebuckling deformations caused by edge constraints in the analysis. It is shown that the use of DISDECO makes it relatively simple to proceed step by step from simple to more complex models and solution procedures. As a final step in the hierarchical analysis approach, the critical buckling load and the estimated imperfection sensitivity of the shell are verified by conducting an analysis of a large finite element model with one of the current generation two-dimensional shell analysis codes with advanced capabilities needed to represent both geometric and material nonlinearities.

Arbocz, Johann↗

On a High-Fidelity Hierarchical Approach to Buckling Load Calculations

As a step towards developing a new design philosophy, one that moves away from the traditional empirical approach used today in design towards a science-based design technology approach, a recent test series of 5 composite shells carried out by Waters at NASA Langley Research Center is used. It is shown how the hierarchical approach to buckling load calculations proposed by Arbocz et al can be used to perform an approach often called "high fidelity analysis", where the uncertainties involved in a design are simulated by refined and accurate numerical methods. The Delft Interactive Shell DEsign COde (short, DISDECO) is employed for this hierarchical analysis to provide an accurate prediction of the critical buckling load of the given shell structure. This value is used later as a reference to establish the accuracy of the Level-3 buckling load predictions. As a final step in the hierarchical analysis approach, the critical buckling load and the estimated imperfection sensitivity of the shell are verified by conducting an analysis using a sufficiently refined finite element model with one of the current generation two-dimensional shell analysis codes with the advanced capabilities needed to represent both geometric and material nonlinearities.

Arbocz, Johann↗

Buckling Load Calculations of the Isotropic Shell A-8 Using a High-Fidelity Hierarchical Approach

As a step towards developing a new design philosophy, one that moves away from the traditional empirical approach used today in design towards a science-based design technology approach, a test series of 7 isotropic shells carried out by Aristocrat and Babcock at Caltech is used. It is shown how the hierarchical approach to buckling load calculations proposed by Arbocz et al can be used to perform an approach often called 'high fidelity analysis', where the uncertainties involved in a design are simulated by refined and accurate numerical methods. The Delft Interactive Shell DEsign COde (short, DISDECO) is employed for this hierarchical analysis to provide an accurate prediction of the critical buckling load of the given shell structure. This value is used later as a reference to establish the accuracy of the Level-3 buckling load predictions. As a final step in the hierarchical analysis approach, the critical buckling load and the estimated imperfection sensitivity of the shell are verified by conducting an analysis using a sufficiently refined finite element model with one of the current generation two-dimensional shell analysis codes with the advanced capabilities needed to represent both geometric and material nonlinearities.

Arbocz, Johann↗

Towards a Global Water Scarcity Risk Assessment Framework: Incorporation of Probability Distributions and Hydro-Climatic Variability

Changing hydro-climatic and socioeconomic conditions increasingly put pressure on fresh water resources and are expected to aggravate water scarcity conditions towards the future. Despite numerous calls for risk-based water scarcity assessments, a global-scale framework that includes UNISDR's definition of risk does not yet exist. This study provides a first step towards such a risk based assessment, applying a Gamma distribution to estimate water scarcity conditions at the global scale under historic and future conditions, using multiple climate change and population growth scenarios. Our study highlights that water scarcity risk, expressed in terms of expected annual exposed population, increases given all future scenarios, up to greater than 56.2% of the global population in 2080. Looking at the drivers of risk, we find that population growth outweigh the impacts of climate change at global and regional scales. Using a risk-based method to assess water scarcity, we show the results to be less sensitive than traditional water scarcity assessments to the use of fixed threshold to represent different levels of water scarcity. This becomes especially important when moving from global to local scales, whereby deviations increase up to 50% of estimated risk levels.

Populations↗

Astrophysics and Technical Study of a Solar Neutrino Spacecraft: NIAC (NASA Innovative Advanced Concepts) Grant - NASA 2018-Phase-1

We report on our study of the design of a neutrino detector, shielding and veto array needed to operate a neutrino detector in space close to the Sun. This study also took into account the expected rates of Galactic gamma and cosmic rays in addition to the particles from the Sun.These preliminary studies show that we can devise a detector such that a small signal of neutrino interactions can be extracted from a large random number of events from the background sources using a double timing method from the conversion electron produced in the neutrino interaction and a secondary delayed signal from the nuclear excited state produced from the initial neutrino interaction; in our case the conversion of Ga 69 or 71 into Ge 69 or 71, but this method could apply to other nuclei with large neutrino cross sections such as Ir 115. Although these types of events need to be above 0.405 megaelectronvolt (MeV) neutrino energy and are only 66 percent of all conversion neutrino interactions on Gallium, this is a small price to pay for an increase of 10,000 by going close to the Sun to enhance the neutrino rate over the background combatorical fake-signal events. The conclusion of this Phase-1 study is very positive in that we can get the backgrounds less than 20 percent fake signals, and in addition to this we have devised another shielding method that makes the Galactic gamma-ray rate a hundred fold less which will make further improvements over these initial estimates. Although these studies are very encouraging it suggests that the next step is a NIAC Phase-II to actually build a test device,measuring basic principles such as light attention within the scintillator with high dopants and to take data in the lab with a cosmic-ray test stand and triggered X-ray source for comparison with simulated expected performance of the detector. This would be the perfect lead into a future proposal beyond a NIAC (NASA Innovative Advanced Concepts) Phase-II for a test flight of a small one-pint detector in orbit of the detector concept beyond Earth outside of the radiation belts.

Astrophysics↗

Residual Stresses Induced due to Curing of the Bulk Matrix in a Simplified Three-Dimensional (3D) Woven Repeating Unit Cell

It has been observed that 3D woven materials are highly susceptible to processes induced defects including tow misalignment, intratow and intertow cracking and/or voids. Adequate process modeling of 3D woven Polymer Matrix Composites (PMCs) is necessary to predict and estimate the effects of the manufacturing process on these defects. A preliminary step towards this is achieved with the commercial Finite Element Analysis (FEA) tool Abaqus complemented with written user subroutines to account for the effect of shrinkage and thermoelastic properties evolution as a function of the degree of cure of the bulk matrix phase in a 3D woven composite repeating unit cell (RUC). As a pathfinder, the process modeling framework is demonstrated on a simplified 3D woven RUC of an AS4/RIMR 135 system. Curing simulations are performed on the 3D woven RUC to observe the processing induced residual stresses. The results of the analysis show that high stress concentrations, present within the 3D intertow matrix bordering the binder tow, coincide with X-Ray computed tomography (X-Ray-CT) scan data depicting cracking in the same location

3D Woven↗

Asymmetry in Subseasonal Surface Air Temperature Forecast Error with Respect to Soil Moisture Initialization

Soil moisture (W) helps control evapotranspiration (ET), and ET variations can in turn have a distinct impact on 2-m air temperature (T2M), given that increases in evaporative cooling encourage reduced temperatures. Soil moisture is accordingly linked to T2M, and realistic soil moisture initialization has, in previous studies, been shown to improve the skill of subseasonal T2M forecasts. The relationship between soil moisture and evapotranspiration, however, is distinctly nonlinear, with ET tending to increase with soil moisture in drier conditions and to be insensitive to soil moisture variations in wetter conditions. Here, through an extensive analysis of subseasonal forecasts produced with a state-of-the-art seasonal forecast system, this nonlinearity is shown to imprint itself on T2M forecast error in the conterminous United States in two unique ways: (i) the T2M forecast bias (relative to independent observations) induced by a negative precipitation bias tends to be larger for dry initializations, and (ii) on average, the unbiased root-mean-square error (ubRMSE) tends to be larger for dry initializations. Such findings can aid in the identification of forecasts of opportunity; taken a step further, they suggest a pathway for improving bias correction and uncertainty estimation in subseasonal T2M forecasts by conditioning each on initial soil moisture state.

Air Temperature↗

T2M Forecasts at Subseasonal Leads: Do Different Soil Moisture Initial States Have Different Impacts?

Soil moisture (W) helps control evapotranspiration (ET), and ET variations can in turn have a distinct impact on 2-m air temperature (T2M), given that increases in evaporative cooling encourage reduced temperatures. Soil moisture is accordingly linked to T2M, and realistic soil moisture initialization has, in previous studies, been shown to improve the skill of subseasonal T2M forecasts. The relationship between soil moisture and evapotranspiration, however, is distinctly nonlinear, with ET tending to increase with soil moisture in drier conditions and to be insensitive to soil moisture variations in wetter conditions. Here, through an extensive analysis of subseasonal forecasts produced with a state-of-the-art seasonal forecast system, this nonlinearity is shown to imprint itself on T2M forecast error in the conterminous United States in two unique ways: (i) the T2M forecast bias (relative to independent observations) induced by a negative precipitation bias tends to be larger for dry initializations, and (ii) on average, the unbiased root-mean-square error (ubRMSE) tends to be larger for dry initializations. Such findings can aid in the identification of forecasts of opportunity; taken a step further, they suggest a pathway for improving bias correction and uncertainty estimation in subseasonal T2M forecasts by conditioning each on initial soil moisture state.

2-meter Temperature↗

Towards an Optimal Noise Versus Resolution Trade-Off in Wind Scatterometry

A scatterometer is a radar that measures the normalized radar cross section sigma(sup 0) of the Earth's surface. Over the ocean this signal is related to the wind via the geophysical model function (GMF). The objective of wind scatterometry is to estimate the wind vector field from sigma(sup 0) measurements; however, there are many subtleties that complicate this problem-making it difficult to obtain a unique wind field estimate. Conventionally, wind estimation is split into two stages: a wind retrieval stage in which several ambiguous solutions are obtained, and an ambiguity removal stage in which ambiguities are chosen to produce an appropriate wind vector field estimate. The most common approach to wind field estimation is to grid the scatterometer swath into wind vector cells and estimate wind vector ambiguities independently for each cell. Then, field wise structure is imposed on the solution by an ambiguity selection routine. Although this approach is simple and practical, it neglects field wise structure in the retrieval step and does not account for the spatial correlation imposed by the sampling. This makes it difficult to develop a theoretically appropriate noise versus resolution trade-off using pointwise retrieval. Fieldwise structure may be imposed in the retrieval step using a model-based approach. However, this approach is generally only practical if a low order wind field model is applied, which may discard more information than is desired. Furthermore, model-based approaches do not account for the structure imposed by the sampling. A more general fieldwise approach is to estimate all the wind vectors for all the WVCs simultaneously from all the measurements. This approach can account for structure of the wind field as well as structure imposed by the sampling in the wind retrieval step. Williams and Long in 2010 developed a fieldwise retrieval method based on maximum a posteriori estimation (MAP). This MAP approach can be extended to perform a noise versus resolution trade-off, and deal with ambiguity selection. This paper extends the fieldwise MAP estimation approach and investigates both the noise versus resolution trade-off as well as ambiguity removal in the fieldwise wind retrieval step. The method is then applied to the Sea Winds scatterometer and the results are analyzed. This paper extends the fieldwise MAP estimation approach and investigates both the noise versus resolution trade-off as well as ambiguity removal in the fieldwise wind retrieval step. The method is then applied to the Sea Winds scatterometer and the results are analyzed.

soil moisture↗

Digital image processing of vascular angiograms

The paper discusses the estimation of the degree of atherosclerosis in the human femoral artery through the use of a digital image processing system for vascular angiograms. The film digitizer uses an electronic image dissector camera to scan the angiogram and convert the recorded optical density information into a numerical format. Another processing step involves locating the vessel edges from the digital image. The computer has been programmed to estimate vessel abnormality through a series of measurements, some derived primarily from the vessel edge information and others from optical density variations within the lumen shadow. These measurements are combined into an atherosclerosis index, which is found in a post-mortem study to correlate well with both visual and chemical estimates of atherosclerotic disease.

Selzer, R. H.↗

Analysis of Pulsed Lidar Measurements of Atmospheric CO2 Column Absorption During the ASCENDS 2009-2011 Airborne Campaigns

We have developed a pulsed lidar technique for measuring the tropospheric CO2 concentrations as a candidate for NASA's ASCENDS mission and have demonstrated the CO2 and O2 measurements from aircraft. Our technique uses two pulsed lasers allowing simultaneous measurement of a single CO2 absorption line near 1572 nm, O2 extinction in the Oxygen A-band, surface height and backscatter profile. The lasers are stepped in wavelength across the CO2 line and an O2 line doublet during the measurement. The column densities for the CO2 and O2 are estimated from the differential optical depths (DOD) of the scanned absorption lines via the IPDA technique. For the 2009 ASCENDS campaign we flew the CO2 lidar on a Lear-25 aircraft, and measured the absorption line shapes of the CO2 line using 20 wavelength samples per scan. Measurements were made at stepped altitudes from 3 to 12.6 km over the Lamont OK, central Illinois, North Carolina, and over the Virginia Eastern Shore. Although the received signal energies were weaker than expected for ASCENDS, clear CO2 line shapes were observed at all altitudes. Most flights had 5-6 altitude steps with 200-300 seconds of recorded measurements per step. We averaged every 10 seconds of measurements and used a cross-correlation approach to estimate the range to the scattering surface and the echo pulse energy at each wavelength. We then solved for the best-fit CO2 absorption line shape, and calculated the DOD of the fitted CO2 line, and computed its statistics at the various altitude steps. We compared them to CO2 optical depths calculated from spectroscopy based on HITRAN 2008 and the column number densities calculated from the airborne in-situ readings. The 2009 measurements have been analyzed and they were similar on all flights. The results show clear CO2 line shape and absorption signals, which follow the expected changes with aircraft altitude from 3 to 13 km. They showed the expected nearly the linear dependence of DOD vs altitude. The measurements showed ~ 1 ppm random errors for 8-10 km altitudes and ~30 sec averaging times. For the 2010 ASCENDS campaigns we flew the CO2lidar on the NASA DC-8 and added an 02lidar channel. During July 2010 we made measurements of CO2 and O2 column absorption during longer flights over Railroad Valley NV, the Pacific Ocean and over Lamont OK. CO2 measurements were made with 30 steps/scan, 300 scans/sec and improved line resolution and receiver sensitivity. Analysis of the 2010 CO2 measurements shows the expected ~linear change of DOD with altitude. For measurements at altitudes> 6 km the random errors were 0.3 ppm for 80 sec averaging times. For the summer 2011 ASCENDS campaigns we made further improvements to the lidar's CO2 line scan and receiver sensitivity. We demonstrated measurements over the California Central Valley, to stratus cloud tops over the Pacific Ocean, over mountain regions with snow, and over several areas with broken clouds. Details of the lidar measurements and their analysis will be described in the presentation.

Abshire, J. B.↗

Pulsed Lidar Measurements of Atmospheric CO2 Column Absorption and Range During the ASCENDS 2009-2011 Airborne Campaigns

We have developed a pulsed lidar technique for measuring the tropospheric CO2 concentrations as a candidate for NASA's ASCENDS mission and have demonstrated the CO2 and O2 measurements from aircraft. Our technique uses two pulsed lasers allowing simultaneous measurement of a single CO2 absorption line near 1572 nm, O2 extinction in the Oxygen A-band, surface height and backscatter profile. The lasers are stepped in wavelength across the CO2 line and an O2 line doublet during the measurement. The column densities for the CO2 and O2 are estimated from the differential optical depths (DOD) of the scanned absorption lines via the IPDA technique. For the 2009 ASCENDS campaign we flew the CO2 lidar only on a Lear-25 aircraft, and measured the absorption line shapes of the CO2 line using 20 wavelength samples per scan. Measurements were made at stepped altitudes from 3 to 12.6 km over the Lamont OK, central Illinois, North Carolina, and over the Virginia Eastern Shore. Although the received signal energies were weaker than expected for ASCENDS, clear C02 line shapes were observed at all altitudes. Most flights had 5-6 altitude steps with 200-300 seconds of recorded measurements per step. We averaged every 10 seconds of measurements and used a cross-correlation approach to estimate the range to the scattering surface and the echo pulse energy at each wavelength. We then solved for the best-fit CO2 absorption line shape, and calculated the DOD of the fitted CO2 line, and computed its statistics at the various altitude steps. We compared them to CO2 optical depths calculated from spectroscopy based on HITRAN 2008 and the column number densities calculated from the airborne in-situ readings. The 2009 measurements have been analyzed in detail and they were similar on all flights. The results show clear CO2 line shape and absorption signals, which follow the expected changes with aircraft altitude from 3 to 13 km. They showed the expected nearly the linear dependence of DOD vs altitude. The measurements showed -1 ppm random errors for 8-10 km altitudes and -30 sec averaging times. For the 2010 ASCENDS campaigns we flew the CO2 lidar on the NASA DC-8 and added an O2 lidar channel. During July 2010 we made measurements of CO2 and O2 column absorption during longer flights over Railroad Valley NV, the Pacific Ocean and over Lamont OK. CO2 measurements were made with 30 steps/scan, 300 scans/sec and improved line resolution and receiver sensitivity. Analysis of the 2010 CO2 measurements shows the expected -linear change of DOD with altitude. For measurements at altitudes> 6 km the random errors were 0.3 ppm for 80 sec averaging times. For the summer 2011 ASCENDS campaigns we made further improvements to the lidar's CO2 line scan and receiver sensitivity. The seven flights in the 2011 Ascends campaign were flown over a wide variety of surface and cloud conditions in the US, which produced a wide variety of lidar signal conditions. Details of the lidar measurements and their analysis will be described in the presentation.

Abshire, J. B.↗

ICESat Observations of Arctic Sea Ice: A First Look

Analysis of near-coincident ICESat and RADARSAT imagery shows that the retrieved elevations from the laser altimeter are sensitive to new openings (containing thin ice or open water) in the sea ice cover as well as to surface relief of old and first-year ice. The precision of the elevation estimates, measured over relatively flat sea ice, is approx. 2 cm. Using the thickness of thin-ice in recent openings to estimate sea level references, we obtain the sea-ice freeboard along the altimeter tracks. This step is necessitated by the large uncertainties in the sea surface topography compared to that required for accurate determination of freeboard. Unknown snow depth introduces the largest uncertainty in the conversion of freeboard to ice thickness. Surface roughness is also derived, for the first time, from the variability of successive elevation estimates along the altimeter track. Overall, these ICESat measurements provide an unprecedented view of the Arctic Ocean ice cover at length scales at and above the spatial dimension of the altimeter footprint of approx. 70 m.

Arctic Ocean↗

Real-time Detection of Moving Objects from Moving Vehicles Using Dense Stereo and Optical Flow

Dynamic scene perception is very important for autonomous vehicles operating around other moving vehicles and humans. Most work on real-time object tracking from moving platforms has used sparse features or assumed flat scene structures. We have recently extended a real-time. dense stereo system to include realtime. dense optical flow, enabling more comprehensive dynamic scene analysis. We describe algorithms to robustly estimate 6-DOF robot egomotion in the presence of moving objects using dense flow and dense stereo. We then use dense stereo and egomotion estimates to identify other moving objects while the robot itself is moving. We present results showing accurate egomotion estimation and detection of moving people and vehicles under general 6DOF motion of the robot and independently moving objects. The system runs at 18.3 Hz on a 1.4 GHz Pentium M laptop. computing 160x120 disparity maps and optical flow fields, egomotion, and moving object segmentation. We believe this is a significant step toward general unconstrained dynamic scene analysis for mobile robots, as well as for improved position estimation where GPS is unavailable.

moving object detection↗

Potential Predictability of U.S. Summer Climate with "Perfect" Soil Moisture

The potential predictability of surface-air temperature and precipitation over the United States continent was assessed for a GCM forced by observed sea surface temperatures and an estimate of observed ground soil moisture contents. The latter was obtained by substituting the GCM simulated precipitation, which is used to drive the GCM's land-surface component, with observed pentad-mean precipitation at each time step of the model's integration. With this substitution, the simulated soil moisture correlates well with an independent estimate of observed soil moisture in all seasons over the entire US continent. Significant enhancements on the predictability of surface-air temperature and precipitation were found in boreal late spring and summer over the US continent. Anomalous pattern correlations of precipitation and surface-air temperature over the US continent in the June-July-August season averaged for the 1979-2000 period increased from 0.01 and 0.06 for the GCM simulations without precipitation substitution to 0.23 and 0.3 1, respectively, for the simulations with precipitation substitution. Results provide an estimate for the limits of potential predictability if soil moisture variability is to be perfectly predicted. However, this estimate may be model dependent, and needs to be substantiated by other modeling groups.

Yang, Fanglin↗