Search NASA⌕ Search

SEARCH · Search NASA

Results for “random processes”

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 739 records · Page 41

Surrogate Hessian accelerated structural optimization for stochastic electronic structure theories

In this work, we present an efficient energy-based method for structural optimization with stochastic electronic structure theories, such as diffusion quantum Monte Carlo (DMC). This method is based on robust line-search energy minimization in reduced parameter space, exploiting approximate but accurate Hessian information from a surrogate theory, such as density functional theory. The surrogate theory is also used to characterize the potential energy surface, allowing for simple but reliable ways to maximize statistical efficiency while retaining controllable accuracy. We demonstrate the method by finding the minimum DMC energy structures of the selected flake-like aromatic molecules, such as benzene, coronene, and ovalene, represented by 2, 6, and 19 structural parameters, respectively. In each case, the energy minimum is found within two parallel line-search iterations. The method is near-optimal for a line-search technique and suitable for a broad range of applications. It is easily generalized to any electronic structure method where forces and stresses are still under active development and implementation, such as diffusion Monte Carlo, auxiliary-field Monte Carlo, and stochastic configuration interaction, as well as deterministic approaches such as the random-phase approximation. Accurate and efficient means of geometry optimization could shed light on a broad class of materials and molecules, showing high sensitivity of induced properties to structural variables.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ferroelectric to paraelectric phase transition mechanism in poled PVDF-TrFE copolymer films

Direct experimental insights into the structural and dynamical mechanisms for ferroelectric β to paraelectric α phase transition in a poled PVDF-TrFE copolymer is obtained from in situ x-ray diffraction and quasielastic neutron scattering measurements at high temperatures. It is observed that the β-to-α phase transition proceeds through two energetically distinct processes, which are identified here as the nucleation and growth of an intermediate γ phase with random skew linkages followed by a γ-to-α transition. The two energetically distinct microscopic processes can explain the stages of evolution for β-to-α phase transition observed from heat flow measurements.

36 MATERIALS SCIENCE↗

Data-driven Minimum Entropy Control for Stochastic Nonlinear Systems using the Cumulant-Generating Function

Here, we present a novel minimum entropy control algorithm for a class of stochastic nonlinear systems subjected to non-Gaussian noises. The entropy control can be considered as an optimization problem for the system randomness attenuation, but the mean value has to be considered separately. To overcome this disadvantage, a new representation of the system stochastic properties was given using the cumulant-generating function based on the moment-generating function, in which the mean value and the entropy was reflected by the shape of the cumulant-generating function. Based on the samples of the system output and control input, a time-variant linear model was identified, and the minimum entropy optimization was transformed to system stabilization. Then, an optimal control strategy was developed to achieve the randomness attenuation, and the boundedness of the controlled system output was analyzed. The effectiveness of the presented control algorithm was demonstrated by a numerical example. In this paper, a data-driven minimum entropy design is presented without pre-knowledge of the system model; entropy optimization is achieved by the system stabilization approach in which the stochastic distribution control and minimum entropy are unified using the same identified structure; and a potential framework is obtained since all the existing system stabilization methods can be adopted to achieve the minimum entropy objective.

42 ENGINEERING↗

Meso-Structured Polymer Electrolyte Fuel Cell Electrode

Increasing the utilization of Pt and Pt alloy catalysts in polymer electrolyte fuel cell cathodes is critical to improving the high power density operation, particularly at low Pt loadings. State of the art electrodes are fabricated in an ink deposition process that leads to uncontrolled electrode architecture with random aggregates of functional domains (catalyst, ionomer, and pore volume) (1). The randomness in the domains induces high tortuosity transport pathways for ions and fluids, which cause severe transport resistance during high current density operation. Thin ionomer films cause additional transport resistance and poisoning of the Pt catalyst, which becomes more significant at low Pt loadings. Reducing the amount of ionomer in the catalyst domain without affecting the ionic transport resistance is key to improving the utilization of the Pt and reducing the transport resistance at low Pt loading. Rational design of the electrode structure with controlled low tortuous ionic transport pathways could improve performance. The introduction of the ionomer pathways could also enable reduction of the ionomer volume in the catalyst domain, reducing the transport resistance. Middelmen et al. proposed electrode structures consisting of aligned components in a low tortuosity configuration to improve performance (2). In this work, we present an alternative electrode structure based on a vertically aligned array of Nafion pillars in the cathode catalyst layer, as shown in Figure 1a. Figure 1b shows the SEM image of the Nafion pillars. Furthermore, Pt supported on carbon catalyst was deposited on the Nafion pillars to fabricate a meso-structured electrode. Nafion pillars provide high conductive and low tortuous pathways for protons, reducing the effective transport distance, and enabling reduction of the ionomer binder in the catalyst domain.

25 ENERGY STORAGE↗

Using Machine Learning to Discern Eruption in Noisy Environments: A Case Study Using CO2-Driven Cold-Water Geyser in Chimayó, New Mexico

In this paper, we present an approach based on machine learning (ML) to distinguish eruption and precursory signals of Chimayó geyser (New Mexico, U.S.A.) under noisy environmental conditions. This geyser can be considered a natural analog of CO 2 intrusion into shallow water aquifers. By studying this geyser, we can understand upwelling of CO 2 -rich fluids from depth, which has relevance to leak monitoring in a CO 2 sequestration project. ML methods such as random forests (RFs) are known to be robust multiclass classifiers and perform well under unfavorable, noisy conditions. However, the extent of the RF method’s accuracy is poorly understood for this CO 2 -driven geysering application. The current study aims to quantify the performance of RF classifiers to discern the geyser state. Toward this goal, we first present the data collected from the seismometer that is installed near the Chimayó geyser. The seismic signals collected at this site contain different types of noises such as daily temperature variations, animal movement near the geyser, and human activity. First, we filter the signals from these noises by combining the Butterworth high-pass (BH) filter and an autoregressive (AR) method in a multilevel fashion. We show that by combining these filtering techniques in a hierarchical fashion leads to a reduction in noise in the seismic data without removing the precursors and eruption event signals. We then use RF on the filtered data to classify the state of geyser into three classes: remnant noise, precursor, and eruption states. RF classifier is constructed based on the comprehensive features extracted using the Tsfresh Python package. We show that the classification accuracy using RF on the filtered data is greater than 90%. We also evaluate the accuracy of other classical time-series methods such as dynamic time warping (DTW) on filtered data along with RF on partially filtered data in which we remove the daily trends. Classification accuracy shows that DTW performs poorly (44%) and RF on partially filtered data performs decently (87%). Denoising seismic signals from both daily trends and human activity enhances RF classifier performance by 7%. Finally, these aspects make the proposed ML framework attractive for event discrimination and signal enhancement under noisy conditions, with strong potential for application to monitor leaks in CO 2 sequestration.

58 GEOSCIENCES↗

Model-Form Epistemic Uncertainty Quantification for Modeling with Differential Equations: Application to Epidemiology

Modeling real-world phenomena to any degree of accuracy is a challenge that the scientific research community has navigated since its foundation. Lack of information and limited computational and observational resources necessitate modeling assumptions which, when invalid, lead to model-form error (MFE). The work reported herein explored a novel method to represent model-form uncertainty (MFU) that combines Bayesian statistics with the emerging field of universal differential equations (UDEs). The fundamental principle behind UDEs is simple: use known equational forms that govern a dynamical system when you have them; then incorporate data-driven approaches – in this case neural networks (NNs) – embedded within the governing equations to learn the interacting terms that were underrepresented. Utilizing epidemiology as our motivating exemplar, this report will highlight the challenges of modeling novel infectious diseases while introducing ways to incorporate NN approximations to MFE. Prior to embarking on a Bayesian calibration, we first explored methods to augment the standard (non-Bayesian) UDE training procedure to account for uncertainty and increase robustness of training. In addition, it is often the case that uncertainty in observations is significant; this may be due to randomness or lack of precision in the measurement process. This uncertainty typically manifests as “noisy” observations which deviate from a true underlying signal. To account for such variability, the NN approximation to MFE is endowed with a probabilistic representation and is updated using available observational data in a Bayesian framework. By representing the MFU explicitly and deploying an embedded, data-driven model, this approach enables an agile, expressive, and interpretable method for representing MFU. In this report we will provide evidence that Bayesian UDEs show promise as a novel framework for any science-based, data-driven MFU representation; while emphasizing that significant advances must be made in the calibration of Bayesian NNs to ensure a robust calibration procedure.

97 MATHEMATICS AND COMPUTING↗

Dynamic Control of Cutting in Terms of Elastoviscous Perturbations

The simulation of perturbations in machining is considered. The influence of perturbation on the composition of the oscillatory process in the tool and workpiece subsystems is discussed. The composition of the oscillatory process—that is, the proportions of ordered and random components—may be controlled by adjusting the machining conditions.

42 ENGINEERING↗

The bulk-Moon MgO/FeO ratio: A highlands perspective

Compositional data for nonmare (highlands) samples suggest that the Moon's mg ratio (MgO/FeO) is higher than general estimates. Geochemically representative highlands soils have mg ratios of 0.66 (Apollo 16), 0.69 (Luna 20) and 0.73 (ALHA81005). These soils are mixtures of unrelated pristine nonmare rocks, of which there are at least three groups: Mg-rich rocks, ferroan anorthosites, and KREEP. Other than Mg-rich rocks, virtually all pristine rocks have mg 0.65. Thus, assuming the mixing process that sampled Mg-rich materials was random, the average mg of Mg-rich parent magmas was probably at least 0.70. More direct evidence can be derived from the Mg-rich rocks themselves. Nine of them have bulk-rock mg 0.87. Two (15445 A and 67435 PST) contain Fo(92) olivine. Production of melts that crystallized Fo(92) olivine implies that the mg ratios of source regions in lunar mantle were commensurably high. A quantification of this constraint is developed assuming that the parent melts formed by equilibrium (batch) partial melting. Implications of the model are discussed.

Warren, P. H.↗

Optimum structural sizing for gust induced response

The present paper describes an optimization capability for sizing airframe structures that are subjected to a combination of deterministic and random loads. Design constraints are implemented to prevent structural failure in fatigue and due to a single exceedance of the allowable stress. The random load is treated as a stationary, homogeneous process with a Gaussian probability distribution, and a frequency domain method is used for the computation of dynamic response parameters. An equal-probability-of-load-combination criterion is proposed in the formulation of strength constraints. This alleviates problems associated with the incomplete phase information available when using power spectral density methods. The optimization procedure is illustrated by examples using typical built-up structures.

Hajela, P.↗

Particle image velocimetry for the surface tension driven convection experiment using a particle displacement tracking technique

The Surface Tension Driven Convection Experiment (STDCE) is a Space Transportation System flight experiment to study both transient and steady thermocapillary fluid flows aboard the USML-1 Spacelab mission planned for 1992. One of the components of data collected during the experiment is a video record of the flow field. This qualitative data is then quantified using an all electronic, two-dimensional particle image velocimetry technique called particle displacement tracking (PDT) which uses a simple space domain particle tracking algorithm. The PDT system is successful in producing velocity vector fields from the raw video data. Application of the PDT technique to a sample data set yielded 1606 vectors in 30 seconds of processing time. A bottom viewing optical arrangement is used to image the illuminated plane, which causes keystone distortion in the final recorded image. A coordinate transformation was incorporated into the system software to correct this viewing angle distortion. PDT processing produced 1.8 percent false identifications, due to random particle locations. A highly successful routine for removing the false identifications was also incorporated, reducing the number of false identifications to 0.2 percent.

Wernet, Mark P.↗

Mullite Whiskers and Mullite-whisker Felt

The Naval Surface Warfare Center has developed processes for the preparation of mullite (3(Al2O3)(dot)2(SiO2)) whiskers and mullite-whisker felt. Three patents on the technology were issued in 1990. The processes are based on chemical reactions between AlF3, Al2O3, and SiO2. The felt is formed in-situ during the processing of shaped powdered precursors. It consists of randomly oriented whiskers which are mutually intergrown forming a rigid structure. The microstructure and properties of the felt and size of the whiskers can be modified by varying the amount of Al2O3 in the starting mixture. Loose mullite whiskers can be used as a reinforcement for polymer-, metal-, and ceramic-matrix composites. The felt can be used as preforms for fabricating composite materials as well as for thermal insulation and high temperature, chemically stable filters for liquids (melts) and gases.

Talmy, Inna G.↗

Rutger's CAM2000 chip architecture

This report describes the architecture and instruction set of the Rutgers CAM2000 memory chip. The CAM2000 combines features of Associative Processing (AP), Content Addressable Memory (CAM), and Dynamic Random Access Memory (DRAM) in a single chip package that is not only DRAM compatible but capable of applying simple massively parallel operations to memory. This document reflects the current status of the CAM2000 architecture and is continually updated to reflect the current state of the architecture and instruction set.

Smith, Donald E.↗

Particle image velocimetry for the Surface Tension Driven Convection Experiment using a particle displacement tracking technique

The Surface Tension Driven Convection Experiment (STDCE) is a Space Transportation System flight experiment to study both transient and steady thermocapillary fluid flows aboard the USML-1 Spacelab mission planned for 1992. One of the components of data collected during the experiment is a video record of the flow field. This qualitative data is then quantified using an all electronic, two-dimensional particle image velocimetry technique called particle displacement tracking (PDT) which uses a simple space domain particle tracking algorithm. The PDT system is successful in producing velocity vector fields from the raw video data. Application of the PDT technique to a sample data set yielded 1606 vectors in 30 seconds of processing time. A bottom viewing optical arrangement is used to image the illuminated plane, which causes keystone distortion in the final recorded image. A coordinate transformation was incorporated into the system software to correct this viewing angle distortion. PDT processing produced 1.8 percent false identifications, due to random particle locations. A highly successful routine for removing the false identifications was also incorporated, reducing the number of false identifications to 0.2 percent.

Wernet, Mark P.↗

Idealized textile composites for experimental/analytical correlation

Textile composites are fiber reinforced materials produced by weaving, braiding, knitting, or stitching. These materials offer possible reductions in manufacturing costs compared to conventional laminated composites. Thus, they are attractive candidate materials for aircraft structures. To date, numerous experimental studies have been performed to characterize the mechanical performance of specific textile architectures. Since many materials and architectures are of interest, there is a need for analytical models to predict the mechanical properties of a specific textile composite material. Models of varying sophistication have been proposed based on mechanics of materials, classical laminated plate theory, and the finite element method. These modeling approaches assume an idealized textile architecture and generally consider a single unit cell. Due to randomness of the textile architectures produced using conventional processing techniques, experimental data obtained has been of limited use for verifying the accuracy of these analytical approaches. This research is focused on fabricating woven textile composites with highly aligned and accurately placed fiber tows that closely represent the idealized architectures assumed in analytical models. These idealized textile composites have been fabricated with three types of layer nesting configurations: stacked, diagonal, and split-span. Compression testing results have identified strength variations as a function of nesting. Moire interferometry experiments are being used to determine localized deformations for detailed correlation with model predictions.

Adams, Daniel O.↗

Variance Analysis of Unevenly Spaced Time Series Data

We have investigated the effect of uneven data spacing on the computation of delta (sub chi)(gamma). Evenly spaced simulated data sets were generated for noise processes ranging from white phase modulation (PM) to random walk frequency modulation (FM). Delta(sub chi)(gamma) was then calculated for each noise type. Data were subsequently removed from each simulated data set using typical two-way satellite time and frequency transfer (TWSTFT) data patterns to create two unevenly spaced sets with average intervals of 2.8 and 3.6 days. Delta(sub chi)(gamma) was then calculated for each sparse data set using two different approaches. First the missing data points were replaced by linear interpolation and delta (sub chi)(gamma) calculated from this now full data set. The second approach ignored the fact that the data were unevenly spaced and calculated delta(sub chi)(gamma) as if the data were equally spaced with average spacing of 2.8 or 3.6 days. Both approaches have advantages and disadvantages, and techniques are presented for correcting errors caused by uneven data spacing in typical TWSTFT data sets.

Hackman, Christine↗

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↗

Integrated circuit physically unclonable function

An integrated circuit (IC) based physically unclonable function (PUF) that comprises a common source amplifier for generating PUF output voltages, a unity gain, negative feedback operational amplifier for generating bias voltages, a voltage regulator and a bit exclusion circuit that excludes unstable PUF bits. Compensation circuitry built into the IC-PUF provides a high power supply rejection ratio and enables highly reliable operation of the IC-PUF across varying input voltages and operating temperatures. The IC-PUF generates a uniformly random output bit stream by taking advantage of process variations that are inherent to the fabrication of (metal-oxide semiconductor) MOS transistors.

42 ENGINEERING↗

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Implementation and Performance

Developing theoretical understanding of complex reactions and processes at interfaces requires using methods that go beyond semilocal density functional theory to accurately describe the interactions between solvent, reactants and substrates. Methods based on many-body perturbation theory, such as the random phase approximation (RPA), have previously been limited due to their computational complexity. However, this is now a surmountable barrier due to the advances in computational power available, in particular through modern GPU-based supercomputers. In this work, we describe the implementation of RPA calculations within BerkeleyGW and show its favorable computational performance on large complex systems relevant for catalysis and electrochemistry applications. Our implementation builds off of the static subspace approximation which, by employing a compressed representation of the frequency dependent polarizability, enables the evaluation of the RPA correlation energy with significant acceleration and systematically controllable accuracy. We find that the computational cost of calculating the RPA correlation energy scales only linearly with system size for systems containing up to 50 thousand bands, and is expected to scale quadratically thereafter. We also show excellent strong scaling results across several supercomputers, demonstrating the performance and portability of this implementation.

algorithmic development↗