Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “probability generating function”

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 325 records · Page 18

Discovering causal structure with reproducing-kernel Hilbert space ε -machines

We merge computational mechanics’ definition of causal states (predictively equivalent histories) with reproducing-kernel Hilbert space (RKHS) representation inference. The result is a widely applicable method that infers causal structure directly from observations of a system’s behaviors whether they are over discrete or continuous events or time. A structural representation—a finite- or infinite-state kernel ϵ-machine—is extracted by a reduced-dimension transform that gives an efficient representation of causal states and their topology. In this way, the system dynamics are represented by a stochastic (ordinary or partial) differential equation that acts on causal states. We introduce an algorithm to estimate the associated evolution operator. Paralleling the Fokker–Planck equation, it efficiently evolves causal-state distributions and makes predictions in the original data space via an RKHS functional mapping. We demonstrate these techniques, together with their predictive abilities, on discrete-time, discrete-value infinite Markov-order processes generated by finite-state hidden Markov models with (i) finite or (ii) uncountably infinite causal states and (iii) continuous-time, continuous-value processes generated by thermally driven chaotic flows. The method robustly estimates causal structure in the presence of varying external and measurement noise levels and for very high-dimensional data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Evaluation of Refractivity Profiles From Champ and SAC-C GPS Radio Occultation

The CHAMP and SAC-C missions are the first missions to carry a second-generation 'Blackjack' GPS receiver. One of the new features of this receiver is its ability to sense the lower troposphere closer to the surface than the proof-of-concept GPS/MET 1995 experiment. Since their launch, CHAMP and SAC-C have collected thousands of GPS radio occultations, representing a wealth of measurements available for data assimilation in Numerical Weather Prediction (NWP) models. In order to evaluate the refractivity data derived by JPL from raw radio occultation measurements, we use Data Assimilation Office (DAO) shout forecasts as an independent state of the atmosphere. We compare CHAMP and SAC-C refractivity (processed by JPL) with refractivity calculated from the DAO global fields of temperature, water vapor content and humidity. We will show statistics of the differences as well as Probability Density Functions (PDFs) of the differences. Depending upon availability of AIRS data, we plan to show individual profile comparisons between GPS radio occultation and AIRS retrievals.

Poli, Paul↗

Superconducting gap symmetry from Bogoliubov quasiparticle interference analysis on Sr 2 RuO 4

Here, the nature of the superconducting order parameter in Sr 2 RuO 4 has generated intense interest in recent years. Since the superconducting gap is very small, high-resolution methods such as scanning tunneling spectroscopy might be the best chance to directly resolve the gap symmetry. Recently, a Bogoliubov quasiparticle interference imaging (BQPI) experiment has suggested that the d x 2 –y 2 gap symmetry is appropriate for Sr 2 RuO 4 . In this work, we use a material-specific theoretical approach based onWannier functions of the surface of Sr 2 RuO 4 to calculate the continuum density of states as detected in scanning tunneling microscopy experiments. We examine several different proposed gap order parameters and calculate the expected BQPI pattern for each case. Comparing to the available experimental data, our results suggest that a s' + id xy gap order parameter is the most probable state, but the measured BQPI patterns still display features unaccounted for by the theory for any of the states currently under discussion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Evaluation of Refractivity Profiles from CHAMP and SAC-C GPS Radio Occultation

The CHAMP and SAC-C missions are the first missions to carry a second-generation 'Blackjack' GPS receiver. One of the new features of this receiver is its ability to sense the lower troposphere closer to the surface than the proof-of-concept GPS/MET 1995 experiment. Since their launch, CHAMP and SAC-C have collected thousands of GPS radio occultations, representing a wealth of measurements available for data assimilation in Numerical Weather Prediction (NWP) models. In order to evaluate the refractivity data derived by JPL from raw radio occultation measurements, we use Data Assimilation Office (DAO) 6-hour forecasts as an independent state of the atmosphere. We compare CHAMP and SAC-C refractivity (processed by JPL) with refractivity calculated from the DAO global fields of temperature, water vapor content and humidity. We will show statistics of the differences as well as Probability Density Functions (PDFs) of the differences. Depending upon availability of AIRS data, we plan to show individual profile comparisons between GPS radio occultation and AIRS retrievals.

Poli, Paul↗

Density estimation via measure transport: Outlook for applications in the biological sciences

Abstract One among several advantages of measure transport methods is that they allow or a unified framework for processing and analysis of data distributed according to a wide class of probability measures. Within this context, we present results from computational studies aimed at assessing the potential of measure transport techniques, specifically, the use of triangular transport maps, as part of a workflow intended to support research in the biological sciences. Scenarios characterized by the availability of limited amount of sample data, which are common in domains such as radiation biology, are of particular interest. We find that when estimating a distribution density function given limited amount of sample data, adaptive transport maps are advantageous. In particular, statistics gathered from computing series of adaptive transport maps, trained on a series of randomly chosen subsets of the set of available data samples, leads to uncovering information hidden in the data. As a result, in the radiation biology application considered here, this approach provides a tool for generating hypotheses about gene relationships and their dynamics under radiation exposure.

gene expression data↗

Ballistic Performance Model of Crater Formation in Monolithic, Porous Thermal Protection Systems

Porous monolithic ablative systems insulate atmospheric reentry vehicles from reentry plasmas generated by atmospheric braking from orbital and exo-orbital velocities. Due to the necessity that these materials create a temperature gradient up to several thousand Kelvin over their thickness, it is important that these materials are near their pristine state prior to reentry. These materials may also be on exposed surfaces to space environment threats like orbital debris and meteoroids leaving a probability that these exposed surfaces will be below their prescribed values. Owing to the typical small size of impact craters in these materials, the local flow fields over these craters and the ablative process afford some margin in thermal protection designs for these locally reduced performance values. In this work, tests to develop ballistic performance models for thermal protection materials typical of those being used on Orion are discussed. A density profile as a function of depth of a typical monolithic ablator and substructure system is shown in Figure 1a.

Miller, J. E.↗

Pristine Igneous Rocks and the Early Differentiation of Planetary Materials

Our studies are highly interdisciplinary, but are focused on the processes and products of early planetary and asteroidal differentiation, especially the genesis of the ancient lunar crust. Most of the accessible lunar crust consists of materials hybridized by impact-mixing. Rare pristine (unmixed) samples reflect the original genetic diversity of the early crust. We studied the relative importance of internally generated melt (including the putative magma ocean) versus large impact melts in early lunar magmatism, through both sample analysis and physical modeling. Other topics under investigation included: lunar and SNC (martian?) meteorites; igneous meteorites in general; impact breccias, especially metal-rich Apollo samples and polymict eucrites; effects of regolith/megaregolith insulation on thermal evolution and geochronology; and planetary bulk compositions and origins. We investigated the theoretical petrology of impact melts, especially those formed in large masses, such as the unejected parts of the melts of the largest lunar and terrestrial impact basins. We developed constraints on several key effects that variations in melting/displacement ratio (a strong function of both crater size and planetary g) have on impact melt petrology. Modeling results indicate that the impact melt-derived rock in the sampled, megaregolith part of the Moon is probably material that was ejected from deeper average levels than the non-impact-melted material (fragmental breccias and unbrecciated pristine rocks). In the largest lunar impacts, most of the impact melt is of mantle origin and avoids ejection from the crater, while most of the crust, and virtually all of the impact-melted crust, in the area of the crater is ejected. We investigated numerous extraordinary meteorites and Apollo rocks, emphasizing pristine rocks, siderophile and volatile trace elements, and the identification of primary partial melts, as opposed to partial cumulates. Apollo 15 sample 15434,28 is an extraodinarily large glass spherule, nearly if not entirely free of meteoritic contamination, and provides insight into the diversity of mare basalts in the Hadley-Apennine region. Apollo 14 sample 14434 is in many respects a new rock type, intermediate between nonmare gabbronorites and mare basalts. We helped to both plan and implement a consortium to study the Yamato-793605 SNC/martian meteorite.

Warren, Paul H.↗

Extending Newtonian Dynamics to Include Stochastic Processes

A paper presents further results of continuing research reported in several previous NASA Tech Briefs articles, the two most recent being Stochastic Representations of Chaos Using Terminal Attractors (NPO-41519), [Vol. 30, No. 5 (May 2006), page 57] and Physical Principle for Generation of Randomness (NPO-43822) [Vol. 33, No. 5 (May 2009), page 56]. This research focuses upon a mathematical formalism for describing post-instability motions of a dynamical system characterized by exponential divergences of trajectories leading to chaos (including turbulence as a form of chaos). The formalism involves fictitious control forces that couple the equations of motion of the system with a Liouville equation that describes the evolution of the probability density of errors in initial conditions. These stabilizing forces create a powerful terminal attractor in probability space that corresponds to occurrence of a target trajectory with probability one. The effect in configuration space (ordinary three-dimensional space as commonly perceived) is to suppress exponential divergences of neighboring trajectories without affecting the target trajectory. As a result, the post-instability motion is represented by a set of functions describing the evolution of such statistical quantities as expectations and higher moments, and this representation is stable.

Zak, Michail↗

Bayesian Estimation of Earth’s Undiscovered Mineralogical Diversity Using Noninformative Priors

Recently, statistical distributions have been explored to provide estimates of the mineralogical diversity of Earth, and Earth-like planets. In this paper, a Bayesian approach is introduced to estimate Earth’s undiscovered mineralogical diversity. Samples are generated from a posterior distribution of the model parameters using Markov chain Monte Carlo simulations such that estimates and inference are directly obtained. It was previously shown that the mineral species frequency distribution conforms to a generalized inverse Gauss–Poisson (GIGP) large number of rare events model. Even though the model fit was good, the population size estimate obtained by using this model was found to be unreasonably low by mineralogists. In this paper, several zero-truncated, mixed Poisson distributions are fitted and compared, where the Poisson-lognormal distribution is found to provide the best fit. Subsequently, the population size estimates obtained by Bayesian methods are compared to the empirical Bayes estimates. Species accumulation curves are constructed and employed to estimate the population size as a function of sampling size. Finally, the relative abundances, and hence the occurrence probabilities of species in a random sample, are calculated numerically for all mineral species in Earth’s crust using the Poisson-lognormal distribution. These calculations are connected and compared to the calculations obtained in a previous paper using the GIGP model for which mineralogical criteria of an Earth-like planet were given.

Bayesian statistics↗

Higher-Order Statistical Moments of Predicted Sonic Boom Waveforms Through Turbulence

Sonic booms generated by supersonic aircraft are affected by turbulence in the atmospheric boundary layer through which they propagate. Turbulence effects lead to random variability of the sonic boom waveforms measured on the ground, complicating the prediction of such waveforms. As an initial effort to predict the waveform variability, the solution of the coherent or mean sonic boom waveform has previously been formulated by the author. The current paper extends the formulation to derive an expression of the second-order statistical moment necessary to calculate the variance of the spectral amplitudes. Since the derivation uses a full wave equation, results using the derived expression will be compared with those obtained using a parabolic approximation. In order to fully quantify the uncertainty of our predictions, the higher order statistical moments are also formulated and are shown to be approximately zero. Consequently, the probability density function (pdf) of the spectral amplitudes is predicted to be Gaussian. The formulation is further extended to determine the pdf of the loudness of sonic booms and quantify the uncertainties associated with the loudness prediction. Results from the formulation are compared with available flight test data.

sonic booms↗

Voyager Observations of the Color of Saturn's Ring

Previously unreduced high resolution Voyager 2 images of Saturn's main rings are used to generate reflectivity (I/F) profiles as a function of radius. Ratios of profiles taken from green, violet, orange, and UV filter images are then produced. The I/F ratios are diagnostic of composition, and provide us with information on the rings' present state of compositional evolution. The rings are extremely reddish in color which suggests that they could not be pure water ice. The most likely candidates for the non-icy components are silicates and organics. The sources of these pollutants are of extreme importance in determining the compositional history of the rings. The radial profiles of ring color ratio exhibit several very interesting properties: (a) broad-scale, fairly smooth, color variations which are only weakly correlated with underlying ring structure (optical depth variations) across the outer C ring and inner B ring as well as the Cassini division region. These variations are probably consistent with ballistic transport; (b) fine-scale, noise-Like (but unquestionably real) color variations across at least the outer two-thirds of the B ring. Not only the "redness" but the spectral shape varies. These variations are currently unexplained. Groundbased spectroscopic observations should be pursued to study the implied compositional heterogeneities on at least the larger scales. This data set will be used for modeling of the color and composition of the main rings using ballistic transport and radiative transfer theories.

Estrada, Paul R.↗

The Information Length Concept Applied to Plasma Turbulence

A methodology to study statistical properties of anomalous transport in fusion plasma is investigated. Three time traces generated by the full-f gyrokinetic code GKNET are analyzed for this purpose. The time traces consist of heat flux as a function of the radial position, which is studied in a novel manner using statistical methods. The simulation data exhibit transport processes with both medium and long correlation length along the radius. A typical example of a phenomenon with long correlation length is avalanches. In order to investigate the evolution of the turbulent state, two basic configurations are studied, one flux-driven and one gradient-driven with decaying turbulence. The information length concept in tandem with Boltzmann–Gibbs and Tsallis entropy is used in the investigation. It is found that the dynamical states in both flux-driven and gradient-driven cases are surprisingly similar, but the Tsallis entropy reveals differences between them. This indicates that the types of probability distribution function are nevertheless quite different since the higher moments are significantly different.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Term Matrix: a novel Gene Ontology annotation quality control system based on ontology term co-annotation patterns

Biological processes are accomplished by the coordinated action of gene products. Gene products often participate in multiple processes, and can therefore be annotated to multiple Gene Ontology (GO) terms. Nevertheless, processes that are functionally, temporally and/or spatially distant may have few gene products in common, and co-annotation to unrelated processes probably reflects errors in literature curation, ontology structure or automated annotation pipelines. We have developed an annotation quality control workflow that uses rules based on mutually exclusive processes to detect annotation errors, based on and validated by case studies including the three we present here: fission yeast protein-coding gene annotations over time; annotations for cohesin complex subunits in human and model species; and annotations using a selected set of GO biological process terms in human and five model species. For each case study, we reviewed available GO annotations, identified pairs of biological processes which are unlikely to be correctly co-annotated to the same gene products (e.g. amino acid metabolism and cytokinesis), and traced erroneous annotations to their sources. To date we have generated 107 quality control rules, and corrected 289 manual annotations in eukaryotes and over 52 700 automatically propagated annotations across all taxa.

59 BASIC BIOLOGICAL SCIENCES↗

A reexamination of two-fluid solar wind models

The two-fluid solar-wind equations have been solved by a method which is approximately 50 times faster than any previously developed, through the use of asymptotic expansions which are self-consistently iterated upon to find a solution that passes through the critical point. The energy assumptions in two-fluid solar-wind models are reexamined, and the conclusions are as follows: (1) proton thermal conduction may not be neglected, (2) the Coulomb logarithm must be calculated as a function of radius, and (3) the electron and proton temperatures at the base need not be equal, even when the time scale for energy exchange between the species is an order of magnitude smaller than the expansion time at the base. It is possible to reproduce reasonable quiet-time solar-wind parameters at 1 AU, but only if the proton temperature is approximately twice the electron temperature at 1 solar radius. This may indicate that extended proton heating is important in the outer solar corona. Winds with velocities at 1 AU of 450 km/s are generated without nonthermal energy deposition but require high proton temperatures as well as very low densities at the base. Higher-velocity solutions are not possible in a spherically symmetric geometry for reasonable particle fluxes at 1 AU, and it is suggested that these higher-velocity states probably require additional heating, acceleration mechanisms, or nonradial flow.

Nerney, S.↗

A flight-test and simulation evaluation of the longitudinal final approach and landing performance of an automatic system for a light wing loading STOL aircraft equipped with wing spoilers

As part of a comprehensive flight-test investigation of short takeoff and landing (STOL) operating systems for the terminal systems for the terminal area, an automatic landing system has been developed and evaluated for a light wing-loading turboprop-powered aircraft. An advanced digital avionics system performed display, navigation, guidance, and control functions for the test aircraft. Control signals were generated in order to command powered actuators for all conventional controls and for a set of symmetrically driven wing spoilers. This report describes effects of the spoiler control on longitudinal autoland (automatic landing) performance. Flight-test results, with and without spoiler control, are presented and compared with available (basically, conventional takeoff and landing) performance criteria. These comparisons are augmented by results from a comprehensive simulation of the controlled aircraft that included representations of the microwave landing system navigation errors that were encountered in flight as well as expected variations in atmospheric turbulence and wind shear. Flight-test results show that the addition of spoiler control improves the touchdown performance of the automatic landing system. Spoilers improve longitudinal touchdown and landing pitch-attitude performance, particularly in tailwind conditions. Furthermore, simulation results indicate that performance would probably be satisfactory for a wider range of atmospheric disturbances than those encountered in flight. Flight results also indicate that the addition of spoiler control during the final approach does not result in any measurable change in glidepath track performance, and results in a very small deterioration in airspeed tracking. This difference contrasts with simulations results, which indicate some improvement in glidepath tracking and no appreciable change in airspeed tracking. The modeling problem in the simulation that contributed to this discrepancy with flight was not resolved.

Brown, S. C.↗

Feeling the Strain: Quantifying Ligand Deformation in Photosynthesis

Structural distortion of protein-bound ligands can play a critical role in enzyme function by tuning the electronic and chemical properties of the ligand molecule. However, quantifying these effects is difficult due to the limited resolution of protein structures and the difficulty of generating accurate structural restrains for non-protein ligands. Here, we seek to quantify these effects through a statistical analysis of ligand distortion in Chlorophyll (Chl) proteins (CP), where ring deformation is thought to play a role in energy and electron transfer. To assess the accuracy of ring-deformation estimates from available structural data, we take advantage of the C 2 symmetry of Photosystem II (PSII), comparing ring-deformation estimates for equivalent sites both within and between 113 distinct X-ray and Cryogenic electron microscopy (CryoEM) PSII structures. Significantly, we find that several deformation modes exhibit considerable variability in predictions, even for equivalent monomers, down to 2 Å resolution, to an extent that probably prevents their utilization in optical calculations. We further find that refinement restrains play a critical role in determining deformation values to resolution as low as 2 Å. However, for those modes that are well-resolved in the structural data, ring deformation in PSII is strongly conserved across all species tested, from cyanobacteria to algae. Furthermore, these results highlight both the opportunities and limitations inherent in the structure-based analysis of the bioenergetic and optical properties of CPs and other protein-ligand complexes.

14 SOLAR ENERGY↗

Adding GPU Support to the Markov Chain Monte Carlo Code Catmip

In geophysics, we are confronted with many under-determined inverse problems. For example, all of our observations of earthquakes are made at the Earth’s surface. So, when we try to infer how slip during an earthquake evolves in space and time, we find that there are many potential slip histories that are consistent with our limited observations and our understanding of earthquake physics. One way to approach these problems is with Bayesian analysis which allows us to infer the ensemble of all potential slip models that satisfy the observations and our prior knowledge of earthquake physics. In Bayesian analysis, our prior knowledge is known as the prior probability density function or prior PDF, the fit to the data is known as the data likelihood, and the target PDF that satisfies both the prior PDF and data likelihood is known as the posterior PDF. However, simulating the posterior PDF typically requires using Markov Chain Monte Carlo (MCMC) to draw tens of billions of random realizations of earthquake slip models, which may not be computationally feasible. To make this and similar geophysical inversions computationally tractable, we developed the Cascading Adaptive Transitional Metropolis In Parallel (CATMIP) algorithm. CATMIP is an efficient parallel Markov Chain Monte Carlo (MCMC) sampler that is used for model fitting and uncertainty quantification in geophysics. Example use cases are earthquake rupture modeling, determining mineral composition on Mars, reconstructing the history of ocean salinity, and historical earthquake relocation. CATMIP employs many parallel instances of the Metropolis algorithm for sampling in a transitioning framework. Transitioning is a process in which a set of random samples at equilibrium with a known probability density function (PDF) are used as seeds for the Markov chains to sample successive target PDFs that incrementally move the distribution from the starting seeds to the final desired PDF that describes the relative plausibility of potential values for the model parameters. The algorithm is implemented as a Master-Worker model employing MPI for communication. The worker processes are loosely coupled with global parameters periodically optimized by the master process. This provides a very high amount of parallelism with little communication between updates. During the presentation we will discuss the history of the algorithm and elaborate the earthquake rupture modeling use case for the CATMIP package. Our first step toward GPU optimization was to optimize the code for the CPU. CPU profiling revealed that most of the compute time is spent in calls to level 2 BLAS routines and calls to GSL random number generators. We revised the algorithm to employ level 3 BLAS routines instead. In our presentation we will describe how this was accomplished. Adding GPU support to CATMIP consisted mostly of replacing the calls to GSL with calls to GPU vendor-provided library routines. A small number of loops were directly implemented in CUDA. In the presentation will provide implementation details. Finally, we will discuss methods for profiling and opportunities for further optimizing GPU execution. By creating a code with the flexibility to run on either a CPU or GPU architecture, CATMIP can be used on systems ranging from large CPU-based HPC environments to single servers with GPU acceleration and everything in between.

HECC↗

LeptonInjector and LeptonWeighter: A neutrino event generator and weighter for neutrino observatories

In this work, we present a high-energy neutrino event generator, called LeptonInjector, alongside an event weighter, called LeptonWeighter. Both are designed for large-volume Cherenkov neutrino telescopes such as IceCube. The neutrino event generator allows for quick and flexible simulation of neutrino events within and around the detector volume, and implements the leading Standard Model neutrino interaction processes relevant for neutrino observatories: neutrino-nucleon deep-inelastic scattering and neutrino-electron annihilation. In this paper, we discuss the event generation algorithm, the weighting algorithm, and the main functions of the publicly available code, with examples. Program summary: Program Titles: LeptonInjector and LeptonWeighter CPC Library link to program files: https://doi.org/10.17632/662gkpjfd9.1 Developer's repository links: https://github.com/icecube/LeptonInjector and https://github.com/icecube/LeptonWeighter Licensing provisions: GNU Lesser General Public License, version 3. Programming Language: C++11 External Routines: • Boost • HDF5 • nuflux (https://github.com/icecube/nuflux) • nuSQuIDS (https://github.com/arguelles/nuSQuIDS) • Photospline (https://github.com/icecube/photospline) • SuiteSparse (https://github.com/DrTimothyAldenDavis/SuiteSparse) Nature of problem: LeptonInjector: Generate neutrino interaction events of all possible topologies and energies throughout and around a detector volume. LeptonWeighter: Reweight Monte Carlo events, generated by a set of LeptonInjector Generators, to any desired physical neutrino flux or cross section. Solution method: LeptonInjector: Projected ranges of generated leptons and the extent of the detector, in terms of column depth, are used to inject events in and around the detector volume. Event kinematics follow distributions provided in cross section files. LeptonWeighter: Event generation probabilities are calculated for each Generator, which are then combined into a generation weight and used to calculate an overall event weight.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗