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,261 records · Page 70

Numerical Investigation of Enhanced Dehumidification Processes By Using Dielectrophoresis Principles in Moist Airflows

Dispersed particle-laden flows are encountered in many building and industrial applications, such as flow in a fluidized bed, hydrocarbon transportation in pipelines, and the fouling of air-cooled heat exchangers (Kuruneru et al., 2016; Ray et al., 2019; Wang et al., 2019). Computational fluid dynamic (CFD) models have been developed in recent years to depict particle-fluid and particle-particle interactions in laminar or turbulent flows with increasing accuracy and stability. One particular particle-laden system of interest for moisture control is electrically-enhanced condensation in air and water droplet flows. Electrically-enhanced condensation consists of the use of highly charged water droplets injected in the moist air. The droplets become electric seeds that attract polar water vapor molecules to their surfaces and promote condensation. The nucleation and growth of the charged droplets deplete the vapor phase near a droplet, which is compensated for by the dielectrophoresis flow and diffusion. Dielectrophoresis flow involves surrounding vapor at a distance of about 10 to 100 nm for droplets charged by an electrospray compared to ~2 nm for a single electron charge in a droplet. As the vapor molecules collapse on the surface of the droplets, their initial electrical charge decreases with time due to the neutralization of the ions. While the physics of this phenomena is well known, engineering models for predicting the condensation rates are not available. This work computationally investigates dehumidification of moist airflow in a converging rectangular duct. The objective is to develop an engineering model that predicts water vapor condensation by employing dielectrophoresis principles. We construct a Computational Fluid Dynamics (CFD) model of the duct with electrically-enhanced condensation. The model is implemented in the open-source software OpenFOAM. We utilize the Multi-Phase Particle-In-Cell (MP-PIC) method coupled with a Population Balance Equation (PBE) approach to simulate the particle-laden system. This methodology is an Eulerian-Lagrangian approach used to simulate the droplets' behavior in the humid air. The MP-PIC approach (Andrews and O'Rourke, 1996) mitigates the computational cost by parceling several fundamental particles with similar properties (such as types, sizes, and temperature) into one computational particle. Thus, the billions of particles can be substituted by millions of computational particles without significant loss of information. The PBE was considered with the Lagrangian frame to combine the particle distribution function used in MP-PIC (Kim et al., 2020). This approach preserves mass and energy conservation between the phases in the Eulerian and Lagrangian structures. The PBE in this procedure was directly linked to the discrete parcels, making the simulation of the particle distribution computationally efficient and robust. The MP-PIC-PBE approach used in the present work was applied to the dehumidification of air. Water droplets were injected in the air stream and forced to grow according to experimentally derived correlation. The experiments were conducted on a converging duct with the same geometry and boundary conditions used to build the CFD model. This approach enabled us to approximate the effect of dielectrophoresis phenomena on the droplet and air interface. This presentation will discuss the details of the new CFD model built for the duct, the implementation of the model in OpenFOAM CFD programming language, and the experimental validation of the newly developed model. The results revealed a moderate yet measurable increase in droplet diameter due to water vapor condensation at the vapor-liquid interface of the electrically charged droplets' surface. The seed water droplet particles grew in size by capturing the water vapor in the surrounding air. The OpenFOAM model predicted reductions of humidity in the air from 5 to 10 percent.

Yel Mahi, Maliha↗

Theory of fluctuations of extended air showers

The Monte Carlo method was used to calculate the probability distribution functions of shower characteristics for primary protons at sea level. The calculation was based on the following model of the elementary event: the interaction paths are 90 g/sq cm for nucleons and 120 g/sq cm for pions. The nonelasticity coefficient for nucleons is uniformly distributed between 0.1 and 0.9, and for pions it is equal to 1. Isobaric pions are taken into account. The spectra of secondary particles were determined using Cocconi's approximation formula. The calculation for the nuclei was carried out on the assumption of a breakup of the nucleus into component nucleons. The mean number of particles and the variances of the distributions for electrons when the number of muons was fixed, and for muons when the number of electrons was fixed, were calculated.

Dedenko, L. G.↗

Acoustical theory of turbulence

Fluid functions are expanded in a series of functions constituting a complete and orthogonal system of wave solutions of the linearized Navier-Stokes system. The motion of characteristic waves is described by Hamiltonian equations of ray acoustics. Amplitude of the waves is described with the aid of nonlinear interaction terms. The distribution functions (squares of amplitudes) satisfy Boltzmann-type equations and completely describe the mean turbulence properties.

Kentzer, C. P.↗

Theoretical study of hole initiated impact ionization in bulk silicon and GaAs using a wave-vector-dependent numerical transition rate formulation within an ensemble Monte Carlo calculation

In this paper, calculations of the hole initiated interband impact ionization rate in bulk silicon and GaAs are presented based on an ensemble Monte Carlo simulation with the inclusion of a wave-vector-dependent numerical transition rate formulation. The ionization transition rate is determined for each of the three valence bands, heavy, light, and split-off, using Fermi's golden rule with a two-body, screened Coulomb interaction. The dielectric function used within the calculation is assumed to be wave-vector-dependent. Calculations of the field-dependent impact ionization rate as well as the quantum yield are presented. It is found from both the quantum yield results and examination of the hole distribution function that the effective threshold energy for hole initiated impact ionization is relatively soft, similar to that predicted for the corresponding electron initiated ionization events occur more frequently than either heavy or split-offf initiated ionization events in bulk silicon over the applied electric field strengths examined here, 250-500 kV/cm. Conversely,in GaAs, the vast majority of hole initated ionization events originate from holes within the split-off band.

Oguzman, Ismail H.↗

Plasma Dispersion Function for the Kappa Distribution

The plasma dispersion function is computed for a homogeneous isotropic plasma in which the particle velocities are distributed according to a Kappa distribution. An ordinary differential equation is derived for the plasma dispersion function and it is shown that the solution can be written in terms of Gauss' hypergeometric function. Using the extensive theory of the hypergeometric function, various mathematical properties of the plasma dispersion function are derived including symmetry relations, series expansions, integral representations, and closed form expressions for integer and half-integer values of K.

Podesta, John J.↗

Unpolarized transverse momentum dependent parton distributions of the nucleon from lattice QCD

We present a first lattice QCD calculation of the unpolarized nucleon’s isovector transverse-momentum-dependent parton distribution functions (TMDPDFs), which are essential to predict observables of multiscale, semi-inclusive processes in the standard model. We use a N f = 2 + 1 + 1 MILC ensemble with valence clover fermions on a highly improved staggered quark (HISQ) sea to compute the quark momentum distributions in a large-momentum nucleon on the lattice. The state-of-the-art techniques in renormalization and extrapolation in the correlation distance on the lattice are adopted. The perturbative kernel up to next-to-next-to-leading order is taken into account, and the dependence on the pion mass and the hadron momentum is explored. Our results are qualitatively comparable with phenomenological TMDPDFs, which provide an opportunity to predict high energy scatterings from first principles. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Continued development of a detailed model of arc discharge dynamics

Using a previously developed set of codes (SEMC, CASCAD, ACORN), a parametric study was performed to quantify the parameters which describe the development of a single electron indicated avalanche into a negative tip streamer. The electron distribution function in Teflon is presented for values of the electric field in the range of four-hundred million volts/meter to four billon volts/meter. A formulation of the scattering parameters is developed which shows that the transport can be represented by three independent variables. The distribution of ionization sites is used to indicate an avalanche. The self consistent evolution of the avalanche is computed over the parameter range of scattering set.

Beers, B. L.↗

Theory and application of the vector pair correlation function for real-space crystallographic analysis of order/disorder correlations from STEM images

Deviations of local structure and chemistry from the average crystalline unit cell are increasingly recognized to have a significant influence on the properties of many technologically important materials. Here, we present the vector pair correlation function (vPCF) as a new real-space crystallographic analysis method, which can be applied to atomic-resolution scanning transmission electron microscopy (STEM) images to quantify and analyze structural order/disorder correlations. Our STEM-based vPCFs have several advantages over radial PCFs and/or 3D pair distribution functions from x-ray total scattering: vPCFs explicitly retain crystallographic orientation information, are spatially resolved, can be applied directly on a sublattice basis, and are suitable for any material that can be imaged with STEM. To show the utility of our approach, we measure partial vPCFs in Ba 5 SmSn 3 Nb 7 O 30 (BSSN), a tetragonal tungsten bronze (TTB) structured complex oxide. Many TTBs are known to be classical or relaxor ferroelectrics, and these properties have been correlated with the presence of superlattice ordering. BSSN, specifically, exhibits relaxor behavior and an incommensurate structural modulation. From the vPCF data, we observe that, of the cation sites, only the Ba (A2) sublattice is structurally modulated. We then infer the local modulation vector and reveal a marked anisotropy in its correlation length. Finally, short-range correlated polar displacements on the B2 cation sites are observed. This work introduces the vPCF as a powerful real-space crystallography technique, which enables direct, robust quantification of short-to-long range order on a sublattice-specific basis and is applicable to a wide range of complex material types.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

General mass variable flavor number scheme for 𝑍 boson production in association with a heavy quark at hadron colliders

We present a methodology to streamline implementation of massive-quark radiative contributions in calculations with a variable number of active partons in proton-proton collisions. The methodology introduces subtraction and residual heavy-quark parton distribution functions (PDFs) to implement calculations in the Aivazis–Collins–Olness–Tung (ACOT) factorization scheme and its simplified realization in various processes up to the next-to-the-next-to-leading order in the QCD coupling strength. Interpolation tables for bottom-quark subtraction and residual distributions for CT18 NLO and NNLO PDF ensembles are provided in the common LHAPDF6 format. A numerical calculation of 𝑍-boson production with at least one 𝑏 jet at the Large Hadron Collider beyond the lowest order in QCD is considered for illustration purposes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Characteristics of suprathermal H(+) and He(++) in plasmoids in the distant magnetotail

Rest frame distribution functions were obtained for suprathermal H(+) and He(2+) ions during 20 particle intensity increases that were observed in the distant magnetotail by ISEE-3. All of these increases were characterized by isotropic, approximately 100-keV electrons, preceded by beaming electrons and velocity-dispersed protons distributions whose signatures are typically associated with plasmoids. In each of the plasmoids, both H(+) and He(2+) are found to be convected tailward, with speeds that range from 200 to 1000 km/sec. The acceleration mechanism is found to depend on particle velocity, and it is established that the solar wind is a dominant plasma source. The plasmoids appear to be expanding at about 100 km/sec during their tailward propagation.

Gloeckler, G.↗

Experimenter's data package for the descending layers rocket

In response to a proposal from Science Applications International Corporation (SAIC), NASA Headquarters has approved a sounding rocket mission designed to study the physics of intermediate layers in the Earth's ionosphere at middle latitudes. The experiment will be carried out by a team of scientists and engineers from the NASA Wallops Flight Facility, SAIC, the NASA Goddard Space Flight Center, and the Millstone Hill radar observatory. The mission will involve the launch of an instrumented sounding rocket from the Wallops Island rocket range in the summer of 1994, with the objective of penetrating a descending ionized layer in the E-region between altitudes of 115 and 140 km. Instrumentation aboard the rocket will measure the ion and neutral composition of the layer, its plasma density, driving wind and electric field forces, the thermal ion distribution function, and electron temperature. Depending on payload weight constraints and subject to availability, a particle detector to measure energetic ion and/or electron fluxes near the layer may also be included. This document was prepared as a reference for the NASA payload development and experiment teams, for distribution at the Project Initiation Conference (PIC). The design specifications discussed herein are therefore of a preliminary nature; the intent is to promote open discussions between experimenters and NASA engineers that will lead to a final design capable of achieving the experiment objectives.

Earle, Greg↗

Conditional H I Mass Functions and the H I -to-halo Mass Relation in the Local Universe

We present a new H I mass estimator that relates log 10 (M HI /M * ) to a linear combination of four galaxy properties: stellar surface mass density, color index u - r, stellar mass, and concentration index, with the scatter of individual galaxies around the mean H I mass modeled with a Gaussian distribution function. We calibrate the estimator using the xGASS sample, including both H I detection and nondetection, and constrain the model parameters through Bayesian inferences. Tests with mock catalogs demonstrate that our estimator provides unbiased H I masses for optical samples like SDSS. We apply our estimator to the SDSS spectroscopic sample to estimate the H I mass function (HIMF) of local galaxies, as well as the conditional H I mass function in galaxy groups and the H I –halo mass relation. Our HIMF agrees with the ALFALFA measurements at M HI ≳ 5 × 10 9 M ⊙ , but with higher amplitude and a steeper slope at lower masses. We show that this discrepancy is caused primarily by the cosmic variance, which is corrected for the SDSS sample but not for ALFALFA. The total CHIMFs for all halo masses can be described by a single Schechter function, while those of central galaxies show a double-Gaussian profile. The total H I mass in a group increases monotonically with halo mass, but for central galaxies, the H I mass shows weak dependence on halo mass when M h ≳ 10 12 M ⊙ . The observed H I –halo mass relation is not reproduced by current hydrodynamic simulations and semianalytic models of galaxy formation.

79 ASTRONOMY AND ASTROPHYSICS↗

Novel data interpretation method for DIII-D divertor retarding field energy analyzer with 3-D particle-in-cell simulations

A novel data interpretation process that utilizes comprehensive particle-in-cell (PIC) simulations is developed for the new retarding field energy analyzer (RFEA) currently being constructed at DIII-D for the lower divertor using the Divertor Material Evaluation System. Furthermore, this probe is expected to survive a heat load of up to 100 MW/m 2 for up to 5 s and reliably measure the main ion temperature (T i ) on the divertor target ranging from 10 to 200 eV. These extreme conditions posed significant engineering limitations on the probe geometry, thus extensive validation work has been performed. The conventional fitting method for the RFEA I–V characteristics is based on a simplified 1-D model without considering the ion space charge inside the probe cavity and may not be sufficient for probes designed for the DIII-D divertor environment. In this article, a more realistic description of the particle propagation process within the RFEA cavity is achieved by including both 3-D geometric effects and ion space charge in the PIC simulations, and the capability to reconstruct the ion energy distribution functions is demonstrated with reasonable consistency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Structure of Soot-Containing Laminar Jet Diffusion Flames

The structure and soot properties of nonbuoyant and weakly-buoyant round jet diffusion flames were studied, considering ethylene, propane and acetylene burning in air at pressures of 0.125-2.0 atm. Measurements of flame structure included radiative heat loss fractions, flame shape and temperature distributions in the fuel-lean (overfire) region. These measurements were used to evaluate flame structure predictions based on the conserved-scalar formalism in conjunction with the laminar flamelet concept, finding good agreement betweem predictions and measurements. Soot property measurements included laminar smoke points, soot volume function distributions using laser extinction, and soot structure using thermophoretic sampling and analysis by transmission electron microscopy. Nonbuoyant flames were found to exhibit laminar smoke points like buoyant flames but their properties are very different; in particular, nonbuoyant flames have laminar smoke point flame lengths and residence times that are shorter and longer, respectively, than buoyant flames.

Mortazavi, S.↗

LRAT: Lightning Radiative Transfer

In this report, we extend to cloud physics the work done for single and multiple scattering of electromagnetic waves. We consider the scattering of light, visible or infrared, by a spherical cloud represented by a statistically homogeneous ensemble of configurations of N identical spherical water droplets whose centers are uniformly distributed in its volume V. The ensemble is specified by the average number rho of scatterers in unit volume and by rho f(R) with f(R) as the distribution function for separations R of pairs. The incident light, vector-phi(sub 0) a plane electromagnetic wave with harmonic time dependence, is from outside the cloud. The propagation parameter kappa(sub 0) and the index of refraction eta(sub 0) determine physically the medium outside the distribution of scatterers. We solve the interior problem separately to obtain the bulk parameters for the scatterer equivalent to the ensemble of spherical droplets. With the interior solution or the equivalent medium approach, the multiple scattering problem is reduced to that of an equivalent single scatterer excited from outside illumination. A dispersion relation which determines the bulk propagation parameter K and the bulk index of refraction eta of the cloud is given in terms of the vector equivalent scattering amplitude vector-G and the dyadic scattering amplitude tilde-g of the single object in isolation. Based on this transfer model we will have the ability to consider clouds composed of inhomogeneous distribution of water and/or ice particles and we will be able to take into account particle size distributions within the cloud. We will also be able to study the effects of cloud composition (i.e., particle shape, size, composition, orientation, location) on the polarization of the single or the multiple scattered waves. Finally, this study will provide a new starting point for studying the problem of lightning radiative transfer.

Phanord, Dieudonne D.↗

Processing and mechanical characterization of short carbon fiber-reinforced epoxy composites for material extrusion additive manufacturing

Fiber-reinforced polymer composites have been extensively utilized in recent years as feedstock materials for material extrusion additive manufacturing (AM) processes to improve strength, stiffness, and functionality of printed parts over unfilled printed polymers. However, the widespread adoption of AM of fiber-reinforced polymer composites requires a deeper understanding of the process-structure-property relationships in printed components, and such relationships are not well understood yet. Fiber length is critically important to the mechanical performance of short fiber composites, but very few studies to-date have focused on how the fiber length distribution (FLD) evolves during processing of composite feedstocks and how this evolution affects printing behavior and mechanical properties in 3D-printed composites. Here, FLD is measured for carbon fiber reinforced epoxy composites over a wide range of ink compositions and shear mixing times, and the distributions are fit with a Weibull-type distribution function. The effects of FLD on the tradeoff between ink processability, ink rheology, printing behavior and mechanical properties are investigated. Furthermore, the effects of printing parameters (nozzle size and print speed) on mechanical anisotropy and fiber orientation distribution (FOD) in printed composites are explored. Mechanical properties of printed composites are characterized via 3 pt-flexural testing, and microstructure is investigated using optical and scanning electron microscopy (SEM), and x-ray computed tomography. Finally, the fitted Weibull parameters are fed into a composite model that incorporates FLD and FOD, and model predictions are found to be in excellent agreement with experimental observations.

3D printing↗

Evolution of in-medium baryon-baryon scattering cross sections and stiffness of dense nuclear matter from Bayesian analyses of FOPI proton-flow excitation functions

Within a Bayesian statistical framework using a Gaussian process (GP) emulator for an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model simulator of heavy-ion reactions with momentum-independent Skyrme interactions, we infer from the proton directed and elliptical flow in mid-central Au + Au reactions, measured at beam energies from 150 to 1200 MeV/nucleon by the FOPI Collaboration, the posterior probability distribution functions (PDFs) of the in-medium baryon-baryon scattering cross section (BBSCS) modification factor 𝑋 (with respect to their free-space values) and the stiffness parameter 𝐾 of dense nuclear matter. Here, we find that the most probable value of 𝑋 evolves from around 0.7 to 1.0 as the beam energy 𝐸 beam /𝐴 increases. On the other hand, the posterior PDF⁡(𝐾) may have dual peaks having roughly the same height or extended shoulders at high 𝐾 values. More quantitatively, the posterior PDF⁡(𝐾) changes from having a major peak around 220 MeV, characterizing a soft EOS in the reaction at 𝐸 beam /𝐴 = 150 MeV, to one that peaks around 320 MeV, indicating a stiff EOS in the reactions at 𝐸 beam /𝐴 higher than about 600 MeV. The transition from soft to stiff happens in mid-central Au + Au reactions at beam energies around 250 MeV/nucleon in which 𝐾 = 220 MeV and 𝐾 = 320 MeV are approximately equally probable. Altogether, the FOPI proton-flow excitation function data indicate a gradual hardening of hot and dense nuclear matter as its density and temperature increase in reactions with higher beam energies.

nuclear matter↗

Effects of multiply charged ions on microturbulence-driven electron transport in partially magnetized plasmas

Herein, nonlinear interaction between kinetic instabilities in partially magnetized plasmas in the presence of multiply-charged ion streams is investigated using kinetic simulations. It was observed by Hara and Tsikata [Phys. Rev. E 102, 023202 (2020)] that the axial ion-ion two-stream instability (IITSI) due to singly and doubly charged ion streams, coupled with the azimuthal electron cyclotron drift instability (ECDI), enhances cross-field electron transport. In the present study, it is observed that the addition of triply charged ions (as a third ion species) contributes to damping of the excited modes, leading to a reduction in the cross-field electron transport. The net instability-driven electron transport is shown to be a function not only of the azimuthal modes, such as the ECDI, but of the multiple ion species that dictate the development of additional plasma waves. It is found that trapping of the higher ion charge states within the plasma waves results in a broadening of the ion velocity distribution functions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗