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,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

The Exchange of Mass and Angular Momentum in the Impact Event of Ice Giant Planets: Implications for the Origin of Uranus

Uranus has a tilted rotation axis, which is supposed to have been caused by a giant impact. In general, an impact event also changes the internal compositional distribution and drives mass ejection from the planet, which may provide the origin of satellites. Previous studies of the impact simulation of Uranus investigated the resultant angular momentum and the ejected mass distribution. However, the effect of changing the initial condition of the thermal and compositional structure is not studied. In this paper, we perform hydrodynamics simulations for the impact events of Uranus-size ice giants composed of a water core surrounded by a hydrogen envelope using two variant methods of the smoothed particle hydrodynamics. We find that the higher-entropy target loses its envelope more efficiently than the low-entropy target. However, the higher-entropy target gains more angular momentum than the lower-entropy target since the higher-entropy target has a more expanded envelope. We discuss the efficiency of angular momentum transport and the amount of the ejected mass and find a simple analytical model to roughly reproduce the outcomes of numerical simulations. We suggest the range of possible initial conditions for the giant impact on proto-Uranus that reproduces the present rotation tilt of Uranus and sufficiently provides the total angular momentum of the satellite system that can be created from the fragments from the giant impact.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1999

Experimental Study of Unsupported Nonane fuel Droplet Combustion in Microgravity

Soot formation in droplet flames is the basic component of the particulate emission process that occurs in spray combustion. The complexity of soot formation motivates a one-dimensional transport condition which has obvious advantages in modeling. Recent models of spherically symmetric droplet combustion have made this assumption when incorporating such aspects as detailed chemistry and radiation. Interestingly, spherical symmetry does not necessarily restrict the results because it has been observed that the properties of carbon formed in flames are not strongly affected by the nature of the fuel or flaming configuration. What is affected, however, are the forces acting on the soot aggregates and where they are trapped by a balance of drag and thermophoretic forces. The distribution of these forces depends on the transport conditions of the flame. Prior studies of spherical droplet flames have examined the droplet burning history of alkanes, alcohols and aromatics. Data are typically the evolution of droplet, flame, extinction, and soot shell diameters. These data are only now just beginning to find their way into comprehensive numerical models of droplet combustion to test proposed oxidation schemes for fuels such as methanol and heptane. In the present study, we report new measurements on the burning history of unsupported nonane droplets in a convection-free environment to promote spherical symmetry. The far-field gas is atmospheric pressure air at room temperature. The evolution of droplet diameter was measured using high speed cine photography of a spark-ignited, droplet within a confined volume in a drop tower. The initial droplet diameters varied between 0.5 mm and 0.6 mm. The challenge of unsupported droplets is to form, deploy and ignite them with minimal disturbance, and then to keep them in the camera field of view. Because of the difficulty of this undertaking, more sophisticated diagnostics for studying soot than photographic were not used. Supporting the test droplet by a fiber fixes the droplet position but the fiber can perturb the burning process especially for a sooting fuel. Prior studies on heptane showed little evidence for soot formation due to g-droplets of similar size the relationship between sooting and droplet diameter. For nonane droplets we expect increased sooting due to the greater number of carbon atoms. As a sooting droplet burns and its diameter decreases, proportionally less soot should form. This reduced soot, as well as the influence of soot formed earlier in the burning process which collects in a 'shell', on heat transport to the flame offers the potential for a time-varying burning rate. Such an effect was investigated and revealed in results reported here. Speculation is offered for the cause of this effect and its possible relation to soot formation.

Callahan, B. J.↗
DOE OSTI2015

Phosphorus and boron diffusion paths in polycrystalline silicon gate of a trench-type three-dimensional metal-oxide-semiconductor field effect transistor investigated by atom probe tomography

The dopant (P and B) diffusion path in n- and p-types polycrystalline-Si gates of trench-type three-dimensional (3D) metal-oxide-semiconductor field-effect transistors (MOSFETs) were investigated using atom probe tomography, based on the annealing time dependence of the dopant distribution at 900 °C. Remarkable differences were observed between P and B diffusion behavior. In the initial stage of diffusion, P atoms diffuse into deeper regions from the implanted region along grain boundaries in the n-type polycrystalline-Si gate. With longer annealing times, segregation of P on the grain boundaries was observed; however, few P atoms were observed within the large grains or on the gate/gate oxide interface distant from grain boundaries. These results indicate that P atoms diffuse along grain boundaries much faster than through the bulk or along the gate/gate oxide interface. On the other hand, in the p-type polycrystalline-Si gate, segregation of B was observed only at the initial stage of diffusion. After further annealing, the B atoms became uniformly distributed, and no clear segregation of B was observed. Therefore, B atoms diffuse not only along the grain boundary but also through the bulk. Furthermore, B atoms diffused deeper than P atoms along the grain boundaries under the same annealing conditions. This information on the diffusion behavior of P and B is essential for optimizing annealing conditions in order to control the P and B distributions in the polycrystalline-Si gates of trench-type 3D MOSFETs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Bound-preserving discontinuous Galerkin methods for conservative phase space advection in curvilinear coordinates

In this study, we extend the positivity-preserving method of Zhang & Shu 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 stability-preserving, 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\in[0,1]$. The combination of suitable CFL conditions and the use of the high-order limiter proposed in Zhang & Shu (2010) is sufficient 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 divergence-free property of the phase space flow. 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↗
DOE OSTI2024

An implicit solution for Asay foil trajectories generated by separable, sustained-production ejecta source models

We present a simple implicit solution for the time-dependent trajectory of a thin Asay foil ejecta diagnostic for the general case where the impinging ejecta cloud is generated by a source function characterized by an arbitrary (sustained) time dependence and a time-independent (stationary) particle velocity distribution. In the limit that the source function time dependence becomes a delta function, this solution—which is amenable to rapid numerical calculations of arbitrary accuracy—exactly recovers a previously published solution for the special case of instantaneous ejecta production. We also derive simple expressions for the free-surface arrival (catch-up) time as well as the true ejecta areal mass accumulation on the accelerating foil and place bounds on the level of error incurred when applying instant-production mass solutions to a sustained-production trajectory. We demonstrate these solutions with example calculations for hypothetical source functions spanning a wide range of ejecta production durations, velocity distributions, and temporal behaviors. These calculations demonstrate how the foil trajectory is often insensitive to the temporal dependence of the source function, instead being dominated by the velocity distribution. We quantify this insensitivity using a “compatibility score” metric. Under certain conditions, one may capitalize upon this insensitivity to obtain a good approximation of the second integral of the velocity distribution from the observed foil trajectory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Modeling of H2 Dispersion at ARIES

Hydrogen is a versatile and clean energy carrier that can be produced from various renewable sources such as wind, solar, and hydropower and help decarbonize electricity grids, industry, and transportation. Using the Hydrogen Research Facility under Advanced Research on Integrated Energy Systems (ARIES) at the National Renewable Energy Laboratory's (NREL) Flatirons campus as a test bench, the study examines the feasibility, useability, and value of using computational fluid dynamics (CFD) techniques to model hydrogen dispersion. The ARIES facility was chosen because controlled hydrogen releases can be performed at a rate of 27 kg-H2/hr. Site-specific atmospheric and weather condition data such as wind speed and temperature were used as inputs to the model. The results show statistical distributions and ranges of hydrogen concentrations at locations throughout the domain. Wind conditions are found to significantly impact the release behavior, including the hydrogen cloud's direction and concentrations. At low wind speeds (below 1 mph), hydrogen forms a cloud and at higher wind speeds (> 2-4 mph) hydrogen plume stretches in the direction of wind momentum. From >100 simulations for ARIES site-specific conditions, statistical quantities combined with a clustering algorithm were used to propose sensor location at various elevations from ground.

dispersion↗
DOE OSTI2026

Plutonium Retention by Crystalline Silicotitanate under Hyperalkaline Conditions Relevant to Tank-Side Cesium-Removal at the Hanford Site

Crystalline silicotitanate (CST) is used in Hanford’s Tank-Side Cesium-Removal (TSCR) process to selectively remove Cs-137 from highly caustic, nitrate-rich tank supernatants. Recent testing with actual waste samples suggests that CST can also retain measurable plutonium (Pu), which could affect radiological classification and disposal pathways for spent CST. To quantify this behavior, Pu partitioning to CST was studied under Hanford-relevant conditions using batch-contact experiments in a representative simulant (2 M NaNO3, 0.7 M NaOH). Isotherm data were measured and distribution ratios calculated, with Cs+ uptake used as benchmark. Under low-carbonate conditions, Pu was retained strongly by CST in systems initially contacted with either PuO2 nanoparticles (Pu(IV)) or aqueous Pu(VI), with distribution ratios of ~2,200–3,700 mL/g, generally exceeding those for Cs+ (~400–1,000 mL/g). Increasing carbonate concentration strongly reduced PuO2 nanoparticle retention; at [Na2CO3] = 1 M, distribution ratios decreased by up to one order of magnitude to roughly 100–300 mL/g. Electron microscopy suggests that Pu retention involves a combination of mechanisms such as PuO2 NP aggregation induced by CST leachate components, and association with CST bead surfaces.

Neumann, J.↗
DOE OSTI2017

Interplay of electron and proton instabilities in expanding solar wind

Protons and electrons observed in the solar wind possess temperature anisotropies for which upper and lower bounds appear to be partially regulated by marginal conditions associated with various kinetic plasma instabilities. Such features are most clearly seen when a collection of measurements is plotted as a two-dimensional histogram in (β{sub ∥},T{sub ⊥}/T{sub ∥}) phase space. While the partial outer boundaries of such data distribution may well be explained by various instability threshold conditions, an outstanding issue is that the majority of data points are actually located sufficiently away from the boundaries and reside in near isotropic conditions. This implies that certain processes are operative that counteract the adiabatic effect in the radially expanding solar wind, without which solar wind plasma will inexorably be forced to proceed toward the marginal firehose condition. A number of physical processes have been proposed in the literature to explain such a feature. The present paper suggests yet another mechanism. It considers dynamic electrons and protons in the quasilinear evolution of anisotropy-driven instabilities, which is in contrast to previous studies where either protons or electrons are assumed to be stationary when considering the dynamics of the other particle species. It is shown that the dynamical interplay between the two species during the quasilinear development of parallel electron firehose and proton–cyclotron instabilities leads to a counter-balancing effect, which prevents the uniform progression of the solar wind protons toward the marginal firehose state.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

Origin, transport, and losses of energetic He(+) and He(2+) ions in the magnetosphere of the Earth - AMPTE/CCE observations

Data from the ion charge-energy-mass spectrometer CHEM flown on AMPTE/CCE spacecraft are used to investigate the origin, transport, and losses of energetic He(+) and He(2+) ions in the earth's magnetosphere. The L profiles of the average ion phase space density f were determined as a function of the magnetic momentum. It is shown that the L profiles have an inner part, where f increases with L for both He(+) adn He(2+) and where steady-state conditions are fulfilled. The outer boundary L(lim) of this region is located at a distance that depends on the ion species and the geomagnetic activity level. Steady-state conditions continue outside L(lim) for He(+) ions, while the He(2+) ion distribution outside L(lim) is strongly influenced by ion convection causing a lack of steady-state conditions. It is concluded that solar wind is the origin of the He(2+), while a mixed origin is suggested for the He(+) ions, in which the major contribution is from the solar wind via charge exchange production from the He(2+) ions.

Kremser, G.↗
DOE OSTI2015

Numerical calculation of ion runaway distributions

Ions accelerated by electric fields (so-called runaway ions) in plasmas may explain observations in solar flares and fusion experiments; however, limitations of previous analytic work have prevented definite conclusions. In this work, we describe a numerical solver of the 2D non-relativistic linearized Fokker-Planck equation for ions. It solves the initial value problem in velocity space with a spectral-Eulerian discretization scheme, allowing arbitrary plasma composition and time-varying electric fields and background plasma parameters. The numerical ion distribution function is then used to consider the conditions for runaway ion acceleration in solar flares and tokamak plasmas. Typical time scales and electric fields required for ion acceleration are determined for various plasma compositions, ion species, and temperatures, and the potential for excitation of toroidal Alfvén eigenmodes during tokamak disruptions is considered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Solving high-dimensional inverse problems using amortized likelihood-free inference with noisy and incomplete data

Here, we present a likelihood-free probabilistic inversion method based on normalizing flows for high-dimensional inverse problems. The proposed method is composed of two complementary networks: a summary network for data compression and an inference network for parameter estimation. The summary network encodes raw observations into a fixed-size vector of summary features, while the inference network generates samples of the approximate posterior distribution of the model parameters based on these summary features. The posterior samples are produced in a deep generative fashion by sampling from a latent Gaussian distribution and passing these samples through an invertible transformation. We construct this invertible transformation by sequentially alternating conditional invertible neural network and conditional neural spline flow layers. The summary and inference networks are trained simultaneously. We apply the proposed method to an inversion problem in groundwater hydrology to estimate the posterior distribution of the log-conductivity field conditioned on spatially sparse time-series observations of the system’s hydraulic head responses. The conductivity field is represented with 706 degrees of freedom in the considered problem. Comparison with the likelihood-based iterative ensemble smoother PEST-IES method demonstrates that the proposed method accurately estimates the parameter posterior distribution and the observations’ predictive posterior distribution at a fraction of the inference time of PEST-IES.

conditional invertible neural network↗
DOE OSTI2018

Revisiting alkaline aerobic lignin oxidation

Lignin conversion to renewable chemicals is a promising means to improve the economic viability of lignocellulosic biorefineries. Alkaline aerobic oxidation of lignin has long been employed for production of aromatic compounds such as vanillin and syringaldehyde, but this approach primarily focuses on condensed substrates such as Kraft lignin and lignosulfonates. Conversely, emerging lignocellulosic biorefinery schemes enable the production of more native-like, reactive lignin. Here, we revisit alkaline aerobic oxidation of highly reactive lignin substrates to understand the impact of reaction conditions and catalyst choice on product yield and distribution. The oxidation of native poplar lignin was studied as a function of temperature, NaOH loading, reaction time, and oxygen partial pressure. Besides vanillin and syringaldehyde, other oxidation products include acetosyringone and vanillic, syringic, and p-hydroxybenzoic acids. Reactions with vanillin and syringaldehyde indicated that these compounds are further oxidized to non-aromatic carboxylic acids during alkaline aerobic oxidation, with syringaldehyde being substantially more reactive than vanillin. The production of phenolic compounds from lignin is favored by high NaOH loadings and temperatures, but short reaction times, as the products degrade rapidly, which is further exacerbated by the presence of oxygen. Under optimal conditions, a phenolic monomer yield of 30 wt% was obtained from poplar lignin. Testing a range of catalysts showed that Cu-containing catalysts, such as CuSO 4 and LaMn 0.8 Cu 0.2 O 3 , accelerate product formation; specifically, the catalyst does not increase the maximum yield, but expands the operating window in which high product yields are obtainable. We also demonstrate that other native and isolated lignin substrates that are significantly chemically modified are effectively converted to phenolic compounds. Finally, alkaline aerobic oxidation of native lignins was compared to nitrobenzene oxidation and reductive catalytic fractionation, as these methods constitute suitable benchmarks for lignin depolymerization. While nitrobenzene oxidation achieved a somewhat higher yield, similar monomer yields were obtained through RCF and alkaline aerobic oxidation, especially for lignins with a high guaiacyl- and/or p-hydroxyphenyl-content, as syringyl units are more unstable during oxidation. Overall, this study highlights the potential for aerobic lignin oxidation revisited on native-like lignin substrates.

09 BIOMASS FUELS↗
NASA NTRS2002

Use of Traffic Intent Information by Autonomous Aircraft in Constrained Operations

This paper presents findings of a research study designed to provide insight into the issue of intent information exchange in constrained en-route air-traffic operations and its effect on pilot decision-making and flight performance. The piloted simulation was conducted in the Air Traffic Operations Laboratory at the NASA Langley Research Center. Two operational modes for autonomous flight management were compared under conditions of low and high operational complexity (traffic and airspace hazard density). The tactical mode was characterized primarily by the use of traffic state data for conflict detection and resolution and a manual approach to meeting operational constraints. The strategic mode involved the combined use of traffic state and intent information, provided the pilot an additional level of alerting, and allowed an automated approach to meeting operational constraints. Operational constraints applied in the experiment included separation assurance, schedule adherence, airspace hazard avoidance, flight efficiency, and passenger comfort. The strategic operational mode was found to be effective in reducing unnecessary maneuvering in conflict situations where the intruder's intended maneuvers would resolve the conflict. Conditions of high operational complexity and vertical maneuvering resulted in increased proliferation of conflicts, but both operational modes exhibited characteristics of stability based on observed conflict proliferation rates of less than 30 percent. Scenario case studies illustrated the need for maneuver flight restrictions to prevent the creation of new conflicts through maneuvering and the need for an improved user interface design that appropriately focuses the pilot's attention on conflict prevention information. Pilot real-time assessment of maximum workload indicated minimal sensitivity to operational complexity, providing further evidence that pilot workload is not the limiting factor for feasibility of an en-route distributed traffic management system, even under highly constrained conditions.

Wing, David J.↗
DOE OSTI2017

Relationship between T-invariance and formation mechanisms for T-odd and T-even asymmetries in the angular distributions of fission fragments of oriented nuclei

The T-invariance condition was analyzed for the amplitude T{sub b,a} of multiparticle multistep elastic or inelastic a → b nuclear reactions. This condition leads to the equality of the amplitude T{sub b,a} to the amplitude T̃{sub ā,b̄} of b̄→ā time-reversed reaction, for which the reaction operator T̃ coincides with the inverse-reaction (b → a) operator. It is shown that, in the case where the original, inverse, and time-reversed reactions are governed by a common T-invariant mechanism, the coefficients D of asymmetries appearing in the differential cross sections for these reactions can be represented in the form of a unified scalar (pseudoscalar) function of arguments equal to the momentum and spin vectors of particles of the initial and final channels of the reactions under analysis. It is also shown that the use of the T-invariance condition for the coefficients D of asymmetries in the differential cross section for the original nuclear reaction that are different in P- and T-parity makes it possible to separate mechanisms leading to nonzero coefficients D for a number of the asymmetries under analysis from the remaining mechanisms leading to zero coefficients D of these asymmetries. It is proven that there exist such asymmetries in the differential cross section for the original reaction whose coefficients vanish for all possible T-invariant mechanisms of their appearance, so that, upon proving experimentally the appearance of nonzero coefficients of these asymmetries in the differential cross section for the original reaction, this fact can be used to assess features of T-noninvariant interactions in nuclear systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

An Alternative Frictional Boundary Condition for Computational Fluid Dynamics Simulation of Friction Stir Welding

For better application of numerical simulation in optimization and design of friction stir welding (FSW), this paper presents a new frictional boundary condition at the tool/workpiece interface for computational fluid dynamics (CFD) modeling of FSW. The proposed boundary condition is based on an implementation of the Coulomb friction model. Using the new boundary condition, the CFD simulation yields non-uniform distribution of contact state over the tool/workpiece interface, as validated by the experimental weld macrostructure. It is found that interfacial sticking state is present over large area at the tool-workpiece interface, while significant interfacial sliding occurs at the shoulder periphery, the lower part of pin side, and the periphery of pin bottom. Due to the interfacial sticking, a rotating flow zone is found under the shoulder, in which fast circular motion occurs. The diameter of the rotating flow zone is smaller than the shoulder diameter, which is attributed to the presence of the interfacial sliding at the shoulder periphery. For the simulated welding condition, the heat generation due to friction and plastic deformation makes up 54.4 and 45.6% of the total heat generation rate, respectively. In conclusion, the simulated temperature field is validated by the good agreement to the experimental measurements.

36 MATERIALS SCIENCE↗
NASA NTRS1978

Geoid Anomalies and the Near-Surface Dipole Distribution of Mass

Although geoid or surface gravity anomalies cannot be uniquely related to an interior distribution of mass, they can be related to a surface mass distribution. However, over horizontal distances greater than about 100 km, the condition of isostatic equilibrium above the asthenosphere is a good approximation and the total mass per unit column is zero. Thus the surface distribution of mass is also zero. For this case we show that the surface gravitational potential anomaly can be uniquely related to a surface dipole distribution of mass. Variations in the thickness of the crust and lithosphere can be expected to produce undulations in the geoid.

Turcotte, D. L.↗