Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Probability density 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 757 records · Page 42

Evidence for electron energization accompanying spontaneous formation of ion acceleration regions in expanding plasmas

We report experiments conducted in an expanding argon plasma generated in the inductive mode of a helicon source in the Hot hELIcon eXperiment–Large Experiment on Instabilities and Anisotropies facility. As the neutral gas pressure increases, the supersonic ion acceleration weakens. Increasing neutral pressure also alters the radial profile of electron temperature, density, and plasma potential upstream of the plasma expansion region. Langmuir probe measurements of the electron energy probability function (EEPF) show that heating of electrons at the plasma edge by RF fields diminishes with increasing gas pressure, yielding a plasma with a centrally peaked electron temperature, and flat potential profiles at higher neutral pressures. For neutral pressures at which ion acceleration regions develop in the expanding plasma plume, EEPFs reveal electrons with two temperature components.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Simulation-based inference for parameter estimation of complex watershed simulators

High-resolution, spatially distributed process-based (PB) simulators are widely employed in the study of complex catchment processes and their responses to a changing climate. However, calibrating these PB simulators using observed data remains a significant challenge due to several persistent issues, including the following: (1) intractability stemming from the computational demands and complex responses of simulators, which renders infeasible calculation of the conditional probability of parameters and data, and (2) uncertainty stemming from the choice of simplified representations of complex natural hydrologic processes. Here, we demonstrate how simulation-based inference (SBI) can help address both of these challenges with respect to parameter estimation. SBI uses a learned mapping between the parameter space and observed data to estimate parameters for the generation of calibrated simulations. To demonstrate the potential of SBI in hydrologic modeling, we conduct a set of synthetic experiments to infer two common physical parameters – Manning's coefficient and hydraulic conductivity – using a representation of a snowmelt-dominated catchment in Colorado, USA. We introduce novel deep-learning (DL) components to the SBI approach, including an “emulator” as a surrogate for the PB simulator to rapidly explore parameter responses. We also employ a density-based neural network to represent the joint probability of parameters and data without strong assumptions about its functional form. While addressing intractability, we also show that, if the simulator does not represent the system under study well enough, SBI can yield unreliable parameter estimates. Approaches to adopting the SBI framework for cases in which multiple simulator(s) may be adequate are introduced using a performance-weighting approach. The synthetic experiments presented here test the performance of SBI, using the relationship between the surrogate and PB simulators as a proxy for the real case.

54 ENVIRONMENTAL SCIENCES↗

Random scattering approach to gravitational microlensing

Small random deflections of a narrow beam of radiation due to gravitational scattering by stars randomly distributed within the deflector plane are considered. Using a Fouriere transform method, the probability of scattering is obtained as a function of scattering angle for an arbitrary number of stars with an arbitrary distribution of masses. The probability density, expressed in proper units, depends on one parameter only: the effective number of stars. At small scattering angles the density is a Gaussian, and at large angles it falls off as the scattering angle to the minus fourth power. The probability distribution for scatterings is simply related to the angular distribution of the surface brightness of a macroimage, averaged over many microimages. The isophotes are ellipses, with the ratio of the major axis to the minor axis, determined by the dimensionless surface mass density and the shear of the lensing system. The number of stars that has to be included in the modeling of microlensing is proportional to the amplification due to the macrolens, and to the square of the dimensionless surface mass density.

Katz, N.↗

The influence of nonlinear damping on the random response of panels by time domain simulation

The effects of nonlinear damping on the large-deflection response of symmetrically laminated composite panels is studied by time-domain simulation in an attempt to quantify experimentally observed phenomena. A single-mode deflection approach is used in the equations of motion and in the linear and nonlinear damping models. The probability, peak-probability, and spectral-density data are compiled for both strains and deflection. The rms responses as a function of input excitation level are also obtained and compared to linear, equivalent-linearization, and Fokker-Planck equation solutions.

Robinson, J.↗

Gravitational microlensing of high-redshift supernovae by compact objects

An analysis of the effect of microlensing by a cosmologically dominant density of compact objects is performed, using high-redshift Type Ia supernovae (SN Ia's) as probes. The compact objects are modeled as a three-dimensional distribution of point masses, and Monte Carlo simulations are done to calculate the resulting amplification probability distributions for several column densities and cosmologies. By combining these distributions with the intrinsic SN Ia luminosity function and comparing with the results for a perfectly smooth universe, estimates are made of the number of supernovae that would need to be observed to confirm or rule out this lensing scenario. It is found that about 1000 SN Ia's with redshifts of z = 1 would be needed to perform this test, which is beyond what current searches can hope to accomplish. Observations of many fewer high-redshift supernovae, used merely as standard candles, appears a promising way of distinguishing between different cosmological models.

Rauch, Kevin P.↗

A count probability cookbok: Spurious effects and the scaling model

We study the errors brought by finite volume effects and dilution effects on the practical determination of the count probability distribution function P(sub N)(n,l), which is the probability of having N objects in a cell of volume l cubed for a set of average number density n. Dilution effects are particularly revelant to the so-called sparse sampling strategy. This work is mainly done in the framework of the Bailan & Schaeffer scaling model, which assumes that the Q-body correlation functions obey the scaling relation Xi(sub Q)(lambda r(sub l),....lambda r(sub Q) = lambda(exp -(Q-1)gamma) Xi(sub Q)(r(sub 1),....r(sub Q)). We use three synthetic samples as references to perform our analysis: a fractal generated by a Rayleigh-Levy random walk with approximately 3 x 10(exp 4) objects, a sample dominated by a spherical power-law cluster with approximately 3 x 10(exp 4) objects and a cold dark matter (CDM) universe involving approximately 3 x 10(exp 5) matter particles.

Colombi, S.↗

Pseudospin-doublet bands and Gallagher Moszkowski doublet bands in 100 Y

New transitions in neutron-rich 100 Y have been identified in a 9 Be + 238 U experiment with mass and Z gates to provide full fragment identification. These transitions and high spin levels of 100 Y have been investigated by analyzing the high statistics γ–γ–γ and γ–γ–γ–γ coincidence data from the spontaneous fission of 252 Cf at the Gammasphere detector array. Two new bands, 14 new levels, and 23 new transitions have been identified. The K π = 4 + new band decaying to a 1s isomeric state is assigned to be the high-K Gallagher-Moszkowski (GM) partner of the known K π = 1 + band, with the π5/2[522]Ⓧν3/2[411] configuration. This 4 + band is also proposed to be the pseudospin partner of the new K π = 5 + band with a 5 + π5/2[422] Ⓧ ν5/2[413] configuration, to form a π5/2[422] Ⓧ ν[3125/2,3/2] neutron pseudospin doublet. Here, constrained triaxial covariant density-functional theory and quantal particle rotor model calculations have been applied to interpret the band structure and available electromagnetic transition probabilities and are found to be in good agreement with experimental values.

100Y↗

The track-length extension fitting algorithm for energy measurement of interacting particles in liquid argon TPCs and its performance with ProtoDUNE-SP data

This paper introduces a novel track-length extension fitting algorithm for measuring the kinetic energies of inelastically interacting particles in liquid argon time projection chambers (LArTPCs). The algorithm finds the most probable offset in track length for a track-like object by comparing the measured ionization density as a function of position with a theoretical prediction of the energy loss as a function of the energy, including models of electron recombination and detector response. The algorithm can be used to measure the energies of particles that interact before they stop, such as charged pions that are absorbed by argon nuclei. The algorithm's energy measurement resolutions and fractional biases are presented as functions of particle kinetic energy and number of track hits using samples of stopping secondary charged pions in data collected by the ProtoDUNE-SP detector, and also in a detailed simulation. Additional studies describe the impact of the dE/dx model on energy measurement performance. The method described in this paper to characterize the energy measurement performance can be repeated in any LArTPC experiment using stopping secondary charged pions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A parametric study of the complexity of sequential decoders, volume 2

Plots are presented of complexity versus required signal energy per bit/single-sided noise spectral density to obtain an output probability of error per bit of 0.0001 for an ideal coherent PSK demodulation. Transfer function of convolutional codes and techniques of branch and reference phase synchronization are considered.

Huth, G. K.↗

The effects of transmitter/receiver clock time base instability on coherent communication system performance

The purpose of this paper is to present a model and associated analysis of the deleterious effects which both transmitter and receiver data clock time-base instability have on bit error probability. The types of time-base instability modeled and analyzed include bit jitter (epoch jitter), bit jitter rate (time interval jitter) and data asymmetry. The bit error probability is determined as a function of the statistical characteristics of the time-base instability, the data transition density, the signal-to-noise ratio in the bit rate bandwidth and the mechanizations of both the data bit detector and bit synchronizer. Numerical results are provided.

Chie, C. M.↗

Shield sizing and response equations

A consolidated list is presented of meteoroid debris shield equations which have been given in the referenced memorandums. In some cases, equations have been updated; thus, this memorandum supersedes reference 1. The equations are presented in two parts: (1) shield sizing equations which are used to produce preliminary estimates of shielding weights; and (2) response equations to describe the impact conditions (projectile size as a function of velocity, density, and impact angle) causing failure of a given shield that are to be used for probability analyses (such as in the modified BUMPER program). Specific equations are given that are applicable for the following types of shields: aluminum Whipple shields; Nextel multishock (MS) shields; and mesh double bumper (MDB) shields. These equations will be updated in the future as warranted by the results of additional HVI tests, analyses, and shield modeling.

Christiansen, Eric L.↗

Theoretical calculations of ion acceleration 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 pickup process for comet Giacobini-Zinner is studied in two ways: (1) with 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 and (2) with a quasi-linear diffusion model. The cometary ion distribution function was calculated with the test particle model at several locations upstream of the bow shock and for several types of magnetic fluctuations (or waves). These waves were allowed to propagate in both directions along the magnetic field at the Alfven speed. Both pitch angle scattering and energy diffusion are evident in the derived ion distributions. The monochromatic waves result in less ion acceleration than turbulent fluctuations with about the same amplitude. The calculated ion distribution functions are in reasonable agreement with the distributions measured in the vicinity of comet Giacobini-Zinner in 1985 by particle detectors on the ICE spacecraft when the ratio of power in sunward propagating Alfven waves to the power in antisunward propagating waves is assumed to lie between about 20% and 50%. However, the quasi-linear diffusion model results agree best with the measured distribution functions when the power ratio is only about 3 percent.

Mckenzie, M. L.↗

Quantum information approach to high energy interactions

High energy hadron interactions are commonly described by using a probabilistic parton model that ignores quantum entanglement present in the light-cone wave functions. Here, we argue that since a high energy interaction samples an instant snapshot of the hadron wave function, the phases of different Fock state wave functions cannot be measured—therefore the light-cone density matrix has to be traced over these unobservable phases. Performing this trace with the corresponding U(1) Haar integration measure leads to ‘Haar scrambling’ of the density matrix, and to the emergence of entanglement entropy. This entanglement entropy is determined by the Fock state probability distribution, and is thus directly related to the parton structure functions. As proposed earlier, at large rapidity η the hadron state becomes maximally entangled, and the entanglement entropy is S E ~η according to QCD evolution equations. When the phases of Fock state components are controlled, for example in spin asymmetry measurements, the Haar average cannot be performed, and the probabilistic parton description breaks down.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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↗

New machine learning techniques for simulation-based inference: InferoStatic nets, kernel score estimation, and kernel likelihood ratio estimation

We propose an intuitive, machine-learning approach to multiparameter inference, dubbed the InferoStatic Networks (ISN) method, to model the score and likelihood ratio estimators in cases when the probability density can be sampled but not computed directly. The ISN uses a backend neural network that models a scalar function called the inferostatic potential \varphi φ . In addition, we introduce new strategies, respectively called Kernel Score Estimation (KSE) and Kernel Likelihood Ratio Estimation (KLRE), to learn the score and the likelihood ratio functions from simulated data. We illustrate the new techniques with some toy examples and compare to existing approaches in the literature. We mention en passant some new loss functions that optimally incorporate latent information from simulations into the training procedure.

Kong, Kyoungchul↗

Codebase release 0.1 for infstat

We propose an intuitive, machine-learning approach to multiparameter inference, dubbed the InferoStatic Networks (ISN) method, to model the score and likelihood ratio estimators in cases when the probability density can be sampled but not computed directly. The ISN uses a backend neural network that models a scalar function called the inferostatic potential \varphi φ . In addition, we introduce new strategies, respectively called Kernel Score Estimation (KSE) and Kernel Likelihood Ratio Estimation (KLRE), to learn the score and the likelihood ratio functions from simulated data. We illustrate the new techniques with some toy examples and compare to existing approaches in the literature. We mention en passant some new loss functions that optimally incorporate latent information from simulations into the training procedure.

Kong, Kyoungchul↗

Observation of a driver gas-tangential discontinuity

A complete analysis of an interplanetary disturbance of Nov. 19, 1970 using the Apollo 12-SIDE (Suprathermal Ion Detector Experiment) is presented. The SIDE detectors were pointing at 26.3 degrees from the normal solar-wind direction during the observations. The data were least-squares fitted (using a parabolic hypersurface approximation) to a convected Maxwell-Boltzmann distribution function. The results of the fit combined with two other experiments showed a drastic change in the wind speed (from an steady 352 km/sec down to 219 km/sec), direction, and temperature. Except for a delta-function increase at the onset, the density remained constant. There was a considerable enhancement in the abundance of He and probably of heavier elements. The interplanetary magnetic field exhibited a jump of 21 gamma with a change in latitude from -56 to -76 degrees in solar ecliptic coordinates. It is concluded that the disturbance was due to the driver gas-tangential discontinuity of a solar flare-induced shock wave. The characteristic of the tangential discontinuity fit well with theoretical prediction.

Medrano, R. A.↗

The effects of transmitter/receiver clock time-base instability on coherent communication system performance

The purpose of this paper is to present a model and associated analysis of the deleterious effects which both transmitter and receiver data clock time-base instability have on bit-error probability. The types of time-base instability modeled and analyzed include what is commonly referred to as bit jitter (epoch jitter), bit-jitter rate (time interval jitter), and data asymmetry. The bit-error probability is determined as a function of the statistical characteristics of the time-base instability, the data-transition density, the signal-to-noise ratio in the bit-rate bandwidth and the mechanizations of both the data-bit detector and bit synchronizer. The latter is manifested through the bit synchronizer's natural frequency and loop-damping parameters. Numerical results which indicate the sensitivity of the system bit-error rate to various parameters of the model are provided.

Chie, C. M.↗