Search NASA⌕ Search

SEARCH · Search NASA

Results for “random variables”

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 613 records · Page 34

Heterodyne efficiency for a coherent laser radar with diffuse or aerosol targets

The performance of a Coherent Laser Radar is determined by the statistics of the coherent Doppler signal. The heterodyne efficiency is an excellent indication of performance because it is an absolute measure of beam alignment and is independent of the transmitter power, the target backscatter coefficient, the atmospheric attenuation, and the detector quantum efficiency and gain. The theoretical calculation of heterodyne efficiency for an optimal monostatic lidar with a circular aperture and Gaussian transmit laser is presented including beam misalignment in the far-field and near-field regimes. The statistical behavior of estimates of the heterodyne efficiency using a calibration hard target are considered. For space based applications, a biased estimate of heterodyne efficiency is proposed that removes the variability due to the random surface return but retains the sensitivity to misalignment. Physical insight is provided by simulation of the fields on the detector surface. The required detector calibration is also discussed.

Frehlich, R. G.↗

Approach for Uncertainty Propagation and Robust Design in CFD Using Sensitivity Derivatives

This paper presents an implementation of the approximate statistical moment method for uncertainty propagation and robust optimization for a quasi 1-D Euler CFD (computational fluid dynamics) code. Given uncertainties in statistically independent, random, normally distributed input variables, a first- and second-order statistical moment matching procedure is performed to approximate the uncertainty in the CFD output. Efficient calculation of both first- and second-order sensitivity derivatives is required. In order to assess the validity of the approximations, the moments are compared with statistical moments generated through Monte Carlo simulations. The uncertainties in the CFD input variables are also incorporated into a robust optimization procedure. For this optimization, statistical moments involving first-order sensitivity derivatives appear in the objective function and system constraints. Second-order sensitivity derivatives are used in a gradient-based search to successfully execute a robust optimization. The approximate methods used throughout the analyses are found to be valid when considering robustness about input parameter mean values.

Putko, Michele M.↗

Approach for Input Uncertainty Propagation and Robust Design in CFD Using Sensitivity Derivatives

An implementation of the approximate statistical moment method for uncertainty propagation and robust optimization for quasi 3-D Euler CFD code is presented. Given uncertainties in statistically independent, random, normally distributed input variables, first- and second-order statistical moment procedures are performed to approximate the uncertainty in the CFD output. Efficient calculation of both first- and second-order sensitivity derivatives is required. In order to assess the validity of the approximations, these moments are compared with statistical moments generated through Monte Carlo simulations. The uncertainties in the CFD input variables are also incorporated into a robust optimization procedure. For this optimization, statistical moments involving first-order sensitivity derivatives appear in the objective function and system constraints. Second-order sensitivity derivatives are used in a gradient-based search to successfully execute a robust optimization. The approximate methods used throughout the analyses are found to be valid when considering robustness about input parameter mean values.

Putko, Michele M.↗

The Net Decay Time of Anomalies in Concentrations of Atmospheric Pollutants

This paper introduces a new parameter to characterize the random component in temporal variability of atmospheric pollutants and proposes a simple statistical technique for its evaluation. That parameter is the net decay time (or the time scale) of the local anomalies in concentrations of atmospheric pollutants, rather than the traditionally used chemical lifetimes of total amounts of the species. Using widely available data of hourly multi-year surface trace gas pollutant concentrations we demonstrate a simplified way to estimate the net decay time with an exponential approximation of lag-correlation functions. We assessed the decay times of fluctuations in observations of eight atmospheric pollutants (SO2, NO, NO2, NOy, O3, CO, NH3, and HNO3) at two urban sites and one cleaner rural site in the Eastern US. The time scales of temporal fluctuations (net decay times) vary from about one hour to slightly more than one day. These scales are generally much shorter in urban environments than in remote regions. We also compared day- and night-time observations in warm and cold seasons. At night in the cold season, time scales of fluctuations in atmospheric pollutants are usually the longest. Such estimates should be useful to air quality prediction, public health, and satellite remote sensing research communities.

air quality prediction↗

Probabilistic Look-ahead Contingency Analysis Integration with Commercial Tool and Practical Data

This paper presents an initial effort of integrating a smart sampling-based probabilistic look-ahead contingency analysis algorithm with General Electric (GE) Grid Solutions’ commercial energy management system (EMS) tool as a proof-of-concept for a seamless research tool integration using real world large-scale grid data. With the increasing impact of random forces such as variable generation and load, their stochastic behaviors cannot be ignored. However, the current practices are still dominated by deterministic tools. They are becoming increasingly inadequate for the future grid. The developed look-ahead contingency analysis algorithm incorporates forecast errors of variable energy and load to address the challenges brought by the increasing uncertainty of power system. The algorithm can reveal the potential violations caused by the variance of variable energy and load that are not normally detected by traditional deterministic approaches. To test its performance under practical environments ( real data with real commercial tool), significant efforts have been made to prepare test cases, modify GE EMS tool, and adapt an extreme value distribution algorithm to analyze the GE EMS’s violation-only outputs. The test results clearly demonstrate the effectiveness of the developed algorithm as new transformer violations that were not previously detected have been identified. This performance provides better situational awareness to engineers for their decision-making process under uncertainty. Moreover, with the discussion of computational performance and future work, this paper has shown a clear path for integrating the probabilistic algorithm with commercial tools to make us better equipped for the changing power system.

Modeling and simulation of power systems, constrai↗

NSR&D Program Fiscal Year 2015 Funded Research Stochastic Modeling of Radioactive Material Releases Final Report

Nonreactor nuclear facilities operating under the approval authority of the U.S. Department of Energy use unmitigated hazard evaluations to determine if potential radiological doses associated with design basis events challenge or exceed dose evaluation guidelines. Unmitigated design basis events that sufficiently challenge dose evaluation guidelines or exceed the guidelines for members of the public or workers, merit selection of safety structures, systems, or components or other controls to prevent or mitigate the hazard. Idaho State University, in collaboration with Idaho National Laboratory, has developed a portable and simple to use software application called SODA (Stochastic Objective Decision-Aide) that stochastically calculates the radiation dose distribution associated with hypothetical radiological material release scenarios. Rather than producing a point estimate of the dose, SODA produces a dose distribution result to allow a deeper understanding of the dose potential. SODA allows users to select the distribution type and parameter values for all of the input variables used to perform the dose calculation. Users can also specify custom distributions through a user defined distribution option. SODA then randomly samples each distribution input variable and calculates the overall resulting dose distribution. In cases where an input variable distribution is unknown, a traditional single point value can be used. SODA, developed using the MATLAB coding framework, has a graphical user interface and can be installed on both Windows and Mac computers. SODA is a standalone software application and does not require MATLAB to function. SODA provides improved risk understanding leading to better informed decision making associated with establishing nuclear facility material-at-risk limits and safety structure, system, or component selection. It is important to note that SODA does not replace or compete with codes such as MACCS or RSAC; rather it is viewed as an easy to use supplemental tool to help improve risk understanding and support better informed decisions. The SODA development project was funded through a grant from the DOE Nuclear Safety Research and Development Program.

42 ENGINEERING↗

Closed‐loop hierarchical control strategies for connected and autonomous hybrid electric vehicles with random errors

This study presents a closed‐loop hierarchical control strategy for a group of connected and autonomous hybrid electric vehicles (HEVs) with the purpose of optimising fuel economy while guaranteeing traffic mobility and vehicle safety. In the hierarchical control architecture, a decentralised stochastic model predictive control (SMPC)‐based higher level controller incorporating signal phase and timing information is formulated to optimise the velocity profile of each vehicle, and an adaptive equivalent consumption minimisation strategy based lower level controller is employed for the energy management control of the HEVs. Creatively, random errors of the control variable are considered in the SMPC framework. The errors are discretised and modelled as a Markov process and the state transition matrix of the errors is generated randomly to capture the error transition dynamics. To solve the SMPC problem more efficiently, a scenario‐based SMPC is employed. Moreover, the propulsion and recuperation efficiencies of the lower level controller are calculated at each time step with measurable variables, and fed back to the higher level controller for velocity optimisation in next time step. Simulation results validate the control effectiveness and advantages of the proposed control architecture.

Qiu, Shaolin↗

Secure true random number generation using 1.5-T transistor flash memory

This disclosure relates generally to physically unclonable function (PUF) circuitry along with methods of generating numbers. In one embodiment, the PUF circuitry includes a memory, a memory control circuitry, and whitening circuitry. To reduce or eliminate the systematic bias from the array, whitening circuitry is configured to generate a random number comprising random number bits in response to the memory control circuit implementing at least one sequence of memory cycles on the array of the memory cells in the memory. The whitening circuitry is configured to provide the random number bits of the random number based on the variable bit states stored by the array of the memory cells. On average the whitening circuitry is configured to provide approximately half of the random number bits in the first bit state and half of random number bits in a second bit state.

97 MATHEMATICS AND COMPUTING↗

BROAD ABSORPTION LINE VARIABILITY ON MULTI-YEAR TIMESCALES IN A LARGE QUASAR SAMPLE

We present a detailed investigation of the variability of 428 C IV and 235 Si IV broad absorption line (BAL) troughs identified in multi-epoch observations of 291 quasars by the Sloan Digital Sky Survey-I/II/III. These observations primarily sample rest-frame timescales of 1-3.7 yr over which significant rearrangement of the BAL wind is expected. We derive a number of observational results on, e.g., the frequency of BAL variability, the velocity range over which BAL variability occurs, the primary observed form of BAL-trough variability, the dependence of BAL variability upon timescale, the frequency of BAL strengthening versus weakening, correlations between BAL variability and BAL-trough profiles, relations between C IV and Si IV BAL variability, coordinated multi-trough variability, and BAL variations as a function of quasar properties. We assess implications of these observational results for quasar winds. Our results support models where most BAL absorption is formed within an order-of-magnitude of the wind-launching radius, although a significant minority of BAL troughs may arise on larger scales. We estimate an average lifetime for a BAL trough along our line-of-sight of a few thousand years. BAL disappearance and emergence events appear to be extremes of general BAL variability, rather than being qualitatively distinct phenomena. We derive the parameters of a random-walk model for BAL EW variability, finding that this model can acceptably describe some key aspects of EW variability. The coordinated trough variability of BAL quasars with multiple troughs suggests that changes in 'shielding gas' may play a significant role in driving general BAL variability.

79 ASTRONOMY AND ASTROPHYSICS↗

A Provably Accurate Randomized Sampling Algorithm for Logistic Regression

In statistics and machine learning, logistic regression is a widely-used supervised learning technique primarily employed for binary classification tasks. When the number of observations greatly exceeds the number of predictor variables, we present a simple, randomized sampling-based algorithm for logistic regression problem that guarantees high-quality approximations to both the estimated probabilities and the overall discrepancy of the model. Our analysis builds upon two simple structural conditions that boil down to randomized matrix multiplication, a fundamental and well-understood primitive of randomized numerical linear algebra. We analyze the properties of estimated probabilities of logistic regression when leverage scores are used to sample observations, and prove that accurate approximations can be achieved with a sample whose size is much smaller than the total number of observations. To further validate our theoretical findings, we conduct comprehensive empirical evaluations. Overall, our work sheds light on the potential of using randomized sampling approaches to efficiently approximate the estimated probabilities in logistic regression, offering a practical and computationally efficient solution for large-scale datasets.

Chowdhury, Agniva↗

Supervised detection of anomalous light curves in massive astronomical catalogs

The development of synoptic sky surveys has led to a massive amount of data for which resources needed for analysis are beyond human capabilities. In order to process this information and to extract all possible knowledge, machine learning techniques become necessary. Here we present a new methodology to automatically discover unknown variable objects in large astronomical catalogs. With the aim of taking full advantage of all information we have about known objects, our method is based on a supervised algorithm. In particular, we train a random forest classifier using known variability classes of objects and obtain votes for each of the objects in the training set. We then model this voting distribution with a Bayesian network and obtain the joint voting distribution among the training objects. Consequently, an unknown object is considered as an outlier insofar it has a low joint probability. By leaving out one of the classes on the training set, we perform a validity test and show that when the random forest classifier attempts to classify unknown light curves (the class left out), it votes with an unusual distribution among the classes. This rare voting is detected by the Bayesian network and expressed as a low joint probability. Our method is suitable for exploring massive data sets given that the training process is performed offline. We tested our algorithm on 20 million light curves from the MACHO catalog and generated a list of anomalous candidates. After analysis, we divided the candidates into two main classes of outliers: artifacts and intrinsic outliers. Artifacts were principally due to air mass variation, seasonal variation, bad calibration, or instrumental errors and were consequently removed from our outlier list and added to the training set. After retraining, we selected about 4000 objects, which we passed to a post-analysis stage by performing a cross-match with all publicly available catalogs. Within these candidates we identified certain known but rare objects such as eclipsing Cepheids, blue variables, cataclysmic variables, and X-ray sources. For some outliers there was no additional information. Among them we identified three unknown variability types and a few individual outliers that will be followed up in order to perform a deeper analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

A Universal Power-law Prescription for Variability from Synthetic Images of Black Hole Accretion Flows

We present a framework for characterizing the spatiotemporal power spectrum of the variability expected from the horizon-scale emission structure around supermassive black holes, and we apply this framework to a library of general relativistic magnetohydrodynamic (GRMHD) simulations and associated general relativistic ray-traced images relevant for Event Horizon Telescope (EHT) observations of Sgr A*. We find that the variability power spectrum is generically a red-noise process in both the temporal and spatial dimensions, with the peak in power occurring on the longest timescales and largest spatial scales. When both the time-averaged source structure and the spatially integrated light-curve variability are removed, the residual power spectrum exhibits a universal broken power-law behavior. On small spatial frequencies, the residual power spectrum rises as the square of the spatial frequency and is proportional to the variance in the centroid of emission. Beyond some peak in variability power, the residual power spectrum falls as that of the time-averaged source structure, which is similar across simulations; this behavior can be naturally explained if the variability arises from a multiplicative random field that has a steeper high-frequency power-law index than that of the time-averaged source structure. We briefly explore the ability of power spectral variability studies to constrain physical parameters relevant for the GRMHD simulations, which can be scaled to provide predictions for black holes in a range of systems in the optically thin regime. We present specific expectations for the behavior of the M87* and Sgr A* accretion flows as observed by the EHT.

79 ASTRONOMY AND ASTROPHYSICS↗

Four-D global reference atmosphere technical description, part 1

An empirical atmospheric model was developed which generates values for pressure, density, temperature, and winds from surface levels to orbital altitudes. The output parameters consist of components for: (1) latitude, longitude, and altitude dependent monthly and annual means; (2) quasi-biennial oscillations; and (3) random perturbations to simulate partially the variability due to synoptic, diurnal, planetary wave, and gravity wave variations. Quasi-biennial and random variation perturbations are computed from parameters determined from various empirical studies and are added to the monthly mean values. This model has been developed as a computer program called PROFILE which can be used to generate altitude profiles of atmospheric parameters along any simulated trajectory through the atmosphere. The PROFILE program was developed for design applications in the space shuttle program. Other applications of the model are discussed, such as for global circulation and diffusion studies, and for generating profiles for comparison with other atmospheric measurement techniques, (e.g. satellite measured temperature profiles).

Justus, C. G.↗

Monte Carlo investigation of thrust imbalance of solid rocket motor pairs

The Monte Carlo method of statistical analysis is used to investigate the theoretical thrust imbalance of pairs of solid rocket motors (SRMs) firing in parallel. Sets of the significant variables are selected using a random sampling technique and the imbalance calculated for a large number of motor pairs using a simplified, but comprehensive, model of the internal ballistics. The treatment of burning surface geometry allows for the variations in the ovality and alignment of the motor case and mandrel as well as those arising from differences in the basic size dimensions and propellant properties. The analysis is used to predict the thrust-time characteristics of 130 randomly selected pairs of Titan IIIC SRMs. A statistical comparison of the results with test data for 20 pairs shows the theory underpredicts the standard deviation in maximum thrust imbalance by 20% with variability in burning times matched within 2%. The range in thrust imbalance of Space Shuttle type SRM pairs is also estimated using applicable tolerances and variabilities and a correction factor based on the Titan IIIC analysis.

Sforzini, R. H.↗

Scattering from nonspherical Chebyshev particles. III - Variability in angular scattering patterns

Shape-induced variability in the scattered intensity from randomly oriented nonspherical particles is studied. A comparison is made between bands of shape-induced variability for six size intervals within the size parameter range of 1 to 20 and corresponding spherical intensities. It is found that averaging spherical intensities over narrow size ranges produces effects that are qualitatively similar to a slightly distorted single sphere. The results also show that, for Chebyshev particles which are neither small nor large compared to the wavelength, shape-induced variability is often comparable to the mean.

Mugnai, Alberto↗

The reflection of impulses from a nonlinear random sea

A model of the reflection of radar impulses from the sea at near-vertical incidence is used to account for non-Gaussian ocean waves statistics. The joint probability density function (pdf), of wave height and slope, is calculated according to the theory of Longuet-Higgins (1963) on the distribution of variables in a 'weakly nonlinear' random era. The long-crested approximation is made, a Phillips wave spectrum is assumed, and the Gram-Charlier series is truncated after skewness terms. It is found that the height and height-slope skewness coefficients bear the ratio 1:2 and that the derived impulse response and conditional cross section versus wave height are in excellent agreement with previous observations. Finally, it is suggested that the empirically determined and theoretically predicted sea state bias be corrected for in the routine processing of satellite radar altimeter data.

Jackson, F. C.↗

Missing Data and Multiple Imputation: An Unbiased Approach

The default method of dealing with missing data in statistical analyses is to only use the complete observations (complete case analysis), which can lead to unexpected bias when data do not meet the assumption of missing completely at random (MCAR). For the assumption of MCAR to be met, missingness cannot be related to either the observed or unobserved variables. A less stringent assumption, missing at random (MAR), requires that missingness not be associated with the value of the missing variable itself, but can be associated with the other observed variables. When data are truly MAR as opposed to MCAR, the default complete case analysis method can lead to biased results. There are statistical options available to adjust for data that are MAR, including multiple imputation (MI) which is consistent and efficient at estimating effects. Multiple imputation uses informing variables to determine statistical distributions for each piece of missing data. Then multiple datasets are created by randomly drawing on the distributions for each piece of missing data. Since MI is efficient, only a limited number, usually less than 20, of imputed datasets are required to get stable estimates. Each imputed dataset is analyzed using standard statistical techniques, and then results are combined to get overall estimates of effect. A simulation study will be demonstrated to show the results of using the default complete case analysis, and MI in a linear regression of MCAR and MAR simulated data. Further, MI was successfully applied to the association study of CO2 levels and headaches when initial analysis showed there may be an underlying association between missing CO2 levels and reported headaches. Through MI, we were able to show that there is a strong association between average CO2 levels and the risk of headaches. Each unit increase in CO2 (mmHg) resulted in a doubling in the odds of reported headaches.

Foy, M.↗

VARIABILITY OF OPTICAL COUNTERPARTS IN THE CHANDRA GALACTIC BULGE SURVEY

We present optical light curves of variable stars consistent with the positions of X-ray sources identified with the Chandra X-ray Observatory for the Chandra Galactic Bulge Survey (GBS). Using data from the Mosaic-II instrument on the Blanco 4 m Telescope at CTIO, we gathered time-resolved photometric data on timescales from ∼2 hr to 8 days over the 3/4 of the X-ray survey containing sources from the initial GBS catalog. Among the light curve morphologies we identify are flickering in interacting binaries, eclipsing sources, dwarf nova outbursts, ellipsoidal variations, long period variables, spotted stars, and flare stars. Eighty-seven percent of X-ray sources have at least one potential optical counterpart. Twenty-seven percent of these candidate counterparts are detectably variable; a much greater fraction than expected for randomly selected field stars, which suggests that most of these variables are real counterparts. We discuss individual sources of interest, provide variability information on candidate counterparts, and discuss the characteristics of the variable population.

79 ASTRONOMY AND ASTROPHYSICS↗