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,279 records · Page 71

Land Use in Computable General Equilibrium Models

Constant Elasticity of Transformation (CET) functions are widely used to allocate land across uses in Computable General Equilibrium (CGE) models. These models fail to maintain area of land in balance. Several ways have been used or offered to fix this problem. This paper first examines the extent to which the CET approach may generate imbalances in area of land and identifies determinants of such imbalances. It then evaluates the merits, properties, and welfare impacts of the available approaches that restore area of land in balance. We show that heterogeneity in land prices (land rents) is the main source of imbalance in area of land, not the curvature of the CET function. Furthermore, we show that the available approaches that restore area of land in balance either introduce ad hoc adjustments in land allocation or undermine the conventional welfare assessments of the CET results and ignore the opportunity costs of land transformation. The CGE models which applied stochastic productivity distribution functions to allocate land among uses maintain area of land in balance, but they suffer from important limitations and ignore costs of land transformation.

Taheripour, Farzad↗

The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: exploring the halo occupation distribution model for emission line galaxies

We study the modelling of the halo occupation distribution (HOD) for the eBOSS DR16 emission line galaxies (ELGs). Motivated by previous theoretical and observational studies, we consider different physical effects that can change how ELGs populate haloes. We explore the shape of the average HOD, the fraction of satellite galaxies, their probability distribution function (PDF), and their density and velocity profiles. Our baseline HOD shape was fitted to a semi-analytical model of galaxy formation and evolution, with a decaying occupation of central ELGs at high halo masses. We consider Poisson and sub/super-Poissonian PDFs for satellite assignment. We model both Navarro–Frenk–White and particle profiles for satellite positions, also allowing for decreased concentrations. We model velocities with the virial theorem and particle velocity distributions. Additionally, we introduce a velocity bias and a net infall velocity. We study how these choices impact the clustering statistics while keeping the number density and bias fixed to that from eBOSS ELGs. The projected correlation function, w p , captures most of the effects from the PDF and satellites profile. The quadrupole, ξ 2 , captures most of the effects coming from the velocity profile. We find that the impact of the mean HOD shape is subdominant relative to the rest of choices. Here, we fit the clustering of the eBOSS DR16 ELG data under different combinations of the above assumptions. The catalogues presented here have been analysed in companion papers, showing that eBOSS RSD+BAO measurements are insensitive to the details of galaxy physics considered here. These catalogues are made publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurements of Beam Spin Asymmetries in p+p0 and p´p0 Dihadron Production at CLAS12

Semi-Inclusive Deep Inelastic Scattering (SIDIS) is a powerful experimental tool for studying the internal structure and dynamics of the proton, revealing how quarks and gluons are distributed and interact within it. SIDIS describes a process where an elec tron scatters off one of the constituent quarks within the proton, causing it to undergo hadronization, creating multiple hadrons in the final state. Through factorization, the full process can be split into probabilistic components: one which describes the internal structure of the proton using Parton Distribution Functions (PDFs), and another which describes the hadronization process using Fragmentation Functions (FFs). These functions are non-perturbative quantities of Quantum Chromodynamics (QCD), meaning they cannot be calculated directly from first principles and must instead be extracted from experimental measurements. Acommon approach for accessing PDFs and FFs using SIDIS is to measure asymmetries. In this context, asymmetries correspond to subtle differences in the angular distribution of outgoing particles that arise when the spin orientation of the incoming beam or target is reversed. Because many of these effects only appear when spin is involved, they isolate specific, nuanced properties of the proton’s spin-structure that are otherwise hidden in spin averaged measurements. In practice, they show up as specific azimuthal modulations (e.g., sin ¿R, sin(¿h ´ ¿R)), whose amplitudes isolate convolutions of PDFs and FFs at leading and subleading twist. Non-zero asymmetries of these angular distributions can be traced back to unique combinations of PDFs and FFs, offering a way to probe them directly. In this work, we measure SIDIS by analyzing high energy electron-proton scattering events using the CLAS12 detector at Jefferson Lab. This study focuses on subset of SIDIS referred to as dihadron SIDIS, where pairs of hadrons — here p+p0 and p´p0 — are observed. We analyzed these dihadrons using detector data collected during Fall 2018 and Spring 2019, where longitudinally polarized electrons from the CEBAF accelerator were incident on a liquid hydrogen target. A photon classifier using a Gradient Boosted Trees (GBTs) architecture was trained using Monte Carlo simulations to reduce the amount of iv false combinatorial background p0’s. When deployed on experimental data, the model in creases our dihadron statistics by up to five-fold compared to previous CLAS12 p0 analyses. This work reports the first measurements of beam spin asymmetries for p+p0 and p´p0 dihadron production in SIDIS. The measured asymmetries offer new insights to the spin-dependent structure and dynamics within the proton, as well as the spin-dependent properties of quark fragmentation. Non-zero twist-3 sin¿R amplitudes are observed, pro viding sensitivity to the subleading twist PDF e(x). The PDF e(x) encodes quark-gluon correlations within the proton — a property that is otherwise inaccessible at leading twist. Additionally, this work measured significant twist-2 modulations carried by sin(¿h ´ ¿R) and sin(2¿h ´2¿R), providing experimental access to the helicity dihadron fragmentation function (DiFF) GK 1 . Because there is no equivalent quark helicity-dependent FF in single pion SIDIS, the DiFF GK 1 offers a unique lens into novel spin-dependent fragmentation. For instance, the twist-2 modulations observed in this study are enhanced by vector mesons created during fragmentation — a behavior predicted by phenomenological models. This study broadens our understanding of dihadron fragmentation, revealing new details about the flavor and charge dependence of hadronization.

Matousek, Gregory [Duke Univ., Durham, NC (United ↗

Gaussian Mixture Model Solvers for the Boltzmann Equation

This report documents our experience constructing a numerical method for the collisional Boltzmann equation that is capable of accurately capturing the collisionless through strongly collisional limits. We explore three different functional representations and present a detailed account of a numerical method based on a spatially dependent Gaussian mixture model (GMM). The Kullback-Leibler divergence is used as a closeness measure and various expectation maximization (EM) solution algorithms are implemented to find a compact representation in velocity space for distribution functions that exhibit significant non-Maxwellian character. We discuss issues that appear with this representation over a range of Knudsen numbers for a prototypical test problem and demonstrate that the strongly collisional limit recovers a solution to Euler's equations. Looking forward, this approach is broadly applicable to the non-relativistic and relativistic collisional Vlasov equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kinetic hydromagnetic instabilities due to a spherical shell distribution of pickup ions

The stability of a spherical shell ion distribution function against low-frequency hydromagnetic waves is discussed by relaxing the assumption of parallel propagation which was imposed in the preceding discussions published in the literature. It is found that Alfven waves as well as magnetosonic waves with oblique angles of propagation can be excited by a shell distribution. In the present analysis, the emphasis is placed on the case in which the thickness of the spherical shell is small.

Wu, C. S.↗

Anisotropic jet broadening and jet shape

In this paper, we explore the use of jet substructure as a way of probing phenomena which break the isotropic behavior of jets, such as jet propagation through an anisotropically flowing quark-gluon plasma or spin correlations. We introduce two novel observables for this purpose: the azimuthal-dependent jet broadening and the azimuthal-dependent jet shape, which generalize the traditional isotropic substructure studies. Using Soft-Collinear Effective Theory, we explicitly calculate the jet functions associated with these observables with a standard jet axis and with a Winner-Take-All jet axis in both the resummed and fixed order limits. While our analysis first and foremost establishes the formalism for the azimuthal-dependent jet substructure, it also brings to light new results for jet substructure in the azimuthally integrated case, such as the semi-inclusive jet function and the exclusive jet shape for the Winner-Take-All axis, and the jet broadening in the fixed order region. As an illustrative example for the new formalism we demonstrate that the azimuthal-dependent jet broadening can be used as a direct probe of the transversity parton distribution function in deep inelastic scattering.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

White Dwarf Binaries across the H-R Diagram

We created the APOGEE-GALEX-Gaia catalog to study white dwarf (WD) binaries. This database aims to create a minimally biased sample of WD binary systems identified from a combination of GALEX, Gaia, and APOGEE data to increase the number of WD binaries with orbital parameters and chemical compositions. We identify 3414 sources as WD binary candidates, with nondegenerate companions of spectral types between F and M, including main-sequence stars, main-sequence binaries, subgiants, sub-subgiants, red giants, and red clump stars. Among our findings are (a) a total of 1806 systems having inferred WD radii R < 25 R ⊕ , which constitute a more reliable group of WD binary candidates within the main sample; (b) a difference in the metallicity distribution function between WD binary candidates and the control sample of most luminous giants (M H < -3.0); (c) the existence of a population of sub-subgiants with WD companions; (d) evidence for shorter periods in binaries that contain WDs compared to those that do not, as shown by the cumulative distributions of APOGEE radial velocity shifts; (e) evidence for systemic orbital evolution in a sample of 252 WD binaries with orbital periods, based on differences in the period distribution between systems with red clump, main-sequence binary, and sub-subgiant companions and systems with main-sequence or red giant companions; and (f) evidence for chemical enrichment during common envelope (CE) evolution, shown by lower metallicities in wide WD binary candidates (P > 100 days) compared to post-CE (P < 100 days) WD binary candidates.

79 ASTRONOMY AND ASTROPHYSICS↗

Binder model system to be used for determination of prepolymer functionality

Development of a method for determining the functionality distribution of prepolymers used for rocket binders is discussed. Research has been concerned with accurately determining the gel point of a model polyester system containing a single trifunctional crosslinker, and the application of these methods to more complicated model systems containing a second trifunctional crosslinker, monofunctional ingredients, or a higher functionality crosslinker. Correlations of observed with theoretical gel points for these systems would allow the methods to be applied directly to prepolymers.

Martinelli, F. J.↗

Rough-to-smooth transition of an equilibrium neutral constant stress layer

Purpose of research on rough-to-smooth transition of an equilibrium neutral constant stress layer is to develop a model for low-level atmospheric flow over terrains of abruptly changing roughness, such as those occurring near the windward end of a landing strip, and to use the model to derive functions which define the extent of the region affected by the roughness change and allow adequate prediction of wind and shear stress profiles at all points within the region. A model consisting of two bounding logarithmic layers and an intermediate velocity defect layer is assumed, and dimensionless velocity and stress distribution functions which meet all boundary and matching conditions are hypothesized. The functions are used in an asymptotic form of the equation of motion to derive a relation which governs the growth of the internal boundary layer. The growth relation is used to predict variation of surface shear stress.

Logan, E., Jr.↗

A kinetic cyclotron maser instability associated with a hollow beam of electrons

A kinetic cyclotron maser instability associated with a hollow-beam distribution function is studied. The instability differs from that discussed for the gyrotron device in two respects: in the present case the momentum dispersion is substantial, and furthermore there exists a low-energy background plasma. On the basis of physical arguments it can be demonstrated that the hollow-beam distribution is far more unstable than the loss cone distribution which has been extensively investigated in recent years. A criterion for maximum growth rate is established on the basis of consideration of the resonance ellipse. The validity of this criterion is supported by the results of numerical calculation.

Wu, C. S.↗

Contrasts among bidirectional reflectance of leaves, canopies, and soils

Simple models are presented for predicting the bidirectional reflectance distribution functions (BRDFs) for soils and plant canopies viewed from various directions. BRDFs are predicted for bare soil, individual leaves, and plant canopies, and the results are compared with measurements and a three coefficient empirical equation. BRDF measurements for corn and soybean leaves are presented to contrast with canopy and soil distributions. Estimates of the soil, canopy, and leaf BRDFs are combined into a model called Cupid to predict BRDFs for complex natural surfaces.

Norman, J. M.↗

Propagation and damping of broadband upstream whistlers

Previous studies indicated that damping rates of upstream whistlers strongly depend on the details of the electron distribution function. Moreover, detailed analysis of Doppler-shift and whistler dispersion relation indicated that upstream whistlers propagate obliquely in a broad band. In this paper we present results of a kinetic calculation of damping lengths of wide-band whistlers using the sum of 7-drifting bi-Maxwellian electron distributions as a best fit to the International Sun Earth Explorer (ISEE) 1 electron data. For 2 cases, when upstream whistlers are observed, convective damping lengths derived from ISEE magnetic field and ephemeris data are compared with theoretical results. We find that the calculated convective damping lengths are consistent with the data and that upstream whistlers remain marginally stable. We also show that the slope of plasma frame spectra of upstream whistlers, obtained by direct fitting of the observed spectra is between 5 and 7 with a sharp lower frequency cutoff corresponding to a wavelength of about one ion inertial length. When the solar wind velocity is directed largely along the wave normal of the upstream whistlers the polariztion of the right hand waves becomes reversed and low frequencies are switched to high resulting in a peaked spectrum with a strong high frequency cutoff. The overall spectral, wave and particle characteristics, proximity to the shock as well as propagation and damping properties indicate that these waves cannot be generated locally. Instead the observed upstream whistlers arise in the shock ramp most likely by a variety of cross-field drift and/or anisotropy driven instabilities.

Orlowski, D. S.↗

Consequences of the large ambipolar electric field in the solar wind

The parallel electric field in the solar wind is much smaller than the V x B motional electric field, yet in the proper dimensions units it is very 'large'. At the orbit of earth it is within a few percent of being at the Dreicer limit. This 'large' electric field is required for quasi-neutrality; it will be shown to have interesting consequences for the electron velocity distribution function and the description of transport of heat. Interestingly, a similar dimensionless situation also occurs at the base of the transition region, while below the transition region the dimensionless electric field is very weak. These facts suggest a new way to look at the thermal-suprathermal dichotomy in velocity distributions as the response of a plasma where charge neutrality requires such large dimensionless electric fields.

Scudder, J. D.↗

Origin and relative age of lunar and Mercurian intercrater plains

Differences in the diameter/density distribution of craters in the size range 7-100 km between heavily cratered areas and areas of intercrater plains indicate the deposition of both lunar and Mercurian intercrater plains has preferentially obliterated craters not greater than about 30 km diameter. This obliteration has altered the original lunar population distribution function from about a -1.3 slope to about -0.7 slope. At least a significant portion of both the lunar and Mercurian intercrater plains appear to be volcanic deposits emplaced during the later stages of heavy bombardment on each body rather than basin ejecta or an ancient primordial surface. The very widespread distribution of intercrater plains on Mercury compared to the moon may be related to Mercury's core formation which would have resulted in a large radius increase leading to widespread extensional fracturing to provide egress for the eruption of volcanic deposits on a global scale. This was followed by a radius decrease due to contraction of the lithosphere to produce thrust faulting represented by Mercury's lobate scarps.

Strom, R. G.↗

Simulation of spacecraft charging environments by monoenergetic beams

To simulate the distributed spectra of plasmas in space, the differential current density spectrum of the plasma was divided into a number of energy bands and the beam energy and current were calculated for each band to provide a piecewise reproduction of the distributed spectrum. Beam energies and current densities were chosen to match the velocity moments of the plasma distribution function. The velocity moments are averages related to physical quantities such as particle density, flux, pressure, and energy flux, and have been used extensively to characterize the measured properties of plasmas in space. Combinations of one, two, and three beams were found to match two to six velocity moments of Maxwellian distributions. A computational model was used to compare the charging of a spacecraft by plasmas with distributed spectra and by monoenergetic beams.

Halverson, W.↗

Magnetic field line reconnection experiments. VI - Magnetic turbulence

Extensive statistical analysis of the vector components of fluctuating magnetic fields have been performed in a time dependent neutral magnetic sheet. Cross spectral analysis indicates a variety of wave numbers present for each frequency investigated. Comparison of Fourier components of the cross spectral function with dispersion surfaces in k space demonstrates the waves are large amplitude whistlers; this is verified by polarization analysis, which shows the random waves' magnetic fields to be right-hand circular. Ion acoustic and Langmuir turbulence are also observed along with bursts of microwave radiation. Measurements of the electron distribution function f(v, r, t) and its fluctuations in velocity space relate wave and particle activity.

Gekelman, W.↗

The cause of the coronal temperature inversion of the solar atmosphere and the implications for the solar wind

The energy contained in suprathermal tails at the base of the transition region is shown to be transformed into the rarefied, but hotter, transition region and low corona without any further addition of energy to the gas above the base of the transition region. Possible critical point location and asymptotic wind speed are shown to be controlled by the suprathermal tail strength parameter used to model possible suprathermal velocity distribution functions at the base of the transition region. This process shows promise for producing temperature profiles that peak near, but outside of, the fluid critical point without ad hoc energy deposition. The coronal temperature inversion above the solar photosphere is argued to be a generic feature around all stars with nonthermal distributions at the heights where the atmosphere last becomes mostly ionized.

Scudder, J. D.↗

New description of charged particle propagation in random magnetic fields

When charged particles spiral along a large constant magnetic field, their trajectories are scattered by random components that are superposed on the guiding field. In the simplest analysis of this situation, scattering causes the particles to diffuse parallel to the guiding field. At the next level of approximation, moving pulses that correspond to a coherent mode of propagation are present, but they are represented by delta-functions whose infinitely narrow width makes no sense physically and is inconsistent with the finite duration of coherent pulses observed in solar energetic particle events. To derive a more realistic description, the transport problem is formulated in terms of 4 x 4 matrices, which derive from a representation of the particle distribution function in terms of eigenfunctions of the scattering operator, and which lead to useful approximations that give explicit predictions of the detailed evolution not only of the coherent pulses, but also of the diffusive wake. More specifically, the new description embodies a simple convolution of a narrow Gaussian with the solutions above that involve delta-functions, but with a slightly reduced coherent velocity. The validity of these approximations, which can easily be calculated on a desktop computer, has been exhaustively confirmed by comparison with results of Monte Carlo simulations which kept track of 50 million particles and which were carried out on the Maspar computer at Goddard Space Flight Center.

Earl, James A.↗