Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “probability distribution”

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,657 records · Page 92

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↗

Effect of tectonic processes on biosphere–geosphere feedbacks across a convergent margin

The subsurface is among Earth’s largest biomes, but the extent to which microbial communities vary across tectonic plate boundaries or interact with subduction-scale geological processes remains unknown. Here we compare bacterial community composition with deep-subsurface geochemistry from 21 hot springs across the Costa Rican convergent margin. We find that cation and anion compositions of the springs reflect the dip angle and position of the underlying tectonic structure and also correlate with the bacterial community. Co-occurring microbial cliques related to cultured chemolithoautotrophs that use the reverse tricarboxylic acid cycle (rTCA) as well as abundances of metagenomic rTCA genes correlate with concentrations of slab-volatilized carbon. This, combined with carbon isotope evidence, suggests that fixation of slab-derived CO2 into biomass may support a chemolithoautotrophy-based subsurface ecosystem. We calculate that this forearc subsurface biosphere could sequester 1.4 × 10 9 to 1.4 × 10 10 mol of carbon per year, which would decrease estimates of the total carbon delivered to the mantle by 2 to 22%. Furthermore, based on the observed correlations, we suggest that distribution and composition of the subsurface bacterial community are probably affected by deep tectonic processes across the Costa Rican convergent margin and that, by sequestering carbon volatilized during subduction, these chemolithoautotrophic communities could in turn impact the geosphere.

54 ENVIRONMENTAL SCIENCES↗

Fast and efficient identification of anomalous galaxy spectra with neural density estimation

ABSTRACT Current large-scale astrophysical experiments produce unprecedented amounts of rich and diverse data. This creates a growing need for fast and flexible automated data inspection methods. Deep learning algorithms can capture and pick up subtle variations in rich data sets and are fast to apply once trained. Here, we study the applicability of an unsupervised and probabilistic deep learning framework, the probabilistic auto-encoder, to the detection of peculiar objects in galaxy spectra from the SDSS survey. Different to supervised algorithms, this algorithm is not trained to detect a specific feature or type of anomaly, instead it learns the complex and diverse distribution of galaxy spectra from training data and identifies outliers with respect to the learned distribution. We find that the algorithm assigns consistently lower probabilities (higher anomaly score) to spectra that exhibit unusual features. For example, the majority of outliers among quiescent galaxies are E+A galaxies, whose spectra combine features from old and young stellar population. Other identified outliers include LINERs, supernovae, and overlapping objects. Conditional modelling further allows us to incorporate additional information. Namely, we evaluate the probability of an object being anomalous given a certain spectral class, but other information such as metrics of data quality or estimated redshift could be incorporated as well. We make our code publicly available.

Böhm, Vanessa↗

Generalized approach for rapid entropy calculation of liquids and solids

We build a comprehensive methodology for the fast computation of entropy across both solid and liquid phases. The proposed method utilizes a single trajectory of molecular dynamics (MD) to facilitate the calculation of entropy, which is composed of three components. The electronic entropy is determined through the temporal average acquired from density functional theory MD simulations. The vibrational entropy, typically the predominant contributor to the total entropy, even within the liquid state, is evaluated by computing the phonon density of states via the velocity autocorrelation function. The most arduous component to quantify, the configurational entropy, is assessed by probability analysis of the local structural arrangement and atomic distribution. We illustrate, through a variety of examples, that this method is both a versatile and valid technique for characterizing the thermodynamic states of both solids and liquids. Furthermore, this method is employed to expedite the calculation of melting temperatures, demonstrating its practical utility in computational thermodynamics.

36 MATERIALS SCIENCE↗

Robust VAR Capability Curve of DER with Uncertain Renewable Generation

Active distribution system with high penetration of inverter based distributed energy resources(DER), can be utilized for var-related ancillary services at the transmission side interface. In order to utilize the DER flexibility, transmission system operator must be presented the aggregated DER flexibility of distribution system. However, the uncertainty in renewable generation, questions the credibility of aggregated capability curve in practice. In this paper, we incorporate the uncertainty into aggregation process to develop capability curve while preserving the real physics (unbalance and lossy nature) of distribution system. The Resulting capability curve with the associated probability can be harnessed by the TSO for decision making for both planning and operation.

Kar, Aditya Shankar↗

Proper motions, cluster membership and reddening in NGC 6611

Approaches used in reductions for proper motion are based on methods developed by Vasilevskis et al. (1965). Probabilities of cluster membership were determined by fitting the distribution of stars in the proper motion vector diagram with two bivariate Gaussian distributions for the field and the cluster stars. The spectroscopic data indicate that the stellar content of the cluster is not anomalous. The color-magnitude diagram is considered along with questions regarding member stars lying above the main sequence, the cluster age, and aspects of spatial structure.

Kamp, L. W.↗

Optimal estimation for discrete time jump processes

Optimum estimates of nonobservable random variables or random processes which influence the rate functions of a discrete time jump process (DTJP) are derived. The approach used is based on the a posteriori probability of a nonobservable event expressed in terms of the a priori probability of that event and of the sample function probability of the DTJP. Thus a general representation is obtained for optimum estimates, and recursive equations are derived for minimum mean-squared error (MMSE) estimates. In general, MMSE estimates are nonlinear functions of the observations. The problem is considered of estimating the rate of a DTJP when the rate is a random variable with a beta probability density function and the jump amplitudes are binomially distributed. It is shown that the MMSE estimates are linear. The class of beta density functions is rather rich and explains why there are insignificant differences between optimum unconstrained and linear MMSE estimates in a variety of problems.

Vaca, M. V.↗

Flight and wind-tunnel correlation of boundary-layer transition on the AEDC transition cone

Transition and fluctuating surface pressure data were acquired on a 10 deg included angle cone, using the same instrumentation and technique over a wide range of Mach and Reynolds numbers in 23 wind tunnels and in flight. Transition was detected with a traversing pitot-pressure probe in contact with the surface. The surface pressure fluctuations were measured with microphones set flush in the cone surface. Good correlation of end of transition Reynolds number RE(T) was obtained between data from the lower disturbance wind tunnels and flight up to a boundary layer edge Mach number, M(e) = 1.2. Above M(e) = 1.2, however, this correlation deteriorates, with the flight Re(T) being 25 to 30% higher than the wind tunnel Re(T) at M(e) = 1.6. The end of transition Reynolds number correlated within + or - 20% with the surface pressure fluctuations, according to the equation used. Broad peaks in the power spectral density distributions indicated that Tollmien-Schlichting waves were the probable cause of transition in flight and in some of the wind tunnels.

Fisher, D. L.↗

Flight and wind-tunnel correlation of boundary-layer transition on the AEDC transition cone

Transition and fluctuating surface pressure data were acquired on a 10 degree included angle cone, using the same instrumentation and technique over a wide range of Mach and Reynolds numbers in 23 wind tunnels and in flight. Transition was detected with a traversing pitot pressure probe in contact with the surface. The surface pressure fluctuations were measured with microphones set flush in the cone surface. Good correlation of end of transition Reynolds number Re (sub T) was obtained between data from the lower disturbance wind tunnels and flight up to a boundary layer edge Mach number, M (sub e) = 1.2. Above M (sub e) = 1.2, however, this correlation deteriorates, with the flight Re (sub T) being 25 to 30% higher than the wind tunnel Re (sub T) at M (sub e) = 1.6. The end of transition Reynolds number correlated within + or - 20% with the surface pressure fluctuations. Broad peaks in the power spectral density distributions indicated that Tollmien-Schlichting waves were the probable cause of transition in flight and in some of the wind tunnels.

Fisher, D. F.↗

Spatial structure of high energy photon sources in solar flares

Stereoscopic observations of high energy (greater than about 100 keV) photon emission from five solar flares have been made with the X-ray spectrometers aboard the ISEE-3 and Pioneer Venus Orbiter spacecraft. The observed altitude structure of the photon source and its dependence on the photon energy and time during a flare are compared with the predictions of thermal and non-thermal models of the hard X-ray source. In the case of the impulsive source, it is found that (1) the thermal model with adiabatic compression and expansion of a magnetically confined plasma and the thin target (non-thermal) model are not consistent with the observations; (2) the thick target (non-thermal) model and the dissipative thermal model are partially in agreement with the observations; (3) the emission probably originates in many individual non-thermal sources distributed in altitude, the lower altitude sources being brighter than those at higher altitude. In the case of the gradual source, it is found that (1) models with purely coronal sources are not consistent with the observations; (2) a partial precipitation model with trapped as well as precipitating electrons is consistent with the observations.

Kane, S. R.↗

The spatial distribution of energetic ions and electrons in the magnetotail

The spatial distributions of energetic ion and electron bursts on the IMP 7 and 8 satellites in the earth's magnetotail were studied. It was found that the ion bursts were more frequently detected in the dusk than in the dawn quarter of the neutral sheet whereas the electron bursts, more frequently in the dawn than the dusk quarter. The degree of dawn-dusk asymmetry is however energy dependent; the distribution for higher energy particle bursts exhibits higher degree of asymmetry. The morphologies of the distributions manifest themselves as seasonal variations of the most probable solar ecliptic latitudes at which the ion and electron bursts were observed. The amplitudes of the variations are about 25 deg with the seasonal variation for ions leading that for electrons by about 2 months.

Xu, R.-L.↗

Load sharing in distributed real-time systems with state-change broadcasts

A decentralized dynamic load-sharing (LS) method based on state-change broadcasts is proposed for a distributed real-time system. Whenever the state of a node changes from underloaded to fully loaded and vice versa, the node broadcasts this change to a set of nodes, called a buddy set, in the system. The performance of the method is evaluated with both analytic modeling and simulation. It is modeled first by an embedded Markov chain for which numerical solutions are derived. The model solutions are then used to calculate the distribution of queue lengths at the nodes and the probability of meeting task deadlines. The analytical results show that buddy sets of 10 nodes outperform those of less than 10 nodes, and the incremental benefit gained from increasing the buddy set size beyond 15 nodes is insignificant. These and other analytical results are verified by simulation. The proposed LS method is shown to meet task deadlines with a very high probability.

Shin, Kang G.↗

Test particle calculations of pick-up ions in the vicinity of Comet Giacobini-Zinner

Ionization of cometary neutral molecules produces ions which are picked-up by the solar wind. The cometary ion pick-up process for Comet Giacobini-Zinner is studied using a test particle method in which trajectories are numerically calculated for several thousand ions whose initial locations were chosen randomly with probability proportional to the neutral density. The cometary ion distribution function was calculated at several locations upstream of the bow shock. Several types of magnetic fluctuations were included. Circularly and linearly polarized monochromatic waves result in less pitch-angle scattering of pick-up ions than turbulent fluctuations which have about the same amplitude. Some anisotropy remains in the distribution function just upstream of the shock, even for the turbulent fluctuations.

Cravens, T. E.↗

Comparison of laser transit anemometry data analysis techniques

Two techniques for the extraction of two-dimensional flow information from laser transit anemometry (LTA) data sets are presented and compared via a simulation study and experimental investigation. The methods are a probability density function (PDF) estimation technique and a marginal distribution analysis technique. The simulation study builds on the results of previous work and provides a quantification of the accuracy of both techniques for various LTA data acquisition scenarios. The experimental comparison consists of using an LTA system to survey the flow downstream of a turbulence generator in a small low-speed wind tunnel. The collected data sets are analyzed and compared.

Humphreys, William M., Jr.↗

Low redshift star-forming galaxies: What can they teach us about primeval galaxies?

The analysis of the UV plus optical spectra of three star-forming galaxies, Mrk 496, Mrk 357, TOL1924-416, obtained by matching the size of the optical aperture with that of IUE, has given unexpected results. These can be summarized as follows: (1) the dereddened Ly(alpha)/H(beta) ratios are consistent with the prediction of case B recombination for nebular emission, within the uncertainties; (2) the decrease of the Ly(alpha)/H(beta) ratio with increasing metallicities is not confirmed in our three objects, although the sample is too small to consider this result definitive. The first result is surprising, mainly because at least the two Markarian galaxies have a large enough H1 content to markedly increase the optical depth for the Ly(alpha) photons and to trigger their absorption by dust. This finding can probably be explained as an effect of the inhomogeneous distribution of gas and dust within the galaxies. On the basis of these results, we conclude that the detection of the Ly(alpha) emission line in searching for primeval galaxies (PG's) can be still considered a valid technique.

Calzetti, D.↗

Conflict Probability Estimation for Free Flight

The safety and efficiency of free flight will benefit from automated conflict prediction and resolution advisories. Conflict prediction is based on trajectory prediction and is less certain the farther in advance the prediction, however. An estimate is therefore needed of the probability that a conflict will occur, given a pair of predicted trajectories and their levels of uncertainty. A method is developed in this paper to estimate that conflict probability. The trajectory prediction errors are modeled as normally distributed, and the two error covariances for an aircraft pair are combined into a single equivalent covariance of the relative position. A coordinate transformation is then used to derive an analytical solution. Numerical examples and Monte Carlo validation are presented.

Paielli, Russell A.↗

Occurrency frequency of substorm field and plasma signatures observed near-earth by ISEE-1/2

The onset of the majority of substorms occurs when the tail field stops growing more tail-like and begins to become more dipolar. This corresponds to the onset signatures on the ground and in geosynchronous orbit. The AE indices and the IGS Pi 2 data were used to determine the major substorm onsets of 1978 and 1979. The time delay between successive substorms, the distribution of the substorm growth phase duration and the probability of tailward flows were determined as a function of spacecraft location. About a half of the substorms exhibit a plasma signature including earthward or tailward flows or plasma sheet drop out and recovery. Earthward flows are often seen at substorm onset, and almost always during substorm recovery. Tailward flows are occasionally seen at onset as the spacecraft is close enough to the neutral sheet. The experimental results are compared to predictions based on the neutral line and current sheet disruption models.

Hsu, T. S.↗

"A Richness Study of 14 Distant X-Ray Clusters from the 160 Square Degree Survey"

We have measured the surface density of galaxies toward 14 X-ray-selected cluster candidates at redshifts z(sub i) 0.46, and we show that they are associated with rich galaxy concentrations. These clusters, having X-ray luminosities of Lx(0.5-2 keV) approx. (0.5 - 2.6) x 10(exp 44) ergs/ sec are among the most distant and luminous in our 160 deg(exp 2) ROSAT Position Sensitive Proportional Counter cluster survey. We find that the clusters range between Abell richness classes 0 and 2 and have a most probable richness class of 1. We compare the richness distribution of our distant clusters to those for three samples of nearby clusters with similar X-ray luminosities. We find that the nearby and distant samples have similar richness distributions, which shows that clusters have apparently not evolved substantially in richness since redshift z=0.5. There is, however, a marginal tendency for the distant clusters to be slightly poorer than nearby clusters, although deeper multicolor data for a large sample would be required to confirm this trend. We compare the distribution of distant X-ray clusters in the L(sub X)-richness plane to the distribution of optically selected clusters from the Palomar Distant Cluster Survey. The optically selected clusters appear overly rich for their X-ray luminosities, when compared to X-ray-selected clusters. Apparently, X-ray and optical surveys do not necessarily sample identical mass concentrations at large redshifts. This may indicate the existence of a population of optically rich clusters with anomalously low X-ray emission, More likely, however, it reflects the tendency for optical surveys to select unvirialized mass concentrations, as might be expected when peering along large-scale filaments.

Jones, Christine↗