Engineering Papers⌕ Search

SEARCH · Engineering Papers

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Synthesis of CS/PVA Biodegradable Composite Nanofibers as a Microporous Material with Well Controllable Procedure Through Electrospinning

In this study, we have showed a facile route for fabrication of a novel microporous material based on chitosan (CS) and poly(vinyl alcohol) (PVA) biodegradable nanofibers that have high specific surface area, considerable porosity, and small diameter. Scanning electron microscopy, thermogravimetric analysis, differential scanning calorimetry, fourier transform infrared spectroscopy, Brunauer–Emmett–Teller surface area analysis, and CHNS/O elemental analyser were applied to characterize the fabricated CS/PVA composite nanofibers. Moreover, the influences of spinning conditions including concentration, voltage, electrospinning distance, and flow rate, on size distribution and pore diameter of the final product were systematically studied using 2{sup k−1} factorial design experiments, and the response surface optimization was used for determining the best synthesis parameter. The results obtained from 2{sup K−1} factorial design experiments showed that electrospinning parameters influenced the size distribution and pore diameter of the CS/PVA microporous material. Based on the response surface methodology, the CS/PVA product could be obtained with a high microporous diameter of 1.8 nm and a small diameter distribution of 15.0 nm under optimized conditions. The obtained results showed that the fabricated samples could be utilized in different applications.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Mapping Stellar Surfaces. II. An Interpretable Gaussian Process Model for Light Curves

The use of Gaussian processes (GPs) as models for astronomical time series data sets has recently become almost ubiquitous, given their ease of use and flexibility. In particular, GPs excel at marginalization over the stellar signal when the variability due to starspots is treated as a nuisance, as in exoplanet transit modeling. However, these effective models are less useful in cases where the starspot signal is of primary interest, since it is not obvious how the parameters of the GP relate to physical parameters like the spot size, contrast, and latitudinal distribution. Instead, it is common practice to explicitly model the effect of individual starspots on the light curve and attempt to infer their properties via optimization or posterior inference. Unfortunately, this process is ill-posed and often computationally intractable when applied to stars with more than a few spots and/or to ensembles of many stars. Here we derive a closed-form expression for a GP that describes the light curve of a rotating, evolving stellar surface conditioned on a given distribution of starspot sizes, contrasts, and latitudes. We demonstrate that this model is correctly calibrated, allowing one to robustly infer physical parameters of interest from one or more light curves, including the typical spot radii and latitudes. Our GP has far-ranging implications for understanding the variability and magnetic activity of stars from light curves and radial velocity measurements, as well as for modeling correlated noise in exoplanet searches. Our implementation is efficient, user-friendly, and open-source, available in the package starry-process.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

SU-E-J-141: Activity-Equivalent Path Length Approach for the 3D PET-Based Dose Reconstruction in Proton Therapy

Purpose: Ion beam therapy is sensitive to uncertainties from treatment planning and dose delivery. PET imaging of induced positron emitter distributions is a practical approach for in vivo, in situ verification of ion beam treatments. Treatment verification is usually done by comparing measured activity distributions with reference distributions, evaluated in nominal conditions. Although such comparisons give valuable information on treatment quality, a proper clinical evaluation of the treatment ultimately relies on the knowledge of the actual delivered dose. Analytical deconvolution methods relating activity and dose have been studied in this context, but were not clinically applied. In this work we present a feasibility study of an alternative approach for dose reconstruction from activity data, which is based on relating variations in accumulated activity to tissue density variations. Methods: First, reference distributions of dose and activity were calculated from the treatment plan and CT data. Then, the actual measured activity data were cumulatively matched with the reference activity distributions to obtain a set of activity-equivalent path lengths (AEPLs) along the rays of the pencil beams. Finally, these AEPLs were used to deform the original dose distribution, yielding the actual delivered dose. The method was tested by simulating a proton therapy treatment plan delivering 2 Gy on a homogeneous water phantom (the reference), which was compared with the same plan delivered on a phantom containing inhomogeneities. Activity and dose distributions were were calculated by means of the FLUKA Monte Carlo toolkit. Results: The main features of the observed dose distribution in the inhomogeneous situation were reproduced using the AEPL approach. Variations in particle range were reproduced and the positions, where these deviations originated, were properly identified. Conclusions: For a simple inhomogeneous phantom the 3D dose reconstruction from PET-activity induced by proton beams was shown to be feasible.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Comparative Cooling Season Performance of Air Distribution Systems in Multistory Townhomes

IBACOS investigated the performance of a small-diameter high velocity heat pump system compared to a conventional system in a new construction triplex townhouse. A ductless heat pump system also was installed for comparison, but the homebuyer backed out because of aesthetic concerns about that system. In total, two buildings, having identical solar orientation and comprised of six townhomes, were monitored for comfort and energy performance. Results show that the small-diameter system provides more uniform temperatures from floor to floor in the three-story townhome. No clear energy consumption benefit was observed from either system. The builder is continuing to explore the small-diameter system as its new standard system to provide better comfort and indoor air quality. The homebuilder also explored the possibility of shifting its townhome product to meet the U.S. Department of Energy Challenge Home National Program Requirements. Ultimately, the builder decided that adoption of these practices would be too disruptive midstream in the construction cycle. However, the townhomes met the ENERGY STAR Version 3.0 program requirements.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

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

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↗
NASA NTRS1989

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.↗
NASA NTRS2017

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.↗
DOE OSTI2021

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

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

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↗
NASA NTRS1988

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.↗
NASA NTRS2005

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.↗
DOE OSTI2016

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↗
NASA NTRS2002

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↗
NASA NTRS1992

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

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.↗
DOE OSTI2019

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↗