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,567 records · Page 87

Kinetic Plasma Simulation Capabilities in the MOOSE Framework: Verification of Particle-Particle Collisions

High-fidelity simulations of complex plasma systems allow researchers to gain key insights into and understanding of these systems. To facilitate massively parallel high-fidelity plasma simulations, finite-element-based particle-in-cell capabilities are being developed within the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) based framework called Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering & Research (SALAMANDER). While SALAMANDER’s primary objective is modeling edge plasmas and plasma-facing components in fusion devices, the particle-in-cell capabilities being developed are general and will support modeling low-temperature plasmas as well. Previously, collisionless magnetostatic simulation capabilities have been verified with the two-stream and Dorey-Guest-Harris instabilities, and single particle motion. Collisions were implemented using the direct simulation Monte Carlo method, and verification of this capability will be presented here several verification problems: relaxation of a randomly initialized gas to a Maxwellian distribution, Fourier heat flow, and comparison of reaction rates to both analytic calculations and those calculated using a multi-term Boltzmann solver.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Overview of Methods for Deriving the Radiative Transfer Theory from the Maxwell Equations. I: Approach Based on the Far-Field Foldy Equations

In this paper, we revisit, with further enhancements and clarifications, the self-consistent first-principles approach developed previously for deriving the vector radiative transfer theory for a discrete random medium with a sparse concentration of particles. We specifically consider the case of a plane-parallel particulate layer embedded in an otherwise homogeneous unbounded medium. The solution method is based on the far-field Foldy equations, an order-of-scattering expansion for the total field derived under the Twersky approximation, the computation of the coherent field by assuming that the positions of the particles are uncorrelated, and the ladder approximation for the coherency dyadic. The latter yields an integral equation for the diffuse specific coherency dyadic, defined through an angular spectrum repre- sentation for the coherency dyadic, which in turn, gives the vector radiative transfer equation for the diffuse specific intensity column vector. We analyze specifically the computation of the coherent field for inhomogeneous particulate media and multiple species of particles, the continuous extension of the far-field representation to the near field, the Foldy approximation, and the Foldy integral equation for the coherent field. Finally, we discuss the transition from the vector to the scalar radiative transfer equation.

Doicu, Adrian↗

Delineating Facies Spatial Distribution by Integrating Ensemble Data Assimilation and Indicator Geostatistics with Level Set Transformation.

A new approach is developed to delineate the spatial distribution of discrete facies (geological units that have unique distributions of hydraulic, physical, and/or chemical properties) conditioned not only on direct data (measurements directly related to facies properties, e.g., grain size distribution obtained from borehole samples) but also on indirect data (observations indirectly related to facies distribution, e.g., hydraulic head and tracer concentration). Our method integrates for the first time ensemble data assimilation with traditional transition probability-based geostatistics. The concept of level set is introduced to build shape parameterization that allows transformation between discrete facies indicators and continuous random variables. The spatial structure of different facies is simulated by indicator models using conditioning points selected adaptively during the iterative process of data assimilation. To evaluate the new method, a two-dimensional semi-synthetic example is designed to estimate the spatial distribution and permeability of two distinct facies from transient head data induced by pumping tests. The example demonstrates that our new method adequately captures the spatial pattern of facies distribution by imposing spatial continuity through conditioning points. The new method also reproduces the overall response in hydraulic head field with better accuracy compared to data assimilation with no constraints on spatial continuity on facies.

97 MATHEMATICS AND COMPUTING↗

Growth and magnetoelectric properties of (00l)-oriented La{sub 0.67}Sr{sub 0.33}MnO{sub 3}/PbZr{sub 0.52}Ti{sub 0.48}O{sub 3} heterostructure films

C-axis oriented La{sub 0.67}Sr{sub 0.33}MnO{sub 3}(LSMO)/PbZr{sub 0.52}Ti{sub 0.48}O{sub 3}(PZT) films are fabricated successfully by sol-gel method on LaAlO{sub 3} (00l) substrates. The structure, composition and morphology of the films are investigated by X-ray diffractometer (XRD, θ-2θ scan, ω-scan and ϕ-scan), X-ray photoelectron spectroscope (XPS), field emission scanning electron microscope (FESEM) and high resolution transmission electron microscope (HRTEM). The electric and magnetic properties of randomly and c-axis oriented LSMO/PZT films are studied comparably using ferroelectric testing apparatus and physical property measurement system (PPMS). It is found that the epitaxial LSMO/PZT composite films show well controlled growth along c-axis, and much better magnetoelectric properties than the randomly oriented ones. The ME voltage coefficient increases from 23 mV cm{sup −1} Oe{sup −1} for the randomly oriented LSMO/PZT composite films to 52 mV cm{sup −1} Oe{sup −1} for c-axis oriented ones prepared using the low cost sol-gel method presented in this study, which shows high potential in promising applications. - Highlights: •Epitaxial LSMO/PZT films were fabricated successfully by sol-gel method on LAO (00l) substrate. •The prepared films exhibit well-defined multiferroic properties for the epitaxial LSMO/PZT films. •Epitaxial LSMO/PZT films show superior magnetoelectric properties to the randomly oriented ones.

36 MATERIALS SCIENCE↗

Optical devices and methods employing nanoparticles, microcavities, and semicontinuous metal films

An optical sensing enhancing material (and corresponding method of making) comprising: a medium, the medium comprising a plurality of aggregated nanoparticles comprising fractals; and a microcavity, wherein the medium is located in a vicinity of the microcavity. Also an optical sensor and sensing method comprising: providing a doped medium, the medium comprising a plurality of aggregated nanoparticles comprising fractals, with the material; locating the doped medium in the vicinity of a microcavity; exciting the doped medium with a light source; and detecting light reflected from the doped medium. Also an optical sensing enhancing material comprising a medium, the medium comprising a semicontinuous metal film of randomly distributed metal particles and their clusters at approximately their percolation threshold. The medium preferably additionally comprises a microcavity/microresonator. Also devices and methods employing such material.

Armstrong, Robert L.↗

Large gradients via correlation in random parameterized quantum circuits

Scaling of variational quantum algorithms to large problem sizes requires efficient optimization of random parameterized quantum circuits. For such circuits with uncorrelated parameters, the presence of exponentially vanishing gradients in cost function landscapes is an obstacle to optimization by gradient descent methods. In this work, we prove that reducing the dimensionality of the parameter space by utilizing circuit modules containing spatially or temporally correlated gate layers can allow one to circumvent the vanishing gradient phenomenon. Here, examples are drawn from random separable circuits and asymptotically optimal variational versions of Grover's algorithm based on the quantum alternating operator ansatz. In the latter scenario, our bounds on cost function variation imply a transition between vanishing gradients and efficient trainability as the number of layers is increased toward $\mathcal{O}\left({2}^{n/2}\right)$, the optimal oracle complexity of quantum unstructured search.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Validation of a Deterministic Vibroacoustic Response Prediction Model

This report documents the recently completed effort involving validation of a deterministic theory for the random vibration problem of predicting the response of launch pad structures in the low-frequency range (0 to 50 hertz). Use of the Statistical Energy Analysis (SEA) methods is not suitable in this range. Measurements of launch-induced acoustic loads and subsequent structural response were made on a cantilever beam structure placed in close proximity (200 feet) to the launch pad. Innovative ways of characterizing random, nonstationary, non-Gaussian acoustics are used for the development of a structure's excitation model. Extremely good correlation was obtained between analytically computed responses and those measured on the cantilever beam. Additional tests are recommended to bound the problem to account for variations in launch trajectory and inclination.

Caimi, Raoul E.↗

Sequential Image Recovery Using Joint Hierarchical Bayesian Learning

Abstract Recovering temporal image sequences (videos) based on indirect, noisy, or incomplete data is an essential yet challenging task. We specifically consider the case where each data set is missing vital information, which prevents the accurate recovery of the individual images. Although some recent (variational) methods have demonstrated high-resolution image recovery based on jointly recovering sequential images, there remain robustness issues due to parameter tuning and restrictions on the type of sequential images. Here, we present a method based on hierarchical Bayesian learning for the joint recovery of sequential images that incorporates prior intra- and inter-image information. Our method restores the missing information in each image by “borrowing” it from the other images. More precisely, we couple sequential images by penalizing their pixel-wise difference. The corresponding penalty terms (one for each pixel and pair of subsequent images) are treated as weakly-informative random variables that favor small pixel-wise differences but allow occasional outliers. As a result, all of the individual reconstructions yield improved accuracy. Our method can be used for various data acquisitions and allows for uncertainty quantification. Some preliminary results indicate its potential use for sequential deblurring and magnetic resonance imaging.

Xiao, Yao↗

Hybrid Ultrasonic TOFD Imaging of Weld Flaws Using Wavelet Transforms and Image Registration

Ultrasonic time of flight diffraction (TOFD) is an effective weld crack inspection technique. Due to the intensity of diffraction wave is rather weak compared with the lateral wave and the bottom echo wave, thus the signal-to-ratio (SNR) of TOFD image is low. A new dichotomous method is comprised of two steps that contains wavelet shrinkage and image registration is proposed in this paper to reduce the noise and improve the resolution of TOFD images as well. In order to evaluate the reliability of our proposed method in this paper, we have established the experiment system, and sampled a number of TOFD data with random distribution of noise characteristics. We adopted one-dimension wavelet transform and two-dimension wavelet transform in the very beginning of the first step of the proposed algorithm respectively. The SNR of the result obtained in this step is improved significantly compared with the classic algorithms. Next, the image registration is applied. After the registered images have been added to form a new one, then it comes to the final result that shows not only the SNR but also the definition of the image is enhanced effectively.

42 ENGINEERING↗

A data-driven framework for sparsity-enhanced surrogates with arbitrary mutually dependent randomness

The challenge of quantifying uncertainty propagation in real-world systems is rooted in the high-dimensionality of the stochastic input and the frequent lack of explicit knowledge of its probability distribution. Traditional approaches show limitations for such problems, especially when the size of the training data is limited. To address these diculties, we have developed a general framework of constructing surrogate models on spaces of stochastic input with arbitrary probability measure irrespective of the mutual dependencies between individual components of the random inputs and the analytical form. The present Data-driven Sparsity-enhancing Rotation for Arbitrary Randomness (DSRAR) framework includes a data-driven construction of multivariate polynomial basis for arbitrary mutually dependent probability measure and a sparsity enhancement rotation procedure. This sparsity-enhancing rotation method was initially proposed in our previous work [1] for Gaussian density distributions, which may not be feasible for non-Gaussian distributions due to the loss of orthogonality after the rotation. To remedy such diculties, we developed a new data-driven approach to construct orthonormal polynomials for polynomials for arbitrary mutually dependent (amdP) randomness, ensuring the constructed basis maintains the orthogonality/near-orthogonality with respect to the density of the rotated random vector, where directly applying the regular polynomial chaos including arbitrary polynomial chaos (aPC) [2] shows limitations due to the assumption of the mutual independence between the components of the random inputs. The developed DSRAR framework leads to accurate recovery, with only limited training data, of a sparse representation of the target functions. The eectiveness of our method is demonstrated in challenging problems such as PDEs and realistic molecular systems within high-dimensional conformational space (O(10)) where the underlying density is implicitly represented by a large collection of sample data, as well as systems with explicitly given non-Gaussian probabilistic measures.

Lei, Huan↗

Layered Electrode Architectures Enabled by Holey Graphene Hosts

Battery electrodes are typically prepared using a slurry-based method, where high boiling point organic solvents are typically used with a polymer binder to disperse active electrode materials and conductive carbon additives followed by casting into current collector sheets. In this process, the use of organic solvents are environmentally hazardous, energy intensive, and time-consuming, while polymeric binders add weight and can undergo parasitic reactions during battery cycling. Here we discuss the use of a unique conductive carbon nanomaterial, namely holey graphene, as a multifunctional host that can be processed via a facile, solvent-free, binder-free, cold pressing method to prepare electrodes with unique architectures. This approach can accommodate a broad range of battery chemistries, such as lithium oxygen/carbon dioxide (Li-O2/CO2), lithium sulfur/selenium (Li-S/Se), and various lithium ion chemistries. The successful use of dry, cold compression allows the convenient modification of electrode mass loading with unique architectures containing active materials randomly mixed or layered with single or multiple layers of one or more active materials. In the conventional slurry-based electrode fabrication method, it is difficult, if not impossible, to achieve such layered architectures. In particular, the layered electrode architectures from the dry, cold press approach are easy to fabricate without significantly sacrificing device performance, and are amenable to continuous roll-to-roll dry processing. In addition, these architectures offer unprecedented perspectives on the battery material loading vs. utilization efficiency, as well as enabling unique characterization opportunities to advance mechanistic understanding of various battery chemistries. Potential applications of this unique electrode fabrication approach for solid-state batteries will also be discussed.

holey graphene, batteries, electrode architecture,↗

Analysis and prediction of intersection traffic violations using automated enforcement system data

We report that the automated enforcement system (AES) is an effective way of supplementing traditional traffic enforcement, and the traffic violation data from AES can also be effectively used for safety research. In this study, traffic violation data were used to analyze the influencing factors associated with traffic violations and to predict the probability of violations at intersections. The potential factors influencing violations include 24 independent factors related to time, space, traffic and weather. Results from a logistic model showed that the midday period, weekends, residential districts, collector roads, congested traffic conditions, high traffic flow, lower wind speed and low temperature would increase the probability of traffic violations. The probability of violations was predicted by the random forest algorithm, which was proven to be the best traffic violation prediction model among logistic regression, Gaussian naive Bayes, and support vector machine. Moreover, the proximity weighted synthetic oversampling technique (ProWSyn) method was applied to reduce the impact of the imbalance ratio (IR) and improve the model’s prediction performance. The receiver operating characteristics (ROC) curves and Precision-Recall (PR) curves illustrated that the random forest algorithm using oversampling data had the best classifier prediction performance than undersampling data. The area under curve (AUC) and out-of-bag (OOB) error with IR = 1 reached 0.914 and 0.0787, which showed the better performance of the random forest algorithm using ProWSyn in dealing with imbalanced traffic violation data.

42 ENGINEERING↗

Generation of Tunable Stochastic Sequences Using the Insulator–Metal Transition

Probabilistic computing is a paradigm in which data are not represented by stable bits, but rather by the probability of a metastable bit to be in a particular state. The development of this technology has been hindered by the availability of hardware capable of generating stochastic and tunable sequences of “1s” and “0s”. The options are currently limited to complex CMOS circuitry and, recently, magnetic tunnel junctions. Here, we demonstrate that metal–insulator transitions can also be used for this purpose. We use an electrical pump/probe protocol and take advantage of the stochastic relaxation dynamics in VO 2 to induce random metallization events. A simple latch circuit converts the metallization sequence into a random stream of 1s and 0s. The resetting pulse in between probes decorrelates successive events, providing a true stochastic digital sequence.

97 MATHEMATICS AND COMPUTING↗

Effects of Dose Error and Sample Size on Sonic Boom Dose-Response Curves

NASA will soon be collecting noise-annoyance community survey data as the X-59 aircraft flies supersonically over several communities in the USA. Sparse measurements of the X-59 sonic thumps will be used together with physics-based simulations to estimate noise doses at survey participant locations. These dose estimates have associated error that affects the accuracy of modeled dose-response curves, which can result in misestimation of annoyance. The precision in dose-response curves is also a consideration in selecting the number of survey participants. To enable pretest studies of dose error and precision, simulated dose-response data were generated based on NASA’s Quiet Supersonic Flights 2018 test. The data included various degrees of dose error and sample size. Frequentist multilevel logistic regression models were fit to the true and perturbed dose-response data. Simple proportional relationships were identified between the model parameters and the perturbation standard deviation. The summary dose-response curves illustrate the impact on accuracy if dose error is not accounted for in the model. The precision in the dose-response curves is also shown as the number of participants and degree of participation is varied. Finally, sampling variability is illustrated by showing the dose-response curves for several replicates with random draws of participants and errors.

X-59↗

Inhibition of Notch signaling pathway using γ-secretase inhibitor delivered by a low dose of Triton-X100 in cultured oral cancer cells

Highlights: • A novel method is developed to deliver γ-secretase inhibitor DAPT using a low dose of Triton-X100. • This approach effectively blocks the cleavage of Notch1 receptor. • The elimination of Notch1 intracellular domain from cell nuclei is confirmed. • Notch1 signaling blockade in mRNA level is demonstrated upon applying this delivery system. How to effectively delivering therapeutic agents, including γ-secretase inhibitors (GSIs), into live cells, remains a significant challenge. This study assessed the effect of Notch signaling inhibition by examining levels of the Notch1 intracellular domain (N1ICD) in cultured oral cancer cells analyzed with random stitched images (2D) and 3D visualizations using confocal microscopy and quantitative gene analysis. Substantially, we have developed a novel method to assist the delivery of γ-secretase inhibitor, DAPT, into live cells in the presence of an effective minimum concentration of Triton-X100 (0.001%) without damaging cell activity and membrane integrity assessed with cell proliferation assays. The images obtained in this study showed that DAPT alone could not block the γ-secretase inhibitor despite inhibiting cell growth. Further analysis of quantitative gene expressions of Notch signaling canonical pathway to verify the effectiveness of the novel method for delivering inhibitor into live cells, displayed deregulation of Notch1, Delta-like ligand 1 (DLL1) and hairy and enhancer of split 1 (Hes1). Our data suggest that Notch1/Hes1 signaling pathway is deactivated using DAPT with a low dose of Triton-X100 in this cancer cells. And the finding also suggests that Notch1 could be engaged by DLL1 to promote differentiation in oral cancer cells. Using this approach, we demonstrate that Triton-X100 is a promising and effective permeabilization agent to deliver γ-secretase inhibitor DAPT into live oral epithelial cells. This strategy has the potential to implicate in the treatment of cancer diseases.

60 APPLIED LIFE SCIENCES↗

Faster Randomized Dynamical Decoupling

We present a randomized dynamical decoupling (DD) protocol that can substantially improve the performance of any given deterministic DD scheme for suppressing coherent noise by using no more than two additional pulses. Our construction is implemented by probabilistically applying sequences of pulses, which, when combined, effectively eliminate the error terms that scale linearly with the system-environment coupling strength. As a result, we show that a randomized protocol using a few pulses can outperform deterministic DD protocols that require considerably more pulses. Furthermore, we prove that the randomized protocol provides an improvement compared to deterministic DD sequences that aim to reduce the error in the system’s Hilbert space, such as Uhrig DD, which had been previously regarded to be optimal. To rigorously evaluate the performance, we introduce new analytical methods suitable for analyzing higher-order DD protocols that might be of independent interest. Here, we also present numerical simulations confirming the significant advantage of using randomized protocols compared to widely used deterministic protocols.

Quantum algorithms & computation↗

An uncertainty quantification method relevant to material test reactors

Within material test reactor calculations, energy dependent flux and reaction rate uncertainties are typically not quantified when performing as-run analyses to determine the neutron field experienced by the experiment. When high fidelity Monte-Carlo codes are used in such analyses, straight forward methods to calculate output uncertainties are not available, instead expert opinion is used to postulate computational uncertainties. New methods to propagate uncertainties through these high fidelity simulations are available when sufficient computational power is available. A tool is developed here for sampling any part of an MCNP input from random distributions to determine output uncertainties based on those inputs. Another tool is developed to sample nuclear data cross-section in ACE format using multi-group nuclear data covariances. The Total Monte-Carlo Method and Gesellschaft für Anlagen-und Reaktorsicherheit method (GRS) are implemented and compared to one another as well as MCNP sensitivity and uncertainty calculations. The methods were applied to the Godiva critical sphere k-eigenvalue, the UAM pin-cell benchmark energy dependent flux and reaction rates, and the Advanced Test Reactor energy dependent flux within an experimental location. Furthermore, the two methods agree well, with GRS allowing for an order of magnitude speedup for reaction rate uncertainty calculations and several orders of magnitude for eigenvalue uncertainty calculations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗