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 1,675 records · Page 93

Parametric Study of Pulse-Combustor-Driven Ejectors at High-Pressure

Pulse-combustor configurations developed in recent studies have demonstrated performance levels at high-pressure operating conditions comparable to those observed at atmospheric conditions. However, problems related to the way fuel was being distributed within the pulse combustor were still limiting performance. In the first part of this study, new configurations are investigated computationally aimed at improving the fuel distribution and performance of the pulse-combustor. Subsequent sections investigate the performance of various pulse-combustor driven ejector configurations operating at high pressure conditions, focusing on the effects of fuel equivalence ratio and ejector throat area. The goal is to design pulse-combustor-ejector configurations that maximize pressure gain while achieving a thermal environment acceptable to a turbine, and at the same time maintain acceptable levels of NO(x) emissions and flow non-uniformities. The computations presented here have demonstrated pressure gains of up to 2.8.

Pressure Gain Combustion↗

CaCO 3 Precipitation, Transport and Sensing in Porous Media with In Situ Generation of Reactants

Ureolytically driven calcite precipitation is a promising approach for inducing subsurface mineral precipitation, but engineered application requires the ability to control and predict precipitate distribution. To study the coupling between reactant transport and precipitate distribution, columns with defined zones of immobilized urease were used to examine the distribution of calcium carbonate precipitation along the flow path, at two different initial flow rates. As expected, with slower flow precipitate was concentrated toward the upstream end of the enzyme zone and with higher flow the solid was more uniformly distributed over the enzyme zone. Under constant hydraulic head conditions the flow rate decreased as precipitates decreased porosity and permeability. The hydrolysis/precipitation zone was expected to become compressed in the upstream direction. However, apparent reductions in the urea hydrolysis rate and changes in the distribution of enzyme activity, possibly due to CaCO 3 precipitate hindering urea transport to the enzyme, or enzyme mobilization, mitigated reaction zone compression. Co-injected strontium was expected to be sequestered by coprecipitation with CaCO 3 , but the results suggested that coprecipitation was not an effective sequestration mechanism in this system. In addition, spectral induced polarization (SIP) was used to monitor the spatial and temporal evolution of the reaction zone.

54 ENVIRONMENTAL SCIENCES↗

SPATIAL DISTRIBUTION OF PAIR PRODUCTION OVER THE PULSAR POLAR CAP

Using an analytic, axisymmetric approach that includes general relativity, coupled to a condition for pair production deduced from simulations, we derive general results about the spatial distribution of pair-producing field lines over the pulsar polar cap. In particular, we show that pair production on magnetic field lines operates over only a fraction of the polar cap for an aligned rotator for general magnetic field configurations, assuming the magnetic field varies spatially on a scale that is larger than the size of the polar cap. We compare our result to force-free simulations of a pulsar with a dipole surface field and find excellent agreement. Our work has implications for first-principles simulations of pulsar magnetospheres and for explaining observations of pulsed radio and high-energy emission.

79 ASTRONOMY AND ASTROPHYSICS↗

Nanotomography for Quantitative 3D Particle Reconstruction

Particulates are ubiquitous across fuel cycle operations and carry critical information about particle formation, processing, and potential proliferation-related activities. Traditional analytical techniques, including micro-Raman spectroscopy and standard electron microscopy, are often limited in spatial resolution or dimensionality, particularly when used to examine metallic or submicron-scale features. Understanding particle morphology, phase distribution, and internal porosity is essential for constraining formation conditions, thermodynamic environments, and material transport behavior. In this report, we demonstrate the application of plasma focused ion beam nanotomography to reconstruct micron-scale particulates at nanoscale resolution. Using high-resolution backscattered electron imaging and Avizo software, we obtained 3D reconstructions that enabled quantitative analysis of particle morphology, phase composition, and internal voids. Representative examples include a Ta particle with a large central void and a composite particle with embedded tetrahedral crystalline structures. These reconstructions reveal structural and compositional details that are inaccessible through conventional 2D imaging. The results demonstrate that nanotomography provides both qualitative and quantitative insights into particle formation and behavior. Using nanotomography, porosity and phase distributions can be quantified to inform models of particle density, transport, and solidification conditions. Beyond technical insights, the workflow developed here establishes a transferable capability for analyzing heterogeneous particles and has potential applications in bulk materials studies via x-ray computed tomography or other volumetric imaging modalities. Ongoing efforts are focused on optimizing the workflow to process multiple particles simultaneously, increasing throughput and statistical robustness. Overall, this work illustrates the power of nanotomography as a tool for connecting particulate morphology to formation mechanisms, composition, and transport, thereby strengthening analytical capabilities for nuclear forensics, fuel cycle analysis, and related scientific investigations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Inviscid, nonadiabatic flow fields over blunt, sonic corner bodies for outer planet entry conditions by a method of integral relations

An investigation has been made into the ability of a method of integral relations to calculate inviscid zero degree angle of attack, radiative heating distributions over blunt, sonic corner bodies for some representative outer planet entry conditions is investigated. Comparisons have been made with a more detailed numerical method, a time asymptotic technique, using the same equilibrium chemistry and radiation transport subroutines. An effort to produce a second order approximation (two-strip) method of integral relations code to aid in this investigation is also described and a modified two-strip routine is presented. Results indicate that the one-strip method of integral relations cannot be used to obtain accurate estimates of the radiative heating distribution because of its inability to resolve thermal gradients near the wall. The two-strip method can sometimes be used to improve these estimates; however, the two-strip method has only a small range of conditions over which it will yield significant improvement over the one-strip method.

Gnoffo, P. A.↗

On the Modeling of Wall Heat Flux Partitioning in Subcooled Flow Boiling

Accurate prediction of void distribution in reactor coolant channels at subcooled boiling conditions is an important factor in the analysis of the operating conditions of both pressurized water reactors (PWRs) and boiling water reactors (BWRs). Its importance arises from purely thermal-hydraulic considerations (thermal limits, CRUD formation, etc.) and the impact of void reactivity feedback on core neutronics. A mechanistic multidimensional model of subcooled boiling was originally developed by Kurul and Podowski. Whereas this model has been implemented in several commercial CFD codes and successfully used in computer simulations of flow and phase distributions in the coolant channels of light water reactor, additional new modeling aspects of forced-convection subcooled boiling heat transfer have been investigated in recent studies. The purpose of this paper is to formulate a modified model of wall heat flux partitioning, which allows for capturing and quantifying the effect of simultaneous evaporation and condensation at the heated wall on the net evaporation rate and wall temperature. The new model has been implemented in the NPHASE-CMFD code, numerically tested and validated against experimental data. A mechanistic wall heat flux partitioning model for subcooled flow boiling has been presented. This model has been combined with extensively tested models of interfacial mass, momentum and energy transfer terms. It turns out that the model predictions of the wall temperature and void fraction are in good agreement with the experimental data of reference. Future work will focus on the continuation of model validation against experimental data, and on extending its scope to include the phenomena governing departure from nucleate boiling. (authors)

42 ENGINEERING↗

Bound-Preserving Discontinuous Galerkin Methods for Conservative Phase Space Advection in Curvilinear Coordinates

We extend the positivity-preserving method of Zhang & Shu [49] to simulate the advection of neutral particles in phase space using curvilinear coordinates. The ability to utilize these coordinates is important for non-equilibrium transport problems in general relativity and also in science and engineering applications with specific geometries. The method achieves high-order accuracy using Discontinuous Galerkin (DG) discretization of phase space and strong stabilitypreserving, Runge-Kutta (SSP-RK) time integration. Special care in taken to ensure that the method preserves strict bounds for the phase space distribution function f; i.e., f ϵ [0, 1]. The combination of suitable CFL conditions and the use of the high-order limiter proposed in [49] is su cient to ensure positivity of the distribution function. However, to ensure that the distribution function satisfies the upper bound, the discretization must, in addition, preserve the divergencefree property of the phase space ow. Proofs that highlight the necessary conditions are presented for general curvilinear coordinates, and the details of these conditions are worked out for some commonly used coordinate systems (i.e., spherical polar spatial coordinates in spherical symmetry and cylindrical spatial coordinates in axial symmetry, both with spherical momentum coordinates). Results from numerical experiments - including one example in spherical symmetry adopting the Schwarzschild metric - demonstrate that the method achieves high-order accuracy and that the distribution function satisfies the maximum principle.

97 MATHEMATICS AND COMPUTING↗

Investigating the domain of validity of the Gubser solution to the Boltzmann equation

Here, we study the evolution of the one particle distribution function that solves exactly the relativistic Boltzmann equation within the relaxation time approximation for a conformal system undergoing simultaneously azimuthally symmetric transverse and boost-invariant longitudinal expansion. We show, for arbitrary values of the shear viscosity to entropy density ratio, that the distribution function can become negative in certain kinematic regions of the available phase space depending on the boundary conditions. For thermal equilibrium initial conditions, we determine numerically the physical boundary in phase space where the distribution function is always positive definite. The requirement of positivity of this particular exact solution restricts its domain of validity, and it imposes physical constraints on its applicability.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Iron - Whence it came, where it went

Determinations of the abundances of iron and related elements in the photosphere, chromosphere, and corona of the sun as well as in solar and galactic cosmic rays are reviewed and compared with abundances derived from meteoritic data. Observed solar system abundances are found to be in accord with predictions of nucleosynthesis under either hydrostatic or explosive conditions but cannot yet be used to define these processes uniquely. The distribution of iron among planets and meteorites can probably be adequately modelled by condensation and fragmentation under equilibrium conditions above about 700 K, but below that temperature it is likely that inhibited solid-state diffusion perturbed attainment of equilibrium. Pertinent factors which are presently unknown include the mechanism responsible for metal-silicate fractionation, the grain size achieved by metallic iron in the nebula, and whether iron-bearing silicate formed prior to accretion.

Kerridge, J. F.↗

Quantitative characterization of spatial distribution of particles in materials: Application to materials processing

Most engineering materials contain second phase particles or fibers which serve to reinforce the matrix phase. The effect of reinforcements on material properties is usually analyzed in terms of the average volume fraction and spacing of reinforcements, quantities which are global microstructural characteristics. However, material properties can also depend on local microstructural characteristics; for example, on how uniformly the reinforcing phase is distributed in the material. The analysis method will then be applied to a materials processing problem to discover how processing parameters can be selected to maximize redistribution of the reinforcing phase during processing. Several mathematical analysis methods could be adapted to the problem of characterizing the distribution of particles in materials. A tessellation-based method was selected. In the first phase of the investigation, a software package was written to automate the analysis. Typical results are shown. The analysis technique allows the degree to which particles are clustered together, the size and spacing of particle clusters, and the particle density in clusters to be found. The analysis methods were applied to computer-generated distributions and to a few real particle-containing materials. Methods for analyzing a nonuniform particle distribution in a material can be applied to two broad classes of materials science problems: understanding how the resulting particle distribution affects properties. The analysis method described is applied to a materials processing problem: how to select extrusion conditions to maximize the redistribution of reinforcing particles that are initially nonuniformly distributed. In addition, the tessellation-based method to analyze star distributions in spiral galaxies was adapted, illustrating the diverse types of problems to which the analysis method can be applied.

Parse, J. B.↗

Airplane Spinning

A summary on spinning is presented to point out the state of the art and the most important parameters', and to show the effects of these parameters on the spin and spin-recovery characteristics. The discussion presented applies to the fully developed spin, but does not apply to the spin entry, or incipient spin. The principal factors in spinning are mass distribution, which is by far the most important single parameter, and tail design, which is particularly important for conditions of zero or near-zero loading. By knowing the mass distribution and tail design, it is possible in many cases to predict whether an airplane has satisfactory spin-recovery characteristics. In other cases, however, it is necessary to make spin tests to assure satisfactory recovery.

Bowman, James S.↗

Numerical Modeling of Liquid-Vapor Phase Change

We implemented a two- and three-dimensional finite difference/front tracking technique to solve liquid-vapor phase change problems. The mathematical and the numerical features of the method were explained in great detail in our previous reports, Briefly, we used a single formula representation which incorporated jump conditions into the governing equations. The interfacial terms were distributed as singular terms using delta functions so that the governing equations would be the same as conventional conservation equations away from the interface and in the vicinity of the interface they would provide correct jump conditions. We used a fixed staggered grid to discretize these equations and an unstructured grid to explicitly track the front. While in two dimensions the front was simply a connection of small line segments, in three dimensions it was represented by a connection of small triangular elements. The equations were written in conservative forms and during the course of computations we used regriding to control the size of the elements of the unstructured grid. Moreover, we implemented a coalescence in two dimensions which allowed the merging of different fronts or two segments of the same front when they were sufficiently close. We used our code to study thermocapillary migration of bubbles, burst of bubbles at a free surface, buoyancy-driven interactions of bubbles, evaporation of drops, rapid evaporation of an interface, planar solidification of an undercooled melt, dendritic solidification, and a host of other problems cited in the reference.

Esmaeeli, Asghar↗

Deformation of pores in response to uniaxial and hydrostatic stress cycling in Marcellus Shale: Implications for gas recovery

One of the main challenges during gas production from shale reservoirs is low recovery rate. One contributing factor to this outcome is an insufficient understanding of pore systems, especially pore behaviour following changes in reservoir conditions or resulting from gas production practices. Because the pressure in the producing well can be controlled, understanding the effects of pressure variation on the pore size distribution and methane trapping is necessary to help design optimal conditions to improve the gas recovery rate. This work is the first systematic study of sub-millimeter pore deformation in shale caused by uniaxial and hydrostatic stress up to 100 MPa. Overmature samples from the Middle Devonian Marcellus Shale were analyzed using neutron scattering (SANS and USANS) techniques to interpret the response of nanopores to stress cycling of magnitude and duration compatible with the hydraulic fracturing treatments. Experiments reported here are performed at a series of uniaxial pressures up to 100 MPa and at hydrostatic pressures of deuterated methane 0 and 50 MPa. Since at the original depth of the shale samples’ burial of 2184 m the hydrostatic pressure is approximately 27 MPa and the lithostatic pressure is about 55 MPa, the experimental conditions reasonably well simulate the reservoir pressure regime. Our SANS and USANS results show that different pore sizes are affected by uniaxial stress in different ways. Specifically, in the pore size range from 1 nm to 800 nm, a decrease of pore density with pressure is observed, with the most depleted being mesopores of about 100 nm in diameter. The observed decrease is likely related to deformation of kerogen, followed by a loss of pore nano-volume, as well as methane trapped in the micropores. For pores larger than 5 μm, USANS data suggest that the negative trend is reversed at above 74 MPa, and the number density of large macropores may increase with increased stress even above the original value. The increased number of macropores at high pressure may create new interconnected conduits for gas migration, resulting in a better recovery rate. Another important finding of this study is an irreversible rearrangement of pore size distribution taking place after pressure cycling. Furthermore, this irreversible reorganization of pore size distribution should be taken into account during management of well production to maximize recovery rate.

03 NATURAL GAS↗

State-space receptive fields of semicircular canal afferent neurons in the bullfrog

Receptive fields are commonly used to describe spatial characteristics of sensory neuron responses. They can be extended to characterize temporal or dynamical aspects by mapping neural responses in dynamical state spaces. The state-space receptive field of a neuron is the probability distribution of the dynamical state of the stimulus-generating system conditioned upon the occurrence of a spike. We have computed state-space receptive fields for semicircular canal afferent neurons in the bullfrog (Rana catesbeiana). We recorded spike times during broad-band Gaussian noise rotational velocity stimuli, computed the frequency distribution of head states at spike times, and normalized these to obtain conditional pdfs for the state. These state-space receptive fields quantify what the brain can deduce about the dynamical state of the head when a single spike arrives from the periphery. c2001 Elsevier Science B.V. All rights reserved.

Non-NASA Center↗

Determining Water Content and Distribution in PEMFCs to Predict Aging While in Storage

Proton membrane exchange fuel cells (PEMFCs) have the potential to be long term backup power sources with a startup time on the order of seconds. Water management is the key issue in being able to successfully store PEMFCs for extended periods of time. In this work custom made PEMFCs were humidified at various relative humidities (%RH) and subsequently stored for different lengths of time. The fuel cell’s water content was then imaged at the National Institute of Standards and Technology (NIST) neutron imaging facility. In conclusion, the cells’ startup performances were measured simulating quick startup conditions to define the effect of different water distributions.

25 ENERGY STORAGE↗

Predicting a contact's sensitivity to initial conditions using metrics of frictional coupling

This paper presents a method for predicting how sensitive a frictional contact’s steady-state behavior is to its initial conditions. Previous research has proven that if a contact is uncoupled , i.e. if slip displacements do not influence the contact pressure distribution, then its steady-state response is independent of initial conditions, but if the contact is coupled, the steady-state response depends on initial conditions. In this paper, two metrics for quantifying coupling in discrete frictional systems are examined. These metrics suggest that coupling is dominated by material dissimilarity due to Dundurs’ composite material parameter β when β ≥ 0.2, but geometric mismatch becomes the dominant source of coupling for smaller values of β. Based on a large set of numerical simulations with different contact geometries, material combinations, and friction coefficients, a contact’s sensitivity to initial conditions is found to be correlated with the product of the coupling metric and the friction coefficient. For cyclic shear loading, this correlation is maintained for simulations with different contact geometries, material combinations, and friction coefficients. Furthermore, for cyclic bulk loading, the correlation is only maintained when the contact edge angle is held constant.

42 ENGINEERING↗

ERS-1 Investigations of Southern Ocean Sea Ice Geophysics Using Combined Scatterometer and SAR Images

Coregistered ERS-1 SAR and Scatterometer data are presented for the Weddell Sea, Antarctica. Calibrated image backscatter statistics are extracted from data acquired in regions where surface measurements were made during two extensive international Weddell Sea experiments in 1992. Changes in summer ice-surface conditions, due to temperature and wind, are shown to have a large impact on observed microwave backscatter values. Winter calibrated backscatter distributions are also investigated as a way of describing ice thickness conditions in different location during winter. Coregistered SAR and EScat data over a manned drifting ice station are used to illustrate the seasonal signature changes occurring during the fall freeze-up transition.

sea-ice Weddell Sea ice formation↗