Search NASA⌕ Search

SEARCH · Search NASA

Results for “Randomized methods”

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,441 records · Page 80

Experimental Design for the Evaluation of Detection Techniques of Hidden Corrosion Beneath the Thermal Protective System of the Space Shuttle Orbiter

The detection of corrosion beneath Space Shuttle Orbiter thermal protective system is traditionally accomplished by removing the Reusable Surface Insulation tiles and performing a visual inspection of the aluminum substrate and corrosion protection system. This process is time consuming and has the potential to damage high cost tiles. To evaluate non-intrusive NDE methods, a Proof of Concept (PoC) experiment was designed and test panels were manufactured. The objective of the test plan was three-fold: establish the ability to detect corrosion hidden from view by tiles; determine the key factor affecting detectability; roughly quantify the detection threshold. The plan consisted of artificially inducing dimensionally controlled corrosion spots in two panels and rebonding tile over the spots to model the thermal protective system of the orbiter. The corrosion spot diameter ranged from 0.100" to 0.600" inches and the depth ranged from 0.003" to 0.020". One panel consisted of a complete factorial array of corrosion spots with and without tile coverage. The second panel consisted of randomized factorial points replicated and hidden by tile. Conventional methods such as ultrasonics, infrared, eddy current and microwave methods have shortcomings. Ultrasonics and IR cannot sufficiently penetrate the tiles, while eddy current and microwaves have inadequate resolution. As such, the panels were interrogated using Backscatter Radiography and Terahertz Imaging. The terahertz system successfully detected artificially induced corrosion spots under orbiter tile and functional testing is in-work in preparation for implementation.

Kemmerer, Catherine C.↗

GIFTS SM EDU Level 1B Algorithms

The Geosynchronous Imaging Fourier Transform Spectrometer (GIFTS) SensorModule (SM) Engineering Demonstration Unit (EDU) is a high resolution spectral imager designed to measure infrared (IR) radiances using a Fourier transform spectrometer (FTS). The GIFTS instrument employs three focal plane arrays (FPAs), which gather measurements across the long-wave IR (LWIR), short/mid-wave IR (SMWIR), and visible spectral bands. The raw interferogram measurements are radiometrically and spectrally calibrated to produce radiance spectra, which are further processed to obtain atmospheric profiles via retrieval algorithms. This paper describes the GIFTS SM EDU Level 1B algorithms involved in the calibration. The GIFTS Level 1B calibration procedures can be subdivided into four blocks. In the first block, the measured raw interferograms are first corrected for the detector nonlinearity distortion, followed by the complex filtering and decimation procedure. In the second block, a phase correction algorithm is applied to the filtered and decimated complex interferograms. The resulting imaginary part of the spectrum contains only the noise component of the uncorrected spectrum. Additional random noise reduction can be accomplished by applying a spectral smoothing routine to the phase-corrected spectrum. The phase correction and spectral smoothing operations are performed on a set of interferogram scans for both ambient and hot blackbody references. To continue with the calibration, we compute the spectral responsivity based on the previous results, from which, the calibrated ambient blackbody (ABB), hot blackbody (HBB), and scene spectra can be obtained. We now can estimate the noise equivalent spectral radiance (NESR) from the calibrated ABB and HBB spectra. The correction schemes that compensate for the fore-optics offsets and off-axis effects are also implemented. In the third block, we developed an efficient method of generating pixel performance assessments. In addition, a random pixel selection scheme is designed based on the pixel performance evaluation. Finally, in the fourth block, the single pixel algorithms are applied to the entire FPA.

Tian, Jialin↗

GIFTS SM EDU Data Processing and Algorithms

The Geosynchronous Imaging Fourier Transform Spectrometer (GIFTS) Sensor Module (SM) Engineering Demonstration Unit (EDU) is a high resolution spectral imager designed to measure infrared (IR) radiances using a Fourier transform spectrometer (FTS). The GIFTS instrument employs three Focal Plane Arrays (FPAs), which gather measurements across the long-wave IR (LWIR), short/mid-wave IR (SMWIR), and visible spectral bands. The raw interferogram measurements are radiometrically and spectrally calibrated to produce radiance spectra, which are further processed to obtain atmospheric profiles via retrieval algorithms. This paper describes the processing algorithms involved in the calibration stage. The calibration procedures can be subdivided into three stages. In the pre-calibration stage, a phase correction algorithm is applied to the decimated and filtered complex interferogram. The resulting imaginary part of the spectrum contains only the noise component of the uncorrected spectrum. Additional random noise reduction can be accomplished by applying a spectral smoothing routine to the phase-corrected blackbody reference spectra. In the radiometric calibration stage, we first compute the spectral responsivity based on the previous results, from which, the calibrated ambient blackbody (ABB), hot blackbody (HBB), and scene spectra can be obtained. During the post-processing stage, we estimate the noise equivalent spectral radiance (NESR) from the calibrated ABB and HBB spectra. We then implement a correction scheme that compensates for the effect of fore-optics offsets. Finally, for off-axis pixels, the FPA off-axis effects correction is performed. To estimate the performance of the entire FPA, we developed an efficient method of generating pixel performance assessments. In addition, a random pixel selection scheme is designed based on the pixel performance evaluation.

Tian, Jialin↗

The effect of random Alfven waves on the propagation of hydromagnetic waves in a finite-beta plasma

Using first-order smoothing theory, Fourier analysis and perturbation methods, the evolution equation of the wave spectrum as well as the nonlinear forces generated by random Alfven waves in a finite-beta plasma with phenomenological Landau-damping effects are obtained. The effect of microscale random Alfven waves on the propagation of large-scale hydromagnetic waves is also investigated by solving the mean-field equations. It is shown that parallel-propagating random Alfven waves are modulationally stable and that obliquely propagating random Alfven waves can be modulationally unstable when the energy of random waves is converted to slow magnetoacoustic waves that can be Landau-damped, providing a dissipation mechanism for the Alfven waves.

Hamabata, Hiromitsu↗

Determining the Mass Density Along Magnetic Field Lines from Toroida Eigenfrequencies

Toroidal eigenfrequencies can be used to remotely sense the equatorial mass density rho(sub eq) and the density dependence along a magnetic field line. Here we present improvements to the method of Schulz [1996], which allows rho(sub eq) and the power law index alpha (for mass density along a field line proportional to R(sup -alpha), where R is the radial distance from the center of the Earth) to be determined from the y intercept and slope of a plot of toroidal frequency versus toroidal harmonic number n. Our modifications include a model form for eigenfrequencies with a fractional precision of 0.0005 for -6 less than or = alpha less than or = 6 and 2 less than or = L less than or = 8 (accuracy is doubtful beyond L = 5) and an iterative procedure for getting more accurate results than those found using Schulz's method. In addition, we do an analysis of the effect of random measurement errors. Observed frequencies need to be accurate to approx. 6% (3%) of the fundamental frequency in order to determine rho(sub eq) (alpha) to a precision of 30% (unity). We then apply our method to data generated using the Global Core Plasma Model for plasmaspheric mass density; our analysis demonstrates clearly bow the alpha index represents the mass density dependence on the outer part of the field line (R/(LR(sub E)) greater than or approx. 2/3).

Denton, R.↗

Evaluation of Time-Series Gap-Filling Methods for Solar Irradiance Applications

A complete solar resource data set is essential for any stage of a solar energy project - from feasibility studies to daily operations. But measured or modeled solar resource data are prone to data gaps and data quality issues. To mitigate these issues, a data imputation process should be implemented to obtain a complete and reliable temporal and spatial data series. This study focused on imputing temporal scales by applying random and artificial data gaps and then implementing eight imputation methods, including the Kalman filtering and smoothing and stine interpolations. These methods were implemented on 1-minute to half hourly irradiance data for 1 year using a few locations from the National Solar Radiation Database (NSRDB) and ground measurement data set. The results demonstrated that some of the simpler methods, such as the stine and linear interpolation methods, were the relatively best models based on the statistical metrics for imputing NSRDB and ground measurement data, respectively.

14 SOLAR ENERGY↗

Improved statistical models for limited datasets in uncertainty quantification using stochastic collocation

This paper presents a data-driven framework for performing uncertainty quantification (UQ) by choosing a stochastic model that accurately describes the sources of uncertainty in a system. This model is propagated through an appropriate response surface function that approximates the behavior of this system using stochastic collocation. Given a sample of data describing the uncertainty in the inputs, our goal is to estimate a probability density function (PDF) using the kernel moment matching (KMM) method so that this PDF can be used to accurately reproduce statistics like mean and variance of the response surface function. Instead of constraining the PDF to be optimal for a particular response function, we show that we can use the properties of stochastic collocation to make the estimated PDF optimal for a wide variety of response functions. We contrast this method with other traditional procedures that rely on the Maximum Likelihood approach, like kernel density estimation (KDE) and its adaptive modification (AKDE). We argue that this modified KMM method tries to preserve what is known from the given data and is the better approach when the available data is limited in quantity. We test the performance of these methods for both univariate and multivariate density estimation by sampling random datasets from known PDFs and then measuring the accuracy of the estimated PDFs, using the known PDF as a reference. Comparing the output mean and variance estimated with the empirical moments using the raw data sample as well as the actual moments using the known PDF, we show that the KMM method performs better than KDE and AKDE in predicting these moments with greater accuracy. This improvement in accuracy is also demonstrated for the case of UQ in electrostatic and electrothermomechanical microactuators. We show how our framework results in the accurate computation of statistics in micromechanical systems.

97 MATHEMATICS AND COMPUTING↗

Packetized energy management control systems and methods of using the same

Aspects of the present disclosure include anonymous, asynchronous, and randomized control schemes for distributed energy resources (DERs). Such control schemes may include packetized energy management (PEM) control schemes for managing DERs that may provide near-optimal tracking performance under imperfect information and consumer quality of service (QoS) constraints.

Frolik, Jeffrey↗

Estimating rates and patterns of diversification with incomplete sampling: a case study in the rosids

Premise Recent advances in generating large‐scale phylogenies enable broad‐scale estimation of species diversification. These now common approaches typically are characterized by (1) incomplete species coverage without explicit sampling methodologies and/or (2) sparse backbone representation, and usually rely on presumed phylogenetic placements to account for species without molecular data. We used empirical examples to examine the effects of incomplete sampling on diversification estimation and provide constructive suggestions to ecologists and evolutionary biologists based on those results. Methods We used a supermatrix for rosids and one well‐sampled subclade (Cucurbitaceae) as empirical case studies. We compared results using these large phylogenies with those based on a previously inferred, smaller supermatrix and on a synthetic tree resource with complete taxonomic coverage. Finally, we simulated random and representative taxon sampling and explored the impact of sampling on three commonly used methods, both parametric (RPANDA and BAMM) and semiparametric (DR). Results We found that the impact of sampling on diversification estimates was idiosyncratic and often strong. Compared to full empirical sampling, representative and random sampling schemes either depressed or inflated speciation rates, depending on methods and sampling schemes. No method was entirely robust to poor sampling, but BAMM was least sensitive to moderate levels of missing taxa. Conclusions We suggest caution against uncritical modeling of missing taxa using taxonomic data for poorly sampled trees and in the use of summary backbone trees and other data sets with high representative bias, and we stress the importance of explicit sampling methodologies in macroevolutionary studies.

59 BASIC BIOLOGICAL SCIENCES↗

Switching characteristics for ferroelectric random access memory based on RC model in poly(vinylidene fluoride-trifluoroethylene) ultrathin films

The switching characteristic of the poly(vinylidene fluoride-trifluoroethlene) (P(VDF-TrFE)) films have been studied at different ranges of applied electric field. It is suggest that the increase of the switching speed upon nucleation protocol and the deceleration of switching could be related to the presence of a non-ferroelectric layer. Remarkably, a capacitor and resistor (RC) links model plays significant roles in the polarization switching dynamics of the thin films. For P(VDF-TrFE) ultrathin films with electroactive interlayer, it is found that the switching dynamic characteristics are strongly affected by the contributions of resistor and non-ferroelectric (non-FE) interface factors. A corresponding experiment is designed using poly(3,4-ethylene dioxythiophene):poly(styrene sulfonic) (PEDOT-PSSH) as interlayer with different proton concentrations, and the testing results show that the robust switching is determined by the proton concentration in interlayer and lower leakage current in circuit to reliable applications of such polymer films. These findings provide a new feasible method to enhance the polarization switching for the ferroelectric random access memory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Validation of Force Limited Vibration Testing at NASA Langley Research Center

Vibration tests were performed to develop and validate the forced limited vibration testing capability at the NASA Langley Research Center. The force limited vibration test technique has been utilized at the Jet Propulsion Laboratory and other NASA centers to provide more realistic vibration test environments for aerospace flight hardware. In standard random vibration tests, the payload is mounted to a rigid fixture and the interface acceleration is controlled to a specified level based on a conservative estimate of the expected flight environment. In force limited vibration tests, both the acceleration and force are controlled at the mounting interface to compensate for differences between the flexible flight mounting and rigid test fixture. This minimizes the over test at the payload natural frequencies and results in more realistic forces being transmitted at the mounting interface. Force and acceleration response data was provided by NASA Goddard Space Flight Center for a test article that was flown in 1998 on a Black Brant sounding rocket. The measured flight interface acceleration data was used as the reference acceleration spectrum. Using this acceleration spectrum, three analytical methods were used to estimate the force limits. Standard random and force limited vibration tests were performed and the results are compared with the flight data.

Rice, Chad↗

Milky Way mass with K giants and BHB stars using LAMOST, SDSS/SEGUE, and Gaia : 3D spherical Jeans equation and tracer mass estimator

ABSTRACT We measure the enclosed Milky Way mass profile to Galactocentric distances of ∼70 and ∼50 kpc using the smooth, diffuse stellar halo samples of Bird et al. The samples are Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) and Sloan Digital Sky Survey/Sloan Extension for Galactic Understanding and Exploration (SDSS/SEGUE) K giants (KG) and SDSS/SEGUE blue horizontal branch (BHB) stars with accurate metallicities. The 3D kinematics are available through LAMOST and SDSS/SEGUE distances and radial velocities and Gaia DR2 proper motions. Two methods are used to estimate the enclosed mass: 3D spherical Jeans equation and Evans et al. tracer mass estimator (TME). We remove substructure via the Xue et al. method based on integrals of motion. We evaluate the uncertainties on our estimates due to random sampling noise, systematic distance errors, the adopted density profile, and non-virialization and non-spherical effects of the halo. The tracer density profile remains a limiting systematic in our mass estimates, although within these limits we find reasonable agreement across the different samples and the methods applied. Out to ∼70 and ∼50 kpc, the Jeans method yields total enclosed masses of 4.3 ± 0.95 (random) ±0.6 (systematic) × 1011 M⊙ and 4.1 ± 1.2 (random) ±0.6 (systematic) × 1011 M⊙ for the KG and BHB stars, respectively. For the KG and BHB samples, we find a dark matter virial mass of $M_{200}=0.55^{+0.15}_{-0.11}$ (random) ±0.083 (systematic) × 1012 M⊙ and $M_{200}=1.00^{+0.67}_{-0.33}$ (random) ±0.15 (systematic) × 1012 M⊙, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-E-J-37: Combining Proton Radiography and X-Ray CT Information to Better Estimate Relative Proton Stopping Power in a Clinical Environment

Purpose: In standard proton therapy clinical practice, proton stopping power uncertainties are in the order of 3.5%, which affects the ability of placing the proton Bragg peak at the edge of the tumor. The innovating idea of this project is to approach the uncertainty problem in RSP by using combined information from X-ray CT and proton radiography along a few beam angles. In addition, this project aims to quantify the systematic error introduced by the theoretical models (Janni, ICRU49, Bischel) for proton stopping power in media. Methods: A 3D phantom of 36 cm3 composed of 9 materials randomly placed is created. Measured RSP values are obtained using a Gammex phantom with a proton beam. Theoretical RSP values are calculated with Beth-Block equation in combination with three databases (Janni, ICRU49 and Bischel). Clinical RSP errors are simulated by introducing a systematic (1.5%, 2.5%, 3.5%) and a random error (+/−0.5%) to the theoretical RSP. A ray-tracing algorithm uses each of these RSP tables to calculate energy loss for proton crossing the phantom through various directions. For each direction, gradient descent (GD) method is done on the clinical RSP table to minimize the residual energy difference between the simulation with clinical RSP and with theoretical RSP. The possibility of a systematic material dependent error is investigated by comparing measured RSP to theoretical RSP as calculated from the three models. Results: Using 10,000 iterations on GD algorithm, RSP differences between theoretical values and clinical RSP have converged (<1%) for each error introduced. Results produced with ICRU49 have the smallest average difference (0.021%) to the measured RSP. Janni (1.168%) and Bischel (−0.372%) database shows larger systematic errors. Conclusion: Based on these results, ray-tracing optimisation using information from proton radiography and X-ray CT demonstrates a potential to improve the proton range accuracy in a clinical environment.

07 ISOTOPE AND RADIATION SOURCES↗

SATYA: Defending Against Adversarial Attacks Using Statistical HypothesisTesting

The paper presents a new defense against adversarial attacks for deep neural networks. We demonstrate the effectiveness of our approach against the popular adversarial image generation method DeepFool. Our approach uses Wald’s Sequential Probability Ratio Test to sufficiently sample a carefully chosen neighborhood around an input image to determine the correct label of the image. On a benchmark of 50,000 randomly chosen adversarial images generated by DeepFool we demonstrate that our method SATYA is able to recover the correct labels for 95.76% of the images for CaffeNet and 97.43% of the correct label for GoogLeNet.

Raj, Sunny↗

Analysis of the Uncertainty in Wind Measurements from the Atmospheric Radiation Measurement Doppler Lidar during XPIA: Field Campaign Report

In March and April of 2015, the ARM Doppler lidar that was formerly operated at the Tropical Western Pacific site in Darwin, Australia (S/N 0710-08) was deployed to the Boulder Atmospheric Observatory (BAO) for the eXperimental Planetary boundary-layer Instrument Assessment (XPIA) field campaign. The goal of the XPIA field campaign was to investigate methods of using multiple Doppler lidars to obtain high-resolution three-dimensional measurements of winds and turbulence in the atmospheric boundary layer, and to characterize the uncertainties in these measurements. The ARM Doppler lidar was one of many Doppler lidar systems that participated in this study. During XPIA the 300-m tower at the BAO site was instrumented with well-calibrated sonic anemometers at six levels. These sonic anemometers provided highly accurate reference measurements against which the lidars could be compared. Thus, the deployment of the ARM Doppler lidar during XPIA offered a rare opportunity for the ARM program to characterize the uncertainties in their lidar wind measurements. Results of the lidar-tower comparison indicate that the lidar wind speed measurements are essentially unbiased (~1cm s-1), with a random error of approximately 50 cm s-1. Two methods of uncertainty estimation were tested. The first method was found to produce uncertainties that were too low. The second method produced estimates that were more accurate and better indicators of data quality. As of December 2015, the first method is being used by the ARM Doppler lidar wind value-added product (VAP). One outcome of this work will be to update this VAP to use the second method for uncertainty estimation.

Atmospheric boundary layer, lidar-derived wind spe↗

A Gain Scheduling Wide-Area Damping Controller for the Efficient Integration of Photovoltaic Plant

With photovoltaic generation increasing, the wide-area oscillation has a vital influence on the long-distance power transmission of photovoltaic power. For an efficient integration of photovoltaic generation, a gain scheduling wide-area damping controller (WADC) is designed to guarantee system damping performance under a wide-area signal delay and a changing operating condition, which could be aggravated by the fluctuant and uncertain photovoltaic generation. A linear parameter-varying (LPV) model is proposed to represent all possible operating conditions within varying ranges of the photovoltaic generation and the load. Then, a tensor product (TP) model transformation is applied to transform the LPV model into a TP model. Based on the TP model, a parallel distributed compensation controller (PDC) framework is adopted to design the gain scheduling WADC, which is solved by linear matrix inequalities (LMIs) considering the accommodation to the random interval delay. In addition, a staircase basis function method is proposed to reduce the conservatism of LMIs within the PDC framework. The gain of the proposed WADC will be adjusted according to the operating condition acknowledged through the real-time scheduling and control center, once a wide-area oscillation is detected. Lastly, the performance of the proposed WADC is verified by a 2-area 4-machine system and a 68-bus 16-machine system.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

The {sup 150}Nd({sup 3}He,t) and {sup 150}Sm(t,{sup 3}He) reactions with applications to {beta}{beta} decay of {sup 150}Nd

The {sup 150}Nd({sup 3}He,t) reaction at 140 MeV/u and {sup 150}Sm(t,{sup 3}He) reaction at 115 MeV/u were measured, populating excited states in {sup 150}Pm. The transitions studied populate intermediate states of importance for the (neutrinoless) {beta}{beta} decay of {sup 150}Nd to {sup 150}Sm. Monopole and dipole contributions to the measured excitation-energy spectra were extracted by using multipole decomposition analyses. The experimental results were compared with theoretical calculations obtained within the framework of the quasiparticle random-phase approximation, which is one of the main methods employed for estimating the half-life of the neutrinoless {beta}{beta} decay (0{nu}{beta}{beta}) of {sup 150}Nd. The present results thus provide useful information on the neutrino responses for evaluating the 0{nu}{beta}{beta} and 2{nu}{beta}{beta} matrix elements. The 2{nu}{beta}{beta} matrix element calculated from the Gamow-Teller transitions through the lowest 1{sup +} state in the intermediate nucleus is maximally about half that deduced from the half-life measured in 2{nu}{beta}{beta} direct counting experiments, and at least several transitions through 1{sup +} intermediate states in {sup 150}Pm are required to explain the 2{nu}{beta}{beta} half-life. Because Gamow-Teller transitions in the {sup 150}Sm(t,{sup 3}He) experiment are strongly Pauli blocked, the extraction of Gamow-Teller strengths was complicated by the excitation of the 2({h_bar}/2{pi}){omega}, {Delta}L=0, {Delta}S=1 isovector spin-flip giant monopole resonance (IVSGMR). However, the near absence of Gamow-Teller transition strength made it possible to cleanly identify this resonance, and the strength observed is consistent with the full exhaustion of the non-energy-weighted sum rule for the IVSGMR.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗