Search NASA⌕ Search

SEARCH · Search NASA

Results for “Conditional distributions”

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 595 records · Page 33

The distribution of neutral hydrogen in the colour–magnitude plane of galaxies

ABSTRACT We present the conditional H i (neutral hydrogen) mass function (HIMF) conditioned on observed optical properties, Mr (r-band absolute magnitude), and Cur (u − r colour), for a sample of 7709 galaxies from Arecibo Legacy Fast ALFA (40 per cent data release – α.40) which overlaps with a common volume in SDSS DR7. Based on the conditional HIMF, we find that the luminous red, luminous blue, and faint blue populations dominate the total HIMF at the high-mass end, knee, and the low-mass end, respectively. We use the conditional HIMF to derive the underlying distribution function of ΩH i (H i density parameter), p(ΩH i), in the colour–magnitude plane of galaxies. The distribution, p(ΩH i), peaks in the blue cloud at $\mathit{ M_{{r}}}^{\text{max}}=-19.25, \mathit{ C_{{ur}}}^{\text{max}}=1.44$ but is skewed. It has a long tail towards faint blue galaxies, and luminous red galaxies. We argue that p(ΩH i) can be used to reveal the underlying relation between cold gas, stellar mass, and the star formation rate in an unbiased way, that is, the derived relation does not suffer from survey or sample selection.

Dutta, Saili (ORCID:0000000288585845)↗

Learning the temporal evolution of multivariate densities via normalizing flows

In this work, we propose a method to learn multivariate probability distributions using sample path data from stochastic differential equations. Specifically, we consider temporally evolving probability distributions (e.g., those produced by integrating local or nonlocal Fokker–Planck equations). Here, we analyze this evolution through machine learning assisted construction of a time-dependent mapping that takes a reference distribution (say, a Gaussian) to each and every instance of our evolving distribution. If the reference distribution is the initial condition of a Fokker–Planck equation, what we learn is the time-T map of the corresponding solution. Specifically, the learned map is a multivariate normalizing flow that deforms the support of the reference density to the support of each and every density snapshot in time. We demonstrate that this approach can approximate probability density function evolutions in time from observed sampled data for systems driven by both Brownian and Lévy noise. We present examples with two- and three-dimensional, uni- and multimodal distributions to validate the method.

97 MATHEMATICS AND COMPUTING↗

Communication requirements of sparse Cholesky factorization with nested dissection ordering

Load distribution schemes for minimizing the communication requirements of the Cholesky factorization of dense and sparse, symmetric, positive definite matrices on multiprocessor systems are presented. The total data traffic in factoring an n x n sparse symmetric positive definite matrix representing an n-vertex regular two-dimensional grid graph using n exp alpha, alpha not greater than 1, processors are shown to be O(n exp 1 + alpha/2). It is O(n), when n exp alpha, alpha not smaller than 1, processors are used. Under the conditions of uniform load distribution, these results are shown to be asymptotically optimal.

Naik, Vijay K.↗

Numerical Simulation of Illumination and Thermal Conditions at the Lunar Poles Using LOLA DTMs

We are interested in illumination conditions and the temperature distribution within the upper two meters of regolith near the lunar poles. Here, areas exist receiving almost constant illumination near areas in permanent shadow, which were identified as potential exploration sites for future missions. For our study a numerical simulation of the illumination and thermal environment for lunar near-polar regions is needed. Our study is based on high-resolution, twenty meters per pixel and 400 x 400 km large polar Digital Terrain Models (DTMs), which were derived from Lunar Orbiter Laser Altimeter (LOLA) data. Illumination conditions were simulated by synthetically illuminating the LOLA DTMs using the horizon method considering the Sun as an extended source. We model polar illumination for the central 50 x 50 km subset and use it as an input at each time-step (2 h) to evaluate the heating of the lunar surface and subsequent conduction in the sub-surface. At surface level we balance the incoming insolation with the subsurface conduction and radiation into space, whereas in the sub-surface we consider conduction with an additional constant radiogenic heat source at the bottom of our two-meter layer. Density is modeled as depth-dependent, the specific heat parameter as temperature-dependent and the thermal conductivity as depth- and temperature-dependent. We implemented a fully implicit finite-volume method in space and backward Euler scheme in time to solve the one-dimensional heat equation at each pixel in our 50 x 50 km DTM. Due to the non-linear dependencies of the parameters mentioned above, Newton's method is employed as the non-linear solver together with the Gauss-Seidel method as the iterative linear solver in each Newton iteration. The software is written in OpenCL and runs in parallel on the GPU cores, which allows for fast computation of large areas and long time scales.

Glaser, P.↗

Experimental decoy-state Bennett-Brassard 1984 quantum key distribution through a turbulent channel

In free-space quantum key distribution (QKD) in turbulent conditions, scattering and beam wandering cause intensity fluctuations which decrease the detected signal-to-noise ratio. This effect can be mitigated by rejecting received bits when the channel's transmittance is below a threshold. Thus, the overall error rate is reduced and the secure key rate increases despite the deletion of bits. In this work, we implement recently proposed selection methods focusing on the prefixed-threshold real-time selection (P-RTS) where a cutoff can be chosen prior to data collection and independently of the transmittance distribution. We perform finite-size decoy-state Bennett-Brassard 1984 QKD in a laboratory setting where we simulate the atmospheric turbulence using an acousto-optical modulator. We show that P-RTS can yield considerably higher secure key rates for a wide range of the atmospheric channel parameters. In addition, we evaluate the performance of the P-RTS method for a realistically finite sample size. We demonstrate that a near-optimal selection threshold can be predetermined even with imperfect knowledge of the channel transmittance distribution parameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On the Morphodynamics of a Wide Class of Large‐Scale Meandering Rivers: Insights Gained by Coupling LES With Sediment‐Dynamics

Abstract In meandering rivers, interactions between flow, sediment transport, and bed topography affect diverse processes, including bedform development and channel migration. Predicting how these interactions affect the spatial patterns and magnitudes of bed deformation in meandering rivers is essential for various river engineering and geoscience problems. Computational fluid dynamics simulations can predict river morphodynamics at fine temporal and spatial scales but have traditionally been challenged by the large scale of natural rivers. We conducted coupled large‐eddy simulation and bed morphodynamics simulations to create a unique database of hydro‐morphodynamic data sets for 42 meandering rivers with a variety of planform shapes and large‐scale geometrical features that mimic natural meanders. For each simulated river, the database includes (a) bed morphology, (b) three‐dimensional mean velocity field, and (c) bed shear stress distribution under bankfull flow conditions. The calculated morphodynamics results at dynamic equilibrium revealed the formation of scour and deposition patterns near the outer and inner banks, respectively, while the location of point bars and scour regions around the apexes of the meander bends is found to vary as a function of the radius of curvature of the bends to the width ratio. A new mechanism is proposed that explains this seemingly paradoxical finding. The high‐fidelity simulation results generated in this work provide researchers and scientists with a rich numerical database for morphodynamics and bed shear stress distributions in large‐scale meandering rivers to enable systematic investigation of the underlying phenomena and support a range of river engineering applications.

54 ENVIRONMENTAL SCIENCES↗

Extension of moment projection method to the fragmentation process

The method of moments is a simple but efficient method of solving the population balance equation which describes particle dynamics. Recently, the moment projection method (MPM) was proposed and validated for particle inception, coagulation, growth and, more importantly, shrinkage; here the method is extended to include the fragmentation process. The performance of MPM is tested for 13 different test cases for different fragmentation kernels, fragment distribution functions and initial conditions. Comparisons are made with the quadrature method of moments (QMOM), hybrid method of moments (HMOM) and a high-precision stochastic solution calculated using the established direct simulation algorithm (DSA) and advantages of MPM are drawn.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Data traffic reduction schemes for sparse Cholesky factorizations

Load distribution schemes are presented which minimize the total data traffic in the Cholesky factorization of dense and sparse, symmetric, positive definite matrices on multiprocessor systems with local and shared memory. The total data traffic in factoring an n x n sparse, symmetric, positive definite matrix representing an n-vertex regular 2-D grid graph using n (sup alpha), alpha is equal to or less than 1, processors are shown to be O(n(sup 1 + alpha/2)). It is O(n(sup 3/2)), when n (sup alpha), alpha is equal to or greater than 1, processors are used. Under the conditions of uniform load distribution, these results are shown to be asymptotically optimal. The schemes allow efficient use of up to O(n) processors before the total data traffic reaches the maximum value of O(n(sup 3/2)). The partitioning employed within the scheme, allows a better utilization of the data accessed from shared memory than those of previously published methods.

Naik, Vijay K.↗

Constraints on Grain Formation Around Carbon Stars from Laboratory Studies of Presolar Graphite

We report the results of an investigation into the physical conditions in the mass outflows of asymptotic giant branch (AGB) carbon stars that are required for the formation of micron-sized presolar graphite grains, either with or without internal crystals of titanium carbide (TiC). In addition to providing detailed information about stellar nucleosynthesis, the structure and composition of presolar grains give unique information about the conditions of grain formation. In the present work we use laboratory observations of presolar graphite to gain insight into the physical conditions in circumstellar outflows from carbon AGB stars. The periodic pulsation of AGB stars enhances the gas density through shocks in the stellar atmosphere above the photosphere, promoting the condensation of dust grains. Copious mass outflow occurs largely because grains are coupled to the radiation field of the star, which accelerates them by radiation pressure; momentum is in turn transferred to gas molecules by collisions with grains. The dust/gas mixture is effectively a two-component fluid whose motion depends on atmospheric structure and which, in turn, influences that structure. In particular, the radiation pressure on the grains determines the velocity field of the outflow and thus the density distribution, while the density distribution itself determines the conditions of radiative transfer within the outflow and thus the effective radiation pressure.

Bernatowicz, T. J.↗

Wave-optical assessment of alignment tolerances in nano-focusing with ellipsoidal mirror

High-precision ellipsoidal mirrors, which can efficiently focus X-rays to the nanometer dimension with a mirror, have not been realized because of the difficulties in the fabrication process. The purpose of our study was to develop nano-focusing ellipsoidal mirrors in the hard X-ray region. We developed a wave-optical focusing simulator for investigating alignment tolerances in nano-focusing with a designed ellipsoidal mirror, which produce a diffraction-limited focus size of 30 × 35 nm{sup 2} in full width at half maximum at an X-ray energy of 7 keV. The simulator can calculate focusing intensity distributions around the focal point under conditions of misalignment. The wave-optical simulator enabled the calculation of interference intensity distributions, which cannot be predicted by the conventional ray-trace method. The alignment conditions with a focal length error of ≲ ±10 µm, incident angle error of ≲ ±0.5 µrad, and in-plane rotation angle error of ≲ ±0.25 µrad must be satisfied for nano-focusing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On the Exit Boundary Condition for One-Dimensional Calculations of Pulsed Detonation Engine Performance

In one-dimensional calculations of pulsed detonation engine (PDE) performance, the exit boundary condition is frequently taken to be a constant static pressure. In reality, for an isolated detonation tube, after the detonation wave arrives at the exit plane, there will be a region of high pressure, which will gradually return to ambient pressure as an almost spherical shock wave expands away from the exit, and weakens. Initially, the flow is supersonic, unaffected by external pressure, but later becomes subsonic. Previous authors have accounted for this situation either by assuming the subsonic pressure decay to be a relaxation phenomenon, or by running a two-dimensional calculation first, including a domain external to the detonation tube, and using the resulting exit pressure temporal distribution as the boundary condition for one-dimensional calculations. These calculations show that the increased pressure does affect the PDE performance. In the present work, a simple model of the exit process is used to estimate the pressure decay time. The planar shock wave emerging from the tube is assumed to transform into a spherical shock wave. The initial strength of the spherical shock wave is determined from comparison with experimental results. Its subsequent propagation, and resulting pressure at the tube exit, is given by a numerical blast wave calculation. The model agrees reasonably well with other, limited, results. Finally, the model was used as the exit boundary condition for a one-dimensional calculation of PDE performance to obtain the thrust wall pressure for a hydrogen-air detonation in tubes of length to diameter ratio (L/D) of 4, and 10, as well as for the original, constant pressure boundary condition. The modified boundary condition had no performance impact for values of L/D > 10, and moderate impact for L/D = 4.

Wilson, Jack↗

Diffusive crack growth at a bimaterial interface

The diffusional microcrack growth behavior in a bimaterial system is investigated with an aim at estimating service life of advanced ceramic composites under creep-rupture conditions. The crack is assumed to grow via a coupled surface and grain-boundary diffusion under steady state conditions. The tensile stress distribution along the interface ahead of the moving crack tip is solved, and it is found that a new length parameter exists as a scaling factor for which the solution becomes identical to the single phase case when plotted on the nondimensional physical plane. In contrast to the elastic stress solution which shows singularity at the tip, together with oscillatory character away from the tip, the creep stresses have a peak value away from the tip due to a wedging effect and interfacial sliding eliminates stress oscillation resulting in a decoupling between mode I and mode II loading. This solution ties the far-field loading parameter to the crack tip conditions in terms of the unknown crack velocity to give a specific V-K(1) relationship. It is shown that an exponent of 12 in the conventional crack growth power law emerges at the higher applied stress range.

Chuang, Tze-Jer↗

Resiliency of Distribution Systems Incorporating Asynchronous Information for System Restoration

Resiliency of distribution systems under extreme operating conditions is critical, especially when the utility is not available. With the large-scale deployment of distributed resources, it becomes possible to restore critical loads with local non-utility resources. Distribution system operators (DSOs) need to determine the critical loads to be restored, considering limited resources and distribution facilities. Several studies on resiliency have been conducted for the restoration of distribution systems. However, the inherent asynchronous characteristic on the information availability has not been incorporated. With incomplete and asynchronous information, decisions may be made that result in underutilization of generation resources. In this paper, a new distribution system restoration approach is proposed, considering uncertain devices and associated asynchronous information. It uses a two-module architecture that efficiently optimizes restoration actions using a binary linear programming model and evaluates their feasibility with unbalanced optimal power flow. Networked microgrids are included in the model. The IEEE 123-node test feeder is used for validation. Results show that asynchronous messages may affect the restoration actions significantly and the impacts can be mitigated by the proposed decision support tool for the DSOs.

Bedoya, Juan C.↗

Electron Energy Partition across Interplanetary Shocks. I. Methodology and Data Product

Analyses of 15,314 electron velocity distribution functions (VDFs) within ±2 hr of 52 interplanetary (IP) shocks observed by the Wind spacecraft near 1 au are introduced. The electron VDFs are fit to the sum of three model functions for the cold dense core, hot tenuous halo, and field-aligned beam/strahl component. The best results were found by modeling the core as either a bi-kappa or a symmetric (or asymmetric) bi-self-similar VDF, while both the halo and beam/strahl components were best fit to bi-kappa VDF. This is the first statistical study to show that the core electron distribution is better fit to a self-similar VDF than a bi-Maxwellian under all conditions. The self-similar distribution deviation from a Maxwellian is a measure of inelasticity in particle scattering from waves and/or turbulence. The ranges of values defined by the lower and upper quartiles for the kappa exponents are κ ec ~ 5.40–10.2 for the core, κ eh ~ 3.58–5.34 for the halo, and κ eb ~ 3.40–5.16 for the beam/strahl. The lower-to-upper quartile range of symmetric bi-self-similar core exponents is s ec ~ 2.00–2.04, and those of asymmetric bi-self-similar core exponents are p ec ~ 2.20–4.00 for the parallel exponent and q ec ~ 2.00–2.46 for the perpendicular exponent. The nuanced details of the fit procedure and description of resulting data product are also presented. The statistics and detailed analysis of the results are presented in Paper II and Paper III of this three-part study.

79 ASTRONOMY AND ASTROPHYSICS↗

Evolution of deep-bed filtration of engine exhaust particulates with trapped mass

Micro-scale filtration experiments were performed on cordierite filter samples using particulate matter (PM) generated by a spark-ignition direct-injection (SIDI) engine fueled with tier II EEE certification gasoline. Size-resolved mass and number concentrations were obtained from several engine operating conditions. The resultant mass-mobility relationships showed weak dependence on the operating condition. An integrated particle size distribution (IPSD) method was used estimate the PM mass concentration in the exhaust stream from the SIDI engine and a heavy duty diesel (HDD) engine. The average estimated mass concentration between all conditions was ~77****** % of the gravimetric measurements performed on Teflon filters. Despite the relatively low elemental carbon fraction (~0.4 to 0.7), the IPSD mass for stoichiometric SIDI exhaust was ~83±38 % of the gravimetric measurement. Identical cordierite filter samples with properties representative of diesel particulate filters were sequentially loaded with PM from the different SIDI engine operating conditions, in order of increasing PM mass concentration. Simultaneous particle size distribution measurements upstream and downstream of the filter sample were used to evaluate filter performance evolution and the instantaneous trapped mass within the filter for two different filter face velocities. The evolution of filtration performance for the different samples was sensitive only to trapped mass, despite using PM from a wide range of operating conditions. Higher filtration velocity resulted in a more rapid shift in the most penetrating particle size towards smaller mobility diameters.

33 ADVANCED PROPULSION SYSTEMS↗

Generation of Kappa Distributions in Solar Wind at 1 au

We examine the generation of kappa distributions in the solar wind plasma near 1 au. Several mechanisms are mentioned in the literature, each characterized by a specific relationship between the solar wind plasma features, the interplanetary magnetic field (IMF), and the kappa index-the parameter that governs the kappa distributions. This relationship serves as a signature condition that helps the identification of the mechanism in the plasma. In general, a mechanism that generates kappa distributions involves a single or a series of stochastic or physical processes that induces local correlations among particles. We identify three fundamental solar wind plasma conditions that can generate kappa distributions, noted as (i) Debye shielding, (ii) frozen IMF, and (iii) temperature fluctuations, each one prevailing in different scales of solar wind plasma and magnetic field properties. Moreover, our findings show that the kappa distributions, and thus, their generating mechanisms, vary significantly with solar wind features: (i) the kappa index has different dependence on the solar wind speed for slow and fast modes, i.e., slow wind is characterized by a quasi-constant kappa index, K (is) approx. 4.3 +/- 0.7, while fast wind exhibits kappa indices that increase with bulk speed; (ii) the dispersion of magnetosonic waves is more effective for lower kappa indices (i.e., further from thermal equilibrium); and (iii) the kappa and polytropic indices are positively correlated, as it was anticipated by the theory.

Livadiotis, G.↗

Negative-ion formation in the explosives RDX, PETN, and TNT using the Reversal Electron Attachment Detection (READ) technique

In the search for high sensitivity and direct atmospheric sampling of trace species, techniques have been developed such as atmospheric-sampling, glow-discharge ionization (ASGDI), corona discharge, atmospheric pressure ionization (API), electron-capture detection (ECD), and negative-ion chemical ionization (NICI) that are capable of detecting parts-per-billion to parts-per-trillion concentrations of trace species. These techniques are based on positive- or negative-ion formation via charge-transfer to the target, or electron capture under multiple-collision conditions in a Maxwellian distribution of electron energies at the source temperature. One drawback of the high-pressure, corona- or glow-discharge devices is that they are susceptible to interferences either through indistinguishable product masses, or through undesired ion-molecule reactions. The ASGDI technique is relatively immune from such interferences, since at target concentrations of less than 1 ppm the majority of negative ions arises via electron capture rather than through ion-molecule chemistry. A drawback of the conventional ECD, and possibly of the ASGDI, is that they exhibit vanishingly small densities of electrons with energies in the range 0-10 millielectron volts (meV), as can be seen from a typical Maxwellian electron energy distribution function at T = 300 K. Slowing the electrons to these subthermal (less than 10 meV) energies is crucial, since the cross section for attachment of several large classes of molecules is known to increase to values larger than 10(exp -12) sq cm at near-zero electron energies. In the limit of zero energy these cross sections are predicted to diverge as epsilon(exp -1/2), where epsilon is the electron energy. In order to provide a better 'match' between the electron energy distribution function and attachment cross section, a new concept of attachment in an electrostatic mirror was developed. In this scheme, electrons are brought to a momentary halt by reversing their direction with electrostatic fields. At this turning point the electrons have zero or near-zero energy. A beam of target molecules is introduced, and the resultant negative ions extracted. This basic idea has been recently improved to allow for better reversal geometry, higher electron currents, lower backgrounds, and increased negative-ion extraction efficiency. We present herein application of the so-called reversal electron attachment detector (READ) to the study of negative-ion formation in the explosives molecules RDX, PETN, and TNT under single-collision conditions.

Chutijian, Ara↗

PLTEMP/ANL Verification for Assemblies Having Five-Layer Fuel Tubes

PLTEMP/ANL v4.2 is a steady-state single-phase (light or heavy water) code used by ANL and by research laboratories in several countries since 1980 in the analysis of various research reactors whose cores contain assemblies made of multiple fuel plates or nested fuel tubes. This analysis is typically needed in the conversion of cores from highly-enriched uranium to low-enriched uranium. PLTEMP calculates temperature distribution and safety margins, i.e., onset of nucleate boiling ratio (ONBR), onset of Ledinegg flow instability ratio (OFIR), and critical heat flux ratio (CHFR). The analytical method in the code for calculating radial and axial temperature distributions in fuel assemblies having multiple nested coaxial fuel tubes was recently extended and verified to allow each tube to have five layers instead of three. The method also works for fuel assemblies having flat fuel plates by setting the radii of curvature of plates to large values ({approx}50 m). The method models heat sources in all five layers and in the coolant channels. To verify the extended method, an Excel spreadsheet was developed using an alternative solution method to calculate the temperature distribution in a reactor of two fuel assemblies, each made of two five-layer flat fuel plates. The code and the spreadsheet solution methods are described and their solutions are compared.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗