SEARCH · Engineering Papers

Results for “Distribution 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,171 records · Page 65

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Separation of variables solution for non-linear radiative cooling

A separation of variables solution has been obtained for transient radiative cooling of an absorbing-scattering plane layer. The solution applies after an initial transient period required for adjustment of the temperature and scattering source function distributions. The layer emittance, equal to the instantaneous heat loss divided by the fourth power of the instantaneous mean temperature, becomes constant. This emittance is a function of only the optical thickness of the layer and the scattering albedo; its behavior as a function of these quantities is considerably different than for a layer at constant temperature.

Siegel, Robert↗

NASA NTRS1988

Modeling and measurement of error propagation in a multimodule computing system

An error propagation model has been developed for multimodule computing systems in which the main parameters are the distribution functions of error propagation times. A digraph model is used to represent a multimodule computing system, and error propagation in the system is modeled by general distributions of error propagation times between all pairs of modules. Two algorithms are developed to compute systematically and efficiently the distributions of error propagation times. Experiments are also conducted to measure the distributions of error propagation times with the fault-tolerant microprocessor (FTMP). Statistical analysis of experimental data shows that the error propagation times in FTMP do not follow a well-known distribution, thus justifying the use of general distributions in the present model.

Shin, Kang G.↗

NASA NTRS2008

A Probabilistic Wake Vortex Lateral Transport Model Using Data from SFO and DEN

In a previous report, we considered the behavior of the lateral position of vortices as a function of time after vortex formation for Out of Ground Effects (OGE) data for aircraft landing at San Francisco International Airport (SFO). We quantified the spread in lateral position as a function of time and examined how predictable lateral position is under a variety of assumptions. The combination of spread and predictability allowed us to derive probability distribution functions (PDFs) for lateral position given observed crosswind (CW) velocities. In this study, we examine the portability of these PDFs with respect to other landing sites. To this end, we consider OGE data obtained by the Federal Aviation Administration for landings at Denver International Airport (DEN) between 04/05/2006 and 06/03/2006. We consider vortices from both B733 (Boeing 737 models 200-500) and B757 (Boeing 757) aircraft. The data set contains 635 B733 landings and 506 B757 landings. The glide slope altitude for these measurements was 280 m, determined by the average initial vortex observation adjusted for a 3-second delay in the initial observation. The comparable SFO altitude was 158 m. We note that the principal mechanism for lateral transport in the OGE regime is advection by the ambient wind. This implies that a simple crosswind correction may be effective in explaining much of the variation in the lateral transport data. In this study, we again consider the use of ASOS data and average Lidar crosswind data over the vortex altitude range to predict vortex location as a function of time.

Mellman, George R.↗

NASA NTRS2022

Spectral Power-law Formation by Sequential Particle Acceleration in Multiple Flare Magnetic Islands

We present a first-principles model of pitch-angle and energy distribution function evolution as particles are sequentially accelerated by multiple flare magnetic islands. Data from magnetohydrodynamic (MHD) simulations of an eruptive flare/coronal mass ejection provide ambient conditions for the evolving particle distributions. Magnetic islands, which are created by sporadic reconnection at the self-consistently formed flare current sheet, contract and accelerate the particles. The particle distributions are evolved using rules derived in our previous work. In this investigation, we assume that a prescribed fraction of particles sequentially "hops" to another accelerator and receives an additional boost in energy and anisotropy. This sequential process generates particle number spectra that obey an approximate power law at mid-range energies and presents low- and high-energy breaks. We analyze these spectral regions as functions of the model parameters. We also present a fully analytic method for forming and interpreting such spectra, independent of the sequential acceleration model. The method requires only a few constrained physical parameters, such as the percentage of particles transferred between accelerators, the energy gain in each accelerator, and the number of accelerators visited. Our investigation seeks to bridge the gap between MHD and kinetic regimes by combining global simulations and analytic kinetic theory. The model reproduces and explains key characteristics of observed flare hard X-ray spectra as well as the underlying properties of the accelerated particles. Our analytic model provides tools to interpret high-energy observations for missions and telescopes, such as RHESSI, FOXSI, NuSTAR, Solar Orbiter, EOVSA, and future high-energy missions.

magnetic reconnection↗

DOE OSTI2025

Predicting the Galactic population of free-floating planets from realistic initial conditions

ABSTRACT We present the first prediction for the mass distribution function of Galactic free-floating planets (FFPs) that aims to accurately include the relative contributions of multiple formation pathways and stellar populations. We derive our predicted distribution from dedicated simulations of planet birth, growth, migration, and ejection around circumbinary systems and extend these results to also include the contributions from single and wide binary systems. Our resulting FFP mass distribution shows several distinct features, including a strong peak at ${\sim} 8\, {\rm M}_{\oplus }$ arising from the transition between pebble and gas accretion regimes and a trough at ${\sim} 1\, {\rm M}_{\oplus }$ due to the shift in the dominant ejection process from planet–planet scattering to ejection through interactions with stars in circumbinary systems. We find that interactions with the central binary in close circumbinary systems are likely the dominant progenitor for FFPs more massive than Earth, leading to a steep power-law dependence in mass that agrees well with existing observations. In contrast, we find planet–planet scattering events in single and wide binary systems likely produce the majority of planets at Mars mass and below, resulting in a shallower power-law dependence. Our results suggest that existing extrapolations into the sub-terrestrial mass range may significantly overestimate the true FFP abundance. The features we predict in the mass distribution of FFPs will be detectable by upcoming space-based microlensing surveys and, if observed, will provide key insight into the origins of FFPs and the environments in which they form.

Coleman, Gavin A. L. (ORCID:0000000151118963)↗

DOE OSTI2021

Deep-inelastic scattering with positron beams

Here we discuss several advantages of using positron beams in addition to electron beams for the study of parton distribution functions (PDFs) of the nucleon in deep-inelastic scattering. Specific topics include the flavor separation of the unpolarized PDFs as well as the spin-dependent distributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

fp-tools

The output of a molecular dynamics simulation is a data file containing properties of the motion of the simulated particles as a function of time (e.g., positions, velocities). From this output, scientists can extract various properties which give insight into the physics of the system; for example, the radial distribution function can be computed from the particle positions, which gives insight into how close the system is to melting or freezing. Many publicly available codes to perform molecular dynamics simulations exist and have been well adopted by the scientific community (e.g., LAMMPS, VASP). However, it remains standard practice for scientists to write their own post-processing scripts to extract various properties of interest from the molecular dynamics simulation data. Since in many cases scientists are interested in the same set of “textbook” properties, there is a large duplication of effort in writing these post-processing scripts. Our proposed code, FP-Tools, is a C++ toolkit designed to help scientists extract and analyze quantities of interest from the output of a molecular dynamics simulation, eliminating the need for them to write these codes themselves and thus reducing this duplication of effort. The quantities which are computed by our code are well documented in the literature (either in publications or textbooks), and the algorithms we implement to compute these properties are also well known in the field. We are not introducing new science or methods here; rather, our primary goal is to provide a useful tool to the community in the form of a well-documented, easy to use package.

Hartman, Leah↗

NASA NTRS1985

On the determination of the depth of EAS development maximum using the lateral distribution of Cerenkov light at distances 150 m from EAS axis

The Samarkand extensive air showers (EAS) array was used to measure the mean and individual lateral distribution functions (LDF) of EAS Cerenkov light. The analysis of the individual parameters b showed that the mean depth of EAS maximum and the variance of the depth distribution of maxima of EAS with energies of approx. 2x10 to the 15th power eV can properly be described in terms of Kaidalov-Martirosyan quark-gluon string model (QGSM).

Aliev, N.↗

NASA NTRS2019

Consistency Between NASS Surveyed Soil Moisture Conditions and SMAP Soil Moisture Observations

The USDA National Agricultural Statistics Survey (NASS) collects and publishes crop growth status and soil moisture conditions in major US agricultural regions. The operationally-produced weekly reports are based on survey information. The surveys are based on visual assessments and ? in the case of soil moisture ? report soil moisture levels in one of four categories (Very Short, Short, Adequate and Surplus). In this study, we show that these reports have remarkable correspondence with the NASA Soil Moisture Active Passive (SMAP) Level-4 Soil Moisture (L4SM) product. This consistency allows the combining the two distinct types of data to produce a value-added combination, which is mapped fields rather than State-by-State tables and it is refreshed daily rather than weekly. In this study classification thresholds are derived for L4SM by mapping cumulative distribution functions of L4SM surface and root-zone SM to the categorical NASS SM conditions. The results show that, year-over-year, the SMAP cumulative SM distributions are consistent with the NASS SM conditions and, furthermore, that the temporal evolution of the SMAP-derived thresholds is consistent with the seasonal crop growth cycles from year to year. The results signify that the SMAP SM retrievals are relatable to SM estimation conducted in agriculture by land managers and farmers, which underlines the general applicability of the SMAP data.

Andreas Colliander↗

NASA NTRS2006

The Solar Wind and Geomagnetic Activity as a Function of Time Relative to Corotating Interaction Regions

Corotating interaction regions during the declining phase of the solar cycle are the cause of recurrent geomagnetic storms and are responsible for the generation of high fluxes of relativistic electrons. These regions are produced by the collision of a high-speed stream of solar wind with a slow-speed stream. The interface between the two streams is easily identified with plasma and field data from a solar wind monitor upstream of the Earth. The properties of the solar wind and interplanetary magnetic field are systematic functions of time relative to the stream interface. Consequently the coupling of the solar wind to the Earth's magnetosphere produces a predictable sequence of events. Because the streams persist for many solar rotations it should be possible to use terrestrial observations of past magnetic activity to predict future activity. Also the high-speed streams are produced by large unipolar magnetic regions on the Sun so that empirical models can be used to predict the velocity profile of a stream expected at the Earth. In either case knowledge of the statistical properties of the solar wind and geomagnetic activity as a function of time relative to a stream interface provides the basis for medium term forecasting of geomagnetic activity. In this report we use lists of stream interfaces identified in solar wind data during the years 1995 and 2004 to develop probability distribution functions for a variety of different variables as a function of time relative to the interface. The results are presented as temporal profiles of the quartiles of the cumulative probability distributions of these variables. We demonstrate that the storms produced by these interaction regions are generally very weak. Despite this the fluxes of relativistic electrons produced during those storms are the highest seen in the solar cycle. We attribute this to the specific sequence of events produced by the organization of the solar wind relative to the stream interfaces. We also show that there are large quantitative differences in various parameters between the two cycles.

McPherron, Robert L.↗

DOE OSTI2022

Numerical schemes for a multi-species BGK model with velocity-dependent collision frequency

Here, we consider a kinetic description of multi-species gas mixture modeled with Bhatnagar-Gross-Krook (BGK) collision operators, in which the collision frequency varies not only in time and space but also with the microscopic velocity. In this model, the Maxwellians typically used in standard BGK operators are replaced by a generalization of such target functions, which are defined by a variational procedure. In this paper we present a numerical method for simulating this model, which uses an Implicit-Explicit (IMEX) scheme to minimize a certain potential function, mimicking the Lagrange functional that appears in the theoretical derivation. We show that theoretical properties such as conservation of mass, total momentum and total energy as well as positivity of the distribution functions are preserved by the numerical method, and illustrate its usefulness and effectiveness with numerical examples.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Theory and applications of parton pseudodistributions

We review the basic theory of the parton pseudodistributions approach and its applications to lattice extractions of parton distribution functions. The crucial idea of the approach is the realization that the correlator $M(z,p)$ of the parton fields is a function $M(ν,-z^2)$ of Lorentz invariants $ν=-(zp)$, the Ioffe time, and the invariant interval $z^2$. This observation allows to extract the Ioffe-time distribution $M(ν,-z^2)$ from Euclidean separations z accessible on the lattice. Another basic feature is the use of the ratio $M(ν,-z^2)≡M(ν,-z^2)/M(0,-z^2)$, that allows to eliminate artificial ultraviolet divergence generated by the gauge link for spacelike intervals. The remaining $z^2$-dependence of the reduced Ioffe-time distribution $M(ν,-z^2)$ corresponds to perturbative evolution, and can be converted into the scale-dependence of parton distributions $f(x,μ^2)$ using matching relations. The ν-dependence of $M(ν,-z^2)$ governs the x-dependence of parton densities $f(x,μ^2)$. The perturbative evolution was successfully observed in exploratory quenched lattice calculation. The analysis of its precise data provides a framework for extraction of parton densities using the pseudodistributions approach. It was used in the recently performed calculations of the nucleon and pion valence quark distributions. We also discuss matching conditions for the pion distribution amplitude and generalized parton distributions, the lattice studies of which are now in progress.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Magnetic Energy Dissipation and γ-Ray Emission in Energetic Pulsars

Abstract Some of the most energetic pulsars exhibit rotation-modulated γ -ray emission in the 0.1–100 GeV band. The luminosity of this emission is typically 0.1%–10% of the pulsar spin-down power ( γ -ray efficiency), implying that a significant fraction of the available electromagnetic energy is dissipated in the magnetosphere and reradiated as high-energy photons. To investigate this phenomenon we model a pulsar magnetosphere using 3D particle-in-cell simulations with strong synchrotron cooling. We particularly focus on the dynamics of the equatorial current sheet where magnetic reconnection and energy dissipation take place. Our simulations demonstrate that a fraction of the spin-down power dissipated in the magnetospheric current sheet is controlled by the rate of magnetic reconnection at microphysical plasma scales and only depends on the pulsar inclination angle. We demonstrate that the maximum energy and the distribution function of accelerated pairs is controlled by the available magnetic energy per particle near the current sheet, the magnetization parameter. The shape and the extent of the plasma distribution is imprinted in the observed synchrotron emission, in particular, in the peak and the cutoff of the observed spectrum. We study how the strength of synchrotron cooling affects the observed variety of spectral shapes. Our conclusions naturally explain why pulsars with higher spin-down power have wider spectral shapes and, as a result, lower γ -ray efficiency.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Neutron Structure Functions at Large Bjorken X

The E12-10-002 (F2) experiment ran in Hall C at JLab to extract the F2 structure functions of proton and neutron from inclusive H(e, e0) and D(e, e0) reactions in the high Bjorken x region. Those extractions put constraints on the Parton Distribution Functions PDFs (especially at large Bjorken x), and facilitate the study of Quark Hadron Duality. The F2 experiment measurements covers a large kinematic range in x from 0.2 to 1.0, and in Q2 from 4 to 16 GeV2 . While the proton F2 structure function has been studied extensively through inelastic electron-proton scattering, much less is known about neutron structure due to the unavailability of high density, free neutron targets. The BONuS12 experiment was proposed to measure the neutron F2 on a nearly free neutron within a weakly bound deuteron target via the spectator tagging method. Tagging the slow backward moving spectator protons minimizes both off-shell and final-state interaction effects, and the measured proton momentum is used to correct for the initial-state momentum of the struck neutron. The recoil detector is used to detect spectator protons with momenta 70 < ps < 150 MeV/c. This dissertation outlines the design, construction and testing of the recoil detector, Radial Time Projection Chamber (RTPC) and discusses the results of the H(e, e0) and D(e, e0) cross sections and the F2 structure functions obtained from the Hall C F2 data analysis.

Nadeeshani, Sooriyaarachchilage↗

DOE OSTI2021

Gluon distributions and their applications to Ioffe-time distributions

We investigate unpolarized and polarized gluon distributions and their applications to the Ioffe-time distributions, which are related to lattice QCD calculations of parton distribution functions. Guided by the counting rules based on the perturbative QCD at large momentum fraction x and the color coherence of gluon couplings at small x, we parametrize gluon distributions in the helicity basis. By fitting the unpolarized gluon distribution, the inferred polarized gluon distribution from our parametrization agrees with the one from global analysis. A simultaneous fit to both unpolarized and polarized gluon distributions is also performed to explore the model uncertainty. The agreement with the global analysis supports the (1 - x) power suppression of the helicity-antialigned distribution relative to the helicity-aligned distribution. The corresponding Ioffe-time distributions and their asymptotic expansions are calculated from the gluon distributions. Our results of the Ioffe-time distributions can provide guidance to the extrapolation of lattice QCD data to the region lacking precise gluonic matrix elements. Therefore, they can help regulate the ill-posed inverse problem associated with extracting the gluon distributions from discrete data from first-principle calculations, which are available in a limited range of the nucleon momentum and the spatial separation between the gluonic currents. Given various limitations in obtaining lattice QCD data at large Ioffe time, phenomenological approaches can provide important complementary information to extract the gluon distributions in the entire momentum fraction region, especially at small x. The possibility of investigating higher-twist effects and other systematic uncertainties in the contemporary first-principle calculations of parton distributions from phenomenologically well-determined Ioffe-time distributions in the large Ioffe-time region is also discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1993

Comparison of hydrodynamic and semi-kinetic treatments for plasma flow along closed field lines

Hydrodynamic and semi-kinetic treatments of plasma flow along closed geomagnetic field lines are compared. The hydrodynamic treatment is based on a simplified 16-moment set of transport equations as the equations for the heat flows are not solved; the heat flows are treated heuristically. The semi-kinetic treatment is based on a particle code. The comparison deals with the distributions of the plasma density, flow velocity, and parallel and perpendicular temperatures as obtained from the two treatments during the various stages of the flow. In the kinetic treatment, the appropriate boundary condition is the prescription of the velocity distribution functions for the particles entering the flux tubes at the ionospheric boundaries; those particles leaving the system are determined by the processes occurring in the flux tube. The prescribed distributions are half-Maxwellian with temperature T(sub 0) and density n(sub 0). In the hydrodynamic model, the prescribed boundary conditions are on density (n(sub 0)), flow velocity (V(sub 0)) and temperature (T(sub 0). It was found that results from the hydrodynamic treatment critically depend on V(sub 0); for early stages of the flow this treatment yields results in good agreement with those from the kinetic treatment, when V(sub 0) = square root of (kT(sub 0)/2 (pi)m), which is the average velocity of particles moving in a given direction for a Maxwellian distribution. During this early stage, the flows developing form the conjugate ionospheres show some distinct transitions. For the first hour or so, the flows are highly supersonic and penetrate deep into the opposite hemispheres, and both hydrodynamics and kinetic treatments yield almost similar features. It is found that during this period heatflow effects are negligibly small. When a flow penetrates deep into the opposite hemisphere, the kinetic treatment predicts reflection and setting up of counterstreaming. In contrast, the hydrodynamic treatment yields a shock in the flow. The reasons for this difference in the two treatments is discussed, showing that in view of the relatively warm ions, the coupling of ion beams and the consequent shock formation in the offequatorial region are not likely due to the enhancements in the beam temperatures. The counterstreaming in the kinetic treatment and the shock in the hydrodynamic treatment first advance upward to the equator and then downward to the ionospheric boundary from where the flow originated. The transit time for this advancement is found to be about 1 hour for the respective models. After 2 hours or so, both models predict that the flows from the ionospheric boundaries are generally subsonic with respect to the local ion-sound speed. At late stages of the flow, when a substantial fraction of ions entering the flux tube begin to return back in the kinetic treatment, the hydrodynamic treatment with the boundary condition V(sub 0) = square root of (kT(sub 0)/2(pi)m) yields an over-refilling, and the choice of V(sub 0) becomes uncertain.

Singh, Nagendra↗

NASA NTRS2013

Correlated Uncertainties in Radiation Shielding Effectiveness

The space radiation environment is composed of energetic particles which can deliver harmful doses of radiation that may lead to acute radiation sickness, cancer, and even death for insufficiently shielded crew members. Spacecraft shielding must provide structural integrity and minimize the risk associated with radiation exposure. The risk of radiation exposure induced death (REID) is a measure of the risk of dying from cancer induced by radiation exposure. Uncertainties in the risk projection model, quality factor, and spectral fluence are folded into the calculation of the REID by sampling from probability distribution functions. Consequently, determining optimal shielding materials that reduce the REID in a statistically significant manner has been found to be difficult. In this work, the difference of the REID distributions for different materials is used to study the effect of composition on shielding effectiveness. It is shown that the use of correlated uncertainties allows for the determination of statistically significant differences between materials despite the large uncertainties in the quality factor. This is in contrast to previous methods where uncertainties have been generally treated as uncorrelated. It is concluded that the use of correlated quality factor uncertainties greatly reduces the uncertainty in the assessment of shielding effectiveness for the mitigation of radiation exposure.

Werneth, Charles M.↗

DOE OSTI2024

Exploring the non-Gaussianity of the cosmic infrared background and its weak gravitational lensing

ABSTRACT Gravitational lensing deflects the paths of photons, altering the statistics of cosmic backgrounds and distorting their information content. We take the cosmic infrared background (CIB), which provides plentiful information about galaxy formation and evolution, as an example to probe the effect of lensing on non-Gaussian statistics. Using the Websky simulations, we first quantify the non-Gaussianity of the CIB, revealing additional detail on top of its well-measured power spectrum. To achieve this, we use needlet-like multipole-band filters to calculate the variance and higher-point correlations. Using our simulations, we show the two-, three- and four-point spectra, and compare our calculated power spectra and bispectra to Planck values. We then lens the CIB, shell-by-shell with corresponding convergence maps, to capture the broad redshift extent of both the CIB and its lensing convergence. The lensing of the CIB changes the three- and four-point functions by a few tens of per cent at large scales, unlike with the power spectrum, which changes by less than two per cent. We expand our analyses to encompass the full intensity probability distribution functions (PDFs) involving all n-point correlations as a function of scale. In particular, we use the relative entropy between lensed and unlensed PDFs to create a spectrum of templates that can allow estimation of lensing. The underlying CIB model is missing the important role of star bursting, which we test by adding a stochastic lognormal term to the intensity distributions. The novel aspects of our filtering and lensing pipeline should prove useful for any radiant background, including line intensity maps.

79 ASTRONOMY AND ASTROPHYSICS↗