Engineering Papers⌕ Search

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.

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

At least 1,063 records · Page 59

Disordering and H-embrittlement of Pb(110) surface using embedded atom method and molecular dynamics

The dynamics and structure of the crystal melt interface are of great importance in crystal growth, melting, and other properties of the solid-liquid interface. The Embedded Atom Method (EAM), a semiempirical method developed by Daw and Baskes for determining the energetics of atoms in a bulk environment, is a useful method in the computer simulation. The EAM functions of PB in conjunction with the MD code are employed to study melting and disordering of the Pb(100) surface. The EAM functions are fitted to the bulk experimental data at zero temperature, and it is not clear if it will behave correctly at higher temperatures. In fact, it is shown that use of a similar type of potential at higher temperatures might result in a negative coefficient of thermal expansion. The primary purpose here is to perform MD simulation combined with the EAM functions of PB to study melting of the Pb(100) surface. This provides an accurate test of the EAM functions at higher temperatures. In particular, we calculate the two dimensional structure factor, the two dimensional distribution function, and density profile as a function of temperature up to the melting point of lead.

Karimi, Majid↗

Generalized probability model for calculation of interference to the Deep Space Network due to circularly Earth-orbiting satellites

The probability of exceeding interference power levels and the duration of interference at the Deep Space Network (DSN) antenna is calculated parametrically when the state vector of an Earth-orbiting satellite over the DSN station view area is not known. A conditional probability distribution function is derived, transformed, and then convolved with the interference signal uncertainties to yield the probability distribution of interference at any given instant during the orbiter's mission period. The analysis is applicable to orbiting satellites having circular orbits with known altitude and inclination angle.

Ruggier, C. J.↗

Classification of multifrequency multilook synthetic aperture radar data

A technique for segmentation of multifrequency multilook intensity synthetic aperture radar (SAR) data into regions of homogeneous and similar backscatter characteristics is presented. Two statistical models, one for the multifrequency multilook SAR intensities and the other for the distribution of the region labels, are combined to obtain the a posteriori probability distribution function of the region labels given the multifrequency speckled intensities. As the maximization of the posterior distribution is computationally intensive, a suboptimal technique for region labeling is proposed. Several examples using both simulated and real multifrequency multilook imagery are given to illustrate the performance of the algorithm.

Rignot, Eric↗

Plasma heating by collisionless magnetic reconnection - Analysis and computation

Analytic and numerical results on particle acceleration in 2D collisionless magnetic reconnection are presented. The particles are followed until they reach an outgoing flux surface at the same distance from the origin as the starting surface. The magnetic moment is not conserved for particles passing through the unmagnetized region around the X line at the origin. Other particles cross the separatrix without passing near the X line. The magnetic moment of the first class of outgoing particles is randomized, whereas it can be considered for the second class. The analytic model is based upon the observation of the final kinetic energy as a function of the initial conditions. Analytic results are shown to predict a Maxwellian tail for the distribution function in the perpendicular kinetic energy, with this energy much greater than the parallel kinetic energy. Numerical results showing that the predicted tail temperature agrees with the numerically computed temperature to within 10 percent over 4 orders of magnitude in the electric field are presented.

Moses, R. W.↗

A comparative study of turbulence decay using Navier-Stokes and a discrete particle simulation

A comparative study of the two dimensional temporal decay of an initial turbulent state of flow is presented using a direct Navier-Stokes simulation and a particle method, ranging from the near continuum to more rarefied regimes. Various topics related to matching the initial conditions between the two simulations are considered. The determination of the initial velocity distribution function in the particle method was found to play an important role in the comparison. This distribution was first developed by matching the initial Navier-Stokes state of stress, but was found to be inadequate beyond the near continuum regime. An alternative approach of using the Lees two-sided Maxwellian to match the initial strain-rate is discussed. Results of the comparison of the temporal decay of mean kinetic energy are presented for a range of Knudsen numbers. As expected, good agreement was observed for the near continuum regime, but the differences found for the more rarefied conditions were unexpectedly small.

Goswami, A.↗

Research on Spectroscopy, Opacity, and Atmospheres

With this funding I produced a web site kurucz.harvard.edu that can also be accessed by FTP. it has a 73GB disk that holds all of my atomic and diatomic molecular data, my tables of distribution function opacities, my grids of model atmospheres, colors, fluxes, etc., my programs that are ready for distribution, and most of my recent papers. Atlases and computed spectra will be added as they are completed. New atomic and molecular calculations will be added as they are completed.

Kurucz, Robert L.↗

TROPOMI Geometry-dependent Lambertian-Equivalent surface Reflectivity (GLER) Product for Improved Trace-Gas Retrieval

Accurate information about the reflectivity of the Earth's surface is required for most satellite retrievals of atmospheric composition, and this information is generally taken from monthly surface reflectivity climatology that neglects angular dependence. Previously we introduced Geometry-dependent Lambertian-equivalent surface reflectivity (GLER) which captures solar and satellite viewing angle dependence as well as daily and seasonal changes. GLER is calculated from simulations of Rayleigh-only top-of-atmosphere (TOA) radiances over non-Lambertian surfaces. We use NASA's Moderate Resolution Imaging Spectroradiometer (MODIS) bidirectional reflectance distribution function (BRDF) product over land and the wind-dependent Cox–Munk wave-facet slope distribution including water-leaving radiance over water to accounts for surface BRDF. We have developed global GLER product, previously for the Ozone Monitoring Instrument (OMI) and recently for Sentinel-5 Precursor (S5P) TROPOspheric Monitoring Instrument (TROPOMI) with several new improvements and updates. We have implemented the near real time daily V006 MODIS MCD43C1 BRDF data and gap-filled with a daily BRDF coefficient climatology created from 2002-2017 V006 MCD43GF data. The NASA’s Global Modelling Initiative hourly 0.25 x 0.25 deg Replay simulations are used for more accurate determination of pixel specific terrain pressure. To improve detection of seasonal snow/ice scenes, we use the 4-km snow cover product from the Interactive Multi-sensor Snow and Ice Mapping System (IMS). Finally, we use an improved version of the vector linearized discrete ordinate radiative transfer (VLIDORT) for update of the top-of-atmosphere (TOA) radiance look-up-tables (LUTs). We demonstrate how the use of GLER is beneficial to TROPOMI’s high spatial resolution (up to 3.5 km x 3.5 km) measurements to monitor atmospheric trace gas pollutants down to the sub-city scale.

TROPOMI↗

F2 extraction from Inclusive Cross Section Data at Large Bjorken x

While we have made significant progress in probing nuclear structure at low Bjorken x and high Q 2 , we still have gaps in our knowledge regarding the high Bjorken x and intermediate Q 2 kinematic region. This is not an accident, but instead due to the diffculty in assessing the non-perturbative region of nuclear physics. For this reason, the E12-10-002 experiment measured the H(e,e') and D(e,e') inclusive cross section in the resonance region. This will allow us to study both perturbative and non-perturbative physics in a kinematic region that lacks precision measurements. Measurements were made up to a Bjorken x of 0.99 and Q 2 up to 17 GeV 2 . These measurements are made using the independent HMS and SHMS spectrometers at the upgraded beam energy of 11 GeV in Hall C at Jefferson Lab. Further, we extract F 2 structure functions from the measured data and study the effect of their inclusion in a global Parton Distribution Function analysis.

Araiza Gonzalez, Fernando↗

Studies in Astronomical Time Series Analysis: VII. An Enquiry Concerning Non-Linearity, the RMS-Mean Flux Relation, and log-Normal Flux Distributions

A broad and widely used class of stationary, linear, additive time series models can have statistical properties which many authors have asserted imply that the underlying process must be non-linear, non-stationary, multiplicative, or inconsistent with shot noise. This result is demonstrated with exact and numerical evaluation of the model flux distribution function and dependence of flux standard deviation on mean flux (here and in the literature called the rms-flux relation). These models can: (1) exhibit normal, log-normal or other flux distributions; (2) show linear or slightly non-linear rms-mean flux dependencies; as well as (3) match arbitrary second order statistics of the time series data. Accordingly the above assertions cannot be made on the basis of statistical time series analysis alone. Also idscussed are ambiguities in the meaning of terms relevant to this study – linear, stationary and multiplicative – and functions that can transform observed fluxes to a normal distribution as well or better than the logarithm.

Jeffrey D Scargle↗

A Study of Cirrus Ice Particle Size Distribution Using TC4 Observations

An analysis of two days of in situ observations of ice particle size spectra, in convectively generated cirrus, obtained during NASA s Tropical Composition, Cloud, and Climate Coupling (TC4) mission is presented. The observed spectra are examined for their fit to the exponential, gamma, and lognormal function distributions. Characteristic particle size and concentration density scales are determined using two (for the exponential) or three (for the gamma and lognormal functions) moments of the spectra. It is shown that transformed exponential, gamma, and lognormal distributions should collapse onto standard curves. An examination of the transformed spectra, and of deviations of the transformed spectra from the standard curves, shows that the lognormal function provides a better fit to the observed spectra.

Tian, Lin↗

Global QCD analysis of spin PDFs in the proton with high-𝑥 and lattice constraints

We perform a comprehensive global QCD analysis of spin-dependent parton distribution functions (PDFs), combining all available data on inclusive and semi-inclusive deep-inelastic scattering (DIS), as well as inclusive weak boson and jet production in polarized 𝑝⁢𝑝 collisions, simultaneously extracting spin-averaged PDFs and fragmentation functions. Including recent Jefferson Lab DIS data at high 𝑥, together with subleading power corrections to the leading-twist framework, allows us to verify the stability of the PDFs for 𝑊 2 ≥ 4 GeV 2 and quantify the uncertainties on the spin structure functions more reliably. We explore the use of new lattice QCD data on gluonic pseudo-Ioffe time distributions, which, together with jet production and high-𝑥 DIS data, improve the constraints on the polarized gluon PDF. The expanded kinematic reach afforded by the data into the high-𝑥 region allows us to refine the bounds on higher-twist contributions to the spin structure functions, and test the validity of the Bjorken sum rule.

Perturbative QCD↗

Numerical models for Halley dust environments

The algebra, parameter values, and results of a reference model for the dust environment of Comet P/Halley are presented. The model employs simplified treatments of the icy conglomerate nucleus and of the hydrodynamic outflow of gas and dust. The model includes a gas production rate derived from observations of Halley's visual magnitude in 1910, gas and dust flux nonuniformly distributed over the nuclear surface, a dust distribution function consistent with large particle inferences from comet antitails, dust density variable with mass, and an approximation to the distribution of particles near the dust envelope. Many parameter values are derived by analogy with the results for other comets. The features are incorporated in a computer program which calculates the concentration, flux, and fluence of dust particles as functions of time, position and particle mass.

Divine, N.↗

Particlization in fluid dynamical simulations of heavy-ion collisions: The IS3D module

The IS3D particlization module simulates the emission of hadrons from heavy-ion collisions via Monte-Carlo sampling of the Cooper–Frye formula which converts fluid dynamical information into local phase-space distributions for hadrons. The code package includes multiple choices for the non-equilibrium correction to these distribution functions: the 14-moment approximation, first-order Chapman–Enskog expansion, and two types of modified equilibrium distributions. This makes it possible to explore to what extent heavy-ion experimental data are sensitive to different choices for δf n , presently the main source of theoretical uncertainty in the particlization stage. Here, we validate our particle sampler with a high degree of precision by generating several million hadron emission events from a longitudinally boost-invariant hypersurface and comparing the event-averaged particle spectra and space–time distributions to the Cooper–Frye formula.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Conservative DG method for the micro-macro decomposition of the Vlasov–Poisson–Lenard–Bernstein model

The micro-macro (mM) decomposition approach is considered for the numerical solution of the Vlasov– Poisson–Lenard–Bernstein (VPLB) system, which is relevant for plasma physics applications. In the mM approach, the kinetic distribution function is decomposed as $f$ = $\mathscr{E}$[$ρ_f$] + g , where $\mathscr{E}$ is a local equilibrium distribution, depending on the macroscopic moments $ρ_f$ = $∫_{\mathbb{R}}$ $efdv$ = $\langle$$ef$$\rangle$$_{\mathbb{R}}$, where e = (1, $v$, $\frac{1}{2}$ $v^2$) T , and $g$, the microscopic distribution, is defined such that $\langle$$eg$$\rangle$$_{\mathbb{R}}$ = 0. We aim to design numerical methods for the mM decomposition of the VPLB system, which consists of coupled equations for $ρ_f$ and $g$. To this end, we use the discontinuous Galerkin (DG) method for phase-space discretization, and implicit-explicit (IMEX) time integration, where the phase-space advection terms are integrated explicitly and the collision operator is integrated implicitly. We give special consideration to ensure that the resulting mM method maintains the $\langle$$eg$$\rangle$$_{\mathbb{R}}$ = 0 constraint, which may be necessary for obtaining (i) satisfactory results in the collision dominated regime with coarse velocity resolution, and (ii) unambiguous conservation properties. The constraint-preserving property is achieved through a consistent discretization of the equations governing the micro and macro components. Here, we present numerical results that demonstrate the performance of the mM method. The mM method is also compared against a corresponding DG-IMEX method solving directly for $f$.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Electron heating in kinetic-Alfvén-wave turbulence

We report analytical and numerical investigations of subion-scale turbulence in low-beta plasmas using a rigorous reduced kinetic model. We show that efficient electron heating occurs and is primarily due to Landau damping of kinetic Alfvén waves, as opposed to Ohmic dissipation. This collisionless damping is facilitated by the local weakening of advective nonlinearities and the ensuing unimpeded phase mixing near intermittent current sheets, where free energy concentrates. The linearly damped energy of electromagnetic fluctuations at each scale explains the steepening of their energy spectrum with respect to a fluid model where such damping is excluded (i.e., a model that imposes an isothermal electron closure). In conclusion, the use of a Hermite polynomial representation to express the velocity-space dependence of the electron distribution function enables us to obtain an analytical, lowest-order solution for the Hermite moments of the distribution, which is borne out by numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Models of the quasar population. II - The effects of dust obscuration

A model for dust obscuration is proposed to account for the apparent lack of high-redshift quasars which are assumed to be obscured by dust in intervening galaxies. The distribution function for the optical depth as a function of redshift is constructed and used to determine the expected properties of quasars. The apparent magnitude distribution at three points is used to fit the model. Optical results are given for the comoving number density, redshift distribution, and mean color excess. If the theory is correct, evidence for dust may be seen in quasar spectra and a correlation between quasar continuum slopes and strength of metal line absorption systems should be observed.

Heisler, Julia↗

Electrostatic charging of spacecraft in response to electron beam injection

Electron beam injections from spacecraft now constitute a major activity in space research. Here, the charging level of a conducting surface when an electron beam is injected from it is investigated. Injections into both vacuum and an ambient plasma are considered. When a Maxwellian beam is injected into vacuum, the surface changes to a potential much greater than the average beam energy. The dependence of this excess is examined by considering beams with water-bag types of velocity distribution functions in which no electron has a velocity V(max) above a certain value. The electric field distribution in the electron sheath near the surface is determined by the pressure distribution. Thus, the surface potential is determined not only by V(max) but by all the beam parameters. The ambient plasma reduces the charging level and causes an oscillation in the surface potential. The oscillation frequency is the electron-plasma frequency associated with the ambient plasma.

Singh, Nagendra↗

'Self-consistent' production of ion conics on return current region auroral field lines - A time-dependent, semi-kinetic model

We describe initial results from a time-dependent, semi-kinetic model of plasma outflow incorporating wave-particle interactions along current-carrying auroral field lines. Electrostatic waves are generated by the current driven ion cyclotron instability (CDICI), causing perpendicular velocity diffusion of ions plus electron heating via anomalous resistivity when and where the relative drift between electrons and ions exceeds certain critical velocities. Using the local bulk parameters we calculate these critical velocities, and so are able to self-consistently switch on and off the heating of the various particle species. Due to the dependence of these critical velocities on the bulk parameters of the species the heating effects exhibit quite complex spatial and temporal variations. A wide range of ion distribution functions are observed in these simulations, including conics with energies of a few electron volts and 'ring' distributions. The rings are seen to be a natural result of transverse heating and velocity filter effects and do not require coherent acceleration processes. We also observe the formation of a density depletion in hydrogen and enhanced oxygen densities at high altitudes.

Brown, David G.↗