Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “exponential probability density”

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 37 records · Page 2

A Numerical Approach to Estimate the Ballistic Coefficient of Space Debris from TLE Orbital Data

Low Earth Orbit (LEO) is full of space debris, which consist of spent rocket stages, old satellites and fragments from explosions and collisions. As of 2009, more than 21,000 orbital debris larger than 10 cm are known to exist], and while it is hard to track anything smaller than that, the estimated population of particles between 1 and 10 cm in diameter is approximately 500,000, whereas small as 1 cm exceeds 100 million. These objects orbit Earth with huge kinetic energies speeds usually exceed 7 kms. The shape of their orbit varies from almost circular to highly elliptical and covers all LEO, a region in space between 160 and 2,000 km above sea level. Unfortunately, LEO is also the place where most of our active satellites are situated, as well as, International Space Station (ISS) and Hubble Space Telescope, whose orbits are around 400 and 550 km above sea level, respectively.This poses a real threat as debris can collide with satellites and deal substantial damage or even destroy them.Collisions between two or more debris create clouds of smaller debris, which are harder to track and increase overall object density and collision probability. At some point, the debris density couldthen reach a critical value, which would start a chain reaction and the number of space debris would grow exponentially. This phenomenon was first described by Kessler in 1978 and he concluded that it would lead to creation of debris belt, which would vastly complicate satellite operations in LEO. The debris density is already relatively high, as seen from several necessary debris avoidance maneuvers done by Shuttle, before it was discontinued, and ISS. But not all satellites have a propulsion system to avoid collision, hence different methods need to be applied. One of the proposed collision avoidance concepts is called LightForce and it suggests using photon pressure to induce small orbital corrections to deflect debris from colliding. This method is very efficient as seen from theoretical simulations, even few continuous mode 10 kW ground-based lasers, focused by 1.5 m telescopes with adaptive optics, were enough to prevent significant amount of the debris collisions. Simulations were done by propagating all space objects in LEO by 1 year into the future and checking whether the probability of collision was high. For those space objects different ground-based lasers were used to divert them, afterwards collision probabilities were reevaluated. However, the actual accuracy of the LightForce software, which has been developed at NASA AmesResearch Center, depends on the veracity of the input parameters, one of which is the objects ballistic coefficient. It is a measure of bodys ability to overcome air resistance, which has a significant impact on the debris in LEO, and thus it is responsible for the shape of the trajectory of the debris. Having the exact values of the ballistic coefficient would make significantly better collision predictions, unfortunately, we do not know what are the values for most of the objects.In this research, we were working with part of LightForce code, which estimates the ballistic coefficient from ephemerides. Previously used method gave highly inaccurate values, when compared to known objects, and it needed to be changed. The goal of this work was to try out a different method of estimating the ballistic coefficient and to check whether or not it gives noticeable improvements.

Collisions↗

Lagrangian conditional statistics and flow topology in edge plasma turbulence

Lagrangian statistics and particle transport in edge plasma turbulence are investigated using the Hasegawa–Wakatani model and its modified version. The latter shows the emergence of pronounced zonal flows. Different values of the adiabaticity parameter are considered. The main goal is to characterize the role of coherent structures, i.e., vortices and zonal flows, and their impact on the Lagrangian statistics of particles. Computationally intensive long time simulations following ensembles of test particles over hundreds of eddy turnover times are considered in statistically stationary turbulent flows. The flow topology is characterized using the Lagrangian Okubo–Weiss criterion in order to split the flow into topologically different domains. In elliptic and hyperbolic regions, the probability density functions (PDFs) of the residence time have self-similar algebraic decaying tails. However, in the intermediate regions, the PDFs exhibit exponentially decaying tails. Topologically conditioned PDFs of the Lagrangian velocity, and acceleration and density fluctuations are likewise computed. The differences between the classical Hasegawa–Wakatani system and its modified version are assessed, and the role of zonal flows is highlighted. Here, the density flux spectrum, which characterizes the contributions of different length scales, is studied, and its inertial scaling is found to be in agreement with predictions based on dimensional arguments. Analyzing the angular change of particle tracers at different time scales, corresponding to coarse grained curvature, completes the study, and these multiscale geometric statistics quantify the directional properties of the particle motion in different flow regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Radar scattering laws for the lunar surface

Bistatic radar observations of quasi-specular scatter from the lunar surface at 13 and 116 cm wavelengths have been expressed in terms of probability densities of the underlying surface slope distributions. These show best agreement with a Hagfors scattering function, although other forms having predominantly Gaussian or exponential characteristics also occur. Estimates of root mean square (rms) surface slope derived from these function range from 4 deg rms in maria to at least 8 deg rms in highlands. These are values appropriate to 25 m horizontal scales and represent averages over tens of square kilometers of surface area. The effective scattering roughness varies with wavelength, most strongly in maria, where its behavior is consistent with that expected from an immature, i.e., unsaturated cratered surface. There is very little variation in scattering roughness with wavelength in highlands area.

Simpson, R. A.↗

Modeling gap probability in discontinuous vegetation canopies

In the present model for the gap probability of a discontinuous vegetation canopy, the assumption of a negative exponential attenuation within individual plant canopies will yield a problem involving the distribution distances within canopies through which a ray will pass. If, however, the canopies intersect and/or overlap, so that foliage density remains constant within the overlap area, the problem can be approached with two types of approximations. Attention is presently given to the case of a comparison of modeled gap probabilities with those observed for a stand of Maryland pine, which shows good agreement for zenith angles of illumination up to about 45 deg.

Li, Xiaowen↗

Generalized master equation for particle transport in binary random media with renewal statistics

Particle transport in binary stochastic mixtures is classically modeled assuming Markovian or exponential mixing statistics but in many applications material memory invalidates the Markov assumption. For non-Markovian mixing characterized by alternating renewal processes, a transport-theoretic framework is presented that provides an exact description of transport in nonscattering random binary media with general non-exponential statistics. Our approach is to Markovianize the problem by augmenting the {material type, particle flux} state space with the age or distance from the last interface. A Chapman-Kolmogorov equation is formulated for the joint probability density of the material type, particle flux, and age, and subsequently reduced to a generalized Master equation (GME) in differential form. This constitutes the primary result of this work. A state-updating Monte Carlo algorithm consistent with the GME is developed and benchmarked against analytical solutions for multiple chord-length laws. For purely absorbing renewal statistical media, the GME reproduces analytical benchmarks for the equilibrium age distribution, interior mean/variance of material-conditioned fluxes, and boundary transmittance. Simulations further demonstrate that a Markov (exponential) approximation of non-exponential statistics can introduce large errors in transmittance and interior flux profiles. Lastly, the reintroduction of memory due to scattering is briefly addressed through heuristic considerations.

Fluctuations & noise↗

Comparison between mirror Langmuir probe and gas-puff imaging measurements of intermittent fluctuations in the Alcator C-Mod scrape-off layer

Statistical properties of the scrape-off layer plasma fluctuations are studied in ohmically heated plasmas in the Alcator C-Mod tokamak. For the first time, plasma fluctuations as well as parameters that describe the fluctuations are compared across measurements from a mirror Langmuir probe (MLP) and from gas-puff imaging (GPI) that sample the same plasma discharge. This comparison is complemented by an analysis of line emission time-series data, synthesized from the MLP electron density and temperature measurements. The fluctuations observed by the MLP and GPI typically display relative fluctuation amplitudes of order unity together with positively skewed and flattened probability density functions. Such data time series are well described by an established stochastic framework that models the data as a superposition of uncorrelated, two-sided exponential pulses. The most important parameter of the process is the intermittency parameter, $γ=τ_d/τ_w$, where τd denotes the duration time of a single pulse and $τ_w$ gives the average waiting time between consecutive pulses. In this study we show, using a new deconvolution method, that these parameters can be consistently estimated from different statistics of the data. We also show that the statistical properties of the data sampled by the MLP and GPI diagnostic are very similar. Finally, a synthetic GPI signal using only plasma parameters sampled by the MLP shows qualitatively different fluctuation statistics from the measured GPI signal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Interrelated structure of high altitude atmospheric profiles

A preliminary development of a mathematical model to compute probabilities of thermodynamic profiles is presented. The model assumes an exponential expression for pressure and utilizes the hydrostatic law and equation of state in the determination of density and temperature. It is shown that each thermodynamic variable can be factored into the produce of steady state and perturbation functions. The steady state functions have profiles similar to those of the 1962 standard atmosphere while the perturbation functions oscillate about 1. Limitations of the model and recommendations for future work are presented.

Engler, N. A.↗

A Simple, Velocity Dependent, Collision Probability Algorithm for Small Combined Hard Body Radii Close Approach Events

Over the last few years progress has been made in the velocity dependent collision probability problem. Recent progress in this problem has been approached via integration of the probability flux into the surface of the combined hard body object. The algorithm presented uses the surface flux approach and is designed to compute the collision probability rate between two objects given their time dependent states and state error covariance matrices. With regards to the computation of the collision probability (Pc) rate, two major differences exist between the development of this algorithm and the usual Pc algorithm. First, the shape of the at-risk volume is assumed to be a cube rather than a sphere. The size of the cube is chosen so that it circumscribes the usual hard body sphere chosen for the spherical Pc problem. This will result in half the length of a side of the cube being equal to the hard body radius (HBR) of the sphere. Second, it is assumed that the HBR of the cube is much smaller than the smallest combined position uncertainty (σ(sub min) > 5*HBR) at each time point of evaluation. This is necessary as the actual collision probability calculation is based on a first order, small variable expansion in the position components of the Gaussian probability density function. This leads to the collision probability rate being to second order in the combined HBR. The collision probability rate is in a concise, closed form containing exponential and error functions.

Frisbee, Joe↗

Simulation of flight maneuver-load distributions by utilizing stationary, non-Gaussian random load histories

Random numbers were generated with the aid of a digital computer and transformed such that the probability density function of a discrete random load history composed of these random numbers had one of the following non-Gaussian distributions: Poisson, binomial, log-normal, Weibull, and exponential. The resulting random load histories were analyzed to determine their peak statistics and were compared with cumulative peak maneuver-load distributions for fighter and transport aircraft in flight.

Leybold, H. A.↗

Eulerian Mapping Closure Approach for Probability Density Function of Concentration in Shear Flows

The Eulerian mapping closure approach is developed for uncertainty propagation in computational fluid mechanics. The approach is used to study the Probability Density Function (PDF) for the concentration of species advected by a random shear flow. An analytical argument shows that fluctuation of the concentration field at one point in space is non-Gaussian and exhibits stretched exponential form. An Eulerian mapping approach provides an appropriate approximation to both convection and diffusion terms and leads to a closed mapping equation. The results obtained describe the evolution of the initial Gaussian field, which is in agreement with direct numerical simulations.

He, Guowei↗

Highly ionized interstellar gas located in the Galactic disk and halo

High-resolution IUE absorption line spectra have been obtained for 40 distant stars in order to study the distribution of interstellar H I, Si IV, C IV, and N V in the Galactic disk and lower halo. Respective midplane densities of 2 x 10 to the -9th, 7 x 10 to the -9th, and 3 x 10 to the -9th are found for Si IV, C IV, and Ni V. Both column density and velocity data indicate that the highly ionized gas (HIG) is considerably more extended in directions away from the Galactic plane than is H I or Si II. The absorption-line velocities for the halo HIG are consistent with the notion that halo gas in the inner Galaxy rotates more slowly than gas in the underlying disk. The derived column densities suggest an exponential scale height for the HIG of about 3 kpc; however, a simple exponential distribution is a poor representation of the distribution of the gas. It is concluded that a full explanation of the origin of the halo HIG will probably require a blending of ideas from the Galactic fountain and the photoionized halo models.

Savage, Blair D.↗

Dispersion and Lifetime of the SO2 Cloud from the August 2008 Kasatochi Eruption

Hemispherical dispersion of the SO2 cloud from the August 2008 Kasatochi eruption is analyzed using satellite data from the Ozone Monitoring Instrument (OMI) and the Goddard Trajectory Model (GTM). The operational OMI retrievals underestimate the total SO2 mass by 20-30% on 8-11 August, as compared with more accurate offline Extended Iterative Spectral Fit (EISF) retrievals, but the error decreases with time due to plume dispersion and a drop in peak SO2 column densities. The GTM runs were initialized with and compared to the operational OMI SO2 data during early plume dispersion to constrain SO2 plume heights and eruption times. The most probable SO2 heights during initial dispersion are estimated to be 10-12 km, in agreement with direct height retrievals using EISF algorithm and IR measurements. Using these height constraints a forward GTM run was initialized on 11 August to compare with the month-long Kasatochi SO2 cloud dispersion patterns. Predicted volcanic cloud locations generally agree with OMI observations, although some discrepancies were observed. Operational OMI SO2 burdens were refined using GTM-predicted mass-weighted probability density height distributions. The total refined SO2 mass was integrated over the Northern Hemisphere to place empirical constraints on the SO2 chemical decay rate. The resulting lower limit of the Kasatochi SO2 e-folding time is approx.8-9 days. Extrapolation of the exponential decay back in time yields an initial erupted SO2 mass of approx.2.2 Tg on 8 August, twice as much as the measured mass on that day.

Krotkov, N. A.↗

Validation of the two-region Rossi-alpha model for reflected assemblies

Nuclear criticality safety, nonproliferation and safeguards, emergency response, and stockpile stewardship utilize estimates of the k eff multiplication factor. The value of k eff cannot be directly measured, but it can be inferred from the prompt neutron period. One modality of measuring the prompt neutron period is the Rossi-alpha method, which is predicated on fitting a histogram of detection time differences due to the nonrandom temporal distribution of same-fission-chain neutrons. Recent works have developed the motivation and theory to expand traditional one-region point kinetic models resulting in one-exponential histogram fits to two-region models that result in two-exponential fits. This paper validates the new two-region model using organic scintillator measurements of copper-reflected weapons-grade plutonium (0:83 ≤ k eff ≤ 0:94) and high-density- polyethylene-reflected highly enriched uranium (0:73 ≤ k eff ≤ 0:95). Furthermore, the results show that more thermal systems have shorter prompt neutron periods in the core region due to increased induced-fission probabilities for moderated neutrons. Here, a new parameter introduced by the two- region model is also shown to be correlated to the amount of reflection, and may be used to infer assembly properties such as type and amount of reflector.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Exponential Escape Rate of Filamentary Incubation in Mott Spiking Neurons

Mott materials such as vanadium oxides, when subject to a strong applied voltage, present an inhomogeneous insulator-to-metal transition with formation of metallic filaments within the insulating bulk. This property is enabling the development of compact and power-efficient neuromorphic devices known as Mott neurons. However, the nature of the transition has not been fully understood yet, as it may be attributed to different effects, including Joule self-heating and hot-carrier injection. Moreover, the experimental determination of the threshold voltage needed to induce the transition has proven to be challenging, as the transition becomes increasingly unpredictable when the threshold is approached. The physical understanding of these issues would not only deepen our understanding of Mott insulators, but would also be an important step toward the realization of neuromorphic devices based on such materials. In this work we use numerical simulations based on the Mott resistor network model to study the nature of the filament incubation and formation process. Here, we show that both electronic and thermal effects, in the form of current density focusing and Joule self-heating, respectively, contribute to the filamentary incubation and growth. Remarkably, we find that the percolation of the metallic filaments near the threshold is intrinsically stochastic, qualitatively similar to the familiar Arrhenius activated behavior and to the stochastic firing of biological neurons. More precisely, we characterize the filament percolation as a Poisson point process, which has the same probability distribution as mathematical models of neuronal firing with an exponential escape rate. Finally, we support the numerical simulation results by performing experiments in VO 2 that are in agreement with the exponential escape rate behavior. Thus, we establish a functionality of Mott insulators that opens a path toward implementing neuromorphic hardware with quantum materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Results of critical velocity experiments with barium, strontium, and calcium releases from CRRES satellite

As part of the NASA Combined Release and Radiation Effects Satellite (CRRES) chemical release program in September 1990, two Ba and also one each Sr and Ca canisters of a boron-titanium thermite mixture, which vaporizes the element on ignition, were released near perigee after dusk in the South Pacific to study the critical velocity effect proposed by Alfven. The critical velocities of these three elements are 2.7, 3.5, and 5.4 km/s respectively, all well below the orbital velocity of 9.4 km/s. On September 10, 1990, a Sr and Ba pair (G-13, or critical ionization velocity (CIV) I) was released near Rarotonga at approximately 515 km altitude in a background electron density of 3.4 x 10(exp 6)/cu cm. On September 14, 1990, G-14 or CIV II released a Ca and Ba pair west of New Caledonia near 595 km at an electron density of 1.5 x 10(exp 6)/cu cm. Ions of all three elements were observed with low-light level imagers from two aircraft after they had transited up the magnetic field lines into the sunlight. Emissions from the spherically expanding neutral gas shells below the solar terminator, observed with cameras filtered for the Ba(+) ion line at 4554 A and also in unfiltered imagers for approximately 15 s after release, are probably due to excitation by hot electrons created in the CIV process. The ions created clearly lost much of their energy, which we now show can be explained by elastic collisions: Ba(+) + O. Inventories of the observed ions indicate yields of 0.15% and 1.84% for Ba in the first and second experiments, 0.02% for Sr and 0.27% for Ca. Ionization from all the releases continued along the satellite trajectory much longer (greater than 45 s) than expected for a CIV process. The ion production along the satellite track versus time typically shows a rapid rise to a peak in a few seconds followed by an exponential decrease to a level essentially constant rate. The characteristic distances for CIV I and II are 47 and 62 km, respectively. We interpret the early time rise and exponential fall to be due to CIV ionization, of 0.014% (CIV I) and 0.40% (CIV II) for the Ba releases. The later ions produced at a constant rate probably have origins from other such processes as stripping and associative ionization collisions with atmospheric constituents primarily O, and charge exchange with O(+), He(+), and H(+). We suggest that the much larger Ba ionization rate in CIV II than CIV I is due to the fact that the release occurred in the peak Ca density where hot electrons were already present.

Wescott, E. M.↗

Coronal gas in the Galactic halo: ORFEUS observations of NGC 346 #1

We present high resolution (lambda/delta lambda = 3000) observations of the Small Magellanic Cloud (SMC) star NGC 346 #1 near the O VI resonance doublet (1031.9, 1037.6 A). The data were collected with the Berkeley extreme/far ultraviolet (EUV/FUV) spectrometer in the ORFEUS telescope aboard the space platform Astro-SPAS during the missiom of space shuttle Discovery flown in September 1993. The O VI features shows broad absorption centered near nu(sub LSR) is approximately equal to +140 km/sec. This is almost certainly associated with the OB cluster NGC 346 itself, which also shows deep absorption at Si IV and C IV at comparable velocities. An absorption feature corresponding to O VI near nu(sub LSR) = -230 km/sec is present in our data, but this feature is caused by contamination from a second star, probably an early B-type main sequence star at lower standard of rest (LSR) velocity, which is also present in the 20 min circular aperture. Between these velocity extrema lies one broad feature (b is approximately equal 70 km/sec) or multiple narrower features, with a combined column density of about 1.4 x 10(exp 14)/sq cm. Assuming an exponential distribution of O VI with the in-plane density of Jenkins, we infer a scale height, H(sub 0), of no more than approximately equal 1.2 kpc toward NGC 346 #1. Our O VI column density is about a factor of 4 below the predictions of cooling fountain models, but our low scale height is consistent with more recent predictions.

Hurwitz, M.↗

Nonperturbative and perturbative dynamics of a light QCD axion: Dark matter and the strong 𝐶⁢𝑃 problem

Considerable theoretical efforts have gone into expanding the reach of the quantum chromodynamics (QCD) axion beyond its canonical mass–decay-constant relation. The 𝑍 𝒩 QCD axion model reduces the QCD axion mass naturally, by invoking a discrete 𝑍 𝒩 symmetry through which the axion field is coupled to 𝒩 copies of the Standard Model. Before the QCD phase transition at temperature 𝑇 QCD , the 𝑍 𝒩 potential has a minimum at misalignment angle 𝜃 = 𝜋. At 𝑇 QCD , 𝜃 = 𝜋 becomes a maximum; the axion potential becomes exponentially suppressed and develops 𝒩 minima—only one of which actually solves the strong 𝐶⁢𝑃 problem. Before 𝑇 QCD , 𝜃 relaxes toward 𝜋. After 𝑇 QCD , the axion field starts from around the hilltop and may have sufficient kinetic energy to overcome the newly suppressed potential barriers. Such a field evolution leads to nonperturbative effects via the self-interactions near the hilltop, which can cause the exponential growth of fluctuations and backreaction on the coherent motion. This behavior can influence the relic density of the field and the minimum in which it settles. We conduct the first lattice simulations of the 𝑍 𝒩 QCD axion using 𝒞osmoℒattice to accurately calculate dark matter abundances and find nonperturbative dynamics reduce the abundance by up to a factor of two. We furthermore find that the probability of solving the strong 𝐶⁢𝑃 problem tends to diverge considerably from the naïve expectation of 1/𝒩.

Axions↗

DNTTD (Distributed Non-Negative Tensor Train Decomposition)

The era of exascale computing opens new venues for innovations and discoveries in many scientific, engineering, and commercial fields. However, with the exa flops also come the extra-large high-dimensional data generated by high performance computing. High-dimensional data is presented as multidimensional arrays, aka tensors. The presence of latent (not directly observable) structures in the tensor allows a unique representation and compression of the data by classical tensor factorization techniques. However, the classical tensor methods are not always stable or they can be exponential in their memory requirements, which makes them not suitable for high-dimensional tensors. Tensor train (TT) is a state-of-the-art tensor network introduced for factorization of high-dimensional tensors. TT transforms the initial high-dimensional tensor in a network of three-dimensional tensors that requires only a linear storage. Many real-world data, such as, density, temperature, population, probability, etc., are non-negative and for an easy interpretation, the algorithms preserving non-negativity are preferred. Here, we introduce a distributed non-negative tensor-train and demonstrate its scalability and the compression on synthetic and real world big datasets.

Bhattarai, Manish↗