Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stochastic inference”

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 145 records · Page 8

Dynamical constraints on the formation and evolution of planetary bodies

The present investigation is concerned with a number of inferences as to the origin of planetary bodies, taking into account the present dynamical state of the solar system and some of the limitations which apply to the considered conclusions. Attention is given to the dynamical processes, specifically those processes which may have influenced the orbital or rotational properties of the planets and satellites. Collisional processes are explored, taking into consideration orbital spacing, planetary rotation, and stochastic effects. In connection with a discussion of the evolution of rotational motion, spin state evolution is investigated along with spin axis precession and resonance variation, and the Cassini states. The evolution of planetary orbits is also studied. The subjects considered are related to tides, secular resonances, disk dynamics, and disk-satellite interactions.

Harris, A. W.↗

Shape‐Dependent Multi‐Weight Magnetic Artificial Synapses for Neuromorphic Computing

Abstract In neuromorphic computing, artificial synapses provide a multi‐weight (MW) conductance state that is set based on inputs from neurons, analogous to the brain. Herein, artificial synapses based on magnetic materials that use a magnetic tunnel junction (MTJ) and a magnetic domain wall (DW) are explored. By fabricating lithographic notches in a DW track underneath a single MTJ, 3–5 stable resistance states that can be repeatably controlled electrically using spin‐orbit torque are achieved. The effect of geometry on the synapse behavior is explored, showing that a trapezoidal device has asymmetric weight updates with high controllability, while a rectangular device has higher stochasticity, but with stable resistance levels. The device data is input into neuromorphic computing simulators to show the usefulness of application‐specific synaptic functions. Implementing an artificial neural network (NN) applied to streamed Fashion‐MNIST data, the trapezoidal magnetic synapse can be used as a metaplastic function for efficient online learning. Implementing a convolutional NN for CIFAR‐100 image recognition, the rectangular magnetic synapse achieves near‐ideal inference accuracy, due to the stability of its resistance levels. This work shows MW magnetic synapses are a feasible technology for neuromorphic computing and provides design guidelines for emerging artificial synapse technologies.

Leonard, Thomas↗

Connection between ambient density fluctuations and clumpy Langmuir waves in type III radio sources

A recent stochastic-growth theory of clumpy Langmuir waves in type III sources is shown to imply that the clumps will have the same size distribution as the ambient low-frequency density fluctuations in the solar wind. Spectral analysis of Langmuir-wave time series from the ISEE 3 plasma wave instrument confirms this prediction to within the uncertainties in the spectra. The smallest Langmuir clump size is inferred to be in the range 0.4-30 km in general, and 2-30 km for beam-resonant waves, and it is concluded that the diffusion of waves in the source is anomalous.

Robinson, P. A.↗

Stochastic modeling and statistical calibration with model error and scarce data

This paper introduces a procedure to assess the predictive accuracy of stochastic models subject to model error and sparse data. Model error is introduced as uncertainty on the coefficients of appropriate polynomial chaos expansions (PCE). The error associated with finite sample size allows us to conceive of these coefficients as statistics of the data that we describe as random variables whose influence on output quantities of interest is evaluated through the extended polynomial chaos expansion (EPCE). A Bayesian data assimilation scheme is introduced to update these expansions by considering the resulting nested chaos expansion as a hierarchical probabilistic model. Stochastic models of quantities of interest (QoI) are thus constructed and efficiently evaluated. Here, the Metropolis–Hastings Markov chain Monte Carlo procedure is used to sample the posterior. Two illustrative analytical and numerical problems are used to demonstrate the proposed approach.

Bayesian inference↗

Molecular phylogenetic trees - On the validity of the Goodman-Moore augmentation algorithm

A response is made to the reply of Nei and Tateno (1979) to the letter of Holmquist (1978) supporting the validity of the augmentation algorithm of Moore (1977) in reconstructions of nucleotide substitutions by means of the maximum parsimony principle. It is argued that the overestimation of the augmented numbers of nucleotide substitutions (augmented distances) found by Tateno and Nei (1978) is due to an unrepresentative data sample and that it is only necessary that evolution be stochastically uniform in different regions of the phylogenetic network for the augmentation method to be useful. The importance of the average value of the true distance over all links is explained, and the relative variances of the true and augmented distances are calculated to be almost identical. The effects of topological changes in the phylogenetic tree on the augmented distance and the question of the correctness of ancestral sequences inferred by the method of parsimony are also clarified.

Holmquist, R.↗

Stochastic acceleration of solar flare protons

The acceleration of solar flare protons by the cyclotron damping of intense Alfven wave turbulence in a magnetic trap is considered. The energy diffusion coefficient is computed for a near-isotropic distribution of super-Alfvenic protons, and a steady-state solution for the particle spectrum is found for both transit-time and diffusive losses out of the ends of the trap. The acceleration time to a characteristic energy of about 20 MeV per nucleon can be as short as 10 sec. On the basis of phenomenological arguments, it is inferred that the Alfven wave spectrum typically has an inverse square frequency dependence and that the correlation time of the turbulence lies in the range between 0.5 and 50 msec.

Barbosa, D. D.↗

Experimental and calculational investigation of laser-heated additive manufactured foams

Foam materials are starting to find application in laser-heated Hohlraums used to drive inertial confinement fusion implosions. Foams made using additive manufacturing (AM) techniques are now available and may have advantages over traditional chemical (aerogel) foams. Here, we present new experimental data on laser-heated AM foams. Samples of four different types of printed AM foams were heated using a single 527 nm laser beam at the Jupiter Laser Facility. The laser pulse was ∼180 J square pulse with an FWHM of 1.6 ns and a peak intensity of 3–4 × 1014 W/cm2. The foam densities ranged from 12 to 93 mg/cc (all supercritical for 527 nm light). We measured the backscattered light (power and spectrum), the transmitted light, side-on x-ray images, and the Ti K-shell emission that was used to infer the time-integrated temperature. The fraction of backscattered light was 6%–15% of the input laser energy. The pure carbon foam sample had less backscatter than a C8H9O3 foam of similar density, which was consistent with multi-fluid calculations that predicted less ion heating for the C8H9O3 foam. The level of backscatter and the thermal front speeds for the AM foams were similar to values measured for stochastic (aerogel) foams under similar conditions.

36 MATERIALS SCIENCE↗

Neural network-based model of galaxy power spectrum: fast full-shape galaxy power spectrum analysis

ABSTRACT We present a neural network-based emulator for the galaxy redshift-space power spectrum that enables several orders of magnitude acceleration in the galaxy clustering parameter inference, while preserving 3$\sigma$ accuracy better than 0.5 per cent up to $k_{\mathrm{max}}$ = 0.25 $\, h\text{Mpc}^{-1}$ within Lambda-cold dark matter ($\Lambda$CDM) and around 0.5 per cent $w_0$–$w_a$CDM. Our surrogate model only emulates the galaxy bias-invariant terms of one-loop perturbation theory predictions, these terms are then combined analytically with galaxy bias terms, counter-terms, and stochastic terms in order to obtain the non-linear redshift-space galaxy power spectrum. This allows us to avoid any galaxy bias prescription in the training of the emulator, which makes it more flexible. Moreover, we include the redshift $z \in [0,1.4]$ in the training which further avoids the need for re-training the emulator. We showcase the performance of the emulator in recovering the cosmological parameters of $\Lambda$CDM by analysing the suite of 25 AbacusSummit simulations that mimic the Dark Energy Spectroscopic Instrument luminous red galaxies at $z=0.5$ and 0.8, together as the emission line galaxies at $z=0.8$. We obtain similar performance in all cases, demonstrating the reliability of the emulator for any galaxy sample at any redshift in $0 \lt z \lt 1.4$. We will make our emulator public at github repository.

Trusov, Svyatoslav (ORCID:0000000224146720)↗

Stochastic acceleration of solar flare protons

The acceleration of solar flare protons is considered by cyclotron damping of intense Alfven wave turbulence in a magnetic trap. The energy diffusion coefficient is computed for a near-isotropic distribution of super-Alfvenic protons and a steady-state solution for the particle spectrum is found for both transit-time and diffusive losses out of the ends of the trap. The acceleration time to a characteristic energy approximately 20 Mev/nucl can be as short as 10 sec. On the basis of phenomenological arguments an omega/2 frequency dependence for the Alfven wave spectrum is inferred. The correlation time of the turbulence lies in the range .0005 less than tau/corr less than .05s.

Barbosa, D. D.↗

Competition between roughness and strength for scale-dependent surfaces

Rocks famously have scale-dependent strength, yet the actual dependence is notoriously hard to measure or incorporate into any theoretical framework. Natural rough surfaces present an opportunity to solve the problem. Surfaces sliding in shear evolve as protrusions collide. These asperities can deform or break, thus creating a new surface shape. In particular, natural surfaces have roughness at all scales as well as scale-dependent strength. Based on a scaling analysis, we have previously suggested that the scale-dependent aspect ratio of steady-state surfaces should be proportional to the scale-dependent shear strain at yield. If true, scale-dependent strength could easily be inferred from natural surfaces. Thus, moving beyond the scaling argument to a rigorous treatment of scale-dependent strength for multiscale rough surfaces in shear is important. However, analytic frameworks for analyzing multiscale problems are challenging, as conventional continuum mechanics typically involves a single value for a material property across scales. Here, in this work, we build on the formalism of Persson (2001) that presents a method to compute contact area for rough surfaces with a prescribed topographic spectrum using a stochastic differential equation. The Persson formalism allows for plastic yield under normal loading of otherwise elastic materials and leaves open the possibility of scale-dependent yield stress. In this study, we pursue this route to develop a theory and numerical results for the yielding of a rough, elastoplastic surface with scale-dependent yield stress. Here, we examine surfaces for which the power spectrum of the topography 𝐶 and yield stress 𝑌 follow power laws as a function of scale 𝜆, such that 𝐶∼𝜆 −𝑚 and 𝑌∼𝜆 −𝑛 , respectively. In this formal treatment of the problem, we focus on surfaces in contact and the resulting yield and do not impose shear. Numerical solutions show that the deviation from the elastic scaling solution is bounded as expected by the prior 1D heuristic scaling argument that anticipates the Hurst exponent as 1−𝑛. We also show that the plasticity is expected to erode the contacts if 𝑚 is lower than 𝑛−3, which corresponds to a Hurst exponent lower than 1−𝑛/2. This result is rigorously sound for 2D, i.e., realistic surfaces, and quantitatively different than the prior scaling argument. The theory now permits a correspondingly quantitative approach to interpreting natural surfaces.

elasticity↗

Meteorological factors in Earth-Satellite propagation

A 5-meter paraboloidal antenna operated at 28 GHz showed gain changes of 2 dB due to rain. While precise estimation of the corresponding angle of arrival changes is difficult, they appear to have been on the order of 0.02 degrees. The attenuation at 28.6 GHz inferred from radiometry agreed well with that measured simultaneously over a satellite/Earth link at the same frequency. The radiometers so calibrated have been used to add to the available site diversity data base using a 9 km baseline. An improved empirical model of site diversity gain was obtained by applying regression techniques to available published data. A brief review of the literature has led to suggestions for two experimental programs, one dealing with multifrequency radiometry and the other with the effects of the stochastic properties of precipitation on wideband data transmission.

Levis, C. A.↗

Loss of Control Detection for Commercial Transports Using Conditional Variational Autoencoders

This work describes a detector for the loss of control of a commercial transport in flight. The detector has a belief state defined by the latent variable stochastic modeling of a conditional variational autoencoder (CVAE) constructed with bidirectional recurrent layers. In 2000, the Boeing Company and the NASA Langley Research Center jointly developed a quantitative set of metrics for defining loss-of-control (LOC) for a commercial transport. We use the thresholds for these quantitative metrics to define a condition vector for training the CVAE. We demonstrate through experimentation that reconstruction probability is an accurate indicator that the vehicle has shifted to an LOC state. Second, we introduce a technique for inferring that the vehicle is experiencing a flight state change is approaching by measuring a shift in the sampling Gaussian distributions of the latent space. We provide an analysis of its applicability to flight data from a NASA generic commercial transport-type aircraft.

Newton H Campbell↗

Edge supercurrent reveals competition between condensates in a Weyl superconductor

In topological materials, the edge states are readily distinguished from the bulk states. The situation where a topological semimetal becomes superconducting so that Cooper pairs occupy both the bulk and the edge states is not well understood. In particular, we do not know if we can force their pairing symmetries to be different. Here we show that, when supercurrent is injected into the superconducting Weyl semimetal MoTe 2 from Nb contacts, the invasive s-wave pairing potential from Nb is incompatible with the intrinsic Cooper pair condensate in MoTe 2 . This incompatibility leads to strong stochasticity in the switching current and an unusual anti-hysteretic behaviour in the current–voltage loops. There is also an asymmetry in the edge oscillations where, as the magnetic field crosses zero, the phase noise switches from one with a noisy spectrum to one that is noise free. Furthermore, using the noise spectrum as a guide, we track the anomalous features to field-induced switching of the device gap function between s-wave symmetry and the unconventional symmetry intrinsic to MoTe 2 . We infer that the behaviour of the gap function along the edges is different from that in the bulk.

Superconducting properties and materials↗

Filament-induced breakdown spectroscopy with structured beams

Filament-induced ablation represents an attractive scheme for long-range material identification via optical spectroscopy. However, the delivery of laser energy to the target can be severely hindered by the stochastic nature of multiple-filamentation, ionization of ambient gas, and atmospheric turbulence. In order to mitigate some of these adverse effects, we examine the utility of beam shaping for femtosecond filament-induced breakdown spectroscopy with Gaussian and structured (Laguerre-Gaussian, Airy, and Bessel-Gaussian) beams in the nonlinear regime. Interaction of filaments with copper, zinc, and brass targets was studied by recording axially-resolved broadband emission from the filament-induced plasma. The laser-solid coupling efficacy was assessed by inferring thermodynamic parameters such as excitation temperature and electron density. While under our experimental conditions the ablation rate with Gaussian- and Laguerre-Gaussian beams is found to be similar, the Airy and Bessel-Gaussian beams offer the advantage of longitudinally extended working zones. These results provide insights into potential benefits of structuring ultrafast laser beams for standoff sensing applications.

47 OTHER INSTRUMENTATION↗

Loss of Control Detection for Commercial Transport Aircraft Using Conditional Variational Autoencoders

This work describes a detector for the loss of control condition of a commercial transport in flight. The detector has a belief state defined by the latent variable stochastic modeling of a conditional variational autoencoder (CVAE) constructed with bidirectional recurrent layers. In 2000, the Boeing Company and the NASA Langley Research Center jointly developed a quantitative set of metrics for defining loss-of-control (LOC) for a commercial transport. We use the thresholds for these quantitative metrics to define a condition vector for training the CVAE. First, we demonstrate through experimentation that reconstruction probability is an accurate indicator that the vehicle has shifted to an LOC state. Second, we introduce a technique for inferring that the vehicle is approaching a flight state change by measuring a shift in the sampling distributions of the CVAE latent space. The sampling distributions for flight observations that are approaching envelope limits are localized to external areas of the latent space. We provide an analysis of its applicability to flight data from NASA’s dynamically-scaled generic transport model (GTM) aircraft.

Loss of control↗

NESSUS/expert and NESSUS/FPI in the Probabilistic Structural Analysis Methods (PSAM) program

The Numerical Evaluation of Stochastic Structures under Stress (NESSUS) is the primary computer code being developed in the NASA Probabilistic Structural Analysis Methods (PSAM) project. It consists of four modules NESSUS/EXPERT, NESSUS/FPI, NESSUS/PRE and NESSUS/FEM. This presentation concentrates on EXPERT and FPI. To provide an effective interface between NESSUS and the user, an expert system module called NESSUS/EXPERT is being developed. That system uses the CLIPS artificial intelligence code developed to NASA-JSC. The code is compatible with FORTRAN, the standard language for codes in PSAM. The user interacts with the CLIPS inference engine, which is linked to the knowledge database. The perturbation database generated by NESSUS/FEM and managed in EXPERT is used to develop the so-called response or performance model in the random variables. Two independent probabilistic methods are available in PSAM for the computation of the probabilistic structural response. These are the Fast Probability Integration (FPI) method and Monte Carlo simulation. FPI is classified as an advanced reliability method and has been developed over the past ten years by researchers addressing the reliability of civil engineering structures. Monte Carlo is a well-established technique for computing probabilities by conducting a number of deterministic analyses with specified input distributional information.

Burnside, O. H.↗

Statistical Treatment of Inverse Problems Constrained by Differential Equations-based Models with Stochastic Terms

This paper introduces a statistical treatment of inverse problems constrained by models with stochastic terms. The solution of the forward problem is given by a distribution represented numerically by an ensemble of simulations. The goal is to formulate the inverse problem, in particular the objective function, to find the closest forward distribution (i.e., the output of the stochastic forward problem) that best explains the distribution of the observations in a certain metric. We use proper scoring rules, a concept employed in statistical forecast verification, namely energy, variogram, and hybrid (i.e., combination of the two) scores. We study the performance of the proposed formulation in the context of two applications: a coefficient field inversion for subsurface flow governed by an elliptic partial differential equation with a stochastic source and a parameter inversion for power grid governed by differential-algebraic equations. In both cases we show that the variogram and the hybrid scores produce better parameter inversion results than does the energy score, whereas the energy score leads to better probabilistic predictions.

97 MATHEMATICS AND COMPUTING↗

Extracting scene feature vectors through modeling, volume 3

The remote estimation of the leaf area index of winter wheat at Finney County, Kansas was studied. The procedure developed consists of three activities: (1) field measurements; (2) model simulations; and (3) response classifications. The first activity is designed to identify model input parameters and develop a model evaluation data set. A stochastic plant canopy reflectance model is employed to simulate reflectance in the LANDSAT bands as a function of leaf area index for two phenological stages. An atmospheric model is used to translate these surface reflectances into simulated satellite radiance. A divergence classifier determines the relative similarity between model derived spectral responses and those of areas with unknown leaf area index. The unknown areas are assigned the index associated with the closest model response. This research demonstrated that the SRVC canopy reflectance model is appropriate for wheat scenes and that broad categories of leaf area index can be inferred from the procedure developed.

Berry, J. K.↗