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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

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↗
DOE OSTI2022

Efficient reconstruction and validation of heterogeneous microstructures for energy applications

The digital reconstruction of microstructures is necessary for simulations in fields ranging from geology to electrochemistry, but the state-of-the-art digital reconstruction techniques often compromise between resolution and field of view. It is challenging to retain detailed microstructure information in large-scale reconstructions. Here, this study investigates different aspects of the Yeong-Torquato algorithm based on correlation functions to make it more efficient. We achieve this goal by reducing the computational complexity of the chord-length distribution function and the two-point correlation function, applying the random sphere-packing method as the initial condition, and restricting potential voxel swaps to interfaces. In addition, a novel superposition parallel scheme is introduced to aid in searching for potential voxel swaps. The algorithm proposed is validated by comparing the pore-size distributions of reconstructed 3D custom battery electrodes from a sample dataset obtained from transmission X-ray microscopy. From a sample image with 200 x 200 pixels, the code can reconstruct a 300 x 300 x 300 structure in under 22 h and reconstruct a 400 x 400 x 400 structure in 43 h with eight cores.

42 ENGINEERING↗
DOE OSTI2021

Photometric redshifts for the S-PLUS Survey: Is machine learning up to the task?

The Southern Photometric Local Universe Survey (S-PLUS) is a novel project that aims to map the Southern Hemisphere using a twelve filter system, comprising five broad-band SDSS-like filters and seven narrow-band filters optimized for important stellar features in the local universe. In this paper we use the photometry and morphological information from the first S-PLUS data release (S-PLUS DR1) cross-matched to unWISE data and spectroscopic redshifts from Sloan Digital Sky Survey DR15. We explore three different machine learning methods (Gaussian Processes with GPz and two Deep Learning models made with TensorFlow) and compare them with the currently used template-fitting method in the S-PLUS DR1 to address whether machine learning methods can take advantage of the twelve filter system for photometric redshift prediction. Using tests for accuracy for both single-point estimates such as the calculation of the scatter, bias, and outlier fraction, and probability distribution functions (PDFs) such as the Probability Integral Transform (PIT), the Continuous Ranked Probability Score (CRPS) and the Odds distribution, we conclude that a deep-learning method using a combination of a Bayesian Neural Network and a Mixture Density Network offers the most accurate photometric redshifts for the current test sample. In conclusion, it achieves single-point photometric redshifts with scatter (σ NMAD ) of 0.023, normalized bias of -0.001, and outlier fraction of 0.64% for galaxies with r_auto magnitudes between 16 and 21.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Modeling the Behavior of Complex Aqueous Electrolytes Using Machine Learning Interatomic Potentials: The Case of Sodium Sulfate

Understanding the structure and thermodynamics of solvated ions is essential for advancing applications in electrochemistry, water treatment, and energy storage. While ab initio molecular dynamics methods are highly accurate, they are limited by short accessible time and length scales whereas classical force fields struggle with accuracy. Herein, we explore the structure and thermodynamics of complex monovalent-divalent ion pairs using Na 2 SO 4 (aq) as a case study by applying a machine learning interatomic potential (MLIP) trained on density functional theory (DFT) data. Our MLIP-based approach reproduces key bulk properties such as density and radial distribution functions of water. We provide the hydration structure of the sodium and sulfate ions in the 0.1–2 M concentration range and the one-dimensional and two-dimensional potentials of mean force for the sodium–sulfate ion pairing at the low concentration limit (0.1 M), which are inaccessible to DFT. At low concentrations, the sulfate ion is strongly solvated, leading to the stabilization of solvent-separated ion pairs over contact ion pairs. Minimum energy pathway analysis revealed that coordinating two sodium ions with a sulfate ion is a multistep process whereby the sodium ions coordinate to the sulfate ion sequentially. Finally, we demonstrate that MLIPs allow the study of solvated ions beyond simple monovalent pairs with DFT-level accuracy in their low concentration limit (0.1 M) via statistically converged properties from ns-long simulations.

anions↗
DOE OSTI2022

Substructure in the stellar halo near the Sun: II. Characterisation of independent structures

In an accompanying paper, we present a data-driven method for clustering in ‘integrals of motion’ space and apply it to a large sample of nearby halo stars with 6D phase-space information. The algorithm identified a large number of clusters, many of which could tentatively be merged into larger groups. The goal here is to establish the reality of the clusters and groups through a combined study of their stellar populations (average age, metallicity, and chemical and dynamical properties) to gain more insights into the accretion history of the Milky Way. To this end, we developed a procedure that quantifies the similarity of clusters based on the Kolmogorov–Smirnov test using their metallicity distribution functions, and an isochrone fitting method to determine their average age, which is also used to compare the distribution of stars in the colour–absolute magnitude diagram. Also taking into consideration how the clusters are distributed in integrals of motion space allows us to group clusters into substructures and to compare substructures with one another. We find that the 67 clusters identified by our algorithm can be merged into 12 extended substructures and 8 small clusters that remain as such. The large substructures include the previously known Gaia-Enceladus, Helmi streams, Sequoia, and Thamnos 1 and 2. We identify a few over-densities that can be associated with the hot thick disc and host a small metal-poor population. Especially notable is the largest (by number of member stars) substructure in our sample which, although peaking at the metallicity characteristic of the thick disc, has a very well populated metal-poor component, and dynamics intermediate between the hot thick disc and the halo. We also identify additional debris in the region occupied by Sequoia with clearly distinct kinematics, likely remnants of three different accretion events with progenitors of similar masses. Although only a small subset of the stars in our sample have chemical abundance information, we are able to identify different trends of [Mg/Fe] versus [Fe/H] for the various substructures, confirming our dissection of the nearby halo. We find that at least 20% of the halo near the Sun is associated to substructures. When comparing their global properties, we note that those substructures on retrograde orbits are not only more metal-poor on average but are also older. We provide a table summarising the properties of the substructures, as well as a membership list that can be used for follow-up chemical abundance studies for example.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The Milky Way’s bulge star formation history as constrained from its bimodal chemical abundance distribution

ABSTRACT We conduct a quantitative analysis of the star formation history (SFH) of the Milky Way’s (MW) bulge by exploiting the constraining power of its stellar [Fe/H] and [Mg/Fe] distribution functions. Using Apache Point Observatory Galactic Evolution Experiment survey data, we confirm the previously established bimodal [Mg/Fe]–[Fe/H] distribution within 3 kpc of the inner Galaxy. To fit the chemical bimodal distribution, we use a simple but flexible star formation framework, which assumes two distinct stages of gas accretion and star formation, and systematically evaluate a wide multidimensional parameter space. We find that the data favour a three-phase SFH that consists of an initial starburst, followed by a rapid star formation quenching episode, and a lengthy, quiescent secular evolution phase. The metal-poor, high-α bulge stars ([Fe/H] < 0.0 and [Mg/Fe] > 0.15) are formed rapidly (<2 Gyr) during the early starburst. The density gap between the high- and low-α sequences is due to the quenching process. The metal-rich, low-α population ([Fe/H] > 0.0 and [Mg/Fe] < 0.15) then accumulates gradually through inefficient star formation during the secular phase. This is qualitatively consistent with the early SFH of the inner disc. Given this scenario, a notable fraction of young stars (age <5 Gyr) is expected to persist in the bulge. Combined with extragalactic observations, these results suggest that a rapid star formation quenching process is responsible for bimodal distributions in both the MW’s stellar populations and in the general galaxy population and thus plays a critical role in galaxy evolution.

Lian, Jianhui↗
DOE OSTI2021

Quark sivers function at small x: spin-dependent odderon and the sub-eikonal evolution

We apply the formalism developed earlier for studying transverse momentum dependent parton distribution functions (TMDs) at small Bjorken x to construct the small- x asymptotics of the quark Sivers function. First, we explicitly construct the complete fundamental “polarized Wilson line” operator to sub-sub-eikonal order: this object can be used to study a variety of quark TMDs at small x . We then express the quark Sivers function in terms of dipole scattering amplitudes containing various components of the “polarized Wilson line” and show that the dominant (eikonal) term which contributes to the quark Sivers function at small x is the spin-dependent odderon, confirming the re- cent results of Dong, Zheng and Zhou. Our conclusion is also similar to the case of the gluon Sivers function derived by Boer, Echevarria, Mulders and Zhou. We also analyze the sub-eikonal corrections to the quark Sivers function using the constructed “polarized Wilson line” operator. We derive new small- x evolution equations re-summing double-logarithmic powers of α s ln 2 (1 /x ) with α s the strong coupling constant. We solve the corresponding novel evolution equations in the large- N c limit, obtaining a sub-eikonal correction to the spin-dependent odderon contribution. We conclude that the quark Sivers function at small x receives contributions from two terms and is given by ${f}_{1T}^{\perp q}\left(x,{k}_T^2\right)={C}_O\left(x,{k}_T^2\right)\frac{1}{x}+{C}_1\left({k}_T^2\right){\left(\frac{1}{x}\right)}^0+\cdots$ with the function C O ( x, ${k}_T^2$) varying slowly with x and the ellipsis denoting the subasymptotic and sub-sub-eikonal (order- x ) corrections.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Detection of spatial clustering in the 1000 richest SDSS DR8 redMaPPer clusters with nearest neighbor distributions

ABSTRACT Distances to the k-nearest-neighbor (kNN) data points from volume-filling query points are a sensitive probe of spatial clustering. Here, we present the first application of kNN summary statistics to observational clustering measurement, using the 1000 richest redMaPPer clusters (0.1 ≤ z ≤ 0.3) from the SDSS DR8 catalog. A clustering signal is defined as a difference in the cumulative distribution functions (CDFs) of kNN distances from fixed query points to the observed clusters versus a set of unclustered random points. We find that the k = 1, 2-NN CDFs of redMaPPer deviate significantly from the randoms’ across scales of 35 to 155 Mpc, which is a robust signature of clustering. In addition to kNN, we also measure the two-point correlation function for the same set of redMaPPer clusters versus random points, which shows a noisier and less significant clustering signal within the same radial scales. Quantitatively, the χ2 distribution for both the kNN-CDFs and the two-point correlation function measured on the randoms peak at χ2 ∼ 50 (null hypothesis), whereas the kNN-CDFs (χ2 ∼ 300, p = 1.54 × 10−36) pick up a much more significant clustering signal than the two-point function (χ2 ∼ 100, p = 1.16 × 10−6) when measured on redMaPPer. Finally, the measured 3NN and 4NN CDFs deviate from the predicted k = 3, 4-NN CDFs assuming an ideal Gaussian field, indicating a non-Gaussian clustering signal for redMaPPer clusters, although its origin might not be cosmological due to observational systematics. Therefore, kNN serves as a more sensitive probe of clustering complementary to the two point correlation function, providing a novel approach for constraining cosmology and galaxy–halo connection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Residual mean field model of valence quarks in the nucleon

Abstract We develop a non-perturbative model for valence parton distribution functions (PDFs) based on the mean field interactions of valence quarks in the nucleonic interior. The main motivation for the model is to obtain a mean field description of the valence quarks as a baseline to study the short range quark–quark interactions that generate the high x tail of PDFs. The model is based on the separation of the valence three-quark cluster and residual system in the nucleon. Then the nucleon structure function is calculated within the effective light-front diagrammatic approach introducing nonperturbative light-front valence quark and residual wave functions. Within the model a new relation is obtained between the position, $$x_p$$ x p , of the peak of $$xq_V(x)$$ x q V ( x ) distribution of the valence quark and the effective mass of the residual system, $$m_R$$ m R , in the form: $$x_{p} \approx {1\over 4} (1-{m_R\over m_N})$$ x p ≈ 1 4 ( 1 - m R m N ) at starting $$Q^2$$ Q 2 . This relation explains the difference in the peak positions for d- and u-quarks through the expected difference of residual masses for valence d- and u-quark distributions. The parameters of the model are fixed by fitting the calculated valence quark distributions to the phenomenological PDFs. This allowed us to estimate the overall mean field contribution in baryonic and momentum sum rules for valence d- and u-quarks. Finally, the evaluated parameters of the non-perturbative wave functions of valence 3q-cluster and residual system can be used in calculation of other quantities such as nucleon form factors, generalized partonic and transverse momentum distributions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Highly efficient semi-hydrogenation in strained ultrathin PdCu shell and the atomic deciphering for the unlocking of activity-selectivity

Excellent ethylene selectivity in acetylene semi-hydrogenation is often obtained at the expense of activity. To break the activity-selectivity trade-off, precise control and in-depth understanding of the three-dimensional atomic structure of surfacial active sites are crucial. Here, we designed a novel Au@PdCu core–shell nanocatalyst featuring diluted and stretched Pd sites on the ultrathin shell (1.6 nm), which showed excellent reactivity and selectivity, with 100% acetylene conversion and 92.4% ethylene selectivity at 122 °C, and the corresponding activity was 3.3 times higher than that of the PdCu alloy. The atomic three-dimensional decoding for the activity-selectivity balance was revealed by combining pair distribution function (PDF) and reverse Monte Carlo simulation (RMC). The results demonstrate that a large number of active sites with a low coordination number of Pd–Pd pairs and an average 3.25% tensile strain are distributed on the surface of the nanocatalyst, which perform a pivotal function in the simultaneous improvement of hydrogenation activity and ethylene selectivity. Our work not only develops a novel strategy for unlocking the linear scaling relation in heterogeneous catalysis but also provides a paradigm for atomic 3D understanding of lattice strain in core–shell nanocatalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Computational insights into the interaction of water with the UiO-66 metal–organic framework and its functionalized derivatives

The UiO-66 metal–organic framework (MOF) has been identified as a promising hydrophilic material for water harvesting. Recent studies show that its water uptake ability at low relative humidity (RH) can be improved by incorporating hydrophilic functional groups into the framework. In this work, we provide computational insights into the adsorption of water in UiO-66 and its functionalized derivatives to reveal the role played by different adsorption sites and functional groups in the adsorption mechanism. We started by developing molecular models for UiO-66, UiO-66-NH 2 , UiO-66-OH, and UiO-66-(OH) 2 compatible with the MB-pol data-driven many-body potential of water. We then benchmarked these models against ab initio data. Here we used these models to perform molecular dynamics simulations and calculate radial distribution functions, IR spectra, and two-dimensional density distribution maps for water in the MOFs. These results consistently show that the μ 3 -OH sites are the preferential interaction sites for water in UiO-66 and all its variants, and the formation of localised water clusters inside the octahedral pores is responsible for the abrupt step in the experimental adsorption isotherms. Furthermore, the presence of functional groups in the framework allows water to cluster in the octahedral pores at lower RH, thus making the MOF a more efficient water harvester. Overall, this study provides molecular-level insights into the pore filling process of UiO-66 and its functionalized derivatives, which are needed for the design of efficient water harvesting materials based on MOFs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Blocky Selective Postpolymerization C–H Functionalization of Polyolefins

C–H functionalization of commodity polyolefins affords functional materials derived from a high-volume, low-cost resource. However, current postpolymerization modification strategies result in randomly distributed functionalization along the length of the polymer backbone, which has a negative impact on the crystallinity of the resultant polymers, and thus the thermomechanical properties. Here, we demonstrate an amidyl radical mediated C–H functionalization of polyolefins to access blocky microstructures, which exhibit a higher crystalline fraction, larger crystallite size, and improved mechanical properties compared to their randomly functionalized analogues. Taking inspiration from the site-selective C–H functionalization of small molecules, we leverage the steric protection provided by crystallites and target polymer functionalization to amorphous domains in a semicrystalline polyolefin gel. The beneficial outcomes of blocky functionalization are independent of the identity of the pendant functional group that is installed through functionalization. Here, the decoupling of functional group incorporation and crystallinity highlights the promise in accessing nonrandom microstructures through selective functionalization to circumvent traditional tradeoffs in postpolymerization modification, with potential impact in advanced materials and upcycling plastic waste.

Neidhart, Eliza K. [The University of North Caroli↗
DOE OSTI2021

Radial structure and formation of the Milky Way disc

The formation of the Galactic disc is an enthusiastically debated issue. Numerous studies and models seek to identify the dominant physical process(es) that shaped its observed properties; for example, satellite accretion, starburst, quenching, gas infall, and stellar radial migration. Taking advantage of the improved coverage of the inner Milky Way provided by the SDSS DR16 APOGEE catalogue and of the ages published in the APOGEE-AstroNN Value Added Catalogue (VAC), we examined the radial evolution of the chemical and age properties of the Galactic stellar disc with the aim of better constraining its formation. Using a sample of 199 307 giant stars with precise APOGEE abundances and APOGEE-AstroNN ages, selected in a ±2 kpc layer around the galactic plane, we assessed the dependency as a function of guiding radius of (i) the median metallicity, (ii) the ridge lines of the [Fe/H] – [Mg/Fe] and age–[Mg/Fe] distributions, and (iii) the age distribution function (ADF). The giant star sample allows us to probe the radial behaviour of the Galactic disc from R g = 0 to 14–16 kpc. The thick disc [Fe/H] – [Mg/Fe] ridge lines follow closely grouped parallel paths, supporting the idea that the thick disc did form from a well-mixed medium. However, the ridge lines present a small drift in [Mg/Fe], which decreases with increasing guiding radius. At sub-solar metallicity, the intermediate and outer thin disc [Fe/H] – [Mg/Fe] ridge lines follow parallel sequences shifted to lower metallicity as the guiding radius increases. We interpret this pattern as the signature of a dilution of the interstellar medium from R g ~ 6 kpc to the outskirts of the disc, which occurred before the onset of the thin disc formation. The APOGEE-AstroNN VAC provides stellar ages for statistically significant samples of thin disc stars from the Galactic centre up to R g ~ 14 kpc. An important result provided by this dataset is that the thin disc presents evidence of an inside-out formation up to R g ~ 10 – 12 kpc. Moreover, about ~7 Gyr ago, the [Mg/Fe] ratio in the outer thin disc (R g > 10 kpc) was higher by about ~0.03–0.05 dex than in the more internal regions of the thin disc. This could be the fossil record of a pollution of the outer disc gas reservoir by the thick disc during its starburst phase.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Observation of Kolmogorov turbulence due to multiscale vortices in dusty plasma experiments

We report the experimental observation of fully developed Kolmogorov turbulence originating from self-excited vortex flows in a three-dimensional (3D) dust cloud. The characteristic -5/3 scaling of 3D Kolmogorov turbulence is consistent in both the spatial and temporal energy spectra within a statistical variation of experimental data. Additionally, the 2/3 scaling in the second-order structure function further supports the presence of Kolmogorov turbulence. We also identified a slight deviation in the tails of the probability distribution functions for velocity gradients, a reflection of intermittency. The experiment showed the formation of a dust cloud in the diffused plasma region away from the electrodes. The dust rotation was observed in multiple experimental campaigns under different discharge conditions at different spatial locations and background plasma environments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Measurements of the $W$ boson rapidity, helicity, double-differential cross sections, and charge asymmetry in $pp$ collisions at $\sqrt {s}$ = 13 TeV

The differential cross section and charge asymmetry for inclusive W boson production at s=13 TeV is measured for the two transverse polarization states as a function of the W boson absolute rapidity. The measurement uses events in which a W boson decays to a neutrino and either a muon or an electron. The data sample of proton-proton collisions recorded with the CMS detector at the LHC in 2016 corresponds to an integrated luminosity of 35.9 fb-1. The differential cross section and its value normalized to the total inclusive W boson production cross section are measured over the rapidity range |yW|<2.5. In addition to the total fiducial cross section, the W boson double-differential cross section, d2σ/dpTℓd|ηℓ|, and the charge asymmetry are measured as functions of the charged lepton transverse momentum and pseudorapidity. The precision of these measurements is used to constrain the parton distribution functions of the proton using the next-to-leading order NNPDF3.0 set.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Phase-space sensitivity (weight functions) of 3 MeV proton diagnostics

Abstract The weight function W ( X ) is a diagnostic sensitivity to phase-space variables X that relates the measured signal C to the distribution function F ( X ) through the equation C = ∫ W ( X ) F ( X ) d X . In the present work, an algorithm to calculate W for a diagnostic that measures 3 MeV protons produced in d(d,p)t fusion reactions between a fast ion and a thermal deuteron is developed. The emitted protons escape the tokamak on curved orbits and are detected. These curved orbits constitute effective diagnostic ‘sightlines’. The presented algorithm accounts for the complications associated with these curved sightlines. An initial calculation of time-reversed proton orbits computes effective solid angles and sightlines for the relevant range of incident proton velocity vectors. These precomputed orbits are then used within the framework of FIDASIM (2020 Plasma Phys. Cont. Fusion 62 105008) to calculate the reactivity averaged over the thermal distribution of the ‘target’ deuterons and the probability that a fast ion of specified energy and pitch has a gyroangle that is consistent with the kinematic equations along each of the sightlines. Comparisons with analytic formulas and with independent calculations for the Mega Amp Spherical Tokamak 3 MeV proton diagnostic verify the algorithm.

Physics↗
DOE OSTI2020

Spatial Engineering Direct Cooperativity between Binding Sites for Uranium Sequestration

Preorganization is a basic design principle used by nature that allows for synergistic pathways to be expressed. Herein, a full account of the conceptual and experimental development from randomly distributed functionalities to a convergent arrangement that facilitates cooperative binding is given, thus conferring exceptional affinity toward the analyte of interest. The resulting material with chelating groups populated adjacently in a spatially locked manner displays up to two orders of magnitude improvement compared to a random and isolated manner using uranium sequestration as a model application. This adsorbent shows exceptional extraction efficiencies, capable of reducing the uranium concentration from 5 ppm to less than 1 ppb within 10 min, even though the system is permeated with high concentrations of competing ions. The efficiency is further supported by its ability to extract uranium from seawater with an uptake capability of 5.01 mg g -1 , placing it among the highest-capacity seawater uranium extraction materials described to date. The concept presented here uncovers a new paradigm in the design of efficient sorbent materials by manipulating the spatial distribution to amplify the cooperation of functions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Evolution of the age-included nearest pair distribution in disperse multiphase flows

The age of the nearest particle pair is introduced as the difference between the current time and the most recent time when the nearest particle pair was formed. The evolution equation for the age-included nearest pair distribution function is derived. With the assumption of random destruction of the nearest particle pairs, the evolution equation predicts the exponential probability distribution of the ages of the nearest particle pairs. Particle-resolved numerical simulations with moving particles are performed to verify this prediction. The equation is then used to derive the evolution equation for the particle–fluid–particle (PFP) stress, which is known to be related to hyperbolicity of the two-fluid equations. It is found that the relaxation time of the age probability distribution is also the relaxation time for the PFP stress. Guided by the closure terms in the PFP stress evolution equation, we study kinematics of the nearest particle pairs in the particle-resolved simulations for flows caused by sedimentation of the particles with initially isotropic and homogeneous particle distributions. At the steady states, the particle Reynolds numbers are around 20. Anisotropy and inhomogeneity of particle distributions are seen to develop in these flows. The mean distances to the nearest particles and evolution of the distribution of the Voronoi cell volumes are studied. We also found the PFP stress is closely related to the changes in these inter-particle scale quantities.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗