Search NASA⌕ Search

SEARCH · Search NASA

Results for “random parameters”

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 415 records · Page 23

Accurate Magnetometer/Gyroscope Attitudes Using a Filter with Correlated Sensor Noise

Magnetometers and gyroscopes have been shown to provide very accurate attitudes for a variety of spacecraft. These results have been obtained, however, using a batch-least-squares algorithm and long periods of data. For use in onboard applications, attitudes are best determined using sequential estimators such as the Kalman filter. When a filter is used to determine attitudes using magnetometer and gyroscope data for input, the resulting accuracy is limited by both the sensor accuracies and errors inherent in the Earth magnetic field model. The Kalman filter accounts for the random component by modeling the magnetometer and gyroscope errors as white noise processes. However, even when these tuning parameters are physically realistic, the rate biases (included in the state vector) have been found to show systematic oscillations. These are attributed to the field model errors. If the gyroscope noise is sufficiently small, the tuned filter 'memory' will be long compared to the orbital period. In this case, the variations in the rate bias induced by field model errors are substantially reduced. Mistuning the filter to have a short memory time leads to strongly oscillating rate biases and increased attitude errors. To reduce the effect of the magnetic field model errors, these errors are estimated within the filter and used to correct the reference model. An exponentially-correlated noise model is used to represent the filter estimate of the systematic error. Results from several test cases using in-flight data from the Compton Gamma Ray Observatory are presented. These tests emphasize magnetometer errors, but the method is generally applicable to any sensor subject to a combination of random and systematic noise.

Sedlak, J.↗

Energy Infrastructure Futures: A Multiscale Evaluation of Projected Power Plant Siting Across the Western Interconnection

Energy Infrastructure Futures: A Multiscale Evaluation of Projected Power Plant Siting Across the Western Interconnection Description This dataset contains input and output data for the manuscript Mongird, K. et al. (under review) titled "Energy Infrastructure Futures: A Multiscale Evaluation of Projected Power Plant Siting Across the Western Interconnection". Input data corresponds to gridded spatial siting attributes that are necessary to conduct a random forest machine learning analysis of siting feature importance. Output data includes SHAP feature analysis outputs, and classification report values. For data on power plant siting results referred to in the manuscript, please refer to the CERF: IM3 Projected Western US Power Plant Locations data download page. The downloadable data includes values for eight different future scenarios for the Western US. The scenarios include combinations of two Shared Socioeconomic Pathways (SSP3 and SSP5) with four high-resolution climate projections specific to the United States (see, https://tgw-data.msdlive.org/). These climate projections include "hotter" and "cooler" variants for two Representative Concentration Pathways (RCP4.5 and RCP8.5). The resulting eight simulations are: rcp45cooler_ssp3 rcp45cooler_ssp5 rcp45hotter_ssp3 rcp45hotter_ssp5 rcp85cooler_ssp3 rcp85cooler_ssp5 rcp85hotter_ssp3 rcp85hotter_ssp5 Technical Information The dataset includes two sets of data files: (1) CERF gridded siting parameters and (2) Feature analysis outputs and classification reports. All downloadable data is in csv file format. Files with x/y coordinate information use the Albers Equal Area Conic projection (ESRI:102003). 1. CERF Gridded Siting Parameters This directory provides a balanced sample of gridded CERF siting parameters data for eight different scenarios for the Western US through 2055, seven different technologies, and eight timesteps. This data serves as input to the feature analysis. It contains the following parameters. region_name - name of region (i.e., state) sited - binary value representing whether the grid cell received a siting of that technology type (1=True) rcp - binary value representing scenario resource concentration pathway (0 = RCP4.5, 1 = RCP8.5) ssp - binary value representing scenario shared socioeconomic pathway (0 = SSP3, 1 = SSP5) climate - binary value representing cooler (0) or hotter (1) GCM forcing tech_name - generation technology name sited_year - year that values correspond to transmission_cost - cost of transmission interconnection pipeline_cost - cost of natural gas pipeline interconnection interconnection_cost - total interconnection cost (sum of transmission cost and gas pipeline cost) lmp - associated locational marginal value ($/MWh) associated with the grid cell, timestep, scenario, and technology xcoord - x-coordinate of location ycoord - y-coordinate of location 2a. Feature Analysis Output The dataset includes the feature analysis shap output for locational marginal price and interconnection cost. It contains the following parameters. technology - generator technology name scenario - name of scenario feature - name of feature, either locational_marginal_price or interconnection_cost value - the mean of absolute value of SHAP values for given feature 2b. Feature Analysis Classification Report This download includes the classification report associated with each random forest model. The dataset contains the following parameters. technology - generation technology name scenario - name of scenario test - one of precision (the proportion of predicted positives that are actually correct), recall (the proportion of actual positives that were correctly identified), f1-score (the harmonic mean of precision and recall) 0.0 - value of test for classification of 0 (grid cell not chosen for siting) 1.0 - value of test for classification of 1 (grid cell chosen for siting) accuracy - accuracy of model (i.e., fraction of all predictions that were right) macro avg - Simple average of test values for all classes weighted avg - Weighted average of test values for all classes, weighted based on Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program. License This data is made available under a CCBY4 License Disclaimer This material was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor the United States Department of Energy, nor the Contractor, nor any or their employees, nor any jurisdiction or organization that has cooperated in the development of these materials, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness or any information, apparatus, product, software, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof, or Battelle Memorial Institute. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. PACIFIC NORTHWEST NATIONAL LABORATORYoperated byBATTELLEfor theUNITED STATES DEPARTMENT OF ENERGYunder Contract DE-AC05-76RL01830

Mongird, Kendall [Pacific Northwest National Labor↗

Exploring the Connection Between Sampling Problems in Bayesian Inference and Statistical Mechanics

The Bayesian and statistical mechanical communities often share the same objective in their work - estimating and integrating probability distribution functions (pdfs) describing stochastic systems, models or processes. Frequently, these pdfs are complex functions of random variables exhibiting multiple, well separated local minima. Conventional strategies for sampling such pdfs are inefficient, sometimes leading to an apparent non-ergodic behavior. Several recently developed techniques for handling this problem have been successfully applied in statistical mechanics. In the multicanonical and Wang-Landau Monte Carlo (MC) methods, the correct pdfs are recovered from uniform sampling of the parameter space by iteratively establishing proper weighting factors connecting these distributions. Trivial generalizations allow for sampling from any chosen pdf. The closely related transition matrix method relies on estimating transition probabilities between different states. All these methods proved to generate estimates of pdfs with high statistical accuracy. In another MC technique, parallel tempering, several random walks, each corresponding to a different value of a parameter (e.g. "temperature"), are generated and occasionally exchanged using the Metropolis criterion. This method can be considered as a statistically correct version of simulated annealing. An alternative approach is to represent the set of independent variables as a Hamiltonian system. Considerab!e progress has been made in understanding how to ensure that the system obeys the equipartition theorem or, equivalently, that coupling between the variables is correctly described. Then a host of techniques developed for dynamical systems can be used. Among them, probably the most powerful is the Adaptive Biasing Force method, in which thermodynamic integration and biased sampling are combined to yield very efficient estimates of pdfs. The third class of methods deals with transitions between states described by rate constants. These problems are isomorphic with chemical kinetics problems. Recently, several efficient techniques for this purpose have been developed based on the approach originally proposed by Gillespie. Although the utility of the techniques mentioned above for Bayesian problems has not been determined, further research along these lines is warranted

Pohorille, Andrew↗

How Well do Multisatellite Products Capture the Space-Time Dynamics of Precipitation? Part II: Building an Error Model Through Spectral System Identification

Satellite precipitation products, as all quantitative estimates, come with some inherent degree of uncertainty. To associate a quantitative value of the uncertainty to each individual estimate, error modeling is necessary. Most of the error models proposed so far compute the uncertainty as a function of precipitation intensity only, and only at one specific spatio-temporal scale. We propose a spectral error model which accounts for the neighboring space-time dynamics of precipitation into the uncertainty quantification. Systematic distortions of the precipitation signal and random errors are characterized distinctively in every frequency-wavenumber band in the Fourier domain, to accurately characterize error across scales. The systematic distortions are represented as a deterministic space-time linear filtering term. The random errors are represented as a non-stationary additive noise. The spectral error model is applied to the IMERG multi satellite precipitation product and its parameters are estimated empirically through a system identification approach using the GV-MRMS gauge-radar measurements as reference (“truth”) over the eastern United States. The filtering term is found to be essentially low-pass. While traditional error models attribute most of the error variance to random errors, it is found here that the systematic filtering term explains 48% of the error variance at the native resolution of IMERG. This fact confirms that, at high resolution, filtering effects in satellite precipitation products cannot be ignored, and that the error cannot be represented as a purely random additive or multiplicative term. An important consequence is that precipitation estimates derived from totally different sources shall not be expected to automatically have statistically independent errors.

Precipitation↗

Self‐Potential Tomography Preconditioned by Particle Swarm Optimization—Application to Monitoring Hyporheic Exchange in a Bedrock River

Abstract A self‐potential (SP) data‐inversion algorithm was developed and tested on an analytical model of electrical‐potential profile data attributed to single and multiple polarized electrical sources. The developed algorithm was then validated by an application to SP‐monitoring field data measured on the floodplain of East Fork Poplar Creek, Oak Ridge, Tennessee, to image electrical sources in areas conducive to preferential flow into the flood plain from the bedrock‐lined riverbed. The algorithm combined stochastic source‐localization by particle‐swarm‐optimization (PSO) of electrical sources characterized by simplified geometries with source tomography by regularized weighted least‐squares minimization of a quadratic objective function. Prior information was incorporated by preconditioning the tomography algorithm by PSO results. Variable percentages of random noise were added to analytical‐model data to evaluate the algorithm performance. Results indicated that true parameters of single‐source models were inverted and approximated with small residual error, whereas inversion of analytical‐model data representing multiple electrical sources accurately approximated the locations of the sources but miscalculated some parameters because of the non‐uniqueness of the inverse‐model solution. Source tomography applied to analytical model data during testing produced a spatially continuous parameter field that identified the locations of point‐scale synthetic dipole sources of electrical current flow with varying degrees of accuracy depending on the prior information incorporated into the tomography. When applied to SP‐monitoring field data, the algorithm imaged electrical sources within a known fault that intersects the bedrock riverbed and flood plain of East Fork Poplar Creek and depicted dynamic electrical conditions attributed to hyporheic exchange.

54 ENVIRONMENTAL SCIENCES↗

The design of a device for hearer and feeler differentiation, part A

A speech modulated white noise device is reported that gives the rhythmic characteristics of a speech signal for intelligible reception by deaf persons. The signal is composed of random amplitudes and frequencies as modulated by the speech envelope characteristics of rhythm and stress. Time intensity parameters of speech are conveyed through the vibro-tactile sensation stimuli.

Creecy, R.↗

Methodology for error analysis and simulation of lidar aerosol measurements

A methodology is presented for objective and automated determination of the uncertainty in lidar aerosol measurements. This methodology is based on standard error-propagation procedures, a large data base on atmospheric behavior, and long experience in lidar data processing. Algebraic expressions for probable error are derived as a function of the relevant parameters. The validity of these expressions is then tested by making simulated measurements and analyses in which random errors of appropriate size are injected at proper steps of the measurement and analysis process. An illustrative example is given where the methodology is applied to a new lidar system now being used for airborne measurements of the stratospheric aerosol.

Russell, P. B.↗

The Analysis of Spaced Sensor Records by Correlation Techniques

The use of cross correlation and autocorrelation techniques to determine drift velocity and other parameters of moving patterns is examined. Background on the nature of moving patterns is provided. Stationary randomly-changing patterns and moving randomly-changing patterns are considered.

Briggs, B. H.↗

The D/H ratio in Saturn's atmosphere from Voyager IRIS spectra

A selection of 1656 spectra recorded between 180 and 2300 mm in the 20 to 40 N latitude range by the Voyager interferometer IRIS were analyzed. The CH3D/H2 and CH4/H2 ratios are determined from a best fit of the radiances measured in the nu6 and nu4 fundamental bands, respectively centered at 1161 and 1304 mm. Results are CH3D/H2 = 3.9 + or - 2.5 times 10 to the minus 7th power and CH4/H2 = 4.5 + 2.4, -1.9 x 0.001. The uncertainty includes random errors due to instrumental noise and possible systematic errors in the assumed temperature profile and spectroscopic parameters. The D/H isotopic ratio is determined from abundance ratios from the expression D/H = 1/4f (CH3D/H2)/(CH4/H2) where f accounts for equilibrium deuterium fractionation between the molecular hydrogen and methane phases. Beer and Taylor (1973, 1978) estimated f = 1.37 + or - 0.07 for Jupiter assuming a wide range of convective velocities and neglecting catalytic effects from dust particles. The same value was adopted for Saturn and yields D/H = 1.6 +1.3, -1.2 times 10 to the minus 5th power.

Courtin, R.↗

On-Orbit System Identification

Information derived from accelerometer readings benefits important engineering and control functions. Report discusses methodology for detection, identification, and analysis of motions within space station. Techniques of vibration and rotation analyses, control theory, statistics, filter theory, and transform methods integrated to form system for generating models and model parameters that characterize total motion of complicated space station, with respect to both control-induced and random mechanical disturbances.

Mettler, E.↗

Turbulent and directed plasma motions in solar flares

An improved method for fitting asymmetric soft X-ray line profiles from solar flares is presented. A two-component model is used where one component represents the total emission from directed upflow plasma and the other the emission from the plasma at rest. Unlike previous methods, the width of the moving component is independent from that of the stationary component. Time variations of flare plasma characteristics (i.e., temperature, emission measure of moving and stationary plasma, upflow and turbulent velocities) are derived from the Ca XIX and Fe XXV spectra recorded by the Bent Crystal Spectrometer on the Solar Maximum Mission. The fitting technique provides a statistical estimation for the uncertainties in the fitting parameters. The relationship between the directed and turbulent motions has been studied, and a correlation of the random and directed motions has been found in some flares with intensive plasma upflows. Mean temperatures of the upflowing and stationary plasmas are compared for the first time from ratios of calcium to iron X-ray line intensities. Finally, evidence for turbulent motions and the possibility of plasma upflow late into the decay phase is presented and discussed.

Fludra, A.↗

On the RNG theory of turbulence

The Yakhot and Orszag (1986) renormalization group (RNG) theory of turbulence has generated a number of scaling law constants in reasonable quantitative agreement with experiments. The theory itself is highly mathematical, and its assumptions and approximations are not easily appreciated. The present paper reviews the RNG theory and recasts it in more conventional terms using a distinctly different viewpoint. A new formulation based on an alternative interpretation of the origin of the random force is presented, showing that the artificially introduced epsilon in the original theory is an adjustable parameter, thus offering a plausible explanation for the remarkable record of quantitative success of the so-called epsilon-expansion procedure.

Lam, S. H.↗

Radiative transfer modelling of discontinuous tree canopies at microwave frequencies

Development of a first-order radiative transfer model for predicting backscatter from tree canopies has been underway at the University of Michigan Radiation Laboratory for some time. This model is known as the Michigan Microwave Canopy Scattering (MIMICS) model. This article presents the second-generation MIMICS model (MIMICS II) which accounts for canopies with discontinuous (open) crown layer geometries. MIMICS II models open crown layers by treating the location, size, and shape of the individual tree crowns as random variables. The backscattering coefficients for the canopy are then determined by introducing statistics derived from these parameters into the radiative transfer solution. Application of the radiative transfer equations to the discontinuous canopy geometry is presented. The resulting model is a robust fully polarimetric solution that is applicable over a wide variety of canopy architectures. Model simulations are compared to results generated with the continuous canopy model. The effect of the open crown geometry is found to be most significant at shallow incidence angles and at high frequencies for trees with well-developed crowns. Under these conditions, the gaps in the crown layer give rise to a notable increase in crown layer transmissivity which allows the radar to see through to the lower layers of the canopy more easily, thereby directly affecting the backscatter contribution of the trunks and ground.

Mcdonald, K. C.↗

The Dripping Handrail Model: Transient Chaos in Accretion Systems

We define and study a simple dynamical model for accretion systems, the "dripping handrail" (DHR). The time evolution of this spatially extended system is a mixture of periodic and apparently random (but actually deterministic) behavior. The nature of this mixture depends on the values of its physical parameters - the accretion rate, diffusion coefficient, and density threshold. The aperiodic component is a special kind of deterministic chaos called transient chaos. The model can simultaneously exhibit both the quasiperiodic oscillations (QPO) and very low frequency noise (VLFN) that characterize the power spectra of fluctuations of several classes of accretion systems in astronomy. For this reason, our model may be relevant to many such astrophysical systems, including binary stars with accretion onto a compact object - white dwarf, neutron star, or black hole - as well as active galactic nuclei. We describe the systematics of the DHR's temporal behavior, by exploring its physical parameter space using several diagnostics: power spectra, wavelet "scalegrams," and Lyapunov exponents. In addition, we note that for large accretion rates the DHR has periodic modes; the effective pulse shapes for these modes - evaluated by folding the time series at the known period - bear a resemblance to the similarly- determined shapes for some x-ray pulsars. The pulsing observed in some of these systems may be such periodic-mode accretion, and not due to pure rotation as in the standard pulsar model.

Young, Karl↗

Remote Sensing/gis Integration for Site Planning and Resource Management

The development of an interactive/batch gridded information system (array of cells georeferenced to USGS quad sheets) and interfacing application programs (e.g., hydrologic models) is discussed. This system allows non-programer users to request any data set(s) stored in the data base by inputing any random polygon's (watershed, political zone) boundary points. The data base information contained within this polygon can be used to produce maps, statistics, and define model parameters for the area. Present/proposed conditions for the area may be compared by inputing future usage (land cover, soils, slope, etc.). This system, known as the Hydrologic Analysis Program (HAP), is especially effective in the real time analysis of proposed land cover changes on runoff hydrographs and graphics/statistics resource inventories of random study area/watersheds.

Fellows, J. D.↗

Discrete Random Media Techniques for Microwave Modeling of Vegetated Terrain

Microwave remote sensing of agricultural crops and forested regions is studied. Long term goals of the research involve modeling vegetation so that radar signatures can be used to infer the parameters which characterize the vegetation and underlying ground. Vegetation is modeled by discrete scatterers viz, leaves, stems, branches and trunks. These are replaced by glossy dielectric discs and cylinders. Rough surfaces are represented by their mean and spectral characteristics. Average scattered power is then calculated by employing discrete random media methodology such as the distorted Born approximation or transport theory. Both coherent and incoherent multiple scattering techniques are explored. Once direct methods are developed, inversion techniques can be investigated.

Lang, R. H.↗

On the VLBI-satellite laser ranging Iron Triangle intercomparison experiment

Simulations were performed to determine the suitability of the proposed station locations. The criterion was a comparison among the possible station configuration of the standard deviations of baseline and Earth rotation parameters estimated from a least squares covariance analysis. Only the relative magnitudes of the standard deviations were addressed in the analysis. Only random errors were assumed and no provision was made for systematic effects. The Iron Triangle, consisting of the stations at Westford (Haystack), Massachusetts, Owens Valley, California, and Ft. Davis, Texas, was regarded as the structure of the proposed network with options to incorporate either the Goldstone, California, or Green Bank, West Virginia, station or both. It was decided to include the Richmond, Florida, station in the analysis since it offered more North-South separation and therefore could strengthen the geometry of the network especially in the recovery of Earth rotation parameters.

Bock, Y.↗

Picasso: Memory-Efficient Graph Coloring Using Palettes With Applications in Quantum Computing

A coloring of a graph is an assignment of colors to vertices such that no two neighboring vertices have the same color. The need for memory-efficient coloring algorithms is motivated by their application in computing clique partitions of graphs arising in quantum computations where the objective is to map a large set of Pauli strings into a compact set of unitaries. We present Picasso, a randomized memory-efficient iterative parallel graph coloring algorithm with theoretical sublinear space guarantees under practical assumptions. The parameters of our algorithm provide a trade-off between coloring quality and resource consumption. To assist the user, we also propose a machine learning model to predict the coloring algorithm’s parameters considering these trade-offs. We provide a sequential and a parallel implementation of the proposed algorithm. We perform an experimental evaluation on a 64-core AMD CPU equipped with 512 GB of memory and an Nvidia A100 GPU with 40GB of memory. For a small dataset where existing coloring algorithms can be executed within the 512 GB memory budget, we show up to 68× memory savings. On massive datasets we demonstrate that GPU-accelerated Picasso can process inputs with 49.5× more Pauli strings (vertex set in our graph) and 2,478× more edges than state-of-the-art parallel approaches.

artificial intelligence, quantum computing↗